We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the extended stochastic Liouville--von Neumann equation. Our approach generalizes previous work based on Caldeira-Leggett models and a partitioned initial density matrix. This provides a simple, yet exact, closed-form description for the evolution of open systems from equilibriated initial conditions. The applicability of this model and the potential for numerical implementations are also discussed.
We show that all nonrelativistic quantum processes, whether open or closed, are either unitary or probabilistic unitary, i.e., probabilistic combinations of unitary evolutions. This means that, for open quantum systems, their continuous dynamics can always be described by the Lindblad master equation with all jump operators being unitary. We call this formalism the probabilistic unitary formulation of open quantum system dynamics. This formalism is shown to be exact under all cases, and does not rely on any assumptions other than the continuity and differentiability of the density matrix. Moreover, it requires as few as $d\ensuremath{-}1$ jump operators, instead of ${d}^{2}\ensuremath{-}1$, to describe the open dynamics in the most general case, where
d is the dimension of Hilbert space of the system. Importantly, different from the conventional Lindblad master equation, this formalism is state dependent, meaning that the Hamiltonian, jump operators, and rates in general all depend on the current state of the density matrix. Hence one needs to know the explicit expression of the density matrix in order to write down the probabilistic unitary master equation explicitly. Experimentally, the formalism provides a scheme to control a quantum state to evolve along designed nonunitary quantum trajectories, and can be potentially useful in quantum computing and quantum control scenes since only unitary resources are needed for implementation.
We extend the coupled-cluster method to correlated quantum dynamics of both closed and open systems at finite temperatures using the thermofield formalism. The approach expresses the time-dependent density matrix in an exponential ansatz and describes time-evolution along the Keldysh path contour. A distinct advantage of the approach is exact trace-preservation as a function of time, ensuring conservation of probability and particle number. Furthermore, the method avoids the computation of correlated bra-states, simplifying the computational implementation. We develop the method in a thermal quasiparticle representation, which allows seamless connection to the projection method and diagrammatic techniques of the traditional coupled-cluster formalism. For comparison, we also apply the thermofield framework to the density-matrix renormalization-group method to obtain reference results for closed and open systems at finite temperature. We test the singles and doubles approximation to the density-matrix coupled-cluster method on the correlated electronic dynamics of the single-impurity Anderson model, demonstrating that the new method successfully captures the correlated dynamics of both closed systems at finite temperature and driven-dissipative open systems. This encouraging performance motivates future applications to nonequilibrium quantum many-body dynamics in realistic systems.