Stochastic simulation of anharmonic dissipation. I. Linear response regime

Hierarchical Equations-of-Motion Method for Momentum System-Bath Coupling

For a broad class of quantum models of practical interest, we demonstrate that the Hamiltonian of the system nonlinearly coupled to a harmonic bath through the system and bath coordinates can be equivalently mapped into the Hamiltonian of the system bilinearly coupled to the bath through the system and bath momenta. We show that the Hamiltonian with bilinear system-bath momentum coupling can be treated by the hierarchical equations-of-motion (HEOM) method and present the corresponding proof-of-principle simulations. The developed methodology creates the opportunity to scrutinize a new family of nonlinear quantum systems by the numerically accurate HEOM method.

Stochastic simulation of anharmonic dissipation. II. Harmonic bath potentials with quadratic couplings

The workhorse simulating the dissipative dynamics is mainly based on the harmonic bath potentials together with linear system-bath couplings, but a realistic bath always assumes anharmonicity. In this work, we extend the linear dissipation model to include quadratic couplings and suggest a stochastic simulation scheme for the anharmonic dissipation. We show that the non-Gaussian noises induced by the anharmonic bath can be rigorously constructed, and the resulting stochastic Liouville equation has the same form as that for the linear dissipation model. As a preliminary application, we use this stochastic method to investigate the vibration-induced symmetry breaking in two-level electronic systems and find that the characteristic function of the non-Gaussian noises determines the absorption and fluorescence spectra.