Dynamics and spectroscopy of molecular ensembles in a lossy microcavity

Theory and quantum dynamics simulations of exciton-polariton motional narrowing

The motional narrowing effect has been extensively studied for cavity exciton-polariton systems in recent decades both experimentally and theoretically, which is featured by (1) the subaverage behavior and (2) the asymmetric linewidths for the upper polariton and the lower polariton. However, a minimal theoretical model that is clear and adequate to address all these effects as well as the linewidth scaling relations remains missing. In this work, based on the single mode 1D Holstein-Tavis-Cummings (HTC) model, we studied the motional narrowing effect of the polariton linear absorption spectra via both semi-analytic derivations and numerically exact quantum dynamics simulations using the hierarchical equations of motion approach. The results reveal that under collective light-matter coupling between a cavity mode and N molecules, the polariton linewidth scales as 1/N under the slow limit, while scales as 1/N under the fast limit, due to the polaron decoupling effect. Furthermore, by varying the detunings, the polariton linewidths exhibit significant motional narrowing, covering both characters mentioned above. Our analytic linewidth expressions [Eqs. (34) and (35)] agree well with the numerical exact simulations in all the parameter regimes we explored. These results indicate that the physics of motional narrowing is adequately accounted for by the single-mode 1D HTC model. We envision that both the numerical results and the analytic polariton linewidths expression presented in this work will offer great theoretical value for providing a better understanding of the exciton-polariton motional narrowing based on the HTC model.

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.

Non-Hermitian molecular dynamics simulations of exciton-polaritons in lossy cavities

The observation that materials can change their properties when placed inside or near an optical resonator has sparked a fervid interest in understanding the effects of strong light-matter coupling on molecular dynamics, and several approaches have been proposed to extend the methods of computational chemistry into this regime. Whereas the majority of these approaches have focused on modeling a single molecule coupled to a single cavity mode, changes to chemistry have so far only been observed experimentally when very many molecules are coupled collectively to multiple modes with short lifetimes. While atomistic simulations of many molecules coupled to multiple cavity modes have been performed with semi-classical molecular dynamics, an explicit description of cavity losses has so far been restricted to simulations in which only a very few molecular degrees of freedom were considered. Here, we have implemented an effective non-Hermitian Hamiltonian to explicitly treat cavity losses in large-scale semi-classical molecular dynamics simulations of organic polaritons and used it to perform both mean-field and surface hopping simulations of polariton relaxation, propagation, and energy transfer.

A Quantum Chemistry Approach to Linear Vibro-Polaritonic Infrared Spectra with Perturbative Electron-Photon Correlation

In the vibrational strong coupling (VSC) regime, molecular vibrations and resonant low-frequency cavity modes form light-matter hybrid states, vibrational polaritons, with characteristic infrared (IR) spectroscopic signatures. Here, we introduce a molecular quantum chemistry-based computational scheme for linear IR spectra of vibrational polaritons in polyatomic molecules, which perturbatively accounts for nonresonant electron-photon interactions under VSC. Specifically, we formulate a cavity Born-Oppenheimer perturbation theory (CBO-PT) linear response approach, which provides an approximate but systematic description of such electron-photon correlation effects in VSC scenarios while relying on molecular ab initio quantum chemistry methods. We identify relevant electron-photon correlation effects at the second order of CBO-PT, which manifest as static polarizability-dependent Hessian corrections and an emerging polarizability-dependent cavity intensity component providing access to transmission spectra commonly measured in vibro-polaritonic chemistry. Illustratively, we address electron-photon correlation effects perturbatively in IR spectra of CO2 and Fe(CO)5 vibro-polaritonic models in sound agreement with nonperturbative CBO linear response theory.

Finite temperature dynamics of the Holstein-Tavis-Cummings model

By employing the numerically accurate multiple Davydov Ansatz (mDA) formalism in combination with the thermo-field dynamics (TFD) representation of quantum mechanics, we systematically explore the influence of three parameters-temperature, photonic-mode detuning, and qubit-phonon coupling-on population dynamics and absorption spectra of the Holstein-Tavis-Cummings (HTC) model. It is found that elevated qubit-phonon couplings and/or temperatures have a similar impact on all dynamic observables: they suppress the amplitudes of Rabi oscillations in photonic populations as well as broaden the peaks and decrease their intensities in the absorption spectra. Our results unequivocally demonstrate that the HTC dynamics is very sensitive to the concerted variation of the three aforementioned parameters, and this finding can be used for fine-tuning polaritonic transport. The developed mDA-TFD methodology can be efficiently applied for modeling, predicting, optimizing, and comprehensively understanding dynamic and spectroscopic responses of actual molecular systems in microcavities.

Shaping the laser control landscape of a hydrogen transfer reaction by vibrational strong coupling. A direct optimal control approach

Controlling molecular reactivity by shaped laser pulses is a long-standing goal in chemistry. Here, we suggest a direct optimal control approach that combines external pulse optimization with other control parameters arising in the upcoming field of vibro-polaritonic chemistry for enhanced controllability. The direct optimal control approach is characterized by a simultaneous simulation and optimization paradigm, meaning that the equations of motion are discretized and converted into a set of holonomic constraints for a nonlinear optimization problem given by the control functional. Compared with indirect optimal control, this procedure offers great flexibility, such as final time or Hamiltonian parameter optimization. A simultaneous direct optimal control theory will be applied to a model system describing H-atom transfer in a lossy Fabry-Pérot cavity under vibrational strong coupling conditions. Specifically, optimization of the cavity coupling strength and, thus, of the control landscape will be demonstrated.

Coupling Molecular Systems with Plasmonic Nanocavities: A Quantum Dynamics Approach

Plasmonic nanoparticles have the capacity to confine electromagnetic fields to the subwavelength regime and provide strong coupling with few or even a single emitter at room temperature. The photophysical properties of the emitters are highly dependent on the relative distance and orientation between them and the nanocavity. Therefore, there is a need for accurate and general light-matter interaction models capable of guiding their design in application-oriented devices. In this work, we present a Hermitian formalism within the framework of quantum dynamics and based on first-principles electronic structure calculations. Our vibronic approach considers the quantum nature of the plasmonic excitations and the dynamics of nonradiative channels to model plasmonic nanocavities and their dipolar coupling to molecular electronic states. Thus, the quantized and dissipative nature of the nanocavity is fully addressed.

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.

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.

The role of dephasing for dark state coupling in a molecular Tavis-Cummings model

The collective coupling of an ensemble of molecules to a light field is commonly described by the Tavis-Cummings model. This model includes numerous eigenstates that are optically decoupled from the optically bright polariton states. Accessing these dark states requires breaking the symmetry in the corresponding Hamiltonian. In this paper, we investigate the influence of non-unitary processes on the dark state dynamics in the molecular Tavis-Cummings model. The system is modeled with a Lindblad equation that includes pure dephasing, as it would be caused by weak interactions with an environment, and photon decay. Our simulations show that the rate of pure dephasing, as well as the number of two-level systems, has a significant influence on the dark state population.

Certifying Multimode Light-Matter Interaction in Lossy Resonators

Quantum models based on few-mode master equations have been a central tool in the study of resonator quantum electrodynamics, extending the seminal single-mode Jaynes-Cummings model to include loss and multiple modes. Despite their broad application range, previous approaches within this framework have either relied on a Markov approximation or a fitting procedure. By combining ideas from pseudomode and quasinormal mode theory, we develop a certification criterion for multi-mode effects in lossy resonators. It is based on a witness observable, and neither requires a fitting procedure nor a Markov approximation. Using the resulting criterion, we demonstrate that such multi-mode effects are important for understanding previous experiments in x-ray cavity QED with Mössbauer nuclei and that they allow one to tune the nuclear ensemble properties.

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.

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.

Multidimensional Quantum Dynamical Simulation of Infrared Spectra under Polaritonic Vibrational Strong Coupling

Recent experimental and theoretical studies demonstrate that the chemical reactivity of molecules can be modified inside an optical cavity. Here, we provide a theoretical framework for conducting multidimensional quantum simulations of the infrared (IR) spectra for molecules interacting with cavity modes. A single water molecule under polaritonic vibrational strong coupling serves as an illustrative example. Combined with accurate potential energy and dipole moment surfaces, our cavity vibrational self-consistent field/virtual state configuration interaction (cav-VSCF/VCI) approach can predict the IR spectra when the molecule is inside or outside the cavity. The spectral signatures of Rabi splittings and shifts of certain bands are found to be strongly dependent on the frequency and polarization direction of the cavity modes. Analyses of the simulated spectra show that polaritonic vibrational strong coupling can induce unconventional couplings among the molecule's vibrational modes, suggesting that intramolecular vibrational energy transfer can be significantly accelerated by the cavity.

Suppressing non-radiative decay of photochromic organic molecular systems in the strong coupling regime

The lifetimes of electronic excited states have a strong influence on the efficiency of organic solar cells. However, in some molecular systems a given excited state lifetime is reduced due to the non-radiative decay through conical intersections. Several strategies may be used to suppress this decay channel. The use of the strong light-matter coupling provided in optical nano-cavities is the focus of this paper. Here, we consider the meso-tert-butyl-4,4-difluoro-4-bora-3a,4a-diaza-s-indacene molecule (meso-tert-butyl-BODIPY) as a showcase of how strong and ultrastrong coupling might help in the development of organic solar cells. The meso-tert-butyl-BODIPY is known for its low fluorescence yield caused by the non-radiative decay through a conical intersection. However, we show here that, by considering this system within a cavity, the strong coupling can lead to significant changes in the multidimensional landscape of the potential energy surfaces of meso-tert-butyl-BODIPY, suppressing almost completely the decay of the excited state wave packet back to the ground state. By means of multi configuration electronic structure calculations and nuclear wave packet dynamics, the coupling with the cavity is analyzed in-depth to provide further insight of the interaction. By fine-tuning the cavity field strength and resonance frequency, we show that one can change the nuclear dynamics in the excited state, and control the non-radiative decay. This may lead to a faster and more efficient population transfer or the suppression of it.

Cavity-induced non-adiabatic dynamics and spectroscopy of molecular rovibrational polaritons studied by multi-mode quantum models

We study theoretically the quantum dynamics and spectroscopy of rovibrational polaritons formed in a model system composed of a single rovibrating diatomic molecule, which interacts with two degenerate, orthogonally polarized modes of an optical Fabry–Pérot cavity. We employ an effective rovibrational Pauli–Fierz Hamiltonian in length gauge representation and identify three-state vibro-polaritonic conical intersections (VPCIs) between singly excited vibro-polaritonic states in a two-dimensional angular coordinate branching space. The lower and upper vibrational polaritons are of mixed light–matter hybrid character, whereas the intermediate state is purely photonic in nature. The VPCIs provide effective population transfer channels between singly excited vibrational polaritons, which manifest in rich interference patterns in rotational densities. Spectroscopically, three bright singly excited states are identified when an external infrared laser field couples to both a molecular and a cavity mode. The non-trivial VPCI topology manifests as pronounced multi-peak progression in the spectral region of the upper vibrational polariton, which is traced back to the emergence of rovibro-polaritonic light–matter hybrid states. Experimentally, ubiquitous spontaneous emission from cavity modes induces a dissipative reduction of intensity and peak broadening, which mainly influences the purely photonic intermediate state peak as well as the rovibro-polaritonic progression.

Incorporating Lindblad Decay Dynamics into Mixed Quantum-Classical Simulations

We derive the $\mathcal{L}$-MFE method to incorporate Lindblad jump operator dynamics into the mean-field Ehrenfest (MFE) approach. We map the density matrix evolution of Lindblad dynamics onto pure state coefficients using trajectory averages. We use simple assumptions to construct the $\mathcal{L}$-MFE method that satisfies this exact mapping. This establishes a method that uses independent trajectories which exactly reproduces Lindblad decay dynamics using a wavefunction description, with deterministic changes of the magnitudes of the quantum expansion coefficients, while only adding on a stochastic phase. We further demonstrate that when including nuclei in the Ehrenfest dynamics, the $\mathcal{L}$-MFE method gives semi-quantitatively accurate results, with the accuracy limited by the accuracy of the approximations present in the semiclassical MFE approach. This work provides a general framework to incorporate Lindblad dynamics into semiclassical or mixed quantum-classical simulations.

Accurate Simulation of Spectroscopic Signatures of Cavity-Assisted, Conical-Intersection-Controlled Singlet Fission Processes

A numerically accurate, fully quantum methodology has been developed for the simulation of the dynamics and nonlinear spectroscopic signals of cavity-assisted, conical-intersection-controlled singlet fission systems. The methodology is capable of handling several molecular systems strongly coupled to the photonic mode of the cavity and treats the intrinsic conical intersection and cavity-induced polaritonic conical intersections in a numerically exact manner. Contributions of higher-lying molecular electronic states are accounted for comprehensively. The intriguing process of cavity-modified fission dynamics, including all of its electronic, vibrational, and photonic degrees of freedom, together with its two-dimensional spectroscopic manifestation, is simulated for two rubrene dimers strongly coupled to the cavity mode.

Theoretical Challenges in Polaritonic Chemistry

Polaritonic chemistry exploits strong light-matter coupling between molecules and confined electromagnetic field modes to enable new chemical reactivities. In systems displaying this functionality, the choice of the cavity determines both the confinement of the electromagnetic field and the number of molecules that are involved in the process. While in wavelength-scale optical cavities the light-matter interaction is ruled by collective effects, plasmonic subwavelength nanocavities allow even single molecules to reach strong coupling. Due to these very distinct situations, a multiscale theoretical toolbox is then required to explore the rich phenomenology of polaritonic chemistry. Within this framework, each component of the system (molecules and electromagnetic modes) needs to be treated in sufficient detail to obtain reliable results. Starting from the very general aspects of light-molecule interactions in typical experimental setups, we underline the basic concepts that should be taken into account when operating in this new area of research. Building on these considerations, we then provide a map of the theoretical tools already available to tackle chemical applications of molecular polaritons at different scales. Throughout the discussion, we draw attention to both the successes and the challenges still ahead in the theoretical description of polaritonic chemistry.

Not dark yet for strong light-matter coupling to accelerate singlet fission dynamics

Polaritons are unique hybrid light-matter states that offer an alternative way to manipulate chemical processes. In this work, we show that singlet fission dynamics can be accelerated under strong light-matter coupling. For superexchange-mediated singlet fission, state mixing speeds up the dynamics in cavities when the lower polariton is close in energy to the multiexcitonic state. This effect is more pronounced in non-conventional singlet fission materials in which the energy gap between the bright singlet exciton and the multiexcitonic state is large ( > 0.1 eV). In this case, the dynamics is dominated by the polaritonic modes and not by the bare-molecule-like dark states, and, additionally, the resonant enhancement due to strong coupling is robust even for energetically broad molecular states. The present results provide a new strategy to expand the range of suitable materials for efficient singlet fission by making use of strong light-matter coupling.

Non-Hermitian Cavity Quantum Electrodynamics - Configuration Interaction Singles Approach for Polaritonic Structure with ab initio Molecular Hamiltonians

We combine ab initio molecular electronic Hamiltonians with a cavity quantum electrodynamics model for dissipative photonic modes and apply mean-field theories to the ground- and excited-states of resulting polaritonic systems. In particular, we develop a non-Hermitian configuration interaction singles theory for mean-field ground- and excited-states of the molecular system strongly interacting with a photonic mode, and apply these methods to elucidating the phenomenology of paradigmatic polaritonic systems. We leverage the Psi4Numpy framework to yield open-source and accessible reference implementations of these methods.

Simple but accurate estimation of light-matter coupling strength and optical loss for a molecular emitter coupled with photonic modes

Light-matter coupling strength and optical loss are two key physical quantities in cavity quantum electrodynamics (CQED), and their interplay determines whether light-matter hybrid states can be formed or not in chemical systems. In this study, by using macroscopic quantum electrodynamics (MQED) combined with a pseudomode approach, we present a simple but accurate method, which allows us to quickly estimate the light-matter coupling strength and optical loss without free parameters. Moreover, for a molecular emitter coupled with photonic modes (including cavity modes and plasmon polariton modes), we analytically and numerically prove that the dynamics derived from the MQED-based wavefunction approach is mathematically equivalent to the dynamics governed by the CQED-based Lindblad master equation when the Purcell factor behaves like Lorentzian functions.

Ground state properties and infrared spectra of anharmonic vibrational polaritons of small molecules in cavities

Recent experiments and theory suggest that ground state properties and reactivity of molecules can be modified when placed inside a nanoscale cavity, giving rise to strong coupling between vibrational modes and the quantized cavity field. This is commonly thought to be caused either by a cavity-distorted Born-Oppenheimer ground state potential or by the formation of light-matter hybrid states, vibrational polaritons. Here, we systematically study the effect of a cavity on ground state properties and infrared spectra of single molecules, considering vibration-cavity coupling strengths from zero up to the vibrational ultrastrong coupling regime. Using single-mode models for Li-H and O-H stretch modes and for the NH₃ inversion mode, respectively, a single cavity mode in resonance with vibrational transitions is coupled to position-dependent molecular dipole functions. We address the influence of the cavity mode on polariton ground state energies, equilibrium bond lengths, dissociation energies, activation energies for isomerization, and on vibro-polaritonic infrared spectra. In agreement with earlier work, we observe all mentioned properties being strongly affected by the cavity, but only if the dipole self-energy contribution in the interaction Hamiltonian is neglected. When this term is included, these properties do not depend significantly on the coupling anymore. Vibro-polaritonic infrared spectra, in contrast, are always affected by the cavity mode due to the formation of excited vibrational polaritons. It is argued that the quantized nature of vibrational polaritons is key to not only interpreting molecular spectra in cavities but also understanding the experimentally observed modification of molecular reactivity in cavities.

Polar diatomic molecules in optical cavities: Photon scaling, rotational effects, and comparison with classical fields

We address topics related to molecules coupled to quantum radiation. The formalism of light-matter interaction is different for classical and quantum fields, but some analogies remain, such as the formation of light induced crossings. We show that under particular circumstances, the molecular dynamics under quantum or classical fields produce similar results, as long as the radiation is prepared as a Fock state and far from ultra-strong coupling regimes. At this point, the choice of specific initial Fock states is irrelevant since the dynamics scales. However, in realistic multistate molecular systems, radiative scaling may fail due to the presence of simultaneous efficient non-radiative couplings in the dynamics. Polar molecules have permanent dipoles, and within the context of the full quantum Rabi model with a Pauli-Fierz Hamiltonian, they play a crucial role in the polaritonic dynamics since both permanent dipole moments and self-energy terms produce drastic changes on the undressed potential energy surfaces at high coupling strengths. We also gauge the effect of including rotational degrees of freedom in cavity molecular photodynamics. For diatomic molecules, the addition of rotation amounts to transform (both with classical or quantum fields) a light induced crossing into a light induced conical intersection. However, we show that conical intersections due to molecular rotation do not represent the standard properties of well-known efficient intrinsic conical intersections inasmuch they do not enhance the quantum transition rates.

Optical-Cavity Manipulation of Conical Intersections and Singlet Fission in Pentacene Dimers

We demonstrate how the singlet fission process in pentacene dimers mediated by a conical intersection is controlled by coupling the molecule to a confined optical cavity photon mode. By following the polariton quantum dynamics of a conical intersection coupled to a cavity mode taking into account vibrational relaxation and cavity loss, we find that the singlet fission can be significantly suppressed because the polaritonic conical intersection is pushed away from the initial Franck-Condon excitation region.

Nonequilibrium effects of cavity leakage and vibrational dissipation in thermally activated polariton chemistry

In vibrational strong coupling (VSC), molecular vibrations strongly interact with the modes of an optical cavity to form hybrid light-matter states known as vibrational polaritons. Experiments show that the kinetics of thermally activated chemical reactions can be modified by VSC. Transition-state theory, which assumes that internal thermalization is fast compared to reactive transitions, has been unable to explain the observed findings. Here, we carry out kinetic simulations to understand how dissipative processes, namely, those introduced by VSC to the chemical system, affect reactions where internal thermalization and reactive transitions occur on similar timescales. Using the Marcus-Levich-Jortner type of electron transfer as a model reaction, we show that such dissipation can change reactivity by accelerating internal thermalization, thereby suppressing nonequilibrium effects that occur in the reaction outside the cavity. This phenomenon is attributed mainly to cavity decay (i.e., photon leakage), but a supporting role is played by the relaxation between polaritons and dark states. When nonequilibrium effects are already suppressed in the bare reaction (the reactive species are essentially at internal thermal equilibrium throughout the reaction), we find that reactivity does not change significantly under VSC. Connections are made between our results and experimental observations.

Controlling the Photostability of Pyrrole with Optical Nanocavities

Strong light-matter coupling provides a new strategy to manipulate the non-adiabatic dynamics of molecules by modifying potential energy surfaces. The vacuum field of nanocavities can couple strongly with the molecular degrees of freedom and form hybrid light-matter states, termed as polaritons or dressed states. The photochemistry of molecules possessing intrinsic conical intersections can be significantly altered by introducing cavity couplings to create new conical intersections or avoided crossings. Here, we explore the effects of optical cavities on the photo-induced hydrogen elimination reaction of pyrrole. Wave packet dynamics simulations have been performed on the two-state, two-mode model of pyrrole, combined with the cavity photon mode. Our results show how the optical cavities assist in controlling the photostability of pyrrole and influence the reaction mechanism by providing alternative dissociation pathways. The cavity effects have been found to be intensely dependent on the resonance frequency. We further demonstrate the importance of the vibrational cavity couplings and dipole-self interaction terms in describing the cavity-modified non-adiabatic dynamics.

Many-photon excitation of organic molecules in a cavity-Superradiance as a measure of coherence

Coherent excitation of a molecular ensemble coupled to a common radiation mode can lead to the collective emission of radiation known as superradiance. This collective emission only occurs if there is an entanglement between the molecules in their ground and excited state and can, therefore, serve as a macroscopic measure of coherence in the ensemble. Reported here are wave packet propagations for various pyrazine models of increasing complexity and molecular ensembles thereof. We show that ensemble coherence upon photoexcitation can prevail up to relatively long time scales although the effect can diminish quickly with increasing ensemble size. Coherence can also build up over time and even reemerge after the molecules have passed through a conical intersection. The effect of the pump pulse characteristics on the collective response of the molecular ensemble is also studied. A broadband pulse imprints a large amount of initial coherence to the system, as compared to a longer pulse with a smaller spread in the frequency domain. However, the differential effects arising from a different pulse duration and coherent bandwidth become less prominent if the emission of light from the ensemble takes place after a non-adiabatic decay process.

Simulating photodissociation reactions in bad cavities with the Lindblad equation

Optical cavities, e.g., as used in organic polariton experiments, often employ low finesse mirrors or plasmonic structures. The photon lifetime in these setups is comparable to the timescale of the nuclear dynamics governing the photochemistry. This highlights the need for including the effect of dissipation in the molecular simulations. In this study, we perform wave packet dynamics with the Lindblad master equation to study the effect of a finite photon lifetime on the dissociation of the MgH⁺ molecule model system. Photon lifetimes of several different orders of magnitude are considered to encompass an ample range of effects inherent to lossy cavities.

Role of Cavity Losses on Nonadiabatic Couplings and Dynamics in Polaritonic Chemistry

We present a non-Hermitian formulation of the polaritonic structure of azobenzene strongly coupled to a photonic mode that explicitly accounts for the fleeting nature of the photon-molecule interaction. This formalism reveals that the polaritonic nonadiabatic couplings that facilitate cis-trans isomerization can be dramatically modified by photonic dissipation. We perform Fewest-Switches Surface Hopping dynamics on the surfaces that derive from our non-Hermitian formalism and find that the polaritonic isomerization yields are strongly suppressed for moderate dissipation rates and that cavity-free isomerization dynamics are recovered under large dissipation rates. These findings highlight the important role that the finite lifetime of photonic degrees of freedom play in polaritonic chemistry.

Photoprotecting Uracil by Coupling with Lossy Nanocavities

We analyze how the photorelaxation dynamics of a molecule can be controlled by modifying its electromagnetic environment using a nanocavity mode. In particular, we consider the photorelaxation of the RNA nucleobase uracil, which is the natural mechanism to prevent photodamage. In our theoretical work, we identify the operative conditions in which strong coupling with the cavity mode can open an efficient photoprotective channel, resulting in a relaxation dynamics twice as fast as the natural one. We rely on a state-of-the-art chemically detailed molecular model and a non-Hermitian Hamiltonian propagation approach to perform full-quantum simulations of the system dissipative dynamics. By focusing on the photon decay, our analysis unveils the active role played by cavity-induced dissipative processes in modifying chemical reaction rates, in the context of molecular polaritonics. Remarkably, we find that the photorelaxation efficiency is maximized when an optimal trade-off between light-matter coupling strength and photon decay rate is satisfied. This result is in contrast with the common intuition that increasing the quality factor of nanocavities and plasmonic devices improves their performance. Finally, we use a detailed model of a metal nanoparticle to show that the speedup of the uracil relaxation could be observed via coupling with a nanosphere pseudomode, without requiring the implementation of complex nanophotonic structures.