Energy Transport Velocity in Photonic Time Crystals

Steep or near-vertical Floquet dispersion in photonic time crystals (PTCs) is often read as fast, even apparently superluminal, transport. Here, we demonstrate that this anomaly arises from modulation-driven geometric drift, not energy flow. By deriving a Maxwell-flux Hellmann-Feynman relation, we prove that the cycle-averaged energy velocity remains strictly bounded. We further establish a universal velocity-product law conserved throughout the passband, $ v_E v_g=\langle v_{\rm ph}^2\rangle_T $, fixing transport solely by the temporal average of the inverse permittivity. The divergent group velocity is then traced to a mismatch between electric and magnetic geometric phase connections, revealing apparent superluminality as a geometric effect of temporal modulation.

💡 Research Summary

The manuscript investigates the true nature of energy transport in photonic time crystals (PTCs), a class of temporally periodic, spatially homogeneous media whose permittivity ε(t) repeats with period T. While recent experiments and theoretical works have highlighted the appearance of extremely steep, almost vertical Floquet dispersion curves ω(k) in such systems—often interpreted as “super‑luminal” group velocities—the authors demonstrate that these features do not correspond to actual energy flow.

Starting from Maxwell’s equations in a one‑dimensional, source‑free medium with constant permeability μ>0 and a strictly positive, real, T‑periodic permittivity ε(t), the authors define the instantaneous electromagnetic energy density u(z,t) and Poynting flux S_z(z,t). The Poynting theorem acquires an extra term −½ \dot{ε}(t)|E|², which represents the exchange of energy between the field and the temporal modulation. By averaging over one modulation period, they introduce the cycle‑averaged energy‑transport velocity

 v_E ≡ ⟨S_z⟩_T / ⟨u⟩_T,

and note the trivial bound |v_E| ≤ v_ph,max ≤ c, where v_ph,max = 1/√(μ ε_min).

To connect v_E with the Floquet band structure, the authors expand the fields in a plane‑wave ansatz E(z,t)=E(t)e^{ikz}, H(z,t)=H(t)e^{ikz} and Fourier‑decompose both the fields and ε(t). This yields a non‑Hermitian eigenvalue problem in Sambe space, M(ω,k)|R⟩=0, where |R⟩ = (E,H)^T. Conventional Hellmann–Feynman formulas involve the left null vector ⟨L|, but the electromagnetic problem possesses a natural symplectic (flux) structure encoded in the matrix

 J =