[论文解读] Finite-Difference Time-Domain (FDTD) Modelling of Space-Time Modulated Metasurfaces

Scott A. Stewart, T. Smy+1 more|arXiv (Cornell University)|Dec 6, 2016

Metamaterials and Metasurfaces Applications被引用 9

一句话总结

本文提出一种有限-difference time-domain (FDTD) 方法，用于利用广义面过渡条件 (GSTCs) 对有限尺寸、零厚度的时空调制赫兹超表面进行建模。通过将全介质单元的时空介电调制映射为洛伦兹形式的等效电与磁极化率，作者推导出具有非恒定系数的二阶微分方程，并通过显式 FDTD 方法求解，以模拟平面波与高斯波激励下产生的角向分离衍射与谐波生成。

ABSTRACT

A finite-difference time-domain (FDTD) modelling of finite-size zero thickness space-time modulated Huygens' metasurfaces based on Generalized Sheet Transition Conditions (GSTCs), is proposed and numerically demonstrated. A typical all-dielectric Huygens' unit cell is taken as an example and its material permittivity is modulated in both space and time, to emulate a travelling-type spatio-temporal perturbation on the metasurface. By mapping the permittivity variation onto the parameters of the equivalent Lorentzian electric and magnetic susceptibility densities, $\ ilde{\\chi}_\ ext{ee}$ and $\ ilde{\\chi}_\ ext{mm}$, the problem is formulated into a set of second-order differential equations in time with non-constant coefficients. The resulting field solutions are then conveniently solved using an explicit finite-difference technique and integrated with a Yee-cell based propagation region to visualize the scattered fields taking into account the various diffractive effects from the metasurface of finite size. Several examples are shown for both linear and space-time varying metasurfaces which are excited with normally incident plane and Gaussian beams, showing detailed scattering field solutions. While the time-modulated metasurface leads to the generation of new collinearly propagating temporal harmonics, these harmonics are angularly separated in space, when an additional space modulation is introduced in the metasurface.

研究动机与目标

• 开发一种严格的时域 FDTD 框架，用于对具有精确表面不连续性处理的有限尺寸时空调制超表面进行仿真。
• 利用从物理全介质单元单元推导出的等效洛伦兹形式电与磁极化率密度来建模赫兹超表面。
• 展示在同时存在空间与时间调制的超表面中，谐波生成与空间-频谱分解现象。
• 将基于 GSTC 的场方程与基于 Yee 单元的 FDTD 求解器相结合，以准确可视化散射场与衍射效应。
• 通过平面波与高斯波正入射下的数值实例验证该方法，展示其在线性与时空调制响应下的表现。

提出的方法

• 利用广义面过渡条件 (GSTCs) 将超表面建模为零厚度边界，以强制满足切向场连续性及表面电流的跃变条件。
• 从实际的全介质赫兹单元中推导出电与磁表面极化率 $\chi_{\text{ee}}$ 与 $\chi_{\text{mm}}$，并以洛伦兹形式表达，以捕捉其频率相关的色散特性。
• 将单元的介电常数的空间与时间调制映射为洛伦兹极化率模型中的时变与空间变参数，从而在控制方程中引入非恒定系数。
• 将所得具有非恒定系数的二阶微分方程组重新表述为适用于显式时域有限差分积分的矩阵形式。
• 采用 Yee 非结构化网格实现 FDTD 算法，以求解传播区域中的电场与磁场，边界条件通过 GSTC 矩阵公式施加。
• 该方法可实现对散射场的全波仿真，包括由时空调制引起的谐波生成与衍射级次的角向分离。

实验结果

研究问题

• RQ1 如何利用严格的电磁框架，在时域中精确建模时空调制赫兹超表面？
• RQ2 同时存在空间与时间调制对散射场的角向分布与频谱内容有何影响？
• RQ3 采用 GSTC 与可变极化率模型的 FDTD 方法，如何捕捉有限尺寸超表面中的谐波生成与波前整形？
• RQ4 有效极化率的洛伦兹色散模型在多大程度上能准确表征真实全介质单元在动态调制下的电磁响应？
• RQ5 在平面波与高斯波激励下，线性（时不变）与时空调制超表面在散射场方面有何关键差异？

主要发现

• 时变调制的超表面由于在入射波时间周期尺度上的周期性调制，会产生新的共线传播的时域谐波。
• 当引入空间调制与时间调制并存时，所产生的谐波在空间上实现角向分离，从而实现空间-频谱复用。
• FDTD 方法成功捕捉了有限尺寸超表面的复杂衍射效应，包括波束转向与谐波束分裂。
• 仿真结果表明，谐波的角向分离程度随调制包络的空间频率增加而增大。
• 该方法准确再现了赫兹超表面在动态调制下预期的波前整形行为，验证了等效极化率模型的有效性。
• 针对高斯波与平面波激励的数值结果证实了 FDTD 公式在建模时空调制超表面瞬态与稳态响应方面的鲁棒性。