Atoms and molecules in cavities, from weak to strong coupling in quantum-electrodynamics (QED) chemistry

Nonadiabatic Field: A Conceptually Novel Approach for Nonadiabatic Quantum Molecular Dynamics

Reliable trajectory-based nonadiabatic quantum dynamics methods at the atomic/molecular level are critical for the practical understanding and rational design of many important processes in real large/complex systems, where the quantum dynamical behavior of electrons and that of nuclei are coupled. The paper reports latest progress of nonadiabatic field (NaF), a conceptually novel approach for nonadiabatic quantum dynamics with independent trajectories. Substantially different from the mainstreams of Ehrenfest-like dynamics and surface hopping methods, the nuclear force in NaF involves the nonadiabatic force arising from the nonadiabatic coupling between different electronic states, in addition to the adiabatic force contributed by a single adiabatic electronic state. NaF is capable of faithfully describing the interplay between electronic and nuclear motion in a broad regime, which covers where the relevant electronic states keep coupled in a wide range or all the time and where the bifurcation characteristic of nuclear motion is essential. NaF is derived from the exact generalized phase space formulation with coordinate-momentum variables, where constraint phase space (CPS) is employed for discrete electronic-state degrees of freedom (DOFs) and infinite Wigner phase space is used for continuous nuclear DOFs. We propose efficient integrators for the equations of motion of NaF in both adiabatic and diabatic representations. Since the formalism in the CPS formulation is not unique, NaF can in principle be implemented with various phase space representations of the time correlation function (TCF) for the time-dependent property. They are applied to a suite of representative gas-phase and condensed-phase benchmark models where numerically exact results are available for comparison. It is shown that NaF is relatively insensitive to the phase space representation of the electronic TCF and will be a potential tool for practical and reliable simulations of the quantum mechanical behavior of both electronic and nuclear dynamics of nonadiabatic transition processes in real systems.

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.

Polariton-induced Purcell effects via a reduced semiclassical electrodynamics approach

Recent experiments have demonstrated that polariton formation provides a novel strategy for modifying local molecular processes when a large ensemble of molecules is confined within an optical cavity. Herein, a numerical strategy based on coupled Maxwell-Schrödinger equations is examined for simulating local molecular processes in a realistic cavity structure under collective strong coupling. In this approach, only a few molecules, referred to as quantum impurities, are treated quantum mechanically, while the remaining macroscopic molecular layer and the cavity structure are modeled using dielectric functions. When a single electronic two-level system embedded in a Lorentz medium is confined in a two-dimensional Bragg resonator, our numerical simulations reveal a polariton-induced Purcell effect: the radiative decay rate of the quantum impurity is significantly enhanced by the cavity when the impurity frequency matches the polariton frequency, while the rate can sometimes be greatly suppressed when the impurity is near resonance with the bulk molecules forming strong coupling. In addition, this approach demonstrates that the cavity absorption of light exhibits Rabi-splitting-dependent suppression due to the inclusion of a realistic cavity structure. Our simulations also identify a fundamental limitation of this approach-an inaccurate description of polariton dephasing rates into dark modes. This arises because the dark-mode degrees of freedom are not explicitly included when most molecules are modeled using simple dielectric functions. As the polariton-induced Purcell effect alters molecular radiative decay differently from the Purcell effect under weak coupling, this polariton-induced effect may facilitate understanding the origin of polariton-modified photochemistry under electronic strong coupling.

Trajectory-based non-adiabatic simulations of the polariton relaxation dynamics

We benchmark the accuracy of various trajectory-based non-adiabatic methods in simulating the polariton relaxation dynamics under the collective coupling regime. The Holstein-Tavis-Cummings Hamiltonian is used to describe the hybrid light-matter system of N molecules coupled to a single cavity mode. We apply various recently developed trajectory-based methods to simulate the population relaxation dynamics by initially exciting the upper polariton state and benchmark the results against populations computed from exact quantum dynamical propagation using the hierarchical equations of motion approach. In these benchmarks, we have systematically varied the number of molecules N, light-matter detunings, and the light-matter coupling strengths. Our results demonstrate that the symmetrical quasi-classical method with γ correction and spin-mapping linearized semi-classical approaches yield more accurate polariton population dynamics than traditional mixed quantum-classical methods, such as the Ehrenfest and surface hopping techniques.

Stochastic resonance in vibrational polariton chemistry

In this work, we systematically investigate the impact of ambient noise intensity on the rate modifications of ground-state chemical reactions in an optical cavity under vibrational strong-coupling conditions. To achieve this, we utilize a numerically exact open quantum system approach-the hierarchical equations of motion in twin space, combined with a flexible tree tensor network state solver. Our findings reveal a stochastic resonance phenomenon in cavity-modified chemical reactivities: an optimal reaction rate enhancement occurs at an intermediate noise level. In other words, this enhancement diminishes if ambient noise, sensed by the cavity-molecule system through cavity leakage, is either too weak or excessively strong. In the collective coupling regime, when the cavity is weakly damped, rate enhancement strengthens as more molecules couple to the cavity. In contrast, under strong cavity damping, reaction rates decline as the number of molecules grows.

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.

Non-adiabatic quantum electrodynamic effects on electron-nucleus-photon systems: Single photonic mode vs infinite photonic modes

The quantum-electrodynamic non-adiabatic emission (QED-NAE) is a type of radiatively assisted vibronic de-excitation due to electromagnetic vacuum fluctuations on non-adiabatic processes. Building on our previous work [Tsai et al., J. Phys. Chem. Lett. 14, 5924 (2023)], we extend the theory of the QED-NAE rate from a single cavity photonic mode to infinite photonic modes and calculate the QED-NAE rates of 9-cyanoanthracene at the first-principles level. To avoid the confusion, the quantum electrodynamic internal conversion process is renamed as "QED-NAE" in our present work. According to our theory, we identify three key factors influencing the QED-NAE processes: light-matter coupling strength (mode volume), mass-weighted orientation factor, and photonic density of states. The mode volume is the primary factor causing rate differences between the two scenarios. In a single cavity with a small mode volume, strong light-matter coupling strength boosts QED-NAE rates. In contrast, in free space with infinite photonic modes, weak coupling strength significantly reduces these rates. From a single cavity photonic mode to infinite photonic modes, the mass-weighted orientation factor only causes an 8π/3-fold increase in the QED-NAE rate. In free space, the photonic density of state exhibits a flat and quadratic distribution, which slightly reduces the QED-NAE rate. Our study shows that cavities can significantly enhance non-adiabatic QED effects while providing a robust analysis demonstrating that QED vibronic effects can be safely ignored in free space.

Analytical derivative approaches for vibro-polaritonic structures and properties. I. Formalism and implementation

Vibro-polaritons are hybrid light-matter states that arise from the strong coupling between the molecular vibrational transitions and the photons in an optical cavity. Developing theoretical and computational methods to describe and predict the unique properties of vibro-polaritons is of great significance for guiding the design of new materials and experiments. Here, we present the ab initio cavity Born-Oppenheimer density functional theory (CBO-DFT) and formulate the analytic energy gradient and Hessian as well as the nuclear and photonic derivatives of dipole and polarizability within the framework of CBO-DFT to efficiently calculate the harmonic vibrational frequencies, infrared absorption, and Raman scattering spectra of vibro-polaritons as well as to explore the critical points on the cavity potential energy surface. The implementation of analytic derivatives into the electronic structure package is validated by a comparison with the finite-difference method and with other reported computational results. By adopting appropriate exchange-correlation functionals, CBO-DFT can better describe the structure and properties of molecules in the cavity than CBO-Hartree-Fock method. It is expected that CBO-DFT is a useful tool for studying the polaritonic structures and properties.

Toward Collective Chemistry under Strong Light-Matter Coupling

Collective strong light-matter coupling provides a versatile means to manipulate physicochemical properties of molecules and materials. Understanding collective polaritonic dynamics is hindered by the macroscopic number of molecules interacting collectively with photonic modes. We develop a many-body theory to investigate the spectroscopy and dynamics of a molecular ensemble embedded in an optical cavity in the collective strong coupling regime. This theory is constructed by a pseudoparticle representation of the molecular Hamiltonian, which maps the polaritonic Hamiltonian into a coupled fermion-boson model under particle number constraints. The mapped model is then analyzed using the nonequilibrium Green's function theory with the self-energy diagrams identified through a large N expansion. We demonstrate that in the thermodynamic limit, the necessary condition to have any collective effects is to have a macroscopic cavity field. Numerical illustrations are shown for the driven Tavis-Cummings model, which shows an excellent agreement with exact results.

Chiral polaritonics: cavity-mediated enantioselective excitation condensation

Separation of the two mirror images of a chiral molecule, the enantiomers, is a historically complicated problem of major relevance for biological systems. Since chiral molecules are optically active, it has been speculated that strong coupling to circularly polarized fields may be used as a general procedure to unlock enantiospecific reactions. In this work, we focus on how chiral cavities can be used to drive asymmetry in the photochemistry of chiral molecular systems. We first show that strong coupling to circularly polarized fields leads to enantiospecific Rabi splittings, an effect that displays a collective behavior in line with other strong coupling phenomena. Additionally, entanglement with circularly polarized light generates an asymmetry in the enantiomer population of the polaritons, leading to a condensation of the excitation on a preferred molecular configuration. These results confirm that chiral cavities represent a tantalizing opportunity to drive asymmetric photochemistry in enantiomeric mixtures.

Multiple interacting photonic modes in strongly coupled organic microcavities

Room-temperature cavity quantum electrodynamics with molecular materials in optical cavities offers exciting prospects for controlling electronic, nuclear and photonic degrees of freedom for applications in physics, chemistry and materials science. However, achieving strong coupling with molecular ensembles typically requires high molecular densities and substantial electromagnetic-field confinement. These conditions usually involve a significant degree of molecular disorder and a highly structured photonic density of states. It remains unclear to what extent these additional complexities modify the usual physical picture of strong coupling developed for atoms and inorganic semiconductors. Using a microscopic quantum description of molecular ensembles in realistic multimode optical resonators, we show that the emergence of vacuum Rabi splitting in linear spectroscopy is a necessary but not sufficient metric of coherent admixing between light and matter. In low-finesse multi-mode situations, we find that molecular dipoles can be partially hybridized with photonic dissipation channels associated with off-resonant cavity modes. These vacuum-induced dissipative processes ultimately limit the extent of light-matter coherence that the system can sustain.This article is part of the theme issue 'The quantum theory of light'.

Theory of molecular emission power spectra. III. Non-Hermitian interactions in multichromophoric systems coupled with polaritons

Based on our previous study [Wang et al., J. Chem. Phys. 153, 184102 (2020)], we generalize the theory of molecular emission power spectra (EPS) from one molecule to multichromophoric systems in the framework of macroscopic quantum electrodynamics. This generalized theory is applicable to ensembles of molecules, providing a comprehensive description of the molecular spontaneous emission spectrum in arbitrary inhomogeneous, dispersive, and absorbing media. In the far-field region, the analytical formula of EPS can be expressed as the product of a lineshape function (LF) and an electromagnetic environment factor (EEF). To demonstrate the polaritonic effect on multichromophoric systems, we simulate the LF and EEF for one to three molecules weakly coupled to surface plasmon polaritons above a silver surface. Our analytical expressions show that the peak broadening originates from not only the spontaneous emission rates but also the imaginary part of resonant dipole-dipole interactions (non-Hermitian interactions), which is associated with the superradiance of molecular aggregates, indicating that the superradiance rate can be controlled through an intermolecular distance and the design of dielectric environments. This study presents an alternative approach to directly analyze the hybrid-state dynamics of multichromophoric systems coupled with polaritons.

Cavity-enhanced superconductivity in MgB2 from first-principles quantum electrodynamics (QEDFT)

Strong laser pulses can control superconductivity, inducing nonequilibrium transient pairing by leveraging strong-light matter interaction. Here, we demonstrate theoretically that equilibrium ground-state phonon-mediated superconductive pairing can be affected through the vacuum fluctuating electromagnetic field in a cavity. Using the recently developed ab initio quantum electrodynamical density-functional theory approximation, we specifically investigate the phonon-mediated superconductive behavior of MgB[Formula: see text] under different cavity setups and find that in the strong light-matter coupling regime its superconducting transition temperature T[Formula: see text] can be enhanced at most by [Formula: see text]10% in an in-plane (or out-of-plane) polarized and realistic cavity via photon vacuum fluctuations. The results highlight that strong light-matter coupling in extended systems can profoundly alter material properties in a nonperturbative way by modifying their electronic structure and phononic dispersion at the same time. Our findings indicate a pathway to the experimental realization of light-controlled superconductivity in solid-state materials at equilibrium via cavity materials engineering.

Polaritonic Chemistry Using the Density Matrix Renormalization Group Method

The emerging field of polaritonic chemistry explores the behavior of molecules under strong coupling with cavity modes. Despite recent developments in ab initio polaritonic methods for simulating polaritonic chemistry under electronic strong coupling, their capabilities are limited, especially in cases where the molecule also features strong electronic correlation. To bridge this gap, we have developed a novel method for cavity QED calculations utilizing the Density Matrix Renormalization Group (DMRG) algorithm in conjunction with the Pauli-Fierz Hamiltonian. Our approach is applied to investigate the effect of the cavity on the S0-S1 transition of n-oligoacenes, with n ranging from 2 to 5, encompassing 22 fully correlated π orbitals in the largest pentacene molecule. Our findings indicate that the influence of the cavity intensifies with larger acenes. Additionally, we demonstrate that, unlike the full determinantal representation, DMRG efficiently optimizes and eliminates excess photonic degrees of freedom, resulting in an asymptotically constant computational cost as the photonic basis increases.

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.

Strong-Coupling Modification of Singlet-Fission Dynamical Pathways

We investigate theoretically the influence of strong light-matter coupling on the initial steps of the phototriggered singlet-fission process. In particular we focus on intramolecular singlet fission in a TIPS-pentacene dimer derivative described by a vibronic Hamiltonian including the optically active singlet excited states, doubly excited and charge transfer states, as well as the final triplet-triplet pair state. Quantum dynamics simulations of up to four dimers in the cavity indicate that the modified resonance condition imposed by the cavity strongly quenches the passage through the intermediate charge transfer and double-excitation states, thus largely reducing the triplet-triplet yield in the bare system. Subsequently, we modify the system parameters and construct a model Hamiltonian where the optically active singlet excitation lies below the final triplet-triplet state such that the yield of the bare system becomes insignificant. In this case we find that using the upper polariton as the doorway state for photoexcitation can lead to a much enhanced yield. This pathway is operative provided that the system is sufficiently rigid to prevent vibronic losses from the upper polariton to the dark-states manifold.

Simulating Cavity-Modified Electron Transfer Dynamics on NISQ Computers

We present an algorithm based on the quantum-mechanically exact tensor-train thermo-field dynamics (TT-TFD) method for simulating cavity-modified electron transfer dynamics on noisy intermediate-scale quantum (NISQ) computers. The utility and accuracy of the proposed methodology is demonstrated on a model for the photoinduced intramolecular electron transfer reaction within the carotenoid-porphyrin-C60 molecular triad in tetrahydrofuran (THF) solution. The electron transfer rate is found to increase significantly with increasing coupling strength between the molecular system and the cavity. The rate process is also seen to shift from overdamped monotonic decay to under-damped oscillatory dynamics. The electron transfer rate is seen to be highly sensitive to the cavity frequency, with the emergence of a resonance cavity frequency for which the effect of coupling to the cavity is maximal. Finally, an implementation of the algorithm on the IBM Osaka quantum computer is used to demonstrate how TT-TFD-based electron transfer dynamics can be simulated accurately on NISQ computers.

Modulating molecular plasmons in naphthalene via intermolecular interactions and strong light-matter coupling

We conducted a theoretical investigation on the modulation of plasmon-like resonances in naphthalene - the so-called molecular plasmons - through intermolecular interactions and strong light-matter coupling. The configuration interaction with single excitations (CIS) approach and its quantum electrodynamics extension (QED-CIS-1) are used to describe the molecular plasmon states under these interactions. We detail the effects of changing intermolecular distances of the naphthalene dimer and incorporating the naphthalene molecule into optical cavities, both allowing for precise control of naphthalene's plasmonic responses. Our results show significant shifts of the plasmon peak in the absorption spectra of naphthalene, depending on the spatial configuration of the dimer and cavity parameters such as polarization, frequency, and coupling strength. Further investigation of the naphthalene dimer in a cavity reveals a synergistic effect on the plasmon peak when the two types of interactions are combined. This research provides insights into the plasmonic behavior of simple polyacenes like naphthalene and opens up possibilities for plasmon modulation in more complex systems.

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.

Electromagnetic Enantiomer: Chiral Nanophotonic Cavities for Inducing Chemical Asymmetry

Chiral molecules, a cornerstone of chemical sciences with applications ranging from pharmaceuticals to molecular electronics, come in mirror-image pairs called enantiomers. However, their synthesis often requires complex control of their molecular geometry. We propose a strategy called "electromagnetic enantiomers" for inducing chirality in molecules located within engineered nanocavities using light, eliminating the need for intricate molecular design. This approach works by exploiting the strong coupling between a nonchiral molecule and a chiral mode within a nanocavity. We provide evidence for this strong coupling through angular emission patterns verified by numerical simulations and with complementary evidence provided by luminescence lifetime measurements. In simpler terms, our hypothesis suggests that chiral properties can be conveyed on to a molecule with a suitable chromophore by placing it within a specially designed chiral nanocavity that is significantly larger (hundreds of nanometers) than the molecule itself. To demonstrate this concept, we showcase an application in display technology, achieving efficient emission of circularly polarized light from a nonchiral molecule. The electromagnetic enantiomer concept offers a simpler approach to chiral control, potentially opening doors for asymmetric synthesis.

i-PI 3.0: A flexible and efficient framework for advanced atomistic simulations

Atomic-scale simulations have progressed tremendously over the past decade, largely thanks to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler-Parinello, DeePMD, and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows for deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities.

Kadanoff-Baym Equations for Interacting Systems with Dissipative Lindbladian Dynamics

The extraordinary quantum properties of nonequilibrium systems governed by dissipative dynamics have become a focal point in contemporary scientific inquiry. The nonequilibrium Green's functions (NEGF) theory provides a versatile method for addressing driven nondissipative systems, utilizing the powerful diagrammatic technique to incorporate correlation effects. We here present a second-quantization approach to the dissipative NEGF theory, reformulating Keldysh ideas to accommodate Lindbladian dynamics and extending the Kadanoff-Baym equations accordingly. Generalizing diagrammatic perturbation theory for many-body Lindblad operators, the formalism enables correlated and dissipative real-time simulations for the exploration of transient and steady-state changes in the electronic, transport, and optical properties of materials.

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.

Thermal disorder prevents the suppression of ultra-fast photochemistry in the strong light-matter coupling regime

Strong coupling between molecules and confined light modes of optical cavities to form polaritons can alter photochemistry, but the origin of this effect remains largely unknown. While theoretical models suggest a suppression of photochemistry due to the formation of new polaritonic potential energy surfaces, many of these models do not account for the energetic disorder among the molecules, which is unavoidable at ambient conditions. Here, we combine simulations and experiments to show that for an ultra-fast photochemical reaction such thermal disorder prevents the modification of the potential energy surface and that suppression is due to radiative decay of the lossy cavity modes. We also show that the excitation spectrum under strong coupling is a product of the excitation spectrum of the bare molecules and the absorption spectrum of the molecule-cavity system, suggesting that polaritons can act as gateways for channeling an excitation into a molecule, which then reacts normally. Our results therefore imply that strong coupling provides a means to tune the action spectrum of a molecule, rather than to change the reaction.

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.

cclib 2.0: An updated architecture for interoperable computational chemistry

Interoperability in computational chemistry is elusive, impeded by the independent development of software packages and idiosyncratic nature of their output files. The cclib library was introduced in 2006 as an attempt to improve this situation by providing a consistent interface to the results of various quantum chemistry programs. The shared API across programs enabled by cclib has allowed users to focus on results as opposed to output and to combine data from multiple programs or develop generic downstream tools. Initial development, however, did not anticipate the rapid progress of computational capabilities, novel methods, and new programs; nor did it foresee the growing need for customizability. Here, we recount this history and present cclib 2, focused on extensibility and modularity. We also introduce recent design pivots-the formalization of cclib's intermediate data representation as a tree-based structure, a new combinator-based parser organization, and parsed chemical properties as extensible objects.

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.

Mesoscale Molecular Simulations of Fabry-Pérot Vibrational Strong Coupling

Developing theoretical frameworks for vibrational strong coupling (VSC) beyond the single-mode approximation is crucial for a comprehensive understanding of experiments with planar Fabry-Pérot cavities. Herein, a generalized cavity molecular dynamics (CavMD) scheme is developed to simulate VSC of a large ensemble of realistic molecules coupled to an arbitrary 1D or 2D photonic environment. This approach is built upon the Power-Zienau-Woolley Hamiltonian in the normal mode basis and uses a grid representation of the molecular ensembles to reduce the computational cost. When simulating the polariton dispersion relation for a homogeneous distribution of molecules in planar Fabry-Pérot cavities, our data highlight the importance of preserving the in-plane translational symmetry of the molecular distribution. In this homogeneous limit, CavMD yields the consistent polariton dispersion relation as an analytic theory, i.e., incorporating many cavity modes with varying in-plane wave vectors (k∥) produces the same spectrum as the system with a single cavity mode. Furthermore, CavMD reveals that the validity of the single-mode approximation is challenged when nonequilibrium polariton dynamics are considered, as polariton-polariton scattering occurs between modes with the nearest neighbor k∥. The procedure for numerically approaching the macroscopic limit is also demonstrated with CavMD by increasing the system size. Looking forward, our generalized CavMD approach may facilitate understanding vibrational polariton transport and condensation.

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.

Modeling the electroluminescence of atomic wires from quantum dynamics simulations

Static and time-dependent quantum-mechanical approaches have been employed in the literature to characterize the physics of light-emitting molecules and nanostructures. However, the electromagnetic emission induced by an input current has remained beyond the realm of molecular simulations. This is the challenge addressed here with the help of an equation of motion for the density matrix coupled to a photon bath based on a Redfield formulation. This equation is evolved within the framework of the driven-Liouville von Neumann approach, which incorporates open boundaries by introducing an applied bias and a circulating current. The dissipated electromagnetic power can be computed in this context from the time derivative of the energy. This scheme is applied in combination with a self-consistent tight-binding Hamiltonian to investigate the effects of bias and molecular size on the electroluminescence of metallic and semiconducting chains. For the latter, a complex interplay between bias and molecular length is observed: there is an optimal number of atoms that maximizes the emitted power at high voltages but not at low ones. This unanticipated behavior can be understood in terms of the band bending produced along the semiconducting chain, a phenomenon that is captured by the self-consistency of the method. A simple analytical model is proposed that explains the main features revealed by the simulations. The methodology, applied here at a self-consistent tight-binding level but extendable to more sophisticated Hamiltonians such as density functional tight binding and time dependent density functional theory, promises to be helpful for quantifying the power and quantum efficiency of nanoscale electroluminescent devices.

Fully Quantum Simulation of Polaritonic Vibrational Spectra of Large Cavity-Molecule System

The formation of molecular vibrational polaritons, arising from the interplay between molecular vibrations and infrared cavity modes, is a quantum phenomenon necessitating accurate quantum dynamical simulations. Here, we introduce the cavity vibrational self-consistent field/virtual state configuration interaction method, enabling quantum simulation of the vibrational spectra of many-molecule systems within the optical cavity. Focusing on the representative (H2O)21 system, we showcase this parameter-free quantum approach's ability to capture both linear and nonlinear vibrational spectral features. Our findings highlight the growing prominence of molecular couplings among OH stretches and bending excited bands with increased light-matter interaction, revealing distinctive nonlinear spectral features induced by vibrational strong coupling.

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.

Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states

A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.

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.

High-quality nanocavities through multimodal confinement of hyperbolic polaritons in hexagonal boron nitride

Compressing light into nanocavities substantially enhances light-matter interactions, which has been a major driver for nanostructured materials research. However, extreme confinement generally comes at the cost of absorption and low resonator quality factors. Here we suggest an alternative optical multimodal confinement mechanism, unlocking the potential of hyperbolic phonon polaritons in isotopically pure hexagonal boron nitride. We produce deep-subwavelength cavities and demonstrate several orders of magnitude improvement in confinement, with estimated Purcell factors exceeding 108 and quality factors in the 50-480 range, values approaching the intrinsic quality factor of hexagonal boron nitride polaritons. Intriguingly, the quality factors we obtain exceed the maximum predicted by impedance-mismatch considerations, indicating that confinement is boosted by higher-order modes. We expect that our multimodal approach to nanoscale polariton manipulation will have far-reaching implications for ultrastrong light-matter interactions, mid-infrared nonlinear optics and nanoscale sensors.

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.

Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics

The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. In this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.

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.

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.

Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction Theory

Polariton chemistry has attracted great attention as a potential route to modify chemical structure, properties, and reactivity through strong interactions among molecular electronic, vibrational, or rovibrational degrees of freedom. A rigorous theoretical treatment of molecular polaritons requires the treatment of matter and photon degrees of freedom on equal quantum mechanical footing. In the limit of molecular electronic strong or ultrastrong coupling to one or a few molecules, it is desirable to treat the molecular electronic degrees of freedom using the tools of ab initio quantum chemistry, yielding an approach we refer to as ab initio cavity quantum electrodynamics, where the photon degrees of freedom are treated at the level of cavity quantum electrodynamics. Here, we present an approach called Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction theory to provide ground- and excited-state polaritonic surfaces with a balanced description of strong correlation effects among electronic and photonic degrees of freedom. This method provides a platform for ab initio cavity quantum electrodynamics when both strong electron correlation and strong light-matter coupling are important and is an important step toward computational approaches that yield multiple polaritonic potential energy surfaces and couplings that can be leveraged for ab initio molecular dynamics simulations of polariton chemistry.

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.

Machine Learning Framework for Modeling Exciton Polaritons in Molecular Materials

A light-matter hybrid quasiparticle, called a polariton, is formed when molecules are strongly coupled to an optical cavity. Recent experiments have shown that polariton chemistry can manipulate chemical reactions. Polariton chemistry is a collective phenomenon, and its effects increase with the number of molecules in a cavity. However, simulating an ensemble of molecules in the excited state coupled to a cavity mode is theoretically and computationally challenging. Recent advances in machine learning (ML) techniques have shown promising capabilities in modeling ground-state chemical systems. This work presents a general protocol to predict excited-state properties, such as energies, transition dipoles, and nonadiabatic coupling vectors with the hierarchically interacting particle neural network. ML predictions are then applied to compute the potential energy surfaces and electronic spectra of a prototype azomethane molecule in the collective coupling scenario. These computational tools provide a much-needed framework to model and understand many molecules' emerging excited-state polariton chemistry.

Simulating the Excited-State Dynamics of Polaritons with Ab Initio Multiple Spawning

Over the past decade, there has been a growth of interest in polaritonic chemistry, where the formation of hybrid light-matter states (polaritons) can alter the course of photochemical reactions. These hybrid states are created by strong coupling between molecules and photons in resonant optical cavities and can even occur in the absence of light when the molecule is strongly coupled with the electromagnetic fluctuations of the vacuum field. We present a first-principles model to simulate nonadiabatic dynamics of such polaritonic states inside optical cavities by leveraging graphical processing units (GPUs). Our first implementation of this model is specialized for a single molecule coupled to a single-photon mode confined inside the optical cavity but with any number of excited states computed using complete active space configuration interaction (CASCI) and a Jaynes-Cummings-type Hamiltonian. Using this model, we have simulated the excited-state dynamics of a single salicylideneaniline (SA) molecule strongly coupled to a cavity photon with the ab initio multiple spawning (AIMS) method. We demonstrate how the branching ratios of the photodeactivation pathways for this molecule can be manipulated by coupling to the cavity. We also show how one can stop the photoreaction from happening inside of an optical cavity. Finally, we also investigate cavity-based control of the ordering of two excited states (one optically bright and the other optically dark) inside a cavity for a set of molecules, where the dark and bright states are close in energy.

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.

Coupled cluster cavity Born-Oppenheimer approximation for electronic strong coupling

Chemical and photochemical reactivity, as well as supramolecular organization and several other molecular properties, can be modified by strong interactions between light and matter. Theoretical studies of these phenomena require the separation of the Schrödinger equation into different degrees of freedom as in the Born-Oppenheimer approximation. In this paper, we analyze the electron-photon Hamiltonian within the cavity Born-Oppenheimer approximation (CBOA), where the electronic problem is solved for fixed nuclear positions and photonic parameters. In particular, we focus on intermolecular interactions in representative dimer complexes. The CBOA potential energy surfaces are compared with those obtained using a polaritonic approach, where the photonic and electronic degrees of freedom are treated at the same level. This allows us to assess the role of electron-photon correlation and the accuracy of CBOA.

Vacuum-Induced Symmetry Breaking of Chiral Enantiomer Formation in Chemical Reactions

A material with symmetry breaking inside can transmit the symmetry breaking to its vicinity by vacuum electromagnetic fluctuations. Here, we show that vacuum quantum fluctuations proximate to a parity-symmetry-broken material can induce a chirality-dependent spectral shift of chiral molecules, resulting in a chemical reaction process that favors producing one chirality over the other. We calculate concrete examples and evaluate the chirality production rate with experimentally realizable parameters, showing the promise of selecting chirality with symmetry-broken vacuum quantum fluctuations.

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.