The celebrated Lindblad equation governs the nonunitary time evolution of density operators used in the description of open quantum systems. It is usually derived from the von Neumann equation for a large system, at given physical conditions, when a small subsystem is explicitly singled out and the rest of the system acts as an environment whose degrees of freedom are traced out. In the specific case of a subsystem with variable particle number, the equilibrium density operator is given by the well-known grand canonical Gibbs state. Consequently, solving the Lindblad equation in this case should automatically yield, without any additional assumptions, the corresponding density operator in the limiting case of statistical equilibrium. Current studies of the Lindblad equation with varying particle number assume, however, the grand canonical Gibbs state a priori: the chemical potential is externally imposed rather than derived from first principles, and hence the corresponding density operator is not obtained as a natural solution of the equation. In this work, we investigate the compatibility of grand canonical statistical mechanics with the derivation of the Lindblad equation. We propose an alternative and complementary approach to the current literature that consists in using a generalized system Hamiltonian which includes a term μN. In a previous paper, this empirically well-known term has been formally derived from the von Neumann equation for the specific case of equilibrium. Including μN in the system Hamiltonian leads to a modified Lindblad equation which yields the grand canonical state as a natural solution, meaning that all the quantities involved are obtained from the physics of the system without any external assumptions.
The celebrated Lindblad equation governs the nonunitary time evolution of density operators used in the description of open quantum systems. It is usually derived from the von Neumann equation for a large system, at given physical conditions, when a small subsystem is explicitly singled out and the rest of the system acts as an environment whose degrees of freedom are traced out. In the specific case of a subsystem with variable particle number, the equilibrium density operator is given by the well-known grand canonical Gibbs state. Consequently, solving the Lindblad equation in this case should automatically yield, without any additional assumptions, the corresponding density operator in the limiting case of statistical equilibrium. Current studies of the Lindblad equation with varying particle number assume, however, the grand canonical Gibbs state a priori: the chemical potential is externally imposed rather than derived from first principles, and hence the corresponding density operator is not obtained as a natural solution of the equation. In this work, we investigate the compatibility of grand canonical statistical mechanics with the derivation of the Lindblad equation. We propose an alternative and complementary approach to the current literature that consists in using a generalized system Hamiltonian which includes a term μN. In a previous paper, this empirically well-known term has been formally derived from the von Neumann equation for the specific case of equilibrium. Including μN in the system Hamiltonian leads to a modified Lindblad equation which yields the grand canonical state as a natural solution, meaning that all the quantities involved are obtained from the physics of the system without any external assumptions.
Abstract We discuss the derivation of an effective Hamiltonian for open quantum many-particle systems. The aim is to define an operator that can be used for (molecular) simulations where, through the exchange of energy and matter with the surrounding environment (reservoir), the number of particles, n, becomes a variable of the problem. The Hamiltonian is formally derived from the Von Neumann equation; specifically, we derive an n-hierarchy of equations for the density matrix, ρ ^ n , for near equilibrium situations. Such a hierarchy, in case of stationary equilibrium, delivers the standard grand canonical density matrix as it would be expected. We report that a similar Hamiltonian was conjectured, from empirical considerations, in the field of superconductivity. Thus our result also provide a formal basis for this long-standing hypothesis. Finally, an application is discussed for Path Integral simulations of molecular systems.
The Lindblad equation determines the time evolution of the density operator of open quantum systems. While valid for any system size, its use is, in practice, restricted to prototype/surrogate models with the aim of tackling specific aspects of the overall quantum complexity of a multi-atomic system. Path integral molecular dynamics (PIMD) instead provides static and dynamical quantum statistical averages of physical observables for systems in equilibrium composed of up to thousands of atoms over timescales up to nanoseconds, under the condition that short-time quantum coherence is not relevant for the properties of interest. PIMD relies on the well-established technique of molecular dynamics (MD) with its associated classical trajectories. However, it cannot describe a direct time evolution of a system and its convergence to a stationary state in situations out of equilibrium. In this work, we analyze the link between the Lindblad equation and PIMD; specifically, we will discuss how PIMD can actually be used to calculate the time evolution of ensemble-averaged physical observables and their convergence to a stationary state for situations out of equilibrium, bypassing the need of explicitly solving the Lindblad equation. Yet, at the same time, the Lindblad equation and PIMD are linked to one another through a formal relation of equivalence, which provides an argument for the consistency of PIMD results, namely the positivity of the density operator at any time. A numerical study of a prototype system, which is of interest in chemical physics, will be used to showcase the method.
The Lindblad (GKLS) master equation, which represents the mathematical form for the general evolution of a density matrix, is a versatile and widely-used tool in open quantum systems. In contrast with the typical approach of imposing additional conditions on the system, such as weak coupling or energy conservation, we explore the structure of the equation with no assumptions. We demonstrate that general quantum dynamics can be expressed through a combination of free evolution, exchanges of some physical quantities (generalised charges), not necessarily commuting with the Hamiltonian, between the system and the bath, and pure dephasing. This result comprises a novel perspective on quantum master equations, employing physical processes as elemental parts. We use it to explore the dynamics and stationary states of a two-level system and show that strong coupling, particle exchange, and non-Abelian effects all share the same physical origin. Moreover, we demonstrate that the generalised Gibbs state for all three cases contains a non-commutation term, which has not been previously considered.
We present a derivation of the Lindblad equation -an important tool for the treatment of nonunitary evolutions -that is accessible to undergraduate students in physics or mathematics with a basic background on quantum mechanics. We consider a specific case, corresponding to a very simple situation, where a primary system interacts with a bath of harmonic oscillators at zero temperature, with an interaction Hamiltonian that resembles the Jaynes-Cummings formato We start with the Born-Markov equation and, tracing out the bath degrees of freedom, we obtain an equation in the Lindblad formo The specific situation is very instructive, for it makes it easy to realize that the Lindblads represent the effect on the main system caused by the interaction with the bath, and that the Markov approximation is a fundamental condition for the emergence of the Lindbladian operator. The formal derivation of the Lindblad equation for a more general case requires the use of quantum dynamical semi-groups and broader considerations regarding the environment and temperature than we have considered in the particular case treated here.
In many cases, measurements are performed on “mixed” ensembles of realizations of a system, which cannot be associated with a single vector in its Hilbert space. In these cases the state of the system is represented by a proper “density operator.” It is instructive to associate density operators with state vectors in a vector space, termed Liouville’s space, where the Schrödinger equation is reformulated as the Liouville–Von Neumann equation. An important consequence of this equation is that the density operator of a system at equilibrium must commute with its Hamiltonian. Namely, the matrix representation of the density operator in the basis of Hamiltonian eigenstates is diagonal, where the diagonal elements are the relative populations of the system’s Hamiltonian eigenstates. The equilibrium populations are revealed by maximizing the (Von Neumann) entropy, subject to given constraints. We derive the equilibrium density operator explicitly for the cases of canonical and grand canonical ensembles.
Generalized Gibbs ensembles have been used as powerful tools to describe the steady state of integrable many-particle quantum systems after a sudden change of the Hamiltonian. Here, we demonstrate numerically that they can be used for a much broader class of problems. We consider integrable systems in the presence of weak perturbations which break both integrability and drive the system to a state far from equilibrium. Under these conditions, we show that the steady state and the time evolution on long timescales can be accurately described by a (truncated) generalized Gibbs ensemble with time-dependent Lagrange parameters, determined from simple rate equations. We compare the numerically exact time evolutions of density matrices for small systems with a theory based on block-diagonal density matrices (diagonal ensemble) and a time-dependent generalized Gibbs ensemble containing only a small number of approximately conserved quantities, using the one-dimensional Heisenberg model with perturbations described by Lindblad operators as an example.