Abstract
Van der Waals materials and heterostructures that manifest strongly bound exciton states at room temperature also exhibit emergent physical phenomena and are of great promise for optoelectronic applications. Here, we demonstrate that nanostructured, multilayer transition metal dichalcogenides (TMDCs) by themselves provide an ideal platform for excitation and control of excitonic modes, paving the way to exciton-photonics. Hence, we show that by patterning the TMDCs into nanoresonators, strong dispersion and avoided crossing of exciton, cavity photons and plasmon polaritons with effective separation energy exceeding 410âmeV can be controlled with great precision. We further observe that inherently strong TMDC exciton absorption resonances may be completely suppressed due to excitation of hybrid light-matter states and their interference. Our work paves the way to the next generation of integrated exciton optoelectronic nano-devices and applications in light generation, computing, and sensing.
Similar content being viewed by others
Introduction
Control of light-wave dispersion at the nanoscale is of great importance for applications ranging from lasing to sensing to computing. Plasmonic and high contrast dielectric nanostructures1,2,3,4 are a frontier approach for dispersion engineering (Fig. 1a, b) at the dimensions comparable to and smaller than the light wavelength. Materials with resonant quantum confined states provide an alternative route for controlling light propagation and interaction. This strategy has been effectively employed in a wide variety of materials ranging from inorganic IIIâV epitaxial quantum wells to organic small molecules and polymers, to even carbon nanotubes, and most recently, two-dimensional (2D) transition metal dichalcogenides (TMDCs) (Fig. 1c). Due to weak intermolecular bonding or physical confinement, excitons in this class of materials exhibit strong binding energies and dominate optical responses even at room temperatures5. Among these excitonic materials, the 2D TMDCs of Mo and W are of particular interest since the strong exciton binding manifests in high refractive indices and extinction coefficients (k) near the exciton resonance. Lightâexciton interaction has been extensively studied in monolayer direct-band gap TMDC films6,7,8,9,10,11. These studies have involved coupling the monolayers to extrinsic plasmonic or dielectric optical meta-elements or cavities or metasurfaces12,13,14. Strong coupling regimes and families of hybrid exciton modes have been observed. Nevertheless, light-interaction with monolayer films is challenged by the very large disparity between optical wavelength and film thickness. Recent demonstration of enhanced resonant light interaction with few-layer thick TMDCs suggests that thin systems may be well suited for light dispersion5,15,16,17.
Light interaction with nanostructured metals (a) and dielectrics (b) gives rise to the excitation of plasmonic and dielectric-cavity resonances, respectively. Dispersion of these resonances is determined mainly by the geometry of the system. c Reduced dimensional materials provide a different approach to strong lightâmaterials interaction. Here, sharp optical resonances manifest in the excitation of long-lived bound electron-hole pairs, with layered two-dimensional transition metal dichalcogenides (TMDCs) as prominent examples. d By mixing engineered geometric dispersion in metal-dielectric nanostructures with intrinsic excitonic resonances in TMDCs, a multifaceted, strongly-interacting interplay of photonic and electronic states can be achieved.
Here, we show that nanostructured multilayer TMDCs provide a natural platform for taming lightâexciton interaction. Specifically, within a single integrated platform we observe multipartite interaction and hybridization, leading to strong coupling between three resonant excitations with an avoided crossing of 410âmeV at room temperature. We further reveal a regime of light-TMDC interaction associated with suppression of inherent exciton absorption resonances due to interference of hybridized states at the nanoscale.
Results
Nanopatterned multi-layer WS2 grating resonators on Au
We begin by systematically exploring optical reflectance spectroscopy, first in unpatterned (Fig. 2a) and then in patterned cases (Fig. 2b) of varying TMDC thickness. Our choice of the TMDC is WS2, but our results are generalizable to other room temperature excitonic materials (see Supplementary Note 1). Micromechanically exfoliated WS2 flakes were obtained on a template-stripped Au substrate by mechanical exfoliation from bulk crystals using a tape18. The flakes were directly transferred using the tape onto the template stripped Au substrates. Additional details of sample preparation are provided in the âMethodsâ section. Upon exfoliation, these flakes exhibit thickness-dependent colors resulting from FabryâPerot-like resonances16, resulting in translucent red, purple, blue, and white regions of crystal (Fig. 2c). Periodic one-dimensional (1D) gratings were etched in single WS2 flakes to increase the lightâmaterials interaction strength and engineer optical dispersion. Additional fabrication details and measurement methods are provided in Supplementary Note 2. Structural characterization of grating etch patterns by atomic force microscopy (AFM) and scanning electron microscopy (SEM) (Fig. 2c, d) shows deep-subwavelength WS2 flake thicknesses, along with gratings etched through the full flake. The etching process results in some edge roughness, but these non-perturbative structural deviations do not modify the optically critical parameter, which is the grating period. Detailed optical micrographs, AFM profiles and SEM images are shown in Supplementary Note 2.
a Schematic of an unpatterned multi-layer WS2 flake on an Au substrate. b Schematic of a multi-layer WS2 grating structure on a gold substrate. Red arrows denote the incident light polarization: TM polarization is defined such that the electric field is perpendicular to the grating, while the TE polarization electric field is parallel to the grating. Grate width (w), period (p), and thickness (t) are defined. c Unpolarized white light micrograph of WS2 flake on Au substrate, showing patterned and unpatterned regions. Different colors appear due to different absorption spectra for different material thicknesses; 15ânm thick WS2 exhibits a purple hue while 25ânm thick WS2 appears pale blue. d Atomic force microscopy (AFM) image of the grating topography corresponding to dashed box of (c). e False color scanning electron microscopy (SEM) image of the boxed region in panel d acquired at 70° tilted angle. The grating structure with dimensions wâ=â300ânm and pâ=â600ânm is clearly seen.
Resonant absorption in nanopatterned WS2 on Au
Reflectance spectra for unpatterned WS2 of varying thickness on an Au back reflector (Fig. 3a, b) exhibit two reflectance dips (absorption resonance peaks). The primary exciton mode of WS2 at 2.0âeV demonstrates high absorption independent of sample thickness and is therefore the focus of this study. While excitons which possess out-of-plane dipoles have also been predicted and observed in monolayer TMDCs19 that have a direct band-gap, the multilayer samples in this study combined with normal incidence illumination prevents any in-coupling of light in reflectance measurements or out-coupling of light in photoluminescence measurements (see Supplementary Note 3). The simulated magnetic field profile indicates strong exciton driven light absorption, suggesting BeerâLambert like absorbance, typical for bulk samples. Additional details on simulation methods are provided in the âMethodsâ section. In contrast, the magnetic field profile representing the FabryâPerot resonance indicates a low-Q cavity mode formed by the combination of metal and high-index dielectric in the heterostructure16,20. The cavity mode is observed for deep-subwavelength thicknesses, owing to the lossy dielectric of WS2 on a reflective substrate facilitating near-unity absorptive resonances. The cavity mode couples to the primary exciton mode at the critical thickness ~15ânm, resulting in overall absorption enhancement. This coupling results in hybridization of the cavity photon and exciton modes into exciton-polariton modes with a characteristic avoided crossing in the reflectance spectra. As the cavity mode is tuned through the exciton energy, the upper exciton-polariton (UEP) branch and lower exciton-polariton branch (LEP) are split21. At the point of avoided crossing of the exciton and the cavity photon, the splitting energy is evaluated to be ħΩâ=â170âmeV by fitting the simulation data to a coupled oscillator model (Supplementary Note 4). Experimental reflectance spectra with varying thicknesses of WS2 show remarkably good qualitative and quantitative agreement with the simulations (Fig. 3b). Both the primary exciton peak and the weak B exciton peak (λâ~â515ânm) locations match well. Note that this coupling between a weak cavity mode and an exciton mode in the deep subwavelength thickness regime can be extended to other TMDC materials and geometries. The coupling also extends to higher order cavity modes for larger thicknesses of WS2 (see Supplementary Note 5). Notably, such behavior cannot be observed in IIIâV semiconductors (e.g., GaP or GaAs) that possess similar band gaps but lack room-temperature excitonic features in their dielectric functions (see Supplementary Note 1). We further note that the optical constants used in our simulations are assumed constant as a function of thickness and are derived from bulk optical constants of WS25. This is because the phenomena reported in this study are focused in the range of 8â25ânm thickness which represents the electronic bulk in WS222 as the effects of electronic quantum confinement become negligible at ~5ânm (6â7 layer) thickness.
a Plot of simulated reflectance spectra of unpatterned WS2 on Au with varying WS2 thickness on the y axis. The primary exciton mode of WS2 can be seen at ~2.0âeV which couples with the cavity mode, splitting the exciton into the upper exciton-polariton (UEP) and the lower exciton-polariton (LEP). A low Q FabryâPerot-like cavity is formed even in deep-subwavelength thickness regime as a result of the thin-film interference. Insets show the HY field profiles corresponding to the exciton mode (left) and the cavity mode (right). Reflectance color scale is the same for (a, c, f). Dashed lines represent a model of two coupled oscillators, while circles represent experimentally measured peak positions shown in (b). b Experimental unpolarized reflectance spectra of unpatterned WS2 on Au with varying thickness. The experimentally measured thicknesses are 5, 10, 14, 16, 21, and 25ânm. The UEP and LEP peaks are emphasized with circles and are superimposed on panel a to show strong matching with simulation. c Simulated reflectance spectra under TM polarization of a 1D WS2 grating on an Au substrate with varying thickness. The structure has a fixed width and period of 300ânm and 500ânm, respectively. The non-dispersive mode induced by dielectric grating emerges as the highest-energy resonance (~500ânm) while plasmon resonances of varying orders emerge at lower energies. For higher thickness regimes, an additional plasmon resonance emerges which red-shifts with increasing thickness. d Experimental TM reflectance spectra of the 1D WS2 grating on the Au substrate with wâ=â300ânm and pâ=â500ânm. The peaks corresponding to various modes are coded with matching symbols on (c). e EZ and HY mode profiles for the various modes matching with the corresponding symbols in (c, d). f Simulated reflectance spectra under TE polarization of a 1D WS2 grating on an Au substrate as thickness is swept.
Periodic patterning provides an additional degree of control for lightâexciton interaction and dispersion engineering (Fig. 3c, d). Under TM polarization, we observe that several additional absorption resonances emerge for a periodically patterned TMDC. Primarily, a resonant mode associated with a dielectric grating emerges at the highest-energies (~500ânm) while several plasmonic modes are excited at lower energies due to the change in sign of the real permittivity across the WS2/Au interface (Fig. 3e)23,24. Characteristic mode profiles are shown in Fig. 3e. Different types of plasmonic mode localization and confinement are observed. For higher thicknesses (more than ~15ânm thick), these plasmonic modes further red-shift with increasing thickness. Detailed discussion on E and H field profiles of the modes is provided below. Our simulation results closely match experimental reflectance spectra for three discrete thicknesses of 10, 14, and 20ânm (Fig. 3d). Under TE polarization (Fig. 3f), the simulated reflectance spectra behave similarly to the unpatterned case, implying that grating and plasmon modes are not excited, since there is no breaking of symmetry along TE polarization upon etching. Experimentally measured reflectance spectra also match the simulations for the TE case, and the spectra, as well as mode profiles, are shown in detail in Supplementary Note 6. All simulations and experimental measurements are for TM polarization for the remainder of this work. We would like to note that the observed modes in this case are non-propagating unlike the prior predicted25 and observed26 guided modes in ultrathin TMDCs. This is also evident from simulations of our field profiles that do not show evidence of any guided-mode resonances.
Plasmonic modes above critical-thickness of WS2
To understand further how the lateral patterning affects mode dispersion, we have investigated the dependence of reflection spectra on grating finger width. At the critical thickness (~15ânm), a plasmon-mode branch emerges from the LEP branch which suggests it has a hybrid character (Fig. 4a). The experimental reflectance spectra for fixed tâ=â15ânm and pâ=â500ânm, and widths of 200ânm and 300ânm match well with the simulation. This along with the non-dispersive grating induced mode and a shallow exciton-polariton mode around 600ânm are all observable in the experimental spectra (Fig. 4b). Beyond the critical thickness regime, more (higher order) plasmon modes emerge from the hybrid LEP mode (Fig. 4c) which is once again verified in the experimental spectra, (Fig. 4d) where a second peak emerges at 695ânm besides the first peak at 740ânm. Grating structures also offer a possibility of tuning frequency of these excited plasmons by changing grating geometrical parameters i.e. width. (Fig. 4e). The grating with its own wave-vector presents momentum matching conditions to excite plasmons in integer multiples of this wavevector. Higher-order surface plasmons (A,B,C in Fig. 4e) emerge in the reflectance spectra as the width is increased, while the plasmon mode indicated by D (Fig. 4d) retains its place as the lowest energy branch. The Ez and Hy field profiles show plasmons of different orders trapped between the WS2/Au interface (Fig. 4f) verifying our assumption. Another interesting observation that we make in our system is the occurrence of plasmonic modes emerging from the interaction of two parent modes27,28. At the crossing point of the plasmon mode labeled D in Fig. 4e with one labeled B, there is a minor dip in intensity followed by re-emergence of the linear dispersion, which is slightly offset from the original plasmon line. We designate these as hybrid mode G retaining character of both the B and D plasmon branches. This can be clearly verified by the mode profile in Fig. 4f where both the Hy and Ez field profiles of mode G are a combination of the field profiles of modes B and D. Extensive width dependence of these modes is shown in detail as a set of figures in Supplementary Note 7 and Supplementary video/animation 1. We observe two different types of plasmon polaritons in our system. Mode D is representative of localized resonance which is invariant of grating width whereas A, B and C have linear dependence on resonator width and hence can be referred to as surface plasmon modes. This suggests that modes A, B and C are excited by the dielectric discontinuity created by the WS2 (εâ>â0) resonator at the Au (εâ<â0)/air surface and are confined at the interface between Au and WS2. Hence their resonances occur at widths equal to odd integers of half of resonance wavelengths as also evident from their field profiles in Fig. 4f. Note that the plasmons observed in our system may be correlated to different types of plasmon modes (e.g., localized and surface plasmons)1. However, here due to the complex nature and multipartite interactions involved, we avoid making any such categorical classifications.
a Width dependence of simulated TM reflectance spectra. Period is fixed at 500ânm, and thickness is fixed at 15ânm (critical thickness). The plasmonic peak emerges from the lower exciton-polariton (LEP) of the unpatterned case indicating a hybrid nature of the mode. Reflectance color scale is the same for (a, c, e). b Experimental reflectance spectra of the WS2 grating structure with tâ=â15ânm and pâ=â500ânm, and widths of 200ânm and 300ânm. Corresponding mode peaks are coded with matching symbols in (a). Inset is the EZ field profile of the first order plasmon mode. The symbols on peaks measured in experimental spectra are marked on the simulated plots in (a). c Width dependence of simulated TM reflectance spectra. Period is fixed at 500ânm, and thickness is fixed at 20ânm (above critical thickness i.e. uncoupled regime). Since the thickness is sufficient to support cavity modes in the unpatterned case, plasmon modes (marked with square symbols) start emerging from the cavity mode. d Experimental reflectance spectra of the WS2 grating structure with thickness and period fixed at 20ânm and 500ânm respectively. Spectra for widths of 200ânm and 300ânm are compared. Inset are the Ez field profiles of the two plasmonic modes. The symbols on peaks measured in experimental spectra are marked on the simulated plots in (c). e Width dependence of simulated TM reflectance spectra. Period is fixed at 500ânm, and thickness is fixed at 30ânm, sufficiently above the critical thickness regime. For a photon wavelength of 750ânm, the resonant widths are equal to 100ânm (A), 300ânm (B), and 470ânm (C), forming 1st, 2nd, and 3rd order plasmon modes. The plasmon mode (D) stands at the lowest-energy resonance. The mode G marked at the crossing point of the 2nd order plasmon mode (B) and (D) indicates a hybrid mode with characters of both. f EZ, HY field profiles of the plasmonic modes labeled with corresponding points in (e).
Strong three-oscillator coupling below critical thickness
The above discussion summarizes the emergence of various plasmonic modes above critical thickness and their hybridization with each other. Thus far there has been no mention of the dielectric grating mode and its interaction with the exciton. This interaction is strong and notably observed for thicknesses below the critical value. The incident wave couples to the guided mode, and the dielectric-grating mode is formed under the momentum matching condition between the incident light-frequency (Ï) and in-plane wave momentum (kx) component generated by the 1D dielectric (WS2) gratings24,26. The dielectric grating induced mode couples the incident light in-plane which allows us to estimate a momentum vs frequency dispersion of the resonances. Below the critical thickness (tâ=â15ânm, or in the strong-coupling regime), the original primary exciton absorption peak is split into LEP (~630ânm) and UEP (~550ânm) modes as a result of strong coupling between the exciton and the cavity mode, as shown in Fig. 5a for the case of 10ânm thick WS2. Both the UEP and the LEP can further couple to the dielectric grating induced mode, forming a three-oscillator coupling29. The three-oscillator plasmon-exciton polaritons coupling has also been reported in several previous works30,31,32. In a structure with a fixed width of 300ânm and thickness of 10ânm WS2, the grating mode can be tuned by varying the period to be in resonance with both UEP and LEP modes. It is worth noting that the UEP and LEP, while coupled as seen in the Fig. 3a, b of the unpatterned case, are period independent unless they interact with the linearly dispersive grating mode. Upon the interaction of grating modes with the two coupled branches of polaritons, strong-coupling occurs among all three modes, leading to redshift of both UEP and LEP modes. The grating mode merges into the non-dispersive upper polariton branch (UPB), while the lower polariton branch (LPB) switches from the non-dispersive LEP mode to plasmonic mode with increasing period. The UEP mode forms the middle polariton branch (MPB), which has a cut-off in the period range of 500ânm to 650ânm. A three-oscillator model is applied to fit the system governed by Eq. (1) below. (Fig. 5c). In this equation, Ed, EUEP and ELEP represent the resonant energies of the dielectric-grating mode, the split upper-exciton-polariton branch and the lower-exciton-polariton branch respectively. g represents the coupling strength between two branches (e.g., gd-UEP represents the coupling strength between the dielectric-grating branch and the UEP branch).
a Period dependence of simulated TM reflectance spectra. Thickness and width are fixed at 10ânm and 300ânm, respectively. Period is swept from the unpatterned case (wâ=âpâ=â300ânm) to a period of 750ânm. The dashed lines represent the calculated absorption spectra based on the three-oscillator model in (c), while the circles represent experimentally measured peak positions from (b). b Experimental reflectance spectra of the WS2 grating structure with a fixed thickness of 10ânm and a width of 300ânm. Spectra for the patterned case, pâ=â500ânm, 600ânm, 650ânm, 700ânm, and 750ânm are plotted and offset for clarity. The marked peaks are also plotted in (a), for comparison with simulation. c The calculated energy dispersion relation based on the three-oscillator model. The three oscillators represent the grating mode, LEP mode, and UEP mode. The bare modes are represented by dashed lines in the figure. By fitting the oscillator model to the simulation and experiment data, we obtain an effective mode separation energy of 410âmeV. d Comparison of the field-line profiles of the UP, MP cut-off, and LP regions marked with triangles, diamonds and crosses, respectively in (a, b). At the point marked with a diamond, almost all the electric field lines are repelled out of the resonator rendering the highly lossy WS2 resonator invisible in the far-field. e Simulated reflectance spectra of grating resonators showing the evolution of the strong-coupling-induced separation as a function of width in the wavelength range of ~550â650ânm. Inset: the EZ field profile of the plasmonic mode for tâ=â10ânm. f Period dependence of reflectance of the WS2 grating structure for tâ=â30ânm and wâ=â300ânm. The three dashed lines representing the bare modes are also shown in the figure, indicating the exciton mode at λâ=â610ânm almost vanishes. The dielectric grating induced mode is unable to couple with the pure exciton without the exciton being hybridized with the cavity photon.
(See Supplementary Note 8 for in-depth discussion of the model and parameters). We would like to note that a similar fitting-based phenomenological description of complex multi-particle interactions has been conducted in many other experimental nanophotonic systems33,34. Because polaritonic coupling dominates the response of this system, we assume gdâLEP to be zero so the effective separation energy can be expressed in terms of the coupling strengths gdâUEP and gLEPâUEP. From the experimental spectrum corresponding to the zero-detuning point (pâ=â600ânm), gdâUEP and gLEPâUEP are taken to be 116âmeV and 89âmeV, respectively (Fig. 5b). When the dielectric grating induced mode and the plasmonic mode are closest in energy, and the MPB is cut-off in between them, the giant effective separation energy which is the maximum experimentally observed separation is estimated to be \(\hbar \Omega = 2\left( {g_{{\mathrm{d}} - {\mathrm{UEP}}} + g_{{\mathrm{LEP}} - {\mathrm{UEP}}}} \right) = 410\,{\rm{meV}}\). It is worth noting that for these nanopatterned grating samples, the linewidth of the dielectric grating mode and the exciton-cavity hybrid mode are γâ=â133âmeV and κâ=â50âmeV, respectively. In comparison, the effective mode separation is 2gâ=â410âmeV which is much larger than (γâ+âκ)/2, further supporting our evidence for strong-coupling35. Interestingly, the MPB cut-off is only replicated in the oscillator model when it is assumed that the MPB oscillator, i.e. the UEP mode, is the only one driven. This suggests the UEP mode is the primary excitation mechanism of this system. This cut-off of the MPB upon strong interaction between these three modes has not been observed in earlier reports of strong lightâmatter coupling of quasiparticles in TMDCs36. From the electromagnetic simulation, we observe that a weak additional plasmon mode begins to split off from the primary plasmon mode at periods near 600ânm (Fig. 5a, left of the black X). Our coupled oscillator model assumes a fixed number of resonators. Hence the emergence of this weak fourth mode is not effectively captured in our model. Based on the electromagnetic simulation and experimental spectra (diamond symbol in Fig. 5a, b), this cut-off of the middle branch, where the absorption reaches almost zero, suggests that the strong mode coupling results in an interference condition that prevents any electromagnetic field interaction with the lossy dielectric part. This is evident from the field line profile (Fig. 5d). As compared with the UPB and LPB of absorption resonances where all the electric-field lines are recirculating or focused in the WS2 resonators (Fig. 5d triangle and cross), in the cut-off middle branch region (absorption dip), most electric-field lines are pushed out of or take the shortest part of interaction through the resonator, resulting in an electromagnetically âinvisibleâ or transparent resonator, qualitatively analogous to atomic systems37,38 (Fig. 5d diamond). Such electromagnetic âinvisibilityâ has not been observed in any lossy photonic system let alone an excitonic, external cavity-less system. Lossy grating structures have been heavily investigated for Si, Ge, and lossy metals based structures39,40,41,42,43 for the purposes of light trapping and structural coloring. In all cases, electromagnetic wave simulations suggest field lines concentrating (focused) inside the lossy dielectric part44 in contrast to our excitonic grating system where the field lines are pushed out of the lossy resonator part (Fig. 5d).
Discussion
The evolution of the strong coupling with thickness (Supplementary Note 9) and resonator width (Fig. 5e) shows that there exists a critical range of resonator widths (50â330ânm) and a critical thickness (15ânm) under which this three-oscillator coupling resulting in cut-off of the MPB is supported. The field profile of the high-Q plasmon mode (Fig. 5e inset) shows significant leakage out of the TMDC, indicating sensitivity to the dielectric environment, with potential applications in chemical and biomolecular sensing45 as well as coupling/tuning quantum emitters in other TMDCs or h-BN46,47,48. Because the grating modeâs field enhancement is localized at the top-corners of the resonator, the grating width is particularly critical for strong-coupling. For pâ=â600ânm and widths larger than 330ânm, since the location of the two modes are too separated, the grating mode lacks the momentum to strongly interact with the exciton, which can explain the fact that no more strong interaction occurs at the point where the grate width is larger than 330ânm. Additionally, the strength of the polariton can be reduced when one dimension of the grating is limited. However, the LEP oscillator strength is much more robust compared to the pure excitons at widths lower than 300ânm (Fig. 5e, f). In comparison to the strongly coupled case (tâ<â15ânm), the pure exciton absorption almost vanishes in the uncoupled regime (tâ>â15ânm) (Fig. 5f). For these thicknesses of WS2, the exciton does not couple to the cavity mode to form hybrid polaritons. The exciton mode at λâ~â600ânm almost vanishes in strength before it can even interact with the grating mode due to its instability in the discontinuous WS2 structure (Fig. 5f). Modes B and Câ correspond to the 2nd and 3rd order plasmon modes while mode D corresponds to the localized plasmonic mode. This indicates that the exciton-polaritons (UEP and LEP modes) or hybrid particles have a much greater robustness and stability to strongly couple to a second photon as opposed to pure excitons. We further stress that in our experiments and observation we observe coupling between purely optical as well as optical and matter resonances. While the latter case is clear and well-studied in both atomic, molecular, and condensed matter systems, the former case has also been a subject of active study49,50. For our unpatterned case, we observe coupling between optical and matter resonances as discussed above. However, as we modify the geometry of the system (patterned case), the geometric dispersion plays an important role causing optical fields to mix and couple further which results in the suppression of the mode vs non-interaction as seen in Fig. 5a vs 5f, respectively.
Our results highlight a previously unexplored aspect of TMDC excitons and optics, and opens a regime in exploring strong lightâmatter interactions where the TMDC serves as a host for both excitonic and cavity modes, without the need for an external cavity. Prior studies have focused on either monolayers or very thick (>100ânm) TMDCs15. Our results show that in this intermediate thickness regime, there is possibility of combining multiple classes of lightâmatter interactions (plasmons, excitons and cavity photons) resulting in strong-coupling phenomena involving two photons or three oscillators resulting in observation of optical invisibility conditions51. Further, our results showing a large effective separation energy (~410âmeV) of the modes at room temperatures can help serve as a platform for exploring a rich diversity of polaritonic phenomena. While multilayer TMDCs are not emissive (see Supplementary Note 4), they can serve as underlying couplers for quantum point defects or other emitting gain media on top, potentially paving a way to exploring interaction between multiple quantum emitters, polariton condensation in molecular or quantum dot assemblies52,53 as well as low-threshold lasers at room temperatures48,52,54,55.
Methods
Sample fabrication
Measurements were performed on a template stripped Au substrate16. To ensure a flat surface of the gold back reflector rather than the rough surface of the evaporated gold, an epoxy-based peeling procedure is applied to the 100ânm-thick Au film evaporated on a clean polished Si wafer (parent wafer) using an e-beam evaporator (Kurt J. Lesker PVD 75). A piece of silicon wafer (transfer wafer) is glued to the noble metal film using a thin layer of thermal epoxy (Epo-Tek 375, Epoxy Technology) and then peeled upwards after the epoxy layer achieves its final hardness, resulting in stripping of the gold film. WS2 was mechanically exfoliated from bulk crystal (HQ-graphene) using Scotch Tape and transferred onto the Au substrate. Flakes with different thicknesses (6â30ânm) were transferred onto the substrate and identified using an optical microscope under reflected white light illumination.
Patterning of the WS2 flakes was achieved by a combination of electron-beam lithography (Elionix ELS-7500EX), exposure of poly methyl metha acrylate (PMMA A4) and dry etching process. Dry-etching was performed using a combination of XeF2 and Ar for 210âs, which is suitable for all the samples with varying thickness from 6ânm to 20ânm to be fully etched. Since XeF2 is highly selective, the etching process does not affect the Au substrate. Finally, the samples were cleaned in acetone to remove the 250ânm-thick PMMA resist. The grating region exhibits a clear change in color compared to the unpatterned region when viewed under an optical microscope in reflection mode under white-light illumination.
Optical characterization
The WS2 reflectance spectra were obtained using a microscope with external white-light illumination (AvaLight-HAL) under normal incidence. The reflected light signals were collected from the microscope objective (Olympus SLMPLN 50X N.A.â=â0.35) and analyzed using a grating spectrometer coupled to a Si focal plane array (FPA). All these instruments are integrated in the LabRAM HR Evolution Confocal Microscope. In order to avoid the signals from the system background and the gold bulk plasmon resonances, the reflectance spectra are normalized by the reflectance of a silver mirror (Thorlabs PF30-03-P01) which can be considered a spectrally uniform high reflector. While the reflectance from the Ag reference mirror is not unity, the minor (few-percent) absorption by the reference mirror will translate to a similarly minor uncertainty in absolute reflectance calibrations. We assume minimum scattering of the incident light due to normal incidence illumination, small sample roughness and very thin features compared to wavelengths. The separated TM and TE polarization signals are obtained using a linear-polarizer (Edmund Linear Glass Polarizing Filter #43-783) for the incident white light. The incident white-light spot size is ~4âμm in diameter (shown in Supplementary Fig. 2a).
Simulations
Numerical simulations of the electromagnetic response of the WS2 gratings were achieved with the Lumerical FDTD and COMSOL Multiphysics solvers. 1-dimensional heterostructures were solved using the transfer matrix method. The reflectance spectra and field distributions were calculated by simulating incident plane waves with TM or TE polarizations. The anisotropic permittivity of bulk WS2 reported in ref. 5,15 was adapted for our simulations.
Data availability
The data that support the conclusions of this study are available from the corresponding author upon request.
Code availability
The Python scripts used in this study for plotting and modeling are available from the corresponding author upon request.
References
Schuller, J. A. et al. Plasmonics for extreme light concentration and manipulation. Nat. Mater. 9, 193â204 (2010).
Lukâyanchuk, B. et al. The Fano resonance in plasmonic nanostructures and metamaterials. Nat. Mater. 9, 707 (2010).
Kuznetsov, A. I., Miroshnichenko, A. E., Brongersma, M. L., Kivshar, Y. S. & Lukâyanchuk, B. Optically resonant dielectric nanostructures. Science 354, aag2472 (2016).
Zhang, L., Mei, S., Huang, K. & Qiu, C. W. Advances in full control of electromagnetic waves with metasurfaces. Adv. Optical Mater. 4, 818â833 (2016).
Li, Y. et al. Measurement of the optical dielectric function of monolayer transition-metal dichalcogenides: MoS2, MoSe2, WS2, and WSe2. Phys. Rev. B 90, 205422 (2014).
Krasnok, A., Lepeshov, S. & Alú, A. Nanophotonics with 2D transition metal dichalcogenides. Opt. Express 26, 15972â15994 (2018).
Zheng, D. et al. Manipulating coherent plasmonâexciton interaction in a single silver nanorod on monolayer WSe2. Nano Lett. 17, 3809â3814 (2017).
Hu, T. et al. Strong coupling between Tamm plasmon polariton and two dimensional semiconductor excitons. Appl. Phys. Lett. 110, 051101 (2017).
Wen, J. et al. Room-temperature strong lightâmatter interaction with active control in single plasmonic nanorod coupled with two-dimensional atomic crystals. Nano Lett. 17, 4689â4697 (2017).
Wang, M. et al. Tunable fano resonance and plasmonâexciton coupling in single au nanotriangles on monolayer WS2 at room temperature. Adv. Mater. 30, 1705779 (2018).
Schneider, C., Glazov, M. M., Korn, T., Höfling, S. & Urbaszek, B. Two-dimensional semiconductors in the regime of strong light-matter coupling. Nat. Commun. 9, 1â9 (2018).
Jahani, S. & Jacob, Z. All-dielectric metamaterials. Nat. Nanotechnol. 11, 23 (2016).
Lundt, N. et al. Room-temperature Tamm-plasmon exciton-polaritons with a WSe2 monolayer. Nat. Commun. 7, 1â6 (2016).
Kleemann, M.-E. et al. Strong-coupling of WSe2 in ultra-compact plasmonic nanocavities at room temperature. Nat. Commun. 8, 1â7 (2017).
Verre, R. et al. Transition metal dichalcogenide nanodisks as high-index dielectric Mie nanoresonators. Nature Nanotechnol. 14, 679â683 (2019).
Jariwala, D. et al. Near-unity absorption in van der Waals semiconductors for ultrathin optoelectronics. Nano Lett. 16, 5482â5487 (2016).
Zhang, L., Gogna, R., Burg, W., Tutuc, E. & Deng, H. Photonic-crystal exciton-polaritons in monolayer semiconductors. Nat. Commun. 9, 713 (2018).
Novoselov, K. S. et al. Two-dimensional atomic crystals. Proc. Natl Acad. Sci. USA 102, 10451â10453 (2005).
Loh, K. P. Brightening the dark excitons. Nat. Nanotechnol. 12, 837â838 (2017).
Kats, M. A., Blanchard, R., Genevet, P. & Capasso, F. Nanometre optical coatings based on strong interference effects in highly absorbing media. Nat. Mater. 12, 20 (2013).
Flatten, L. C. et al. Room-temperature exciton-polaritons with two-dimensional WS2. Sci. Rep. 6, 33134 (2016).
Zhao, W. et al. Evolution of electronic structure in atomically thin sheets of WS2 and WSe2. ACS Nano 7, 791â797 (2013).
Pi, S. et al. Dielectric-grating-coupled surface plasmon resonance from the back side of the metal film for ultrasensitive sensing. IEEE Photonics J. 8, 1â7 (2015).
Schröter, U. & Heitmann, D. Grating couplers for surface plasmons excited on thin metal films in the Kretschmann-Raether configuration. Phys. Rev. B 60, 4992 (1999).
Khurgin, J. B. Two-dimensional excitonâpolaritonâlight guiding by transition metal dichalcogenide monolayers. Optica 2, 740â742 (2015).
Zhang, X. et al. Guiding of visible photons at the Ã¥ngström thickness limit. Nature Nanotechnol. 14, 844â850 (2019).
Prodan, E., Radloff, C., Halas, N. J. & Nordlander, P. A hybridization model for the plasmon response of complex nanostructures. science 302, 419â422 (2003).
Fofang, N. T. et al. Plexcitonic nanoparticles: plasmonâexciton coupling in nanoshell-J-aggregate complexes. Nano Lett. 8, 3481â3487 (2008).
Mukherjee, S. et al. Fanoshells: nanoparticles with built-in Fano resonances. Nano Lett. 10, 2694â2701 (2010).
Rodriguez, S. & Rivas, J. G. Surface lattice resonances strongly coupled to Rhodamine 6G excitons: tuning the plasmon-exciton-polariton mass and composition. Opt. Express 21, 27411â27421 (2013).
Väkeväinen, A. et al. Plasmonic surface lattice resonances at the strong coupling regime. Nano Lett. 14, 1721â1727 (2014).
Todisco, F. et al. Toward cavity quantum electrodynamics with hybrid photon gap-plasmon states. ACS Nano 10, 11360â11368 (2016).
Artar, A., Yanik, A. A. & Altug, H. Multispectral plasmon induced transparency in coupled meta-atoms. Nano Lett. 11, 1685â1689 (2011).
Artar, A., Yanik, A. A. & Altug, H. Directional double Fano resonances in plasmonic hetero-oligomers. Nano Lett. 11, 3694â3700 (2011).
Khitrova, G., Gibbs, H., Kira, M., Koch, S. W. & Scherer, A. Vacuum Rabi splitting in semiconductors. Nat. Phys. 2, 81â90 (2006).
Li, B. et al. Large Rabi splitting obtained in Ag-WS2 strong-coupling heterostructure with optical microcavity at room temperature. Opto-Electron. Adv. 2, 190008 (2019).
Boller, K.-J., ImamoÄlu, A. & Harris, S. E. Observation of electromagnetically induced transparency. Phys. Rev. Lett. 66, 2593 (1991).
Fleischhauer, M. & Lukin, M. D. Dark-state polaritons in electromagnetically induced transparency. Phys. Rev. Lett. 84, 5094 (2000).
Pala, R. A. et al. Optimization of non-periodic plasmonic light-trapping layers for thin-film solar cells. Nat. Commun. 4, 2095 (2013).
Li, Z., Palacios, E., Butun, S. & Aydin, K. Visible-frequency metasurfaces for broadband anomalous reflection and high-efficiency spectrum splitting. Nano Lett. 15, 1615â1621 (2015).
Aydin, K., Ferry, V. E., Briggs, R. M. & Atwater, H. A. Broadband polarization-independent resonant light absorption using ultrathin plasmonic super absorbers. Nat. Commun. 2, 517 (2011).
Cao, L., Fan, P., Barnard, E. S., Brown, A. M. & Brongersma, M. L. Tuning the color of silicon nanostructures. Nano Lett. 10, 2649â2654 (2010).
Gu, Y., Zhang, L., Yang, J. K., Yeo, S. P. & Qiu, C.-W. Color generation via subwavelength plasmonic nanostructures. Nanoscale 7, 6409â6419 (2015).
Kim, S. J., Fan, P., Kang, J.-H. & Brongersma, M. L. Creating semiconductor metafilms with designer absorption spectra. Nat. Commun. 6, 7591 (2015).
Rodrigo, D. et al. Mid-infrared plasmonic biosensing with graphene. Science 349, 165â168 (2015).
Ren, S., Tan, Q. & Zhang, J. Review on the quantum emitters in two-dimensional materials. J. Semiconductors 40, 071903 (2019).
He, Y.-M. et al. Single quantum emitters in monolayer semiconductors. Nat. Nanotechnol. 10, 497 (2015).
Wu, S. et al. Monolayer semiconductor nanocavity lasers with ultralow thresholds. Nature 520, 69 (2015).
Yang, Y., Kravchenko, I. I., Briggs, D. P. & Valentine, J. All-dielectric metasurface analogue of electromagnetically induced transparency. Nat. Commun. 5, 1â7 (2014).
Staude, I. et al. Tailoring directional scattering through magnetic and electric resonances in subwavelength silicon nanodisks. ACS Nano 7, 7824â7832 (2013).
Papasimakis, N., Fedotov, V. A., Zheludev, N. & Prosvirnin, S. Metamaterial analog of electromagnetically induced transparency. Phys. Rev. Lett. 101, 253903 (2008).
Deng, H., Haug, H. & Yamamoto, Y. Exciton-polariton bose-einstein condensation. Rev. Mod. Phys. 82, 1489 (2010).
Balili, R., Hartwell, V., Snoke, D., Pfeiffer, L. & West, K. Bose-Einstein condensation of microcavity polaritons in a trap. Science 316, 1007â1010 (2007).
Deng, H., Weihs, G., Snoke, D., Bloch, J. & Yamamoto, Y. Polariton lasing vs. photon lasing in a semiconductor microcavity. Proc. Natl Acad. Sci. USA 100, 15318â15323 (2003).
Zhou, W. et al. Lasing action in strongly coupled plasmonic nanocavity arrays. Nat. Nano 8, 506â511 (2013).
Acknowledgements
D.J. acknowledges primary support for this work by the U.S. Army Research Office under contract number W911NF-19-1-0109. D.J. also acknowledges support from Penn Engineering Start-up funds, National Science Foundation (DMR-1905853) and University of Pennsylvania Materials Research Science and Engineering Center (MRSEC) (DMR-1720530) and Northrop Grumman Corporation. B.A. acknowledges support from Vagelos Integrated Program in Energy Research at Penn. D.J. and S.R. acknowledge support from Center for Undergraduate Research and Fellowships at Penn. A.R.D. acknowledges start-up funds from UCLA. This work was carried out in its entirety at the Singh Center for Nanotechnology at the University of Pennsylvania which is supported by the National Science Foundation (NSF) National Nanotechnology Coordinated Infrastructure Program grant NNCI-1542153. The authors acknowledge assistance from Xingdu Qiao for SEM imaging.
Author information
Authors and Affiliations
Contributions
D.J. and H.Z. conceived the project. H.Z. made the samples, performed microfabrication, reflectance measurements, and electron microscopy. B.A. performed all simulations and modeling with help from A.R.D., M.W.K., Q.Z., H.Z., and S.R. J.M. assisted with sample preparation and microfabrication. K.J. performed atomic force microscopy characterization. D.J., H.Z., B.A., and A.R.D. analyzed and interpreted the optical spectroscopy and simulation data. D.J. supervised the study. All others contributed to writing of the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Seyedeh Mehri Hamidi, Jun Song, and the other, anonymous reviewer(s) for their contribution to the peer review of this work.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhang, H., Abhiraman, B., Zhang, Q. et al. Hybrid exciton-plasmon-polaritons in van der Waals semiconductor gratings. Nat Commun 11, 3552 (2020). https://doi.org/10.1038/s41467-020-17313-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-17313-2
This article is cited by
-
Extinction Behavior and Near-Field Enhancement of 3D Hybrid Nanostructures Consisting of Ag-Coated SiO2 Nanorods
Plasmonics (2024)
-
Recent progress of exciton transport in two-dimensional semiconductors
Nano Convergence (2023)
-
Efficient Manipulation of Strong Coupling Towards the Deep Ultraviolet Region
Plasmonics (2023)
-
Ultra-confined Propagating ExcitonâPlasmon Polaritons Enabled by Cavity-Free Strong Coupling: Beating Plasmonic Trade-Offs
Nanoscale Research Letters (2022)
-
Cavity-enhanced linear dichroism in a van der Waals antiferromagnet
Nature Photonics (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.