[Paper Review] 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 Applications38 references4 citations

TL;DR

This paper presents a rigorous semi-analytical Floquet analysis for zero-thickness space-time modulated Huygens' metasurfaces using Generalized Sheet Transition Conditions (GSTCs) and Lorentzian surface susceptibilities to model temporal dispersion. The method computes scattered field harmonics via numerically solved linear equations, enabling efficient steady-state analysis of oblique plane wave and Gaussian beam excitations, with demonstrated non-reciprocity under traveling wave modulation.

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.

Motivation & Objective

• To develop an efficient, semi-analytical method for computing the steady-state response of space-time modulated Huygens’ metasurfaces under periodic excitation.
• To incorporate physically motivated Lorentzian dispersion in surface electric and magnetic susceptibilities to accurately model resonant behavior under time and space modulation.
• To validate the method against FDFD and FDTD simulations for pure-space and pure-time modulation cases.
• To demonstrate the breaking of Lorentz reciprocity in space-time modulated metasurfaces using traveling wave perturbations.
• To extend the method to arbitrary incident fields, such as Gaussian beams, via Fourier decomposition into plane wave components.

Proposed method

• The metasurface is modeled as a zero-thickness discontinuity using Generalized Sheet Transition Conditions (GSTCs), relating tangential field discontinuities to surface electric and magnetic currents.
• Surface susceptibilities χₑ and χₘ are described by Lorentzian profiles with time- and space-varying parameters (e.g., resonant frequency, damping), enabling accurate modeling of temporal dispersion.
• The scattered fields are expanded into space-time Floquet harmonics, with amplitudes determined by solving a large system of linear equations derived from the GSTC formulation.
• The method maps harmonic indices (m, n) to a single index p to form a finite matrix system, enabling numerical solution of the unknown field amplitudes.
• The matrix formulation includes terms for electric and magnetic susceptibilities, incident field components, and coupling between harmonics, with proper handling of phase and impedance matching.
• The approach is extended to Gaussian beams by decomposing them into plane wave spectra, each of which is analyzed using the same Floquet-based framework.

Experimental results

Research questions

• RQ1 How can the steady-state scattered fields of a space-time modulated Huygens’ metasurface be efficiently computed when the modulation is periodic in both space and time?
• RQ2 To what extent does incorporating Lorentzian dispersion in the surface susceptibilities improve the accuracy of field predictions in time-varying metasurfaces?
• RQ3 Can the proposed method accurately reproduce results from established FDFD and FDTD simulations for pure-space and pure-time modulation cases?
• RQ4 Does a space-time modulated Huygens’ metasurface exhibit non-reciprocal behavior under traveling wave modulation, and how is this captured in the Floquet analysis?
• RQ5 How can the method be generalized to handle arbitrary incident fields such as Gaussian beams while preserving computational efficiency?

Key findings

• The proposed method successfully validates against FDFD simulations for pure-space modulation and FDTD simulations for pure-time modulation, confirming accuracy and consistency.
• A standing wave perturbation in the metasurface parameters results in a reciprocal response, as expected from time-reversal symmetry.
• A traveling wave perturbation breaks Lorentz reciprocity, demonstrated by asymmetric scattering patterns between forward and backward propagation directions.
• The method enables efficient computation of scattered fields across multiple Floquet harmonics, with the strongest harmonic components observed at frequencies ω₀ ± ωₚ and 2ω₀ ± ωₚ.
• The scattered field distribution for a Gaussian beam incident on a modulated metasurface shows significant energy transfer to higher-order harmonics, with peak amplitudes at 1.0ω₀, 1.1ω₀, and 1.2ω₀ for different harmonic orders.
• The matrix formulation allows for scalable and numerically stable solution of the harmonic amplitudes, with convergence observed across multiple simulation cases using M=8 and N=8 harmonic modes.