Subdiffractional focusing and guiding of polaritonic rays in a natural hyperbolic material

Subdiffractional focusing and guiding of polaritonic rays in a natural hyperbolic material

Abstract

concepts↗content_tags↗content_tags↗Uniaxial materials whose axial and tangential permittivities have opposite signs are referred to as indefinite or hyperbolic media. In such materials, light propagation is unusual leading to novel and often non-intuitive optical phenomena. object_of_studymethod↗sensesHere we report infrared nano-imaging experiments demonstrating that crystals of hexagonal boron nitride, a natural mid-infrared hyperbolic material, can act as a ‘hyper-focusing lens' and as a multi-mode waveguide. phenomenoncontent_tagsThe lensing is manifested by subdiffractional focusing of phonon–polaritons launched by metallic disks underneath the hexagonal boron nitride crystal. The waveguiding is revealed through the modal analysis of the periodic patterns observed around such launchers and near the sample edges. Our work opens new opportunities for anisotropic layered insulators in infrared nanophotonics complementing and potentially surpassing concurrent artificial hyperbolic materials with lower losses and higher optical localization.

One of the primary goals of nanophotonics is concentration of light on scales shorter than the free-space wavelength λ . According to the general principles of Fourier optics, this is only possible provided electromagnetic modes of large tangential momenta k t > ω /(2 π ), normally evanescent, are nonetheless able to reach the focal plane (the x – y plane). Here ω=λ −1 is the measure of frequency common in spectroscopy and . In devices known as superlenses 1 2 3 4 5 6 , this requirement is realized via resonant tunneling between the opposite sides of the structure. However, the tunneling is very sensitive to damping, for example, the magnitude of the imaginary part of the permittivity ɛ of the superlens material 7 . The largest characteristic momentum that can pass through a superlens of thickness d can be found from the relation . concepts↗concept_edgesIn this regard, hyperbolic media (HM) 8 9 promise a significant advantage as they support large- k hyperbolic polaritons that remain propagating rather than evanescent, so that the condition on damping is much softer (see below). The unusual properties of hyperbolic polaritons in HM 8 9 10 11 12 13 14 15 16 17 18 19 20 stem from the dispersion of these modes that is described by the equation of a hyperboloid:

where ɛ z and ɛ t ≡ ɛ x = ɛ y are the axial and tangential permittivities, respectively. The hyperboloid is single-sheeted if ɛ z >0, ɛ t <0 (type II) and two-sheeted if ɛ z <0, ɛ t >0 (type I), see Fig 1a,b , respectively. In both cases, the slope of the propagation (group velocity) direction, which is orthogonal to the dispersion surface, asymptotically approaches

The condition for achieving super-resolution is Im k z d =( k t d )Im tan θ ∼1. Hence, admissible Im ɛ z , Im ɛ t scale algebraically rather than exponentially with the resolution .

Directional propagation of hyperbolic polaritons along ‘resonance cones' of apex angle θ has been observed in a magnetized plasma 21 22 , which behaves as a natural HM in the microwave domain. A major resurgence of interest to HM was prompted by their discussion in the context of artificial materials (metamaterials) 23 24 . Examples of such hyperbolic metamaterials include microstrip arrays, where directional propagation and focusing of hyperbolic polaritons have been experimentally observed 25 26 . Directional optical beams have been studied in planar 25 26 27 28 and curved 12 29 metamaterials made of alternating layers of metals and semiconductors. The work on non-planar structures 12 29 was motivated by theoretical proposals of a hyperlens 30 31 32 , a device in which directional beams outgoing from a subdiffractional source enable optical magnification. research_questionHowever, improvement over the diffraction limit has so far been severely impeded by losses in constituent metals and imperfections of nanofabrication.

concepts↗concepts↗concept_edgesRecent work 33 34 identified hexagonal boron nitride (hBN) as a low-loss natural HM in the mid-infrared domain. This layered insulator has emerged as a premier substrate or a spacer for van der Waals heterostructures 35 36 . concepts↗concepts↗concept_edgesLight atomic masses, strong anisotropy and the polar band between B and N yield prominent optical phonon modes that create two widely separated stop-bands—spectral intervals where one of the principal values of the dielectric tensor is negative 33 34 37 . The upper band comprises ω =1,370–1,610 cm −1 where the real part of ɛ t (the in-plane permittivity) is negative while that of ɛ z is positive. In the lower band spanning ω =746–819 cm −1 , the signs of the permittivity components are reversed. concepts↗concepts↗concept_edgesThus, the out-of-plane crystal vibrations enable the type I hyperbolic response, whereas the in-plane ones accounts for the type II behaviour. The momentum-frequency dispersion surface for the hyperbolic polaritons of the upper band resembles a ‘butterfly' ( Fig. 1c ) composed of individual hyperbolas sketched in Fig. 1a . It can be contrasted with the flat dispersion surfaces of longitudinal phonons typical for isotropic materials. Effectively, in hBN the longitudinal phonons are hybridized with the transverse ones by quasi-static Coulomb interaction mediated by large- k photons 38 . Because the hyperbolic response in hBN originates from the anisotropic phonons, in the following, the large- k hyperbolic polaritons are referred to as hyperbolic phonon polaritons (HP 2 ).

Results

Subdiffractional focusing and imaging through hBN

In our experiments, efficient excitation and detection of HP 2 in hBN are accomplished with the help of optical antenna structures 39 40 . The antennas concentrate electric field and bridge the large momentum mismatch between the free-space photons and the HP 2 . In our previous work 33 , we used for this purpose a sharp tip of an atomic force microscope (AFM) incorporated in our scattering-type scanning near-field optical microscopy (s-SNOM) apparatus (Methods). Here, we additionally demonstrate the antenna and polariton-launching capabilities of Au disks patterned on a SiO 2 substrate. The AFM topography image in Fig. 2a depicts Au disks of diameters (top to bottom) 1,000, 500 and 200 nm and thickness of about 50 nm. After the subsequent deposition of hBN crystals of thickness d =100–1,060 nm and lateral sizes up to 10 μm, these Au disks become encapsulated between hBN and SiO 2 . The hBN crystal remains essentially flat, as verified by AFM. Below we present experimental results demonstrating that interaction of these disks with an incident infrared beam excites polaritons that travel across hBN and produce specific contrast patterns at the other surface. We show that the observed dependence of the near-field images on the frequency and hBN thickness is the result of directional propagation of the polaritons along conical surfaces with frequency-tunable apex angle given by Equation (2) . Thus, hBN may emerge as a new standard bearer for mid-infrared nanophotonics by enabling devices for deeply subdiffractional propagation, focusing and imaging with tunable characteristics.

Representative s-SNOM imaging data are shown in Fig. 2 . Figure 2b depicts an s-SNOM scan taken at the top surface of hBN of thickness d =395 nm at frequency ω =1,515 cm −1 ( λ =6.6 μm). Here we plot the third harmonics of the scattering amplitude s ( ω ) (Methods). In this image, each Au disk is surrounded by a series of concentric ‘hot rings' of strongly enhanced nano-infrared contrast. The diameters of all the disks are much smaller than λ (see also Fig. 2a ), the smallest one being 200 nm= λ /33. The diameters of the hot rings can be larger, smaller or equal to those of the disks. The spacing between adjacent hot rings in the same sample increases with the infrared frequency but decreases with the sample thickness. We stress that images displayed in Fig. 2b could only be detected if the infrared wavelength falls inside the hyperbolic spectral regions. Outside of the hBN stop bands, no hot rings can be identified by the s-SNOM. In fact, the entire image is homogeneous, comprised of nothing but random noise, as illustrated by Fig. 2d for ω =1,740 cm −1 ( λ =5.7 μm).

We now elaborate on the formation of images in Fig. 2 recorded with our s-SNOM apparatus with the help of a model of HP 2 propagation through a slab of hBN ( Fig. 3a,c,d ). Consider a perfectly thin metallic disk sandwiched between a slab of a HM of thickness d and an isotropic dielectric substrate. The system is subject to a uniform electric field of frequency ω and amplitude E 0 in the x direction. An approximate solution for the total field in this system can be found analytically ( Supplementary Note 1 ). The corresponding distributions of the z -component of the field E z ( x , y , z ) in the two cross-sections, y =0 (the vertical symmetry plane) and z = d −0 (just below the top surface of the hBN slab), are illustrated in Fig. 3c . These plots are computed for three representative radii of the disk using permittivity values at ω =1,515 cm −1 . The plots demonstrate a series of concentric high-intensity rings on the top surface, very similar to the data in Fig. 2b . The interpretation ( Fig. 3a,c ) is straightforward: the external field polarizes the disk, which perturbs the adjacent HM (hBN in our case) and launches polaritons. The HP 2 emission occurs predominantly at disk edges due to the high concentration of electric field therein. Polaritonic rays propagate across the slab, maintaining a fixed angle θ with respect to the z axis: the ‘resonance cone' direction 18 21 22 25 26 27 28 . Upon reaching the other slab surface, they undergo a total internal reflection with the reflected cone extending toward the bottom surface. The process repeats until eventually the field vanishes because of radial spreading and/or damping. The role of the s-SNOM tip in imaging experiments in Fig. 2 is to out-couple HP 2 fields at the top surface ( Fig. 3a ). The observed s-SNOM signal is roughly proportional to the amplitude of the electric field immediately above the slab E z ( z = d =0). (Note that it is related to the field just inside the slab by a constant factor, E z ( z = d +0)= ɛ z ( ω ) E z ( z = d −0).)

The above model of image formations via HP 2 yields a number of quantitative predictions that are in accord with our observations. The scenario of oblique propagation implies that upon each roundtrip across the slab, the excitation front returns to the same surface displaced radially by the distance

Accordingly, the radii of the ‘hot rings' at the top surface of the slab are given by

where a is the disk radius. The intensity of the rings is expected to decrease with | n |. Consistent with this formula, the smallest rings in Fig. 3c have the radius r 0 =| a −|δ|/2|. Particularly interesting is the case where the innermost ring shrinks to a single bright spot, r 0 =0. Experimentally, we observed spots of diameter 200 nm (the full width at half maximum, see Supplementary Note 1 ), which corresponds to λ /33 for Fig. 2b (top). Focal spots of similar size 185–210 nm were observed in all other hBN crystals, with the thickness up to 1,050 nm ( Supplementary Fig. 5 ).

A proposal for focusing of electromagnetic radiation via resonance-cone propagation in hyperbolic media was theoretically discussed in the context of magneto-plasmas 21 . Experimental confirmation of this idea in an artificial hyperbolic multi-layer was reported where λ /6 focusing was deduced from examining the pattern of a polymerized photoresist behind a two-slit polaritonic launcher 26 . Here, using a natural hyperbolic slab (hBN crystal), we demonstrated the λ /33 focusing in both spatial directions via out-coupling of polaritons with the infrared nano-probe. We stress that a distinction should be made between ‘focusing' and ‘imaging.' Focusing devices can be of both imaging and non-imaging type 41 and both are important in applications. Our hBN device ( Fig. 3a ) is an example of the latter.

Continuing with the verifiable predictions of our model, we note that Equations (3) and (4) indicate that the slope tan θ of the resonance cone is uniquely related to the radii of the hot rings ( Fig. 3a ). To test this prediction, we analysed the images collected from samples of different hBN thicknesses and different Au disk diameters. For each of these, we determined the radius r 1 of the first-order ring and computed |tan θ |=( r 1 − a )/ d as a function of the infrared frequency ( Fig. 3a ). As shown in Fig. 3b , all the data collapse toward a single smooth curve computed from Equation (2) using optical constants of hBN from ref. 33 . Yet another prediction of the model: the polaritonic rays travel along the z axis provided that ɛ t ( ω ) and therefore θ ( ω ) are vanishingly small. This condition is satisfied at ω =1,610 cm −1 ( Fig. 2c ) where we observe almost 1:1 images of Au disks. Similar behaviour was observed when instead of the disks more complicated metallic shapes were imaged ( Supplementary Fig. 5 ). Thus, the totality of our data establishes the notion of directional propagation of HP 2 in hBN over macroscopic distances with a frequency-tunable slope ( Fig. 3b ).

Real-space imaging of multiple guided polaritons in hBN

The outlined real-space picture has a counterpart in its conjugate momentum space. Mathematically, the resonance cones in the real space are coherent superpositions of an infinite number of polariton modes of a slab. Such modes are characterized by quantized momenta, k z , l =( π / d )( l + α ), labelled by integer index l 33 . Here α ∼1 (in general, ω -dependent) quantifies the phase shift acquired at the total internal reflection from the slab surfaces. Per Equation (1) , the tangential momenta of these modes are also quantized,

In the last step, we have applied Equation (3) . For illustration, the dispersion curves of such guided modes in the upper stop-band of hBN of thickness 105 nm are shown in Fig. 1c , where they are overlaid on the dispersion surface of bulk hBN. The same curves are replotted as ω vs k t in Fig. 4a . In Fig. 4b the dispersion curves of the guided modes of lower stop-band are shown. An intriguing aspect of these curves is that their slope ∂ ω /∂ k t is positive (negative) in the upper (lower) band. This sign difference is a consequence of the opposite direction of the group velocity vector for the type I and type II cases, cf. Fig. 1a,b . Central to the connection between the resonance cones in the real space and the quantized momenta in the k -space is that these momenta form an equidistant sequence of period Δ k t = k t, l +1 − k t, l =2 π / δ . Therefore, if several guided modes are excited simultaneously by a source, their superposition would produce beats with period 2 π /Δ k t in real space. This is precisely the spacing δ between periodic revivals of the ‘hot rings' ( Equation (3) and Fig. 2 ). Thus, the multi-ring images and the existence of higher-order guided modes are complementary manifestations of the same fundamental physics. In our previous work 33 , we reported nano-imaging and nano-spectroscopic study of the lowest-momentum guided mode l =0 in hBN crystals. Below we present new results documenting the first observation of the higher-order (up to three) guided modes in such materials by direct nano-infrared imaging.

To map the dispersion of HP 2 , we utilized hBN crystals on SiO 2 substrate without any intervening metallic disks (Methods). Here the sharp tip of the s-SNOM serves as both the emitter and the detector of the polariton waves on the open surface of the hBN. As the tip is scanned toward the sample edge, distinct variations in the detected scattering amplitude s ( ω ) are observed. Such variations are caused by passing over minima and maxima of the standing waves created by interference of the polaritons launched by the tip and their reflections off the sample edges ( Fig. 5a ). Representative data for the upper stop-band (the type II hyperbolic region) are shown in Fig. 5b–f , where we plot s ( ω ) at various infrared frequencies. Specifically, the image presented in Fig. 5b exhibits oscillations with the period ∼1 μ m extending parallel to the edge of a 31-nm-thick hBN crystal. While these oscillations are similar to those reported previously 33 , a high-resolution scan performed very close to the edge (the olive square) reveals additional oscillations occurring on a considerably shorter scale: down to hundreds of nanometres ( Fig. 5c–e ). Similar results have been obtained using many other samples. For example, Fig. 5f also shows short-scale oscillations near the edges co-existing with longer-range oscillations further away from the edge in the data collected for a thicker hBN crystal ( d =105 nm).

To analyse the harmonic content of the measured s ( ω ) quantitatively, we used the spatial Fourier transform (FT). An example shown in Fig. 5h is the FT of the line trace α from Fig. 5g . The three dominant peaks in the FT are marked with β' (blue), γ' (magenta) and ζ' (olive). These peaks have been deemed statistically significant and their positions k β , k γ and k ζ have been recorded for each of the traces studied. We reasoned that including additional weaker peaks into consideration may be unwarranted at this stage. Indeed, the gross features in the real-space trace α exceeding the noise level of ∼1 a.u. are accounted for by oscillations in the three partial traces β, γ and ζ, which are obtained by the inverse FT of the shaded regions in Fig. 5h .

The remaining step in the analysis is to establish the connection of thus determined momenta k β , k γ and k ζ and the momenta k t ,l of the guided modes, Equation (5) . This requires more care than in prior studies of single-mode waves in two-dimensional materials 33 42 43 44 . The interference patterns near the edge can be created by various combinations of the tip-launched waves (labeled by l ) and edge-reflected waves (labeled by r ). The total momentum of a particular combination is k t, l + k t, r . If the mode index is conserved, l = r , the set of possible periods narrows down to 2 k t, l . This is consistent with our data obtained for several infrared frequencies ( Fig. 4a ), where the symbols indicate k β , k γ and k ζ . These data are in a quantitative agreement with the calculated dispersion curves for the l =0, 1 and 2 polariton-guided waves in the upper stop-band. The analysis of polariton propagation length 33 shows that the loss factor is as low as γ ∼0.03 ( Supplementary Fig. 4 ). Dispersion mapping in the lower band (746–819 cm −1 ) where no monochromatic lasers are available is discussed in Supplementary Fig. 3 . Broad-band lasers used in an independent study by Li et al . have allowed to demonstrate focusing behaviour of hBN in this challenging frequency region 45 .

Discussion

Data presented in Figs 2 , 3 , 4 , 5 demonstrate launching, long-distance waveguiding transport and focusing of electromagnetic energy in thin crystals of hBN. These phenomena are enabled by directional propagation of large-momentum polariton beams in this natural hyperbolic material. contributionclaim_typecontent_tagsThe sharpness of the attained focusing, λ /33 at distances up to λ /6 ( Supplementary Fig. 5 ), in units of the free-space wavelength, surpasses all prior realizations of superlenses and hyperlenses. Remarkably, a simple addition of a circular metallic launcher transforms an hBN crystal into a powerful focusing 19 device! The analysis presented in Supplementary Note 1 ( Supplementary Equation 10 ) limitationsindicates that the size of the focal spot in our system is limited by the finite thickness ∼50 nm of Au disks. By using thinner disks, say 20 nm thick, one should be able to achieve focal spots as small as ∼ λ /10 2 , comparable to the spatial resolution of our nano-infrared apparatus. A fundamental advantage of using natural rather than artificial hyperbolic materials is the magnitude of the upper momentum cutoff. In a natural material such as hBN, this cutoff is ultimately set by interatomic spacing thus immensely enhancing the spatial resolution. content_tagsIn addition, we have shown that hBN can serve as a multi-mode waveguide for polaritons with excellent figure of merit: loss factor as small as γ ∼0.03. These characteristics exceed the benchmarks 46 47 48 of current metal-based plasmonics and metamaterials. mechanisms↗The physics behind this fundamental advantage of phonon polaritons over plasmons in conducting media is in the absence of electronic losses in insulators. Applications of hBN for non-imaging focusing devices 41 , subdiffractional waveguides and nanoresonators 34 readily suggest themselves 45 . Combining such elements together may lead to development of sophisticated nanopolaritonic circuits.

Methods

Experimental setup

The nano-imaging and nano-FTIR experiments described in the main text were performed at UCSD using a commercial s-SNOM ( www.neaspec.com ). The s-SNOM is based on a tapping-mode AFM illuminated by monochromatic quantum cascade lasers (QCLs) ( www.daylightsolutions.com ) and a broad-band laser source utilizing the difference frequency generation ( www.lasnix.com ) 49 . Together, these lasers cover a frequency range of 700–2,300 cm −1 in the mid-infrared. The nanoscale near-field images were registered by pseudo-heterodyne interferometric detection module with AFM tapping frequency and amplitude around 250 kHz and 60 nm, respectively. To obtain the background-free images, the s-SNOM output signal used in this work is the scattering amplitude s ( ω ) demodulated at the n th harmonics of the tapping frequency. We chose n =3 in this work.

Sample fabrication

Silicon wafers with 300-nm-thick SiO 2 top layer were used as substrates for all the samples. The Au patterns of various lateral shapes and 50-nm thickness were fabricated on these wafers by electron beam lithography. The hBN microcrystals of various thicknesses were exfoliated from bulk samples synthesized under high pressure 50 . Such microcrystals were subsequently mechanically transferred onto either patterned or unpatterned parts of the substrates.

Author contributions

All the authors were involved in designing the research, performing the research and writing the paper.

Additional information

How to cite this article : Dai, S. et al . Subdiffractional focusing and guiding of polaritonic rays in a natural hyperbolic material. Nat. Commun . 6:6963 doi: 10.1038/ncomms7963 (2015).

Supplementary Material

Supplementary Figures 1-5, Supplementary Notes 1-2 and Supplementary References.