Objective Two-dimensional (2D) van der Waals (vdW) crystals like graphene and a-MoO3 can support polaritons in the spectral range from terahertz to mid-/far-infrared regime, enabling nanoscale confining, focusing, and controlling of the electromagnetic fields. The hybridization between different polaritons can further enrich the properties of polaritons and bring more degrees of freedom for the regulation of electromagnetic fields at the nanoscale. In this paper, we studied the hybridization of plasmon polaritons and phonon polaritons in a heterostructure composed of an a-MoO3 vdW thin lamina stacking onto a monolayer graphene. An analytical waveguide model was developed to calculate the polariton propagation characteristics in the vdW heterostructure. The dispersion contours, dispersion relations, and localized electromagnetic field distributions of the hybridized polariton waveguide modes were derived. The theoretical results were then verified by real-space optical nano-imaging and numerical simulations. Our study can provide a quantitative model for the calculation of the hybridized polariton waveguide modes in vdW heterostructures, which can help further exploring the interactions between different types of polaritons in 2D vdW crystals.MethodsIn our theoretical model, the vdW heterostructure is treated as a 2D infinite waveguide supported onto SiO2 substrate,which consists of a monolayer graphene of 0.5-nm thickness and an a-MoO3 lamina of 115-nm thickness (Fig.1). The dielectric functions of the graphene and a-MoO3 are described using the Drude model and Lorentz model, respectively. Because in a-MoO3 the polaritons can approximately be treated as transverse magnetic (TM) mode, the electromagnetic modes and the associated dispersion relations of the heterostructure can then be obtained by solving the Maxwell's equations upon the continuities of the electric and magnetic fields at interfaces. Monolayer graphene was grown by the chemical vapor deposition (CVD) method. Microfabrication technique combining electron beam lithography (EBL) and reactive ion etching (RIE) was employed to pattern the graphene into microstructures. The vdW a-MoO3 laminas were grown using a physical vapor deposition method. Dry transfer method was utilized to prepare the a-MoO3/graphene (a-MoO3/Gr) heterostructures, where the monolayer graphene microstructures were covered with a-MoO3 laminas of different thicknesses. Real-space nano-imaging was conducted using a scattering-type scanning near-field optical microscope (NeaSNOM, Neaspec GmbH, Germany). In a specific measurement, a metal-coated tip (Arrow-IrPt, Nanoworld,Switzerland) was illuminated using a mid-infrared laser (Access Laser, USA) with a wavelength range of 9.20-10.70 mu m (934.5-1087.0 cm-1). The tip was vibrated vertically with a frequency of about 280 kHz. The backscattered light from the tip was detected in a pseudo-heterodyne interferometric manner, where the scattered light was demodulated at the fourth harmonic of the tip vibration frequency. The optical and morphological images of the sample can be simultaneously obtained by scanning the heterostructure underneath the tip. For the numerical study, the real-space polariton waves were manifested as the real-part of the z-component of the electric field, Re(Ez), on the surface of the SiO2 substrate. They were calculated using the finite element method (FEM) simulations (COMSOL Multiphysics). A vertically-polarized electric dipole source was fixed above the a-MoO3 with a separation of 50 nm (Fig.S5, Supporting materials). The thicknesses of the air,a-MoO3, graphene and SiO2 layers were 500 nm,115 nm,0.5 nm and 500 nm, respectively. The permittivities along the three principle axes were calculated according to Eq.(1). The anisotropic dielectric tensors of the a-MoO3 layer are written as epsilon=[epsilon x,0,0;0,epsilon y,0;0,0,epsilon z] and epsilon=[epsilon xx,epsilon xy,0;epsilon yx,epsilon yy,0;0,0,epsilon zz]. These tensors were imported into the COMSOL package to solve the Maxwell's equations. For the monolayer graphene and a-MoO3lamina, the Re(Ez) was monitored respectively on the planes 5 nm away from their upper surfaces, while for the a-MoO3/Gr heterostructure, the Re(Ez) was monitored on the plane 4 nm away from the graphene upper surface.Results and DiscussionsIn the theoretical model, the vdW heterostructure is modeled as a 2D infinite waveguide (Fig.1). The thicknesses of top (epsilon[1]) and bottom layers (epsilon[2]) are d1 and d2, respectively. It is sandwiched between two semi-infinite plates, which act as the substrate (epsilon s=epsilon[3]) and cover layer (epsilon c=epsilon[0]). The electromagnetic modes [Eqs.(S8)-(S10), Supporting materials] and the associated polariton dispersions [Eq.(4)] are obtained by solving the Maxwell's equations upon the continuities of the electric and magnetic fields at interfaces. With the input of dielectric functions of monolayer graphene [Eq.(2)] and a-MoO3[Eq.(1)], the calculated polariton dispersions and contours of the a-MoO3/Gr heterostructure are shown in Fig.1 and Fig. S4 in the Supporting materials. In comparison with the pristine monolayer graphene and a-MoO3 lamina, due to the isotropic plasmon polariton in the graphene, the dispersion contours of the heterostructure are more complex and distorted along the [100] and [001] directions.Therefore, the hybridized plasmon-phonon polaritons can propagate along these directions that are forbidden respectively for the phonon polaritons in Restrahlen Band 1 and Band 2 in a-MoO3 lamina (Fig.2). The theoretical results are further corroborated by numerical simulations using FEM (Figs.1 and 2) and experimental nano-imaging measurements (Figs.3 and 4). Moreover, the influence of the thickness of a-MoO3 lamina on the polariton hybridizations in a-MoO3/Gr heterostructure is also investigated.Because the polariton fields are of evanescent nature, by reducing the a-MoO3 thickness, the hybridized polaritons converge to the plasmon polariton in monolayer graphene, while with the increase in the a-MoO3 thickness, the polaritons in the heterostructure evolve into the phonon polaritons in the pristine a-MoO3(Fig.4). The results are also corroborated by the nano-imaging measurements (Fig.4).ConclusionsIn conclusion, we have established a theoretical model to investigate the hybridizations of plasmon polaritons and phonon polaritons in a heterostructure consisting of an a-MoO3 lamina covering a monolayer graphene. The propagation characteristics, including the polariton dispersion relation, in-plane dispersion contour, and localized electromagnetic field distribution, were calculated and studied. It is revealed that due to the hybridization effect, the a-MoO3/Gr heterostructure is able to support polariton propagation along the directions that are forbidden for the phonon polaritons in pristine a-MoO3 lamina. Additionally, the influence of the a-MoO3 thickness on the polariton hybridization in the heterostructure was also investigated,indicating that when the a-MoO3 lamina was thinner/thicker, the hybridized polaritons became more plasmon/phonon polariton-like.The theoretical results were corroborated respectively by the numerical simulations and experimental nano-imaging measurements.We strongly believe that the results obtained in our study can on one hand provide a theoretical model for analytically studying the polariton hybridizations in vdW heterostructures, and on the other hand help further our understanding on the polaritonic physics in low-dimensional materials
