Hyperbolic Phonon Polaritons in Natural Crystals

• Hyperbolic phonon polaritons are hybrid electromagnetic–lattice excitations in anisotropic crystals, characterized by opposite-sign dielectric tensor components that allow high in-plane momenta and deep subwavelength confinement.
• They naturally occur in materials such as hexagonal boron nitride and α-MoO₃, enabling applications like nanocavity formation, beam steering, and enhanced local sensing.
• Substrate and geometry engineering dynamically tune their dispersion, confinement, and energy transport, paving the way for advanced mid-IR nanophotonic devices.

Searching arXiv for recent and foundational papers on hyperbolic phonon polaritons to ground the article in the literature.
Hyperbolic phonon polaritons are hybrid electromagnetic–lattice excitations sustained by polar anisotropic crystals in spectral intervals where principal components of the dielectric tensor have opposite signs. In these Reststrahlen bands, the isofrequency surfaces become open hyperboloids rather than closed ellipsoids, enabling deeply subwavelength momenta, directional energy flow, and guided or volume-confined infrared modes with strong field localization. In natural van der Waals crystals such as hexagonal boron nitride and biaxial oxides such as (\alpha)-MoO(_3), these modes appear without artificial metamaterial structuring and have been exploited for nanocavity formation, beam steering, substrate-controlled dispersion engineering, long-range energy transfer, nanofocusing, and subsurface diagnosis [2105.01148], [2604.08996], [2603.06184], [2401.09594], [1808.03660].

1. Definition, dielectric criteria, and hyperbolic dispersion

Phonon polaritons arise from coupling between infrared photons and optical phonons in polar crystals. Within the interval between transverse optical and longitudinal optical phonon frequencies, the real part of the permittivity becomes negative, allowing strongly confined electromagnetic modes. In anisotropic crystals, the permittivity is tensorial, for example
[
\boldsymbol{\varepsilon}(\omega)=
\begin{pmatrix}
\varepsilon_{xx}(\omega) & 0 & 0 \
0 & \varepsilon_{yy}(\omega) & 0 \
0 & 0 & \varepsilon_{zz}(\omega)
\end{pmatrix},
]
and hyperbolicity occurs when the real parts of at least two components have opposite signs, such as
[
\operatorname{Re}[\varepsilon_i(\omega)]\cdot \operatorname{Re}[\varepsilon_j(\omega)] < 0.
]
This sign-indefinite response converts the momentum-space isofrequency surface into a hyperboloid and permits very large in-plane wavevectors relative to the free-space value [2105.01148], [2604.08996].

For uniaxial media such as hBN, extraordinary-wave dispersion is commonly written as
[
\frac{k_x^2 + k_y^2}{\varepsilon_{\parallel}(\omega)} +
\frac{k_z^2}{\varepsilon_{\perp}(\omega)} =
\left(\frac{\omega}{c}\right)^2
]
or, equivalently under the notation used for hBN,
[
\frac{k_x^2 + k_y^2}{\varepsilon_z(\omega)} + \frac{k_z^2}{\varepsilon_t(\omega)} = \left(\frac{\omega}{c}\right)^2.
]
For biaxial systems such as (\alpha)-MoO(_3), in-plane propagation can be described approximately by
[
\frac{k_x^2}{\varepsilon_y(\omega)} + \frac{k_y^2}{\varepsilon_x(\omega)} = \left(\frac{\omega}{c}\right)^2,
]
so the orientation of the open contour depends directly on the signs and magnitudes of (\varepsilon_x) and (\varepsilon_y) [2105.01148], [2103.17016].

Several consequences follow from this dispersion. Hyperbolic phonon polaritons support large in-plane momenta (q), strong confinement, directional or ray-like propagation, and ultra-slow group velocities with long lifetimes because the underlying material excitation is phononic rather than electronic [2105.01148]. In hBN, HPhP wavelengths can reach (1)–(100) nm while free-space wavelengths are of order (10\,\mu\mathrm{m}), evidencing wavelength compression by factors of (100)–(1000) in the relevant bands [1504.05345]. In (\alpha)-MoO(_3), confinement factors around (80) were reported at (990\ \mathrm{cm}^{-1}) in nanobelt cavities [2105.01148], while plasmonic-antenna control produced (\lambda_0/\lambda_p \approx 40) at (946\ \mathrm{cm}^{-1}) [2103.17016].

A common misconception is that hyperbolicity is synonymous with merely having a negative permittivity. The literature distinguishes ordinary phonon polaritons from hyperbolic phonon polaritons by the sign opposition among tensor components, not by negativity alone. Another common simplification is to treat the dielectric tensor as dependent only on electric-field polarization. A more recent analysis argued that once LO–TO splitting is fully included, dielectric functions become wave-vector-direction dependent as well as polarization dependent, improving agreement with measured hBN permittivity and predicting unusual dumbbell-shaped and butterfly-shaped isofrequency curves in some hexagonal crystals [2402.11956].

2. Natural material platforms and Reststrahlen-band taxonomy

Natural HPhP platforms span uniaxial and biaxial polar crystals. Hexagonal boron nitride is the canonical uniaxial system, with dielectric tensor
[
\boldsymbol{\varepsilon}(\omega) =
\begin{pmatrix}
\varepsilon_{\perp}(\omega) & 0 & 0 \
0 & \varepsilon_{\perp}(\omega) & 0 \
0 & 0 & \varepsilon_{\parallel}(\omega)
\end{pmatrix}.
]
Its lower and upper Reststrahlen bands correspond to opposite sign configurations of (\operatorname{Re}\varepsilon_{\perp}) and (\operatorname{Re}\varepsilon_{\parallel}), yielding two hyperbolic windows in the mid-infrared [2604.08996], [1704.01834]. In the terminology of one hBN study, the upper Reststrahlen band (1370\text{–}1610\ \mathrm{cm}^{-1}) is a type II regime with (\varepsilon_t<0\) and \(\varepsilon_z>0) [1704.01834]. Another study reported type I and type II hyperbolic bands for hBN as (750\text{–}797\ \mathrm{cm}^{-1}) and (1349\text{–}1554\ \mathrm{cm}^{-1}), respectively [2402.11956].

(\alpha)-MoO(_3) is a biaxial van der Waals semiconductor with crystal axes (x\parallel[100]), (y\parallel[001]), and (z\parallel[010]). Because (\varepsilon_x(\omega)), (\varepsilon_y(\omega)), and (\varepsilon_z(\omega)) change sign at different phonon resonances, (\alpha)-MoO(_3) hosts multiple Reststrahlen bands across the far-IR and mid-IR. One report describes four distinct bands, RB1–RB4, featuring type I and type II behavior along different axes and frequency-dependent transitions between in-plane and out-of-plane hyperbolicity [2105.01148]. Another gives three Reststrahlen bands in the (545\text{–}1004\ \mathrm{cm}^{-1}) range, with Band 2 characterized by (\operatorname{Re}(\varepsilon_x)<0\), \(\operatorname{Re}(\varepsilon_y),\operatorname{Re}(\varepsilon_z)>0), and Band 3 by (\operatorname{Re}(\varepsilon_z)<0\), \(\operatorname{Re}(\varepsilon_x),\operatorname{Re}(\varepsilon_y)>0) [2104.00429]. These descriptions are consistent in emphasizing that (\alpha)-MoO(_3) supports both in-plane and out-of-plane hyperbolic regimes, with especially strong in-plane anisotropy [2103.17016].

Several materials comparisons recur in the literature. hBN is a natural low-loss hyperbolic crystal with two principal hyperbolic bands [2105.01148]. (\alpha)-MoO(_3) offers more spectral windows, in-plane hyperbolicity, and biaxial control [2105.01148], [2103.17016]. SiC is cited as a bulk polar crystal with a single Reststrahlen band and relatively isotropic behavior unless structured [2105.01148]. A broader survey of hexagonal compounds concluded that, besides hBN, h-AlN exhibits a wide hyperbolic frequency-band range, whereas h-BP, h-AlP, h-GaN, and h-GaP show it scarcely [2402.11956].

This suggests that the width and utility of a hyperbolic band depend not only on anisotropy, but on the separation of infrared-active TO frequencies along different directions, phonon lifetime, and overall phonon energy scale [2402.11956]. That conclusion is drawn explicitly in the h-AlN comparison study rather than inferred from general intuition.

3. Mode structure, confinement, and dispersion orders

In finite slabs, HPhPs are quantized into thickness-dependent branches or dispersion orders. For hBN slabs, one study writes the in-plane momentum of branch (l) as
[
k = k_l' + i k_l'' =
-\left[ \arctan\left(\frac{\varepsilon_t\, \Psi}{\varepsilon_z}\right) + \arctan\left(\varepsilon_s \Psi\right) \right] \frac{l\pi}{d},
]
where (d) is slab thickness, (l) is the mode order, (\varepsilon_s) is the substrate permittivity, and (k_l',k_l'') encode wavelength and damping [2601.16465]. In this framework, the (l=0) mode is the fundamental branch, while (l=1,2,\dots) are higher-order branches with incrementally higher momenta and different parity. Finite-difference frequency-domain calculations in the same work show even and odd (E_x) field profiles across thickness for different (l), establishing a modal taxonomy based on vertical nodal structure [2601.16465].

For thin hBN slabs over metal, gap engineering can transform a fundamental free-standing HPhP continuously into an image HPhP. At (1490\ \mathrm{cm}^{-1}), a (33) nm hBN slab above Au evolves from an “almost suspended” mode at (G=75) nm to a strongly confined image HPhP at (G=0) nm, with the in-plane momentum increasing by more than three times relative to the free-standing case at mid-band [2604.08996]. Experimentally, for the same (33) nm slab, (\operatorname{Re}(q)) increases by about (2.7\times) as the gap changes from (75) nm to (0) nm over (1450\text{–}1510\ \mathrm{cm}^{-1}) [2604.08996].

For (\alpha)-MoO(3), slab-guided HPhPs in lithography-free nanobelts are quantized by the belt width, producing Fabry–Perot modes across several Reststrahlen bands [2105.01148]. The cavity condition is written approximately as
[
2 k{\mathrm{pp}}(\omega) L + \phi_r(\omega) = m\pi,
]
where (L\approx w) is the effective cavity length, (\phi_r) is the edge reflection phase, and (m) is the resonance order [2105.01148]. The experimentally imaged standing-wave maxima correspond to different integer orders (n), and the extracted in-plane momenta lie on the dispersion of the fundamental guided mode (M_0) of an infinite slab of the same thickness [2105.01148].

An important correction to a common oversimplification is that different HPhP orders are not merely higher harmonics of one in-plane wave. They correspond to distinct slab eigenmodes with different field symmetry and thickness quantization [2601.16465]. Another misconception is that higher order always implies stronger environmental sensitivity. A substrate-index study found the opposite in relative terms: the absolute change (\Delta k) under substrate modification is nearly the same across HPhP orders, but the fractional change is smaller for higher orders because they begin at larger baseline (k) [1907.01950].

4. Launching, imaging, and spectroscopy

Because HPhPs possess momenta far larger than those of free-space photons, excitation and detection require momentum-matching strategies. The standard approach is scattering-type scanning near-field optical microscopy, in which a metallic AFM tip launches and detects polaritons locally. This technique has been used extensively in hBN, (\alpha)-MoO(_3), and related systems [1704.01834], [2105.01148], [2103.17016], [2604.08996], [2601.16465].

In s-SNOM, the tip acts as a nanoscale antenna, supplying high-(k) Fourier components to couple incident mid-IR radiation into HPhPs. Reflections from edges, steps, terraces, or internal defects generate interference fringes whose spacing gives the polariton wavelength. For typical tip-launched standing-wave configurations, the fringe spacing is (\lambda_p/2); for edge-launched propagation interfered with a background field, it can equal (\lambda_p) [2604.08996]. In hBN, a representative tip-launched wavelength at (1428\ \mathrm{cm}^{-1}) was reported as (460\pm 5) nm, matching a calculated value of (455) nm [1705.10318].

Several complementary spectroscopic and imaging schemes have expanded the HPhP toolbox. Broadband synchrotron infrared nanospectroscopy combined with pseudo-heterodyne s-SNOM was used to map HPhPs in (\alpha)-MoO(_3) nanobelts from (\approx 320) to (3000\ \mathrm{cm}^{-1}), while tunable QCL-based s-SNOM provided higher-SNR narrowband images for mode assignment [2105.01148]. Mechanical detection via AFM contact-mode cantilever oscillations was shown to image hBN HPhPs without detecting scattered photons, using photothermal expansion induced by locally enhanced lattice vibrations [1704.01834]. That method achieved better than (\lambda/20) resolution and reproduced the same fringe periodicity as s-NSOM on the same flakes [1704.01834].

Far-field access has historically been difficult because of the large momentum mismatch. Two distinct solutions emerge in the literature. One is electrical detection in graphene/hBN heterostructures, where split metallic gates launch HPPs in hBN that propagate as confined rays toward a graphene (pn)-junction, producing photocurrent via hot-carrier photothermoelectric conversion [1705.10318]. The other is the construction of photonic hypercrystals in (\alpha)-MoO(_3): periodic nanohole lattices provide reciprocal lattice vectors that phase-match free-space photons to in-plane HPhPs, enabling direct far-field reflectance signatures of hyperbolic modes and twist-controlled resonance tuning [2309.17146].

A further diagnostic route is purely electronic. In graphene/hBN heterostructures, coupling between graphene plasmons and hBN slab phonon polaritons produces plasmon–phonon polariton branches whose signatures appear in angle-resolved photoemission spectra as satellite features and self-energy peaks [1504.05345]. This does not detect HPhPs optically; rather, it transfers their spectral fingerprints into graphene’s electronic spectral function.

5. Substrate, cavity, and geometry control

Environmental engineering is one of the dominant mechanisms for controlling HPhPs. In hBN on corrugated gold, a sinusoidally varying gap (G(x)) from (0) to (\approx 75) nm produces continuous local tuning of the HPhP wavelength. For (33) nm hBN in the upper Reststrahlen band, (\operatorname{Re}(q)) increases by about (2.7\times) between (G=75) nm and (G=0) nm [2604.08996]. The same structure functions as a “polaritonic taper,” compressing (\lambda_p) from about (451) nm to about (180) nm at (1490\ \mathrm{cm}^{-1}), corresponding to about (2.5\times) lateral nanofocusing [2604.08996].

Substrate refractive index also provides strong control even without metallic image-mode formation. For hBN on suspended, dielectric, and metallic substrates, the thickness-normalized wavevector (kd) was shown to vary strongly with substrate permittivity. The wavevector can be reduced by a factor of (25) simply by changing the substrate from dielectric to metallic behavior [1907.01950]. Incorporating the imaginary part of the substrate dielectric function showed that small local variations in loss or carrier density can dynamically control the wavevector and may spatially separate hyperbolic modes of different orders [1907.01950].

In (\alpha)-MoO(_3), gradient suspension supplies a different environmental degree of freedom. A wedge-shaped air gap under a suspended flake modifies the Fabry–Perot phase condition for guided HPhPs. The dependence of (\lambda_p) on gap thickness is opposite in the lower and upper bands: in the in-plane hyperbolic band with (\operatorname{Re}(\varepsilon_x)<0\), \(\lambda_p\) increases with gap; in the upper band with \(\operatorname{Re}(\varepsilon_x)>0), (\lambda_p) decreases with gap [2104.00429]. The reported tuning reached polariton wavelength elongation up to (160\%) and damping-rate reduction up to (35\%) [2104.00429].

Geometry on the scale of the polariton wavelength is equally important. Naturally grown quasi-1D (\alpha)-MoO(_3) nanobelts form lithography-free Fabry–Perot cavities that support HPhP standing waves in multiple bands [2105.01148]. Geometrically designed Au antennas on (\alpha)-MoO(_3) impose a spatial phase profile on launched in-plane HPhPs, enabling frequency-dependent focusing. At (906\ \mathrm{cm}^{-1}), a convex antenna extremity yields a focal spot of about (1.3\,\mu\mathrm{m}), around (8.5\times) smaller than the free-space wavelength, and at (926\ \mathrm{cm}^{-1}) the lateral confinement reaches about (36) [2103.17016]. The focusing behavior reverses between hyperbolic and elliptic bands, which clarifies that this is not a generic antenna effect but a manifestation of in-plane hyperbolicity [2103.17016].

Mode conversion provides a related geometric control principle. In step-shaped hBN and (\alpha)-MoO(_3) terraces, simple slab edges preserve HPhP order, while asymmetric vdW steps break symmetry and supply additional scattering momentum, enabling (l=0\rightarrow 1) conversion and short-period beating patterns observed by s-SNOM [2601.16465]. In hBN terraces at (1407\ \mathrm{cm}^{-1}), the converted first-order branch produced Fourier peaks near (35\ \mu\mathrm{m}^{-1}), far above the (l=0) peaks near (3) and (6\ \mu\mathrm{m}^{-1}) [2601.16465].

6. Directionality, refraction, focusing, and transport

Directional propagation is one of the defining operational features of HPhPs. In hBN, the propagation angle for type-II hyperbolic modes is expressed in one study as
[
\tan\theta = \sqrt{\frac{\varepsilon_z}{|\varepsilon_t|}},
]
yielding fixed ray trajectories within the slab [1704.01834]. In (\alpha)-MoO(_3), the in-plane propagation direction depends strongly on frequency and band because different tensor components become negative in different spectral intervals [2105.01148], [2103.17016]. Simulations and near-field maps show propagation elongated along ([001]) in one band and along ([100]) in another, directly reflecting the open orientation of the in-plane hyperbola [2105.01148].

This directional control underlies several wave phenomena. Negative refraction was visualized at interfaces between the natural hyperbolic crystals h(^{11})BN and MoO(_3), where type-I and type-II phonon polaritons coexist in the same (740\text{–}822\ \mathrm{cm}^{-1}) spectral window [2209.15155]. At a special frequency (\omega_0), the resulting collimated rays circulate along closed diamond-shaped trajectories, and the eigenmode structure shows regions of both positive and negative dispersion separated by gaps associated with polaritonic level repulsion and strong coupling [2209.15155]. This indicates that interface engineering between different natural hyperbolic crystals can realize ray optics not available in a single medium.

Directional propagation also drives long-range interaction phenomena. A theoretical framework for (\alpha)-MoO(_3) slabs showed that hyperbolic phonon polaritons can mediate dipole–dipole interactions with enhancement factors of order (10^3)–(10^4) along hyperbolic asymptotes, extending significant coupling to (R\sim 50\,\mu\mathrm{m}), about (5\lambda_0), in the mid-IR [2603.06184]. At (900\ \mathrm{cm}^{-1}), enhancement peaks around (4200) were reported along the hyperbolic asymptotes; at (850\ \mathrm{cm}^{-1}), enhancement remained strong at about (2500); and at (825\ \mathrm{cm}^{-1}), loss-induced canalization reduced the peak to about (50) while maximizing propagation length and directional transport [2603.06184]. Twisted bilayers allowed a trade-off between enhancement and directionality, with (\theta\simeq 69.3^\circ) yielding canalization and (\theta=90^\circ) giving weaker, more isotropic coupling [2603.06184].

Thermal transport provides another example of directional and high-speed energy flow. In hBN, pump–probe thermoreflectance measurements showed that volume-confined HPhP modes excited by near-field radiation from hot Au can mediate interfacial heat transfer with an effective thermal boundary conductance of at least (500\ \mathrm{MW\,m^{-2}K^{-1}}), compared to a phonon-only Au/hBN conductance of (12\pm 2\ \mathrm{MW\,m^{-2}K^{-1}}) measured by TDTR [2401.09594]. The same study cites HPhP group velocities approaching (10^{-2}c) and a theoretical upper-limit conductance around (2.5\ \mathrm{GW\,m^{-2}K^{-1}}) [2401.09594]. A plausible implication is that HPhPs can function as an ultrafast photonic heat channel inside solids, not merely as near-field spectral resonances.

7. Diagnostics, applications, and emerging quantum directions

HPhPs have become a metrological tool as much as a subject of fundamental optics. In hBN, buried defects act as internal HPhP reflectors. By analyzing fringe spacing, reflection coefficient, and dispersion near an internal air gap, one study reconstructed a concealed defect at depth (z=111\text{–}113) nm below the top surface with thickness (d=18\text{–}23) nm in a (157) nm slab [1808.03660]. This was achieved from measured polariton wavelengths such as (1215) nm at (1515\ \mathrm{cm}^{-1}) and (1388) nm at (1504\ \mathrm{cm}^{-1}), combined with electromagnetic simulations [1808.03660]. The same work analyzed reflection, transmission, and scattering versus defect size, showing stronger reflection for thicker defects and at higher frequencies because the defect becomes more comparable to the polariton wavelength [1808.03660].

Photodetection offers a direct device-level application. In graphene/hBN heterostructures above split metallic gates, the gate gap launches HPPs that propagate as confined rays and heat graphene at a (pn)-junction, producing a photo-thermoelectric photocurrent [1705.10318]. The external responsivity reached about (1\ \mathrm{V/W}) at room temperature and zero bias, while the internal responsivity was inferred to be about (150\ \mathrm{V/W}) with a noise-equivalent power of about (26\ \mathrm{pW}/\sqrt{\mathrm{Hz}}) [1705.10318]. The spectral peak can be tuned over about (60\ \mathrm{cm}^{-1}) by varying hBN thickness and split-gate gap width [1705.10318].

Far-field optics is becoming increasingly accessible. Patterned (\alpha)-MoO(_3) photonic hypercrystals convert free-space photons into in-plane HPhPs via reciprocal lattice vectors, producing polarization-selective reflectance resonances in both Reststrahlen bands [2309.17146]. Rotating the square hole lattice relative to the (\alpha)-MoO(_3) crystal axes tunes the far-field resonances because the rotated reciprocal vectors sample different directions of the anisotropic HPhP dispersion [2309.17146]. This provides a direct far-field route to probing in-plane hyperbolicity without s-SNOM.

Quantum-optical directions have also emerged. A cavity-QED treatment of hBN color centers proposed two HPP generation schemes: spontaneous emission into the phonon sideband and a stimulated Raman process [2602.05736]. In ultrathin slabs, spontaneous emission becomes effectively single-mode with enhanced decay into a selected HPP branch, while the Raman process provides tunable and narrowband excitation of ray-like HPPs that propagate over micrometer distances [2602.05736]. The work further outlines a two-emitter correlation measurement to test the single-polariton character of these emissions [2602.05736]. This suggests that HPPs may become not only classical nanophotonic carriers but also quantum channels coupling spatially separated emitters.

A recurring misconception is that HPhPs are purely a near-field curiosity. The literature now includes electrically detected HPP-guided photodetectors [1705.10318], far-field hypercrystal coupling [2309.17146], ultrafast thermal transport [2401.09594], and quantum-emitter-based source proposals [2602.05736]. Another misconception is that all HPhP functionality comes from nanofabricated resonators. Naturally grown (\alpha)-MoO(_3) nanobelts already behave as ultrabroadband HPhP nanocavities, demonstrating that crystallographic growth alone can define usable cavity geometries [2105.01148].

Taken together, these results place hyperbolic phonon polaritons at the center of a broad mid-IR and far-IR photonics program: deeply subwavelength waveguiding, cavity formation, local and nonlocal energy transfer, directional beam engineering, environmental sensing, defect tomography, and potentially quantum transduction. The unifying principle is the same in every case: indefinite anisotropic permittivity creates a hyperbolic phase space of high-(k) states, and that phase space can be controlled through thickness, substrate, gap, twist, edge geometry, or heterointerface design [2105.01148], [2604.08996], [2603.06184], [2601.16465].