Molecular van der Waals Fluids in Cavity Quantum Electrodynamics

Diagrams in Polaritonic Coupled Cluster Theory

We present a diagrammatic notation to derive the quantum-electrodynamic coupled cluster (QED-CC) equations needed for the description for polaritonic ground and excited states. Our presented notation is a generalization of the existing diagrammatic notation of standard electronic coupled-cluster theory. It is used to derive the QED-CC and QED-EOM-CC equations for the QED-CCSD-1-SD and QED-CCSD-12-SD truncation schemes. Furthermore, we present the diagrams for the CC Λ-equations and reduced density matrices which are needed for the calculation of molecular properties.

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.

The orientation dependence of cavity-modified chemistry

Recent theoretical studies have explored how ultra-strong light-matter coupling can be used as a handle to control chemical transformations. Ab initio cavity quantum electrodynamics calculations demonstrate that large changes to reaction energies or barrier heights can be realized by coupling electronic degrees of freedom to vacuum fluctuations associated with an optical cavity mode, provided that large enough coupling strengths can be achieved. In many cases, the cavity effects display a pronounced orientational dependence. Here, we highlight the critical role that geometry relaxation can play in such studies. As an example, we consider a recent work [Pavošević et al., Nat. Commun. 14, 2766 (2023)] that explored the influence of an optical cavity on Diels-Alder cycloaddition reactions and reported large changes to reaction enthalpies and barrier heights, as well as the observation that changes in orientation can inhibit the reaction or select for one reaction product or another. Those calculations used fixed molecular geometries optimized in the absence of the cavity and fixed relative orientations of the molecules and the cavity mode polarization axis. Here, we show that when given a chance to relax in the presence of the cavity, the molecular species reorient in a way that eliminates the orientational dependence. Moreover, in this case, we find that qualitatively different conclusions regarding the impact of the cavity on the thermodynamics of the reaction can be drawn from calculations that consider relaxed vs unrelaxed molecular structures.

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.

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.

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.