Epsilon-Near-Zero (ENZ) Modes: Physics & Applications

• Epsilon-near-zero (ENZ) modes are electromagnetic eigenmodes occurring when the real part of a material's permittivity vanishes, leading to uniform phase and dramatic field enhancement.
• ENZ conditions enable ultrafast nonlinear optical effects, such as a sixfold increase in the Kerr nonlinearity and efficient low-power all-optical switching.
• These modes are exploited in hybrid structures and nanodevices to achieve advanced photonic routing, enhanced quantum emission, and precise control over material order parameters.

Epsilon-near-zero (ENZ) modes are electromagnetic eigenmodes occurring in media whose real part of the permittivity vanishes at a particular frequency. They yield uniquely strong and tunable light–matter interaction regimes relevant for nonlinear optics, ultrafast switching, quantum photonics, and active material control. ENZ conditions arise generically from intrinsic phonon resonances, free-carrier plasmas, or engineered cavities, and are characterized by field enhancement, suppressed phase accumulation, slow group velocity, and dramatic sensitivity of material order parameters to optical excitation (Kwaaitaal et al., 2023).

1. Definition and Physical Origin

The ENZ condition is attained when the complex dielectric permittivity $\varepsilon(\omega)$ of a medium satisfies $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$, with $\mathrm{Im}\,\varepsilon(\omega_{ENZ})$ remaining low. In polar ionic crystals, this zero-crossing typically occurs near the longitudinal-optic (LO) phonon due to collective ionic displacement screening, as described by a multi-oscillator Drude–Lorentz model: $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$ Here, $\varepsilon_\infty$ is the high-frequency permittivity, $f_j$ the oscillator strength, $\omega_{p,j}$ the plasma frequency, $\omega_{0,j}$ and $\gamma_j$ the resonance and damping, respectively (Kwaaitaal et al., 2023). As $\omega\rightarrow\omega_{LO,j}$, the material enters the ENZ regime.

In conducting oxides and III–V semiconductors (e.g., ITO, doped CdO, GaN), ENZ behavior is achieved near the screened plasma frequency, while in metamaterial composites, the condition can be engineered geometrically (Inglés-Cerrillo et al., 2 Dec 2025). In two-dimensional stacks or metasurfaces, ENZ regimes can also be realized via the macroscopic Dirac point of coupled plasmonic layers (Mattheakis et al., 2016).

2. Dispersion Relations, Field Enhancement, and Phase Structure

In a homogeneous ENZ medium (relative permeability $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$0), the Maxwellian dispersion reads

$\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$1

As $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$2, $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$3, implying infinite phase velocity, uniform internal phase, and a diverging effective wavelength. For a slab at normal incidence, continuity of the displacement $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$4 enforces an internal electric field enhancement: $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$5 thus, $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$6 as $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$7 (assuming low loss). This field enhancement is responsible for the pronounced nonlinear optical responses and is a robust feature of the ENZ regime (Kwaaitaal et al., 2023, Caspani et al., 2016).

The internal field profile is nearly spatially uniform due to the phase-frozen nature ($\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$8), supporting long-range phase correlations and facilitating phenomena such as tunneling without phase delay and impedance-matched transmission at pseudo-Brewster angles (Ceglia et al., 2013). In finite systems, leaky or radiative ENZ modes can occur, often featuring slow group velocities and strong local density of states enhancement (Liberal et al., 2015).

3. Nonlinear and Ultrafast Light–Matter Interaction

ENZ media exhibit a dramatic amplification of optical nonlinearities. Third-order nonlinear polarization in the bulk, $\mathrm{Re}\,\varepsilon(\omega_{ENZ})=0$9, is strongly boosted by the internal field enhancement at ENZ. The effective nonlinear susceptibility scales as $\mathrm{Im}\,\varepsilon(\omega_{ENZ})$0 (Kwaaitaal et al., 2023, Caspani et al., 2016). This enables orders-of-magnitude increases in the Kerr nonlinearity ($\mathrm{Im}\,\varepsilon(\omega_{ENZ})$1), leading to large, ultrafast refractive index changes at moderate pump intensities (Caspani et al., 2016). Experimentally, sixfold enhancement of $\mathrm{Im}\,\varepsilon(\omega_{ENZ})$2 and ultrafast, unity-level index modulation (with sub-picosecond response) have been demonstrated in ENZ thin films (Caspani et al., 2016).

Because the group velocity vanishes ($\mathrm{Im}\,\varepsilon(\omega_{ENZ})$3), the temporal interaction of light with matter is further extended, reinforcing the efficiency of nonlinear processes. Picosecond and sub-picosecond pump pulses tuned to the ENZ wavelength have been shown to induce ultrafast, transient switching of order parameters such as ferroelectric polarization (Kwaaitaal et al., 2023). The combination of field and density-of-states enhancement is directly exploited in ultrafast optical switches, frequency conversion, and low-threshold bistability (Ceglia et al., 2013).

4. Hybrid and Coupled ENZ Structures

ENZ modes are highly sensitive to strong coupling and hybridization with other resonances, such as cavity modes, surface plasmon polaritons (SPPs), and dielectric antenna Mie resonances.

• Photonic Gap Antennas (PGAs) with embedded ENZ films display hybrid modes resulting from the interaction between the ENZ resonance of the film and the Mie resonances of the high-index dielectric. These hybrids support near-field enhancements up to $\mathrm{Im}\,\varepsilon(\omega_{ENZ})$4 and spontaneous emission (Purcell) factors exceeding $\mathrm{Im}\,\varepsilon(\omega_{ENZ})$5; field symmetry breaking yields unidirectional radiation (Patri et al., 2021, Thouin et al., 3 Feb 2025).
• Metal–Insulator–Metal (MIM) Nanocavities exhibit hybridized ENZ modes with tunable mode splitting and high $\mathrm{Im}\,\varepsilon(\omega_{ENZ})$6-factors, explained via both semi-classical double-well models and coupled-oscillator pictures. Such hybridization can achieve large Rabi-splitting, strong field confinement, and efficient light–matter coupling. Excitation via resonant tunneling obviates the momentum-matching requirement typical for SPPs (Caligiuri et al., 2020).
• Strong Coupling with Optical Cavities (e.g., unpatterned Fabry–Perot cavities) enables ENZ dispersion engineering, transforming flat ENZ bands into dispersive polaritonic branches. Such designs facilitate wide-angle polarizers and actively switchable mid-IR optical elements when integrated with phase-change materials (Johns, 2023).
• Plasmon–ENZ Hybridization in Doped Oxides (e.g., CdO, GaN) combines the field confinement of ENZ with the propagation length of SPPs, opening a route to tunable strong coupling and long-range, subwavelength guided modes in the mid-IR (Runnerstrom et al., 2018, Inglés-Cerrillo et al., 2 Dec 2025).

In ultrathin, transdimensional films, quantum confinement and nonlocal effects (e.g., Keldysh–Rytova interaction) lift plasmon mode degeneracies and further split ENZ modes, allowing control of spontaneous emission rates via film thickness (Bondarev et al., 2019).

5. ENZ-Induced Control of Functional Order Parameters

The amplification of light–matter interaction at ENZ enables not only transient but permanent control of material order parameters:

• Permanent All-Optical Ferroelectric Switching: Experiments on single-domain BaTiO$\mathrm{Im}\,\varepsilon(\omega_{ENZ})$7 (ENZ near its LO phonons) revealed that only under conditions where both $\mathrm{Im}\,\varepsilon(\omega_{ENZ})$8 and $\mathrm{Im}\,\varepsilon(\omega_{ENZ})$9 vanish could ultrafast pump pulses induce stable, persistent $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$0 or $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$1 ferroelectric domain reversal. Off-ENZ excitation yields only transient switching (<10 ms) (Kwaaitaal et al., 2023).
• Mechanism: The pump drives large amplitude LO phonon oscillations (coordinate $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$2) leading to rectified strain $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$3 through nonlinear phonon–strain and piezoelectric coupling, which yields a displacement field $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$4 that can bias and switch the ferroelectric polarization across its energy barrier.
• Universality: This mechanism is independent of specific chemical composition and applies broadly to polar crystals (perovskites, III–V's, oxides), as the ENZ condition generically follows from LO–TO splitting in ionic lattices. Analogous ENZ-driven switching of magnetic order has been reported (Kwaaitaal et al., 2023).

This establishes ENZ-excitation as a generic pathway for ultrafast, low-energy, and in some cases, permanent control over material order (polarization, magnetization).

6. Geometric, Cavity, and Structural Effects

ENZ modes manifest unique geometric and topological electromagnetic features:

• Field Confinement and Nonlocality: In deeply subwavelength films or metacavities, ENZ conditions enable field confinement beyond the classical skin depth, limited only by electronic pressure and nonlocal response. Incorporation of hydrodynamic and quantum nonlocal models is necessary to capture resonance shifts and sharp spectral features observed in ultrathin systems (Thouin et al., 3 Feb 2025, Ceglia et al., 2013, Bondarev et al., 2019).
• Shape-Independent Nonradiating ENZ Modes: ENZ cavities support nonradiating eigenmodes at $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$5 for any cavity shape, due to spatially uniform phase and suppressed emission rates. These can be switched to radiating modes by detuning or by modulating emitter orientation, granting control over light emission properties (Liberal et al., 2015).
• Ferrell–Berreman and Leaky ENZ Modes: In multilayer and finite-thickness structures, ENZ resonances often correspond not just to zero-crossing of the effective permittivity, but to specific poles in the complex propagation constant. These modes typically feature slow-light properties, high field enhancement, and are accessible by conventional free-space excitation (Newman et al., 2015, Ceglia et al., 2013).
• Topologically Distinct Polaritonic Branches: In anisotropic (uniaxial) metasurfaces, ENZ polaritons yield a rich phase diagram with seven distinct regimes, closed and open isofrequency curves, topological transitions, and controllable phase-velocity sign. These features directly impact near-field profiles, local density of states, and directionality of spontaneous emission (Alfaro-Mozaz et al., 6 May 2025).

7. Applications and Design Guidelines

ENZ modes underlie a broad spectrum of technological and scientific applications, with concrete quantitative and practical guidelines:

Application Area	ENZ Mechanism/Feature	Reference
Nonlinear optics and all-optical switching	ENZ-enhanced $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$6, ultrafast Kerr, low-power regime	(Caspani et al., 2016, Kuttruff et al., 2020, Ceglia et al., 2013)
Quantum emission and Purcell enhancement	Local density of states + field enhancement	(Patri et al., 2021, Liberal et al., 2015, Thouin et al., 3 Feb 2025)
Ultrafast and permanent control of order	ENZ field-induced switching of polarization/magnetization	(Kwaaitaal et al., 2023)
Waveguiding and photonic routing	Phase-free propagation, canalization, anisotropic ENZ modes	(Inglés-Cerrillo et al., 2 Dec 2025, Alfaro-Mozaz et al., 6 May 2025, Runnerstrom et al., 2018)
Mid-IR and visible integration	Low-loss operation in GaN, CdO, engineered MIM stacks	(Inglés-Cerrillo et al., 2 Dec 2025, Caligiuri et al., 2020)
ENZ-modified Casimir forces	Dispersionless, repulsive contributions to fluctuation forces	(Gong et al., 2022)
Sensing	Refractive index sensitivity of ENZ confinement in fibers	(Minn et al., 2017)

Design of ENZ-based nanodevices typically involves tuning carrier density or plasma frequency (via doping or composition), control over film thickness (for field overlap and loss management), and selecting suitable cavity or hybridization strategies to achieve desired mode splitting, $$\varepsilon(\omega) = \varepsilon_\infty - \sum_j \frac{f_j\,\omega_{p,j}^2}{\omega^2 - \omega_{0,j}^2 + i\,\gamma_j\,\omega}$$7, and field localization. Losses are minimized by matching the ENZ and near-zero-index (NZI) regimes and carefully engineering material damping parameters (Inglés-Cerrillo et al., 2 Dec 2025). Nonlocal and quantum effects become critical in the ultrathin or high-field limit (Thouin et al., 3 Feb 2025, Bondarev et al., 2019).

• (Kwaaitaal et al., 2023): Epsilon-near-zero regime as the key to ultrafast control of functional properties of solids
• (Patri et al., 2021): Hybrid Epsilon-Near-Zero Modes of Photonic Gap Antennas
• (Caligiuri et al., 2020): Hybridization of Epsilon-Near-Zero Modes via Resonant Tunneling in Layered Metal/Insulator Double Nanocavities
• (Runnerstrom et al., 2018): Polaritonic hybrid-epsilon-near-zero modes: engineering strong optoelectronic coupling and dispersion in doped cadmium oxide bilayers
• (Thouin et al., 3 Feb 2025): Field-enhancement and nonlocal effects in epsilon-near-zero photonic gap antennas
• (Caspani et al., 2016): Enhanced nonlinear refractive index in epsilon-near-zero materials
• (Ceglia et al., 2013): Low-damping epsilon-near-zero slabs: nonlinear and nonlocal optical properties
• (Johns, 2023): Dispersion engineering of infrared epsilon-near-zero modes by strong coupling to optical cavities
• (Liberal et al., 2015): Nonradiating and radiating modes excited by quantum emitters in open epsilon-near-zero cavities
• (Inglés-Cerrillo et al., 2 Dec 2025): GaN mid-IR plasmonics: low-loss epsilon-near-zero modes
• (Gong et al., 2022): The Effect of Epsilon-Near-Zero (ENZ) Modes on the Casimir Interaction between Ultrathin Films
• (Mattheakis et al., 2016): Epsilon-Near-Zero behavior from plasmonic Dirac point: theory and realization using two-dimensional materials
• (Newman et al., 2015): Ferrell-Berreman modes in plasmonic epsilon-near-zero media
• (Bondarev et al., 2019): Transdimensional epsilon-near-zero modes in planar plasmonic nanostructures
• (Minn et al., 2017): Excitation of epsilon-near-zero resonance in ultra-thin indium tin oxide shell embedded nanostructured optical fiber
• (Zhou et al., 2022): Realization of broadband index-near-zero modes in nonreciprocal magneto-optical heterostructures
• (Alfaro-Mozaz et al., 6 May 2025): Generalized Epsilon-Near-Zero Polaritons on Uniaxial Metasurfaces
• (Kuttruff et al., 2020): Ultrafast all-optical switching enabled by epsilon-near-zero modes in metal-insulator nanocavities