The fast pace in progress of visible/near-infrared (IR) plasmonic optics ,  has not been mirrored in the IR part of the electromagnetic (EM) spectrum. However, a plethora of applications would require extraordinary EM control at longer wavelengths throughout the mid-infrared (MIR) spectrum , , between 3 and 30 μm free-space wavelengths, and into the far-IR and terahertz (THz) spectrum  of sub-mm free-space wavelengths. A few characteristic examples of such applications are bio-molecule fingerprinting , , , inspections of energy efficiency of buildings , passive radiative cooling , , as well as diagnostic tools in dentistry and cancer detection . Optical and optoelectronic components, such as sources, detectors, beam guides, beam splitters and beam-steering devices, that function in the mid- and far-IR , especially in integrable chipscale form , , are thus extremely important for studying and advancing the pertinent physics and developing new optical systems at these long wavelengths.
Here, we review the progress and most recent advances of the emerging and fast growing field of “phonon-polaritonics” ,  or “reststrahlen optics” , which has brought about a range of versatile possibilities for harnessing IR light with the so-called phonon-polariton (PhP) or reststrahlen-band materials , , , , , . These materials were initially overlooked for applications requiring manipulation and control of IR light because of their near-perfect reflector properties in bulk form. However, the wide range of works that we review, highlight and discuss here have unveiled their powerful potential for IR photonics. We analyze in the following how these unconventional capabilities for IR light emission, detection, and manipulation emanate from the different possible photonic properties in the vicinity of the reststrahlen band. We will also discuss how some of these extraordinary IR effects are unique to PhP platforms and are not possible with metals.
In particular, this review article is structured as follows. In Section 2, we discuss why traditional metallic-based platforms which have revolutionized light control at optical frequencies do not possess the full range of light manipulation capabilities in the IR. In Section 3, we describe the characteristics of PhP materials, explaining how they obtain their interesting photonic responses from a resonant interaction between their phonon vibrations and the impinging EM field. We present examples of materials in this class, ranging from ionic crystals, semiconductors, as well as two-dimensional (2D) materials, along with the associated frequency spectrum of operation in Section 4. In Section 5 we review the distinct photonic responses that PhPs make available that include exotic indefinite permittivity responses  initially thought of as possible only with metamaterials. In Section 6, we go through the different types of modes supported by PhP media which are the essential building blocks providing the foundation for the progress in the phonon-polaritonics field. In Section 7, we review the capabilities unleashed by popular photonic material structuring routes when applied to PhP materials. In Section 8, we discuss avenues to tune, switch or reconfigure the responses of PhP material. In Section 9, we give examples of how to bring together PhP platforms with other photonic systems for new functionalities in hybrid platforms. Finally, we present our conclusions in Section 10 with a brief outlook on the future of the field.
2 Metallic platforms for IR light: possibilities and limitations
In order to provide a context for what capabilities and directions phonon-polariton platforms bring about in IR photonics we briefly review in this section widely-researched metallic platforms for IR light and discuss their possibilities and limitations. A range of unique IR photonic effects have been reported leveraging metallic platforms with significant relevance to a wide range of IR photonic applications. Examples include IR metasurfaces , resonant or impedance matched platforms for super-absorption1 , , spoof surface plasmons (SPs) , as well as resonant  or non-resonant  IR transparency and extraordinary near-field enhancement in perforated arrays. Actually, most of these platforms operate based on resonances induced by the geometry rather than the metallic material properties of the structural elements , , , , . All share a common feature: there is little overlap between IR light and the metallic matter, and the three-dimensional (3D) localization of the field around the metal constituents is rather weak.
This rather weak overlap between light and metallic matter is inherent to metallic platforms at IR frequencies. It emanates from the relatively small skin depth, δ, when compared to the free-space wavelength λ0 (e.g. see Figure 1).2 We note that despite the ultra-small skin depth with respect to the free-space wavelength in the IR spectrum the dissipative (ohmic) losses inside the metal may not be necessarily negligible; specifically in certain metamaterial-based IR detector architectures metallic losses as high as 20% have been reported . However, such ultra-small relative skin depth most certainly prevents IR metal-optics systems  from reaching the full EM-field-sculpting capabilities of plasmonic optics . As a result, metals in the IR do not support the familiar highly-confined optical states at metallic planar interfaces, – known as SPs or surface-plasmon polaritons (SPPs) ,  –, or around metallic particles, – known as localized surface plasmons resonances (LSPRs) , particle plasmons , or Mie plasmons , .
Indeed, this is attested by the representative examples for SPPs and LSPRs for gold (Au) shown in Figure 2. Specifically, Figure 2A depicts the dispersion relation, i.e. the frequency ω, versus the real part of the lateral wave vector Re(kx), for the SPP modes at the vacuum-Au interface.3 We see that for IR light of wavelength larger than ~3 μm the SPP mode dispersion is very close to the lightline (vertical blue line) signifying the free-space EM wave dispersion. The SPP mode’s field profile, at these long wavelengths, has more in common with a grazing EM wave propagating along the x-direction rather than a surface wave localized around the interface (see panel SP1 of Figure 2B). The situation changes as the frequency increases towards the visible spectrum; the SPP mode moves away from the lightline while becoming more localized around the interface (e.g. see panel SP2 of Figure 2B).
To evaluate IR LSPR modes we calculate the extinction efficiency, Qext (see Figure 2C) accounting both for scattering and absorption for an EM wave launched onto a Au sphere.4 Small, nm-sized, spheres are resonant at the visible spectrum, with absorption dominating their extinction response. On the other hand larger, sub-micron/micron-sized particles show a flat response, in the IR spectrum; the extinction response is absent of clearly defined resonances, and both scattering and absorption make a comparable contribution. Furthermore, the field decays much faster outside the sphere for the nm-sized sphere (mode LSP15; see Figure 2D). In other words, near-resonance visible light can be strongly localized around a nm-sized metal particle (e.g. mode LSP1), while IR EM-fields are only weakly confined around larger micron-sized particles over a spectral region without clearly defined resonances (e.g. modes LSP2 and LSP3) (see also corresponding electric field intensity landscapes in Figure 2E).
The tight electric-field confinement around planar interfaces or particle geometries underpins most observed plasmonic phenomena in the visible spectrum. Moreover, certain types of photonic metamaterials whose design relies on a judicious volumetric mixing between constituents with opposite permittivities, such as epsilon-near-zero ,  or hyperbolic metamaterials ,  also require constituents with photonic responses similar to those yielding plasmonic phenomena. Clearly, then, the aforementioned phenomena and properties would not be transferable to the mid-/far-IR spectrum with traditional metallic structures. A very active research direction to bring the familiar effects of visible/near-IR plasmonics into the mid-/far-IR spectrum involves alternative plasmonic materials for IR light such as graphene , , conductive oxides , , or heavily-doped semiconductors , , . We direct the reader to , ,  reviewing this broad and highly-promising area.
Here, we review a different burgeoning research direction that involves tailoring IR photonic behaviors that do not emanate predominantly from the material free carriers. This research direction involves PhP materials which owe their photonic responses to the strong coupling of light and certain vibrational modes of the material lattice. As we discuss in the sections that follow, these photonic responses not only may mimic those of noble metals at visible frequencies and thus are ideal candidates for enabling “plasmonic effects” at the mid-/far-IR spectrum but also enable unconventional photonic responses and associated extraordinary modes.
3 Beyond metals: IR photonics with PhP materials
3.1 What are PhPs?
The native out-of-phase vibrational modes of the lattice (optical phonons) can strongly couple to impinging IR EM fields (photons) for certain types of solids. Such strong coupling can be triggered only in solids comprising oppositely charged atoms in their lattice basis so that the incident electric field can drive their out-of-phase vibrational modes; also known in the spectroscopy community as IR active phonons . As a result of such strong coupling, an EM wave going through such solid would have mixed wave dispersion characteristics, which are neither “photon-like” nor “phonon-like” . We depict the wave dispersion of these “mixed modes,” known as PhPs, in Figure 3 (red solid line). Polar covalent semiconductors can be seen as effectively comprising two separate oppositely charged ions  and therefore also support IR active phonons and thus PhP modes.
3.2 A microscopic model for the macroscopic PhP photonic responses
To understand the strong phonon-photon coupling phenomenon and its effect on the macroscopic optical response of ionic or ionic-like solids, we consider a basic, simple cubic diatomic lattice model. This comprises two oppositely charged ions with charges +Q, −Q and respective masses M+, M−, as depicted in Figure 4A. These charges correspond to effective charges in real materials as the concept of a perfect ionic medium is an idealized model . We consider an EM wave incident with electric field Einc along the x-direction, launched from atop the structure of Figure 4A along the +y-direction. In order to satisfy the EM boundary conditions at the interface the wave-vector of the coupled radiation field-phonon wave inside the ionic medium, q, must also lie in the y-direction . Also, from Maxwell’s equations it follows that the displacement vector D should be divergence-free inside the ionic material, implying that q·D=0 . Assuming linear and isotropic materials:
where P is the polarization, defined as the dipolar moment per unit volume, ε is the relative permittivity of the material and ε0 is the vacuum permittivity. It follows from all the above, that D, P, and E are along the x-direction inside the ionic material and normal to the wave-vector q of the PhP wave.
In order to understand the phonon-photon strong coupling we must understand how the macroscopic electric field, E, and polarization, P, correlate with the microscopic properties and motion of the atoms, as well as with the microscopic local electric field Eloc. In lattices of tetrahedral symmetry, the local microscopic field at the center of both types of ions ,  is given by:
This local electric field, Eloc, is the superposition of the impinging EM wave, and the long-range electrostatic field of the rest of the ions in the lattice. The macroscopic polarization in the ionic crystal, P, will be a superposition of the polarization resulting from the dipole moments of each kind of the two polarizable ions and the polarization from the dipole moments emanating from the oppositely-charged ions as they move. In particular, for the former contribution, we assume immovable but deformable ions under the influence of the local electric field Eloc. The local electric field shifts the electron cloud around the ionic cores in each ion type (see Figure 4B). The contributing dipole moments of each ion type (p− and p+) yield a total macroscopic polarization of:
where N represents the number of positive-negative ion pairs per unit volume, and α+, α− the respective polarizabilities of the positively-charged and negatively-charged ions. We have set α=α++α− to represent the total polarizability of the oppositely-charged ion pair.
Conversely, for the latter contribution, i.e. for the polarization from the dipole moments emanating from the oppositely-charged ions as they move, we consider the oppositely-charged ions as rigid ions moving under the influence of the local electric field Eloc and the restoring forces from the inter-atomic potential. By rigid ions we mean that the ions move with the electron cloud tight around them , i.e. neither of the ions is polarizable (see Figure 4C). Then a net dipole moment in the elementary cell can be obtained only when there is a net displacement w of the positive ions with respect to the negative ions, yielding a total macroscopic polarization of:
Accounting for both types of contributions to the microscopic dipole moments, the total polarization, P, would be:
As the macroscopic polarization and field are parallel to each other [see Eq. (1)], i.e. E | |P, it follows from Eq. (2) that the microscopic local field, Eloc, is parallel to both of them. Then Eq. (5) implies that the relative displacement between the positive and negative ions, w, must also remain parallel to all fields. So, we have that Einc | |Eloc | |E | |P | |w, with all perpendicular to the wave vector q. In other words, w·q=0, meaning that the impinging EM can drive only the transverse phonon vibrational waves , .
Eqs. (2) and (5) can be also used to express both the local electric field, Eloc, and the polarization, P in term of the macroscopic field E. We get:
So, all displacements and fields are along the x-direction, while the PhP mode propagates along the y-direction. If the incoming EM field is a harmonic field of frequency ω, so will be the macroscopic field inside the ionic medium, which still must satisfy Maxwell’s equations. As Eqs. (2), (6) and (7) must be obeyed at any depth of the lattice, the macroscopic field, E, the local field, Eloc, the polarization, P, and the relative displacement, w, all would follow the same waveform , i.e. we have:
with E0, Eloc,0, P0 and w0 representing their respective amplitudes.
To understand further the motion of the ions in the lattice (in the rigid-ion model of Figure 4C), we take Eloc at a certain depth y to be constant, i.e. we ignore any spatial variations at the inter-ion scale , . Thus, it suffices to focus our attention on two particular ions at certain depth y, one positive and one negative [e.g. see ions designated with (1) and (2) in Figure 4C]. The short-range forces between adjacent ions can be approximated as restoring forces , similar to the elastic force between two masses connected with a spring. A damping force proportional to the velocity of the ion is also taken, to account for dissipation in realistic lattices with the assumption that the more massive ion experiences a proportionally stronger damping force. We take u+, u− to represent the respective displacements along the x-direction of the positive and negative ions in the chain from their equilibrium position. We also take k to represent the elastic constant of the restoring force on the ions, and γ to be a constant related to the strength of the damping force. Then the the equation of motion becomes , :
for the positive ions, and
for the negative ions.
In order to understand better how the incident field drives the motion of the two kinds of ions in the lattice we temporarily ignore the damping term in Eqs. (9) and (10). Then, by adding these and integrating over time (assuming at t=0 the ions are at their equilibrium position) we get that:
Eq. (11) signifies that the motion of the positive and negative ions will be out-of-phase with u+u−<0. Also, if w0 is the maximum value of the relative displacement, w, along the x-direction, i.e. we would have that:
The PhP wave is a slow varying wave in the y-direction. Meaning, the wave would need to cross several atoms typically of the order of a 1000 or more for a full cycle. In Figure 5, we depict a snapshot of the transverse PhP wave.
We follow through to determine the amplitude of the relative ion displacement, w0, needed to determine the macroscopic polarization and thus the effective photonic response of the ionic medium. It can be shown, by subtracting Eq. (10) from Eq. (9), after multiplying them respectively, with M+ and M−, and using Eq. (7) to express the local field in terms of the macroscopic field, that:
with M being the reduced mass of the the ion-pair of Figure 4C given by:
and using Eq. (8) into Eq. (14) we obtain:
We first consider a limiting case: when the driving field is a high frequency field. For such a case, the atomic polarizability from the electron displacement dominates the permittivity response; the fields change too fast for the ions to respond (the atoms having much larger inertia in comparison to electrons). This means that the displacement, w, in Eq. (6) can be taken to be zero. We can see that also from Eq. (17) as w0→0 as ω→∞. Then from Eqs. (1) and (6) we get the familiar Clausius-Mossotti relation  for the high-frequency permittivity, ε∞:
Eq. (18) allows simplifying (16) to:
Now we consider another limiting case: the driving field being a static field, i.e. ω=0. Then from Eqs. (17) along with Eqs. (1) and (6) we get for the static permittivity εs:
Conversely, we get for the permittivity response [with the successful application of Eqs. (1), (6), (17) and (20)]:6
To understand the physical meaning of the frequency ωR in Eq. (21) we assume for a moment an undamped system, i.e. a system with γ=0; then as ω→ωR, Eq. (17) implies that the amplitude of oscillations becomes infinite. Essentially, ωR, represents the frequency of the transverse optical (TO) phonon vibrations that are driven by a static field . We will designate such frequency as ωTO and simply refer to it as TO phonon frequency, i.e. ωR≡ωTO. The TO phonon frequency is somewhat different from the free vibration optical phonon mode, Ωfree, of a diatomic chain which for q=0  is:
Then, with the use of Eqs. (16), (18) and (22), it can be shown that:
which means the TO phonon frequency of modes driven by a static electric field is lower than the free vibrational optical phonon modes of a diatomic chain (both are designated with horizontal long-dashed and short-dashed lines, respectively in Figure 3).
Continuing, we can observe that at the permittivity, ε, becomes zero (for the undamped system). Previously, we stated that only TO phonons couple to the radiation field. Indeed, the divergence-free nature of the displacement field D, mandates εq·E to be zero which requires the wave vector q to be perpendicular to the driving field E, i.e. requires transverse modes. However, when the permittivity, ε, is exactly zero, the divergence-free condition for the displacement field D may be satisfied without requiring transverse modes, meaning longitudinal modes may be possible in this case. This is the reason why the frequency at which the permittivity, ε becomes zero for the undamped system, is called the longitudinal optical (LO) phonon frequency , ωLO. Hence we obtain the famous Lyddane-Sachs-Teller (LST) relation , :
With the use of the LST relation, we are able to obtain from Eq. (21) the familiar expression for the dielectric function of PhP media:
Eq. (25) encompasses the basic photonic responses that result from strong phonon-photon coupling as derived from a classical microscopic picture for the motion of the ions and the associated electron cloud. The particular parameters of a certain material are extracted from experimental reflectivity and/or ellipsometry spectra (see Section 4). Very recently Rivera et al.  showed that a quantum mechanical model leads to an effective permittivity response that is essentially equivalent to that of Eqs. (21) and (25).7 The quantum model , which encompasses treatment of more complex unit cells with multiple-atom basis and tensorial Born effective charges, provides a pathway to first-principle predictions of PhP permittivities, which is particularly important for emergent materials where experimental data may not be available.
3.3 The PhP gap and the reststrahlen band
Figure 3 depicts the coupled phonon-radiation field wave dispersion, ω versus q, i.e. the PhP dispersion (red solid lines). This is determined by where ε(ω) is the permittivity response of Eq. (25). The two mixed phonon-photon mode branches are strongly modified from both the phonon dispersion (long-dashed horizontal line) and the high-frequency photon dispersion (slanted solid line, with ). The phonon-mode dispersion appears flat as it represents phonons with a near-zero phonon wave vector, q. As EM waves are extremely slow varying in the atomic scale it is only the near-zero-q phonons which may couple to EM waves. The frequency region where no mixed phonon-photon modes (yellow-shaded region) is known as the PhP gap signifying forbidden EM propagation; almost all impinging light would get reflected from the material [see Figure 3 (right panel) showing the reflectance from a respective slab with thickness about a half the impinging light’s wavelength].
Alternatively, the PhP gap is called the reststrahlen band; a term coming from German for the spectrum of “residual radiation”. Historically, this term emerged  from experiments where a pulsed signal is reflected back and forth between slabs made of a PhP material. While the radiation comprising most of the pulse’s frequency spectrum transmits through the slabs and escapes the set-up, the part of the pulse with frequencies within the PhP gap nearly-perfectly reflects on the surfaces of these slabs, thereby remaining trapped in the slit area between them. These “residual” radiation/rays  exit eventually through the slit as a narrow-band beam after multiple bounces. Set-ups like these had been used in the past to create narrowband/quasi-monochromatic sources from broadband thermal sources .
The permittivity is negative, i.e. Re(ε)<0, within the PhP gap/reststrahlen band. In the sections that follow we review interesting photonic phenomena and properties within but also around the vicinity of the reststrahlen band; we will be referring to this broader spectral region as the PhP regime.
3.4 Extensions to the basic PhP permittivity model for a broader class of material
The basic permittivity model of Eq. (25) covers successfully the photonic response of some materials but it may be too crude for others where the assumptions made in Section 3.2 fail by a large degree. For example, in the microscopic model the damping was simply associated with a single force proportional to the velocity of the ions. However, there are different distinct mechanisms causing damping of the PhP mode. Mechanisms that can cause PhP damping are scattering at crystal defects , as well as decay of the phonon-part of the PhP into acoustic phonons , , , ; the latter is caused by the presence of anharmonic terms in the interatomic forces . Damping can also occur because of other coupling to low frequency phonon excitations aided by strain or defects , , .
Schwarz and Maier  have asserted that damping caused by the presence of defects and the anharmonicity in the ionic potential can be lumped together in the decay constant γ of Eq. (25). However, according to  coupling to low frequency phonon excitations adds up a number of Lorentzian oscillatory terms in the denominator of Eq. (25) with associated poles at the zeros of where ωn and γn are the resonant frequencies and damping constants of the low-frequency phonon excitations. In materials where multiple of these damping mechanisms are dominant, as for example, in lithium niobate , other authors , ,  have considered a frequency-dependent decay constant in the model of Eq. (25), γ(ω); a causal model  however requires γ(ω) to be itself complex valued . The treatment of  offers more insight to the actual physical mechanisms of damping, beyond considering ad-hoc parameters that improve the agreement between the model-predicted PhPs and experimental results .
Furthermore, depending on the lattice symmetry multiple optical phonon branches may exist also. If these are well-separated in frequency, Eq. (25) serves as a good model in the neighborhood of the TO phonon frequency of interest. By well-separated we mean that there should not be a significant spectral overlap between adjacent TO phonon resonances, when taking into account their respective broadening as determined by the decay constant, γl, associated with the respective phonon resonance, ωTO,l. This however, is not true for many materials for which in order to construct the photonic response N independent harmonic oscillators should be considered. Each oscillator corresponds to one of the relevant TO phonon resonances ωTO,l with damping γl with a corresponding LO phonon frequency ωLO,l, that is independent of the other oscillators . So, the full permittivity response would be in such case:
Not all the oscillators contribute equally; the strength of each oscillator, Sl, is determined by :
implying that a certain oscillator makes a stronger contribution the larger the relative difference between the respective LO and TO phonon frequencies.
Note here that, the independent oscillator model of Eq. (26) would have led to the same decay constant, γl, appearing in both the numerator and denominator of Eq. (28) , . However, Gervais and Piriou  have asserted that a model that permits a different broadening of the LO and TO phonon resonance, – or in other words a different decay constant in the numerator, γLO,l and in the denominator, γTO,l respectively –, better captures the permittivity response for materials where the anharmonic decay channel of the PhP is strong , . There are some constraints in the permissible values for the decay constants, mandated by the passivity condition for the permittivity model which requires the imaginary part of the permittivity function to be positive or zero, i.e. Im[ε(ω)]≥0 . In particular , :
4 Widely-used and emergent PhP materials
4.1 Ionic crystals (insulators)
Ionic crystals are the quintessential example of media with a PhP response. Table 1 outlines the permittivity model parameters of Eq. (26) for some characteristic material of this class. The TO and LO phonon frequencies are accompanied by their corresponding free-space wavelengths, λTO and λLO. The representation of the reststrahlen band in terms of free-space wavelength is helpful for gauging the relevant feature-size ranges for platforms functional in this spectrum. We will therefore include such representation for all materials in the subsequent tables, 1 through 5 that follow.8 While the subject of this review is how to utilize the reststrahlen band for photonic applications, we remark in passing that for some applications the reststrahlen band may need to be avoided. For example, CaF2 and BaF2 are popular as IR transparent materials, in which cases they are utilized for spectral regions away from their reststrahlen band.
4.2 Semiconductors: semi-insulating versus doped materials
4.2.1 Isotropic semiconductors
Heteropolar and polar covalent semiconductors comprise effective cation-anion pairs; so they may also exhibit IR active phonons, and hence PhPs. However, unlike the case of the insulating ionic crystals of Table 1, in general their permittivity response is not fully captured by the models of Section 3; these would be only appropriate for semi-insulating or undoped semiconductors. The full permittivity response of semiconductors must include the contribution of the free carriers typically showing a Drude response  to the macroscopic electric field. It can be reasonably assumed that the free carriers do not contribute to the local field experienced by the ions.9 So the semiconductor’s full permittivity response, εsem, can be considered as the cumulative result of the electric-field induced motion of the ions and their bound electrons as if free carriers were not there and then the electric-field induced motion of the free carriers as if the ions and their bound electrons were immovable. In other words:
where the ionic contribution εion has the same expression as one of the models of Section 3 as appropriate per case; εDrude represents the Drude permittivity response of the free carriers, i.e.:
For generality the contribution of both types of carriers, electrons and holes, is taken  in Eq. (32) with respective carrier densities Ne and Nh and effective masses and e is the electron charge, ε0 the vacuum permittivity and γp the decay-rate/damping constant of the free-carriers motion.11 The Drude response should be suitably adapted to account for cases where contributions from multiple valleys or sub-bands are relevant (e.g. see ).
In Table 2 we show some examples of widely used semiconductor material and the parameters relevant to their PhP permittivity response, εion. The parameters refer to the model of Eq. (26) (with N=1) except for the case of AlN that requires the model of Eq. (28) (with N=1). In Figure 6A, we plot the real part of the ionic permittivity response but without the bound electron contribution, i.e. Re(εion–ε∞), for the materials and model parameters of Table 2. We observe the PhP responses of the different semiconductors cover a wide IR spectral range from about 10 μm to 60 μm. We also show in Figure 6B the related skin depth, δ, scaled with the free-space wavelength λ0 along with the skin-depth value for Au (solid orange line), – multiplied with a factor of a hundred for an easier comparison. We observe high PhP related skin-depth values for all semiconductors, much higher than the skin-depth of Au. Thus, semiconductors in the reststrahlen band are ideal material for strong IR light-matter interactions.
The examples of Table 2 and associated models capture the photonic responses of most semiconductors. Semiconductor alloys tend to show the vibrational resonances of their constituents typically leading to photonic responses that can be described with the multi-pole models of Eqs. (26) and (28) . For small gap or zero-gap semiconductors, as for example, HgTe, the full picture is more involved ,  as the intra- and inter-band transitions spectrally overlap with the PhP regime and contribute to the permittivity response; these fall out of the scope of phenomena and applications in this review. Furthermore, some semiconductors show directional dependent responses that we discuss later.
4.2.2 Anisotropic semiconductors
The permittivity models of Eqs. (26) or (28) may unwittingly imply that the ionic photonic responses are independent of the propagation direction. This is not true; in materials with crystallographic anisotropy the TO phonon will depend on the direction of the electric field inside the material. These anisotropic materials still follow the aforementioned models as appropriate but with different parameters for electric field oriented along different crystallographic directions , . Typically, they are uniaxial media  with their c-axis coinciding with the optical axis.
SiC(6H) shows a small anisotropy per  but has been predominantly treated as an isotropic PhP material , , , , ,  with the parameters of  which are the ones shown in Table 2. On the other hand, SiC (4H) is typically considered as an anisotropic PhP material ,  with the anisotropy coming in large part from the anisotropy of the high-frequency permittivity ε∞ . The nitride semiconductors of Table 2, AlN and GaN, are also anisotropic with the most significant anisotropy contributed by the TO phonon frequency. Table 3 complements Table 2 in capturing their full directional PhP response.
4.3 Other anisotropic PhP material
Aside from the semiconductors of Section 4.2.2, other widely used materials in photonics with an anisotropic PhP response are quartz , sapphire  and vanadium oxide (VO2) , . Their permittivity model parameters are given in Table 4, where we observe these materials possess multiple PhP resonances. The anisotropy in TO phonon frequency underpins the existence of extreme, hyperbolic, photonic responses as we discuss in more detail in Section 5.2 and drives a strong interest towards materials with this property. Furthermore, VO2 has additional interest as a phase-transition material  offering tunability by switching on or off PhP responses (see Section 8.5).
One particular material, hexagonal boron nitride (hBN), has received a recent heightened attention for IR applications because of its anisotropic PhP responses , , , , , , , , , , , , , , , , , . While in principle the exotic hyperbolic modes exist in all uniaxial PhP material, what sets hBN apart is the degree of TO-phonon-frequency anisotropy leading to a much broader bandwidth where these modes are available. hBN is a van der Waals (vdW) crystal (see Figure 7A). The much looser connection between the hexagonal-lattice planar sheets with weaker vdW forces significantly red-shifts the out-of-plane IR-active phonon frequency [see Eq. (22)] and accordingly the TO phonon frequency [see Eq. (23)]. We show hBNs in-plane (⊥ to the c-axis) and out-of-plane (|| to the c-axis) IR active phonons in Figure 7B and C, respectively, designated with Mulliken’s notation  that is common in phonon/Raman spectroscopy communities. In Table 5 we list hBN’s PhP permittivity model parameters. Indeed, we observe a tremendous anisotropy between the in-plane (⊥ to the c-axis) and out-of plane (|| to the c-axis) PhP responses of hBN. To quantify such anisotropy in the observed PhP responses we introduce the PhP anisotropy factor, fanis, defined as:
where ωTO,⊥, ωTO,|| are the respective TO-phonon frequencies for electric field perpendicular and parallel to the c-axis. hBN has a PhP anisotropy factor, fanis, of more than 50% compared to only ~10% and ~5% fanis for AlN and GaN, respectively. Large fanis values indicate a wider spectral band where indefinite uniaxial photonic responses  and associated exotic hyperbolic modes can be available as we discuss in Section 5.2.1.
While the model parameters of the Tables 1–4 are typically obtained by fitting to experimental reflectivity spectra, for hBN this process would be highly non-trivial. Available experimental data highly depend on the polycrystalline nature of the samples, leading to substantial overestimation of the decay constants, as well as on domain off-plane orientation defects , leading to mixing between the on- and off-plane responses . As a result the older experimentally extracted model parameters  may not be appropriate for single crystalline hBN systems of current interest . It appears that similar but different permittivity parameters have been used in different works with no clear consensus. We choose to indicate in Table 5 the model parameters that have been adopted by many current works , , , , following the ones of Dai et al. ; i.e. ωTO, ωLO, from , ε∞ values from the ab-initio calculations of  and γ, from the observed decay of exotic guided modes in hBN by Dai et al. .
Another emergent vdW material for PhP-based photonics is α-MoO3 (see Figure 8A). Due to its crystal structure this material has a very rich vibrational response. We focus our discussion here on the main three IR active modes relevant to photonic applications in the free-space wavelength range around 10 μm–20 μm . The vdW planes connect same-species atoms, O-atoms in this case; so the out-of-plane PhP response comes only from the asymmetric stretching vibration of the Mo-O(1) bond (B2u mode) (see Figure 8B). This is not too different from the B3u asymmetric-stretching mode of the Mo-O(3) bond; so both are spectrally very close. However, the O(2) forms two bonds of significantly unequal length with the adjacent Mo atoms , which significantly red-shifts the associated IR active mode (mode B1u) with respect to the B2u and B3u modes. As a result α-MoO3 has a biaxial PhP response (see also Section 5.2.2) with a large in-plane anisotropy (plane ⊥ to the stacking direction ). To be specific, an analogous PhP anisotropy factor, fanis−IP, can be defined for the in-plane (IP) PhP modes as:
where x and z define the two in-plane symmetry directions of α-MoO3,  and . The in-plane PhP anisotropy factor for α-MoO3, fanis−IP, is about 40% indicating the possibility for broad spectral availability of indefinite biaxial photonic responses  that we discuss in Section 5.2.2.
Thus, hBN and α-MoO3 are PhP material of high current interest because they exhibit an extreme PhP anisotropy, out-of-plane and in-plane, respectively; this anisotropy leads to unusual photonic responses that we discuss in more detail in Section 5.2. But this tells only part of the story. The weaker cross-plane vdW forces also allow for straightforward fabrication of few layer hBN  or α-MoO3 thereby unlocking extreme possibilities for molding EM waves with PhPs in atomic scale materials. In the section that follows, we go beyond the examples of hBN and α-MoO3 and discuss the importance and implications of 2D and quasi-2D PhP materials for phonon-polaritonics.
4.4 PhPs in the “flatland”: 2D and quasi-2D PhP materials
2D materials have broken boundaries in material science and are a subject of intense research efforts both for their appealing electronic  and photonic ,  properties. Going a step further, Geim and Grigorieva  envisioned that 2D material could be used as building-block “legos” that can be joined together plane-by-plane to make materials with new properties. These super-2D materials are structurally stable as the material planes are held together with vdW forces in the absence of dangling bonds; hence their name vdW heterostructures , . When it comes to PhPs these vdW heterostructures can provide an avenue to make atomic-level designer PhP material with tailored functional spectral ranges of their IR-active phonons (see schematic of Figure 9). To the best of our knowledge these atomic-level designer PhP materials, comprising dissimilar adjacent vdW planes, are uncharted territory at present. Their possibility has only been envisioned by Caldwell and Novoselov in  and briefly discussed in , where they were referred to as “crystalline hybrids”.
One could argue that for these atomic-level designer materials the effective in-plane PhP responses might bear some similarities to the case of semiconductor alloys, especially as the intra-plane bonds are much stronger than the inter-plane forces. Nevertheless, interlayer charge transfer effects as well as the charged-phonon effect , , ,  that depend on the particular materials involved may cause completely unexpected responses. The resulting effective permittivity could be very different from the one hypothesized from simpler classical arguments as above. In addition, new IR-active phonons for an electric field perpendicular to the atomic planes could emerge, as a result of possible new IR-activated lattice vibrations in the cross-plane direction. Most certainly, these vdW hetero-composites open exciting bottom-up design routes to new materials with different PhP spectral ranges. vdW hetero-composites could form a very promising research direction for on-demand IR light control with atomically-thin materials.
Although vdW hetero-composites have not been studied in the context of creating new effective photonic responses, vdW crystals supporting IR-active phonons have been interfaced with other 2D material with the purpose of achieving strong coupling between their respective native EM modes (e.g. see ). The underlying physics principles of these systems are to tune/modulate the EM wave propagation properties by a strong interaction between the EM wave modes supported in each sub-system. As these do not relate to changing inherent material properties we will review these phenomena in detail in Section 9 where we discuss hybrid platforms. The vast possibilities of these hybrid platforms can have a strong impact on future photonic applications. This has stimulated a strong interest in the community for realizing 2D/quasi-2D PhP materials.
Many of the research efforts have mainly focused on hBN. Monolayer ,  and few layer hBN , ,  have been reported and studied for its PhP-related hyperbolic responses by various groups [see also Sections 6.1.2 and 6.1.4]. Detailed studies of atomic-scale (few-layer) hBN ,  have shown that exotic hyperbolic-type of guided modes survive even for atomic-scale thickness [albeit with a modified guided wave dispersion, ω(k)]. This fact provides evidence that the strong TO phonon anisotropy, in the respective planar and c-axis directions, is maintained for atomically-thin hBN . Indeed,  reports an excellent agreement between theoretical predictions and experimental results on the PhP guided-mode even for sub-5 nm thicknesses. Furthermore, Raman spectra for the similar in-plane Raman-active phonon mode in monolayer and few-layer hBN have confirmed a small shift only for the optical phonon frequency . These results in combination imply that the TO phonon frequency for the in-plane IR-active mode (Figure 7B) is fairly robust to finite-size effects.
This is likely not the case for the TO phonon mode across vdW planes of hBN, where the atoms oscillate along the c-axis as depicted in Figure 7C. In the extreme case of monolayer hBN, no out-of-plane vibration modes would be possible. Thus, the TO phonon mode along the c-axis completely disappears  in hBN monolayers along with the associated extraordinary hyperbolic photonic responses (see also Section 5.2.1). The significant interest in mid-IR photonics with atomically-thin hBN calls for further experimental evaluation of finite-size effects when thicknesses approach the atomic scale. Characterization of the in-plane and out-of-plane TO-phonon modes and their associated photonic responses for atomically-thin hBN will be extremely important in elucidating such effects.
As an insulator , hBN may have a limited potential for achieving PhP-based structures or devices with actively tunable responses. Semiconductor 2D/quasi-2D PhP materials would be highly desirable for tunable IR devices. Atomically-thin semiconductors that have been most frequently discussed as 2D PhP material, at least theoretically , , , are the transition metal dichalcogenides (TMDs). TMDs are also vdW crystals allowing their realization in a monolayer or few-layer form via mechanical exfoliation . Semiconducting TMD monolayers are typically in the 1H phase (hexagonal/trigonal prismatic arrangements) ,  as depicted in Figure 10A, which is the standard thermodynamically stable phase for sulfides and selenides of group VI TMDs . Although tellurides of group VI TMDs are typically semi-metallic following octahedral (trigonal anti-prismatic) arrangements , a recent work has shown that it is possible to mechanically exfoliate MoTe2 monolayers onto Si/SiO2 substrates in the semiconducting hexagonal phase .
What is remarkable about these 2D semiconducting TMDs is that even in the monolayer form they may support extraordinary hyperbolic modes, in stark contrast to monolayer hBN which loses the hyperbolic photonic response in its monolayer form. The monolayer-TMD hyperbolic photonic response is a direct consequence of the trigonal prismatic atomic arrangement, allowing out-of-phase vibrations between the transition metal and chalcogenide atoms not only along the atomic planes (see Figure 10B) but also in the direction perpendicular to the atomic planes, i.e. along the c-axis (see Figure 10C). The bandwidth of such hyperbolic photonic response depends on the spectral separation of the in- and out-of plane IR-active phonons. To the best of our knowledge, experimental characterization for the permittivity responses of monolayer TMDs in the mid/far-IR is not yet available. However, first-principle calculations based on density functional theory (DFT) can predict the resonant frequencies/wavelengths for the IR-active phonon modes of Figure 10B and C , , , which are related to the TO phonon frequency, ωTO [see Eq. (23)].
We list in Table 6 the IR-active phonon mode wavelengths for certain semiconducting monolayer TMDs (1H-phase), calculated from DFT at the Γ point , , , i.e. at zero phonon wave vector, q~0 (see discussion in Section 3). We observe from Table 6 that the typical anisotropy between the in-plane and out-of-plane IR-active phonon frequencies is about 20%. Although not exactly equal, the anisotropy between the in-plane and out-of-plane IR-active phonon frequencies relates to the PhP anisotropy factor, fanis, of Eq. (33) (see also footnote c in Table 6). Thus, a more moderate fanis of about 20% is expected for TMDs, much smaller than that of hBN, which has an fanis of about 50%, but still much larger than that of the traditional semiconductors of Table 3, which have an fanis of about 10% or less. Such significant PhP anisotropy factor makes hexagonal 2D TMDs excellent media for hyperbolic photonic responses with the added advantage of gate tunability, which is not available in hBN. The PhP modes listed in Table 6 would be of course lossy, just like in hBN and other semiconductors, because of the relaxation to the low-lying phonon bands . The actual PhP resonances in experimental samples may show a slight spectral shift depending on the isotopic purity of the TMD samples which affects the mass of the ionic oscillators .
In reviewing 2D materials as PhP media, special attention should be given to graphene, the first ever 2D material, discovered in 2004 . Since graphene’s discovery there have been numerous reports of its unique electronic properties and vastly tunable photonic response. However, graphene has not received much attention for its potential as a PhP material. Although graphene is comprised of only carbon atoms, it can actually support optical phonon modes . This is because graphene has two atoms within the hexagonal Bravais-lattice unit cell which can vibrate out-of-phase. However, because of the lack of charge asymmetry between the two atoms of the unit-cell basis, it is not possible for EM waves to trigger their out-of-phase vibrational modes. In other words, monolayer graphene supports optical phonon modes, which however are not IR active; hence PhPs cannot be supported on monolayer graphene.
The situation is a bit different in Bernal-stacked bilayer graphene (AB stacking) (see Figure 11). The structural asymmetry of the AB stacking, together with electron hopping between the two layers, breaks the equivalence between the A and B sub-lattice carbon atoms, leading to a dipole moment in each Bravais lattice unit cell. This dipole moment would enable coupling to incident light, i.e. it would make the optical phonon-mode IR active. Indeed, Kuzmenko et al.  observed a PhP response for Bernal-stacked bilayer graphene. However, the experimentally observed PhP showed a much higher oscillator strength, which in turn implied a three-orders-of-magnitude higher dipole moment from what would be expected from classical arguments .12 Also, the observed optical response, as manifest in the measured conductivity, exhibits a spectral profile that is characteristic of Fano resonances .
This IR active response can be strongly tuned by gating the Bernal-stacked bilayer graphene , . These curious results of  were interpreted within the quantum mechanical framework of charged-phonon theory adapted for the graphene-bilayer system , . In particular, based on this theory the “bare” optical phonon mode via interactions with the allowed electron-hole excitations acquires a large dipole moment. Hence, this “dressed” phonon has a much stronger effective dipole moment , , ,  than the one predicted by the asymmetry alone, thus leading to a prominent IR active phonon mode observed in related experiments , , , . In other words, this “dressed” phonon effect gives rise to a PhP with much higher oscillator strength. At the same time, the phonon part of this discrete PhP mode interacts with the electronic continuum via electron-phonon scattering, which leads to a Fano profile for the conductivity response , , , , .
Similar effects can occur also in trilayer graphene  with the ABC-stacking demonstrating a much stronger tunability with the gate voltage. Other explored 2D PhP materials include AB-stacked black-phosphorus-carbide (b-PC) , with IR-active phonons at frequencies similar to the IR-active spectrum graphene. 2D b-PC shows in-plane optical anisotropy and a Fano spectral response; the latter does not appear to be fully understood at present and could be related to similar effects that were observed in bi-/tri-layer graphene.
We discussed previously α-MoO3, a vdW transition metal oxide (TMO) , that also exhibits in-plane anisotropic PhPs , ; α-MoO3 with atomic-scale thickness has been recently realized and studied for its PhP responses . This emergent PhP material is expected to attract increasing interest as an atomically-thin bi-axial material  with extreme in-plane anisotropy and IR activity covering a spectral range between ~10 μm and 20 μm . Finally, phosphorene, i.e. a monolayer of black phosphorus has been shown to support IR-active phonons at wavelengths around 20 μm .
Summarizing, IR-active phonons in 2D materials, and associated PhP photonic responses, cover a large part of the IR spectrum between free-space wavelengths of 6 μm and 45 μm. In Figure 12 we outline such spectral coverage in IR activity with representative 2D materials discussed in this section. The 2D material are arranged from right to left in order of decreasing wavelength of their IR-active optical phonons.
5 Distinct photonic responses in the PhP regime
5.1 Photonic responses in isotropic PhP material
To visualize the photonic responses in the PhP regime we consider a general fictitious material, with permittivity following the model of Eq. (26) that broadly covers many PhP materials. We consider realistic model parameters within the range of the parameters for the semiconductors of Table 2. In particular, we take ωLO=1.3ωTO, ε∞=6, and γ=0.008ωTO (N=1). We plot in Figure 13 the real and imaginary part of the permittivity, versus the frequency, ω, that has been scaled with the TO frequency, ωTO. We can identify three distinct types of photonic responses that underpin distinct design principles for platforms relevant to IR applications.
5.1.1 IR plasmonic-like responses
Within the reststrahlen band, between the TO phonon frequency, ωTO, and LO phonon frequency, ωLO, the photonic response is metallic with Re(ε)<0 (blue-shaded region in Figure 13). The pertinent permittivity range for obtaining powerful light-molding capabilities like that of noble metals in the visible , that we call here plasmonic phenomena, does not have hard limits and depends on the particular geometries and sizes as well as the permittivity of the host media. Surely, large negative permittivity values correspond to smaller relative skin-depths and so do not support a strong light-matter interaction. Generally, plasmonic effects can be expected in the spectral range where the permittivity takes values roughly between ε=–15 and ε=0; note that permittivity values of ε=–1 and ε=–2 represent the quasi-static conditions for surface SPP and particle-plasmon respectively assuming the host medium is vacuum . Now, permittivity values between ε=–15 and ε=0 are available for a large portion of the reststrahlen band wherein IR plasmonic-type effects can be explored. We discuss the basic building blocks of such IR plasmonic-type effects in Section 6. These plasmonic-like responses are also relevant to constructing IR hyperbolic metamaterials that we discuss in Section 7.1.1. We designate the spectral region for IR plasmonic-like responses in Figure 13 with the blue rectangle.
5.1.2 IR high-refractive-index (HRI) responses
We also observe with the example medium of Figure 13 that the permittivity, and thus the refractive index, can reach very high values for frequencies close to but smaller than ωTO. We designate the region of a high-refractive index response with the yellow shaded area. However, there are some differences with pure dielectrics as a significant material loss is present (non-zero imaginary part in the permittivity). We discuss to what extend can this photonic response be utilized for photonic paradigms inspired from all-dielectric photonics in Sections 6.2.1 and 7.2.1. This response is also relevant to constructing IR magnetic metamaterials discussed in Section 7.1.3.
5.1.3 IR ENZ responses
ENZ photonic environments have been associated with enhancing non-linearity ,  and absorption ,  as well as modulating antenna resonances and their collective effects , . Naturally, PhP material have an epsilon-near-zero (ENZ) response near the LO phonon frequency and can thus bring such ENZ phenomena in the IR spectrum. Recent examples include  and  reporting on tailoring of antenna resonances and emission with PhP-enabled ENZ substrates. Also  reports a strong angular and spectrally-selective emission mode of an ultrathin PhP film at the ENZ point associated with the Berreman mode .
We designate in Figure 13 the ENZ point for the example medium we have considered. Generally, the quality of the ENZ response is determined by the accompanying imaginary part where the real-part of the permittivity is near-zero; lower imaginary part of the permittivity, at the ENZ point, generally corresponds to systems with lower losses and responses with near-zero refractive index. In particular, at the ENZ point with Re(ε)~0 the full complex refractive index is:
An ENZ response at the LO phonon frequency, ωLO is attained assuming:
This condition is reasonably satisfied by most single-oscillator PhP materials; in multiple-oscillator materials the off-resonance contributions by other oscillators may overshadow the ENZ response of the corresponding LO phonon frequency. For single-oscillator material the permittivity becomes, at the LO phonon frequency:
Eq. (37) implies that the quality of the ENZ response is better for single-oscillator material with smaller x values as given by Eq. (36).
5.2 Photonic responses in anisotropic PhP material
As we discussed in Section 4 some material, e.g. AlN, GaN, hBN, sapphire, insulating-phase VO2, semiconducting TMDs, α-MoO3 etc., have anisotropic PhP responses, with hBN dominating recent attention , , , , , , , , , , , , , , , , , , . hBN owes such increased attention to its being a vdW crystal that not only leads to an extreme PhP anisotropy but also brings such anisotropy in atomically-thin structures. α-MoO3 is another emergent PhP material attracting significant recent interest , ,  due to its inherent in-plane PhP anisotropy that enables a biaxial response. Also a vdW crystal, α-MoO3 can be made into 2D/quasi-2D structures. Next, we discuss how these anisotropic PhP responses, of uniaxial or biaxial type, translate to distinct types of supported photonic modes by these material.
5.2.1 Uniaxial PhP material
Most anisotropic PhP material would be of uniaxial type, with different permittivity responses, ε|| and ε⊥, for an electric field parallel and perpendicular to the c-axis; thus the c-axis is also the optical axis. Propagating EM modes through uniaxial PhP media have been studied since the 1960s . In analogy to bulk EM modes supported in uniaxial crystals  these modes were called ordinary PhPs (O-PhPs) and extraordinary PhPs (EO-PhPs) with respective wave dispersion relations:
It was not however until the rise of the metamaterials field, and in particular until after the conception of hyperbolic metamaterials , , that the truly “extraordinary” properties of the EO-PhPs were uncovered. Metamaterials spurred the design of man-made EM meta-atoms yielding vast and unusual effective EM responses not readily available in nature, including media with permittivity of opposite sign along different directions (indefinite media) . Later, in 2010, da Silva et al. , noticed that these curious indefinite permittivity responses are actually available in natural quartz by virtue of its anisotropic PhP response. In particular, within a distinct frequency regime, quartz has ε||ε⊥<0, causing a hyperbolic dispersion form for the EO-PhPs (Eq. 39) and in turn negative refraction , , a phenomenon which up to that point had been associated only with man-made artificial metamaterials . An example from  of such EO-PhP-enabled negative refraction in natural crystals is shown in Figure 14; this negative refraction also enables superfocusing effects ,  similar to those reported with artificial hyperbolic metamaterials , .
Generally an arbitrary EM wave incident on an anisotropic uniaxial crystal with a general interface cut couples to both ordinary and extraordinary modes , ; for PhP crystals this would mean both O-PhPs and EO-PhPs would simultaneously couple leading to birefringent PhP propagation . However, in most frequently encountered situations, the interface plane is perpendicular to the c-axis. For this particular interface cut, the O-PhP and EO-PhP modes decouple with O-PhP modes becoming accessible only by S-polarized waves (left panel of Figure 15A) and EO-PhP modes becoming accessible only by P-polarized waves (right panel of Figure 15A) , .
In Figure 15B we summarize the corresponding wave dispersion forms for these cases, determined from Eqs. (36) and (37) as applicable per case; we consider all the possible permittivity value combinations available throughout the PhP regime with widely used materials. We depict in Figure 15C the graphical representations13 of the distinct possible wave dispersion forms. Both types of the hyperbolic dispersion forms, type I and type II (both seen in Figure 15C), are available for different frequency ranges. These are the familiar dispersion forms observed in hyperbolic metamaterials , , , . Aside from negative refraction and superfocusing effects reported in , ,  these PhP-enabled hyperbolic dispersions in natural PhP materials can bring to the IR many of the other exotic effects associated with hyperbolic metamaterials, such us dark-field imaging or a Purcell enhancement ,  from IR quantum dots .
5.2.2 Biaxial PhP material
Certain materials may additionally exhibit in-plane anisotropy, i.e. the PhP response is also anisotropic in the plane perpendicular to the stacking direction, as for example, in α-MoO3 (see also Figure 8 and Table 5). The photonic responses of biaxial media is much richer and a lot more complex than that of uniaxial media. Similar to uniaxial media, the wave dispersion of waves propagating in a non-magnetic biaxial medium, is determined by Maxwell’s equations mandating that the supported modes must satisfy Det[A]=0, with [A] being a 3×3 matrix with elements:
with i=1, 2, 3 and
and i, j=1, 2, 3; δij is Kronecker’s delta and the running indices correspond to the Cartesian directions as: 1↔x, 2↔y, and 3↔z.
However, unlike the case of uniaxial media, the wave dispersion relations, ω(k), resulting from solving the condition above cannot be represented by simple expressions like Eqs. (38) and (39). Given the vast possibilities in biaxial PhP responses, it goes beyond the scope of this review to identify all possible photonic mode forms that are supported. To give the reader a flavor for the unconventional possibilities with PhP biaxial media we showcase an example in Figure 16 that corresponds to εx<0, εy>0, εz>0 with εy≠εz. This case is entirely possible with α-MoO3 in the spectral range λ0 ~ 11–12 μm (free-space wavelength). The corresponding depicted equi-frequency surface (EFS) shows an unusual shape resembling type I hyperbolic dispersion joined with an ellipsoidal dispersion. This unusual EFS has been recently reported for α-MoO3 (see ). The EFS of Figure 16 depicts all available modes in the entire wave vector space for a particular frequency, ω. However, the associated electric field corresponding to a certain wave vector k at a certain frequency, ω, must also be examined to determine whether the mode is or is not accessible by a certain light excitation.
In Section 5 we reviewed the different types of photonic responses enabled by PhPs. In the following section we review the geometries which in synergy with PhP responses lead to confined EM modes which form the backbone of many of the diverse reported IR PhP-based capabilities and phenomena.
6 The essential building blocks of IR phonon-polaritonics
The vast range of amazing phenomena in the field of plasmonic optics all share a common thread: they are build upon the highly-confined optical states at either metal/dielectric planar interfaces or around metallic nanoparticles. Conversely, the essential building blocks relevant to IR photonics would be highly confined photonic modes around planar interfaces or micron/sub-micron particles of pertinent IR material like alternative IR plasmonic materials , ,  or PhP material which is the subject of this review. Therefore, we review in this section the fundamentals of localized modes around PhP material interfaces or PhP particles that set the foundations of the phonon-polaritonics field and present examples of current related research, new phenomena and applications that have emerged. We also discuss the similarities and difference of PhP-based localized modes with their plasmonic counterparts. Furthermore, we highlight new types of localized modes emerging from exotic EM dispersions, such as hyperbolic EM dispersions, which have no counterpart in the optical regime with natural material. The latter have been attracting increasing attention and are expected to play a central role in PhP-based IR photonics.
6.1 Confined PhP modes at planar interfaces and ultrathin films
6.1.1 Ordinary surface phonon-polaritons (O-SPhPs)
Localized quasi-bound EM modes14 may be supported at the surface of PhP materials, with Re(ε)<0. These modes share similar characteristics to SPPs  and are so analogously called surface phonon polaritons (SPhPs) , .15 The dispersion relation for SPhPs, relating their lateral wave vector kx, with their frequency ω, is the same as for SPPs given by :
with c being the speed of light, ε(ω) the complex frequency-dependent permittivity response of an isotropic PhP medium and εd the permittivity of the dielectric material interfacing with the PhP medium.16
SPhPs are supported within the PhP gap/reststrahlen band where Re(ε)<0. The black-solid line in Figure 17A is a graphical representation of the dispersion relation of Eq. (42) for an example SPhP mode supported at the interface between air and the PhP material of Figure 13. The frequency ω is scaled with the TO phonon frequency, ωTO, and the real part of the wave vector is represented with the dimensionless quantity β=Re(kx)c/ω; thus the vertical brown-dashed line at β=1 represents the lightline. Furthermore, the electric-field landscape (x-component) for the mode designated with the green rectangle, shown in the upper part of Figure 17, suggests a well-confined surface mode around the interface of significant propagation range.
The quality of confinement and propagation range are typically in competition for surface waves at interfaces of lossy materials, with Re(ε)<0 , . Different quantifiable measures, i.e. figures of merit (FOM), to capture the performance of these confined surface waves have been hitherto introduced. Here we consider for the SPhP a FOM with the intent to capture simultaneously both the confinement quality and propagation sustainability of the mode, given by:
The defined FOM represents the ratio of the distance the surface wave travels along the interface until its electric field strength drops by a factor of 1/e over the perpendicular-to-the-interface distance at which the electric field strength drops by a factor of 1/e inside the dielectric medium. The FOM value from Eq. (43) should be used in conjunction with the dispersion relation to characterize the overall quality of an SPhP mode. For example, dispersion regions which are very close to or above the lightline certainly signify a weak confinement to the interface; thus any relatively high FOM values in these region only signify a large propagation range. The FOM for the SPhP of Figure 17A is shown in Figure 17B with a black solid line.
To obtain a sense of how SPhPs relate to SPPs, we also include in Figure 17A and B the dispersion relation and FOM-values respectively of an SPP mode supported at the air-Au interface (see also Figure 2A), where realistic optical parameters are used for Au (see Figure 1). To make such comparison, first the SPP’s frequency ω is scaled with the plasma frequency, ωp. As we discussed in Section 2, SPPs in the mid- and longwavelength IR are essentially grazing waves loosely confined at the interface. So we find it useful to compare the IR SPhPs with their SPPs counterparts but around telecomm and visible frequencies. For this reason, the SPP’s ω/ωp value is shown in Figure 17A after being multiplied by a suitable factor yielding, respectively, the red-dashed curves, representing the SPP mode around telecom frequencies, and the cyan curve, representing the SPP mode around visible frequencies.
It appears the SPPs on the surface of Au are generally closer to the lightline in comparison to SPhPs. As with respect to FOM values, the telecom-wavelength SPP modes show a high FOM- but that would mean only a low propagation loss as confinement is ultra-weak for the telecom SPP mode which essentially follows the lightline. The more confined visible SPP mode is clearly outperformed by its SPhP counterpart for most of the reststrahlen band range. This superior performance of the SPhPs modes, in terms of strong confinement and low propagation loss, occurs around the middle of the reststrahlen band away from the TO and LO phonon frequencies. Actually, it is interesting to observe that unlike its SPP counterpart, the SPhP mode runs through a broader range of β values,17 even crossing through the lightline into β values less than 1 at frequencies approaching the LO phonon frequency, ωLO [Foteinopoulou S, unpublished]. It is worth mentioning that such SPhP dispersion crossing into the above-the-lightline region was first experimentally observed in the 1970s for the case of α-quartz , . However, as the FOM shown in Figure 17B attests, high performing SPhPs are associated with the portion of the dispersion curve below the lightline and around the middle of the reststrahlen band.
SPhPs have been utilized to enhance near-field heat transfer , , ,  (see also Figure 18A), as well as to enhance SHG signals  (see Figure 18B). Theoretical calculations have also shown  that SPhPs can enhance the inherent thermal conductivity of a material. SPhPs substrates may enhance coupling with biomolecular vibrational modes and are thus relevant to SEIRA spectroscopy . Moreover, when SPhPs are accessed with an attenuated total internal reflection (ATR)-type prism  ultra-sharp peaks can emerge in the reflectance spectra that can be used for the IR analog  of surface plasmon resonance (SPR) sensors . Pairing up SPhP modes with other physical systems such as phase-change materials , , , 2D materials with gate-tunable properties , , or plasmonic , ,  and dielectric  waveguides, as well as IR active media such as quantum dots  and quantum wells including strained superlattices , unleashes many more possibilities for harnessing IR light pertinent to integrated and/or tunable infrared devices, IR imaging systems ,  and quantum-well IR photodetectors , . We discuss more examples of these hybrid systems in Section 9.
6.1.2 Extraordinary surface phonon-polaritons (EO-SPhPs)
The indefinite anisotropic PhP responses, as discussed in Section 5.2, unlock possibilities for unprecedented types of confined EM modes at PhP material interfaces; this is an unfamiliar territory with no analogue in plasmonics where metals have isotropic negative permittivity responses. The most trivial case involves a uniaxial PhP material in a regime where the permittivity is negative both along and in the normal plane to the c-axis. In such a case, assuming the interface is along the normal to c-axis plane, a minor adjustment to Eq. (42) accounts for such anisotropy and determines the dispersion for the SPhP modes; this can be seen, for example, in , where SPhPs are excited on the surface of 4H-SiC for a frequency within the associated reststrahlen bands of both the parallel and normal to the c-axis directions .
However, the in-plane anisotropy present in certain bi-axial PhP material, such as α-MoO3, introduces entirely new avenues for confining EM waves at PhP-material surfaces beyond the familiar SPPs, such as Dyakonov surface waves . Dyakonov type surface waves are also supported in uniaxial hyperbolic PhP material for an interface cut along a plane that includes the optical axis . Actually, very recently such atypical surface confined waves have been theoretically predicted and experimentally observed on the edge face of an hBN-crystal flake excited by the near-field emanating from a scattering-scanning near-field optical microscopy s-SNOM18 tip (see mode designated as hyperbolic surface polariton (HSP) in Figure 19 after ). Alfaro-Mozaz et al.  also report antenna behavior for an hBN sub-micron strip that is enabled by hybridization of the HSP modes between the edge faces that lead to sharp resonances of a quality factor, Q, of around 100. Moreover, the guided modes on the hBN strips have shown to couple strongly with the molecular vibrations of organic semiconductors making this system ideal for SEIRA spectroscopy .
6.1.3 Ordinary guided PhPs (O-GPhPs): Fuchs-Kliewer interface polaritons
Different types of guided waves may arise at PhP-material thin films. An ultrathin slab comprising an isotropic PhP-material will behave, in the reststrahlen band, analogously to an insulator-metal-insulator waveguide ,  in visible/near-IR frequencies. In particular, the SPhP modes at each interface will hybridize splitting into two bonding/anti-bonding-like branches just like SPs in plasmonic thin films . These coupled SPhP modes at the interfaces of finite slabs were first studied by Fuchs and Kliewer  and are thus sometimes referred to as Fuchs-Kliewer phonons or Fuchs-Kliewer interface polaritons. Ultra-thin PhP-materials slabs supporting such Fuchs-Kliewer polaritons can operate accordingly as IR waveguides. Moreover, the Fuchs-Kliewer polaritons can be utilized for IR dark-field transfer analogous to that in Pendry’s silver lens, also known as the poor man’s superlens . Actually, this phenomenon has been confirmed by Taubner et al. in  where a SiC slab supporting SPhPs functions as an IR superlens.
6.1.4 Extraordinary guided PhPs (EO-GPhPs)
The indefinite permittivity responses  in natural PhP material showcased in Figures 15 and 16 underpin exotic IR guided modes with unconventional dispersion relations. Let us assume a thin slab waveguide with an indefinite uniaxial, i.e. hyperbolic, PhP response with in-plane permittivity, ε⊥<0, and permittivity parallel to the c-axis ε||>0. The dispersion relation for supported guided modes propagating along the plane with a complex wave vector kt has been given as (after , , ):
where εsub is the substrate permittivity and l represents the order of the supported guided mode. An example of such a hyperbolic GPhP is depicted in Figure 19; it is the mode with the fringed pattern on the planar face of the hBN flake designated as a hyperbolic polariton (HP) (after ).
The hyperbolic guided mode of Eq. (44) can be launched onto the PhP-material layer by scattering IR light onto an atomic force microscopy (AFM) tip and was observed for the first time in hBN by the Basov group  with IR imaging from s-SNOM . The Basov group’s pioneering work spurred numerous investigations aiming to further understand such curious guided modes in hBN, e.g. see , , , . Recent reports on this unconventional guided mode include works that show that it can be also launched by scattering light onto the natural defects of the hBN flake such as folds  as well as its realization in suspended hBN as thin as 15 nm . The suspended platform exhibits a clear advantage of reduced damping in comparison to the hBN-on-SiO2 system . The hyperbolic-type GPhPs in hBN stimulated a strong interest to explore these extraordinary guided modes in other PhP material. Very recently, previously unseen guided modes have been reported  in biaxial α-MoO3; α-MoO3 is anisotropic and can be indefinite by virtue of its PhP response on the interface plane (see also Figure 16).
These hyperbolic guided modes in hBN can be utilized for super-focusing/imaging ,  as well as ultra-slow signal propagation , . hBN guided modes have also been shown to boost the sensitivity of graphene photodetectors , an effect that may lead to high resolution IR imaging systems. The quality of the hyperbolic GPhP in hBN can be further tailored by controlling the Boron isotopic purity in the hBN sample. Specifically, processed hBN samples with either 10B or 11B isotopic purity were very recently utilized to achieve a substantially longer propagation range (lower loss) for the hyperbolic guided mode in comparison with the natural species .
6.2 Particle PhPs
Light incident on plasmonic nanoparticles may give rise to highly localized resonant modes that have been termed as particle plasmons , LSPs  or Mie plasmons , , because of the Mie-type light scattering involved . Nanoantennae , high-Q Fano resonances for biosensing , aperture-less near-field optical microscopy  and spasers , etc. form only an indicative list of a plethora of fantastic effects and applications enabled by particle plasmons. Conversely, localized modes “molded” around micron-sized PhP particles can be the basis for shaping IR light for targeted IR applications; analogously we will refer to these modes as particle PhPs.
However, unlike particle plasmons, particle PhPs may arise from distinctively different material responses which are available in the PhP regime; these can range from plasmonic-like responses characteristic of noble metals in the visible spectrum to HRI responses (see Section 5.1 and Figure 13). Furthermore, indefinite uniaxial or biaxial responses (see Section 5.2 and Figures 15 and 16) are also available with PhP material. These responses do not have any counterpart in dielectrics outside the PhP regime or metals. Thus, in the following we will be referring to particle PhPs as ordinary if connected to plasmonic-type or HRI isotropic photonic responses and as extraordinary if connected to indefinite anisotropic  PhP photonic responses as discussed in Section 5.2.
6.2.1 Ordinary particle PhPs
Just like the case of plasmonic nanoparticles, the resonant response of PhP microparticles to incoming light can be captured by the extinction efficiency, Qext , which is the extinction cross-section  normalized to the particle’s cross-sectional area. The extinction efficiency involves the particle’s response both in terms of scattering and absorption and so gives a more complete picture for the resonant response of particles comprising lossy photonic materials, such as PhP materials. The extinction efficiency sums up the effect of all possible resonant orders, that may be excited at the microparticle (multipoles). Assuming a plane wave of angular frequency, ω, is incident on a microsphere of radius R immersed in air (vacuum), the extinction efficiency, Qext, is equal to , , , :
with n=1 representing the dipolar resonant order, n=2 the quadrupolar resonant order and so on.
For each order the mode is a superposition of two characteristically different field distributions each having a strength represented by the coefficients an and bn, determined by the EM boundary conditions at the sphere’s surface.19 Actually, an represents the strength of a field distribution having an electric field with a strong radial component (implying a strong surface charge density) and a magnetic field that is purely tangential. Conversely, the coefficient bn represents the strength of a field distribution having a magnetic field with a strong radial component and an electric field that is purely tangential. Accordingly, an and bn have been associated with an electric- and magnetic-type resonant response, respectively .
We plot in Figure 20A the extinction efficiency for an example case. We do so in order to showcase the existence of distinct types of ordinary particle PhPs emanating from the different isotropic photonic responses in the PhP regime: plasmonic-like and HRI responses, respectively. The particular example we consider here involves a microsphere made of the PhP material of Figure 13, with radius, R, equal to 0.065 λTO, with λTO being the transverse-optical phonon free-space wavelength. We observe two resonances, hence two particle PhPs: one within the plasmonic permittivity regime, that we will be referring to as the type A particle PhP, and one within the HRI regime, that we will be referring to as the type B particle PhP.
The type A ordinary particle PhP is similar to the LSP; and hence is typically referred to as the localized surface phonon polariton (LSPhPs) , , , . Indeed, its dominant response comes from the dipolar electric-type term, a1, yielding an electric field landscape shown in Figure 20C that resembles the familiar dipolar response of plasmonic nanospheres in the visible spectrum (see, e.g. LSP1 mode of Figure 2E). The observed LSPhP resonance is excited near the Fröhlich condition , .20 This is expected given the deep-subwavelength size of the microsphere. Just like in the case of LSPs, as the microsphere becomes larger the electric dipolar resonance would occur increasingly further from the Fröhlich condition and an electric quadrupolar resonance would emerge. Eventually, for microsphere sizes of the order of the wavelength the near-field response to the incoming EM field involves multiple orders and both electric- and magnetic-type field distributions.
The type B ordinary particle PhP shown in Figure 20B resembles the Mie-resonant responses of dielectric particles . For this particular example, the electric field involves contributions both from the electric dipole (ED) term, a1, and the magnetic dipole (MD) term, b1.21 Therefore the type B particle PhPs may be referred to as HRI-Mie PhPs , in order to distinguish these from LSPhPs; however unlike pure dielectric Mie modes, these type B particle PhPs have a lossy character. In fact, the high extinction efficiency that is observed in Figure 20A is in large part a result of a high absorption efficiency, as opposed to a high scattering efficiency. HRI-Mie PhPs have not, thus far, been widely investigated.
Actually, to-date it is mostly type A particle PhPs, i.e. LSPhPs, that have been widely studied. Theoretical reports of LSPhs involve infinite or finite-height cylinders with circular, elliptical and rectangular cross-section cylinders , ,  as well as disks . The existence of LSPhs has been experimentally verified; in fact LSPhPs have been spectroscopically characterized and/or imaged in sub-micron pillars , , microdisks  (see Figure 21A) or nanoparticle clusters . The field enhancements in LSPhP excitations may be utilized for thermal emission , , , , enhancing non-linear phenomena such as second harmonic generation, for example, see  and in sensing platforms like SEIRA . LSPhP properties can be further tuned by harnessing their collective interactions in a 2D array lattice , just like with lattice LSPs . Finally, particle PhPs in SiO2 microspheres have been shown to enable enhanced near-field heat transfer from the microsphere to a SiO2 plane by coupling to the respective plane SPhPs ; likely both LSPhPs and HRI Mie PhPs contribute to this phenomenon in this case.
6.2.2 Extraordinary (EO) particle PhPs
PhPs made available unconventional “metamaterial” photonic responses in natural material. As discussed in Section 5.2, permittivity may be indefinite  in PhP media, i.e. have a different sign in two of the crystallographic directions. The corresponding photonic dispersions can be extraordinary, of hyperbolic type or other curious complex shapes as in the case of the biaxial medium of Figure 16. Accordingly, the discovery of indefinite responses in PhP media ,  brought out an entire exotic garden of related localized particle PhPs which are currently a subject of increased research attention. In this section we use the term of EO-particle PhPs to be inclusive of particle PhPs in any anisotropic PhP medium with corresponding EO, i.e. non-isotropic, photonic dispersions. However, truly unconventional resonant signatures not possible with isotropic particles appear to emerge only in certain instances of microparticles comprising indefinite PhP media.
For example, particle PhPs have been theoretically and experimentally studied in a bow-tie nanoantenna comprising anisotropic 4H-SiC with the purpose of achieving spectral selectivity of the SiC’s emission spectrum. The observed resonant PhP modes in each triangular structure of the bow-tie nanoantenna resemble the multi-resonant response observed at the visible spectrum in triangular plasmonic particles . Moreover hBN-disk dimers , with hBN in the hyperbolic dispersion regime appear to respond somewhat similarly to plasmonic dimers  in the sense of confining and yielding large enhanced fields in the interstitial region between the disks.
On the other hand, particle PhPs observed in hBN nanocones ,  exhibit an extraordinary electric-field landscape with “hot-rings” at deep-subwavelength distances on the nanocone surface (see Figure 21B). These modes are atypical for microparticle resonators comprising isotropic material and originate from the hBN hyperbolic photonic responses that cause unconventional, highly directional, reflections inside the cone structure at its boundary. These peculiar EO-PhPs are promising for shaping resonances with high-quality factors and strong field confinement. Furthermore, very recently, particle PhPs have been also imaged on an α-MoO3 disk exhibiting also a strong deep-subwalength fringe pattern on the disk surface  (see Figure 21C). It can also be seen in Figure 21C that a Fourier-transform of the PhPs curious electric field landscape reveals the underpinning photonic dispersion which emanates from α-MoO3’s indefinite biaxial properties and is similar to the one depicted in Figure 16. Note that, the Cartesian axis system of Figure 21C has the y and z axis interchanged in comparison with the one of Figure 16.
6.3 New physics with particle PhPs: LO-TO phonon mixing and atomic-scale thermal emitters
Both O-PhPs and EO-PhPs are expected to be major protagonists in conceiving new platforms to shape, divert, absorb or emit IR-light as mandated by specific applications. Very importantly, these systems may also unveil new PhP physics. For example, although LO phonon oscillations are known to be inaccessible by IR light in planar geometries  it has been very recently reported  that microparticle resonators comprising certain anisotropic long-cell polytypes of a PhP material, like 4H-SiC, may enable making zone-folded LO phonons accessible with IR light excitations. In turn, this enables new hybridized longitudinal-transverse PhP excitations. This recent study will inspire new investigations for the understanding, control and realization of these hybrid excitations in other PhP material systems. As Gubbin et al. report , these mixed excitations pave the way for the realization of electrically pumped mid-IR emitters. Furthermore, particle PhPs can be brought down to the atomic-scale in low-dimensional particles in which cases atomistic models would be required to capture their PhP responses. Recently, Venkataram et al.  developed an atomistic theory to capture PhP-mediated emission and heat-transfer effects between elongated low-dimensional molecules and macroscopic objects that cannot be interpreted with classical models.
7 Patterned PhP materials: metamaterials and mesophotonic platforms
As we discussed in Section 6, localized PhP modes at planar interfaces or particles can offer a plethora of possibilities in IR light control. There are however certain patterning routes of PhP material that build further from these localized modes by harnessing judiciously their interaction. Moreover, non-resonant PhP responses can be also brought together into platforms with new capabilities. We briefly review below quintessential examples of such patterning routes along with the corresponding characteristic length scales of their feature sizes and discuss how they move forward the already vast potential of the systems of Section 6.
7.1 PhP-based metamaterials
7.1.1 Hyperbolic PhP metamaterials
While PhP media can be “natural metamaterials” with hyperbolic photonic responses of type I or type II, as we discussed in Section 6.2, it is certainly worthwhile to look at the metamaterial patterning route to a hyperbolic response with PhP material. This way a broader parameter space and spectral availability of the hyperbolic photonic dispersion can be accessed and tailored by design at will. The associated building-block units should comprise a PhP material in its reststrahlen band and should be deeply sub-wavelength, with a size of the order of λ0/100, with λ0 being the free-space wavelength , , . When these elemental units are put in periodic arrangements they respond as “meta-atoms” to incoming IR light making a new material with effective photonic properties. The extremely, deep-subwavelength size for the meta-atom is a requirement for a high quality effective hyperbolic response  that sustains up to large lateral wave vectors .
Typical arrangements reported with PhP material involve either periodic-rod , ,  or lamellar structures , ,  (Figure 22A). At such deep-subwavelength dimensions it can be assumed the electric-field is quasi-static within each constituent. Then the respective effective permittivity responses parallel and perpendicular to the optical axis, ε|| and ε⊥, can be obtained by established mixing formulae , , ,  that take into account the EM boundary conditions that apply at the two different electric field orientations. In particular, for the case of periodic rods , :
where ε1, ε2 are the isotropic permittivities of the two material constituents and f is the volumetric ratio of constituent ‘1’ with respect to the volume of the unit meta-atom. The optical-axis direction of the effective metamaterial for both systems has been designated in Figure 22A.
An indefinite response of Re(ε||)· Re(ε⊥)<0 is needed to obtain hyperbolic photonic responses of either type. For this to be possible the constituents should have an opposite permittivity, i.e. Re(ε1)· Re(ε2)<0. One possibility is for one constituent to be a PhP material operating within its reststrahlen band with the other constituent being a dielectric. Alternatively, both constituents can be a PhP material with one operating within and the other outside their respective reststrahlen bands. In principle, in a dielectric-PhP material composite both type I and type II hyperbolic responses are accessible at different spectral bands , . When both constituents are PhP materials, several spectral bands of type I and type II hyperbolic behavior may emerge; this behavior has been reported for self-organized eutectic ionic crystal mixtures such as NaCl-LiF and NaCl-KCl , .
PhP-based hyperbolic metamaterials can be instrumental in enabling new components for table-top IR optics. For example,  reports that PhP-based hyperbolic metamaterials can function as both polarization filters and converters. Furthermore, the IR lensing and beam-steering capabilities reported for natural PhP hyperbolic crystals , , ,  can be extended to PhP hyperbolic metamaterials while adding at-will design flexibility over the operational bandwidth and steering angles. Recently, dark-field imaging in the optical spectrum was shown to be possible with a type-I/typeII hyperbolic metamaterial heterostructure . This capability can be extended to the IR with corresponding judiciously designed heterostructures of either two PhP hyperbolic metamaterials or a PhP hyperbolic metamaterial and a natural PhP hyberbolic crystal and can benefit IR imaging systems , . Finally,  reports that PhP hyperbolic metamaterial may exhibit enhanced super-planckian thermal emission.
7.1.2 Highly tailorable anisotropic PhP metamaterials
The naturally available indefinite responses in PhP materials open-up the possibility of easily implementable and highly tailorable anisotropic photonic metamaterials. Li et al.  very recently reported, with a spectroscopic-based analysis, that patterning hBN as a grating with a deeply sub-wavelength pitch of 100 nm can yield effective indefinite biaxial responses (see Figure 22B) of the kind observed in α-MoO3 (see also Figure 16) but highly tailorable with patterning. Note, that this type of metamaterial is unique to PhP-based systems. It would not be possible to obtain an effective indefinite biaxial behavior with a simple grating patterning in a plasmonic material. In fact, for an analogous biaxial metamaterial to be realized with optical plasmonic systems, a substantially more challenging two-step patterning would be required involving a grating structuring on an already patterned metallodielectric lamellar structure (e.g. see , ).
7.1.3 Magnetic PhP metamaterials
As the meta-atom sizes of periodic-rod structures increases the quasi-static response picture ceases to be valid and along with it the effective permittivity response models of Eqs. (46) and (47). That does not necessarily mean however, that such systems cease to behave as a metamaterial. In fact, for HRI rods the Mie resonances occur at lower frequencies . Effectively this means that HRI-rod periodic structures may support Mie excitations even when building-block units are significantly smaller than the free-space wavelength , . On the other hand, rigorous effective medium tests have shown ,  that periodic dielectric systems with building-block units of the order of ~λ0/10, with λ0 being the free-space wavelength, generally behave as metamaterials. However, their effective photonic responses are generally not captured by the simple standard effective medium theories of Eqs. (46) and (47) . The most general approach to determine the effective EM constitutive properties of periodic PhP-material systems are suitable EM retrieval processes based on their transmission/reflection spectra (see Appendix D of  or  for a retrieval method established especially for these type of systems). HRI-rod Mie resonances and their interaction with the Bragg resonances in a periodic lattice can be tailored to yield an effective magnetic behavior for the periodic-rod system (see ).
PhP media possess a HRI in the vicinity of the TO phonon frequency (see Figure 13). Accordingly, periodic structures with PhP rods, depicted in Figure 22C, are ideal systems to realize the aforementioned magnetic metamaterial for IR frequencies. The first theoretical report of artificial magnetism in PhP-media periodic composites, with an effective negative permeability, μ, came from Huang et al.  in a system utilizing LiTaO3 rods. Following this work, a negative effective μ has been reported in an fcc lattice of LiTaO3  as well as in a periodic SiC-rod system . Furthermore, an effective magnetic response, albeit weaker, was reported in  in a LiF/NaCl composite, a system that is easily realizable with self-assembled eutectic mixtures . Permeability μ negative media (MNG) in heterostructures with permittivity ε negative media (ENG) may act as a photonic analogue of an LC circuit . Accordingly PhP-based magnetic metamaterials can be useful for the construction of resonator elements in integrated IR photonic systems serving as bandpass filters or modulators.
7.2 PhP-based mesophotonic systems
We use the term “mesophotonic systems” here to describe systems in the sub-wavelength regime where geometric optics fails but that may not respond with effective photonic responses as metamaterials. Typical feature/building-block sizes in this class of system are ~λ0/2, with λ0 being the free-space wavelength. We outline below couple of characteristic platforms in this class that have utilized PhP material.
7.2.1 Particle-PhP arrays: metasurfaces
One-dimensional (1D) or 2D arrays of periodic or more general arrangements of particle PhP resonators can be envisioned constituting essentially PhP-based metasurfaces [see, e.g. depiction in Figure 23A]. This line of research has not been widely pursued to the best of our knowledge but is very promising for sculpting, diverting, absorbing/emitting IR light and biosensing methods like SEIRA . For example, one could envision an LSPhP-based platforms for IR beam steering operating analogously to LSP-based ultra-refractive metasurfaces . Then, established design routes of all-dielectric Mie-based metasurfaces , , such as reflection suppression relying on the lattice Kerker effect , can inspire new Mie-PhP-based platforms for near-perfect absorption/emission . Furthermore, new extraordinary particle PhPs set a pathway for unexplored metasurface platforms that have no analogue in IR metallic , plasmonic  or all-dielectric ,  metasurfaces. Collective phenomena analogous to lattice LSPs  will be important to all particle PhP-based metasurfaces.
Currently, specific 1D , ,  and 2D  arrangements of unit or composite particle-PhP resonators have primarily been researched for harnessing enhanced and/or highly-directional thermal emission. As an example we discuss the platform of  seen in Figure 24. LSPhPs in rectangular cross-section rods arranged in a 1D array of micropyramids have been utilized to achieve polarization selective near-perfect absorption that covers about half the SiC reststrahlen band. The micropyramid array responds like moth eyes with a near-zero albedo  and thus near-perfect absorption when light impinges from the tip side of the pyramid, sequentially coupling from the smaller cross-section to the larger cross-section rods. By virtue of the directional Kirchoff’s law  this structure then behaves as a near-unidirectional emitter. The phenomenon can be relevant for highly efficient polarized spectrally selective globar thermal sources which can benefit MIR biomedical imaging . LSPhP-based absorption/emission engineering , ,  is promising for many applications including IR-imaging bolometers , radiative cooling ,  and optothermal devices .
7.2.2 Perforated PhP-material arrays
Holes in a PhP material matrix in the metallic regime, i.e. when Re(ε)<0, would also support site-resonances analogous to holes in metals. Accordingly, the long-range interaction between such site resonances in perforated PhP-material arrays (see depiction in Figure 23A) can yield enhanced optical transmission (EOT) phenomena similar to metal perforated arrays , . EOT in a perforated SiC array has been reported in Refs. , . The interaction between the site resonances at the holes and extended surface-Bloch modes  also enables Fano interference  in these platforms . This class of systems can be useful to modulate/spectrally select emission spectra pertinent to imaging ,  or thermophotovoltaic applications . Also, the stronger electric fields associated with the site resonances may be suitable to enhance molecular absorption fingerprints for SEIRA methods .
7.2.3 Photonic-crystal-based PhP systems
In the photonic crystal (PC) regime, structures do not respond with effective constitutive properties. The behavior of the system is a result of a combination of any localized resonances (e.g. Mie resonances) within the building-block unit and Bragg scattering between the building-block units . PCs with PhP constituents can show extreme field re-configurability within an ultra-narrow spectral window . Specifically, if the photonic band-gap of a periodic-rod PC is aligned with the PhP gap of the PhP-rod constituent material extreme field reconfiguration occurs when crossing the TO phonon frequency: the electric-field is confined within the rods just below the TO phonon frequency and it becomes completely expelled from the rods just above the TO phonon frequency (see indicative structural parameter range for this phenomenon in the left panel of Figure 23B after ). The effect may inspire new switchable/reconfigurable PhP-based IR devices.
Furthermore, PhP-based multilayer media are suitable to harness absorption/emission , ,  relevant to bolometers  and thermo-photovoltaic applications . While PhP media in the reststrahlen band are completely reflective, slicing them up into ultra-thin slices separated by much thicker dielectric spacers can achieve near-100% absorption (see  for the case of SiC). In aperiodic arrangements this effect can become highly broadband, covering the entire reststrahlen band . If the multilayer operates near-the-band edge, it provides near-unity reflectivity, i.e. it becomes a Bragg reflector. However,  reports that an extremely sub-wavelength PhP material layer of judiciously chosen thickness, placed on top of one unit-cell of such a PhP medium-dielectric Bragg reflector, can suppress reflection while concentrating the electric field on the ultrathin PhP layer. In particular, a spectrally-selective absorption as high as 92% occurring in the ultrathin PhP layer, SiC in this case, was theoretically demonstrated (see schematic in the middle panel of Figure 23B).
We note in passing that,  also reports near-perfect absorption on a thin layer of PhP material by utilizing destructive reflection interference in a three-layer system comprising the top PhP material layer, a dielectric spacer and a frequency-selective backreflector. The PhP material layer of  has still a sub-wavelength thickness of ~λ0/140, albeit quite larger in comparison to the PhP material top layer of  which is about ~λ0/1000, with λ0 being the free-space wavelength. Additionally, EM wave-interference phenomena for absorption management in PhP material have also been realized in geometries interfacing a PhP material layer and an all-dielectric 1D photonic crystal (multilayer). In particular,  reports on one-way near-perfect absorption in a λ0/25 thick GaAs layer operating in its reststrahlen band that rests on an all-dielectric Si/SiO2 stack.
Sapphire is also gaining attention as an IR material, given that it also possesses a spectral regime where it behaves as a natural hyperbolic material due to anisotropic PhPs (see also Sections 4.3 and 5.2.1). Narimanov  reports periodic multilayers, i.e. 1D PCs, comprising sapphire in the hyperbolic dispersion regime. These type of PCs with hyperbolic material constituents, termed as photonic hypercrystals [see schematic in Figure 23B (right panel)], may give rise to new types of low-loss ultra-confined surface modes , . These low-loss surface modes can be utilized for new capabilities in longer-wavelength infrared (LWIR) spectroscopic sensing relevant to part of the IR fingerprint spectrum and smaller molecules , .
8 Dynamic PhP material for tunable, switchable and reconfigurable IR platforms
In this section we discuss avenues to tune the natural PhP material responses to enable tuning or extension of the operational bandwidth of pertinent modes and platforms that we discussed in the previous sections. The tunability would permit filling spectral gaps of available PhP material relevant to targeted applications. Furthermore, we discuss possibilities to tune and switch the PhP material responses with external stimuli that would enable new IR tunable, switchable or reconfigurable platforms.
8.1 Doping-tunable PhP permittivity in semiconductors: the LOPC effect
The carriers in semiconductors provide a tuning “knob” for the ionic permittivity response seen in Figure 6 via their Drude contribution to the full permittivity [see Eqs. (30) and (31)]. For a significant carrier concentration, can approach or become larger than the TO phonon frequency, ωTO . While the TO phonon frequency, ωTO remains unaffected by the presence of free carriers, as these do not influence the physical properties of the IR activated lattice vibrations, this is not true for the LO phonon frequency, ωLO. With carriers present, there are two possible frequencies, where the permittivity can become zero thus implying that the original LO optical phonon frequency ωLO splits into two branches: ωLO+ (blue shifted) and ωLO− (red shifted). This splitting of the LO phonon frequency into two branches  has been termed as the longitudinal optical phonon-plasmon coupling (LOPC) effect.
The LOPC effect arises from the ionic and carrier motion being coupled together via the macroscopic electric field. Although known since the 1960s , the significance of the LOPC effect for IR light-matter interactions was realized only recently when utilized in planar  or microstructured platforms ,  for tunable IR light confinement, enhancement and absorption. To get a feeling for the impact of the LOPC effect on the photonic responses of semiconductors around their reststrahlen band we showcase a demonstrating example in Figure 25 depicting the full permittivity, εsem, of InP for different carrier concentrations (assuming room temperature).
Doping may be enabled via photoinjection, i.e carrier generation with optical pulsed or continuous wave fields with frequencies above the semiconductor band gap . This has been utilized by  to dynamically tune particle PhPs in arrays of InP and 4H-SiC sub-micron pillars. In the semimetals and semiconductors of the 2D/quasi-2D PhP material of Section 4.4, tunability via the LOPC can be achieved by applying a gate voltage. Note that some materials demonstrate a stronger phonon-plasmon coupling leading to cross-coupled damping coefficients; the permittivity response of these materials is different from the independent plasmon and phonon responses of Eqs. (30) and (31) (see ).
8.2 Tunable PhPs with strong optical fields: unlocking non-linear phonon-polaritonics
Another interesting optical-fields-based mechanism has been very recently discovered, unlocking new avenues for tunable phenomena in PhP materials with optical pulses . Specifically, Cartella et al.  point out that the linear ionic oscillator model (see Section 3) would be valid only for small-to-moderate displacements. A theoretical model developed by the authors of  indicates that the high-frequency permittivity and the ionic effective charge Q (Born-effective charge) may acquire a term that is quadratic on the displacement, u  under the influence of strong optical fields. Their hypothesis is supported by experimental evidence showing that the PhPs in SiC can act as a type of a parametric oscillator yielding gain in SiC’s reflectance spectrum when excited by a strong pump beam. This very recent study sets the foundations for phonon-enabled IR non-linear photonics opening a whole new range of phenomena and capabilities in IR photonics.
8.3 Electric-field tunable PhP responses: the case of ferroelectrics
Discovered in 1921 , ferroelectrics represent a special class of dielectrics with intricate photonic responses stemming from a special acentric crystal structure that however possesses an axis of symmetry . Ferroelectric materials exhibit a residual polarization even in the absence of an applied electric field; with application of an electric field the electric polarization follows a hysteresis similar to the magnetic polarization hysteresis of ferromagnetic materials . Ferroelectric materials lose their curious responses for temperatures above a critical temperature, TC, known as the Curie temperature , , , , . The Curie point, TC, of transition to the ferroelectric phase, can be tuned with the application of stress  or even material thickness , and can be very high, approaching the melting point for many ferroelectric materials .
The explosion and progress in the field of 2D and quasi-2D materials helped also bring ferroelectrics down to the scale of a few nanometers, despite that being initially thought of as impossible . Examples of 2D/quasi-2D ferroelectrics that have been experimentally realized are SnTe, CuInP2S6 (CIPS) as well as α-In2Se3 . Moreover, first principle atomistic calculations revealed additional ferroelectric material that can be stable in 2D form such as SbN , BiP  and β-GeSe . β-GeSe can be versatile for a range of applications, as it also exhibits strong tunability under strain, which can bring its Curie temperature down to the room temperature range.
Ferroelectrics are very interesting from the point of view of PhPs. Going back to the microscopic model of Section 3.2, we see from Eq. (19) that when:
the PhP frequency ωR=ωTO becomes zero while simultaneously the static dielectric constant, εs, becomes infinite thus still obeying the LST relation [Eq. (24)] , . This unstable PhP behavior is known as polarization catastrophe and signifies the transition to the ferroelectric state , , . Accordingly, at the transition point a small change in the electric field yields a giant change in the polarization . Following thermodynamic arguments (Devonshire model) it can be shown that in the ferroelectric state, i.e. for temperatures below the Curie temperature TC, the low-frequency PhP mode (also known as soft-mode) increases in frequency, i.e. “hardens” with the application of a static electric field. Note, however, that the Curie transition temperature itself may increase with the application of a static electric field .
Clearly, the PhP (ionic) photonic responses of ferroelectric material are strongly tunable with a static external electric field. Thus ferroelectrics are desirable components for IR tunable devices; especially those with a Curie point at or above room temperature. Ferro-electric material with DC controlled ionic permittivity responses have been recently utilized for tuning the defect mode of a 1D photonic crystal cavity . Moreover, tunable near-field radiative heat transfer has been demonstrated  that relies on the electric-field tunability of the associated surface PhP modes (see Section 6.1.1) of the ferroelectric material.
8.4 Magnetic-field tunable photonic responses in the PhP regime: the case of n-InSb
The cyclotron motion of the carrier electrons  under a static magnetic field causes a directional change in the full permittivity response of n-doped semiconductors (see , ). The phenomenon becomes relevant to tuning photonic responses in the PhP regime only when the cyclotron frequency, with B= |B | being the magnetic field strength, is in the range of the TO phonon frequency, ωTO for relatively low magnetic fields. Not many semiconductors would meet that requirement. Given its low electron effective mass, n-doped InSb has been predominantly studied for its magnetic-field tunable photonic responses in its PhP regime from the 1970s  until today (see indicative example in Figure 26A).22 Recent related research includes tunable non-reciprocal plasmonic-type devices functional at THz frequencies for photonic isolators , , THz sensors , THz tunable polarization conversion  as well as magnetic-field tunable near-field radiative heat transfer pertinent to ultrafast control of thermal management , .
8.5 Turning PhP photonic responses on and off by phase transitions
VO2 has been attracting increasing attention as a tunable photonic material owing to its transition from an insulating monoclinic to a metallic rutile phase at the temperature of 68°C , . VO2 exhibits multiple PhP oscillations in its insulating phase (see Table 4). When heated beyond the transition temperature, these PhP responses completely disappear and VO2 responds with a metallic Drude-like permittivity.23
In Figure 26B-I and 26B-II we plot the room-temperature permittivity of VO2 in the PhP regime for electric fields parallel and perpendicular, respectively, to its crystallographic direction, αM (ionic permittivity model parameters are from  shown in Table 4). Indeed, we observe from Figure 26B-III that the permittivity changes dramatically past the phase-transition point at a temperature of 87°C where VO2 is in the rutile metallic phase.24 Based on permittivity values extracted from the Basov-group experimental data , we also considered here a standard fitting Drude model that captures well the permittivity in the metallic rutile phase for the frequency range relevant to the PhP regime of the insulating phase. This Drude model’s parameters are ωp=15.67×1014 rad/s and γp=4.92×1014 rad/s; the associated real (imaginary) part of the permittivity is represented with the black-solid (red-dashed) lines in Figure 26B-III.
The ionic-to-Drude permittivity transition in VO2 has been utilized to tune the near-field radiative heat transfer with temperature  as well as to device a microsecond switchable thermal antenna . Both effects rely on switching on/off surface or particle PhPs in VO2. As other phase-change material, like GST , , are gaining interest for photonics in the last years, new PhP switching capabilities beyond VO2 may come to light. In fact, there is recent on-going research on characterizing IR active phonons in GST and other phase-transitioning alloys , . Temperature enabled volatile PhP control in phase-change material would be highly attractive for applications pertinent to thermal homeostasis capability  or for thermal camouflage . Non-volatile PhP switching is pertinent to reconfigurable elements for integrated IR photonics. Furthermore, recent research that enabled gate-controlled non-volatile phase-transition in VO2  via hydrogen ion infiltration offers potentially exciting pathways towards gate-tunable PhP-based IR photonics with VO2.
Of course phase-transition enabled tunable devices can also be envisioned by interfacing PhP material not undergoing a phase transition with other phase-change material not possessing PhPs in the spectrum of interest. We will review the latter type of systems when we discuss hybrid platforms in the following section.
9 Hybrid platforms with PhPs
Hybrid platforms are systems that bring together, in a synergistic manner, capabilities stemming from distinct physical mechanisms. In this manner, “the best of two worlds” is combined towards breaking performance bottlenecks or enabling new functionalities and thus broadening the applications horizon. Hybrid approaches to designing photonic devices are gaining increased popularity across the EM spectrum. For example, 2D semiconductors can be combined with disk resonators or photonic crystals to enhance their photoluminescence spectrum , ; plasmonic nanoantennae can be combined with dielectric microcavities achieving both high confinement and high quality factor, Q . In this section, we discuss how certain PhP modes reviewed in previous sections can be combined with other types of photonic modes and/or materials to form hybrid systems that advance further extreme IR light control.
In Section 3 we discussed that PhPs are hybridized modes coming from strong coupling of IR light and lattice vibrations. PhP-enabled modes, such as SPhPs or GPhPs, may also themselves hybridize with other types of modes, such as SPPs. Such hybridization manifests itself as a strong modification of their mixed dispersion, ω(k). In fact, the original independent modes have dispersions that cross in the (ω, k) space. As a result of the their strong interaction they may convert into two anti-crossed branches both having a mixed character from the two original independent modes. Such anti-crossing is really the quintessential signature of the strong-coupling effect between two different types of modes . One possibility is SPhPs strongly-coupling with SPPs. For example,  reports on the strong coupling between SPhPs, native of an SiC or oxide layer, with SPPs, native of a plasmonic metallic layer. Such strong coupling is exploited to devise an architecture that provides enhancement of near-field radiative heat transfer while at the same time being tunable by application of a gate voltage. The gate voltage tunability is facilitated by the PhP material layer which in fact serves in a dual role: enhancing the near-field heat transfer and providing an insulating spacer between the plasmonic layer and the metallic contact, needed for tuning the properties of the the plasmonic layer with voltage. Furthermore, strong coupling between SPhPs of a SiO2 layer and Mie plasmons of a gold nanoantenna  as well as SPhPs on a PMMA film and designer SPPs of a plasmonic metamaterial  have also been reported. Such effects can find applications in modulators and filters for IR sources.
Naturally, the current interest in 2D and quasi-2D photonic material , , has also spurred investigations of SPP-SPhP as well as SPP-GPhP strong coupling involving composite systems with either or both constituents being atomically-thin materials. Many related systems have been reported for strong coupling between SPP modes in graphene, and several SPhP modes in SiO2 ,25 SPP modes in graphene and the hyperbolic GPhP in hBN , , ,  (e.g. see Figure 27A) or SPP modes in doped BP and GPhPs in SiO2 . This strong-coupling enables modulation of the composite SPP-SPhP or SPP-GPhP photonic dispersion with gate-tunable capability via gating of the graphene layer. Highly-doped bilayer graphene is a uniquely interesting system because it can by itself support SPP-SPhP hybridization. Although monolayer graphene does not support PhPs, bilayer graphene does (see Section 4.4). Then the native longitudinal (transverse) plasmons can couple to the TO (LO) optical phonons . Such hybridization grows stronger with increased doping and shows the signature of strong coupling: anti-crossing of the mixed-mode dispersion bands. The phenomenon has been experimentally verified in ultra-highly doped bilayer graphene via FeCl3 intercalation . These effects open several possibilities for bringing extreme IR light control with 2D tunable devices.
Furthermore, in a graphene-polar semiconductor heterostructure it was reported that GPhPs supported in an ultrathin polar semiconductor layer may fast decay into graphene SPPs thereby getting completely quenched . Also, hot carriers in graphene  or in highly-biased graphene bilayers  can cool via coupling to the GPhP modes of an adjacent hBN layer. 2D material dielectrics, such as TMDs outside their PhP regime, also can affect the properties of SPhPs in an adjacent polar semiconductor layer; in particular they ultra-squeeze the associated SPhPs into tight confinement . Moreover, gate-tunable hybridized SPP-SPhP modes have been reported in a graphene-LiF heterostructure; these may be tailored to outperform, in terms of combined confinement and propagation length, both the individual SPPs supported by the graphene layer and the SPhPs supported on the LiF surface .
Beyond mode interaction, hybridization and strong coupling, atomically-thin PhP material can be interfaced with other atomically-thin material to tailor their photonic responses or achieve new functionalities. For example, a deep-subwavelength graphene/hBN heterostructure can serve as a gate-tunable meta-atom for a metamaterial with adjustable photonic dispersion. The latter metamaterial’s photonic response can change from ellipsoidal to hyperbolic (see Figure 15) making it an ideal medium to strongly modulate the emission spectrum of single-photon emitters . A similar graphene/polar-dielectric (SiO2) lamellar structure acting as a tunable hyperbolic metamaterial  has been very recently experimentally realized; the structure incorporates additional Al, and Al2O3 layers needed to enable the fabrication of the structure. Furthermore, hBN/BP vertical heterostructures have also been considered operating within the PhP regime of hBN and outside the PhP regime of BP. The strong in-plane anisotropy of the BP dielectric enables GPhPs in hBN that are in-plane anisotropic . In addition, when a lateral heterostructure interface is considered between graphene and hBN, SPPs supported on the graphene side may refract negatively into GPhPs on the hBN side at the graphene-hBN interface  (see also Figure 27B). By virtue of such negative refraction a localized emitter source may be focused across the interface similarly to a point source focusing through Pendry’s negative refraction index lens , , making this an attractive platform for imaging applications. Furthermore, very recently tunable IR superlensing has been theoretically reported in a similar system in an hBN-graphene lateral heterojunction resting on a SiC substrate; this tailorable and tunable superlens phenomenon exploits the strong coupling between the SPhP modes supported by SiC and the hyperbolic PhPs supported in hBN .
Phase-change material (PCM), such as VO2 ,  and GST , , are gaining increased attention as photonic materials. In Section 8.5, we discussed how PhP responses can be turned on and off in VO2 as well as on-going research in uncovering PhP responses in GST. However, even outside their PhP regime these material can serve as a volatile or non-volatile photonic environment to tune or reconfigure the PhP responses of PhP-based platforms interfaced with them. For example,  reports a VO2-hBN-graphene vertical heterostructure with asymmetric light transmission and resonant near-perfect absorption capability that can be tuned both via the VO2 insulator-to-metal transition and via gating of the graphene.
Furthermore, judiciously designed patterns on PCMs that comprise the different material phases can be used for reconfigurable beam-shaping and steering of SPhP- or GPhP-waves supported on an adjacent hBN flake . A similar phenomenon has been used in a GST-quartz heterostructure. The reversible switching of GST between amorphous and crystalline phases allows writing and erasing paths for SPhP guided waves as well LSPhP resonator structures ,  (see also Figure 27C). In addition, SPhPs on SiC interacting with deep-subwavelength metal holes form a high-Q cavity with ultra-small modal volume. A thin GST layer sandwiched between the SiC layer and the perforated gold structure can switch the cavity on and off by modifying GST’s dielectric constant by switching it between the crystalline and amorphous phases .
Finally, there has been recent interest in IR-integrated photonics analogous to integrated optics . Long-range IR waveguides of ultra-confined modes would be key components to advance such area. Natural candidates would be of course SPhP-based or GPhP-based structures. Hybridization of PhP-waveguides with other guiding components enables flexibility and control over their guiding characteristics including propagation range. Different types of such hybrid waveguides have been reported such as 4H-SiC SPhP-ridge Si hybrid waveguide  as well as SPhP-insulator-SPhP waveguide with hBN . Moreover, highly tailorable localized IR photonic modes have been reported in  by enabling hybridization between the modes supported by a Ge grating structure and the SPhP modes supported on the interface of SiC. The PhP-quantum well system may also enable new integrated IR sources as recent research shows the large lateral momenta of any type of polaritons facilitate indirect transitions in semiconductors . Furthermore, Ohtani et al.  very recently have reported an electrically pumped phonon-polariton laser. Specifically, the intersubband gain in a quantum cascade structure leads to emission into PhPs supported by one of the constituents in the superlattice (AlAs). This structure can also be regarded as a type of phonon laser as significant energy is directed into the vibrational motion of the atomic ions.
We have recapped above the most predominant research directions involving hybrid platforms with PhP media. Of course many other hybrid systems of variable degree of complexity bringing together different system types for photonic control may be envisioned. For example,  reports on a PhP slab interfaced with a dielectric grating and a graphene metamaterial aiming to achieve tunable THz deflection and asymmetric transmission.
10 Conclusions and outlook
Phonon-polaritonics is a burgeoning field bringing into the IR domain many of the powerful capabilities that until recently had been only available in the optical domain. We reviewed here key widely-used and emergent PhP material including 2D and quasi-2D material that play a central role in the development of this field. We discussed a microscopic picture explaining how these PhP material obtain their extraordinary responses, recapped the vast range of distinct photonic responses they make available as well as presented avenues by which these can be tuned/reconfigured. We reviewed the different types of modes that can be supported with PhP material and discussed the capabilities they bring out and associated application areas they advance. Finally, we presented many examples that move forward the already versatile potential of PhP material by considering metamaterial and mesophotonic structuring routes as well as interfacing them with other photonic modes or systems, such as SPPs, antenna resonators, ridge waveguides or phase-change material.
Phonon-polaritonics has grown immensely in the last years impacting the advancement of a broad range of high-interest application areas in IR photonics, such as thermal imaging , , thermal photovoltaics , SEIRA , sensors and detectors , heat management and radiation cooling , , , , THz imaging for medical diagnostics  as well as components for integrated IR photonics . The parallel rapid development of the fields of atomically-thin material photonics , , PCM photonics , ,  and thermal management , , , , , ,  serves as a catalyst for the continued progress of the phonon-polaritonics field; PhP photonics is expected to continue to flourish in the years to come bringing out many more exciting new phenomena and applications for IR light.
Phonon-polaritonics is a part of the broader field of polaritonics ,  that involves different types of material excitations that strongly couple with light, such as magnon-polaritons (e.g see ) or exciton-polaritons (e.g. see ). As other polaritons also induce a Lorentzian-resonant response in the material’s constitutive properties, for example, in the permittivity for exciton-polaritons  or in the permeability for magnon-polaritons  the progress in phonon-polaritonics can inspire new photonic platforms with other type of polaritonic materials in the applicable spectal ranges (e.g. see  for exciton polaritons).
We would like to thank the UNM Center for Advanced Research Computing, supported in part by the National Science Foundation, for providing the high-performance computing resources used in this work. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. Department of Energy or the United States Government. D. W. acknowledges the support from NSF grant ECCS-1609912. S. K. acknowledges support from the George R. K. Smith Chair account.
Maier SA. Plasmonics: fundamentals and applications. New York, Springer, 2007. Google Scholar
Law S, Podolskiy V, Wasserman D. Towards nano-scale photonics with micro-scale photons: the opportunities and challenges of mid-infrared plasmonics. Nanophotonics 2013;2:103–30. Google Scholar
Miller DL. Principles of infrared technology: a practical guide to the state of the art. New York, Springer, 1994. Google Scholar
Diakides NA, Bronzino JD. Medical Infrared Imaging. Boca Raton, FL, CRC Press (Taylor & Francis), 2007. Google Scholar
Button KJ. Infrared and millimeter waves Vol. 6: systems and components. New York, Academic Press, 1982. Google Scholar
Wanke MC, Lee M, Nordquist CD, et al. Integrated chip-scale THz technology. Proceedings of the SPIE 8031, Micro- and Nanotechnology Sensors, Systems, and Applications III. 2011:80310E. Google Scholar
Joshua Caldwell D. Mid-IR to THz polaritonics: realizing novel materials for nanophotonics (Conference Presentation). Proceedings of the SPIE 9918, Metamaterials, Metadevices, and Metasystems 2016;2016:991828. Google Scholar
Kittel C. Introduction to solid state physics. Hoboken, NJ, John Wiley & Sons, 2005. Google Scholar
Adachi S. The reststrahlen region. In: Optical properties of crystalline and amorphous semiconductors. New York, Springer, 1999:33–62. Google Scholar
Economou EN. The physics of solids: essentials and beyond. Berlin, Heidelberg, Springer-Verlag, 2010. Google Scholar
Brown FC. The physics of solids: ionic crystals, lattice vibrations and imperfections. New York, W.A. Benjamin Inc., 1967. Google Scholar
Ashcroft NW, David Mermin N. Solid state physics. New York, Holt, Rinehart and Winston, 1976. Google Scholar
Palik ED. Handbook of optical constants of solids, Vols. I and II, III. San Diego, CA, USA; Chestnut Hill, MA, USA; London, UK, Academic Press, 1985, 1991, 1998. Google Scholar
Panofsky WKH, Phillips M. Classical electricity and magnetism. Reading, MA, USA; London, UK, Addison-Wesley Publishing Company Inc., 1962. Google Scholar
Raether H. Surface plasmons on smooth and rough surfaces and on gratings. Berlin, Heidelberg, Springer-Verlag, 1988. Google Scholar
Economou EN. Surface plasmons in thin films. Phys Rev 1969;182, 539. Google Scholar
Stratton JA. Electromagnetic theory. New York, London, McGraw-Hill, 1941. Google Scholar
Low T, Avouris P. Graphene plasmonics for terahertz to mid-infrared applications. ACS Nano 2014;2:1086–101. Google Scholar
Khamh H, Sachet E, Kelly K, Mariab J-P, Franzen S. As good as gold and better: conducting metal oxide materials for mid-infrared plasmonic applications. J Mater Chem C 2018;6: 8326–42. CrossrefGoogle Scholar
Barho FB, Gonzalez-Posada F, Milla M-J, et al. Highly doped semiconductor plasmonic nanoantenna arrays for polarization selective broadband surface-enhanced infrared absorption spectroscopy of vanillin. Nanophotonics 2018;7:507–16. Google Scholar
Zhong Y, Malagari SD, Hamilton T, Wasserman DM. Review of mid-infrared plasmonic materials. J Nanophotonics 2015;9:093971. Google Scholar
Poulet H, Mathieu JP. Vibration spectra and symmetry of crystals. New York, Gordon and Breach, 1976. Google Scholar
Yu P, Cardona M. Fundamentals of semiconductors: physics and material properties, Vol. 3, Berlin; London, Springer, 2010. Google Scholar
Jackson JD. Classical electrodynamics. New York, John Wiley and Sons, 1962. Google Scholar
Huang K. Phenomenological equations of motion for simple ionic lattices. Tech Rep 1950;T339:50–9. Google Scholar
Rivera N, Coulter J, Christensen T, Narang P. Ab initio calculation of phonon polaritons in silicon carbide and boron nitride. 2018:Arxiv: 189.00058. Google Scholar
Houstoun RA. A treatise on light. London; New York, Longmans Green and Co, 1942. Google Scholar
Pahuja OP. Solid state physics. New Delhi, Laxmi Publications, 2005. Google Scholar
Dragoman D, Dragoman M. Optical characterization of solids. Berlin, Heidelberg, New York, Springer-Verlag, 2002. Google Scholar
Srivastava GP. The physics of phonons. New York, Taylor and Francis, 1990. Google Scholar
Bittner B, Scherm M, Schoed T, Tyroller T, Schwarz UT, Maier M. Phonon-polariton damping by low-frequency excitations in lithium tantalate investigated by spontaneous and stimulated Raman scattering. J Phys: Condens Matter 2002;14:9013–28. Google Scholar
Tropf WJ, Thomas ME, Harris TJ. Properties of crystals and glasses. In: Bass M, Van Stryland EW, Williams DR, Wolfe WL, editors. Handbook of optics, vol. 2: devices, measurements and properties. New York, NY, McGraw-Hill, 1995. Google Scholar
Martin BG, Broerman JG. Theory of surface polaritons in a polar zero-gap semiconductor. Phys Rev B 1981;24. Google Scholar
Grynberg M, Le Toullec R, Balkanski M. Dielectric function in HgTe between 8 and 300 K. Phys Rev B 1974;9. Google Scholar
Claus R, Merten L, Brandmuller J. Light scattering by phonon-polaritons. Springer-Verlag Berlin, 1975. Google Scholar
Kong JA. Electromagnetic wave theory. New York, John Wiley & Sons, 1990. Google Scholar
Tiwald TE, Woollam JA, Zollner S, et al. Carrier concentration and lattice absorption in bulk and epitaxial silicon carbide determined using infrared ellipsometry. Phys Rev B 1999;60:11464. CrossrefGoogle Scholar
Le Gall J, Olivier M, Greffet J-J. Experimental and theoretical study of reflection and coherent thermal emission by a SiC grating supporting a surface-phonon-polariton. Phys Rev B 1997;55:10105–14. CrossrefGoogle Scholar
Devarapu GCR, Foteinopoulou S. Compact photonic-crystal superabsorbers from strongly absorbing media. J Appl Phys 2013;114:033504-1–11. Google Scholar
Hong W-K, Cha S, Sohn JI, Kim JM. Metal-insulator phase transition in quasi-one-dimensional VO2 structures. J Nanomater 2015;2015:538954. Google Scholar
Ambrosio A, Tamagnone M, Chaudhary K, et al. Selective excitation and imaging of ultraslow phonon-polaritons in thin hexagonal boron nitride crystals. Light Sci Appl 2018;7:27. PubMedCrossrefGoogle Scholar
Ambrosio A, Tamagnone M, Chaudhary K, et al. imaging of ultra-confined phonon-polaritons in hexagonal boron nitride on gold. In: Conference on Lasers and Electro-Optics, OSA Technical Digest (online) (Washington, DC, Optical Society of America, 2018), paper FTh1K.6. Google Scholar
Gilburd L, Kim KS, Ho K, et al. Hexagonal boron nitride self-launches hyperbolic phonon-polaritons. J Phys Chem Lett 8;2017:2158–62. Google Scholar
Richard Powell C. Symmetry, group theory and the physical properties of crystals. New York, Springer, 2010. Google Scholar
Novoselov KS, Mishchenko A, Carvalho A, Castro AH, Neto, 2D materials and van der Waals heterostructures. Science 2016;353:6298. Google Scholar
Foteinopoulou S, Panoiu NC, Shalaev VM, Subramania GS. Feature issue introduction: Beyond thin films: photonics with ultrathin and atomically thin materials. Opt Mater Express 2019;9:2427–36. CrossrefGoogle Scholar
Cappelluti E, Benfatto L, Kuzmenko AB. Infrared phonon activity and Fano interference in multilayer graphenes. Phys Scr 2014;T62:014018. Google Scholar
Caldwell JD, Vurgaftman I, Tischler JG, Glembocki OJ, Owrutsky JC, Reinecke TL. Atomic-scale photonic hybrids for mid-infrared and terahertz nanophotonics. Nat Nanotech 2016;11:9. CrossrefGoogle Scholar
Qian J, Allen MJ, Yang Y, Dutta M, Stroscio MA. Quantized long-wavelength optical phonon modes in graphene nanoribbon in the elastic continuum model. Superl Micros 2009;46:881–8. CrossrefGoogle Scholar
Warner JH, Schaffel F, Bachmatiuk A, Rummeli AH. Graphene: fundamentals and emergent applications. Elsevier, 2013. Google Scholar
Sun Y-Y, Zhang S. Communication: Effect of accidental mode degeneracy on Raman intensity in 2D materials: hybrid functional study of bilayer phosphorene. J Chem Phys 2016;145:21102. CrossrefGoogle Scholar
Krayer LJ, Kim J, Munday JN. Near-perfect absorption throughout the visible using ultra-thin metal films on index-near-zero substrates [Invited]. Opt Mater Express 2019;9:330–8. CrossrefGoogle Scholar
DeVault CT, Zenin VA, Pors A, et al. Suppression of near-field coupling in plasmonic antennas on epsilon-near-zero substrates. Optica 2018;5:2334–6. Google Scholar
da Silva RR, Silva RM, Dumelow T, da Costa JAP, Honorato SB, Ayala AP. Using phonon resonances as a route to all-angle negative refraction in the far-infrared region: the case of crystal quartz. Phys Rev Lett 2010;105:163903. CrossrefPubMedGoogle Scholar
Macedo R, Rodrigues da Silva R, Dumelow T, da Costa JAP. MgF2 as a material exhibiting all-angle negative refraction and sub-wavelength imaging due to the phonon response in the far infrared. Opt Commun 2014;310:94–9. CrossrefGoogle Scholar
Yang T, Jing H, Liu D. An improved description of Jones vectors of the electric fields of incident and refracted rays in a birefringent plate. J Opt A: Pure Appl Opt 2006;8:295–9. CrossrefGoogle Scholar
Caldwell JD, Lindsay L, Giannini V, et al. Low-loss, infrared and terahertz nanophotonics using surface phonon-polaritons. Nanophotonics 2015;4:44–68. Google Scholar
Schuller E, Falge HJ, Borstel G. Dispersion curves of surface phonon-polaritons with backbending. Phys Lett 1975;54A:317. Google Scholar
Ghashami M, Geng H, Kim T, Iacopino N, Cho SK, Park K. Precision measurement of phonon-polaritonic near-field energy transfer between macroscale planar structures under large thermal gradients. Phys Rev Lett 2018;120:175901. PubMedCrossrefGoogle Scholar
Passler NC, Razdolski I, Gewinner S, Schollkopf W, Wolf M, Paarmann A. Second-harmonic generation from critically coupled surface phonon polaritons. ACS Photonics 2017;4: 1048–53. CrossrefGoogle Scholar
Griffiths PR, Surface-enhanced infrared absorption spectroscopy: principles and applications. In: Yarwood J, Douthwaite R, Duckett S, editors. Spectroscopic properties of inorganic and organometallic compounds: techniques, materials and applications. vol. 44. London, UK, RCS, 2013. Google Scholar
Zhang Q, Zhang Y, Li J, Soref R, Gu T, Hu J. Broadband nonvolatile photonic switching based on optical phase change materials: beyond the classical figure-of-merit. Opt Lett 2018;43:94–7. PubMedCrossrefGoogle Scholar
Peters DW, Goldflam MD, Campione S, et al. Resonant ultrathin infrared detectors enabling high quantum efficiency, 2018 IEEE Research and Applications of Photonics In Defense Conference (RAPID), Miramar Beach, FL, USA, IEEE, 2018:1–3. Google Scholar
Baraniuk RG, Kelly KF, Krishna S, Bridge RF. Compressive sensing architecture advances infrared camera design. Laser Focus World 2011;47:31. Google Scholar
Autore M, Li P, Dolado I, et al. Boron nitride nanoresonators for phonon-enhanced molecular vibrational spectroscopy at the strong coupling limit. Light: Sci Appl 2018;7:17172. PubMedCrossrefGoogle Scholar
Woessner A, Parret R, Davydovskaya D, et al. Electrical detection of hyperbolic phonon-polaritons in heterostructures of graphene and boron nitride. npj 2D Mater Appl 2017;1:25. CrossrefGoogle Scholar
Kalkbrenner T, Ramstein M, Mlynek J, Sandoghdar V. A single gold particle as a probe for apertureless scanning near-field optical microscopy. Microscopy J 2000;202:72–6. Google Scholar
Premaratne M, Stockman MI. Theory and technology of SPASERs. Adv Opt Photonics 2017;9:81. Google Scholar
Bohren CF, Huffman DR. Absorption and scattering of light by small particles. John Wiley and Sons, 1983. Google Scholar
Antonoyiannakis EI. Electromagnetic fields and forces in nanostructures. Imperial College of Science. Technology and Medicine, Ph.D. Thesis, 1998. Google Scholar
Ameen M, Garcia-Etxarri A, Schnell M, Hillenbrand R, Aizpurua J. Infrared phononic nanoantennas: localized surface phonon polaritons in SiC disks. Chinese Sci Bull 2010;55:2625–8. CrossrefGoogle Scholar
Razdolski I, Chen Y, Giles AJ, et al. Resonant enhancement of second-harmonic generation in the mid-infrared using localized surface phonon polaritons in subdiffractional nanostructures. Nano Lett 2016;16:6954–9. PubMedCrossrefGoogle Scholar
Dominguez O, McGinnity TL, Roeder RK, Hoffman AJ. Mid- and far-infrared optical characterization of monoclinic HfO2 nanoparticles and evidence of localized surface phonon polaritons. Proc SPIE 2017;10100:101001G. Google Scholar
Kivshar Y, Miroshnichenko A. Meta-optics with Mie resonances. Opt Photonics News 2017;28:25–31. Google Scholar
Rockstuhl C, Salt MG, Herzig HP. Analysis of the phonon-polariton response of silicon carbide microparticles and nanoparticles by use of the boundary element method. J Opt Soc Am B 2005;22:481. CrossrefGoogle Scholar
Caldwell JD, Glembocki OJ, Francescato Y, et al. Low-loss, extreme subdiffraction photon confinement via silicon carbide localized surface phonon polariton resonators. Nano Lett 2013;13:3690–7. CrossrefPubMedGoogle Scholar
Venkataram PS, Hermann J, Tkatchenko A, Rodriguez AW. Phonon-polariton mediated thermal radiation and heat transfer among molecules and macroscopic bodies: nonlocal electromagnetic response at mesoscopic scales. Phys Rev Lett 2018;121:045901. PubMedCrossrefGoogle Scholar
Narimanov EE. Photonic hypercrystals. Phys Rev X 2014;4:041014. Google Scholar
Foteinopoulou S. Photonic-crystal-based polaritonic metamaterials functional at terahertz. In: 14th International Conference on Transparent Optical Networks (ICTON), We.A6.3, Coventry, UK, IEEE, 2012:1–4. Google Scholar
Schuller JA, Zia R, Taubner T, Brongersma ML. Dielectric metamaterials based on electric and magnetic resonances of silicon carbide particles. Phys Rev Lett 2007;99:107401. CrossrefPubMedGoogle Scholar
Guo T, Zhu L, Chen P-Y, Argyropoulos C. Tunable terahertz amplification based on photoexcited active graphene hyperbolic metamaterials [Invited]. Opt Mater Express 2018;8:3941–52. CrossrefGoogle Scholar
Rybin MV, Filonov DS, Samusev KB, Belov PA, Kivshar YS, Limonov MF. Phase diagram for the transition from photonic crystals to dielectric metamaterials. Nat Commun 2015;6:10102. CrossrefPubMedGoogle Scholar
Yannopapas V, Moroz A. Negative refractive index metamaterials from inherently non-magnetic materials for deep infrared to terahertz frequency ranges. J Phys: Cond Matt 2005;17:3717–34. Google Scholar
Eaton HAC. Infrared imaging bolometers. In: Huber MCE, Pauluhn A, Culhane JL, Timothy JG, Wilhelm K, Zehnder A, editors. Observing photons in space: a guide to experimental space astronomy. 2013;9:515–24. Google Scholar
Korobkin D, Urzhumov YA, Neuner III B, et al. Mid-infrared metamaterial based on perforated SiC membrane: engineering optical response using surface phonon polaritons. Appl Phys A 2007;88:605–9. CrossrefGoogle Scholar
Zhang X, Liu H, Zhong Y. Microscopic analysis of surface Bloch modes on periodically perforated metallic surfaces and their relation to extraordinary optical transmission. Phys Rev B 2014;89:195431. CrossrefGoogle Scholar
Sakakibara R, Stelmakh V, Chan WR, et al. Practical emitters for thermophotovoltaics: a review. J Photonics Energy 2019;9:032713. Google Scholar
Svensson T, Hallberg T. Infrared absorption bands measured with an uncooled interferometric LWIR hyperspectral camera. Proc SPIE 10644, 2018:106440Y. Google Scholar
Whatmore R. Ferroelectric materials. In: Kasap S, Capper P, editors. New York, Springer handbook of electronic and photonic materials, 2017:589–614. Google Scholar
Pertsev NA, Zembilgotov AG, Hoffmann S, Waser R, Tagantsev AK. Ferroelectric thin films grown on tensile substrates: renormalization of the Curie-Weiss law and apparent absence of ferroelectricity. J Appl Phys 1999;85:1698. CrossrefGoogle Scholar
Wang B, Woo CH, Zheng Y. Tunable ferroelectric phase transition, UTAM Symposium on Size Effects on Material and Structural Behavior at Micron- and Nano-Scales. Hong Kong, China, Proceedings of the IUTAM Symposium, 2004:1–12. Google Scholar
Guan S, Liu C, Lu Y, Yao Y, Yang SA. Tunable ferroelectricity and anisotropic electric transport in monolayer β-GeSe. Phys Rev B 2018;97:114104. Google Scholar
Kukreti A, Kumar A, Naithani UC. Condensed matter: electronic structure, electrical, magnetic and optical properties-electric field dependence of Curie temperature in BaxSr1-xTiO3 ferroelectric perovskites. Indian J Pure Appl Phys 2011;49:126–31. Google Scholar
Skoromets V, Nĕmec H, Kadlec C, Fattakhova-Rohlfing D, Kužel P. Electric-field-tunable defect mode in one-dimensional photonic crystal operating in the terahertz range. Appl Phys Lett 2013;102:241106. CrossrefGoogle Scholar
Cardona M. Faraday rotation in semiconductors. In: Advances in Solid State Physics. New York, Springer, 2007:72–88. Google Scholar
Chochol J, Postava K, Cada M, Pistora J. Experimental demonstration of magnetoplasmon polariton at InSb(InAs)/dielectric interface for terahertz sensor application. Sci Rep 2017;7:13117. CrossrefPubMedGoogle Scholar
Moncada-Villa E, Fernandez-Hurtado V, Garcia-Vidal FJ, Garcia-Martin A, Cuevas JC. Magnetic field control of near-field radiative heat transfer and the realization of highly tunable hyperbolic thermal emitters. Phys Rev B 2015;92:125418. CrossrefGoogle Scholar
Abraham Ekeroth RM, Ben-Abdallah P, Cuevas JC, Garcia-Martin A. Anisotropic thermal magnetoresistance for an active control of radiative heat transfer. ACS Photonics 2018;5:705–10. CrossrefGoogle Scholar
Bragaglia V, Holldack K, Boschker JE, et al. Far-infrared and Raman spectroscopy investigation of phonon modes in amorphous and crystalline epitaxial GeTe-Sb2Te3 alloys. Sci Rep 2016;6:28560. PubMedCrossrefGoogle Scholar
Kats MA, Blachard R, Zhang SY, et al. Vanadium dioxide as a natural disordered metamaterial: perfect thermal emission and large broadband negative differential thermal emittance. Phys Rev X 2013;3:041004. Google Scholar
Luxmoore IJ, Gan CH, Qiang Liu P, et al. Strong coupling in the far-infrared between graphene plasmons and the surface optical phonons of silicon dioxide. ACS Photonics 2014;1:1151–5. CrossrefGoogle Scholar
Fischetti MV, Neumayer DA, Cartier EA. Effective electron mobility in Si inversion layers in metal-oxide-semiconductor systems with a high-κ insulator: the role of remote phonon scattering. J Appl Phys 2001;90:4587. CrossrefGoogle Scholar
Jiang Y, Lin X, Low T, Zhang B, Chen H. Group-velocity-controlled and gate-tunable directional excitation of polaritons in graphene-boron nitride heterostructures. Laser Photonics Rev 2018;12:180049. Google Scholar
Lina X, Yangc Y, Rivera N, et al. All-angle negative refraction of highly squeezed plasmon and phonon polaritons in graphene-boron nitride heterostructures. Proc Natl Acad Sci USA 2017;114:6717–21. Google Scholar
Dubrovkin AM, Qiang B, Krishnamoorthy HNS, Zheludev NI, Wang QJ. Ultra-confined surface phonon polaritons in molecular layers of van der Waals dielectrics. Nat Commun 2018;9:1762. CrossrefGoogle Scholar
Hajian H, Serebryannikov AE, Ghobadi A, et al. Tailoring far-infrared surface plasmon polaritons of a single-layer graphene using plasmon-phonon hybridization in graphene-LiF heterostructures. Sci Rep 2018;8:13209. CrossrefPubMedGoogle Scholar
Brouillet J, Papadakis GT, Atwater HA. Experimental demonstration of tunable graphene hyperbolic metamaterial. 2019:arXiv:1906.10663. Google Scholar
Hajian H, Ghobadi A, Serebryannikov AE, et al. VO2-hBN-graphene-based bi-functional metamaterial for mid-infrared bi-tunable asymmetric transmission and nearly perfect resonant absorption. J Opt Soc Am B 2019;36:1607. CrossrefGoogle Scholar
Sumikura H, Wang T, Li P, et al. Highly confined and switchable mid-infrared surface phonon polariton resonances of planar circular cavities with a phase change material. Nano Lett 2019;19:2549–54. PubMedCrossrefGoogle Scholar
Finch MF, Saunders CAB, Premkumar N, Yang YC, Lail BA. A 4H-SiC phonon-polariton enhanced hybrid waveguide. Fajardo, Puerto Rico, IEEE Antennas and Propagation Society International Symposium, 2016:987–8. Google Scholar
Qiang B, Dubrovkin AM, Krishnamoorthy HNS, et al. High Q-factor controllable phononic modes in hybrid phononic-dielectric structures. Adv Photonics 2019;1:026001. Google Scholar
Serebryannikov AE, Hajian H, Beruete M, Ozbay E, Vandenbosch G. Tunable deflection and asymmetric transmission of THz waves using a thin slab of graphene-dielectric metamaterial, with and without ENZ components. Opt Mater Express 2018;8:3887. CrossrefGoogle Scholar
De Silva LMS, Wijewardena Gamalath KAIL. Modelling of exciton-polaritons. World Sci News 2018;106:194–213. Google Scholar
We use the term super-absorption here in line with literature in the area of enhancing absorption with structured materials; for example, see . The term signifies a behavior where a structured platform exhibits a strong absorption, typically more than 80%, which at the same time represents a significant absorption enhancement with respect to the absorption exhibited by the lossy constituent material, when in bulk form.
The skin depth, represents the distance the EM wave travels before its electric field strength attenuates by a factor of 1/e, ; so δ/λ0=1/2πκ, where κ is the imaginary part of the complex refractive index, or otherwise known as extinction coefficient .
This is calculated from the well-known exact Raether expression  in its full complex form with the use of realistic parameters for Au as in Figure 1.
The term mode has been used loosely here for the highly-lossy metallic systems to describe the near-field excitation  and associated extinction that corresponds to a specific frequency, ω; in the strictest definition modes can be associated only with lossless systems.
The (−) sign before the iωγ term is consistent with harmonic fields varying in time as e−iωt. This is the physics convention that will be followed throughout the paper. Note, the engineering convention considers harmonic fields varying in time as eiωt, in which case then a (+) instead of (−) should precede the term iωγ in the related expressions.
The quantum mechanical model captures also non-local permittivity responses dependent on the phonon wave vector q; these however are unimportant for typical photonic set-ups where EM waves interact with the near-zero q phonons.
The free carrier screening length is of the order of several lattice constants and also the electrons move much faster than the more inert ions .
The oscillator strength, defined in Eq. (27), is connected to the static permittivity, εs [see Eq. (24)]. Hence, higher oscillator strength values could imply a higher charge Q and thus a higher dipole moment of the ion pair [see Eq. (20)]
These show a decaying rather than an evanescent wave profile from the interface. Purely bound EM modes with an evanescent wave profile are only an idealized scenario at interfaces between materials of opposite permittivity and no dissipation .
We do not discuss Zenneck waves  here which may be supported for the spectral range with Re(ε)>0; Zenneck waves are also quasi-bound surface waves showing much looser confinement though.
β is essentially the modal refractive index .
s-SNOM is a popular technique for mapping the near-field landscapes in the IR and has been widely used in probing PhP-based confined modes. See related review article of .
Representative values for the carrier concentration and electron effective mass as well as the PhP model parameters have been taken from  (see also Table 2 for the latter).
VO2 undergoes a Mott-type insulator-to-metal transition ; around the transitional temperatures there is coexistence of an insulating phase and a strongly-correlated electron metallic phase that is different from the rutile metallic phase at higher temperatures.
The real (imaginary) parts of the permittivity [green-dotted (blue-dotted) lines of Figure 26B-III] has been calculated from  by taking the experimental properties as depicted in Figure 4 and using relations (1), (C1), and (C2) of this reference.
SiO2, like quartz, possesses multiple PhP oscillations; see .
About the article
Published Online: 2019-10-17
Citation Information: Nanophotonics, Volume 8, Issue 12, Pages 2129–2175, ISSN (Online) 2192-8614, DOI: https://doi.org/10.1515/nanoph-2019-0232.
© 2019 Stavroula Foteinopoulou et al., published by De Gruyter, Berlin/Boston. This work is licensed under the Creative Commons Attribution 4.0 Public License. BY 4.0