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

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.

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.

Simulations of photoinduced processes with the exact factorization: state of the art and perspectives

This perspective offers an overview of the applications of the exact factorization of the electron-nuclear wavefunction to the domain of theoretical photochemistry, where the aim is to gain insights into the ultrafast dynamics of molecular systems via simulations of their excited-state dynamics beyond the Born-Oppenheimer approximation. The exact factorization offers an alternative viewpoint to the Born-Huang representation for the interpretation of dynamical processes involving the electronic ground and excited states as well as their coupling through the nuclear motion. Therefore, the formalism has been used to derive algorithms for quantum molecular-dynamics simulations where the nuclear motion is treated using trajectories and the electrons are treated quantum mechanically. These algorithms have the characteristic features of being based on coupled and on auxiliary trajectories, and have shown excellent performance in describing a variety of excited-state processes, as this perspective illustrates. We conclude with a discussion on the authors' point of view on the future of the exact factorization.

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.

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

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

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.

Direct Observation of Polaritonic Chemistry by Nuclear Magnetic Resonance Spectroscopy

Polaritonic chemistry is emerging as a powerful approach to modifying the properties and reactivity of molecules and materials. However, probing how the electronics and dynamics of molecular systems change under strong coupling has been challenging due to the narrow range of spectroscopic techniques that can be applied in situ. Here we develop microfluidic optical cavities for vibrational strong coupling (VSC) that are compatible with nuclear magnetic resonance (NMR) spectroscopy using standard liquid NMR tubes. VSC is shown to influence the equilibrium between two conformations of a molecular balance sensitive to London dispersion forces, revealing an apparent change in the equilibrium constant under VSC. In all compounds studied, VSC does not induce detectable changes in chemical shifts, J-couplings, or spin-lattice relaxation times. This unexpected finding indicates that VSC does not substantially affect molecular electron density distributions, and in turn has profound implications for the possible mechanisms at play in polaritonic chemistry under VSC and suggests that the emergence of collective behavior is critical.

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

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

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.

Semiclassical Truncated-Wigner-Approximation Theory of Molecular Vibration-Polariton Dynamics in Optical Cavities

It has been experimentally demonstrated that molecular-vibration polaritons formed by strong coupling of a molecular vibration to an infrared cavity mode can significantly modify the physical properties and chemical reactivities of various molecular systems. However, a complete theoretical understanding of the underlying mechanisms of the modifications remains elusive due to the complexity of the hybrid system, especially the collective nature of polaritonic states in systems containing many molecules. We develop here the semiclassical theory of molecular vibration-polariton dynamics based on the truncated Wigner approximation (TWA) that is tractable in large molecular systems and simultaneously captures the quantum character of photons in the optical cavity. The theory is then applied to investigate the nuclear quantum dynamics of a system of identical diatomic molecules having the ground-state Morse potential and being strongly coupled to an infrared cavity mode in the ultrastrong coupling regime. The validity of TWA is examined by comparing it with the full quantum dynamics of a single-molecule system for two different initial states in the dipole and Coulomb gauges. For the initial tensor-product ground state in the dipole gauge, which corresponds to a light-matter entangled state in the Coulomb gauge, the collective and resonance effects of molecular vibration-polariton formation on the nuclear dynamics are observed in a system of many molecules.

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.

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-Mediated Molecular Entanglement and Generation of Non-classical States of Light

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

Ab Initio 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.

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.