# Strong Coupling between Localized Surface Plasmons and Molecules by Coupled Cluster Theory

Jacopo Fregoni, Tor S. Haugland, Silvio Pipolo, Tommaso Giovannini, Henrik

Koch, Stefano Corni

► To cite this version:

Jacopo Fregoni, Tor S. Haugland, Silvio Pipolo, Tommaso Giovannini, Henrik Koch, et al. Strong Coupling between Localized Surface Plasmons and Molecules by Coupled Cluster Theory. Nano Letters, 2021, 21 (15), pp.6664-6670. 10.1021/acs.nanolett.1c02162. hal-04134576

HAL Id: hal-04134576

https://hal.univ-lille.fr/hal-04134576

Submitted on 20 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Strong Coupling between Localized Surface Plasmons and Molecules by Coupled Cluster Theory

Jacopo Fregoni,⁷ Tor S. Haugland,⁷ Silvio Pipolo, Tommaso Giovannini, Henrik Koch,\* and Stefano Corni\*

<figure><img src="image_3.png" /></figure>

Cite This: *Nano Lett.* 2021, 21, 6664–6670

<figure><img src="image_4.png" /></figure>

Read Online

## ACCESS

<figure><img src="image_5.png" /></figure>

Metrics & More

<figure><img src="image_6.png" /></figure>

Article Recommendations

<figure><img src="image_7.png" /></figure>

Supporting Information

**ABSTRACT:** Plasmonic nanocavities enable the confinement of molecules and electromagnetic fields within nanometric volumes. As a consequence, the molecules experience a remarkably strong interaction with the electromagnetic field to such an extent that the quantum states of the system become hybrids between light and matter: polaritons. Here, we present a nonperturbative method to simulate the emerging properties of such polaritons: it combines a high-level quantum chemical description of the molecule with a quantized description of the localized surface plasmons in the nanocavity. We apply the method to molecules of realistic complexity in a typical plasmonic nanocavity, featuring also a subnanometric asperity (picocavity). Our results disclose the effects of the mutual polarization and correlation of plasmons and molecular excitations, disregarded so far. They also quantify to what extent the molecular charge density can be manipulated by nanocavities and stand as benchmarks to guide the development of methods for molecular polaritonics.

<figure><img src="image_8.png" /></figure>

**KEYWORDS:** Plexcitons, Polaritonic Chemistry, Cavity-QED, Quantum Nanoparticles, Quantum Chemistry, Nanoplasmonics, Quantum coupling

Strong coupling between molecules and quantum plasmons¹ in nanocavities² leads to the formation of hybrid plasmon–molecule states: polaritons, or more specifically, plexcitons.³ These new states manifest distinct features compared to the original states,⁴⁻⁷ potentially resulting in modified chemical/photochemical reactivity⁸⁻¹⁰ and relaxation dynamics,¹¹,¹² along with other coherent processes.¹³,¹⁴ Modeling accurately the molecules, nanostructures and their coupling is of the utmost importance to support the experimental advances. Coupling descriptions proper for resonant optical cavities, such as using the molecular dipole only,¹⁵,¹⁶ should be amended¹⁷⁻¹⁹ when the coupling affects the molecule on a submolecular level.¹⁸,²⁰,²¹

There are currently no approaches that treat the coupled plasmon–molecule system nonperturbatively, that is, that include the relaxation of the molecular electronic density upon polariton formation together with plasmon–molecule correlation. In this work we extend the QED-CC method,²²,²³ already applied to resonant optical cavities, to realistic nanoplasmonic cavities. We quantize the plasmonic excitations starting from the classical dielectric description of nanostructures with arbitrary shapes. Each quantized plasmonic mode is associated with a surface charge density (discretized into point charges), analogous of a transition charge density in a molecular system. A similar quantization approach for

plasmons in Drude metals, derived from macroscopic QED,¹⁹,²⁴,²⁵ was adopted in ref 17. The quantization framework we are presenting accounts for the geometrical features of the nanoparticles setups without disregarding the molecular complexity. Indeed, coupled cluster (singles and doubles excitations) is recognized as a highly accurate method in quantum chemistry:²⁶ we extend it to include self-consistent and correlated molecule–plasmon hybridization. With this method, we compute the interaction between a nanocavity realized after ref 18 and two realistic molecules: porphyrin and para-nitroaniline (PNA). For porphyrin, we highlight the role of the geometrical features of the system and the role of electron–plasmon self-consistency and correlation in polaritons, beyond the standard Jaynes-Cumming picture.²⁷ We use PNA to quantify the same effects, already predicted on the ground and excited states electron densities²² for a resonant cavity, in the case of a nanocavity.

Received: June 2, 2021

Revised: July 9, 2021

Published: July 20, 2021

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>
<figure><img src="image_11.png" /><figcaption>Figure 1. Examples of plasmonic modes computed by Q-PCM-NP. (a) Example of plasmonic dipolar mode for a 10 nm sphere. The p index associated with each set of charges does not take into account the 2l + 1 degeneracy for the modes of the sphere. Hence, the dipolar modes are p = 2 to p = 4. (b) Example of a quadrupolar mode for the same sphere. Similar as for the dipolar case, the quadrupolar modes are 2l + 1 degenerate, hence they are from p = 5 to p = 9 and so on. (c) Relevant plasmonic mode for a nanotip on a support. The mode is characterized by a strong charge density concentrated on the tip as in ref 18. The dimensions of the system are about 10 × 10 nm (see Figure S1b)</figcaption></figure>

## ■ RESULTS AND DISCUSSION

**Theory.** The first subsystem we focus on is the nanostructure. We start from a classical dielectric perspective, namely computing the plasmon modes from a continuous medium model. As such, it is not an atomistic description and does not include dielectric nonlocal effects, electron spillout, or chemical bonds. However, the description we adopt is surprisingly accurate:²⁸ it reproduces the inhomogeneous electric field due to atomistic irregularities in the nanostructures²⁹ even at subnanometric level.³⁰ Our goal is to describe the nanostructure at a quantum level as a material system, which differs from other approaches that quantize the electromagnetic field in the presence of dielectric media.²⁴,³¹ By building on polarizable continuum model applied to nanoparticles (PCM-NP) approach,³² we shall achieve quantization (Q-PCM-NP) in the quasi-electrostatic framework, where the retardation effects are disregarded.³³ Similar quantization frameworks for the plasmonic modes were used in different contexts not involving a molecule³⁴ or using a model description of the molecule.³⁵ To classically compute the nanoparticle dielectric response properties, we solve Maxwell's equations for a continuous, frequency-dependent dielectric (nanostructure) under an external electromagnetic perturbation. In the quasi-static limit, the nanoparticle experiences an electrostatic potential at a given frequency, V(ω), which induces polarization charges q(ω) on the nanoparticle surface. The equation defining such response charges is

$$q(\omega) = Q^{\text{IEF}}(\omega)V(\omega) \quad (1)$$

where Q^IEF(ω) is the frequency-dependent response to the external perturbation V(ω). We refer to IEF as to the specific integral equation formalism³⁶ formulation of the PCM problem. In the present work, V(ω) is the potential produced by the molecular transition density. The frequency-dependent linear charge density response function of the nanostructure is taken in a diagonal form³⁷

$$Q^{\text{IEF}}(\omega) = -S^{-1/2}TK(\omega)T^*S^{-1/2} \quad (2)$$

$$K_p(\omega) = \frac{2\pi + \Lambda_p}{2\pi^{\frac{\epsilon(\omega)+1}{\epsilon(\omega)-1}} + \Lambda_p} \quad (3)$$

where K is the diagonal linear-charge-response-matrix derived from eigenvalues Λp of the appropriate IEF matrix,³⁷ ε(ω) is the frequency-dependent dielectric function, and S is the matrix storing the electrostatic potential between discrete points of the dielectric. This is the starting point of the Q-PCM-NP quantization procedure.

**Quantization of the Plasmonic Modes: Q-PCM-NP.** We assume that the metal nanoparticle can be characterized by a Drude-Lorentz dielectric function

$$\epsilon(\omega) = 1 + \frac{\Omega_p^2}{\omega_0^2 - \omega^2 - i\gamma\omega} \quad (4)$$

where Ωp² is the squared plasma frequency of the bulk metal, ω₀ is the natural frequency of the bound oscillators (Lorentz model), and γ is the damping rate. By defining

$\omega_p^2 = \omega_0^2 + \left(1 + \frac{\Lambda_p}{2\pi}\right) \frac{\Omega_p^2}{2}$ (neglecting the second-order term in γ), we retrieve (see Supporting Information) the full form of the response function Q^IEF(ω) from eqs 3 and 4

$$Q_{kj}^{\text{IEF}}(\omega) = -\sum_{p} (S^{-1/2}T)_{k,p} \sqrt{\frac{\omega_p^2 - \omega_0^2}{2\omega_p}} \left( \frac{1}{\omega_p + \omega + i\frac{\gamma}{2}} + \frac{1}{\omega_p - \omega - i\frac{\gamma}{2}} \right) \sqrt{\frac{\omega_p^2 - \omega_0^2}{2\omega_p}} (T^*S^{-1/2})_{p,j} \quad (5)$$

Each matrix element of Q^IEF(ω) is evaluated on representative points of the tesserae, k and j. In eq 5, we recognize the spectral form of a linear response function,³⁸ Q^quant(ω). Our goal now is to identify the quantities characterizing the excited states of the nanostructure (that is, plasmons) relevant to devise the coupling with the molecule. For the quantum description of the nanostructure, the response function Q^quant(ω) is formally written in terms of the surface charge operator ˆq.³⁹ We conveniently write Q^quant(ω) matrix elements in its spectral representation⁴⁰ as

$$Q_{kj}^{\text{quant}}(\omega) = -\sum_{p} \left( \frac{\langle 0l\hat{q}_k|p \rangle \langle pl\hat{q}_j|0 \rangle}{\omega_p + \omega + i\gamma_{p'}} + \frac{\langle pl\hat{q}_k|0 \rangle \langle 0l\hat{q}_j|p \rangle}{\omega_p - \omega - i\gamma_{p'}} \right) \quad (6)$$
<figure><img src="image_12.png" /><figcaption>Figure 2. Interaction between the free-base porphyrin and a nanotip. (a) Setup for porphyrin interacting with a plasmonic nanotip used in panels b and c, where the nanotip has transition energy $\omega = 4.06$ eV. (b) Map of the plasmon-molecule coupling g for two excitations in porphyrin (4.05 eV, 4.07 eV) and the effective coupling $g_{eff} = \sqrt{g_1^2 + g_2^2}$. Axis x and axis y indicate the displacement of the nanotip from the center of the molecule. (c) Oscillator strengths (intensity) for porphyrin displaced 6 Å along x and 6 Å along y (see black square in b)). Intensities are relative to porphyrin's strongest transition. The inset shows the same plot with only porphyrin. JC(g_{eff}) refers to a two state JC model where the coupling and intensities refer to g_{eff}.</figcaption></figure>

By comparing eq 5 and eq 6, we verify that $\omega_p$ represents the excitation frequencies of the plasmonic system. $\gamma_p' = \gamma/2$ for each plasmon state p is a decay rate, while we can identify $\langle 0|\hat{q}_k|p \rangle$ with

$$ \langle 0|\hat{q}_k|p \rangle = (S^{-1/2}T)_{k,p} \sqrt{\frac{\omega_p^2 - \omega_0^2}{2\omega_p}} = q_{p,k} \quad (7) $$

The label $q_{p,k}$ represents the transition charge sitting on the kth tessera associated with the mode p and they are the analogous of transition densities in molecules. Expressed differently, the ensemble of the $q_{p,k}$ set of charges represents the normal mode of the plasmonic system at frequency $\omega_p$ (see Figure 1a-c as an example).

We are now in the position to write the quantum plasmonic Hamiltonian as

$$ \hat{H}_p = \sum_p \omega_p (\hat{b}_p^\dagger \hat{b}_p) \quad (8) $$

where $\omega_p$ is the p-mode frequency and $\hat{b}_p^\dagger, \hat{b}_p$ are the corresponding bosonic creation and annihilation operators.17,25 Assuming the localized surface plasmons to be bosons is not a result of the present derivation (at the linear response level, a Fermionic or bosonic excitation would yield the same $Q^{quant}(\omega)$). It is rather based on the general behavior of plasmons in extended systems.1 In the Supporting Information, we report numerical checks of the present approach against analytical results for a point-dipole molecule interacting with a Drude metal nanosphere.41 We can now write the general quantized plasmon-molecule Hamiltonian as

$$ \hat{H} = \hat{H}_{\text{mol}} + \hat{H}_{p} + \hat{H}_{\text{int}} \quad (9) $$

where $\hat{H}_{mol}$ is the standard molecular Hamiltonian (which we reduce to the electronic Hamiltonian). Making use of the response charges for the quantized plasmon, the interaction term $\hat{H}_{int}$ reads39

$$ \hat{H}_{\text{int}} = \sum_{p,j} q_{p,j} \hat{V}_j (\hat{b}_p^\dagger + \hat{b}_p) \quad (10) $$

where $\hat{V}_j$ is the molecular electrostatic potential operator, evaluated at each tessera representative point j. Our formulation of $\hat{H}_{int}$ presents two major advances with respect to simplified models (see Supporting Information for the formal comparison): we obtain an intuitive representation of the inhomogeneous electromagnetic field in terms of polarization charges and the interaction is general on the molecular side, meaning that it can be interfaced to advanced electronic structure models. Through this interface, our method can describe phenomena occurring in both the weak and strong coupling regime.42,43 Radiative corrections to the quasi-static description, related to radiation reaction,44 can also be included in this model, at the quantum level. The plasmon radiative decay rate can be calculated from the transition dipole moment associated with the charges $q_{p,j}$ and it can then be added to the nonradiative lifetime $\gamma_p$. For the systems treated in the numerical section, this correction is negligible. An estimate of the molecule-plasmon coupling $g_n$ between the plasmonic mode p and the molecular electronic transition from the ground state $S_0$ to the excited state $S_n$ is given by
<figure><img src="image_13.png" /><figcaption>Figure 3. Quantum coupling between a PNA molecule and a nanotip. (a) Sketch of the molecule interacting with the plasmonic nanotip. The parallel orientations identify the maximum plasmon–molecule coupling condition, whereas the perpendicular one corresponds to the minimum coupling. (b) PNA oscillator strength (intensity) at different orientations inside a plasmonic nanotip-cavity. Each calculation is performed at different molecular orientations around the nuclear center of charge. The plasmon is resonant with the bright molecular transition (4.52 eV). The comparison between the realistic structure of the molecule (QED-CCSD) and the two-level JC model is shown by the full colored lines. (c) Plasmon-induced difference in PNA’s ground-state density, when PNA is parallel with the plasmonic mode of frequency ω = 4.52 eV. Black and white isosurfaces indicate increased and reduced density (±10⁻³e⁻/a₀³).</figcaption></figure>

$$g_n = \langle S_0, 1 | \hat{H}_{\text{int}} | S_n, 0 \rangle = \sum_j q_{p,j} V_j^{(0,n)} \quad (11)$$

Here 0 and 1 are the occupation numbers of the plasmonic mode p. The coupling $g_n$ is the central quantity in simplified approaches to strong coupling, such as the Jaynes-Cummings model (JC).

**Extending QED Coupled Cluster (QED-CC) to Plasmon–Molecule Systems.** Determining the eigenfunctions of the plasmon–molecule Hamiltonian in eq 9 requires an accurate description of the molecule. In this section, we outline the extension of the QED-CC treatment (already applied to resonant cavities²²,²³) to the case of plasmons. The detailed algorithm is available in the QED-CC dedicated section of the Supporting Information. The simplest description of the molecule without plasmon interactions is a single Slater determinant, |HF⟩ (Hartree–Fock). To include electron correlation we rely on coupled cluster theory.²⁶ More explicitly, we apply the exponential of the cluster operator to the |HF⟩ determinant

$$|\text{CC}\rangle = e^{\hat{T}} |\text{HF}\rangle \quad (12)$$

The cluster operator $\hat{T}$ is expressed as

$$\hat{T} = \hat{T}_1 + \hat{T}_2 + \dots + \hat{T}_{N_c} \quad (13)$$

where $T_1$ are linear combinations of single excitations, $T_2$ are linear combinations of double excitations and so on. Inspired by the formal similarities of the interaction in molecule–plasmon and molecule–photon systems, we extend the newly developed quantum electrodynamics coupled cluster theory (QED-CC)²² from photons to plasmons. Here the plasmon–molecule interaction is described nonperturbatively, allowing also the ground state to couple to the electromagnetic field.²³ In QED-CC, the cluster operator $\hat{T}$ also includes the bosonic creation and annihilation operators $\hat{b}_p^\dagger, \hat{b}_p$ of the quantized plasmon mode p

$$\hat{T} = \hat{T}_1 + \hat{T}_2 + \hat{S}_1^1 + \hat{S}_2^1 + \hat{\Gamma}^1 \quad (14)$$

Here we identify three contributions: electronic ($\hat{T}_1, \hat{T}_2$), plasmonic ($\hat{T}^1$) and mixed plasmon–molecule excitations ($\hat{S}_1^1, \hat{S}_2^1$). By the effect of such excitations, QED-CCSD-1 (SD denotes single and double excitation of the electronic subspace, whereas 1 is the excitation order of the plasmonic modes) includes the correlation induced on the electronic states by the plasmonic transition. We note that, in the limit where the interaction $\hat{g} \to 0$, CCSD and QED-CCSD-1 are equivalent.

**Calculations with Molecules. Free-Base Porphyrin.** In this section, we exploit Q-PCM-NP coupled to QED-CCSD-1 and EOM-CCSD to simulate plasmon–molecule interactions. Here we consider a free-base porphyrin coupled to the plasmonic nanotip, as sketched in Figure 2a. A similar setup was adopted in refs 17 and 18. According to eq 11, we define the couplings between the ground state $|S_0, 1\rangle$ of the free-base porphyrin and the $|S_1, 0\rangle, |S_2, 0\rangle$ states as $g_1$ and $g_2$ respectively. In Figure 2b, we show the associated interaction maps, where the couplings $g_1$ and $g_2$ are displayed as a function of the molecule–nanotip displacement. The two quasi-degenerate dipolar transitions of free-base porphyrin at 4.05 eV ($S_1$) and 4.07 eV ($S_2$) have different symmetries. Depending on the tip displacement, the interaction element g varies between 0 meV at the center to 36 meV at 6 Å displacement. By displacing the tip accordingly, a polariton is created with either of the two states. When combining the two quasi-degenerate transitions, $g_{\text{eff}} = \sqrt{g_1^2 + g_2^2}$, a doughnut-like shape of the transition density is obtained, in agreement with previous results on zinc phtalocyanine.¹⁷ We stress that this doughnut-like shape is not observable without including the geometry of the nanostructure. In Figure 2c, we show the absorption spectrum of the relaxed tip–molecule complex (QED-CCSD-1) and compare it with the polaritons obtained by the Jaynes-Cummings model (see the Supporting Information for details). The point at which we perform the calculations is highlighted as the black
<figure><img src="image_14.png" /><figcaption>**Figure 4.** Charge displacement analysis of PNA. The molecule is placed parallel with respect to the plasmonic nanotip. The direction of charge transfer (CT) is along the arrow. (a) Plasmon-induced density difference in the ground state. (b) Charge transfer excitation around 4.5 eV. The red and orange areas are upper (UP) and lower (LP) polaritons, respectively. The blue curve is the sum of red and orange.</figcaption></figure>

square in panel *g*<sub>eff</sub> of Figure 2b. In QED-CCSD-1, the light–matter system is allowed to correlate both the ground and the excited states, guaranteeing accurate results with respect to the truncated basis. Even more, the interactions between all singly and doubly excited determinants with and without plasmons are included. Conversely, the JC case is treated by correlating the electronic states only, hence neglecting the mixed molecule–plasmon correlation contributions. As a further difference between JC and QED-CCSD-1, the interaction in the former is a simple norm of the *g*<sub>1</sub> and *g*<sub>2</sub> interactions, and it does not allow for the *S*<sub>1</sub> and *S*<sub>2</sub> states to interact through the plasmon. For QED-CCSD-1, the presence of the polaritonic correlation and several electronic states reveals that the plasmon contribution is distributed over other electronic states displayed as the small peak at 3.88 eV in Figure 2c. At the same time, the Rabi splitting between the two main polaritonic peaks is reduced (46 meV) with respect to JC (70 meV). While the small peak is due to the inclusion of different states, the difference in the Rabi splitting is due to the inclusion of the mixed plasmon–molecule correlation. To corroborate this view, we compare in the Supporting Information an extended JC model using three states with the QED-CCSD-1. Although the description improves with the extended JC model, it still misses the intensities and the positions of the polaritonic peaks, confirming the role of the mixed plasmon–molecule correlation in the description.

**Para-nitroaniline.** Para-nitroaniline (PNA) presents an intense absorption peak in the UV range. We hence expect large Rabi splittings when the PNA is appropriately aligned with the electromagnetic field of the nanotip. The transition density associated with the bright transition of PNA at ~4.5 eV is parallel to the long axis of the molecule, making it an excellent candidate for single-molecule strong coupling.

In Figure 3, we investigate the coupling of the PNA molecule with the nanotip. We identify the maximum (parallel) and the minimum (perpendicular) coupling conditions, as sketched in Figure 3a, and plot the Rabi splitting of the plasmon–molecule coupled system in Figure 3b. When the molecular transition dipole is perpendicular to the plasmon mode, there is little-to-no effect on the spectrum. Instead, when the molecule is gradually rotated to be parallel to the mode, the Rabi splitting increases until the maximum of 179 meV (parallel). By gradually rotating the molecule in the opposite direction (antiparallel) we observe a slightly smaller

splitting (175 meV). The 4 meV difference between the parallel and antiparallel orientations is a signature of the realistic molecular description. Indeed, a dipolar interaction would result in a symmetric spectrum with a Rabi splitting independent of the parallel/antiparallel orientation. In addition, we observe a red-shift of the polaritonic peaks when including the mixed plasmon–molecule correlation energy via QED-CCSD-1.

The Jaynes-Cummings model fails at describing ground-state properties, as no change in the molecular ground state is predicted by the model. By instead using QED-CCSD-1, we show in Figure 3c how the ground-state electronic density changes due to the molecule–plasmon interaction. Since the coupling in the present case is smaller than the one considered in ref 22 the charge localization effects as discussed in the same reference are reduced. We quantify this effect in Figure 4a. We note the opposite trend to what was described in ref 22, namely that here the plasmon induces a small charge transfer from the donor to the acceptor. This inversion is a signature that the molecules interact differently with plasmonic nanotips compared to optical cavities, whereas the Jaynes-Cummings model treats them equivalently. The charge transfer properties of the polariton states are also substantially changed, reminiscent of what is shown in ref 22. However, in the present case the sum of the polaritons’ charge transfer adds up to the charge transfer of the singlet state without the plasmons, as displayed in Figure 4b. This again highlights that interactions with nanoparticle plasmons and optical cavities are not equivalent. The plasmon-induced changes in the electronic density also imply a modification of other ground-state properties as well as the potential energy surfaces. One example is the minor induced change of the dipole moment (−0.03 D). One important limitation in this ground-state study is the inclusion of only a single plasmon mode: higher energy modes are expected to contribute to the ground-state coupling and enhance the minor effects seen here, as the ground-state correlation is not a resonant property.<sup>23</sup>

## CONCLUSIONS

In this work, we have introduced a new methodology, based on PCM-NP<sup>45,46</sup> and QED-CC,<sup>22</sup> to treat the quantum coupling between nanoparticles of arbitrary shapes and molecules. We compared our description of polaritons in free-base porphyrin with the simplified Jaynes-Cummings model, highlighting the
role of the mixed molecule–plasmon correlation. Finally, we
have analyzed the effects on the ground state of PNA in a
realistic plasmonic nanocavity setup, showing that the
modification of ground states density previously investigated
by QED-CC is still present. Our representation of quantum
plasmonic modes as apparent surface charges, in synergy with
coupled cluster, provides a simple-yet-accurate interface to
tackle mixed nanoparticle–molecules systems. In conclusion,
we believe our method provides solid ground to support and
guide experiments investigating and manipulating the chemical
properties at the submolecular level. Accurate calculations of
plexciton-affected reaction-barriers are already possible, as well
as exploring effects related to chiral nanotips. The method is
also prone to various extensions on the nanostructure side: the
inclusion of multiple EM modes or the extension to classical
atomistic models.47–49 Finally, the model may assess the
properties of molecules within plasmonic nanocavities in turn
hosted by optical cavities, a setup whose experimental
exploration has been recently started.50

■ ASSOCIATED CONTENT

● Supporting Information

The Supporting Information is available free of charge at
https://pubs.acs.org/doi/10.1021/acs.nanolett.1c02162.

Details of the theoretical and computational specifics for Q-PCM-NP and QED-CC and further benchmarks and supplementary results; computational details for the reproducibility of the calculations and the codes availability (PDF)

■ AUTHOR INFORMATION

Corresponding Authors

Henrik Koch — *Scuola Normale Superiore, I-56126 Pisa, Italy*; *Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway*;
*orcid.org/0000-0002-8367-8727*; Email: *henrik.koch@sns.it*

Stefano Corni — Dipartimento di Scienze Chimiche and Padua Quantum Technologies Research Center, University of Padova, I-35131 Padova, Italy; Institute of Nanosciences, Consiglio Nazionale delle Ricerche CNR-Nano, I-41125 Modena, Italy; orcid.org/0000-0001-6707-108X;
Email: stefano.corni@unipd.it

Authors

Jacopo Fregoni — Dipartimento di Scienze Chimiche, University of Padova, I-35131 Padova, Italy; Institute of Nanosciences, Consiglio Nazionale delle Ricerche CNR-Nano, I-41125 Modena, Italy; orcid.org/0000-0002-4356-3060

Tor S. Haugland — Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway

Silvio Pipolo — UCCS Unité de Catalyse et Chimie du Solide, Université de Lille, F-59000 Lille, France

Tommaso Giovannini — Scuola Normale Superiore, I-56126 Pisa, Italy; orcid.org/0000-0002-5637-2853

Complete contact information is available at:
https://pubs.acs.org/10.1021/acs.nanolett.1c02162

Author Contributions

J.F. and T.S.H. contributed equally to the realization of the present work.

Notes

The authors declare no competing financial interest.

■ ACKNOWLEDGMENTS

We acknowledge Enrico Ronca's contributions to the charge displacement analysis. The authors also thank Antonio I. Fernandez Dominguez and Johannes Feist for the fruitful scientific discussion on the field quantization. This work has been funded by the European Research Council through Grants ERC-2015-CoG-681285 (J.F., PI S.C.), and the Research Council of Norway through FRINATEK projects 263110 and 275506 (T.S.H., PI H.K.). We acknowledge computing resources through UNINETT Sigma2 - the National Infrastructure for High Performance Computing and Data Storage in Norway, through project number NN2962k. We further acknowledge computer resources from the HPC Center at Scuola Normale Superiore di Pisa.

■ REFERENCES

(1) Tame, M.; Mcenery, K. R.; Ozdemir, S.; Lee, J.; Maier, S.; Kim, M. Quantum plasmonics. *Nat. Phys.* **2013**, *9*, 329–340.

(2) Hugall, J. T.; Singh, A.; van Hulst, N. F. Plasmonic Cavity Coupling. *ACS Photonics* **2018**, *5*, 43–53.

(3) Fofang, N. T.; Grady, N. K.; Fan, Z.; Govorov, A. O.; Halas, N. J. Plexciton Dynamics: Exciton Plasmon Coupling in a J Aggregate Au Nanoshell Complex Provides a Mechanism for Nonlinearity. *Nano Lett.* **2011**, *11*, 1556–1560.

(4) Li, J.-F.; Li, C.-Y.; Aroca, R. F. Plasmon-enhanced fluorescence spectroscopy. *Chem. Soc. Rev.* **2017**, *46*, 3962–3979.

(5) Polak, D.; et al. Manipulating molecules with strong coupling: harvesting triplet excitons in organic exciton microcavities. *Chem. Sci.* **2020**, *11*, 343–354.

(6) Eizner, E.; Martínez-Martínez, L. A.; Yuen-Zhou, J.; Kéna-Cohen, S. Inverting singlet and triplet excited states using strong light-matter coupling. *Sci. Adv.* **2019**, *5*, eaax4482.

(7) Climent, C.; Galego, J.; Garcia-Vidal, F. J.; Feist, J. Plasmonic Nanocavities Enable Self-Induced Electrostatic Catalysis. *Angew. Chem., Int. Ed.* **2019**, *58*, 8698–8702.

(8) Thomas, A.; George, J.; Shalabney, A.; Dryzhakov, M.; Varma, S. J.; Moran, J.; Chervy, T.; Zhong, X.; Devaux, E.; Genet, C.; Hutchison, J. A.; Ebbesen, T. W. Ground-State Chemical Reactivity under Vibrational Coupling to the Vacuum Electromagnetic Field. *Angew. Chem., Int. Ed.* **2016**, *55*, 11462–11466.

(9) Martínez-Martínez, L. A.; Du, M.; Ribeiro, R. F.; Kéna-Cohen, S.; Yuen-Zhou, J. Polariton-Assisted Singlet Fission in Acene Aggregates. *J. Phys. Chem. Lett.* **2018**, *9*, 1951–1957.

(10) Fregoni, J.; Granucci, G.; Persico, M.; Corni, S. Strong Coupling with Light Enhances the Photoisomerization Quantum Yield of Azobenzene. *Chem.* **2020**, *6*, 250–265.

(11) Gu, B.; Mukamel, S. Manipulating nonadiabatic conical intersection dynamics by optical cavities. *Chem. Sci.* **2020**, *11*, 1290–1298.

(12) Antoniou, P.; Suchanek, F.; Varner, J. F.; Foley, J. J. Role of Cavity Losses on Nonadiabatic Couplings and Dynamics in Polaritonic Chemistry. *J. Phys. Chem. Lett.* **2020**, *11*, 9063–9069.

(13) Mauro, L.; Caicedo, K.; Jonusauskas, G.; Avriller, R. Charge-transfer chemical reactions in nanofluidic Fabry-Pérot cavities. *Phys. Rev. B: Condens. Matter Mater. Phys.* **2021**, *103*, 165412.

(14) Coccia, E.; Corni, S. Role of coherence in the plasmonic control of molecular absorption. *J. Chem. Phys.* **2019**, *151*, 044703.

(15) Galego, J.; Garcia-Vidal, F. J.; Feist, J. Cavity-Induced Modifications of Molecular Structure in the Strong-Coupling Regime. *Phys. Rev. X* **2015**, *5*, 041022.

(16) Kowalewski, M.; Bennett, K.; Mukamel, S. Non-adiabatic dynamics of molecules in optical cavities. *J. Chem. Phys.* **2016**, *144*, 054309.
(17) Neuman, T.; Esteban, R.; Casanova, D.; García-Vidal, F. J.; Aizpurua, J. Coupling of Molecular Emitters and Plasmonic Cavities beyond the Point-Dipole Approximation. *Nano Lett.* **2018**, *18*, 2358–2364.

(18) Doppagne, B.; Neuman, T.; Soria-Martinez, R.; López, L. E. P.; Bulou, H.; Romeo, M.; Berciaud, S.; Scheurer, F.; Aizpurua, J.; Schull, G. Single-molecule tautomerization tracking through space- and time-resolved fluorescence spectroscopy. *Nat. Nanotechnol.* **2020**, *15*, 207–211.

(19) Kosik, M.; Burlayenko, O.; Rockstuhl, C.; Fernandez-Corbaton, I.; Slowik, K. Interaction of atomic systems with quantum vacuum beyond electric dipole approximation. *Sci. Rep.* **2020**, *10*, 5879.

(20) Ojambati, O. S.; Chikkaraddy, R.; Deacon, W. D.; Horton, M.; Kos, D.; Turek, V. A.; Keyser, U. F.; Baumberg, J. J. Quantum electrodynamics at room temperature coupling a single vibrating molecule with a plasmonic nanocavity. *Nat. Commun.* **2019**, *10*, 1049.

(21) Chen, X.; Liu, P.; Hu, Z.; Jensen, L. High-resolution tip-enhanced Raman scattering probes sub-molecular density changes. *Nat. Commun.* **2019**, *10*, 2567.

(22) Haugland, T. S.; Ronca, E.; Kjønstad, E. F.; Rubio, A.; Koch, H. Coupled Cluster Theory for Molecular Polaritons: Changing Ground and Excited States. *Phys. Rev. X* **2020**, *10*, 041043.

(23) Haugland, T. S.; Schäfer, C.; Ronca, E.; Rubio, A.; Koch, H. Intermolecular interactions in optical cavities: An ab initio QED study. *J. Chem. Phys.* **2021**, *154*, 094113.

(24) Buhmann, S. Y.; Butcher, D. T.; Scheel, S. Macroscopic quantum electrodynamics in nonlocal and nonreciprocal media. *New J. Phys.* **2012**, *14*, 083034.

(25) Feist, J.; Fernández-Domínguez, A. I.; García-Vidal, F. J. Macroscopic QED for quantum nanophotonics: emitter-centered modes as a minimal basis for multiemitter problems. *Nanophotonics* **2020**, *10*, 477–489.

(26) Helgaker, T.; Jørgensen, P.; Olsen, J. *Molecular Electronic-Structure Theory*; John Wiley & Sons, Ltd: Chichester, U.K., 2000.

(27) Jaynes, E. T.; Cummings, F. W. Comparison of quantum and semiclassical radiation theories with application to the beam maser. *Proc. IEEE* **1963**, *51*, 89.

(28) Sinha-Roy, R.; García-González, P.; Weissker, H.-C.; Rabilloud, F.; Fernández-Domínguez, A. I. Classical and ab Initio Plasmonics Meet at Sub-nanometric Noble Metal Rods. *ACS Photonics* **2017**, *4*, 1484–1493.

(29) Urbieta, M.; Barbry, M.; Zhang, Y.; Koval, P.; Sánchez-Portal, D.; Zabala, N.; Aizpurua, J. Atomic-Scale Lightning Rod Effect in Plasmonic Picocavities: A Classical View to a Quantum Effect. *ACS Nano* **2018**, *12*, 585–595.

(30) Yang, B.; Chen, G.; Ghafoor, A.; Zhang, Y.; Zhang, Y.; Zhang, Y.; Luo, Y.; Yang, J.; Sandoghdar, V.; Aizpurua, J.; Dong, Z.; Hou, J. G. Sub-nanometre resolution in single-molecule photoluminescence imaging. *Nat. Photonics* **2020**, *14*, 693–699.

(31) Judge, A. C.; Steel, M. J.; Sipe, J. E.; de Sterke, C. M. Canonical quantization of macroscopic electrodynamics in a linear, inhomogeneous magnetoelectric medium. *Phys. Rev. A: At., Mol., Opt. Phys.* **2013**, *87*, 033824.

(32) Mennucci, B.; Corni, S. Multiscale modelling of photoinduced processes in composite systems. *Nat. Rev. Chem.* **2019**, *3*, 315–330.

(33) Kelly, K. L.; Coronado, E.; Zhao, L. L.; Schatz, G. C. The Optical Properties of Metal Nanoparticles: The Influence of Size, Shape, and Dielectric Environment. *J. Phys. Chem. B* **2003**, *107*, 668–677.

(34) Cherqui, C.; Thakkar, N.; Li, G.; Camden, J. P.; Masiello, D. J. Characterizing Localized Surface Plasmons Using Electron Energy-Loss Spectroscopy. *Annu. Rev. Phys. Chem.* **2016**, *67*, 331–357.

(35) Trügler, A.; Hohenester, U. Strong coupling between a metallic nanoparticle and a single molecule. *Phys. Rev. B: Condens. Matter Mater. Phys.* **2008**, *77*, 115403.

(36) Cancès, E.; Mennucci, B.; Tomasi, J. A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics. *J. Chem. Phys.* **1997**, *107*, 3032–3041.

(37) Corni, S.; Pipolo, S.; Cammi, R. Equation of Motion for the Solvent Polarization Apparent Charges in the Polarizable Continuum Model: Application to Real-Time TDDFT. *J. Phys. Chem. A* **2015**, *119*, 5405–5416.

(38) Boyd, R. W. *Nonlinear Optics*, 3rd ed.; Academic Press, Inc., 2008; pp 135–171.

(39) Guido, C. A.; Rosa, M.; Cammi, R.; Corni, S. An open quantum system theory for polarizable continuum models. *J. Chem. Phys.* **2020**, *152*, 174114.

(40) Kristensen, K.; Kauczor, J.; Thorvaldsen, A. J.; Jorgensen, P.; Kjaergaard, T.; Rizzo, A. Damped response theory description of two-photon absorption. *J. Chem. Phys.* **2011**, *134*, 214104.

(41) Delga, A.; Feist, J.; Bravo-Abad, J.; Garcia-Vidal, F. J. Quantum Emitters Near a Metal Nanoparticle: Strong Coupling and Quenching. *Phys. Rev. Lett.* **2014**, *112*, 253601.

(42) Flick, J.; Ruggenthaler, M.; Appel, H.; Rubio, A. Atoms and molecules in cavities, from weak to strong coupling in quantum-electrodynamics (QED) chemistry. *Proc. Natl. Acad. Sci. U. S. A.* **2017**, *114*, 3026–3034.

(43) Schäfer, C.; Ruggenthaler, M.; Rubio, A. Ab initio non-relativistic quantum electrodynamics: Bridging quantum chemistry and quantum optics from weak to strong coupling. *Phys. Rev. A: At., Mol., Opt. Phys.* **2018**, *98*, 043801.

(44) Novotny, L.; Hecht, B. *Principles of Nano-Optics*; Cambridge University Press, 2006.

(45) Corni, S.; Tomasi, J. Enhanced response properties of a chromophore physisorbed on a metal particle. *J. Chem. Phys.* **2001**, *114*, 3739–3751.

(46) Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. *Chem. Rev.* **2005**, *105*, 2999–3094.

(47) Jensen, L. L.; Jensen, L. Atomistic Electrodynamics Model for Optical Properties of Silver Nanoclusters. *J. Phys. Chem. C* **2009**, *113*, 15182–15190.

(48) Chen, X.; Liu, P.; Jensen, L. Atomistic electrodynamics simulations of plasmonic nanoparticles. *J. Phys. D: Appl. Phys.* **2019**, *52*, 363002.

(49) Giovannini, T.; Rosa, M.; Corni, S.; Cappelli, C. A classical picture of subnanometer junctions: an atomistic Drude approach to nanoplasmonics. *Nanoscale* **2019**, *11*, 6004–6015.

(50) Baranov, D. G.; Munkhbat, B.; Zhukova, E.; Bisht, A.; Canales, A.; Rousseaux, B.; Johansson, G.; Antosiewicz, T. J.; Shegai, T. Ultrastrong coupling between nanoparticle plasmons and cavity photons at ambient conditions. *Nat. Commun.* **2020**, *11*, 2715.