We develop an exact quantum thermostat that generates the quantum Boltzmann distribution in path integral phase space and so correctly captures the imaginary position–momentum correlations. This extends beyond the capability of existing quantum thermostats such as the PILE thermostat, which can only achieve the marginal of this distribution in position, with the momentum acting purely as a sampling tool and satisfying classical Maxwell–Boltzmann statistics. The dynamics correspond to the treatment of the Caldeira–Leggett system–bath model with Matsubara dynamics, a theory for approximating nuclear quantum dynamics by propagating smooth imaginary-time Feynman paths in real time. A Matsubara generalised Langevin equation that obeys certain quantum fluctuation–dissipation relations is mapped onto a Markovian stochastic process by introducing auxiliary variables; this allows the construction of a corresponding Fokker–Planck equation from which the stationary distribution in the extended variable space is determined. We show that this corresponds to the joint quantum Boltzmann distribution in position and momentum, a finding that holds for any bath spectral density and for anharmonic systems. The sole approximation in Matsubara dynamics is the neglect of real-time quantum coherence; consequently, provided there is sufficient quantum decoherence, the dynamics of the exact quantum thermostat also well approximates the exact quantum dynamics of the system–bath. In the limit of infinite system–bath coupling, the level of real-time coherence becomes negligible and so the Matsubara dynamics of the thermostat tends to the exact dynamics. The quantum thermostat can therefore be used to study this limit exactly and, in so doing, we demonstrate the existence of a quantum mechanical overdamped regime analogous to its classical counterpart that is characterised by a Matsubara Smoluchowski equation. The drawback of the quantum thermostat is that one of its associated quantum fluctuation–dissipation relations is imaginary. This results in complex noise and dynamical instabilities that often preclude the numerical implementation of the thermostat. Instead, it primarily serves as a theoretical framework that can be used as a starting point for introducing approximations that yield more practical methods. One such approximation that results in real noise is considered here and leads to a practical approach for calculating linear time-correlation functions of bound open quantum systems. These exhibit good agreement with exact quantum mechanical calculations, particularly at strong system–bath coupling, where the real-time quantum coherence is increasingly suppressed.
We develop an exact quantum thermostat that generates the quantum Boltzmann distribution in path integral phase space and so correctly captures the imaginary position–momentum correlations. This extends beyond the capability of existing quantum thermostats such as the PILE thermostat, which can only achieve the marginal of this distribution in position, with the momentum acting purely as a sampling tool and satisfying classical Maxwell–Boltzmann statistics. The dynamics correspond to the treatment of the Caldeira–Leggett system–bath model with Matsubara dynamics, a theory for approximating nuclear quantum dynamics by propagating smooth imaginary-time Feynman paths in real time. A Matsubara generalised Langevin equation that obeys certain quantum fluctuation–dissipation relations is mapped onto a Markovian stochastic process by introducing auxiliary variables; this allows the construction of a corresponding Fokker–Planck equation from which the stationary distribution in the extended variable space is determined. We show that this corresponds to the joint quantum Boltzmann distribution in position and momentum, a finding that holds for any bath spectral density and for anharmonic systems. The sole approximation in Matsubara dynamics is the neglect of real-time quantum coherence; consequently, provided there is sufficient quantum decoherence, the dynamics of the exact quantum thermostat also well approximates the exact quantum dynamics of the system–bath. In the limit of infinite system–bath coupling, the level of real-time coherence becomes negligible and so the Matsubara dynamics of the thermostat tends to the exact dynamics. The quantum thermostat can therefore be used to study this limit exactly and, in so doing, we demonstrate the existence of a quantum mechanical overdamped regime analogous to its classical counterpart that is characterised by a Matsubara Smoluchowski equation. The drawback of the quantum thermostat is that one of its associated quantum fluctuation–dissipation relations is imaginary. This results in complex noise and dynamical instabilities that often preclude the numerical implementation of the thermostat. Instead, it primarily serves as a theoretical framework that can be used as a starting point for introducing approximations that yield more practical methods. One such approximation that results in real noise is considered here and leads to a practical approach for calculating linear time-correlation functions of bound open quantum systems. These exhibit good agreement with exact quantum mechanical calculations, particularly at strong system–bath coupling, where the real-time quantum coherence is increasingly suppressed.
We report the first numerical calculations in which converged Matsubara dynamics is compared directly with exact quantum dynamics with no artificial damping of the time-correlation functions (TCFs). The system treated is a Morse oscillator coupled to a harmonic bath. We show that, when the system-bath coupling is sufficiently strong, the Matsubara calculations can be converged by explicitly including up to M = 200 Matsubara modes, with the remaining modes included as a harmonic "tail" correction. The resulting Matsubara TCFs are in near-perfect agreement with the exact quantum TCFs, for non-linear as well as linear operators, at a temperature at which the TCFs are dominated by quantum thermal fluctuations. These results provide compelling evidence that incoherent classical dynamics can arise in the condensed phase at temperatures at which the statistics are dominated by quantum (Boltzmann) effects, as a result of smoothing of imaginary-time Feynman paths. The techniques developed here may also lead to efficient methods for benchmarking system-bath dynamics in the overdamped regime.
Beyond the weak-coupling limit, open quantum systems equilibrate to a highly entangled thermal state. For continuous-variable systems, this state can be written explicitly as an imaginary-time phase-space path integral, in which the positions are directly entangled with the bath, and the momenta are correlated with the positions through a phase term. Here, we ask to what extent this state can be reached by propagating stochastic classical trajectories in path-integral phase space. Surprisingly, we find that the trajectories equilibrate to the exact quantum equilibrium state, recovering the purely imaginary momentum-position correlation in the phase term. The trajectories are generated using a recently derived Matsubara generalized Langevin equation, which produces the imaginary correlations by evolving the stochastic variables into the complex plane. This makes the dynamics numerically unstable, but we are nonetheless able to demonstrate the equilibration of a quartic oscillator coupled to a white-noise bath. These unexpected findings could lead to new approximate methodologies for simulating continuous-variable open quantum systems.
We consider a quantum mechanical system represented in phase space (referred to hereafter as "Wigner space"), coupled to a harmonic oscillator bath. We derive quantum hierarchal Fokker-Planck (QHFP) equations not only in real time but also in imaginary time, which represents an inverse temperature. This is an extension of a previous work, in which we studied a spin-boson system, to a Brownian system. It is shown that the QHFP in real time obtained from a correlated thermal equilibrium state of the total system possesses the same form as those obtained from a factorized initial state. A modified terminator for the hierarchal equations of motion is introduced to treat the non-Markovian case more efficiently. Using the imaginary-time QHFP, numerous thermodynamic quantities, including the free energy, entropy, internal energy, heat capacity, and susceptibility, can be evaluated for any potential. These equations allow us to treat non-Markovian, non-perturbative system-bath interactions at finite temperature. Through numerical integration of the real-time QHFP for a harmonic system, we obtain the equilibrium distributions, the auto-correlation function, and the first- and second-order response functions. These results are compared with analytically exact results for the same quantities. This provides a critical test of the formalism for a non-factorized thermal state and elucidates the roles of fluctuation, dissipation, non-Markovian effects, and system-bath coherence. Employing numerical solutions of the imaginary-time QHFP, we demonstrate the capability of this method to obtain thermodynamic quantities for any potential surface. It is shown that both types of QHFP equations can produce numerical results of any desired accuracy. The FORTRAN source codes that we developed, which allow for the treatment of Wigner space dynamics with any potential form (TanimuranFP15 and ImTanimuranFP15), are provided as the supplementary material.
We establish that the exact quantum dynamics of a Brownian particle in the Caldeira-Leggett model, with at most quadratic external potential, can be mapped, at any temperature, onto a classical, non-Markovian stochastic process in phase space. Starting from a correlated thermal equilibrium state between the particle and bath, we demonstrate that this correspondence is exact for quadratic potentials under arbitrary quantum state preparations of the particle itself. Our approach allows to consider arbitrary initial quantum states - including highly non-classical superpositions - which are incorporated via their Wigner functions, which serve as statistical weights for trajectory ensembles. Furthermore, the formalism naturally accommodates external manipulations and measurements modeled by preparation functions acting at arbitrary times, enabling the simulation of complex driven-dissipative quantum protocols. For more general, smooth potentials, we identify a natural small parameter: the density matrix becomes strongly quasidiagonal in the coordinate representation, with its off-diagonal width shrinking as the bath's spectral cutoff increases, suggesting a controlled parameter for a possible approximation.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.
Conventional derivations of the Wigner phase-space Fokker–Planck equation for open quantum systems often rely on heuristic coarse-graining of system-bath interactions. We present quantum dissipation in a minimally extended state space (QD-MESS), a unified framework that systematically models Gaussian bath dynamics in Wigner phase space with controllable numerical accuracy, yields a transparent quantum-classical picture for arbitrary noise spectra, and recovers the quantum-optical, Caldeira–Leggett, and quantum Smoluchowski limits.