Hyperbolic metamaterials: fusing artificial structures to natural 2D materials

Optical metamaterials have presented an innovative method of manipulating light. Hyperbolic metamaterials have an extremely high anisotropy with a hyperbolic dispersion relation. They are able to support high-k modes and exhibit a high density of states which produce distinctive properties that have been exploited in various applications, such as super-resolution imaging, negative refraction, and enhanced emission control. Here, state-of-the-art hyperbolic metamaterials are reviewed, starting from the fundamental principles to applications of artificially structured hyperbolic media to suggest ways to fuse natural two-dimensional hyperbolic materials. The review concludes by indicating the current challenges and our vision for future applications of hyperbolic metamaterials.

Among various metamaterials that have been realized to date, hyperbolic metamaterials (HMMs) have attracted great interest due to their highly anisotropic characteristics [21][22][23]. The hyperbolic dispersion of HMMs is determined by the effective permittivity tensor where the principal components of the electric or magnetic fields have opposite signs. The relatively easy fabrication of multilayer or nanowire (NW) structures can yield HMMs that have a three-dimensional (3D) bulk response at optical frequencies. Besides, such materials can support the hyperbolic dispersion, where in principle the infinitely large unbounded momentum and thus the very high confinement of light can be allowed. Thus, such unique anisotropic properties have been applied to various distinguished applications including super-resolution imaging [24,25], negative refraction [26,27], and emission engineering [28,29].
The propagation loss is one of the main limitations of bulk HMMs, due to the inevitable Ohmic loss of incorporated plasmonic materials, so methods to reduce it are being actively sought after. As a result, two-dimensional (2D) hyperbolic metasurfaces (HMSs) have been investigated, rendering the hyperbolic propagation and large confinement of light at the surface and not the bulk, and thus are less sensitive the loss and have shown great potentials as planar optical devices. Besides, some pristine 2D materials and bulk crystals can render the natural hyperbolic response, without the need of structuring, and have been explored as a powerful alternative candidate for the realization of HMMs and HMSs. Natural 2D hyperbolic materials have a homogeneous hyperbolic dispersion, so have been explored as a powerful alternative candidate for the realization of HMMs and HMSs. Numerous materials that exhibit hyperbolic dispersion at various wavelengths have been demonstrated, signalling the high potential of HMMs and HMSs for implementation in practical applications [30][31][32][33].
Herein, we first introduce the basic theory for the realization of HMMs. We will present details of hyperbolic dispersion, bulk HMMs, and planar 2D hyperbolic materials. Following the fundamental discussions, we will review the progress of applications of HMMs and HMSs, including the most recently reported papers and newly emerging fields, with emphasis on high-resolution optical imaging and nanoscale lithography, light propagation and manipulation, spontaneous emission engineering, sensors, and absorbers. Then we present the recent tremendous progress of natural hyperbolic materials that are key for further investigation to realize planar hyperbolic applications. Last, we conclude by highlighting the current challenges and suggesting an overview for future research on hyperbolic metamaterials. While following the existing framework, this review contains the latest publications that have not been reviewed elsewhere. We expect that this review will faithfully fill the gap with existing review articles and provide insight and potential of HMMs as a powerful tool for photonic devices.

Dispersion
The concept of HMMs originates from optical crystals, where the effective permittivity ( ̂ ) and permeability ( ̂ ) tensors of the material oriented along the principal axes are described as where, 0 and 0 are the permittivity and permeability of a vacuum, respectively. In general, materials are isotropic, which means that the permittivity and permeability components are equal in all directions ( xx = yy = zz , xx = yy = zz ). If one or more of these principal components is different, the material becomes either a uniaxial or biaxial anisotropic medium. Assuming that the optical axis (OA) is along the z direction, an electric uniaxial medium has xx = yy ≠ zz , whereas a magnetic uniaxial medium has xx = yy ≠ zz ; and an electric biaxial medium has xx ≠ yy ≠ zz , whereas a magnetic biaxial medium has xx ≠ yy ≠ zz . For simplicity, we now consider electric uniaxial material under transverse magnetic (TM) polarized light. The dispersion relation of the light in the medium with effective permittivity in Eq. (1) can be derived from Maxwell's equations as where k x , k y , k z are the wave vector components along the x, y, and z directions, respectively, is the angular frequency, and c is the speed of light. In Eq. (2), the isofrequency contour (IFC) can be obtained at the intersection of the constant frequency plane and the dispersion surface. In general materials, the components of the effective permittivity all have the same sign, so the IFC is closed and forms a sphere or ellipsoid (Fig. 1a). Due to the bounded IFC, waves with large wave vectors become evanescent in those materials. In contrast, under the condition where the signs of two components of the permittivity are opposite ( zz ⋅ xx < 0 ), the IFC becomes an unbounded hyperboloid that can support high-k waves [21-23, 34, 35]. These properties are not commonly obtained in nature [33,36,37], but can be achieved in a fairly straightforward manner using metamaterials [1,[38][39][40]. HMMs are named after the hyperbolic topology of their IFC, from which their unique properties stem from. HMMs can be classified into two types according to the sign of their permittivity. Materials with xx = yy > 0, zz < 0 have a predominantly dielectric nature, and are therefore called dielectric-type or Type-I HMMs (Fig. 1b); materials with xx = yy < 0, zz > 0 have more metallic properties, and are subsequently known as metallic-type or Type-II HMMs (Fig. 1c). Because of their predominant metallic nature, Type-II HMMs show highly absorptive properties, and hence only support high-k waves that have a cut-off spatial frequency. Magnetic HMMs [41][42][43] can also be achieved under transverse electric (TE) polarization by using a magnetic permeability tensor that has opposite signs in one or two directions.

Bulk hyperbolic metamaterials
A few naturally occurring hyperbolic materials have a strong anisotropy [33,36,37]. Examples include the crystal structures of graphite (C) [30], hexagonal boron nitride (hBN) [31], magnesium diboride ( MgB 2 ) [32], calcite ( CaCO 3 ) [44], and crystals with perovskite layers of cuprate and ruthenate [33]. In these structures, anisotropic crystal structures are arrayed in layers, and the natural hyperbolicity results from the two-dimensional conductivity of the parallel atomic plane. The operating wavelength of these materials is intrinsically bounded by their natural crystal structure and can only be slightly tuned by environmental factors. Various naturally occurring hyperbolic materials are summarized in Table 1 according to their operating wavelengths. However, artificially structured HMMs may be beneficial due to the fact that their optical properties can be easily adjusted by changing the constituent meta-atoms. HMMs have both metallic and dielectric properties depending on the direction, so the most straightforward way to realize HMMs is to design structures composed of metals and dielectrics. So far, two structures have been typically used for bulk HMMs: alternating multilayers of sub-wavelength thickness metal-dielectrics (multilayer HMMs) [49,50], and metallic NWs embedded in a dielectric host medium (NW HMMs) [26,51]. The realization of HMMs using these structures relies on effective medium theory (EMT) [52], which describes the macroscopic optical responses of composite materials using homogenized effective parameters. When the size of the constituent materials is sufficiently small, the homogenization is valid and a hyperbolic effective permittivity tensor can be achieved.
Multilayer HMMs are the simplest structures, in which the motion of free electrons is restricted to directions parallel to the interfaces between layers. In sub-wavelength thickness multilayer HMMs, the effective permittivity components parallel ∥ and perpendicular ⟂ to the OA can be approximated as [34] where m and d are the permittivity of the metal and dielectric materials, respectively, and f = t m t m +t d is the volume fraction of metal and dielectric with thicknesses of t m and t d , respectively. The effective permittivity can be customized by the choice of materials and f. The tuning of these two parameters can yield both Type-I and Type-II multilayer HMMs, as well as effective dielectric or metallic properties. Multilayer HMMs fabricated with evaporation or sputtering deposition have generally been used to create horizontal layers. However, to exhibit Type-I HMMs, this structure requires an optical resonance, which results in a narrow working wavelength and high resistive losses. Recently, vertically-alternating multilayer HMMs have been proposed to circumvent these optical resonances and achieve Type-II HMMs [27,53,54].
The use of NW HMMs is another approach to achieve HMMs, where metallic NWs are embedded in a dielectric host medium. In NW HMMs that have sub-wavelength diameters and lattice periods, the effective  permittivity tensor components parallel ∥ and perpendicular ⟂ to the OA are [34] where m and d are the permittivity of the metal wire and dielectric host, respectively, and f is the volume fraction of the metal. The real part of the metal permittivity is negative at frequencies below the plasma frequency, so the NW structures can easily achieve Type-I HMMs without the requirement of an optical resonance. Consequently, they can operate at relatively low loss and high transmission over a large bandwidth. Due to the longwavelength limit condition, structures should have a high aspect ratio with an extremely small diameter compared to the slab thickness. Those structures with high aspect ratio have been mainly fabricated using electrochemical deposition methods in which metals such as silver (Ag) and gold (Au) [21,55] are grown inside an anode alumina template. The electrochemical deposition method allows for efficient large-area fabrication, but limits the choice of materials that can be used.

2D planar hyperbolic materials
The unique hyperbolic IFC of bulk HMMs have demonstrated the capability to control light in unconventional ways, but they generally have high propagation losses and demanding fabrication requirements. Accordingly, attempts to compress bulk HMMs to 2D surfaces have been made to alleviate this problem. Metasurfaces are 2D planar devices that can manipulate light with the advantages of having a compact form factor and lower propagation losses than 3D structures [17,56,57]. In addition, the simplified fabrication and integration processes of metasurfaces facilitate the compatibility of integration with other photonic devices. Metallic metasurfaces exploit the control of the propagation of surface plasmon polaritons (SPPs) along the metasurfaces. For these ultra-thin metasurfaces, the effective surface conductivity tensors, which are analogous to the effective dielectric tensor, Eq. (5) are used to describe the homogenized medium, as where xx , and yy are the x and y components of the surface conductivity. Metasurfaces can support both TE and TM modes, but the confinement of the TE mode is not as strong as the TM mode. Therefore, we only consider TM modes confined to metasurfaces. The dispersion relation of TM modes in metasurfaces is given as [58,59] (4) where 0 is the impedance of free-space. HMSs are extremely anisotropic metasurfaces that exhibit a hyperbolic dispersion of in-plane SPPs [58,60,61]. That is, the sign of the imaginary parts of the effective surface conductivity components are different to each other. Consequently, HMSs exhibit both metallic (inductive, Im( ) > 0 ) and dielectric (capacitive, Im( ) < 0 ) properties depending on the direction of surface conductivity. Due to the unique dispersion relation of in-plane SPPs, HMSs can support extremely large wave vectors, and have a large photonic density of states, so many interesting phenomena can be achieved [61][62][63]. In practice, HMSs have been realized with either metallic gratings [61] or natural 2D hyperbolic materials such as graphene [58,64,65] and hBN [66,67] at high frequencies.
However, at low frequencies, noble metals are regarded as perfect electric conductors, in which SPPs become weakly confined. Instead, structured metal surfaces can confine surface waves called spoof SPPs [68,69], where the structured surfaces allow the support of cut-off frequency modes and sub-wavelength confinement at low frequencies. Spoof SPPs can mimic the novel properties of SPPs, and the dispersion of spoof SPPs can be simply tailored by adjusting the geometrical parameters of the structured surface. Therefore, HMSs that support the propagation of spoof SPPs in the microwave regime have been demonstrated using structured metal resonator structures [70,71].

High-resolution optical imaging and nanoscale lithography
The diffraction limit prevents conventional microscopes from observing sub-wavelength features [74][75][76]. This limit is a result of the depletion of evanescent waves at the image plane in the far-field. The light scattered from an object is composed of various spatial frequency components. The low-k waves contain information about the large geometric features of the object, and are able to propagate in a vacuum. On the contrary, the dispersion relation in a vacuum ensures that the longitudinal component for wave with transverse wave vector higher than k 0 (high-k wave) necessarily becomes imaginary. Consequently, the waves that encode the sub-wavelength features decay exponentially and cannot reach the image plane (inset (i) in Fig. 2a). The hyperbolic dispersion in HMMs supports the propagation of the high-k waves, therefore showing the potential for super-resolution imaging [21,34]. The unbounded IFC of HMMs allows arbitrarily large transverse wave vectors with real values of longitudinal wave vectors. The first theoretical concepts of superresolution imaging have been demonstrated using several HMMs [77][78][79], but the waves become evanescent again as the waves exit the HMMs (inset (ii) in Fig. 2a). A hyperlens is a device designed to overcome this limitation [24,25,80]. Hyperlenses made of HMMs utilize their geometrically-curved shapes to magnify the subwavelength information as well as supporting high-k waves (inset (iii) in Fig. 2a). When the wave propagates in the radial direction of a hyperlens, the transverse wave vectors are compressed in accordance with the law of conservation of angular momentum; i.e., the image is gradually magnified as it is transferred along the radial direction. The magnification ratio of the hyperlens can be estimated from the ratio of the inner and outer radii. Therefore, evanescent waves can be converted to propagating waves by carefully choosing a radius ratio that sufficiently compresses the transverse wave vectors. The converted waves can then be projected onto the image plane in the far-field, and thereby achieve super-resolution. Hyperlenses can be readily integrated into conventional optical systems [73,81] proving their promise for practical applications.
Super-resolution imaging using a hyperlens was first demonstrated in the ultraviolet (UV) regime [24]. A cylindrical hyperlens composed of alternating concentric multilayers of Ag and Al 2 O 3 ( Fig. 2b) was used to achieve a resolution down to 0 ∕3 , where 0 is the wavelength of light. However, the cylindrical shape restricts its magnification to only one direction, so the geometrical shape was then extended to a spherical shape to achieve full-lateral sub-wavelength imaging (Fig. 2c) [25]. Still, the curved inner surfaces of the cylindrical or spherical hyperlens makes sample placement difficult and hinders their practical applications. As alternatives, planar hyperlenses [72,[82][83][84] were developed (Fig. 2d). Planar hyperlenses can be designed using transformation optics [82], where the light path is deliberately modulated to magnify the subwavelength features.
Recent studies on hyperlenses have sought to address some of the practical limitations [73,85]. Hyperlenses are generally inefficient due to high reflection and absorption losses. In particular, the efficiency decreases as the thickness of the hyperlens is increased, which in turn limits the achievable magnification. One factor responsible for the low efficiency is the huge impedance mismatch at the surface of a hyperlens, so to reduce the high reflection, impedance-matched hyperlenses have been proposed [86]. A non-resonant hyperlens with a radial (fan-shaped) configuration has been demonstrated to reduce the absorption caused by resonances [85]. The radially alternating metal-dielectric multilayer allows the realization of a hyperbolic dispersion in a non-resonant way, and hence low losses in the visible range. The correct selection of materials can also avoid the intrinsic losses of metals. Non-metallic materials of transparent conducting oxides [87] have been used to replace lossy metals. To alleviate the sample-positioning problem, a wafer-scale hyperlens device has been proposed ( Fig. 2e) [73]. A scalable hyperlens array was fabricated on a 4-in. wafer, and real-time imaging of biological systems using the hyperlens array was demonstrated as a practical application.
The characteristic of HMM to support sub-diffractionlimited optical patterns has recently been exploited to improve the resolution of conventional microscopy by combining HMMs with conventional structured illumination microscopy [88][89][90]. HMMs are used to generate sub-diffraction-limited structured light patterns, and the scattering fields are collected in the far-field to improve spatial resolution down to ∼ 40 nm [90]. Sub-diffractional focusing and imaging using natural hyperbolic materials have also been demonstrated to avoid intrinsically high optical losses caused by metals in HMMs [31,67,[91][92][93][94][95]. In contrast to HMMs, the upper limit of the propagating wave vector of natural hyperbolic materials can be greatly extended due to their homogenous nature, and it is only limited by the crystal lattice periodicity [96]. In addition, the long scattering lifetime properties of phonon polaritons supported by natural hyperbolic materials make them attractive candidates for imaging. Therefore, near-field focusing and imaging have been actively demonstrated using natural hyperbolic materials such as hBN [31,67,91,92,97], graphene [94,95], and -MoO 3 [93]. Recently, by exploiting ultra-low loss monoisotopic materials of h 11 BN , hyperlens device with significantly improved spatial resolution down to ∕154 has been demonstrated [97]. The unique properties of HMMs beyond the diffraction limit are also helpful in generating deep sub-wavelength patterns in lithography [98][99][100][101]. The reciprocity theorem suggests that a hyperlens can also be used to reduce diffraction-limited patterns to sub-diffraction patterns. Sub-wavelength photolithography using hyperlenses was therefore demonstrated by simply reversing the direction of the incident light [98,100]. As the light travels along the radial direction to the inner radius, the transverse wave vectors are gradually increased, so the pattern size decreases (Fig. 3a). The first experimental demonstration of a de-magnifying hyperlens using an Ag-TiO 2 multilayered structure [53]. A pattern with line-width of 300 nm inscribed on the outer hyperlens surface was reduced to 170 nm (a de-magnification ratio of 1.8). Recently, Fig. 3 Sub-wavelength patterning using HMMs. a Sub-wavelength light focusing using hyperlenses. b A schematic of photolithography with a hyperlens and c its simulation result for demagnifying 1D line patterns. d A schematic of deep sub-wavelength patterning in lithography. e Atomic force microscopy images of a 1D grating structure and f its interference pattern. a Adapted with permission from ref. [98], copyright 2009, AIP. b, c Adapted with permission from ref. [102], copyright 2018, ACS. d-f Adapted with permission from ref. [99], copyright 2015, WILEY a de-magnification ratio of 3.75 was achieved using a hyperlens with a high radius ratio (Fig. 3b, c) [102].
Other approaches that use HMMs to obtain high-resolution patterns have been used in interference lithography [101,103]. Laser interference lithography is an inexpensive large-area nanofabrication method that utilizes the periodicity of the interference pattern between light waves [104]. Therefore, the resolvable pattern size is again restricted by the diffraction limit, so methods to decrease the operating wavelength are being sought [105,106]. HMMs can support high-k waves, and therefore have the ability to generate sub-wavelength patterns [99,101,107]. Numerical simulations suggested that a conventional chromium grating with 320 nm period can generate a 40 nm-sized sub-wavelength pattern [101]. Subsequent rigorous analysis of the mechanism and properties of the SPPs in HMMs have extended the understanding of the formation of sub-wavelength interference patterns [107]. Sub-wavelength interference lithography using HMMs has been experimentally demonstrated using Al-SiO 2 multilayers (Fig. 3d) [99]. The mutual coupling of SPPs between adjacent metal-dielectric interfaces results in bulk plasmon polariton modes, which generate large-area interference patterns. An Au grating with a period of 360 nm (Fig. 3e) was used to generate uniform interference patterns with a linewidth of 45 nm (Fig. 3f ) over a large area of 20 mm × 20 mm. Recently, the use of Type-II epsilon-near-zero HMMs for high transmission has achieved high fidelity deep subwavelength lithography [108].

Light propagation and manipulation using HMMs
The hyperbolic dispersion of HMMs provides a innovative method of manipulating the propagation of light. Among the various light-propagation phenomena that have been achieved using HMMs, negative refraction has been widely discussed. Various metamaterials such as split-ring resonators [112], multilayered fishnet [113], chiral [40,114,115], and asymmetric nanoring [116] structures have been employed to achieve negative refraction. HMMs lead to negative refraction under TM polarized light which is affected by . Therefore, negative refraction using HMMs can be realized by adjusting even if neither nor have negative values that are needed for a negative n. The principle of negative refraction using HMMs can be described using the IFC (Fig. 4a). In conventional isotropic media such as air, the wave vector k and Poynting vector of TM mode light are parallel. However, if TM mode light propagates along the interface between a HMM and air, k and are no longer parallel. If we consider the x-z plane between a Type-II HMM and air, negative refraction can be observed because the tangential components of k and point in opposite directions [80,117]. This negative refraction implies energy-negative refraction only when light passes through the HMM in an isotropic medium, which differs from the case of an isotropic left-handed metamaterial, in which the phase and energy are negative simultaneously. The direction of propagation and the refracted light is anti-parallel [118,119].
Various investigations of negative refraction have been performed using HMMs and HMSs at different wavelengths [26,61,[109][110][111]. At visible wavelengths, negative refraction has been experimentally demonstrated using HMMs composed of sub-wavelength Ag NW structures (Fig. 4b) [26]. The fabricated bulk Ag NW HMM has opposite signs of along the perpendicular and parallel directions of the structures, producing a hyperbolic dispersion to achieve low-loss wide-angle negative refraction at wavelengths from 660 to 780 nm. Negative refraction of SPPs has also been numerically demonstrated using HMSs [109]. The simulations demonstrated that the SPPs travel from the metal surface to the HMS, and negative refraction occurs at the interface. A simple diffraction-grating structure was used to control the dispersion of the SPPs, and broadband wide-angle negative refraction was demonstrated. The dispersion relation of the grating coupled metasurface represents a hyperbolic dispersion that is flat around = 543 nm (Fig. 4c). At this wavelength, the SPPs propagate without diffraction with a group velocity normal to the IFC, so negative refraction of SPPs is observed (Fig. 4d). Negative refraction from metasurfaces has been demonstrated experimentally [61]. Metasurfaces composed of single-crystalline Ag films showed a hyperbolic dispersion below the critical wavelength of = 540 nm ( Fig. 4e). At > 540 nm, the single-crystalline Ag film has an elliptical dispersion. SPPs that propagate along the metal film meet the interface between metal and singlecrystalline Ag film, and are refracted. The propagation of SPPs changes to negative refraction only at the hyperbolic wavelength region below = 540 nm, this result shows that negative refraction in HMSs is dependent on the wavelength of the incident light (Fig. 4f ).
In the near-infrared region, an experimental demonstration of negative refraction using HMM with alternating multilayers of zinc oxide (ZnO) and aluminum-doped zinc oxide (AZO) has been reported [110]. Hyperbolic dispersion was achieved at 1.84 μm ≤ ≤ 2.4 μm (Fig. 4g). Measurements of the blocking of propagating light when negative refraction occurred show a clear dip in the transmittance spectrum to prove the existance of negative refraction (Fig. 4h). In the far-infrared region, wide-angle negative refraction has been achieved using a multilayer of thick semiconductor layers of doped indium gallium arsenide (InGaAs) and intrinsic AllnAs that has low-loss properties ( Fig. 4i) [111]. Negative refraction in the terahertz/gigahertz regions has also been achieved using natural hyperbolic materials, which will be introduced in more detail in Sect. 3.5. All-angle highly-squeezed polariton negative refraction was achieved using HMSs composed of anisotropic or nanostructured 2D materials over a broadband wavelength region. The working bandwidth can also be varied by tuning the chemical potential of the materials [120]. Moreover, experiments have demonstrated HMSs that have an anisotropic magnetic sheet conductivity supporting in-plane highk designer polaritons. This phenomenon overcomes the difficulty of practical realization despite various theoretical proposals of negative in-plane refraction [121].
Negative refraction through multilayered bulk HMMs has been typically achieved using horizontally-stacked multilayered HMMs (hHMMs). A broad range of working wavelengths △ w and low-loss are desirable characteristics, so hHMMs are limited by their use of narrows resonances △ w and high losses. To alleviate these drawbacks, recently, HMMs composed of  [27], copyright 2019, NPG. b Adapted with permission from ref. [26], copyright 2008, AAAS. c, d Adapted with permission from ref. [109], copyright 2013, AIP. e, f Adapted with permission from ref. [61], copyright 2015, NPG. g, h Adapted with permission from ref. [110], copyright 2012, NAS. i Adapted with permission from ref. [111], copyright 2007, NPG vertically-stacked multilayered HMMs (vHMMs) have been demonstrated to realize broadband negative refraction [27]. The calculated effective permittivity of aluminum (Al) and SiO 2 multilayer structures show broader △ w in vHMMs than in hHMMs in the visible range (Fig. 5a). The negative refraction of vHMMs is not dependent on resonances, so they can achieve wide △ w and avoid losses that originate from resonances [124]. Therefore, the transmission in vHMMs can be higher than in hHMMs. Numerical demonstrations using the finite-element method verified the clear negative refraction from vHMMs with higher transmittance than hHMMs (Fig. 5b). A numerical study of the critical layer thickness of vHMMs in the visible region has been reported [125]. A variety of vHMMs composed of various materials have been numerically investigated to determine the critical layer thickness of the potential for negative refraction utilization. Recently, a vHMM composed of Au and copolymer resist (EL8) has been fabricated through two-step focused ion beam milling. This vHMM experimentally demonstrates broadband negative refraction between 450 nm and 550 nm in the visible region [126].
Phase change materials and liquid crystals have also been applied to HMMs to achieve active control of the negative refraction. One such switchable HMM used a chalcogenide glass alloy of Ge 2 Sb 2 Te 5 [122], which undergoes a phase change from amorphous to crystalline. By using a Ge 2 Sb 2 Te 5 nanocomposite as one of the layers in a multilayer structure HMM, the phase change yields a switchable Type-I hyperbolic dispersion and subsequently the negative refraction changes from the visible to near-infrared region (Fig. 5c).
Nanosphere-dispersed liquid crystals (NDLCs) have also been applied to achieve switchable HMMs in the near-infrared region [123]. Appropriate application of an external electric field was shown to produce a switching mechanism between normal refraction, negative refraction, and reflection of TM mode light (Fig. 5d). Negative refraction using 3D assembled nanorod arrays [127] and using 2D materials such as graphene, boron nitride, and organic materials has also been investigated, the latter of which will be introduced in detail in the following section [120,[128][129][130].
In addition to negative refraction, the spin Hall effect of light (SHEL) [133][134][135][136][137] is another interesting application of HMMs. SHEL refers to a transverse and spindependent displacement of light. In the conventional spin Hall effect, spin-up and spin-down electrons are deflected in opposite directions, whereas in the SHEL, refracted or reflected light at an optical interface undergoes spin-dependent splitting (Fig. 6a)  for horizontally polarized incidence and as V = cot i Re(1 − t p ∕t s )∕2 for vertically polarized incidence [135]. The shift is generally much smaller than the wavelength, so effort have been devoted to increasing the shift using metamaterials and metasurfaces [132,138,139]. A straightforward way to increase the shift is to design a medium that has polarization-dependent transmission or reflection characteristics. Hyperbolic metamaterials offer a natural candidate to increase the shift from the SHEL.
Amplification of the SHEL using HMMs has been proposed theoretically using planar [140,141] and curved multilayers [142], while experimental demonstration has been achieved in the visible range using a horizontallystacked multilayer [131]. A large incident-angle-sensitive SHEL shift that reached 165 μm was also achieved experimentally (Fig. 6b). The SHEL may be further amplified by several orders of magnitude through the use of vertically-stacked HMMs instead of horizontally-stacked HMMs [54]. vHMMs have a different OA configuration than horizontal HMMs (Fig. 6c), which results in a large discrepancy between the transmission of s-and p-polarized light (Fig. 6d). A numerical demonstration showed an increase in the shift of several orders of magnitude, and an experimental demonstration using a vHMM with an Au nanograting verified a micrometer scale increase in transverse shift (Fig. 6e).
The analytic formula of the shift implies that the shift diverges as the Fresnel coefficient in the denominator approaches zero. However, the efficiency of the SHEL is defined as the intensity of the spin-separated beam divided by the intensity of the incident light [132] which decreases quickly due to the relationship with the square of the Fresnel coefficient in the denominator. Thus, many previous attempts to increase the SHEL have been marred with low efficiency. A HMM composed of an array of metallic wires in a dielectric host has been suggested to achieve a large SHEL with near-unity efficiency (Fig. 6f ) [132]. Because of the in-plane anisotropy, the wire-based HMM structure shows near-unity transmission for one polarization but extremely low transmission for the other (Fig. 6g). The numerically and experimentally-obtained shifts and efficiencies agreed well and proved a large shift that simultaneously reached 1.57 wavelengths and near 100 % efficiency in the microwave regime.

Spontaneous emission engineering
After the discovery that spontaneous emission can be increased using a resonant cavity, research in photonic and plasmonic fields has verified the increased spontaneous emission rate through the engineering of the local density of states (LDOS) near cavities or nanostructures [143]. HMMs have shown the potential as an alternative to resonant cavities. The resonance in cavities hinders the development of broadband spontaneous emission and highly-confined electromagnetic fields. In contrast, the HMMs can be a promising solution for enhancing the spontaneous emission over the broad bandwidth. From Fermi's golden rule, the radiative spontaneous emission rate can be defined as [144,145] where is a matrix that expresses the transition dipole moment of energy states, (r 0 , 0 ) is the LDOS, r 0 is the position of the emitter, and ℏ 0 is the photon energy. Equation (7) demonstrates that modification of the electromagnetic environment around the emitters, especially by an adjustment of the LDOS, can increase the spontaneous emission rate involving any factors related to the nature of emitters such as transition dipole moment, chemical structure, and energy levels in the transition. In this respect, the HMMs that support the open IFCs and thus diverging LDOS substantially enhance the spontaneous emission rate. Thus, research on HMMs has improved the manipulation of spontaneous emission which can be characterized from the measured emission spectrum and lifetime of the emitters. HMMs have been studied the characteristics of light emitted from a dye or dye-containing layer on a NW or metal-dielectric multilayer-based structure. The effects of HMMs on the increase of spontaneous emission have been studied theoretically and experimentally [28,29,149,[151][152][153][154][155][156]. The enhanced spontaneous emission rate near the HMMs is characterized by the Purcell factor, defined as the total emitted power of the source near the HMMs by that at the air. Furthermore, to maximize the effect of spontaneous emission due to the confinement of higher wave vectors using HMMs, research has been actively conducted to examine the effect of placing the emitter directly inside the HMM [147,157,158]. Although many researchers have used multilayer HMMs to manipulate the spontaneous emission of dyes or emitters, the focus has only been on observing the Purcell factor of emitters located near HMMs. Most of these studies could not overcome the limitation of multilayered thin-film-type HMMs that have non-radiative properties. Trials to increase the radiative emission efficiency from HMMs are required to enhance the spontaneous emission and yield high enough Purcell effect that can be detected even in the far-field. Therefore, outcoupling techniques have been applied to HMMs to effectively control and improve the LDOS and radiative emission [146-148, 157, 159-168]. Recently, photoluminescence has been increased by using HMMs composed of films of deposited perovskites [169] and HMMs composed of coupled perovskite nanocrystals [170]. The advantages of HMMs can be explained by comparing the Purcell factors of single metal layers and HMMs (Fig. 7a, b) [146]. The Purcell factors are calculated for the UV to visible range for each structure, with an emitter 10 nm away from the single Ag film and Ag-Si multilayer. A bulk single Ag layer has a high Purcell factor within a narrow bandwidth at a surface plasmon resonance in the UV region (Fig. 7a). On the other hand, increase the number of alternating layers results in interface plasmon polaritons that have shifted resonance frequency. These multiple bands of plasmon polaritons broaden the bandwidth of the enhanced Purcell factor [34]. Thus, the HMM composed of a Ag-Si multilayer has a high Purcell factor over a broad range of visible frequencies (Fig. 7b). By adjusting the filling ratio of the metal in HMMs, the position of the peak of the Purcell factor can be freely controlled within the target wavelength range (Fig. 7c). In addition, by exploiting the outcoupling effect of a grating or nanopatterns, radiative emission can be increased compared to flat multilayer HMMs. The Purcell factor and radiative enhancement at the emission wavelength can be observed through the comparison of uniform HMMs and nanopatterned multilayer HMMs fabricated by focused ion beam milling. Both uniform and nanopatterned HMMs have a high Purcell factor (Fig. 7d), but the calculated radiative emission results show remarkable differences (Fig. 7e). The radiative enhancement on the nanopatterned HMM at the peak wavelength is almost double; this change increases the external quantum efficiency of the emitters. HMMs can also be outcoupled using bullseye gratings [157], nanoantennas (Fig. 7f) [147], 2D grating patterns, and rolled-up structures (Fig. 7g) [148] for the enhancement of the directional emission of dipoles.
Emission can also be manipulated by controlling the topological transitions in HMMs [149,171]. A transition in the topology of the IFC can be induced by a shift of the wavelength, and as for the HMM composed of Ag-TiO 2 , the IFC may change from closed (elliptical) to open (hyperbolic dispersion) at a specific wavelength (Fig. 7h).
Phase-change materials (PCMs) have been used to develop tunable HMMs to control the spontaneous emission. A topological transition from an elliptical to a hyperbolic dispersion can be achieved for tunable HMMs by incorporating a PCM. The LDOS near a HMM can change according to the phase of the PCM, so the emission is also affected. A tunable HMM composed of titanium nitride (TiN) and the PCM, antimony sulfide ( Sb 2 S 3 ), demonstrated a considerable tunability of the Purcell factor and emission efficiency over the phase change of Sb 2 S 3 from crystalline to amorphous [150]. Another PCM, vanadium dioxide ( VO 2 ), has been i Effective permittivity of vanadium dioxide ( VO 2 ) in metal and semiconductor phases, and the Purcell factor and quantum efficiency in both phases. a-e Adapted with permission from ref. [146], copyright 2014, NPG. f Adapted with permission from ref. [147], copyright 2020, WILEY. g Adapted with permission from ref. [148], copyright 2016, APS. h Adapted with permission from ref. [149], copyright 2012, AAAS. i Adapted with permission from ref. [150], copyright 2020, OSA demonstrated in a tunable HMM that was used to control the LDOS at visible and NIR wavelengths (Fig. 7i).
Besides the optical applications, hyperbolic dispersion in HMMs can be utilized in thermal emission applications [35,[172][173][174]. At thermal equilibrium, the black body emission in the far-field follows Plank's law. However, the high DOS of HMM can achieve super-Plankian thermal emission in the near-field that can overcome the aforementioned limit [172,175]. Although conventional photonic crystals or surface electromagnetic excitations have been reported, the narrow bandwidth has always been the bottleneck. On the other hand, non-resonant media with hyperbolic dispersion have provided thermal emission enhancement in a broad bandwidth [176][177][178], and various applications of emitters and radiative thermal diode have been reported [179][180][181][182]. More applications are expected using artificial HMMs and 2D materials.

HMM-based sensors and absorbers
HMMs have been investigated for the use as ultrasensitive optical sensors and broadband absorbers [23,127,[183][184][185]. The hyperbolic dispersion of HMMs supports large wave vectors that have a large momentum and can propagate inside or on the surface of HMMs. The large wave vector components in HMMs are non-radiative, so a grating or prism can be used as a coupler [50,77]. The bulk plasmon polaritons in HMMs can be excited by the coupler, and the excited mode can stimulate multimode resonances that have high quality factors and can be utilized in ultrasensitive optical sensors [77].
Multilayer or nanorod structures of HMMs have been applied as refractive-index sensors for the detection of biomolecules. Multilayer HMMs composed of Au and Al 2 O 3 with a 2D patterned grating have been used as a biosensor to sense wavelength and angular modulations with high sensitivity in the visible and NIR regions [186]. The measured sensitivity was 30,000 nm/RIU (refractive index unit), which is strong enough to detect small molecules at concentrations less than picomoles. A HMM sensor based on porous nanorods with high sensitivity to refractive-index changes near the HMM medium has also been demonstrated [185]. In addition, HMM sensors have also used chiral structures and PCMs [187,188]. A numerical investigation of a hyperbolic chiral-structure thin film as an angular-modulation sensor showed a sensitivity [187]. A tunable HMM as a reconfigurable plasmonic sensor has been proposed by structuring a hybrid HMM platform with TiN or graphene and Sb 2 S 3 as the PCM [188]. The Sb 2 S 3 -TiN HMMs sensors have demonstrated significant sensitivity to detect extremely low concentrations of molecules below 1 × 10 −15 M.
Broadband absorbers using HMMs have also been proposed. The broadband high absorption of incident light has been achieved using HMM waveguide structures whereby the large attenuation is attributable to slowlight modes [189][190][191]. Typical sawtooth and pyramidal waveguide HMM structures have been extensively studied [192][193][194][195]. In numerical [193] and experimental [195] demonstrations, tapered waveguide HMM structures have demonstrated outstanding absorption of > 95% in the mid-infrared (MIR) region from 3 to 5.5 μm, as well as in the visible and NIR regions. Recently, non-tapered HMM broadband absorbers with a uniform height have been reported [196,197]. A multilayer HMM made up of a Au-zinc sulfide (ZnS) grating structure has been used to demonstrate wide-angle and broadband absorption in the MIR region as a result of multiple slow-light modes. A non-tapered HMM broadband absorber in the MIR region has been demonstrated and applied to radiative cooling [197] to radiate heat through the atmosphere without consuming energy [198][199][200][201].

Natural hyperbolic metamaterials
The homogenized anisotropic responses of metamaterials are obtained by the engineering of artificial structures, such as dielectric and metallic multilayers [49] and periodic nanorod arrays [26]. However, these approaches have some limitations. The first is the homogenization model of metamaterials, i.e., the EMT as discussed in Sects. 2.2 and 2.3. Such simple descriptions overlook nonlocal effects induced by the intrinsic near-field coupling between two neighboring structures, especially close to resonant frequencies. Bragg scattering in inhomogeneous multilayers limit the maximum value of possible wave vectors to the boundary of its first Brillion zones (i.e., /Λ , where Λ is the periodicity of unit cells). In that sense, Λ is required to be as small as possible, which is not practical considering the limited resolution of fabrication capabilities. Moreover, HMMs with plasmonic media have unavoidable Ohmic losses which makes them intrinsically lossy due to electron-electron scattering. To overcome these problems, recent developments in materials science have revealed many natural low-dimensional nanomaterials as alternative candidates to HMM platforms for various applications.
The first class of optical hyperbolicity consists of hyperbolic plasmon polaritons (HPPs) in 2D materials. Polaritons are quasiparticles of the collective excitation of strongly coupled light and matter [202]; plasmon polaritons are due to the coupling between the light and free or quasi-free electrons. Many natural materials can support HPPs, including black phosphorus (BP) [203][204][205][206][207], the semimetal tungsten ditelluride ( WTe 2 ) and others [208,209], which present reduced symmetry (orthorhombic crystal structure for BP and WTe 2 ) in their atomic orientations. When free carriers are injected, the electron motions can differ along the different characteristic axes, i.e., zigzag direction (x axis in the inset of Fig. 8c), and armchair direction (y axis in the inset of Fig. 8c), which induces electromagnetic anisotropy. BP has different band gaps in the two orthogonal directions, so the interband transitions and intraband transitions are correspondingly excited by the light, where the inductive [Im( ) > 0] and capacitive [Im( ) < 0] responses are induced (Fig. 8a), and cause the hyperbolic propagation of plasmons [205][206][207] in the infrared and terahertz ranges (Fig. 8b). Specifically, the surface conductivity can be written as where = x, y, , m , , and s are the relaxation time, effective mass of the carriers, onset frequency of interband transitions, and strength along directions, respectively. Θ( − ) is a step function. The behaviors along the x and y directions differ (Fig. 8a-c). For some particular doping levels, an emerging hyperbolic frequency range with different signs of conductivity along the two directions occur (Fig. 8a), with supported hyperbolic b Adapted with permission from ref. [206], copyright 2016, APS. c Adapted with permission from ref. [205], copyright 2014, APS. e, f Adapted with permission from ref. [210], copyright 2014, AAAS. h, i Adapted with permission from ref. [47], copyright 2018, NPG. k, l Adapted with permission from ref. [211], copyright 2021, NPG in-plane propagation (Fig. 8b). This natural hyperbolicity at the surface is an ideal platform for HMSs. More recently, similar HPPs have been reported and observed in anisotropic Weyl semimetals of WTe 2 [209], while other materials are yet to be explored.
The second class of natural optical hyperbolicity is composed of hyperbolic phonon polaritons (HPhPs) in polar van der Waals (vdW) materials [212,213]. vdW nanomaterials are made of 2D atomic layers. Different interaction mechanisms and consequent interaction strength exist in atomically thin materials, including relatively weak vdW forces between the constituent atomic sheets, and strong Coulomb interactions of in-plane covalent bonds. Even for covalent bonds, the interaction strengths can also be dramatically different, as it depends on the participating atoms and their separation distances. Phonon polaritons (PhPs) are coupled oscillations between the light and ion lattice vibrations in materials that can be polarized (polar dielectrics); [214] rich anisotropic PhP responses can be found in layered vdW materials. The first example is uniaxial hBN, which has a so-called out-of-plane hyperbolicity xx = yy = ⊥ , zz = ∥ , ⊥ ≠ ∥ , Re( ⊥ ) ⋅ Re( ∥ ) < 0 within the Reststrahlen bands. Here, the Reststrahlen bands represent the frequency range between the longitudinal and transverse optical phonons of the materials which, for hBN, contains both Type I ( ∥ < 0 and ⊥ > 0, 746-819 cm −1 ) and Type II ( ∥ > 0 and ⊥ < 0, 1370-1610 cm −1 ) hyperbolicity (Fig. 8d). For a thin hBN slab, HPhPs propagate as in a waveguide mode, with the dispersion given by the constructive interference conditions of highly confined TM modes, which are written as [210] where d is the thickness of the hBN and sup ( sub ) are the dielectric constants of the superstrate and substrate, the integer l is the mode index, tan ( ) = i √ ⊥ √ ∥ , and k = √ k 2 x + k 2 y is the in-plane wave vector. Using recently developed real-space nanoimaging techniques (scattering-type scanning near-field optical microscopy), one can easily map out the spatial distribution of HPhPs (Fig. 8e) and thus their dispersion features (Fig. 8f ). The HPhPs in hBN capture the most of exotic properties of HMMs, such as large field confinement and directional propagation [31,92]. Moreover, the maximum possible wave vector supported in hBN is limited by the size of its atomic lattice, which is several orders of magnitude larger than that of HMMs made of artificial structures. Most importantly, fundamentally distinct from plasmonic systems, HPhPs are, in principle, immune to Ohmic losses due to (9) the absence of electron-electron scattering. Therefore, they provide a promising resource on which to study many applications of HMMs that require low-loss. However, HPhPs in uniaxial hBN do not carry anisotropic responses at the interface because the optical axis of exfoliated hBN is always perpendicular to the interface. To mimic the functionalities of hyperbolic responses over a surface with PhPs, some other materials have recently been discovered. One example is alpha-phase molybdenum trioxide ( -MoO 3 ), a biaxial vdW material, which can be endowed with in-plane hyperbolic responses [Re( xx )⋅Re( yy ) < 0] within certain Reststrahlen bands (for example, around 900 cm −1 as shown in Fig. 8g) [46,47,215,216]. Using real-space nanoimaging techniques, one can probe the in-plane HPhPs distribution at the thin -MoO 3 slab (Fig. 8h) and can experimentally retrieve the dispersion (Fig. 8i), which agrees excellently with the analytical dispersion. Here, the analytical dispersion is obtained by treating the thin slab as the surface with the conductivity tensor connected to the permittivity tensor [217]. Such approaches directly indicate that the thin in-plane anisotripic vdW films can be ideal platforms on which to study infrared HMSs, and the further incorporation of many other possible materials including -V 2 O 5 [48] can add more opportunities to achieve the on-demand manipulation of light in natural hyperbolic materials.
Those HPPs and HPhPs in 2D materials mostly reside in middle-infrared or terahertz frequency ranges, as constrained by either the strength of optical excitations (for plasmon) or the intrinsic lattice vibration resonance (for phonon). For the hyperbolicity in the visible, the exciton polaritons (ExPs) should be visited. Excitons are another elementary material excitations in a form of bound electron-hole pairs by Coulomb interactions, forming the ExPs when coupled to electromagnetic waves. Such ExPs are generic and can exist in most 2D semiconductors, that hold the band gap in the visible and near-infrared, including perovskites, transition metal-dichalcogenides, BP and others. They can also present the extreme anisotropy in the presence of strong atomic lattice asymmetry. For instance, it is revealed that the 2D hybrid perovskites made of alternating layered organic and inorganic sublattices can mimic the artificially metal-dielectric HMMs and thus present the out-of-plane hyperbolic responses [218]. More interestingly, BP can also hold in-plane hyperbolic ExPs (HExPs) [211]. As discussed, BP shows the thickness-dependent band gap, spanning from the visible and near-infrared, and the in-plane distinct lattice arrangements for different principal directions. The excitons therein are hence anisotropic [219] and can be extreme for some specific cases. Figure 8j plots the optical conductivity of monolayer BP in zigzag and armchair directions, showing the HExPs near 1s excitons. Figure 8k shows the loss function, defined by −1∕Im( ( , )) while and represent the wave vector and frequency and = 1 + i ∕ d 0 is the permittivity of BP with thickness d, of monolayer BP along two distinctive directions. One can clearly see the distinctive excitonic resonant loss in zigzag and armchair directions, as the signature of the in-plane hyperbolicity (Fig. 8l). In general, HExPs can host strong light confinement for enhanced lightmatter interactions due to the open IFC and, more importantly, several other unique figures of merits. The first one is the flexible tunability, as the result of sensitive band gaps, taking BP for example, to the conditions such as layer numbers, strain and electrostatic gating. Secondly, excitons can become radiative. The photonic LDOS of those radiative excitons buried in the intrinsic hyperbolic media can be even more enhanced, compared to that of an emitter in the adjacent, promising more efficient energy harvesting and light-emitting applications [218]. Last, more advanced pathway toward the excitonic hyperbolicity can also exist in, for example, Rydberg excitons, which can be programmed by the transient photoexcitations. We note that the HExPs are only recently identified, while more future efforts on its fundamental mechanisms, materials exploration, experimental mappings, optical engineering and applications should definitely be made in the community.
To close this section, we will briefly discuss more possibilities of natural hyperbolic materials, their applications and limitations which may be possible to overcome in combination with HMMs. Firstly, beyond the routes mentioned above, other methods to achieve natural hyperbolicity exist, including photoexcitation of electron-hole pairs in layered transition-metal dichalcogenides for broadband anisotropy [220], and materials with reduced symmetry [221,222], to name just a few. Secondly, since the hyperbolic responses in natural materials have lower loss and promise much larger momentum (meaning higher field confinement), tremendous nanotechnologies are thus possible, such as imaging [31,92,223], hyperlensing [67], biosensing [224,225], waveguide [226], on-chip electro-optic modulator [227], and many others [54,132,228]. Thirdly, we remark that those hyperbolic responses can exist in beyond 2D materials, while the bulk media can also hold the hyperbolicity. For instance, the HPhPs are recently found in calcite, a polar crystal that has been long studied for more than three centuries. The infrared calcite polaritons feature the raylike and ultra-long propagation distance with a titled wavefront even at room-temperature [44]. The uniqueness of 2D materials is their atomic thinness, important for the highly compact and integrated devices, while its challenges of large-scale fabrication is formidable. The bulk media such as calcite can be easily obtained and engineered, offering the readily commercial platform for hyperbolic photonic devices, but may suffer from the large volume and poor compactness. Lastly, those hyperbolicities are intrinsically bounded by pristine forms of material candidates and hence have fixed working frequency ranges, which impose a limit on the working frequency ranges of various exciting applications. Further application of the philosophy of metamaterials, i.e., to mold electromagnetic responses by engineering sub-wavelength structures, could yield architectures of natural hyperbolic materials that induce exotic wave phenomena, such as infrared HMSs in nanostructured hBN [66,67], topological polaritons in twisted ( -MoO 3 ) bilayers [217], stacked multilayers for heterostructures [65,128,[229][230][231][232], and others [233].

Conclusion and outlook
HMMs have extreme anisotropy and can support highk waves; distinctive features that have been exploited to demonstrate various optical applications. In this review, we have focused on the fundamentals of HMMs and HMSs, and reviewed their applications. First, the fundamental background of hyperbolic dispersion in the view of bulk and 2D planar structures was discussed, and then various applications of HMMs were introduced. The wide range of applications of HMMs highlights the potential of functional optical devices and platforms.
The development of hyperbolic media that will lead to the improvement of influential optical devices is more extensive than the contents of this review. For example, photostability, low-loss organic HMMs and super-resolution imaging platform using HMMs are being actively investigated [90,[234][235][236][237][238]. The propagation of light in HMMs can be applied to various optical components for collimation [70,239], splitters [240][241][242][243], airy beams [244,245], and even with absorption and scattering for sensing and cloaking [246]. Furthermore, lasers and hypercrystals that use HMMs are also being actively studied [247][248][249][250], and cavity structures capable of resonance regardless of their size and scale are being investigated [143]. HMMs have also been applied to various physical phenomena such as paritytime symmetry [251], mechanical properties [252], magneto-optical Kerr effect [253], nonlinear responses [254], chiral surface waves [255], and newly-researched substances such as perovskites [256]. Exploration of surface waves at boundaries of HMMs or surface waves with hyperbolic dispersion is also a fascinating area. HMMs support a unique type of directional surface waves with either full or partial hyperbolic dispersion called Dyakonov plasmons [57,61,[257][258][259][260].
Similar surface waves with hyperbolic dispersion also exist between an isotropic and anisotropic media, even when the anisotropic media has an elliptical dispersion [261,262]. Besides, a HMS that supports electric and magnetic surface waves simultaneously has been reported recently [263]. Finally, combination of hyperbolic dispersion with other effective optical properties opens a new path to guide surface waves topologically [264][265][266].
To extend the range of practical applications of HMMs and HMSs, some problems must be overcome. First, the efficiency of HMMs must be increased by minimizing intrinsic Ohmic losses caused by the commonly-used metals [267]. Maximizing the efficiency could be achieved by minimizing the proportion of metal. Active or tunable devices would also be extremely valuable. In addition to the examples introduced above, various materials are expected to be evaluated as candidates for active or switchable devices that have many potential applications. Finally, straightforward and inexpensive fabrication technologies for the realization of large-area hyperbolic media and patterns would also be beneficial [268]. Recently, various designs other than the multilayer and nanowire structures are proposed such as 3D nanorod arrays [127], core-shell nanoparticles [269,270], hyperbolic van der Waals crystals [271], and films of carbon nanotubes [272]. In addition, metasurface studies have included many attempts to achieve the large-area processing of HMMs using nanomanufacturing methods such as nanoimprinting, colloidal lithography, aerosol jets parallel printing, and spark plasma sintering [273][274][275][276]. If successful, these methods could substantially accelerate the use and development of HMMs for practical application in photonic devices.