## Abstract

Ultrastrong coupling is a distinct regime of electromagnetic interaction that enables a rich variety of intriguing physical phenomena. Traditionally, this regime has been reached by coupling intersubband transitions of multiple quantum wells, superconducting artificial atoms, or two-dimensional electron gases to microcavity resonators. However, employing these platforms requires demanding experimental conditions such as cryogenic temperatures, strong magnetic fields, and high vacuum. Here, we use a plasmonic nanorod array positioned at the antinode of a resonant optical Fabry-Pérot microcavity to reach the ultrastrong coupling (USC) regime at ambient conditions and without the use of magnetic fields. From optical measurements we extract the value of the interaction strength over the transition energy as high as *g*/*ω* ~ 0.55, deep in the USC regime, while the nanorod array occupies only ∼4% of the cavity volume. Moreover, by comparing the resonant energies of the coupled and uncoupled systems, we indirectly observe up to ∼10% modification of the ground-state energy, which is a hallmark of USC. Our results suggest that plasmon-microcavity polaritons are a promising platform for room-temperature USC realizations in the optical and infrared ranges, and may lead to the long-sought direct visualization of the vacuum energy modification.

## Introduction

Two coupled harmonic oscillators is one of the most basic physical toy models that can be employed to understand the behavior of various mechanical and electromagnetic systems in simple intuitive terms. Usually, this approach is described by Newton’s equations of motion for the oscillators’ amplitudes *x*_{1} and *x*_{2}, wherein the coupling is mediated by the bi-linear interaction term (∝*gx*_{1}*x*_{2} or *g*(*x*_{1} − *x*_{2})^{2} depending on the system, where *g* is a coupling constant) in the system’s Hamiltonian^{1}. Such a mechanistic approach has been widely used to model the optical response of coupled plasmonic nanoparticles^{2}, exciton–polaritons^{3,4}, plasmon–excitons^{5,6}, and magnon–polaritons^{7,8}. While in the weak or strong coupling regime, *g* ≪ *ω*, the validity of this approach is accepted, it might not provide an adequate description of coupled electromagnetic systems when the coupling constant reaches a considerable fraction of the resonance energy, that is, *g* ~ *ω*.

Ultrastrong coupling (USC) is a regime of light–matter interaction in which the coupling strength, *g*, exceeds about 10% of the transition energy, *ω*^{9,10}. In this regime, the standard quantum optical approximations, such as the commonly made rotating wave approximation (RWA), fail. Thus so-called fast-rotating terms, as well as the quadratic *A*^{2} term must be taken into account in order to correctly describe the system’s behavior^{11,12,13}. The latter arises from the expansion of the minimal coupling Hamiltonian \(\left( {{\mathbf{p}} - \frac{e}{c}{\mathbf{A}}} \right)^2\) (where **p** and **A** are the particle’s momentum and the field’s vector potential, respectively) and is absent in the naive coupled oscillators model. Remarkably, not only quantum two-level systems, but also classical harmonic oscillators in the regime of ultrastrong coupling require description using the full Hamiltonians^{14}. One of the intriguing implications of USC is that the global vacuum energy of the system becomes dependent on the coupling strength^{15}; that is, if the coupling constant is allowed to vary, it will cost a certain amount of energy to adjust the value of the coupling strength. Additionally, the ground state gains a photonic component, that is, it contains a finite amount of virtual photon excitations^{15,16}. This in turn may lead to highly unusual phenomena, such as dynamical Casimir effect^{17,18,19} and single-photon frequency conversion^{20}.

Although the USC domain of light–matter interaction is of significant fundamental interest, it remains largely unexplored experimentally due to technical challenges of its realization. Indeed, so far the record-high realizations (where *g*/*ω* > 1) have been based on Landau polaritons^{21} and superconducting circuits^{22}, which require cryogenic temperatures and high magnetic fields. This specific interaction regime for which *g*/*ω* > 1 is called “deep” strong coupling. Replicating such results under ambient conditions remains a challenge. Room temperature realizations using collective coupling of organic molecules with microcavities have reached *g*/*ω* of “only” ∼0.3^{23,24}, with the recent implementation based on intersubband transitions of doped quantum wells showing *g*/*ω* ~ 0.7^{25}. However, such values are usually achieved by totally saturating the cavity volume with the material. Plasmonic lattices^{26,27} as well as single plasmonic nanorods^{28} have been shown to couple strongly with microcavity modes previously, however, the reported interaction strengths have not reached the level of the USC regime.

Here, we use our recently developed strategy based on plasmon–microcavity polaritons^{29} to achieve considerably higher coupling strengths, well into the USC regime, at room temperature. The plasmon–microcavity polaritons employed in this study consist of densely packed plasmonic nanorod arrays fabricated at the antinode of the Fabry–Pérot microcavity formed by two gold mirrors. By fitting the experimental reflection data by the spectrum of the full Hopfield Hamiltonian, we extract the normalized coupling strength *g*/*ω* as high as 0.55, one of the highest values for room temperature realizations of USC. We stress that such a high value is achieved here for a cavity, whose occupied volume amounts to only about 4%, with just a single layer of plasmonic nanoparticles. This makes a clear difference with respect to organic dyes^{23} and intersubband transitions^{25}, which reach similar numbers only by filling nearly 100% of the cavity interior. Furthermore, the experimental data allow us to indirectly observe the modification of the vacuum energy induced by USC, as well as estimate the photonic occupancy of the new ground state. In contrast to coupling bulk quantum wells or 2D electron gases, this vacuum energy effect can be potentially observed directly via the action of the vacuum energy landscape on a discrete nanoparticle. We thus argue that the large oscillator strength of plasmonic nanoparticle arrays, as well as control over their geometrical parameters and density, makes the plasmon–microcavity polaritons studied here an attractive platform for further investigations of room temperature ultrastrong and deep strong coupling regimes.

## Results

### Ultrastrong coupling in plasmon–microcavity systems

The system under study is illustrated in Fig. 1a. It consists of a sub-diffractive periodic array of gold (Au) nanorods placed at the antinode of the fundamental Fabry–Pérot (FP) microcavity mode formed by two Au mirrors and filled by a SiO_{2} spacer. The nanorod array couples to the vacuum field of the FP microcavity, thus producing plasmon–cavity polaritons manifested as distinct resonant spectral features emerging in transmission, reflection, and absorption spectra of the coupled system.

To provide initial insight into the behavior of the coupled system, we perform numerical finite-difference time-domain (FDTD) simulations (FDTD Solutions, Lumerical). Figure 1b shows a map of absorption spectra of coupled FP–nanorod systems at normal incidence with the electric field parallel to the nanowires as a function of the cavity thickness for nanorod lengths *L* = 300 nm and *d**y* = 30 nm spacing in the *y*-direction. For an easy comparison between these coupled system spectra with the uncoupled elements, we plot the bare FP cavity resonances with curved lines. The vertical dashed line marks the bare plasmon nanorod resonance of the array. A comparison clearly shows a rather complicated picture of new eigenmodes’ dispersion in which the even FP modes are practically unperturbed while the odd FP modes are shifted significantly from the bare cavity positions.

The 1st order FP mode of an empty cavity intersects the bare nanorod array plasmon resonance around 400 nm cavity thickness resulting in a distinct anticrossing (Fig. 1b). The lower polariton (LP) transitions from a plasmon-dominated mode (for a thin cavity) to an FP-dominated mode at large detuning (for a thick cavity). However, the upper polariton (UP) upon acquiring a plasmon-like character at large detuning, crosses the 2nd order FP mode and approaches the spectral position of the 3rd FP cavity mode, which in the coupled system is strongly pushed to the blue due to hybridization with the plasmon. Such qualitative behavior is observed for all the odd FP modes: each odd coupled *i*th mode is pushed to the blue beyond the subsequent even *i* + 1st mode (which is unperturbed) and approaches the *i* + 2nd odd FP mode. In fact, at no point in the spectral analysis do any of the FP–plasmon–polaritons follow the plasmon dispersion. In contrast, the even modes in the coupled system do not significantly interact with the array because they have a node of the electric field in the center of the cavity, where the rods are positioned. These observations suggest that a multimode character of the FP microcavity is important for a detailed interpretation of our results.

Another remarkable feature of the absorption map in Fig. 1b is the dispersion of the lower polariton in the thin cavity limit: for cavities thinner than about 200 nm, the LP dispersion exhibits a back-bending to extremely low energies, Fig. 1b, which is likely related to the near-field interaction of discrete plasmonic nanoparticles with the cavity mirrors in this short-range limit. This behavior is not reproduced by the Hamiltonian modeling, and we will not consider it in detail in the following.

The spatial distributions of the electric field induced by a normally incident plane wave inside the plasmon–cavity system calculated at the resonant energies for a 400 nm thick cavity (Fig. 1c) clearly display the opposite symmetries of the two resonances. While the lower energy mode shows an anti-symmetric combination of cavity and plasmon fields, featuring two saddle points above and below the nanorod, the upper energy mode is a symmetric combination. Such behavior highlights the polaritonic nature of the two resonances of the hybrid system. For a 400 nm thick cavity, corresponding to near-resonant coupling (*ω*_{cav} = *ω*_{pl} ∼ 0.8 eV), the Rabi splitting, Ω_{R}, estimated as the energy difference between the two absorption peaks reaches ~1 eV. Thus, assuming that Ω_{R} = 2*g* on resonance, we estimate the normalized coupling strength of *g*/*ω*_{pl} > 0.5, which clearly indicates the ultrastrong coupling regime in the system. In what follows, we perform a more rigorous estimation of the *g*/*ω*_{pl} values in our systems based on a full Hopfield Hamiltonian.

The same qualitative behavior is observed for diluted arrays and ones with longer nanorods, as shown in Supplementary Fig. 1. Although the coupling strength decreases with a smaller nanorod density, the spectral behavior indicates that all the odd modes mix with the plasmon, yielding a complex polaritonic system. However, for longer rods and/or diluted arrays, in addition to the plasmon, higher-order resonances are observed, which are associated with lattice modes. Additional data, including transmission and reflection spectra and complete data for the incident light polarization perpendicular to the nanorod axis, as well as the uncoupled cavity and array elements, are provided in Supplementary Figs. 1–5.

Samples of coupled plasmon–microcavity systems were fabricated by combination of electron beam evaporation (Au mirrors), plasma-enhanced chemical vapor deposition (dielectric spacers), and electron beam lithography (nanorod arrays) (see “Methods” for details). Figure 2a shows a bright-field optical microscope image of the fabricated nanorod arrays with lengths ranging from 200 to 400 nm with a step of 50 nm. The nanorods have fixed height of *h* = 20 nm and width of *w* = 50 nm. An exemplary scanning electron microscope (SEM) image of gold nanorods array with length of *L*_{rod} = 250 nm is shown in Fig. 2b (see “Methods” for details). Both the figures clearly show high-density plasmonic arrays with an interparticle distance as small as 30 nm, corresponding to the surface filling factor of 60%. More examples are shown in Supplementary Fig. 6.

Next, we proceed to optical measurements of the fabricated plasmon–cavity systems using Fourier transform infrared (FTIR) spectroscopy (see “Methods” for the details of measurements). Figure 3a, b shows exemplary, measured at normal incidence, reflection, and absorption spectra of an empty 400 nm thick cavity, 300 nm long nanorods array, and those of the coupled system (see Supplementary Figs. 7 and 8 for measured reflection and absorption spectra of all uncoupled cavities and nanorods). The uncoupled cavity and array resonances overlap spectrally and, when coupled, unambiguously confirm the realization of a giant Rabi splitting in the spectra of the coupled plasmon–cavity systems. We note that the 2nd order Fabry–Pérot mode redshifts in the hybrid system although it cannot interact with the nanorod array located exactly in the middle of the cavity. This behavior is consistently observed for all measured coupled systems, Supplementary Fig. 9. It could be explained by the hybrid cavities having slightly larger thickness due to the presence of the plasmonic array that could redshift all the uncoupled Fabry–Pérot modes.

Dispersion of measured normal-incidence reflection spectra from coupled systems with 300 nm long nanorods and varying cavity thickness displays a clear anticrossing between the 1st order Fabry–Pérot mode and the plasmon mode of the array, Fig. 3c (see Supplementary Fig. 9 for dispersions of reflection and absorption spectra vs cavity thickness for all nanorod lengths). The spectra also reveal the 2nd order Fabry–Pérot mode (third dip from the left), which does not interact with the nanorods due to the electric field node in the center of the cavity. As revealed by the Hamiltonian analysis of the spectra in the next section, the nanorods array mode additionally redshifts from ~0.8 to ~0.7 eV due to the presence of a dielectric medium around the rods. For the thinnest 100 nm thick plasmon-loaded cavities, neither reflection nor absorption spectra show the exact position of the upper polariton, which was beyond the detection range of the FTIR microscope used for this set of measurements. For this reason, we performed additional reflectivity measurements in the visible range for the 100 nm thick samples using a normal optical microscope to capture the spectral feature of the upper polariton (see Supplementary Fig. 10). Additional reflection spectra at normal incidence in the visible range were collected using a 20× objective (Nikon, NA = 0.45), directed to a fiber-coupled spectrometer and normalized with reflection from a standard dielectric-coated silver mirror. Based on these spectra, the vacuum Rabi splitting taken as the energy difference between the two reflection dips at zero detuning (*ω*_{cav} = *ω*_{pl}, 500 nm thick cavity), reaches ~0.8 eV at the resonant energy of ~0.7 eV, Fig. 3c. Thus, the Rabi splitting in our samples exceeds both the bare cavity and bare plasmon resonance frequencies, indicating that the hybrid plasmon–cavity system is deep into the USC regime.

### Analysis of the ultrastrong coupling using Hopfield Hamiltonian

We now turn to a more thorough analysis of the experimental data. Since a rough estimation already reveals that the Rabi splitting in our system is comparable to the transition energy of uncoupled oscillators, the usual Jaynes–Cummings or Rabi-type coupled Hamiltonians are invalid, and a more general Hamiltonian must be used. Therefore, to analyze our system we employ the full Hopfield Hamiltonian including both the fast-rotating and the quadratic *A*^{2} terms, which capture the essential physical characteristics of an ultrastrongly coupled system^{9}. We will focus on the two lowest modes of the plasmon–cavity structure, hence we will consider only coupling of two oscillators: the 1st order normal incidence FP mode of the cavity with energy \(\hbar \omega_{cav}\), and the collective long-axis plasmon mode of the array with energy \(\hbar \omega_{pl}\). Here, the cavity mode plays the role of the light component of the system, whereas the plasmonic nanorod array mode plays the role of the matter component. The total Hamiltonian thus reads:

where \(\hat a\) and \(\hat b\) are the microcavity and collective plasmon annihilation operators, respectively, and \(\hat H_{int}\) is the interaction Hamiltonian. If we were to consider individual nanoparticle plasmons interacting with each other instead of the collective array mode, the Hamiltonian would also yield additional eigenstates weakly interacting with light^{16}. As long as we work away from the Rayleigh modes of the array^{30}, which is ensured by sub-diffraction periodicity, all the plasmon–plasmon interaction effects can be absorbed into the collective plasmon frequency *ω*_{pl}.

The interaction part can be written differently depending on the gauge in which the electromagnetic field is treated. The two options that are often used are the Coulomb gauge and its dipole representation. The latter can be obtained from the Coulomb gauge by performing the Power–Zienau–Woolley transformation^{31}. When a cavity couples to a two-level system, the two representations are not gauge-invariant because of the two-level approximation^{32,33}. However, since we are considering coupling of two harmonic oscillators, the two pictures provide identical spectra^{14,16}. We will therefore use the Coulomb gauge, in which the single-mode interaction Hamiltonian can be written as^{34,35}:

where \(\hbar g_C = \mu _{pl}\sqrt {a^2\rho } {\cal{E}}_{vac}\frac{{\omega _{pl}}}{{\omega _{cav}}}\) is the coupling strength with *μ*_{pl} being the transition dipole moment of the plasmonic nanorod, *ρ* the plasmonic nanoparticles density per unit area *a*^{2} (\(\sqrt {a^2\rho }\) thus has a familiar \(\sqrt N\) scaling), and \({\cal{E}}_{vac} = \sqrt {\frac{{\hbar \omega _{cav}}}{{2\varepsilon \varepsilon _0a^2L_{\mathrm{eff}}}}}\) the vacuum electric field of the cavity with *L*_{eff} being the effective cavity mode transverse thickness^{16}. The first term in Eq. (2) is the usual Rabi-type interaction including both slow and fast-rotating terms. The second term is the so-called *A*^{2} term, which arises from the expansion of the minimal coupling Hamiltonian \(\left( {{\mathbf{p}} - \frac{e}{c}{\mathbf{A}}} \right)^2\) and “protects” the coupled system from the superradiant phase transition^{12,13}, as well as stabilizes the spectrum against the square-root singularity^{36}. Supplementary Figure 12 shows the spectrum of Hamiltonian (1) for *ω*_{pl} = *ω*_{cav} = 1 eV as a function of the coupling constant *g*_{C}; it also demonstrates that neglecting the *A*^{2} term as well as fast-rotating terms leads to incorrect and unphysical spectra.

We want to emphasize that although our system is essentially classical, we choose to use the quantum Hamiltonian because it provides a convenient description in terms of the modes amplitudes via creation operators from start. Moreover, the use of the quantum Hamiltonian allows us to obtain the characteristics of the ground state of the system in a straightforward way, which we analyze in the following. However, the linear response and the energy spectrum of the system can be equally obtained from a classical description not involving any operator algebra, which is demonstrated above by the FDTD simulations.

In a classical optical experiment, the outcome of which is some response function of the system, such as elastic scattering, reflection, or absorption, one cannot access directly the ground-state energy. However, spectral positions of the resonant features in reflection or absorption spectra reflect approximately the transition energies between the ground and first excited states of the system \(\hbar \omega _ \pm = E_{ \pm 1} - E_0\) . Therefore, to model the system with the Hopfield Hamiltonian framework, we fit the measured dispersions of reflection dips with calculated transition energies \(\hbar \omega _ \pm\) of the Hopfield Hamiltonian^{34}.

The resulting Hamiltonian fit of a coupled system’s resonant transitions as a function of the bare cavity energy is presented in Fig. 4a for *L*_{rod} = 300 nm nanorod arrays. For each cavity thickness, the bare cavity energy was determined from the spectral position of its reflection dip (Supplementary Fig. 7). By assuming that the effective cavity thickness scales as \(L_{\mathrm{eff}} = \frac{{\lambda _{cav}}}{{4n}}\) with *n* being refractive index of the cavity medium, we arrive at the coupling strength in the Coulomb gauge \(g_C = \omega _{pl}\mu _{pl}\sqrt {\frac{{\hbar \rho }}{{\pi \varepsilon _0nc}}}\), which is independent of the cavity's thickness and energy. Hence, we fit the polaritonic dispersion by freely varying plasmon frequency *ω*_{pl} and the coupling strength *g*_{C}. For the *L*_{rod} = 300 nm nanorod arrays, the fitting yields the plasmon frequency of 640 meV and the coupling strength of 300 meV, resulting in Rabi splitting of exactly 2*g*_{C} = 600 meV at resonance, *ω*_{pl} = *ω*_{cav} (see Supplementary Fig. 13 for Hamiltonian fits of other coupled systems, and Supplementary Table I for extracted plasmon energies and coupling strengths). For all five nanorod lengths, we consistently obtain normalized coupling strength values *g*_{C}/*ω*_{pl} in the range from 0.4 to 0.56, Fig. 4b, which unambiguously indicate the USC regime of interaction between the nanorods and the cavity modes^{9}. Furthermore, we notice that the normalized coupling strength \(g_C/\omega _{pl} = \mu _{pl}\sqrt {\frac{{\hbar \rho }}{{\pi \varepsilon _0nc}}}\) is a function of the plasmon transition dipole moment and the particle density only. Therefore, if the product \(\mu _{pl}\sqrt \rho\) grows with increasing nanorod length, we may expect even higher values of *g*_{C}/*ω*_{pl} for longer rods resonating at lower energies.

We also compare the resulting fits with those obtained by applying the multimode Hopfield Hamiltonian accounting for all the normal-incidence modes of a Fabry–Pérot cavity, which can be solved analytically^{35} (Supplementary Note 1). The most prominent difference is that in the multimode picture the upper polariton crosses the bare plasmon energy exactly at the point where it also crosses the even cavity mode (see Supplementary Fig. 14), in perfect agreement with FDTD simulations (Fig. 1b). Nevertheless, the resulting fits and coupling strengths are very close to the results of the single-mode Hamiltonian analysis (see Supplementary Table II), which justifies its validity in our case. Of course, the single-mode Hamiltonian only yields correct spectra as long as all higher-order (odd) cavity modes are detuned from the array mode, which is the case in our experimental configuration. The multimode Hamiltonian picture, however, may become important in other cases, as was also mentioned in the discussion of numerical FDTD results in Fig. 1b.

It is further instructive to compare the obtained values with an estimation for *g*_{C} that can be deduced directly from the geometry of the system. The vacuum electric field can be calculated as \({\cal{E}}_{vac} = \sqrt {\frac{{\hbar \omega _{cav}}}{{2\varepsilon \varepsilon _0a^2L_{\mathrm{eff}}}}}\) where for the effective cavity mode thickness one can use \(L_{\mathrm{eff}} = \frac{{\lambda _{cav}}}{{4n}}\). The plasmon transition dipole moment can be estimated from the scattering cross-section of a single nanorod in free space (see Supplementary Note 2) and by applying the classical Larmor formula for the decay rate of a point dipole^{37}. This yields the value of around 3.4 × 10^{4} Debye for the 400 nm long Au nanorod (corresponding to the radiative decay rate of ~67 meV, see Supplementary Table III). Combining these values with the nanorod density *ρ* = (430 nm⋅80 nm)^{−1}, one obtains the resonant (*ω*_{pl} = *ω*_{cav}) coupling strength of around *g*_{C} ≈ 0.3 eV, which agrees perfectly with the results of fitting.

We further illustrate the importance of keeping the quadratic term by analyzing the data with simpler, albeit a priori incorrect, Hamiltonians. An attempt to fit the experimental with eigenvalues of Hopfield Hamiltonian without the *A*^{2} term does not yield any reasonable result with a region of the energy spectrum becoming imaginary, Supplementary Fig. 17, and slightly overestimated coupling strengths. This imaginary spectrum is a fundamental property of the coupled oscillators Hamiltonian without any kind of quadratic stabilizing term^{36,38}. Fitting the data with no *A*^{2} Hopfield Hamiltonian under RWA (i.e., also without fast-rotating terms), although seems to give a better fit, yields regions with negative LP energy spectrum and largely overestimated coupling strength, Supplementary Fig. 17.

Besides dispersions of polaritonic energies, a remarkable feature of all measured reflection spectra is that the lower (and even upper) polaritons are much narrower than the bare plasmon mode (see, e.g., Fig. 3a and Supplementary Fig. 9). This behavior contrasts the non-Hermitian Jaynes–Cummings and Rabi Hamiltonians, where the decay rates of the photonic mode and electronic transition are “shared” equally at zero detuning between lower and upper polaritons: Im*ω*_{±} = −(*γ*_{cav} + *γ*_{x})/4, where *γ*_{x} is the linewidth of the electronic transition^{3,4}. A similar suppression of the polariton linewidth has been observed with a cyclotron resonance coupled to a Fabry–Pérot cavity^{39}. While we cannot directly introduce a non-Hermitian part into Hamiltonian (1), since it will render the ground-state energy complex-valued, we can qualitatively describe the decay processes in the coupled system. As we showed above, the linewidth of a single plasmonic nanorod outside the cavity is largely determined by its radiative loss. For nanorod arrays outside the cavity the radiative decay is even more dominant, reaching ~95% of the total linewidth (see Supplementary Note 2). However, when the dipolar oscillator, such as our plasmonic array, is placed between the mirrors, its radiation does not instantaneously leave the system—instead, it first bounces between the mirrors, and leaves the system only at the cavity’s leakage rate *γ*_{cav}. Therefore, fast radiative decay of the plasmon array mode becomes irrelevant, and the polariton linewidths are mostly determined by the total cavity’s and non-radiative array’s decay rates. This qualitative argument explains the apparent narrowness of the observed polaritonic bands. For a more rigorous description of the polaritons linewidth in the ultrastrongly coupled system (which is outside the scope of this work), the master equation approach might be needed^{40}.

### Ground-state energy and photonic occupancy

Having performed the fitting of the experimental data, we can analyze how the ground state of the system \(\left. {|G} \right\rangle\) is modified by the ultrastrong coupling. Again, we will restrict ourselves to the single-mode model of the system. In the uncoupled case, the global ground state is a direct product of the zero-photon and zero-plasmon states \({ { {|G} \rangle = |0_{cav}} \rangle \otimes |0_{pl}} \rangle\), and the energy of this state is \(E_G = \langle {0|H_{cav} + H_{pl}|0} \rangle = \frac{{\hbar} }{2}(\omega _{cav} + \omega _{pl})\), correspondingly. The USC modifies the global ground state \(\left. {|\tilde G} \right\rangle\) by admixing the states with different number of excitations, i.e., the global ground state with the higher excited states^{15}, thus modifying the ground-state energy. Since after diagonalization, the coupled system comprises two new harmonic oscillators, its ground-state energy is \(\tilde E_G = \frac{\hbar }{2}(\omega _ + + \omega _ - )\).

The *A*^{2} term modifies not only the energies of the polaritonic excited states^{41}, but also the ground-state energy of the system. Since the ground-state energy of a harmonic oscillator (or a set thereof) is half the transition energy (sum of those), its modification can be calculated as \(\delta E_G = \tilde E_G - E_G = \frac{\hbar }{2} (\omega _ + + \omega _ - - \omega _{cav} - \omega _{pl})\). By expanding the solution of Eq. (6) (see “Methods”) into a Taylor series near *g*_{C} = 0, the ground-state energy modification can be approximated by:

where we notated \(A=\omega_{cav}^2+\omega_{pl}^2\) and \(B=|\omega_{cav}^2-\omega_{pl}^2|\). At zero detuning (*ω*_{pl} = *ω*_{cav}) this expression yields \(\delta E_G = \frac{{g_C^2}}{{2\omega _{cav}}} + O\left( {g_C^4} \right)\).

The ground-state energy change at zero cavity–plasmon detuning can be estimated as \(\delta E_G \approx \frac{{g_C^2}}{{2\omega _{cav}}}\), which for *g*_{C}/*ω*_{pl} ≈ 0.5 becomes roughly *δE*_{G} ≈ *g*_{C}/4 ≈ 75 meV accounting for about 12% of the unperturbed ground-state energy *E*_{G}. This value of the relative energy modification is smaller than what could be obtained in the system studied in ref. ^{25} with \(\frac{{g_C}}{{\omega _0}} \approx 0.73\). However, our absolute value is much greater because our system exhibits interacting resonances in the near-IR to visible range with resonant energies around 1 eV, while the characteristic energies of the system in ref. ^{25} lie 1 order of magnitude lower at around 100 meV (additionally the plasmonic nanoparticle array in our case fills only about 4% of the cavity interior, as opposed to ref. ^{25} where the active material fully saturates the mode volume). Thus, the absolute ground-state energy change in our system is several times greater than *k*_{B}*T* at room temperature. This, in turn, implies that such ground-state energy modification might be important in practice and may show up in realistic USC-related effects even at room temperature.

The normalized ground-state energy variation \(\frac{{\delta E_G}}{{E_G}} = \frac{{\tilde E_G - E_G}}{{E_G}} = \frac{{\omega _ + + \omega _ - }}{{\omega _{cav} + \omega _{pl}}} - 1\) calculated using the obtained coupling strengths and analytical expressions for polariton energies *ω*_{±}, Fig. 5a, predicts up to ∼10% modification of the ground-state energy for normal incidence Fabry–Pérot mode upon coupling with the plasmonic array (see Supplementary Fig. 18 for the ground-state energy modification for other nanorod lengths). We want to emphasize that this vacuum energy is not an arbitrary reference level for all the higher energy states of the system. If the coupling constant can freely vary, for example, by allowing a single particle to move across a landscape with varying coupling strength, it will come to a state with the lowest vacuum energy even if the system is not coherently or thermally excited.

These theoretical values follow the experimentally obtained trend (circles), which was obtained using the measured cavity and polariton energies with only the bare plasmon frequency *ω*_{pl} adopted from the fitting. The theory predicts a relatively slow dependence of the normalized ground-state energy change on the detuning, whereas the experiment is more sensitive to that. This can be explained by the error in the determination of polaritons energies: in particular, the UP energy has been extracted with the error of up to ±0.1 eV, which already constitutes a few percent of the total ground state energy. Calculating the difference of two close values \(\tilde E_G - E_G\) makes the relative error even worse. An additional possible source of disagreement is that the true polariton energies do not exactly correspond to the extrema of a response function, such as transmission or reflection, but lie rather close to them due to the multi-mode nature of the system and the Fano resonance mechanism^{6}. Despite the non-ideal agreement, however, we stress that both theoretical predictions and experimental reflectivity data signal the ground-state energy modification of the order of 10% in our plasmon–microcavity systems. Such a modification is a clear hallmark of ultrastrong coupling, since in the conventional strong coupling picture, where *g* ≪ *ω*, the additive coupled and uncoupled energies are exactly the same, i.e., *ω*_{+} + *ω*_{−} = *ω*_{cav} + *ω*_{pl}, as can be seen from the Jaynes–Cummings model.

Lastly, we study the photonic occupancy \(\tilde n_{\mathrm{phot}} = \left\langle {\tilde G|\hat a^\dagger \hat a|\tilde G} \right\rangle\) of the modified ground state \(\left. {|\tilde G} \right\rangle\) (the plasmon occupancy of the ground state \(\left\langle {\tilde G|\hat b^\dagger \hat b|\tilde G} \right\rangle\) equals the photonic one^{15}). In the USC regime, the ground state of the system acquires a non-zero photonic component due to the aforementioned admixing of states with different excitation numbers^{15}. The photonic occupancy calculated with the use of the extracted coupling strength for the *L*_{rod} = 300 nm coupled systems, shown in Fig. 5b, suggests that the ground state of the ultrastrongly coupled system may contain up to 0.06 bare cavity photons for cavities resonant with the nanorod array (*ω*_{pl} = *ω*_{cav} ∼ 0.5 eV). This is smaller than 0.37 photons estimated for Landau polaritons in the THz range^{21}, but it is still a feasible number for converting to real photons by fast modulation of the coupling strength. The photonic occupancies calculated for other plasmonic nanorods predict almost identical values, Supplementary Fig. 19.

## Discussion

Above we have presented the ground-state modification taking into consideration only the normal incidence (*k*_{∥} = 0) mode of the cavity, whereas in reality all cavity modes having various in-plane momenta *k*_{∥} as well as TM and TE polarizations will couple to the nanorod array. Due to the periodicity of the system, modes with different *k*_{∥} do not interact and can be treated with independent Hamiltonians. The full vacuum energy per unit area of the cavity therefore can be calculated by integrating the vacuum energy over the entire *k*-space of the system. However, such an integration will diverge due to the asymptotic growth of the Fabry–Pérot modes energy at large *k*_{∥}. A regularization scheme will likely be needed to obtain a finite value similarly to the well-known result of Casimir^{42}; these calculations will be considered elsewhere.

Our plasmon–microcavity system offers a number of interesting perspectives. First, we have studied coupled systems with only one layer of plasmonic nanoparticles that occupies only ∼4% of the cavity interior. However, one can readily scale up the process and place several plasmonic layers in the center of the cavity close to the electric field anti-node. For example, placing four identical layers, assuming they all interact with the maximal electric field, will double the coupling strength and enable deep ultrastrong coupling with *g*_{C}/*ω*_{pl} > 1. One can also note from Fig. 5b that the extracted value of *g*_{C}/*ω*_{pl} monotonically increases with the nanorod length in the range of studied parameters. It is therefore interesting whether the normalized coupling strength can be further boosted by increasing the nanorods length, and at which rod length the maximal *g*_{C}/*ω*_{pl} ratio can be expected? As we showed previously, *g*_{C}/*ω*_{pl} in our system scales as \(\mu _{pl}\sqrt \rho\), which likely has an optimum. Furthermore, by precisely controlling the nanoparticles density, our system allows creating a vacuum energy gradient in the lateral direction. Lastly, the nanoparticles can be made chiral^{43}, opening the opportunities to create chiral vacuum states with various vacuum energies depending on the handedness of the chiral meta-atom.

To conclude, we have demonstrated a room-temperature ultrastrong coupling between two optical harmonic oscillators: a Fabry–Pérot microcavity and an array of plasmonic nanorods. The coupling strength reaches more than half of the cavity transition energy, thus unambiguously indicating the USC regime and one of the highest values for room-temperature implementations of *g*_{C}/*ω*_{pl} > 0.55. Importantly, this high value is achieved by filling only about 4% of the mode interior by plasmonic nanostructures, as opposed to alternative room temperature realization such as organic molecules and intersubband polaritons, which reach USC at 100% cavity filling factor. Analysis of the experimental data with the use of the Hopfield Hamiltonian reveals significant deviation of the coupled system’s eigenenergies from those predicted by the naive coupled oscillators model. Remarkably, the naive models fail to describe our system despite its obvious classical nature—both system’s components, plasmonic arrays, and Au mirrors contain millions of electrons and thus can be treated as classical harmonic oscillators. Furthermore, we indirectly observed a modification of the ground-state energy (up to 10%) and associated with that finite photonic occupancy induced by the ultrastrong coupling. We expect this vacuum energy change to be observable directly by the action of the cavity vacuum potential on a single nanoparticle, or via the dynamical Casimir effect. Our findings thus introduce a promising platform for studies of USC and related phenomena in the optical and infrared range at ambient conditions.

## Methods

### Samples fabrication

All samples were prepared on thin microscope glass (170 µm) coverslips. The glass coverslips were cleaned in acetone and isopropanol at 60 °C in ultrasonicator, dried with N_{2} blow, followed by oxygen plasma cleaning. Subsequently, 10 nm of gold (Au) mirror was prepared by e-beam evaporator with adhesion layer of chromium (2 nm) to form a bottom mirror. Then, various thicknesses of SiO_{2} layer for half-cavities were deposited by plasma-enhanced chemical vapor deposition (PECVD at 300 °C) on top of a freshly-prepared bottom gold mirror.

To fabricate a coupled system, lattice arrays of gold nanorods with various sizes and densities were fabricated on top of the half-cavities using a standard e-beam lithography. Then, the top-half SiO_{2} layers with the same thicknesses as the bottom SiO_{2} half cavities were deposited using PECVD. Finally, the coupled samples were completed by a deposition of 10 nm gold film as a top mirror for Fabry–Pérot cavity. Bare nanorod samples were prepared directly on top of glass substrates as a reference sample. To perform further SEM characterization, the samples were coated by a thin layer of conductive polymer (E-spacer). Morphology of the samples was characterized using a Zeiss (Germany) scanning electron microscope (SEM ULTRA 55 FEG).

### Optical measurements

Infrared optical measurements were performed with a Bruker Hyperion 2000 IR microscope (Schwarzschild-objective with 15× magnification, NA = 0.4) coupled to a Fourier-transform Bruker Vertex 80v spectrometer with a liquid-nitrogen-cooled mercury cadmium telluride detector. Reflection and transmission spectra were collected at normal incidence from a sample area of about 80 × 80 μm^{2} with 2 cm^{−1} resolution. All spectra were obtained with CaF_{2} IR polarizer in two principle orientations with the electric field polarization parallel and perpendicular to the nanorods long axis. A plane gold mirror was used as a reference in the reflection configuration experiment. Broad band absorption spectra were calculated from the measured reflection and transmission spectra. Reflection spectra in visible spectrum range were collected at normal incidence using a 20× magnification objective (Nikon, NA = 0.45), directed to a fiber-coupled spectrometer and normalized with reflection from a standard dielectric-coated silver mirror.

### FDTD simulations

FDTD simulations of the electromagnetic response of the coupled plasmon–cavity system were performed using commercial software (FDTD Solutions, Lumerical, Inc., Canada). Transmission and absorption spectra, as well as electromagnetic field distributions, were obtained with the use of a linearly polarized normally incident plane wave source and periodic boundary conditions with symmetries. The plane wave was polarized either along the nanorods or perpendicular to them. The permittivity of gold was approximated by interpolating the experimental data from Palik in the range 600–8000 nm. The simulation volume was discretized into a Δ*r* = 4 nm mesh with further refinement of 2 nm around the metal structures (nanorod and both mirrors).

### Hopfield Hamiltonian diagonalization

Spectrum of transition energies of Hamiltonian (1) with the interaction part (2) can be obtained as solutions of the following eigenproblem^{34}:

where \(\hat P = \alpha \hat a + \beta \hat b + \gamma \hat a^\dagger + \delta \hat b^\dagger\) is the polariton operator. Rewriting the eigenproblem in the basis of \(\hat a\), \(\hat b\), \(\hat a^\dagger\), and \(\hat b^\dagger\), solutions can be found as eigenvalues of the Hopfield matrix:

Two eigenvalues *ω*_{±} of the above matrix are given by the positive solutions of the bi-quadratic equation:

Thanks to the harmonicity of the coupled system, its entire energy ladder can be restored by collecting all possible values *ω*_{n,m} = *ω*_{0} + *nω*_{+} + *mω*_{−} where *n* and *m* are non-negative integers and \(\omega _0 = \frac{{\omega _ + + \omega _ - }}{2}\).

## Data availability

The data supporting the findings of this study, including experimental reflection and absorption spectra of samples, extracted polariton energies, the Mathematica code used for the data analysis, and simulated electric field plots can be downloaded at https://doi.org/10.5281/zenodo.3727173. Additional data are provided in Supplementary Information and are available from the authors upon reasonable request.

## References

- 1.
Landau, L. D. & Lifshitz, E. M.

*Mechanics*(Butterworth-Heinemann, 1976). - 2.
Prodan, E., Radloff, C., Halas, N. J. & Nordlander, P. A hybridization model for the plasmon response of complex nanostructures.

*Science***302**, 419–422 (2003). - 3.
Törma, P. & Barnes, W. L. Strong coupling between surface plasmon polaritons and emitters: a review.

*Rep. Prog. Phys.***78**, 13901 (2015). - 4.
Khitrova, G., Gibbs, H. M., Kira, M., Koch, S. W. & Scherer, A. Vacuum Rabi splitting in semiconductors.

*Nat. Phys.***2**, 81–90 (2006). - 5.
Wu, X., Gray, S. K. & Pelton, M. Quantum-dot-induced transparency in a nanoscale plasmonic resonator.

*Opt. Express***18**, 23633 (2010). - 6.
Pelton, M., Storm, S. D. & Leng, H. Strong coupling of emitters to single plasmonic nanoparticles: exciton-induced transparency and Rabi splitting.

*Nanoscale***11**, 14540–14552 (2019). - 7.
Zhang, X., Zou, C. L., Jiang, L. & Tang, H. X. Strongly coupled magnons and cavity microwave photons.

*Phys. Rev. Lett.***113**, 156401 (2014). - 8.
Zare Rameshti, B. & Bauer, G. E. W. Indirect coupling of magnons by cavity photons.

*Phys. Rev. B***97**, 014419 (2018). - 9.
Frisk Kockum, A. et al. Ultrastrong coupling between light and matter.

*Nat. Rev. Phys.***1**, 19–40 (2019). - 10.
Forn-Díaz, P., Lamata, L., Rico, E., Kono, J. & Solano, E. Ultrastrong coupling regimes of light–matter interaction.

*Rev. Mod. Phys.***91**, 025005 (2019). - 11.
Huppert, S., Vasanelli, A., Pegolotti, G., Todorov, Y. & Sirtori, C. Strong and ultrastrong coupling with free-space radiation.

*Phys. Rev. B***94**, 155418 (2016). - 12.
Nataf, P. & Ciuti, C. No-go theorem for superradiant quantum phase transitions in cavity QED and counter-example in circuit QED.

*Nat. Commun.***1**, 1069 (2010). - 13.
Rzaewski, K., Wódkiewicz, K. & Zakowicz, W. Phase transitions, two-level atoms, and the A2 term.

*Phys. Rev. Lett.***35**, 432–434 (1975). - 14.
Todorov, Y. & Sirtori, C. Few-electron ultrastrong light–matter coupling in a quantum LC circuit.

*Phys. Rev. X***4**, 041031 (2014). - 15.
Ciuti, C., Bastard, G. & Carusotto, I. Quantum vacuum properties of the intersubband cavity polariton field.

*Phys. Rev. B***72**, 115303 (2005). - 16.
Todorov, Y. & Sirtori, C. Intersubband polaritons in the electrical dipole gauge.

*Phys. Rev. B Condens. Matter Mater. Phys.***85**, 045304 (2012). - 17.
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**, 103602 (2007). - 18.
De Liberato, S., Gerace, D., Carusotto, I. & Ciuti, C. Extracavity quantum vacuum radiation from a single qubit.

*Phys. Rev. A Mol. Opt. Phys.***80**, 053810 (2009). - 19.
Stassi, R., Ridolfo, A., Di Stefano, O., Hartmann, M. J. & Savasta, S. Spontaneous conversion from virtual to real photons in the ultrastrong-coupling regime.

*Phys. Rev. Lett.***110**, 243601 (2013). - 20.
Sanchez-Burillo, E., Zueco, D., Garcia-Ripoll, J. J. & Martin-Moreno, L. Scattering in the ultrastrong regime: nonlinear optics with one photon.

*Phys. Rev. Lett.***113**, 263604 (2014). - 21.
Bayer, A. et al. Terahertz light–matter interaction beyond unity coupling strength.

*Nano Lett.***17**, 6340–6344 (2017). - 22.
Yoshihara, F. et al. Superconducting qubit-oscillator circuit beyond the ultrastrong-coupling regime.

*Nat. Phys.***13**, 44–47 (2017). - 23.
Gambino, S. et al. Exploring light–matter interaction phenomena under ultrastrong coupling regime.

*ACS Photonics***1**, 1042–1048 (2014). - 24.
Barachati, F. et al. Tunable third-harmonic generation from polaritons in the ultrastrong coupling regime.

*ACS Photonics***5**, 119–125 (2018). - 25.
Askenazi, B. et al. Ultra-strong light-matter coupling for designer Reststrahlen band.

*N. J. Phys.***16**, 043029 (2014). - 26.
Ameling, R. & Giessen, H. Cavity plasmonics: large normal mode splitting of electric and magnetic particle plasmons induced by a photonic microcavity.

*Nano Lett.***10**, 4394–4398 (2010). - 27.
Ameling, R. & Giessen, H. Microcavity plasmonics: strong coupling of photonic cavities and plasmons.

*Laser Photonics Rev.***7**, 141–169 (2013). - 28.
Konrad, A., Kern, A. M., Brecht, M. & Meixner, A. J. Strong and coherent coupling of a plasmonic nanoparticle to a subwavelength Fabry–Pérot resonator.

*Nano Lett.***15**, 4423–4428 (2015). - 29.
Bisht, A. et al. Collective strong light–matter coupling in hierarchical microcavity-plasmon-exciton systems.

*Nano Lett.***19**, 189–196 (2019). - 30.
De Abajo, F. J. G. Colloquium: Light scattering by particle and hole arrays.

*Rev. Mod. Phys.***79**, 1267–1290 (2007). - 31.
Woolley, R. G. Molecular quantum electrodynamics.

*Proc. R. Soc. A Math. Phys. Eng. Sci.***321**, 557–572 (1971). - 32.
De Bernardis, D., Pilar, P., Jaako, T., De Liberato, S. & Rabl, P. Breakdown of gauge invariance in ultrastrong-coupling cavity QED.

*Phys. Rev. A***98**, 053819 (2018). - 33.
Di Stefano, O. et al. Resolution of gauge ambiguities in ultrastrong-coupling cavity quantum electrodynamics.

*Nat. Phys.***15**, 803–808 (2019). - 34.
Hopfield, J. J. Theory of the contribution of excitons to the complex dielectric constant of crystals.

*Phys. Rev.***112**, 1555–1567 (1958). - 35.
De Liberato, S. Light–matter decoupling in the deep strong coupling regime: the breakdown of the purcell effect.

*Phys. Rev. Lett.***112**, 016401 (2014). - 36.
Bender, C. M., Felski, A., Hassanpour, N., Klevansky, S. P. & Beygi, A. Analytic structure of eigenvalues of coupled quantum systems.

*Phys. Scr.***92**, 015201 (2017). - 37.
Novotny, L. & Hecht, B.

*Principles of Nano-Optics*(Cambridge University Press, 2006). - 38.
Rokaj, V., Welakuh, D. M., Ruggenthaler, M. & Rubio, A. Light–matter interaction in the long-wavelength limit: no ground-state without dipole self-energy.

*J. Phys. B Mol. Opt. Phys.***51**, 034005 (2018). - 39.
Zhang, Q. et al. Collective non-perturbative coupling of 2D electrons with high-quality-factor terahertz cavity photons.

*Nat. Phys.***12**, 1005–1011 (2016). - 40.
Ciuti, C. & Carusotto, I. Input–output theory of cavities in the ultrastrong coupling regime: the case of time-independent cavity parameters.

*Phys. Rev. A Mol. Opt. Phys.***74**, 033811 (2006). - 41.
Anappara, A. A. et al. Signatures of the ultrastrong light–matter coupling regime.

*Phys. Rev. B***79**, 201303(R) (2009). - 42.
Casimir, H. B. G. & Polder, D. The influence of retardation on the London-van der Waals Forces.

*Phys. Rev.***73**, 360–372 (1948). - 43.
Baranov, D. G. et al. Circular dichroism mode splitting and bounds to its enhancement with cavity-plasmon-polaritons.

*Nanophotonics***9**, 283 (2020).

## Acknowledgements

D.G.B., B.M., A.C., B.R., G.J., and T.S. acknowledge financial support from Swedish Research Council (VR Miljö grant 2016-06059 and VR project grant 2017-04545). E.Z. acknowledges financial support from RSF (Optical measurements of the structures were performed with the financial support of the Grant No. 17-79-20418). T.J.A. acknowledges support from the Polish National Science Center via the project 2017/25/B/ST3/00744. Open access funding provided by Chalmers University of Technology.

## Author information

### Affiliations

### Contributions

D.G.B. performed Hopfield Hamiltonian analysis. B.M. prepared the samples. E.Z. and B.M. performed FTIR and optical reflectivity measurements. A.B. and A.C. prepared initial samples and performed initial simulations. T.J.A. and D.G.B. performed numerical FDTD simulations. D.G.B. and B.R. analyzed the data. D.G.B., T.J.A., and T.S. wrote the manuscript with input from all co-authors. G.J. and T.S. supervised the project. All authors discussed the results.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

**Peer review information** *Nature Communications* thanks Michael Rugentaller and Salvatore Savasta for their contribution to the peer review of this work. Peer reviewer reports are available.

**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

Baranov, D.G., Munkhbat, B., Zhukova, E. *et al.* Ultrastrong coupling between nanoparticle plasmons and cavity photons at ambient conditions.
*Nat Commun* **11, **2715 (2020). https://doi.org/10.1038/s41467-020-16524-x

Received:

Accepted:

Published:

## 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.