We build an exact framework to evaluate heat, energy, and particle transport between Gaussian reservoirs mediated by a quadratic quantum system. By combining full counting statistics with newly developed non-Markovian master equation approaches, we introduce an effective master equation whose solution can be used to generate arbitrary moments of the heat statistics for any number of reservoirs. This theory applies equally to fermionic and bosonic systems, holds at arbitrarily strong coupling, and resolves out-of-equilibrium transient dynamics determined by the system's initial state. In the steady-state, weak-coupling limit, we recover results analogous to those of the well-known Landauer-Büttiker formalism. We conclude our discussion by demonstrating an application of the method to a prototypical fermionic system. Our results uncover a regime of transient negative heat conductance contingent upon the initial system preparation, providing a clear signature of non-trivial out-of-equilibrium dynamics.
Abstract A system of <?CDATA
N≫1?> N ≫ 1 interacting spinless quantum particles, described by a statistical operator F(t), is considered. A time-dependent projection operator formalism for a family of projectors, which select a statistical operator <?CDATA
FS(t)?> F S ( t ) for a group of S < N relevant particles by integration of the variables of the irrelevant N − S ‘environment’ particles, is presented. This formalism results in a nonlinear version of the inhomogeneous Nakajima–Zwanzig generalized master equation (GME) for the relevant part of the statistical operator F(t), which contains an undesirable initial condition term related to the irrelevant part of the statistical operator depending on N variables. By introducing an additional operator identity, the obtained nonlinear inhomogeneous Nakajima–Zwanzig equation is exactly converted into a homogeneous nonlinear equation, accounting for initial correlations in the kernel governing its evolution. Both of these equations are equivalent to the corresponding evolution equations for <?CDATA
FS(t)?> F S ( t ) and take into account the dynamics of the environment particles. The equations for <?CDATA
FS(t)?> F S ( t ) are specialized in the linear approximation for the particle density n and analyzed for a one-particle statistical operator <?CDATA
F1(t)?> F 1 ( t ) . It turns out that, in addition to the correlations caused by interparticle interaction, the irrelevant initial condition term in the inhomogeneous equation for <?CDATA
F1(t)?> F 1 ( t ) , defined by the two-particle correlations, also contains quantum correlations conditioned by the quantum particle statistics. The latter do not vanish with time and cannot be disregarded in the conventional way by applying e.g. Bogoliubov’s principle of weakening of initial correlations. In order to take initial quantum and other correlations (commonly discarded in an unconvincing way) into consideration, the obtained nonlinear homogeneous GME has been applied. Exact in the first approximation in n, we obtain a new homogeneous equation for <?CDATA
F1(t)?> F 1 ( t ) , which accounts for correlations and is valid for any time moment t. It is shown how this reversible-in-time equation changes at different timescales and converts into the irreversible quantum Boltzmann equation on the macroscopic timescale without the use of the conventional approximations.
We develop an exact framework for describing the non-Markovian dynamics of an open quantum system interacting with an environment modeled by a generalized spectral density function.The approach relies on mapping the initial system onto an auxiliary configuration, comprising the original open system coupled to a small number of discrete modes, which in turn are each coupled to an independent Markovian reservoir.Based on the connection between the discrete modes and the poles of the spectral density function, we show how expanding the system using the discrete modes allows for the full inclusion non-Markovian effects within an enlarged open system whose dynamics is governed by an exact Lindblad master equation.Initially we apply this result to obtain a generalization of the pseudomode method [B.M. Garraway, Phys. Rev. A. 55, 2290 (1997)] in cases where the spectral density function has a Lorentzian structure.For many other types of spectral density function, we extend our proof to show that an open system dynamics may be modeled physically using discrete modes which admit a non-Hermitian coupling to the system and for such cases determine the equivalent master equation to no longer be of Lindblad form.For applications involving two discrete modes, we demonstrate how to convert between pathological and Lindblad forms of the master equation using the techniques of the pseudomode method.
The Markovian dynamics of open quantum many-body systems are typically governed by the Lindblad master equation, yet obtaining the reduced density matrix for bosonic and fermionic systems remains a formidable numerical challenge due to the exponential growth of the Hilbert space. Here, we introduce a computationally efficient framework for constructing the reduced density matrix by modeling the environment as a copy of the primary system with a monotonically decaying coupling that enforces unidirectional energy flow. Our method directly yields exact solutions that rigorously satisfy the Lindblad master equation for both bosonic and fermionic cases, bypassing the need for costly Liouvillian diagonalization. We validate our approach through applications to paradigmatic models, demonstrating accurate reproduction of dissipative dynamics across a broad parameter regime. This work provides a straightforward and powerful tool for simulating Markovian open-system evolution, with immediate applicability to quantum transport and control problems in many-body physics.
In this paper we generalize the quantum Brownian motion (QBM) to include the momentum-dependent system-environment couplings. The resulting Hamiltonian in the particle number representation is given by ${H}_{t}ot=\ensuremath{\hbar}{\ensuremath{\omega}}_{\mathrm{S}}{a}^{\ifmmode\dagger\else\textdagger\fi{}}a+{\ensuremath{\sum}}_{k}\ensuremath{\hbar}{\ensuremath{\omega}}_{k}{b}_{k}^{\ifmmode\dagger\else\textdagger\fi{}}{b}_{k}+{\ensuremath{\sum}}_{k}\ensuremath{\hbar}({V}_{k}{a}^{\ifmmode\dagger\else\textdagger\fi{}}{b}_{k}+{V}_{k}^{*}{b}_{k}^{\ifmmode\dagger\else\textdagger\fi{}}a)+{\ensuremath{\sum}}_{k}\ensuremath{\hbar}({W}_{k}{a}^{\ifmmode\dagger\else\textdagger\fi{}}{b}_{k}^{\ifmmode\dagger\else\textdagger\fi{}}+{W}_{k}^{*}{b}_{k}a)$. The conventional QBM model corresponds to the special case
Wk=Vk. The generalized QBM can capture more physical applications because the single-particle transition and the two-particle pair production between the system and the environment represent two very different physical processes that usually do not have the same coupling strengths. We discuss the physical realizations of the generalized QBM in different physical systems, and derive its exact master equation for both the initial decoupled states and initial correlated states. The Hu-Paz-Zhang master equation of the conventional QBM model is reproduced as a special case. We find that the renormalized system Hamiltonian after tracing out all the environmental states induced naturally a momentum-dependent potential, which shows the consistency of including the momentum-dependent coupling in the QBM Hamiltonian. In the Hu-Paz-Zhang master equation, such a renormalized potential is misplaced so that the correct renormalization Hamiltonian has not been found. With the exact master equation for both the initial decoupled and initial correlated states, the issues about the initial jolt are also reexamined. We find that the so-called ``initial jolt,'' which has been thought to be an artificial effect due to the use of the initial decoupled system-environment states, has nothing to do with the initial decoupled state. The exact master equation for the generalized QBM also has potential applications to photonics quantum computing.
The non-Markovian master equation for open quantum systems is obtained by generalization of the ordinary Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations for the reduced density matrices of the bath
ϱB(t) and the system
ϱS(t) are written down. A formal solution of the equation for
ϱB(t), having been inserted in the equation for the reduced density matrix of the system, in the 2-nd approximation in interaction yields a very specific extra term in the generalized master equation. This term, being nonlinear in
ϱS(t), is related to the intrinsic bath dynamics and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we applied the generalized ZN projection scheme to a simple dephasing model. It is shown that consideration of the lowest order in interaction is insufficient to describe time evolution of the system coherence adequately. We explain this fact by analyzing the exact and approximate forms of
ϱB(t) and give some hints how to take the dynamic correlations (which originate from the spin-bath coupling) into account.