Search for: All records

Abstract Recent experiments demonstrate polaritons under the vibrational strong coupling (VSC) regime can modify chemical reactivity. Here, we present a complete theory of VSC-modified rate constants when coupling a single molecule to an optical cavity, where the role of photonic mode lifetime is understood. The analytic expression exhibits a sharp resonance behavior, where the maximum rate constant is reached when the cavity frequency matches the vibration frequency. The theory explains why VSC rate constant modification closely resembles the optical spectra of the vibration outside the cavity. Further, we discussed the temperature dependence of the VSC-modified rate constants. The analytic theory agrees well with the numerically exact hierarchical equations of motion (HEOM) simulations for all explored regimes. Finally, we discussed the resonance condition at the normal incidence when considering in-plane momentum inside a Fabry-Pérot cavity. 

Abstract We present a theory that explains the resonance effect of the vibrational strong coupling (VSC) modified reaction rate constant at the normal incidence of a Fabry–Pérot (FP) cavity. This analytic theory is based on a mechanistic hypothesis that cavity modes promote the transition from the ground state to the vibrational excited state of the reactant, which is the rate-limiting step of the reaction. This mechanism for a single molecule coupled to a single-mode cavity has been confirmed by numerically exact simulations in our recent work in [J. Chem. Phys. 159, 084104 (2023)]. Using Fermi’s golden rule (FGR), we formulate this rate constant for many molecules coupled to many cavity modes inside a FP microcavity. The theory provides a possible explanation for the resonance condition of the observed VSC effect and a plausible explanation of why only at the normal incident angle there is the resonance effect, whereas, for an oblique incidence, there is no apparent VSC effect for the rate constant even though both cases generate Rabi splitting and forming polariton states. On the other hand, the current theory cannot explain the collective effect when a large number of molecules are collectively coupled to the cavity, and future work is required to build a complete microscopic theory to explain all observed phenomena in VSC. 

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. 

When matter is strongly coupled to an optical cavity, new hybrid light–matter states are formed, the so-called polariton states. These polaritons can qualitatively change the physical properties of the matter coupled to the cavity by completely altering its energy eigenspectrum. Fueled by experimental innovations in recent years, much progress has been made in simulating the intrinsic quantum behavior of these hybrid states. At the heart of each simulation is the choice of Hamiltonian to represent the total light–matter hybrid system. Even at this fundamental level, there has been significant progress in developing new gauges and representations for this Hamiltonian, whether exact or under approximations. As such, this review aims to discuss several different forms of Hamiltonians for the researcher trying to enter this field by clearly and concisely deriving each different representation from the fundamental Minimal Coupling Hamiltonian. In addition, this review provides commentary on the optimal usage and extent of approximations for each individual representation to assist the reader in choosing the appropriate Hamiltonian for their work. 

We perform on-the-fly non-adiabatic molecular dynamics simulations using the recently developed spin-mapping formalism. Two quantum dynamics approaches based on this mapping formalism, (i) the fully linearized Spin-LSC and (ii) the partially linearized Spin-PLDM, are explored using the quasi-diabatic propagation scheme. We have performed dynamics simulations in four ab initio molecular models for which benchmark ab initio multiple spawning (AIMS) data have been published. We find that the spin-LSC and the previously reported symmetric quasi-classical (SQC) approaches provide nearly equivalent population dynamics. While we expected the more involved spin-PLDM method to provide superior accuracy compared to the other mapping-based approaches, SQC and spin-LSC, we found that it performed with equivalent accuracy compared to the AIMS benchmark results. We further explore the underpinnings of the spin-PLDM correlation function by decomposing its N2 density matrix-focused initial conditions, where N is the number of states in the quantum subsystem. Finally, we found an approximate form of the spin-PLDM correlation function, which simplifies the simulation and reduces the computational costs from N2 to N. 

In our previous work [Mondal et al., J. Chem. Phys. 162, 014114 (2025)], we developed several efficient computational approaches to simulate exciton–polariton dynamics described by the Holstein–Tavis–Cummings (HTC) Hamiltonian under the collective coupling regime. Here, we incorporated these strategies into the previously developed Lindblad-partially linearized density matrix (L-PLDM) approach for simulating 2D electronic spectroscopy (2DES) of exciton–polariton under the collective coupling regime. In particular, we apply the efficient quantum dynamics propagation scheme developed in Paper I to both the forward and the backward propagations in the PLDM and develop an efficient importance sampling scheme and graphics processing unit vectorization scheme that allow us to reduce the computational costs from O(K2)O(T3) to O(K)O(T0) for the 2DES simulation, where K is the number of states and T is the number of time steps of propagation. We further simulated the 2DES for an HTC Hamiltonian under the collective coupling regime and analyzed the signal from both rephasing and non-rephasing contributions of the ground state bleaching, excited state emission, and stimulated emission pathways. 

We outline two general theoretical techniques to simulate polariton quantum dynamics and optical spectra under the collective coupling regimes described by a Holstein–Tavis–Cummings (HTC) model Hamiltonian. The first one takes advantage of sparsity of the HTC Hamiltonian, which allows one to reduce the cost of acting polariton Hamiltonian onto a state vector to the linear order of the number of states, instead of the quadratic order. The second one is applying the well-known Chebyshev series expansion approach for quantum dynamics propagation and to simulate the polariton dynamics in the HTC system; this approach allows us to use a much larger time step for propagation and only requires a few recursive operations of the polariton Hamiltonian acting on state vectors. These two theoretical approaches are general and can be applied to any trajectory-based non-adiabatic quantum dynamics methods. We apply these two techniques with our previously developed Lindblad-partially linearized density matrix approach to simulate the linear absorption spectra of the HTC model system, with both inhomogeneous site energy disorders and dipolar orientational disorders. Our numerical results agree well with the previous analytic and numerical work. 

We investigate the quantum dynamics of a spin coupling to a bath of independent spins via the dissipaton equation of motion (DEOM) approach. The bath, characterized by a continuous spectral density function, is composed of spins that are independent level systems described by the su(2) Lie algebra, representing an environment with a large magnitude of anharmonicity. Based on the previous work by Suarez and Silbey [J. Chem. Phys. 95, 9115 (1991)] and by Makri [J. Chem. Phys. 111, 6164 (1999)] that the spin bath can be mapped to a Gaussian environment under its linear response limit, we use the time-domain Prony fitting decomposition scheme to the bare–bath time correlation function (TCF) given by the bosonic fluctuation–dissipation theorem to generate the exponential decay basis (or pseudo modes) for DEOM construction. The accuracy and efficiency of this strategy have been explored by a variety of numerical results. We envision that this work provides new insights into extending the hierarchical equations of motion and DEOM approach to certain types of anharmonic environments with arbitrary TCF or spectral density.