[论文解读] Floquet Analysis of Space-Time Modulated Huygens' Metasurfaces with Lorentz Dispersion

Ville Tiukuvaara, T. Smy+1 more|arXiv (Cornell University)|Jul 14, 2020

Metamaterials and Metasurfaces Applications参考文献 38被引用 4

一句话总结

本文提出了一种针对零厚度时空调制赫兹 metasurface 的严格半解析 Floquet 分析方法，采用广义面过渡条件（GSTCs）和 Lorentz 表面极化率来建模时间色散。该方法通过数值求解线性方程组计算散射场谐波，实现了对斜入射平面波和高斯光束激励下稳态响应的高效分析，并在行波调制下展示了非互易性。

ABSTRACT

A rigorous semi-analytical Floquet analysis is proposed for a zero-thickness space-time modulated Huygens' metasurface to model and determine the strengths of the new harmonic components of the scattered fields. The proposed method is based on Generalized Sheet Transition Conditions (GSTCs) treating a metasurface as a spatial discontinuity. The metasurface is described in terms of Lorentzian electric and magnetic surface susceptibilities, $χ_ ext{e}$ and $χ_ ext{m}$, respectively, with parameters (e.g. resonant frequency) that are periodically modulated in both space and time. The unknown scattered fields are expressed in terms of Floquet harmonics, for which the amplitudes can be found by numerically solving a set of linear equations, leading to the total scattered fields. Using existing computational techniques, the method is validated using several examples of pure-space and pure-time modulation with different modulation strengths and pumping frequencies. Finally, two cases of spacetime modulation (standing wave perturbation and a traveling wave perturbation) are presented to demonstrate the breaking of Lorentz reciprocity. The proposed method is simple and versatile and able to determine the steady-state response of a space-time modulated Huygen's metasurface that is excited with an oblique plane wave, or a general incident field such as a Gaussian beam.

研究动机与目标

• 开发一种高效、半解析的方法，用于计算在周期性激励下时空调制赫兹 metasurface 的稳态响应。
• 在表面电极化率和磁极化率中引入物理解释的 Lorentz 表型，以准确建模时空调制下谐振行为。
• 通过与 FDFD 和 FDTD 模拟对比，验证该方法在纯空间调制和纯时间调制情况下的准确性。
• 通过行波扰动，展示时空调制 metasurface 中 Lorentz 互易性被打破的现象。
• 通过将任意入射场（如高斯光束）分解为平面波分量，将该方法推广至任意入射场，同时保持计算效率。

提出的方法

• 将 metasurface 建模为使用广义面过渡条件（GSTCs）的零厚度不连续界面，将切向场不连续性与表面电流传导和磁流传导关联。
• 表面极化率 χₑ 和 χₘ 采用随时间和空间变化的 Lorentz 表型参数（如谐振频率、阻尼），以准确建模时间色散。
• 将散射场展开为时空 Floquet 谐波，其振幅通过从 GSTC 公式推导出的大规模线性方程组求解确定。
• 将谐波指数 (m, n) 映射为单一索引 p，形成有限矩阵系统，从而实现对未知场振幅的数值求解。
• 矩阵公式包含电极化率和磁极化率项、入射场分量以及谐波间的耦合项，并正确处理相位和阻抗匹配。
• 通过将高斯光束分解为平面波谱，每一分量均采用相同的基于 Floquet 的分析框架，将该方法扩展至高斯光束。

实验结果

研究问题

• RQ1 当调制在时空上均为周期性时，如何高效计算时空调制赫兹 metasurface 的稳态散射场？
• RQ2 在表面极化率中引入 Lorentz 色散在多大程度上提升了对时变 metasurface 场预测的准确性？
• RQ3 所提出的方法能否准确再现纯空间调制和纯时间调制情况下已有 FDFD 和 FDTD 模拟的结果？
• RQ4 在行波调制下，时空调制赫兹 metasurface 是否表现出非互易行为？该现象如何在 Floquet 分析中体现？
• RQ5 如何将该方法推广至处理任意入射场（如高斯光束）的同时保持计算效率？

主要发现

• 所提出的方法在纯空间调制情况下与 FDFD 模拟结果一致，在纯时间调制情况下与 FDTD 模拟结果一致，验证了其准确性和一致性。
• 在 metasurface 参数中引入驻波扰动后，响应保持互易，符合时间反演对称性预期。
• 引入行波扰动后，Lorentz 互易性被打破，表现为正向与反向传播方向的散射图案不对称。
• 该方法可高效计算多个 Floquet 谐波的散射场，最强的谐波分量出现在频率 ω₀ ± ωₚ 和 2ω₀ ± ωₚ 处。
• 高斯光束入射到调制 metasurface 上时，散射场分布显示显著的能量转移到高阶谐波，不同阶次的谐波峰值振幅分别出现在 1.0ω₀、1.1ω₀ 和 1.2ω₀ 处。
• 矩阵公式支持谐波振幅的可扩展且数值稳定的求解，在使用 M=8 和 N=8 谐波模式的多个仿真案例中均观察到收敛。