
==== Front
Nanophotonics
Nanophotonics
nanoph
nanoph
Nanophotonics
2192-8606
2192-8614
De Gruyter

38681678
nanoph-2023-0604
10.1515/nanoph-2023-0604
Research Article
Sculpting ultrastrong light–matter coupling through spatial matter structuring
https://orcid.org/0000-0001-5984-3190
Mornhinweg Joshua joshua.mornhinweg@ur.de

Diebel Laura Laura-Katharina.Diebel@physik.uni-regensburg.de

Halbhuber Maike maike.halbhuber@ur.de

Riepl Josef josef1.riepl@ur.de

Cortese Erika E.Cortese@soton.ac.uk

https://orcid.org/0000-0002-4851-2633
De Liberato Simone S.De-Liberato@soton.ac.uk

Bougeard Dominique Dominique.Bougeard@physik.uni-regensburg.de

https://orcid.org/0000-0001-6617-9283
Huber Rupert rupert.huber@ur.de

https://orcid.org/0000-0002-2134-6612
Lange Christoph christoph.lange@tu-dortmund.de

Department of Physics, University of Regensburg, 93040 Regensburg, Germany
School of Physics and Astronomy, University of Southampton, Southampton, SO17 1BJ, UK
IFN – Istituto di Fotonica e Nanotecnologie, CNR, I-20133 Milan, Italy
Department of Physics, TU Dortmund University, 44227 Dortmund, Germany
Corresponding authors: Joshua Mornhinweg, Department of Physics, University of Regensburg, 93040 Regensburg, Germany; and Department of Physics, TU Dortmund University, 44227 Dortmund, Germany, E-mail: joshua.mornhinweg@ur.de; and Christoph Lange, Department of Physics, TU Dortmund University, 44227 Dortmund, Germany, E-mail: christoph.lange@tu-dortmund.de
11 1 2024
4 2024
13 10 19091915
18 9 2023
10 12 2023
© 2024 the author(s), published by De Gruyter, Berlin/Boston
2024
the author(s), published by De Gruyter, Berlin/Boston GmbH, Berlin/Boston
https://creativecommons.org/licenses/by/4.0/ This work is licensed under the Creative Commons Attribution 4.0 International License.
Abstract

The central theme of cavity quantum electrodynamics is the coupling of a single optical mode with a single matter excitation, leading to a doublet of cavity polaritons which govern the optical properties of the coupled structure. Especially in the ultrastrong coupling regime, where the ratio of the vacuum Rabi frequency and the quasi-resonant carrier frequency of light, ΩR/ω c, approaches unity, the polariton doublet bridges a large spectral bandwidth 2ΩR, and further interactions with off-resonant light and matter modes may occur. The resulting multi-mode coupling has recently attracted attention owing to the additional degrees of freedom for designing light–matter coupled resonances, despite added complexity. Here, we experimentally implement a novel strategy to sculpt ultrastrong multi-mode coupling by tailoring the spatial overlap of multiple modes of planar metallic THz resonators and the cyclotron resonances of Landau-quantized two-dimensional electrons, on subwavelength scales. We show that similarly to the selection rules of classical optics, this allows us to suppress or enhance certain coupling pathways and to control the number of light–matter coupled modes, their octave-spanning frequency spectra, and their response to magnetic tuning. This offers novel pathways for controlling dissipation, tailoring quantum light sources, nonlinearities, correlations as well as entanglement in quantum information processing.

Keywords

cavity-QED
THz metasurfaces
multi-mode coupling
polaritonic control
Leverhulme TrustRPG-2022-037 Future and Emerging Technologies (FET)FET 737017 (MIR-BOSE) Deutsche Forschungsgemeinschaft422 31469 5032-SFB1277 (Subproject A01) BO 3140/3-2 HU 1598/8 LA 3307/1-2
==== Body
pmc1 Introduction

Large light–matter coupling strengths have opened up entirely new perspectives for cavity quantum electrodynamics (c-QED): In ultrastrongly coupled structures, exciting quantum phenomena such as the vacuum Bloch–Siegert shift [1], cavity quantum chemistry [2], [3], cavity-controlled electronic transport [4], [5] or the creation of photon-bound excitons [6] have been observed and applications in nanophotonics including novel light sources [7], [8], [9], single-photon manipulation [10] and squeezed quantum states of light [11] are being explored. Harnessing such effects requires control of the polariton frequencies, line shapes and oscillator strengths. These properties ultimately result from the underlying exotic ground state of ultrastrong coupling that is characterized by a non-vanishing expectation value of vacuum excitations as well as correlations and entanglement of optical and electronic modes [12], [13]. Landau cavity polaritons combining planar metal cavities with strongly subwavelength field confinement and the large dipole moments of cyclotron resonances (CR) of two-dimensional electron gases in semiconductor quantum wells (QWs) [1], [4], [14]–[19] have been a highly successful route. The intrinsically large dipole moments of CRs have enabled record coupling strengths of up to ΩR/ω c = 3.19 [19] as well as ultrastrong coupling in few-electron systems using nanogap cavities [18], [20] or ultrastrongly coupled single meta-atoms [21]. Generally, the optical and electronic properties of a light–matter coupled system are defined by the overlap integral of the electric field of the resonator mode and the polarization field of the electronic resonance, each with their specific subwavelength structure [22].

Whereas most investigations have thus far focused on maximising the light–matter coupling strength of a single pair of light–matter coupled modes, or controlling their far-field properties such as polariton frequencies, the large design space accessible by subwavelength control of the shape and overlap of multiple modes has been largely neglected. First explorations include the very-strong coupling regime [6], [23]–[30] or nanophotonic polariton control in equilibrium [31]. The large frequency range of 2ΩR ∼ ω c covered by polaritons of ultrastrongly coupled structures moreover necessitates control of the coupling mechanism over a correspondingly broad spectral range, which can include resonances more than one optical octave above the fundamental cavity mode [19]. Strong nonlinearities and nonlinear polariton mixing were observed in structures featuring multiple polariton modes and corresponding generalized vacuum Rabi frequencies [32]. Moreover, subcycle control of ΩR was demonstrated by femtosecond reshaping and switch-off of cavity modes [33], predicted to release the virtual photon population of the exotic vacuum ground state of ultrastrong coupling [33]–[36]. These examples underscore the importance of the spatial structure of the coupled fields for controlling ultrastrong coupling and its vacuum ground state.

Here, we experimentally explore this design space by tailoring the spatial mode overlap of multiple light and matter modes by structuring the matter component. In similarity to selection rules of classical optics, this new degree of freedom enables us to select particular coupling pathways: nanostructuring controls the extent to which each mode contributes to light–matter hybridization, boosts or suppresses coupling of certain modes, and designs the spectral shape of multiple polariton resonances within a large frequency range. Such an ability is particularly important for large coupling strengths ΩR/ω c, where coupling to multiple modes over an octave-spanning spectral range becomes inevitable and leads to significant resonant and non-resonant interactions [19], [24], [25], [26], [28], [29], [30], [37]. As a result, multiple coupling pathways influence the non-trivial structure of the vacuum state of ultrastrong coupling at the same time. We demonstrate this concept for two structures which ultrastrongly couple the light field of planar metallic THz resonators to the CR of Landau-quantized two-dimensional electron gases hosted in semiconductor QWs. The reference sample implements a planar QW and is referred to as unstructured, while the QW of the second sample, referred to as structured, is laterally etched to control the overlap with the resonator modes. Our analysis of experimental data and theory calculations reveals that even without modifying the cavity, this sub-wavelength structuring of the electronic resonance allows for controlling the frequencies, magnetic field dependencies, light–matter coupling strengths, and the number of light–matter coupled modes, and correspondingly, the exotic vacuum ground state. In the future, this idea may increase the level of flexibility for designing multi-mode polaritonic systems and their light–matter couplings and dispersion over multiple octaves in the deep-strong [17], [19], [33] or very-strong coupling regime [6], [23]–[30].

2 Results and discussion

Our resonators (Figure 1a) feature rectangular double gaps in the center that provide a capacitive part and are connected to inductive loops defining the outer rectangular perimeter [17]. Within the spectral range of interest here, this design supports a fundamental mode which is characterized by periodic energy exchange between the capacitive and the inductive elements at a frequency of ν 1 = 0.8 THz (LC mode, Figure 1b), and correspondingly, a strong near-field enhancement in the center region, which exhibits a phase difference of ∼π relative to the outer parts of the resonator (Figure 2b). A second mode (dipolar mode) with a resonance frequency of ν 2 = 1.6 THz (Figure 1b) partially bypasses the current path through the capacitance and features a more delocalized near-field enhancement which extends to the outer corners of the structure and lacks the phase difference found for the LC mode (Figure 2c). These characteristic spatial field distributions are the key to designing the overlap of resonator modes and electronic polarization by structuring the matter system. In particular, they enable control of coupling pathways by continuously tuning the overlap of the cavity modes from zero to unity (Figure 2d–f). Our semiconductor heterostructures consist of 3 QWs of a thickness of 20 nm, each n-doped with an electron density of ρ QW = 1.25 × 1012 cm−2, and separated by AlGaAs barriers of a thickness of 25 nm. Structuring of the QWs is performed by first covering the heterostructure with resist which is subsequently patterned by electron-beam lithography to form a temporary mask for lateral wet-etching, resulting in a periodic array of rectangular QW patches each measuring 15 µm by 15 µm (Figure 2a, red region). Subsequently, the resonators are fabricated by electron-beam lithography, metal evaporation and a lift-off process, whereby the field-concentrating central gap region is aligned with the etched patches. For the reference sample, the resonators are fabricated on an unstructured QW heterostructure.

Figure 1: Ultrastrongly coupled sample structure. (a) Schematic of the structure consisting of an array of THz resonators (gold shapes), fabricated on top of the GaAs substrate (dark brown). The transmission is probed by broadband, few-cycle THz transients (blue waveform). (b) Simulated transmission spectrum of the bare resonator structure. LC: fundamental mode with a frequency of ν 1 = 0.8 THz. DP: higher-order mode with a frequency of ν 2 = 1.6 THz.

Figure 2: Mode-selective ultrastrong coupling by nanostructured matter resonances. (a) Schematic of a single THz resonator (gold shape) fabricated on top of the GaAs substrate (dark brown), and the cyclotron resonances hosted in spatially confined, high-quality GaAs quantum well structures (bright red patch). (b) Simulated x-polarized near-field distribution of the LC mode of the bare resonator, and (c) the dipolar mode. (d) Schematic of the coupling with zero mode overlap, (e) partial overlap and (f), full overlap of the 1st (LC, ν = 1) and 2nd (DP, ν = 2) cavity modes.

Transmission measurements are performed with linearly polarized, single-cycle THz waveforms generated and detected in ZnTe crystals of a thickness of 0.5 mm and 1 mm, respectively. The corresponding amplitude spectra are obtained by Fourier transform and feature a spectral bandwidth ranging from ∼200 GHz up to >2 THz. The samples are kept at cryogenic temperatures in a magnet cryostat which provides a static magnetic bias field B of up to 5.5 T applied perpendicularly to the QW plane, introducing Landau quantization of the electrons. The resulting cyclotron resonances are tuneable in frequency, ν c = eB/2πm *, where e is the elementary charge and m * = 0.07 m e and m e are the electron effective and free masses, respectively.

We first investigate the unstructured sample implementing the planar QWs, in which the CR couples to both resonator modes without spatial selectivity. In Figure 3a, the spectra are plotted as a function of ν c as a color plot and reveal five distinct resonances. The diagonal feature at ν = ν c (dashed red line) originates from the CR in uncoupled areas of the structure between resonators. The remaining four dominant modes originate from the light–matter coupling. The lower (LP1) and upper polariton (UP1) modes associated with the first resonator mode emerge at frequencies of 0 THz and 0.9 THz, respectively, for ν c = 0 THz, and increase in frequency with opposite curvatures as ν c is increased, forming the typical anti-crossing shape. Similarly, the corresponding LP2 mode of the second resonator mode branches off the CR near ν = ν c ≈ 1.2 THz, further increasing in frequency with a declining slope as ν c increases. The associated UP2 starts at ν = 1.75 THz, for ν c = 0 THz, whereby its initially vanishing slope increases with increasing ν c.

Figure 3: Multi-mode coupling with an unstructured QW film. (a) Experimental transmission spectra of the unstructured sample as a function of the cyclotron frequency ν c. The semi-transparent curves represent the eigenmodes calculated by theoretical fit. (b) FEFD simulation of the transmission spectra of the structure as a function of the cyclotron frequency ν c, including the eigenmodes of the theoretical fit in panel (a). The two arrows indicate the frequencies of the uncoupled LC and DP modes obtained by our fit. Outset: Schematic of the resonator on top of an unstructured QW (bright red area).

While coupling to the CR requires a vanishing in-plane momentum of the light field, the cavity furthermore generates non-zero in-plane wave vectors which are discretized in energy, resulting from the confinement of the near-field to the central gap region and its periodic structure given by the resonator array, respectively [19]. These components couple to magnetoplasmon modes, resulting in additional light–matter coupled magnetoplasmon-polaritons, in general. Here, we observe only one additional resonance of sizeable oscillator strength, at approximately 0.3 THz above the UP1 mode. Since this feature does not significantly influence the dependence of polariton formation on the spatial structure of the matter component, we restrict our analysis to the LP and UP modes discussed above.

Next, we measure the transmission of the structured sample (Figure 4a) where the QW polarization can couple to the near-field only in the central gap area (cf. Figure 2a) and hence, the CR signal observed in Figure 3 is absent. The transmission spectra reveal a fundamentally different polariton spectrum characterized by three rather than four coupled modes. The resonances associated with UP1 and LP2 for the first structure here merge into a single coupled resonance. At ν c = 0, this mode has a frequency of 0.92 THz which remains approximately constant up to ν c = 0.7 THz, from where the frequency increases to the inflection point at ν c = 1.25 THz. As ν c is further increased, the frequency converges towards ν = 1.6 THz. The magnetic tuning function ννc of this novel feature resembles an S-shape without anti-crossing as theoretically predicted in ref. [31]. In addition, also the frequency response of the other coupled resonances differs from the situation in the unstructured sample, while the absorption strength of the polariton resonances remains comparable (see Supplementary Material). In more universal terms, the very large coupling strengths of, e.g., the deep-strong coupling regime, inevitably cause all fundamental resonances to influence each other even across an ultrabroadband spectral range. Tailoring the spatial overlap of these modes by near-field structuring allows us to take control of the of individual coupling constants that link pairs of optical and electronic modes of a light–matter coupled structure, in analogy to selection rules of classical optics (Figure 2d–f).

Figure 4: Multi-mode coupling with a structured QW film. (a) Experimental transmission spectra of the unstructured sample, each normalized to the average transmission Tavg for each frequency ν, as a function of the cyclotron frequency ν c. The semi-transparent curves represent the eigenmodes calculated by theoretical fit. (b) FEFD simulation of the transmission spectra of the structure, each normalized to the average transmission Tavg for each frequency ν as a function of the cyclotron frequency ν c, including the eigenmodes of the theoretical fit in panel (a). The two arrows indicate the frequencies of the uncoupled LC and DP modes obtained by our fit. Outset: Schematic of the resonator on top of a structured QW patch (bright red patch).

In a first approach, we describe our data by a classical electrodynamical formalism employing finite-element frequency-domain (FEFD) simulations [17], which reproduce the transmission spectra and the frequency response of all modes with high accuracy, without free fit parameters (Figures 3b and 4b). These calculations furthermore deliver the spatial profiles of uncoupled (see Figure 2) and coupled modes, enabling the identification of the polariton resonances. To explain the observed differences between both structures and extract the coupling and overlap parameters, we developed a novel quantum mechanical formalism reported in detail elsewhere [31], which extends the established Hopfield model to non-orthogonal modes in order to take the fractional overlap of light and matter modes of our structured sample into account. The bosonic Hamiltonian reads: (1) H^=∑νℏωνa^ν†a^ν+∑νℏωcb^ν†b^ν+∑ν∑μ≤νℏΩR,ν,μb^μ+ΩR,ν,μ*b^μ†a^ν†+a^ν+∑ν,μhν,μa^ν†+a^νa^μ†+a^μ,

where the first two sums describe the bare cavity and matter systems. Here, the bosonic annihilation operator a^ν is associated with the ν-th cavity mode with a frequency ω ν = 2πν ν, and b^μ is the collective bosonic matter operator which describes the μ-th degenerate Landau excitation with frequency ω c = 2πν c. The third double sum describes coupling pathways between all light and matter modes. Here, (2) hν,μ=∑γ≤ν,μΩR,ν,γΩR,μ,γωc

contains the coupling parameters ΩR,ν,γ , i.e., the vacuum Rabi frequencies quantifying the spatial overlap and thus the mutual interaction between modes ν, γ, whereby ΩR,ν,γωc quantifies the associated normalized coupling strength. Moreover, this term accounts for interactions caused by the non-orthogonality of modes in the limited subspace of interaction, i.e., the QW plane [31], as compared to the full three-dimensional space, where orthogonality holds. Note that within a local response approximation, we consider the matter modes as degenerate since they exhibit different real-space wavefunctions, whereas in systems with ultra-narrow features one would be obliged to take the dispersion of the magneto-polaritons into account [30]. The last sum of Eq. (1) is the diamagnetic term that describes the cavity blue shifts for each resonator mode arising from self-interaction. Owing to the bosonic nature of the Hamiltonian, the standard Hopfield diagonalization can be performed.

As shown in greater detail in [31], for two photonic modes we can express the coupling strengths in terms of a single overlap parameter η 2,1 which we calculate based on the electric near-field distribution of the resonator modes (see Supplementary Material) and which assumes values between 0 and 1: (3) ΩR,1,1,ΩR,2,1=Ω~R,2η2,1,ΩR,2,2=Ω~R,2⋅1−η2,12.

For the matter excitations, we chose a basis for which the photonic mode ν = 1 is always exclusively coupled to the matter mode μ = 1, implying that ΩR,1,2 = 0, while the coupling of the second mode, Ω~R,2 , is redistributed between both matter modes for a non-zero overlap parameter η 2,1. While higher-order photonic modes are present in our structures, their coupling strength and influence is negligible (see Supplementary Material).

We employ our multi-mode formalism to fit the resonances of the spectra calculated as a function of ν c (Figures 3b, and 4b, black curves). The coupling parameters ΩR,1,1 and Ω~R,2 serve as fitting parameters. Comparing these spectra for the two configurations enables us to determine the key differences. First, only the fit in Figure 3b exhibits the S-shaped magnetic tuning curve. This distinct feature can be understood by comparing the overlap parameter which assumes η 2,1 = 0.15 for the unstructured system and η 2,1 ≅ 1, for the structured sample. Considering that in principle, 0 < η 2,1 < 1, these non-zero overlap factors show that interactions of the photonic modes are present also in the case of the unstructured sample. Whereas the global orthogonality of the resonator modes rules out interactions, the restriction of the light–matter interaction volume to a finite sub-domain represented either by the full QW plane or by the etched QW patch lifts this restriction, allowing us to tailor the interaction pathways to a large degree. Correspondingly, the two values of η 2,1 lead to strongly different polariton spectra, where the S-shaped spectral signature is a hallmark of nearly full mode overlap, η 2,1 ≈ 1. Our model shows that this spectral shape indicates a change of character of the light–matter coupling mechanism which manifests in a reduction of the number of participating modes from 4 to 3 as η 2,1 approaches unity. This can be formally seen from Eq. (3), where ΩR,2,2 = 0 for η 2,1 = 1, i.e., the interaction between the second photon mode and the second electronic mode (orthogonal to and degenerate with the first one), vanishes, and the cross-interaction term proportional to ΩR,2,1 becomes dominant. Therefore, the same electronic mode interacts with both photon modes ν = 1 and ν = 2, which confine the resulting coupled mode in the spectral range enclosed by them, force the formation of an inflection point, and mark the asymptotic frequencies attained for ν c → 0 and ν c → 2 THz, respectively. The uncoupled electronic mode follows ν = ν c and thus produces a diagonal spectral signature identical to that of the CR in areas between resonator structures. As an additional result, we observe a general reduction of both the polariton splitting as well as the two modes’ diamagnetic blue shifts, resulting in reduced normalized coupling strengths of ΩR,1,1ω1=0.28 , ΩR.2,2ω2=6×10−3 , and ΩR,2,1ω1=0.27 , as compared to the unstructured case where ΩR,1,1ω1=0.37 , ΩR.2,2ω2=0.21 , and ΩR,2,1ω1=0.07 . While overall, the coupling strengths decrease when structuring the QWs, the relative weight of the DP mode, indicated by the ratio ΩR,2,1/ΩR,1,1, increases from 0.19 to 0.96 upon introducing the QW patch. A simulation systematically sweeping the patch size shows how the parameters change from minimum to maximum mode overlap (see Supplementary Material).

3 Conclusions

In conclusion, the concept of lateral confinement of the quantum wells adds a previously unexplored parameter space for tailoring ultrastrong light–matter coupling by controlling the spatial mode overlap in a situation of multiple interacting light and matter modes, across several optical octaves. As we have experimentally demonstrated, the character of the resulting polariton modes is strongly controlled by the overlap, including polariton frequencies, magnetic field dependencies, and spatial field distributions which are relevant for nanophotonic applications. Since for deep-strong coupling conditions, coupling of multiple modes over an ultrabroadband range is almost inevitable, future implementations following our approach may benefit from the ability to selectively boost or suppress certain coupling pathways, in similarity to tailoring of selection rules in classical optics. This novel design approach increases the flexibility for sculpting polaritonic structures by exploiting mutual couplings, controlling modal dispersion, or tailoring resonator structures with strongly non-overlapping optical modes.

Supplementary Material

Supplementary Material Details

Acknowledgments

The authors thank Dieter Schuh and Imke Gronwald for valuable discussions and technical support.

Supplementary Material

This article contains supplementary material (https://doi.org/10.1515/nanoph-2023-0604).

Research funding: We gratefully acknowledge support by the Deutsche Forschungsgemeinschaft (dx.doi.org/10.13039/501100001659) through Project IDs 422 31469 5032-SFB1277 (Subproject A01), grants no. BO 3140/3-2, LA 3307/1-2, and HU 1598/8, as well as by the European Research Council (ERC) (dx.doi.org/10.13039/100010663) through Future and Emerging Technologies (FET) grant no. 737017 (MIR-BOSE) and by the Leverhulme Trust thought the grant RPG-2022-037.

Author contributions: J.M., D.B., R.H. and C.L. conceived the study. L.D., M.H. and D.B. designed, realized, and characterized the semiconductor heterostructures. J.M. and L.D. modelled the metasurfaces and fabricated the samples with support from M.H., J.R., D.B. and C.L. J.M., M.H., L.D. and J.R. carried out the experiments with support from R.H. and C.L. The theoretical modelling was carried out by E.C. and S.D.L., S.D.L., D.B., R.H. and C.L. supervised the study. All authors analyzed the data and discussed the results. J.M., E.C., S.D.L., R.H. and C.L. wrote the manuscript with contributions from all authors. All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

Conflict of interest: Authors state no conflicts of interest.

Data availability: The datasets generated and/or analyzed during the current study are available from the corresponding author upon reasonable request.
==== Refs
References

[1] Li X. Vacuum Bloch–Siegert shift in Landau polaritons with ultra-high cooperativity Nat. Photonics 12 6 324 329 2018 10.1038/s41566-018-0153-0
[2] Thomas A. Tilting a ground-state reactivity landscape by vibrational strong coupling Science 363 6427 615 619 2019 10.1126/science.aau7742 30733414
[3] Dunkelberger A. D. Simpkins B. S. Vurgaftman I. Owrutsky J. C. Vibration-cavity polariton chemistry and dynamics Annu. Rev. Phys. Chem. 73 429 451 2022 10.1146/annurev-physchem-082620-014627 35081324
[4] Paravicini-Bagliani G. L. Magneto-transport controlled by Landau polariton states Nat. Phys. 15 2 186 190 2019 10.1038/s41567-018-0346-y
[5] Valmorra F. Vacuum-field-induced THz transport gap in a carbon nanotube quantum dot Nat. Commun. 12 1 5490 2021 10.1038/s41467-021-25733-x 34531384
[6] Cortese E. Excitons bound by photon exchange Nat. Phys. 17 1 31 35 2021 10.1038/s41567-020-0994-6
[7] Kockum A. F. Miranowicz A. Macrì V. Savasta S. Nori F. Deterministic quantum nonlinear optics with single atoms and virtual photons Phys. Rev. A 95 6 063849 2017 10.1103/PhysRevA.95.063849
[8] Garziano L. Ridolfo A. De Liberato S. Savasta S. Cavity QED in the ultrastrong coupling regime: photon bunching from the emission of individual dressed qubits ACS Photonics 4 9 2345 2351 2017 10.1021/acsphotonics.7b00635
[9] Kockum A. F. Miranowicz A. De Liberato S. Savasta S. Nori F. Ultrastrong coupling between light and matter Nat. Rev. Phys. 1 1 19 40 2019 10.1038/s42254-018-0006-2
[10] Ridolfo A. Savasta S. Hartmann M. J. Nonclassical radiation from thermal cavities in the ultrastrong coupling regime Phys. Rev. Lett. 110 16 163601 2013 10.1103/PhysRevLett.110.163601 23679600
[11] Fedortchenko S. Output squeezed radiation from dispersive ultrastrong light-matter coupling Phys. Rev. A 94 1 13821 2016 10.1103/PhysRevA.94.013821
[12] Carusotto I. Ciuti C. Quantum fluids of light Rev. Mod. Phys. 85 1 299 366 2013 10.1103/RevModPhys.85.299
[13] Forn-Díaz P. Lamata L. Rico E. Kono J. Solano E. Ultrastrong coupling regimes of light-matter interaction Rev. Mod. Phys. 91 2 025005 2019 10.1103/RevModPhys.91.025005
[14] Hagenmüller D. De Liberato S. Ciuti C. Ultrastrong coupling between a cavity resonator and the cyclotron transition of a two-dimensional electron gas in the case of an integer filling factor Phys. Rev. B 81 23 2010 235303 10.1103/PhysRevB.81.235303
[15] Scalari G. Ultrastrong coupling of the cyclotron transition of a 2D electron gas to a THz metamaterial Science 335 6074 1323 1326 2012 10.1126/science.1216022 22422976
[16] Maissen C. Ultrastrong coupling in the near field of complementary split-ring resonators Phys. Rev. B 90 20 205309 2014 10.1103/PhysRevB.90.205309
[17] Bayer A. Terahertz light-matter interaction beyond unity coupling strength Nano Lett. 17 10 6340 6344 2017 10.1021/acs.nanolett.7b03103 28937772
[18] Keller J. Few-electron ultrastrong light-matter coupling at 300 GHz with nanogap hybrid LC microcavities Nano Lett. 17 12 7410 7415 2017 10.1021/acs.nanolett.7b03228 29172537
[19] Mornhinweg J. Mode-multiplexing deep-strong light-matter coupling 2023 arXiv:2309.06915v1 .
[20] Todorov Y. Sirtori C. Few-electron ultrastrong light-matter coupling in a quantum LC circuit Phys. Rev. X 4 4 041031 2014 10.1103/PhysRevX.4.041031
[21] Rajabali S. An ultrastrongly coupled single terahertz meta-atom Nat. Commun. 13 1 2528 2022 10.1038/s41467-022-29974-2 35534466
[22] Gubbin C. R. Maier S. A. De Liberato S. Real-space Hopfield diagonalization of inhomogeneous dispersive media Phys. Rev. B 94 20 205301 2016 10.1103/PhysRevB.94.205301
[23] Khurgin J. B. Excitonic radius in the cavity polariton in the regime of very strong coupling Solid State Commun. 117 5 307 310 2001 10.1016/S0038-1098(00)00469-5
[24] Meiser D. Meystre P. Superstrong coupling regime of cavity quantum electrodynamics Phys. Rev. A 74 6 65801 2006 10.1103/PhysRevA.74.065801
[25] Załużny M. Zietkowski W. Intersubband cavity polaritons: the role of higher photonic modes Phys. Rev. B 80 24 245301 2009 10.1103/PhysRevB.80.245301
[26] Sundaresan N. M. Beyond strong coupling in a multimode cavity Phys. Rev. X 5 2 21035 2015 10.1103/PhysRevX.5.021035
[27] Brodbeck S. Experimental verification of the very strong coupling regime in a GaAs quantum well microcavity Phys. Rev. Lett. 119 2 27401 2017 10.1103/PhysRevLett.119.027401
[28] Moores B. A. Sletten L. R. Viennot J. J. Lehnert K. W. Cavity quantum acoustic device in the multimode strong coupling regime Phys. Rev. Lett. 120 22 227701 2018 10.1103/PhysRevLett.120.227701 29906138
[29] Balasubrahmaniyam M. Genet C. Schwartz T. Coupling and decoupling of polaritonic states in multimode cavities Phys. Rev. B 103 24 L241407 2021 10.1103/PhysRevB.103.L241407
[30] Rajabali S. Cortese E. Beck M. De Liberato S. Faist J. Scalari G. Polaritonic nonlocality in light–matter interaction Nat. Photonics 15 9 690 695 2021 10.1038/s41566-021-00854-3
[31] Cortese E. Mornhinweg J. Huber R. Lange C. De Liberato S. Real-space nanophotonic field manipulation using non-perturbative light–matter coupling Optica 10 1 11 2023 10.1364/OPTICA.473085
[32] Mornhinweg J. Halbhuber M. Ciuti C. Bougeard D. Huber R. Lange C. Tailored subcycle nonlinearities of ultrastrong light-matter coupling Phys. Rev. Lett. 126 17 177404 2021 10.1103/PhysRevLett.126.177404 33988443
[33] Halbhuber M. Non-adiabatic stripping of a cavity field from deep-strongly coupled electrons Nat. Photonics 14 11 675 679 2020 10.1038/s41566-020-0673-2 34221109
[34] De Liberato S. Ciuti C. Carusotto I. Quantum vacuum radiation spectra from a semiconductor microcavity with a time-modulated vacuum Rabi frequency Phys. Rev. Lett. 98 10 103602 2007 10.1103/PhysRevLett.98.103602 17358533
[35] Günter G. Sub-cycle switch-on of ultrastrong light-matter interaction Nature 458 7235 178 181 2009 10.1038/nature07838 19279631
[36] Garziano L. Ridolfo A. Stassi R. Di Stefano O. Savasta S. Switching on and off of ultrastrong light-matter interaction: photon statistics of quantum vacuum radiation Phys. Rev. A 88 6 63829 2013 10.1103/PhysRevA.88.063829
[37] Załużny M. Coupling between intersubband modes in wide quantum wells Phys. Status Solidi B 161 1 K19 K23 1990 10.1002/pssb.2221610149
