Search for: All records

We derive a rigorous nuclear gradient for a molecule-cavity hybrid system using the quantum electrodynamics Hamiltonian. We treat the electronic–photonic degrees of freedom (DOFs) as the quantum subsystem and the nuclei as the classical subsystem. Using the adiabatic basis for the electronic DOF and the Fock basis for the photonic DOF and requiring the total energy conservation of this mixed quantum–classical (MQC) system, we derived the rigorous nuclear gradient for the molecule–cavity hybrid system, which is naturally connected to the approximate gradient under the Jaynes–Cummings approximation. The nuclear gradient expression can be readily used in any MQC simulations and will allow one to perform the non-adiabatic on-the-fly simulation of polariton quantum dynamics. The theoretical developments in this work could significantly benefit the polariton quantum dynamics community with a rigorous nuclear gradient of the molecule–cavity hybrid system and have a broad impact on the future non-adiabatic simulations of polariton quantum dynamics. 

We present the rigorous theoretical framework of the generalized spin mapping representation for non-adiabatic dynamics. Our work is based upon a new mapping formalism recently introduced by Runeson and Richardson [J. Chem. Phys. 152, 084110 (2020)], which uses the generators of the [Formula: see text] Lie algebra to represent N discrete electronic states, thus preserving the size of the original Hilbert space. Following this interesting idea, the Stratonovich–Weyl transform is used to map an operator in the Hilbert space to a continuous function on the SU( N) Lie group, i.e., a smooth manifold which is a phase space of continuous variables. We further use the Wigner representation to describe the nuclear degrees of freedom and derive an exact expression of the time-correlation function as well as the exact quantum Liouvillian for the non-adiabatic system. Making the linearization approximation, this exact Liouvillian is reduced to the Liouvillian of several recently proposed methods, and the performance of this linearized method is tested using non-adiabatic models. We envision that the theoretical work presented here provides a rigorous and unified framework to formally derive non-adiabatic quantum dynamics approaches with continuous variables and connects the previous methods in a clear and concise manner. 

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

This work provides the fundamental theoretical framework for few-mode cavity quantum electrodynamics by resolving the gauge ambiguities between the Coulomb gauge and the dipole gauge Hamiltonians under the photonic mode truncation. We first propose a general framework to resolve ambiguities for an arbitrary truncation in a given gauge. Then, we specifically consider the case of mode truncation, deriving gauge invariant expressions for both the Coulomb and dipole gauge Hamiltonians that naturally reduce to the commonly used single-mode Hamiltonians when considering a single-mode truncation. We finally provide the analytical and numerical results of both atomic and molecular model systems coupled to the cavity to demonstrate the validity of our theory.

 

We investigate the Polariton induced conical intersection (PICI) created from coupling a diatomic molecule with the quantized photon mode inside an optical cavity, and the corresponding Berry Phase effects. We use the rigorous Pauli–Fierz Hamiltonian to describe the quantum light-matter interactions between a LiF molecule and the cavity, and use the exact quantum propagation to investigate the polariton quantum dynamics. The molecular rotations relative to the cavity polarization direction play a role as the tuning mode of the PICI, resulting in an effective CI even within a diatomic molecule. To clearly demonstrate the dynamical effects of the Berry phase, we construct two additional models that have the same Born–Oppenheimer surface, but the effects of the geometric phase are removed. We find that when the initial wavefunction is placed in the lower polaritonic surface, the Berry phase causes a π phase-shift in the wavefunction after the encirclement around the CI, indicated from the nuclear probability distribution. On the other hand, when the initial wavefunction is placed in the upper polaritonic surface, the geometric phase significantly influences the couplings between polaritonic states and therefore, the population dynamics between them. These BP effects are further demonstrated through the photo-fragment angular distribution. PICI created from the quantized radiation field has the promise to open up new possibilities to modulate photochemical reactivities.