Relevance of the Quadratic Diamagnetic and Self-Polarization Terms in Cavity Quantum Electrodynamics

Molecular Polarizability under Vibrational Strong Coupling

Polaritonic chemistry offers the possibility of modifying molecular properties and even influencing chemical reactivity through strong coupling between vibrational transitions and confined light modes in optical cavities. Despite considerable theoretical progress, and due to the complexity of the coupled light-matter system, the fundamental mechanism of how and if collective strong coupling can induce local changes in individual molecules is still unclear. We derive an analytical formulation of static polarizabilities within linear-response theory for molecules under strong coupling using the cavity Born-Oppenheimer Hartree-Fock ansatz. This ab-initio method consistently describes vibrational strong coupling and electron-photon interactions even for ensembles of molecules. For different types of molecular ensembles, we observed local changes in the polarizabilities and dipole moments that are induced by collective strong coupling. Furthermore, we used the polarizabilities to calculate vibro-polaritonic Raman spectra in the harmonic approximation. This allows us to comprehensively compare the effect of vibrational strong coupling on IR and Raman spectra on an equal footing.

Strong Coupling Quantum Electrodynamics Hartree-Fock Response Theory

The development of reliable ab initio methods for light-matter strong coupling is necessary for a deeper understanding of molecular polaritons. The recently developed strong coupling quantum electrodynamics Hartree-Fock model (SC-QED-HF) provides cavity-consistent molecular orbitals, overcoming several difficulties related to the simpler QED-HF wave function. In this paper, we further develop this method by implementing the response theory for SC-QED-HF. We compare the derived linear response equations with the time-dependent QED-HF theory and discuss the validity of equivalence relations connecting matter and electromagnetic observables. Our results show that electron-photon correlation induces an excitation redshift compared to the time-dependent QED-HF energies, and we discuss the effect of the dipole self-energy on the ground and excited state properties with different basis sets.

Controlling Molecular Dynamics by Exciting Atoms in a Cavity

Placing an atom and a molecule in a cavity opens the door to initialize molecular dynamics by exciting a level of the atom. This approach enlarges the range of choosing the light source to trigger molecular dynamics substantially. The interplay of the atomic, molecular, and photonic populations gives rise to rich dynamics. The cavity photon plays the role of a mediator between the atom and the molecule and it is found that the photonic population is rather low throughout and its evolution follows that of the molecule. Cavities are known to be subject to losses. In spite of the losses it is demonstrated that the presence of the atom gives rise to a long-lived dynamics that should be of relevance for experimental investigations. The presence of more atoms and molecules is expected to further enrich the dynamics.

Density-Functional Theory for the Dicke Hamiltonian

A detailed analysis of density-functional theory for quantum-electrodynamical model systems is provided. In particular, the quantum Rabi model, the Dicke model, and a generalization of the latter to multiple modes are considered. We prove a Hohenberg-Kohn theorem that manifests the magnetization and displacement as internal variables, along with several representability results. The constrained-search functionals for pure states and ensembles are introduced and analyzed. We find the optimizers for the pure-state constrained-search functional to be low-lying eigenstates of the Hamiltonian and, based on the properties of the optimizers, we formulate an adiabatic-connection formula. In the reduced case of the Rabi model we can even show differentiability of the universal density functional, which amounts to unique pure-state v-representability.

Spontaneous single-molecule dissociation in infrared nanocavities

Ultrastrong light-matter interaction with molecular vibrations in infrared cavities has emerged as a tool for manipulating and controlling chemical reactivity. By studying the wavepacket dynamics of an individual polar diatomic molecule in a quantized infrared electromagnetic environment, we show that chemical bonds can efficiently dissociate in the absence of additional thermal or coherent energy sources, provided that the coupled system is prepared in a suitable diabatic state. Using hydrogen fluoride as a case study, we predict dissociation probabilities of up to 35% in less than 200 fs for a vibration-cavity system that is rapidly initialized with a low number of bare vibrational and cavity excitations. We develop a simple and general analytical model based on the multipolar formulation of quantum electrodynamics to show that the Bloch-Seigert shift of the bare vibrational ground state is a predictor of a threshold coupling strength below which no spontaneous dissociation is expected. The role of state-dependent permanent dipole moments in the light-matter interaction process is clarified. Our work paves the way toward the development of vacuum-assisted chemical reactors powered by ultrastrong light-matter interaction at the single-molecule level.

On-demand heralded MIR single-photon source using a cascaded quantum system

We propose a mechanism for generating single photons in the mid-infrared (MIR) using a solid-state or molecular quantum emitter. The scheme uses cavity quantum electrodynamics (QED) effects to selectively enhance a Frank-Condon transition, deterministically preparing a single Fock state of a polar phonon mode. By coupling the phonon mode to an antenna, the resulting excitation is then radiated to the far field as a single photon with a frequency matching the phonon mode. By combining macroscopic QED calculations with methods from open quantum system theory, we show that optimal parameters to generate these MIR photons occur for modest light-matter coupling strengths, which are achievable with state-of-the-art technologies. Combined, the cascaded system we propose provides a quasi-deterministic source of heralded single photons in a regime of the electromagnetic spectrum where this previously was not possible.

Quantum-Electrodynamical Density-Functional Theory Exemplified by the Quantum Rabi Model

The key features of density-functional theory (DFT) within a minimalistic implementation of quantum electrodynamics are demonstrated, thus allowing to study elementary properties of quantum-electrodynamical density-functional theory (QEDFT). We primarily employ the quantum Rabi model that describes a two-level system coupled to a single photon mode and also discuss the Dicke model, where multiple two-level systems couple to the same photon mode. In these settings, the density variables of the system are the polarization and the displacement of the photon field. We give analytical expressions for the constrained-search functional and the exchange-correlation potential and compare them to established results from QEDFT. We further derive a form for the adiabatic connection that is almost explicit in the density variables, up to only a nonexplicit correlation term that gets bounded both analytically and numerically. This allows several key features of DFT to be studied without approximations.

Impact of Dipole Self-Energy on Cavity-Induced Nonadiabatic Dynamics

The coupling of matter to the quantized electromagnetic field of a plasmonic or optical cavity can be harnessed to modify and control chemical and physical properties of molecules. In optical cavities, a term known as the dipole self-energy (DSE) appears in the Hamiltonian to ensure gauge invariance. The aim of this work is twofold. First, we introduce a method, which has its own merits and complements existing methods, to compute the DSE. Second, we study the impact of the DSE on cavity-induced nonadiabatic dynamics in a realistic system. For that purpose, various matrix elements of the DSE are computed as functions of the nuclear coordinates and the dynamics of the system after laser excitation is investigated. The cavity is known to induce conical intersections between polaritons, which gives rise to substantial nonadiabatic effects. The DSE is shown to slightly affect these light-induced conical intersections and, in particular, break their symmetry.

Disentangling collective coupling in vibrational polaritons with double quantum coherence spectroscopy

Vibrational polaritons are formed by strong coupling of molecular vibrations and photon modes in an optical cavity. Experiments have demonstrated that vibrational strong coupling can change molecular properties and even affect chemical reactivity. However, the interactions in a molecular ensemble are complex, and the exact mechanisms that lead to modifications are not fully understood yet. We simulate two-dimensional infrared spectra of molecular vibrational polaritons based on the double quantum coherence technique to gain further insight into the complex many-body structure of these hybrid light-matter states. Double quantum coherence uniquely resolves the excitation of hybrid light-matter polaritons and allows one to directly probe the anharmonicities of the resulting states. By combining the cavity Born-Oppenheimer Hartree-Fock ansatz with a full quantum dynamics simulation of the corresponding eigenstates, we go beyond simplified model systems. This allows us to study the influence of self-polarization and the response of the electronic structure to the cavity interaction on the spectral features even beyond the single-molecule case.

Cavity-modified local and non-local electronic interactions in molecular ensembles under vibrational strong coupling

Resonant vibrational strong coupling (VSC) between molecular vibrations and quantized field modes of low-frequency optical cavities constitutes the conceptual cornerstone of vibro-polaritonic chemistry. In this work, we theoretically investigate the role of complementary nonresonant electron-photon interactions in the cavity Born-Oppenheimer (CBO) approximation. In particular, we study cavity-induced modifications of local and non-local electronic interactions in dipole-coupled molecular ensembles under VSC. Methodologically, we combine CBO perturbation theory (CBO-PT) [E. W. Fischer and P. Saalfrank, J. Chem. Theory Comput. 19, 7215 (2023)] with non-perturbative CBO Hartree-Fock (HF) and coupled cluster (CC) theories. In a first step, we derive up to second-order CBO-PT cavity potential energy surfaces, which reveal non-trivial intra- and inter-molecular corrections induced by the cavity. We then introduce the concept of a cavity reaction potential (CRP), minimizing the electronic energy in the cavity subspace to discuss vibro-polaritonic reaction mechanisms. We present reformulations of CBO-HF and CBO-CC approaches for CRPs and derive second-order approximate CRPs from CBO-PT for unimolecular and bimolecular scenarios. In the unimolecular case, we find small local modifications of molecular potential energy surfaces for selected isomerization reactions dominantly captured by the first-order dipole fluctuation correction. Excellent agreement between CBO-PT and non-perturbative wave function results indicates minor VSC-induced state relaxation effects in the single-molecule limit. In the bimolecular scenario, CBO-PT reveals an explicit coupling of interacting dimers to cavity modes besides cavity-polarization dependent dipole-induced dipole and van der Waals interactions with enhanced long-range character. An illustrative CBO-coupled cluster theory with singles and doubles-based numerical analysis of selected molecular dimer models provides a complementary non-perturbative perspective on cavity-modified intermolecular interactions under VSC.

Analytical Evaluation of Ground State Gradients in Quantum Electrodynamics Coupled Cluster Theory

Analytical gradients of potential energy surfaces play a central role in quantum chemistry, allowing for molecular geometry optimizations and molecular dynamics simulations. In strong coupling conditions, potential energy surfaces can account for strong interactions between matter and the quantized electromagnetic field. In this paper, we derive expressions for the ground state analytical gradients in quantum electrodynamics coupled cluster theory. We also present a Cholesky-based implementation for the coupled cluster singles and doubles model. We report timings to show the performance of the implementation and present optimized geometries to highlight cavity-induced molecular orientation effects in strong coupling conditions.

Quantized embedding approaches for collective strong coupling-Connecting ab initio and macroscopic QED to simple models in polaritonics

Collective light-matter interactions have been used to control chemistry and energy transfer, yet accessible approaches that combine ab initio methodology with large many-body quantum optical systems are missing due to the fast increase in computational cost for explicit simulations. We introduce an accessible ab initio quantum embedding concept for many-body quantum optical systems that allows us to treat the collective coupling of molecular many-body systems effectively in the spirit of macroscopic quantum electrodynamics while keeping the rigor of ab initio quantum chemistry for the molecular structure. Our approach fully includes the quantum fluctuations of the polaritonic field and yet remains much simpler and more intuitive than complex embedding approaches such as dynamical mean-field theory. We illustrate the underlying assumptions by comparison to the Tavis-Cummings model. The intuitive application of the quantized embedding approach and its transparent limitations offer a practical framework for the field of ab initio polaritonic chemistry to describe collective effects in realistic molecular ensembles.

Controlling Plasmonic Catalysis via Strong Coupling with Electromagnetic Resonators

Plasmonic excitations decay within femtoseconds, leaving nonthermal (often referred to as "hot") charge carriers behind that can be injected into molecular structures to trigger chemical reactions that are otherwise out of reach─a process known as plasmonic catalysis. In this Letter, we demonstrate that strong coupling between resonator structures and plasmonic nanoparticles can be used to control the spectral overlap between the plasmonic excitation energy and the charge injection energy into nearby molecules. Our atomistic description couples real-time density-functional theory self-consistently to an electromagnetic resonator structure via the radiation-reaction potential. Control over the resonator provides then an additional knob for nonintrusively enhancing plasmonic catalysis, here more than 6-fold, and dynamically reacting to deterioration of the catalyst─a new facet of modern catalysis.

Unraveling abnormal collective effects via the non-monotonic number dependence of electron transfer in confined electromagnetic fields

Strong light-matter coupling within an optical cavity leverages the collective interactions of molecules and confined electromagnetic fields, giving rise to the possibilities of modifying chemical reactivity and molecular properties. While collective optical responses, such as enhanced Rabi splitting, are often observed, the overall effect of the cavity on molecular systems remains ambiguous for a large number of molecules. In this paper, we investigate the non-adiabatic electron transfer process in electron donor-acceptor pairs influenced by collective excitation and local molecular dynamics. Using the timescale difference between reorganization and thermal fluctuations, we derive analytical formulas for the electron transfer rate constant and the polariton relaxation rate. These formulas apply to any number of molecules (N) and account for the collective effect as induced by cavity photon coupling. Our findings reveal a non-monotonic dependence of the rate constant on N, which can be understood by the interplay between electron transfer and polariton relaxation. As a result, the cavity-induced quantum yield increases linearly with N for small N (as predicted by a simple Dicke model) but shows a turnover and suppression for large N. We also interrelate the thermal bath frequency and the number of molecules, suggesting the optimal number for maximizing enhancement. The analysis provides an analytical insight for understanding the collective excitation of light and electron transfer, helping to predict the optimal condition for effective cavity-controlled chemical reactivity.

Exact factorization of the photon-electron-nuclear wavefunction: Formulation and coupled-trajectory dynamics

We employ the exact-factorization formalism to study the coupled dynamics of photons, electrons, and nuclei at the quantum mechanical level, proposing illustrative examples of model situations of nonadiabatic dynamics and spontaneous emission of electron-nuclear systems in the regime of strong light-matter coupling. We make a particular choice of factorization for such a multi-component system, where the full wavefunction is factored as a conditional electronic amplitude and a marginal photon-nuclear amplitude. Then, we apply the coupled-trajectory mixed quantum-classical (CTMQC) algorithm to perform trajectory-based simulations, by treating photonic and nuclear degrees of freedom on equal footing in terms of classical-like trajectories. The analysis of the time-dependent potentials of the theory along with the assessment of the performance of CTMQC allows us to point out some limitations of the current approximations used in CTMQC. Meanwhile, comparing CTMQC with other trajectory-based algorithms, namely multi-trajectory Ehrenfest and Tully surface hopping, demonstrates the better quality of CTMQC predictions.

Quantum nature of reactivity modification in vibrational polariton chemistry

In this work, we present a mixed quantum-classical open quantum system dynamics method for studying rate modifications of ground-state chemical reactions in an optical cavity under vibrational strong-coupling conditions. In this approach, the cavity radiation mode is treated classically with a mean-field nuclear force averaging over the remaining degrees of freedom, both within the system and the environment, which are handled quantum mechanically within the hierarchical equations of motion framework. Using this approach, we conduct a comparative analysis by juxtaposing the mixed quantum-classical results with fully quantum-mechanical simulations. After eliminating spurious peaks that can occur when not using the rigorous definition of the rate constant, we confirm the crucial role of the quantum nature of the cavity radiation mode in reproducing the resonant peak observed in the cavity frequency-dependent rate profile. In other words, it appears necessary to explicitly consider the quantized photonic states in studying reactivity modification in vibrational polariton chemistry (at least for the model systems studied in this work), as these phenomena stem from cavity-induced reaction pathways involving resonant energy exchanges between photons and molecular vibrational transitions.

Exciton energy transfer inside cavity-A benchmark study of polaritonic dynamics using the surface hopping method

Strong coupling between the molecular system and photon inside the cavity generates polaritons, which can alter reaction rates by orders of magnitude. In this work, we benchmark the surface hopping method to simulate non-adiabatic dynamics in a cavity. The comparison is made against a numerically exact method (the hierarchical equations of motion) for a model system investigating excitonic energy transfer for a broad range of parameters. Surface hopping captures the effects of the radiation mode well, both at resonance and off-resonance. We have further investigated parameters that can increase or decrease the rate of population transfer, and we find that surface hopping in general can capture both effects well. Finally, we show that the dipole self-energy term within our parameter regime does not significantly affect the system's dynamics.

Do Molecular Geometries Change Under Vibrational Strong Coupling?

As pioneering experiments have shown, strong coupling between molecular vibrations and light modes in an optical cavity can significantly alter molecular properties and even affect chemical reactivity. However, the current theoretical description is limited and far from complete. To explore the origin of this exciting observation, we investigate how the molecular structure changes under strong light-matter coupling using an ab initio method based on the cavity Born-Oppenheimer Hartree-Fock ansatz. By optimizing H2O and H2O2 resonantly coupled to cavity modes, we study the importance of reorientation and geometric relaxation. In addition, we show that the inclusion of one or two cavity modes can change the observed results. On the basis of our findings, we derive a simple concept to estimate the effect of the cavity interaction on the molecular geometry using the molecular polarizability and the dipole moments.

Tailoring of the polarization-resolved second harmonic generation in two-dimensional semiconductors

Second harmonic generation is a non-linear optical phenomenon in which coherent radiation with frequency ω interacts with a non-centrosymmetric material and produces coherent radiation at frequency 2ω. Owing to the exciting physical phenomena that take place during the non-linear optical excitation at the nanoscale, there is currently extensive research in the non-linear optical responses of nanomaterials, particularly in low-dimensional materials. Here, we review recent advancements in the polarization-resolved second harmonic generation propertied from atomically thin two-dimensional (2D) crystals and present a unified theoretical framework to account for their nonlinear optical response. Two major classes of 2D materials are particularly investigated, namely metal chalcogenides and perovskites. The first attempts to tune and control the second harmonic generation properties of such materials via the application of specific nanophotonic schemes are additionally demonstrated and discussed. Besides presenting recent advances in the field, this work also delineates existing limitations and highlights emerging possibilities and future prospects in this field.

Extending the Tavis-Cummings model for molecular ensembles-Exploring the effects of dipole self-energies and static dipole moments

Strong coupling of organic molecules to the vacuum field of a nanoscale cavity can be used to modify their chemical and physical properties. We extend the Tavis-Cummings model for molecular ensembles and show that the often neglected interaction terms arising from the static dipole moment and the dipole self-energy are essential for a correct description of the light-matter interaction in polaritonic chemistry. On the basis of a full quantum description, we simulate the excited-state dynamics and spectroscopy of MgH+ molecules resonantly coupled to an optical cavity. We show that the inclusion of static dipole moments and the dipole self-energy is necessary to obtain a consistent model. We construct an efficient two-level system approach that reproduces the main features of the real molecular system and may be used to simulate larger molecular ensembles.

Insights into the mechanisms of optical cavity-modified ground-state chemical reactions

In this work, we systematically investigate the mechanisms underlying the rate modification of ground-state chemical reactions in an optical cavity under vibrational strong-coupling conditions. We employ a symmetric double-well description of the molecular potential energy surface and a numerically exact open quantum system approach-the hierarchical equations of motion in twin space with a matrix product state solver. Our results predict the existence of multiple peaks in the photon frequency-dependent rate profile for a strongly anharmonic molecular system with multiple vibrational transition energies. The emergence of a new peak in the rate profile is attributed to the opening of an intramolecular reaction pathway, energetically fueled by the cavity photon bath through a resonant cavity mode. The peak intensity is determined jointly by kinetic factors. Going beyond the single-molecule limit, we examine the effects of the collective coupling of two molecules to the cavity. We find that when two identical molecules are simultaneously coupled to the same resonant cavity mode, the reaction rate is further increased. This additional increase is associated with the activation of a cavity-induced intermolecular reaction channel. Furthermore, the rate modification due to these cavity-promoted reaction pathways remains unaffected, regardless of whether the molecular dipole moments are aligned in the same or opposite direction as the light polarization.

Cavity Quantum Electrodynamics Enables para- and ortho-Selective Electrophilic Bromination of Nitrobenzene

Coupling molecules to a quantized radiation field inside an optical cavity has shown great promise to modify chemical reactivity. In this work, we show that the ground-state selectivity of the electrophilic bromination of nitrobenzene can be fundamentally changed by strongly coupling the reaction to the cavity, generating ortho- or para-substituted products instead of the meta product. Importantly, these are products that are not obtained from the same reaction outside the cavity. A recently developed ab initio approach was used to theoretically compute the relative energies of the cationic Wheland intermediates, which indicate the kinetically preferred bromination site for all products. Performing an analysis of the ground-state electron density for the Wheland intermediates inside and outside the cavity, we demonstrate how strong coupling induces reorganization of the molecular charge distribution, which in turn leads to different bromination sites directly dependent on the cavity conditions. Overall, the results presented here can be used to understand cavity induced changes to ground-state chemical reactivity from a mechanistic perspective as well as to directly connect frontier theoretical simulations to state-of-the-art, but realistic, experimental cavity conditions.

Coupling polyatomic molecules to lossy nanocavities: Lindblad vs Schrödinger description

The use of cavities to impact molecular structure and dynamics has become popular. As cavities, in particular plasmonic nanocavities, are lossy and the lifetime of their modes can be very short, their lossy nature must be incorporated into the calculations. The Lindblad master equation is commonly considered an appropriate tool to describe this lossy nature. This approach requires the dynamics of the density operator and is thus substantially more costly than approaches employing the Schrödinger equation for the quantum wave function when several or many nuclear degrees of freedom are involved. In this work, we compare numerically the Lindblad and Schrödinger descriptions discussed in the literature for a molecular example where the cavity is pumped by a laser. The laser and cavity properties are varied over a range of parameters. It is found that the Schrödinger description adequately describes the dynamics of the polaritons and emission signal as long as the laser intensity is moderate and the pump time is not much longer than the lifetime of the cavity mode. Otherwise, it is demonstrated that the Schrödinger description gradually fails. We also show that the failure of the Schrödinger description can often be remedied by renormalizing the wave function at every step of time propagation. The results are discussed and analyzed.

Extracting kinetic information from short-time trajectories: relaxation and disorder of lossy cavity polaritons

The emerging field of molecular cavity polaritons has stimulated a surge of experimental and theoretical activities and presents a unique opportunity to develop the many-body simulation methodology. This paper presents a numerical scheme for the extraction of key kinetic information of lossy cavity polaritons based on the transfer tensor method (TTM). Steady state, relaxation timescales, and oscillatory phenomena can all be deduced directly from a set of transfer tensors without the need for long-time simulation. Moreover, we generalize TTM to disordered systems by sampling dynamical maps and achieve fast convergence to disordered-averaged dynamics using a small set of realizations. Together, these techniques provide a toolbox for characterizing the interplay of cavity loss, disorder, and cooperativity in polariton relaxation and allow us to predict unusual dependences on the initial excitation state, photon decay rate, strength of disorder, and the type of cavity models. Thus, using the example of cavity polaritons, we have demonstrated significant potential in the use of the TTM toward both the efficient computation of long-time polariton dynamics and the extraction of crucial kinetic information about polariton relaxation from a small set of short-time trajectories.

Investigating the collective nature of cavity-modified chemical kinetics under vibrational strong coupling

In this paper, we develop quantum dynamical methods capable of treating the dynamics of chemically reacting systems in an optical cavity in the vibrationally strong-coupling (VSC) limit at finite temperatures and in the presence of a dissipative solvent in both the few and many molecule limits. In the context of two simple models, we demonstrate how reactivity in the collective VSC regime does not exhibit altered rate behavior in equilibrium but may exhibit resonant cavity modification of reactivity when the system is explicitly out of equilibrium. Our results suggest experimental protocols that may be used to modify reactivity in the collective regime and point to features not included in the models studied, which demand further scrutiny.

Unraveling a Cavity-Induced Molecular Polarization Mechanism from Collective Vibrational Strong Coupling

We demonstrate that collective vibrational strong coupling of molecules in thermal equilibrium can give rise to significant local electronic polarizations in the thermodynamic limit. We do so by first showing that the full nonrelativistic Pauli-Fierz problem of an ensemble of strongly coupled molecules in the dilute-gas limit reduces in the cavity Born-Oppenheimer approximation to a cavity-Hartree equation for the electronic structure. Consequently, each individual molecule experiences a self-consistent coupling to the dipoles of all other molecules, which amount to non-negligible values in the thermodynamic limit (large ensembles). Thus, collective vibrational strong coupling can alter individual molecules strongly for localized "hotspots" within the ensemble. Moreover, the discovered cavity-induced polarization pattern possesses a zero net polarization, which resembles a continuous form of a spin glass (or better polarization glass). Our findings suggest that the thorough understanding of polaritonic chemistry, requires a self-consistent treatment of dressed electronic structure, which can give rise to numerous, so far overlooked, physical mechanisms.

Impact of Cavity on Molecular Ionization Spectra

Ionization phenomena have been widely studied for decades. With the advent of cavity technology, the question arises how quantum light affects molecular ionization. As the ionization spectrum is recorded from the neutral ground state, it is usually possible to choose cavities which exert negligible effect on the neutral ground state, but have significant impact on the ion and the ionization spectrum. Particularly interesting are cases where the ion exhibits conical intersections between close-lying electronic states, which gives rise to substantial nonadiabatic effects. Assuming single-molecule strong coupling, we demonstrate that vibrational modes irrelevant in the absence of a cavity play a decisive role when the molecule is in the cavity. Here, dynamical symmetry breaking is responsible for the ion-cavity coupling and high symmetry enables control of the coupling via molecular orientation relative to the cavity field polarization. Significant impact on the spectrum by the cavity is found and shown to even substantially increase for less symmetric molecules.

Generalized Born-Huang expansion under macroscopic quantum electrodynamics framework

Born-Huang expansion is the cornerstone for studying potential energy surfaces and non-adiabatic couplings (NACs) in molecular systems. However, the traditional approach is insufficient to describe the molecular system, which strongly interacts with quantum light. Inspired by the work by Schäfer et al., we develop the generalized Born-Huang expansion theory within a macroscopic quantum electrodynamics (QED) framework. The theory we present allows us to describe electromagnetic vacuum fluctuations in dielectric media and incorporate the effects of dressed photons (or polaritons) into NACs. With the help of the generalized Born-Huang expansion, we clearly classify electronic nuclear NACs, polaritonic nuclear NACs, and polaritonic electronic NACs. Furthermore, to demonstrate the advantage of the macroscopic QED framework, we estimate polaritonic electronic NACs without any free parameter, such as the effective mode volume, and demonstrate the distance dependence of the polaritonic electronic NACs in a silver planar system. In addition, we take a hydrogen atom in free space as an example and derive spontaneous emission rates from photonic electronic NACs (polaritonic electronic NACs are reduced to photonic electronic NACs). We believe that this work not only provides an avenue for the theoretical exploration of NACs in a nucleus-electron-polariton coupled system but also offers a more comprehensive understanding for molecules coupled with quantum light.

Machine Learning for Polaritonic Chemistry: Accessing Chemical Kinetics

Altering chemical reactivity and material structure in confined optical environments is on the rise, and yet, a conclusive understanding of the microscopic mechanisms remains elusive. This originates mostly from the fact that accurately predicting vibrational and reactive dynamics for soluted ensembles of realistic molecules is no small endeavor, and adding (collective) strong light-matter interaction does not simplify matters. Here, we establish a framework based on a combination of machine learning (ML) models, trained using density-functional theory calculations and molecular dynamics to accelerate such simulations. We then apply this approach to evaluate strong coupling, changes in reaction rate constant, and their influence on enthalpy and entropy for the deprotection reaction of 1-phenyl-2-trimethylsilylacetylene, which has been studied previously both experimentally and using ab initio simulations. While we find qualitative agreement with critical experimental observations, especially with regard to the changes in kinetics, we also find differences in comparison with previous theoretical predictions. The features for which the ML-accelerated and ab initio simulations agree show the experimentally estimated kinetic behavior. Conflicting features indicate that a contribution of dynamic electronic polarization to the reaction process is more relevant than currently believed. Our work demonstrates the practical use of ML for polaritonic chemistry, discusses limitations of common approximations, and paves the way for a more holistic description of polaritonic chemistry.

Variational Lang-Firsov Approach Plus Møller-Plesset Perturbation Theory with Applications to Ab Initio Polariton Chemistry

We apply the Lang-Firsov (LF) transformation to electron-boson coupled Hamiltonians and variationally optimize the transformation parameters and molecular orbital coefficients to determine the ground state. Møller-Plesset (MP-n, with n = 2 and 4) perturbation theory is then applied on top of the optimized LF mean-field state to improve the description of electron-electron and electron-boson correlations. The method (LF-MP) is applied to several electron-boson coupled systems, including the Hubbard-Holstein model, diatomic molecule dissociation (H2, HF), and the modification of proton transfer reactions (malonaldehyde and aminopropenal) via the formation of polaritons in an optical cavity. We show that with a correction for the electron-electron correlation, the method gives quantitatively accurate energies comparable to that by exact diagonalization or coupled-cluster theory. The effects of multiple photon modes, spin polarization, and the comparison to the coherent state MP theory are also discussed.

Cavity-Mediated Molecular Entanglement and Generation of Non-classical States of Light

The generation and control of entanglement in a quantum mechanical system are critical elements of nearly all quantum applications. Molecular systems are promising candidates, with numerous degrees of freedom able to be targeted. However, knowledge of intersystem entanglement mechanisms in such systems is limited. In this work, we demonstrate the generation of entanglement between vibrational degrees of freedom in molecules via strong coupling to a cavity mode driven by a weak coherent field. In a bimolecular system, we show that entanglement can be generated not only between the cavity and molecular system but also between molecules. This process also results in the generation of nonclassical states of light, providing potential pathways for harnessing entanglement in molecular systems.

Ab Initio Calculations of Quantum Light-Matter Interactions in General Electromagnetic Environments

The emerging field of strongly coupled light-matter systems has drawn significant attention in recent years because of the prospect of altering both the physical and chemical properties of molecules and materials. Because this emerging field draws on ideas from both condensed-matter physics and quantum optics, it has attracted the attention of theoreticians from both fields. While the former often employ accurate descriptions of the electronic structure of the matter, the description of the electromagnetic environment is often oversimplified. In contrast, the latter often employs sophisticated descriptions of the electromagnetic environment while using oversimplified few-level approximations of the electronic structure. Both approaches are problematic because the oversimplified descriptions of the electronic system are incapable of describing effects such as light-induced structural changes in the electronic system, while the oversimplified descriptions of the electromagnetic environments can lead to unphysical predictions because the light-matter interactions strengths are misrepresented. In this work, we overcome these shortcomings and present the first method which can quantitatively describe both the electronic system and general electromagnetic environments from first principles. We realize this by combining macroscopic QED (MQED) with Quantum Electrodynamical Density-Functional Theory. To exemplify this approach, we consider the example of an absorbing spherical cavity and study the impact of different parameters of both the environment and the electronic system on the transition from weak-to-strong coupling for different aromatic molecules. As part of this work, we also provide an easy-to-use tool to calculate the cavity coupling strengths for simple cavity setups. Our work is a significant step toward parameter-free ab initio calculations for strongly coupled quantum light-matter systems and will help bridge the gap between theoretical methods and experiments in the field.

Quantum to Classical Cavity Chemistry Electrodynamics

Polaritonic chemistry has ushered in new avenues for controlling molecular dynamics. However, two key questions remain: (i) Can classical light sources elicit the same effects as certain quantum light sources on molecular systems? (ii) Can semiclassical treatments of light-matter interactions capture nontrivial quantum effects observed in molecular dynamics? This work presents a quantum-classical approach addressing issues of realizing cavity chemistry effects without actual cavities. It also highlights the limitations of the standard semiclassical light-matter interaction. It is demonstrated that classical light sources can mimic quantum effects up to the second order of light-matter interaction provided that the mean-field contribution, the symmetrized two-time correlation function, and the linear response function are the same in both situations. Numerical simulations show that the quantum-classical method aligns more closely with exact quantum molecular-only dynamics for quantum light states such as Fock states, superpositions of Fock states, and vacuum squeezed states than does the conventional semiclassical approach.

Ab Initio Vibro-Polaritonic Spectra in Strongly Coupled Cavity-Molecule Systems

Recent experiments have revealed the profound effect of strong light-matter interactions in optical cavities on the electronic ground state of molecular systems. This phenomenon, known as vibrational strong coupling, can modify reaction rates and induce the formation of molecular vibrational polaritons, hybrid states involving both photon modes, and vibrational modes of molecules. We present an ab initio methodology based on the cavity Born-Oppenheimer Hartree-Fock ansatz, which is specifically powerful for ensembles of molecules, to calculate vibro-polaritonic IR spectra. This method allows for a comprehensive analysis of these hybrid states. Our semiclassical approach, validated against full quantum simulations, reproduces key features of the vibro-polaritonic spectra. The underlying analytic gradients also allow for optimization of cavity-coupled molecular systems and performing semiclassical dynamics simulations.

Numerically Exact Solution for a Real Polaritonic System under Vibrational Strong Coupling in Thermodynamic Equilibrium: Loss of Light-Matter Entanglement and Enhanced Fluctuations

The first numerically exact simulation of a full ab initio molecular quantum system (HD+) under strong ro-vibrational coupling to a quantized optical cavity mode in thermal equilibrium is presented. Theoretical challenges in describing strongly coupled systems of mixed quantum statistics (bosons and Fermions) are discussed and circumvented by the specific choice of our molecular system. Our numerically exact simulations highlight the absence of zero temperature for the strongly coupled matter and light subsystems, due to cavity-induced noncanonical conditions. Furthermore, we explore the temperature dependency of light-matter quantum entanglement, which emerges for the ground state but is quickly lost already in the deep cryogenic regime. This is in contrast to predictions from the Jaynes-Cummings model, which is the standard starting point to model collective strong-coupling chemistry phenomenologically. Moreover, we find that the fluctuations of matter remain modified by the quantum nature of the thermal and vacuum-field fluctuations for significant temperatures, e.g., at ambient conditions. These observations (loss of entanglement and coupling to quantum fluctuations) have implications for the understanding and control of polaritonic chemistry and materials science, since a semiclassical theoretical description of light-matter interaction becomes reasonable, but the typical (classical) canonical equilibrium assumption for the nuclear subsystem remains violated. This opens the door for quantum fluctuation-induced stochastic resonance phenomena under vibrational strong coupling, which have been suggested as a plausible theoretical mechanism to explain the experimentally observed resonance phenomena in the absence of periodic driving that has not yet been fully understood.

Cavity-Modified Chemiluminescent Reaction of Dioxetane

Chemiluminescence is a thermally activated chemical process that emits a photon of light by forming a fraction of products in the electronic excited state. A well-known example of this spectacular phenomenon is the emission of light in the firefly beetle, where the formation of a four-membered cyclic peroxide compound and subsequent dissociation produce a light-emitting product. The smallest cyclic peroxide, dioxetane, also exhibits chemiluminescence but with a low quantum yield as compared to that of firefly dioxetane. Employing the strong light-matter coupling has recently been found to be an alternative strategy to modify the chemical reactivity. In the presence of an optical cavity, the molecular degrees of freedom greatly mix with the cavity mode to form hybrid cavity-matter states called polaritons. These newly generated hybrid light-matter states manipulate the potential energy surfaces and significantly change the reaction dynamics. Here, we theoretically investigate the effects of a strong light-matter interaction on the chemiluminescent reaction of dioxetane using the extended Jaynes-Cummings model. The cavity couplings corresponding to the electronic and vibrational degrees of freedom have been included in the interaction Hamiltonian. We explore how the cavity alters the ground- and excited-state path energy barriers and reaction rates. Our results demonstrate that the formation of excited-state products in the dioxetane decomposition process can be either accelerated or suppressed, depending on the molecular orientation with respect to the cavity polarization.

Beyond Cavity Born-Oppenheimer: On Nonadiabatic Coupling and Effective Ground State Hamiltonians in Vibro-Polaritonic Chemistry

The emerging field of vibro-polaritonic chemistry studies the impact of light-matter hybrid states known as vibrational polaritons on chemical reactivity and molecular properties. Here, we discuss vibro-polaritonic chemistry from a quantum chemical perspective beyond the cavity Born-Oppenheimer (CBO) approximation and examine the role of electron-photon correlation in effective ground state Hamiltonians. We first quantitatively review ab initio vibro-polaritonic chemistry based on the molecular Pauli-Fierz Hamiltonian in dipole approximation and a vibrational strong coupling (VSC) Born-Huang expansion. We then derive nonadiabatic coupling elements arising from both "slow" nuclei and cavity modes compared to "fast" electrons via the generalized Hellmann-Feynman theorem, discuss their properties, and reevaluate the CBO approximation. In the second part, we introduce a crude VSC Born-Huang expansion based on adiabatic electronic states, which provides a foundation for widely employed effective Pauli-Fierz Hamiltonians in ground state vibro-polaritonic chemistry. Those do not strictly respect the CBO approximation but an alternative scheme, which we name crude CBO approximation. We argue that the crude CBO ground state misses electron-photon correlation relative to the CBO ground state due to neglected cavity-induced nonadiabatic transition dipole couplings to excited states. A perturbative connection between both ground state approximations is proposed, which identifies the crude CBO ground state as a first-order approximation to its CBO counterpart. We provide an illustrative numerical analysis of the cavity Shin-Metiu model with a focus on nonadiabatic coupling under VSC and electron-photon correlation effects on classical activation barriers. We finally discuss the potential shortcomings of the electron-polariton Hamiltonian when employed in the VSC regime.

Understanding Polaritonic Chemistry from Ab Initio Quantum Electrodynamics

In this review, we present the theoretical foundations and first-principles frameworks to describe quantum matter within quantum electrodynamics (QED) in the low-energy regime, with a focus on polaritonic chemistry. By starting from fundamental physical and mathematical principles, we first review in great detail ab initio nonrelativistic QED. The resulting Pauli-Fierz quantum field theory serves as a cornerstone for the development of (in principle exact but in practice) approximate computational methods such as quantum-electrodynamical density functional theory, QED coupled cluster, or cavity Born-Oppenheimer molecular dynamics. These methods treat light and matter on equal footing and, at the same time, have the same level of accuracy and reliability as established methods of computational chemistry and electronic structure theory. After an overview of the key ideas behind those ab initio QED methods, we highlight their benefits for understanding photon-induced changes of chemical properties and reactions. Based on results obtained by ab initio QED methods, we identify open theoretical questions and how a so far missing detailed understanding of polaritonic chemistry can be established. We finally give an outlook on future directions within polaritonic chemistry and first-principles QED.

Cavity Born-Oppenheimer Hartree-Fock Ansatz: Light-Matter Properties of Strongly Coupled Molecular Ensembles

Experimental studies indicate that optical cavities can affect chemical reactions through either vibrational or electronic strong coupling and the quantized cavity modes. However, the current understanding of the interplay between molecules and confined light modes is incomplete. Accurate theoretical models that take into account intermolecular interactions to describe ensembles are therefore essential to understand the mechanisms governing polaritonic chemistry. We present an ab initio Hartree-Fock ansatz in the framework of the cavity Born-Oppenheimer approximation and study molecules strongly interacting with an optical cavity. This ansatz provides a nonperturbative, self-consistent description of strongly coupled molecular ensembles, taking into account the cavity-mediated dipole self-energy contributions. To demonstrate the capability of the cavity Born-Oppenheimer Hartree-Fock ansatz, we study the collective effects in ensembles of strongly coupled diatomic hydrogen fluoride molecules. Our results highlight the importance of the cavity-mediated intermolecular dipole-dipole interactions, which lead to energetic changes of individual molecules in the coupled ensemble.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

QED Theory for Controlling the Molecule-Cavity Interaction: From Solvable Analytical Models to Realistic Ones

The study of the interactions of chemical systems in a cavity and the ability to control the reactions inside the cavities become an evolving and hot field of research. Despite that, there is still a significant gap between experiment and theory. Herein, we aim to bridge this gap by starting with the analysis of solvable analytical models for reactions inside a cavity, then continuing to realistic models for many molecules inside a single mode and in a multimode cavity. In addition, we investigate different ways to control the strength of the molecule-cavity coupling term, which in turn allows controlling chemical reactions. Our analysis can benefit the development of ab initio computational methods to simulate molecular systems in polariton cavities; in addition, we show how to parameterize the model Hamiltonians in order to simulate a specific molecular system. Finally, we demonstrate the possibility of achieving isomerization, in case it is prohibited out of the cavity, by placing the reaction inside a cavity.

Theory for Cavity-Modified Ground-State Reactivities via Electron-Photon Interactions

We provide a simple and intuitive theory to explain how coupling a molecule to an optical cavity can modify ground-state chemical reactivity by exploiting intrinsic quantum behaviors of light-matter interactions. Using the recently developed polarized Fock states representation, we demonstrate that the change of the ground-state potential is achieved due to the scaling of diabatic electronic couplings with the overlap of the polarized Fock states. Our theory predicts that for a proton-transfer model system, the ground-state barrier height can be modified through light-matter interactions when the cavity frequency is in the electronic excitation range. Our simple theory explains several recent computational investigations that discovered the same effect. We further demonstrate that under the deep strong coupling limit of the light and matter, the polaritonic ground and first excited eigenstates become the Mulliken-Hush diabatic states, which are the eigenstates of the dipole operator. This work provides a simple but powerful theoretical framework to understand how strong coupling between the molecule and the cavity can modify ground-state reactivities.

Computational study on the catalytic control of endo/exo Diels-Alder reactions by cavity quantum vacuum fluctuations

Achieving control over chemical reaction's rate and stereoselectivity realizes one of the Holy Grails in chemistry that can revolutionize chemical and pharmaceutical industries. Strong light-matter interaction in optical or nanoplasmonic cavities might provide the knob to reach such control. In this work, we demonstrate the catalytic and selectivity control of an optical cavity for two selected Diels-Alder cycloaddition reactions using the quantum electrodynamics coupled cluster (QED-CC) method. Herein, we find that by changing the molecular orientation with respect to the polarization of the cavity mode the reactions can be significantly inhibited or selectively enhanced to produce major endo or exo products on demand. This work highlights the potential of utilizing quantum vacuum fluctuations of an optical cavity to modulate the rate of Diels-Alder cycloaddition reactions and to achieve stereoselectivity in a practical and non-intrusive way. We expect that the present findings will be applicable to a larger set of relevant reactions, including the click chemical reactions.

Ab Initio Molecular Cavity Quantum Electrodynamics Simulations Using Machine Learning Models

We present a mixed quantum-classical simulation of polariton dynamics for molecule-cavity hybrid systems. In particular, we treat the coupled electronic-photonic degrees of freedom (DOFs) as the quantum subsystem and the nuclear DOFs as the classical subsystem and use the trajectory surface hopping approach to simulate non-adiabatic dynamics among the polariton states due to the coupled motion of nuclei. We use the accurate nuclear gradient expression derived from the Pauli-Fierz quantum electrodynamics Hamiltonian without making further approximations. The energies, gradients, and derivative couplings of the molecular systems are obtained from the on-the-fly simulations at the level of complete active space self-consistent field (CASSCF), which are used to compute the polariton energies and nuclear gradients. The derivatives of dipoles are also necessary ingredients in the polariton nuclear gradient expression but are often not readily available in electronic structure methods. To address this challenge, we use a machine learning model with the Kernel ridge regression method to construct the dipoles and further obtain their derivatives, at the same level as the CASSCF theory. The cavity loss process is modeled with the Lindblad jump superoperator on the reduced density of the electronic-photonic quantum subsystem. We investigate the azomethane molecule and its photoinduced isomerization dynamics inside the cavity. Our results show the accuracy of the machine-learned dipoles and their usage in simulating polariton dynamics. Our polariton dynamics results also demonstrate the isomerization reaction of azomethane can be effectively tuned by coupling to an optical cavity and by changing the light-matter coupling strength and the cavity loss rate.

Chiral Polaritonics: Analytical Solutions, Intuition, and Use

Preferential selection of a given enantiomer over its chiral counterpart has become increasingly relevant in the advent of the next era of medical drug design. In parallel, cavity quantum electrodynamics has grown into a solid framework to control energy transfer and chemical reactivity, the latter requiring strong coupling. In this work, we derive an analytical solution to a system of many chiral emitters interacting with a chiral cavity similar to the widely used Tavis-Cummings and Hopfield models of quantum optics. We are able to estimate the discriminating strength of chiral polaritonics, discuss possible future development directions and exciting applications such as elucidating homochirality, and deliver much needed intuition to foster the newly flourishing field of chiral polaritonics.

Polaritonic Huang-Rhys Factor: Basic Concepts and Quantifying Light-Matter Interactions in Media

The Huang-Rhys (HR) factor, a dimensionless factor that characterizes electron-phonon (vibronic) coupling, has been extensively employed to investigate a variety of material properties. In the same spirit, we propose a quantity called the polaritonic HR factor to quantitatively describe the effects of (i) light-matter coupling induced by permanent dipoles and (ii) dipole self-energy. The former leads to polaritonic displacements, while the latter is associated with the electronic coupling shift named reorganization dipole self-coupling. In the framework of macroscopic quantum electrodynamics, our theory can evaluate the polaritonic HR factor, reorganization dipole self-coupling, and modified light-matter coupling strength in an arbitrary dielectric environment without free parameters, whose magnitudes are in good agreement with the previous experimental results. We believe that this study provides a useful perspective on understanding and quantifying light-matter interactions in media.

Vibrational Energy Redistribution and Polaritonic Fermi Resonances in the Strong Coupling Regime

Intramolecular vibrational energy redistribution (IVR) plays a significant role in cavity-modified chemical reaction rates. As such, understanding the fundamental mechanisms by which the cavity modifies the IVR pathways is a fundamental step toward engineering the effect of the confined electromagnetic modes on the outcome of chemical processes. Here we consider an ensemble of M two-mode molecules with intramolecular anharmonic couplings interacting with an infrared cavity mode and consider their quantum dynamics and infrared spectra. Polaritonic Fermi resonances involving fundamental and overtone states of the polaritonic subsystem mediate efficient energy transfer pathways between otherwise off-resonant molecular states. These pathways are of collective nature, yet enabled by the intramolecular anharmonic couplings. Hence, through polaritonic Fermi resonances, cavity excitation can efficiently spread toward low-frequency modes while becoming delocalized over several molecules.

Nonadiabatic Wave Packet Dynamics with Ab Initio Cavity-Born-Oppenheimer Potential Energy Surfaces

Strong coupling of molecules with quantized electromagnetic fields can reshape their potential energy surfaces by forming dressed states. In such a scenario, it is possible to manipulate the dynamics of the molecule and open new photochemical reaction pathways. A theoretical approach to describe such coupled molecular-photon systems is the Cavity-Born-Oppenheimer (CBO) approximation. Similarly to the standard Born-Oppenheimer (BO) approximation, the system is partitioned and the electronic part of the system is treated quantum mechanically. This separation leads to CBO surfaces that depend on both nuclear and photonic coordinates. In this work, we demonstrated, for two molecular examples, how the concept of the CBO approximation can be used to perform nonadiabatic wave packet dynamics of a coupled molecular-cavity system. The light-matter interaction is incorporated in the CBO surfaces and the associated nonadiabatic coupling elements. We show that molecular and cavity contributions can be treated on the same numerical footing. This approach gives a new perspective on the description of light-matter coupling in molecular systems.

Dissociation dynamics of a diatomic molecule in an optical cavity

We study the dissociation dynamics of a diatomic molecule, modeled as a Morse oscillator, coupled to an optical cavity. A marked suppression of the dissociation probability, both classical and quantum, is observed for cavity frequencies significantly below the fundamental transition frequency of the molecule. We show that the suppression in the probability is due to the nonlinearity of the dipole function. The effect can be rationalized entirely in terms of the structures in the classical phase space of the model system.

Shining light on the microscopic resonant mechanism responsible for cavity-mediated chemical reactivity

Strong light-matter interaction in cavity environments is emerging as a promising approach to control chemical reactions in a non-intrusive and efficient manner. The underlying mechanism that distinguishes between steering, accelerating, or decelerating a chemical reaction has, however, remained unclear, hampering progress in this frontier area of research. We leverage quantum-electrodynamical density-functional theory to unveil the microscopic mechanism behind the experimentally observed reduced reaction rate under cavity induced resonant vibrational strong light-matter coupling. We observe multiple resonances and obtain the thus far theoretically elusive but experimentally critical resonant feature for a single strongly coupled molecule undergoing the reaction. While we describe only a single mode and do not explicitly account for collective coupling or intermolecular interactions, the qualitative agreement with experimental measurements suggests that our conclusions can be largely abstracted towards the experimental realization. Specifically, we find that the cavity mode acts as mediator between different vibrational modes. In effect, vibrational energy localized in single bonds that are critical for the reaction is redistributed differently which ultimately inhibits the reaction.