Local quantum master equations provide a simple description of interacting subsystems coupled to different reservoirs. They have been widely used to study nonequilibrium critical phenomena in open quantum systems. We here investigate the validity of such a local approach by analyzing a paradigmatic system made of two harmonic oscillators each in contact with a heat bath. We evaluate the steady-state mean occupation number for varying temperature differences and find that local master equations generally fail to reproduce the results of an exact quantum-Langevin-equation description. We relate this property to the inability of the local scheme to properly characterize intersystem correlations, which we quantify with the quantum mutual information.
Local quantum master equations provide a simple description of interacting subsystems coupled to different reservoirs. They have been widely used to study nonequilibrium critical phenomena in open quantum systems. We here investigate the validity of such a local approach by analyzing a paradigmatic system made of two harmonic oscillators each in contact with a heat bath. We evaluate the steady-state mean occupation number for varying temperature differences and find that local master equations generally fail to reproduce the results of an exact quantum-Langevin-equation description. We relate this property to the inability of the local scheme to properly characterize intersystem correlations, which we quantify with the quantum mutual information.
We investigate how to model Markovian evolution of coupled harmonic oscillators, each of them interacting with a local environment. When the coupling between the oscillators is weak, dissipation may be modeled using local Lindblad terms for each of the oscillators in the master equation, as is commonly done. When the coupling between oscillators is strong, this model may become invalid. We derive a master equation for two coupled harmonic oscillators that are subject to individual heat baths modeled by a collection of harmonic oscillators and show that this master equation in general contains nonlocal Lindblad terms. We compare the resulting time evolution with that obtained for dissipation through local Lindblad terms for each individual oscillator and show that the evolution is different in the two cases. In particular, the two descriptions give different predictions for the steady state and for the entanglement between strongly coupled oscillators. This shows that when describing strongly coupled harmonic oscillators, one must take great care in how dissipation is modeled and that a description using local Lindblad terms may fail. This may be particularly relevant when attempting to generate entangled states of strongly coupled quantum systems.
When deriving a master equation for a multipartite weakly-interacting open quantum systems, dissipation is often addressed locally on each component, i.e. ignoring the coherent couplings, which are later added ‘by hand’. Although simple, the resulting local master equation (LME) is known to be thermodynamically inconsistent. Otherwise, one may always obtain a consistent global master equation (GME) by working on the energy basis of the full interacting Hamiltonian. Here, we consider a two-node ‘quantum wire’ connected to two heat baths. The stationary solution of the LME and GME are obtained and benchmarked against the exact result. Importantly, in our model, the validity of the GME is constrained by the underlying secular approximation. Whenever this breaks down (for resonant weakly-coupled nodes), we observe that the LME, in spite of being thermodynamically flawed: (a) predicts the correct steady state, (b) yields with the exact asymptotic heat currents, and (c) reliably reflects the correlations between the nodes. In contrast, the GME fails at all three tasks. Nonetheless, as the inter-node coupling grows, the LME breaks down whilst the GME becomes correct. Hence, the global and local approach may be viewed as complementary tools, best suited to different parameter regimes.
Identifying which master equation is preferable for the description of a multipartite open quantum system is not trivial and has led in recent years to the local vs global debate in the context of Markovian dissipation. We treat here a paradigmatic scenario in which the system is composed of two interacting harmonic oscillators A and B, with only A interacting with a thermal bath---collection of other harmonic oscillators---and we study the equilibration process of the system initially in the ground state with the bath finite temperature. We show that the completely positive version of the Redfield equation obtained using coarse-grain and an appropriate time-dependent convex mixture of the local and global solutions give rise to the most accurate approximations of the whole exact system dynamics, i.e., both at short and at long timescales, outperforming the local and global approaches.
In the following, we study the dissipative time evolution of a quantum chain consisting of three coupled harmonic oscillators, the first and third of which weakly interact quadratically with two independent thermal baths in equilibrium at different temperatures. Due to the quadratic form of the total Hamiltonian, the uni- tary dynamics of the compound system is formally analytically solvable and defines a one-parameter group of Gaussian maps which enables us to solve the exact dynamics of the chain numerically. Following the Gorini-Kossakowski-Sudarshan-Lindblad approach to open quantum systems, one can perform the rotating-wave approximation with respect to the interacting or noninteracting chain Hamiltonian and respectively derive the so-called global and local master equations. The solutions of the ensuing different master equations can then be compared with the exact one, possibly sorting out the two approaches in correspondence to different timescales of the system. We derive the steady states of the open chain quantum dynamics in the two approaches and show that the behavior of fidelity between them versus interoscillator coupling depends on the two bath temperatures, revealing the existence of a temperature-dependent critical interoscillator coupling strength that determines the domain of validity of each approach. When the newly found coupling is less than this critical value, the local approach outperforms the global approach, whereas for larger interoscillator coupling, the global approach is a better approximation of the exact evolution. This critical value of interoscillator coupling depends on the two bath temperatures, which then play a crucial role in deciding the best possible approximating open dynamics.
In the following, we study the dissipative time-evolution of a quantum chain consisting of three coupled harmonic oscillators, the first and third of which weakly interact quadratically with two independent thermal baths in equilibrium at different temperatures. Due to the quadratic form of the total Hamiltonian, the unitary dynamics of the compound system is formally analytically solvable and defines a one-parameter group of Gaussian maps which enables us to solve the exact dynamics of the chain numerically. Following the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) approach to open quantum systems, one can perform the rotating wave approximation with respect to the interacting, or non-interacting chain Hamiltonian and respectively derive the so-called global and local master equations. The solutions of the ensuing different master equations can then be compared with the exact one, possibly sorting out the two approaches in correspondence to different time-scales of the system. We derive the steady states of the open chain quantum dynamics in the two approaches and show that the behaviour of fidelity between them versus inter-oscillator coupling depends on the two bath temperatures, revealing the existence of a temperature-dependent critical inter-oscillator coupling strength that determines the domain of validity of each approach. When the newly found coupling is less than this critical value, the local approach outperforms the global approach, whereas for larger inter-oscillator coupling, the global approach is a better approximation of the exact evolution. This critical value of inter-oscillator coupling depends on the two bath temperatures, which then play a crucial role in deciding the best possible approximating open dynamics.
We develop an exact analytical framework for studying the nonequilibrium dynamics of interacting quantum harmonic oscillators coupled to a modelled thermal reservoir and driven by external classical fields. The total Hamiltonian is diagonalized via successive Bogoliubov transformations, enabling a nonperturbative treatment of dissipation and driving. For two coupled oscillators with an external source applied to the first, we derive explicit energy expressions and analyze their dependence on coupling and bath parameters. The phase-space structure is characterized through Husimi Q-functions for initial separable coherent and number states. We obtain the reduced density matrix elements in the number-state basis and demonstrate that, in the absence of driving and at zero temperature, the reduced density matrix of the main system satisfies the Lindblad master equation. A generating function for the Husimi function is introduced to facilitate computation of higher-order correlations. The analysis is then generalized to n interacting oscillators, with analytical energy expressions provided and the case n=3 examined in detail. Our results establish a versatile and exact framework for driven-dissipative quantum systems, with potential applications in quantum thermodynamics, open quantum systems, and many-body physics.