Nonreciprocal Negative Refraction Enabled by Photonic Time Crystals

TL;DR

The paper tackles magnet-free, nonreciprocal negative refraction by embedding hyperbolic media between time-varying interfaces, which induce Floquet sidebands and break reciprocity while preserving negative refraction. It develops a general time-crystal framework and demonstrates two realizations: optical 3D time crystals using permittivity-modulated dielectric slabs around an AZO/ZnO hyperbolic stack, and microwave 2D time crystals using conductance-modulated metasurfaces around a wire medium. The key contributions include a unified theory based on Floquet-harmonic expansions, a practical optical implementation achieving >$46$ dB isolation, a microwave counterpart achieving ~$11$ dB isolation, and a Harmonic-Balance FEM simulation methodology that validates the concept across frequency regimes. The work significantly broadens the design space for magnet-free nonreciprocal components by showing how temporal modulation and hyperbolic dispersion can be combined to control forward and backward propagation in a frequency-agnostic manner. The findings point to future on-chip isolators and asymmetric beam routing devices that exploit time-varying metasurfaces and photonic time crystals across optical and microwave domains.

Abstract

We propose and theoretically demonstrate nonreciprocal negative refraction enabled by time-varying photonic structures. By engineering temporal modulations at the interfaces of hyperbolic media, we achieve isolation between forward and backward beams while preserving the hallmark property of negative refraction. Two complementary approaches are developed: in the optical regime, a multilayer AZO/ZnO hyperbolic slab is sandwiched between permittivity-modulated dielectric layers (3D time crystals); in the microwave regime, a wire medium is sandwiched between time-modulated resistive metasurfaces (2D time crystals). Both designs exploit Floquet harmonic expansions and are validated with a customized harmonic-balance finite-element solver. We report isolation exceeding 46 dB in the optical device and ~11 dB in the microwave counterpart. This work introduces a general framework for nonreciprocal negative refraction across frequency regimes, expanding the design space of time-varying metasurfaces and photonic time crystals.

Figures (5)

• Figure 1: Modeling time-varying interfaces through negative refraction. (a) Schematic of a multilayer structure where a TM-polarized wave impinges from air (blue) onto a time-varying (TV) dielectric slab of thickness $d$ (red) adjacent to a hyperbolic medium (green). The incident beam (thick white arrow) excites reflected and transmitted harmonics (black arrows). Inside the TV slab, the field is a superposition of forward and backward multi-harmonic modes, $\psi_m^{+}(t)$ and $\psi_m^{-}(t)$, with amplitudes $a_m^{\text{TV}}$ and $b_m^{\text{TV}}$. Different harmonics of these modes are shown by the red arrows. Black and red wavy arrows with varying oscillation periods denote distinct harmonics. The structure is homogeneous along the $y$-axis. (b) Equivalent circuit model of (a). Each harmonic in the air and hyperbolic regions is represented by a transmission line with propagation constant $k_{z,m}$ and wave admittance $Y_m$. The incident excitation couples into the zeroth-order line of the air region. The TV slab is modeled (red box) as a linear combination of forward and backward harmonics, with voltage-controlled current sources mediating inter-harmonic coupling. Three dots before and after the lines indicate that only a subset of the infinite harmonic spectrum is shown. (c) Schematic of a configuration with a time-modulated conductive sheet (red line) at the air--hyperbolic interface. Unlike the TV slab in (a), the sheet generates multi-harmonic surface currents that couple the incident excitation to reflected and transmitted harmonics. Similar to (a), the structure is homogeneous along the $y$-axis. (d) Equivalent circuit model of (c). Transmission lines again represent harmonics in the air and hyperbolic regions. The time-modulated sheet is modeled (red box) by shunt admittances realized through voltage-controlled current sources, which directly couple harmonics across the interface.
• Figure 2: Working principle and typical response of the proposed nonreciprocal negative refraction in the optical range. (a) A schematic of the structure composed of an alternating layer stack of thin plasmonic and dielectric films (AZO/ZnO), sandwiched between two dielectric slabs with modulated permittivity. Besides the subwavelength features of the hyperbolic layer, the system is mainly characterized by the hyperbolic slab thickness, $L = 2.8 \ {\rm{ \mu m}}$, and the incident angle, $\theta_\mathrm{inc} = 40^\circ$. The permittivities of the two modulated slabs are varied with the same frequency and modulation depth but have different phases, i.e., $\phi_1$ and $\phi_2$. The forward and backward oblique beams exhibit significantly different transmissions, resulting in nonreciprocal behavior. The red light beams schematically depict diffraction (due to harmonic modulations) and refraction through the structure. (b) The real and imaginary parts of the effective perpendicular and transverse permittivities of the thin film stack, used as the hyperbolic medium to realize negative refraction. The periodicity of the thin films along the structure is subwavelength, with a metal filling factor of $f = t_m / (t_m + t_d) = 0.3$, where $t_m$ and $t_d$ are the thicknesses of the metal and dielectric layers, respectively. (c) Forward and backward transmission and reflection field amplitudes, obtained through simulation and analytical calculation, versus harmonic order, $m$, with each frequency harmonic $\omega_m$ being $\omega_0 + m\Omega$. (d) The magnetic field profile ($Hy$) when the structure is illuminated from the top (forward) and bottom (backward). Both surface plots use the same linear scale, shown via the colorbar on the right. (e) Power vector fields corresponding to the field profile in (d). The vector fields are normalized, and the arrow colors indicate power amplitude. Both vector field plots share the same linear scale, represented by the colorbar on the right.
• Figure 3: Performance analysis of the optical design. (a) Defined figure of merit (FoM) as a function of the hyperbolic region thickness (normalized to free-space wavelength by multiplication with $k_0=\omega_0/c$) and the modulation frequency (normalized to the optical frequency $\omega_0$). The FoM (color scale) is expressed in dB. The black circle marks the configuration with the maximum FoM. (b) FoM and zeroth-order transmission as a function of frequency detuning, both in dB scale. The isolation is quantified by $I = 46.4$ dB.
• Figure 4: Conceptual illustration and representative response of the proposed approach for achieving nonreciprocal negative refraction in the microwave domain. (a) Schematic of the structure made by a wire medium sandwiched by two sheets whose conductance is modulated. Apart from the subwavelength structural parameters of the hyperbolic layer, this system is defined by only two geometrical parameters: hyperbolic medium thickness, $L$, and the incident angle $\theta_{\rm{inc}} = 40^\circ$. The driving circuit, consisting of a source with a modulation frequency of $\Omega$, a power divider (PD), and a phase shifter (PS), provides temporal modulations with different phases. The forward and backward oblique beams exhibit very different transmissions, resulting in nonreciprocity. The red beams (wavy arrows) qualitatively indicate diffraction induced by temporal harmonics and refraction through the structure. (b) Extracted real and imaginary components of the effective transverse permittivity of the wire medium, serving as the hyperbolic layer to enable negative refraction. The lattice constant is $a = 4$ mm along both x and y, with wire radius $r = 0.2$ mm. (c) Simulated and analytically computed amplitudes of the transmitted and reflected fields for both illumination directions, plotted against harmonic index $m$, corresponding to frequency components $\omega_m = \omega_0 + m\Omega$. (d) Spatial profiles of the magnetic field component $H_y$ for forward (top) and backward (bottom) illumination scenarios. Both plots share a common linear color scale shown on the right. (e) Time-averaged Poynting vector fields corresponding to the field maps in panel (d), where arrow direction and color denote energy flow direction and magnitude, respectively. The same logarithmic color scale is used for both cases and displayed on the right.
• Figure 5: Performance analysis of the microwave design. (a) Figure of merit (FoM) as a function of the hyperbolic region thickness, normalized through $k_0=\omega_0/c$, and the modulation frequency, normalized to the wave frequency $\omega_0$. The FoM, shown on a dB scale, highlights the parameter space leading to the optimal performance. (b) Frequency detuning dependence of the FoM and zeroth-order transmission, both in dB. The design achieves an isolation of $I = 11.4$ dB, arising from the two time-modulated sheets.