Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2112.08344 · Journal of Statistical Mechanics Theory and Experiment

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

Thomas Barthel, Yikang Zhang

Published Nov 1, 2022• 81 citations
arXivHTMLPDF

Abstract

Abstract The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k -point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

DOI
10.1088/1742-5468/ac8e5c

Related papers

Criticality and Phase Classification for Quadratic Open Quantum Many-Body Systems

Yikang Zhang, Thomas Barthel

We study the steady states of translation-invariant open quantum many-body systems governed by Lindblad master equations, where the Hamiltonian is quadratic in the ladder operators, and the Lindblad operators are either linear or quadratic and Hermitian. These systems are called quasifree and quadratic, respectively. We find that steady states of one-dimensional systems with finite-range interactions necessarily have exponentially decaying Green's functions. For the quasifree case without quadratic Lindblad operators, we show that fermionic systems with finite-range interactions are noncritical for any number of spatial dimensions and provide bounds on the correlation lengths. Quasifree bosonic systems can be critical in D>1 dimensions. Last, we address the question of phase transitions in quadratic systems and find that, without symmetry constraints beyond invariance under single-particle basis and particle-hole transformations, all gapped Liouvillians belong to the same phase.

arXiv:2204.05346•Physical Review Letters•Sep 12, 2022•29 citations

Third quantization for bosons: symplectic diagonalization, non-Hermitian Hamiltonian, and symmetries

Steven Kim, Fabian Hassler

Abstract Open quantum systems that interact with a Markovian environment can be described by a Lindblad master equation. The generator of time-translation is given by a Liouvillian superoperator acting on the density matrix of the system. As the Fock space for a single bosonic mode is already infinite-dimensional, the diagonalization of the Liouvillian has to be done on the creation- and annihilation-superoperators, a process called ‘third quantization’. We propose a method to solve the Liouvillian for quadratic systems using a single symplectic transformation. We show that the non-Hermitian effective Hamiltonian of the system, next to incorporating the dynamics of the system, is a tool to analyze its symmetries. As an example, we use the effective Hamiltonian to formulate a -‘symmetry’ of an open system. We describe how the inclusion of source terms allows us to obtain the cumulant generating function for observables such as the photon current.

arXiv:2304.02367•Journal of Physics A Mathematical and Theoretical•Aug 17, 2023•7 citations

Time-dependent correlation functions in open quadratic fermionic systems

Pavel Kos, Tomaž Prosen

Abstract We formulate and discuss the explicit computation of dynamic correlation functions in open quadratic fermionic systems which are driven and dissipated by Lindblad jump processes that are linear in canonical fermionic operators. Dynamic correlators are interpreted in terms of a local quantum quench where the pre-quench state is a non-equilibrium steady state, i.e. a fixed point of the Liouvillian. As an example, we study the XY spin 1/2 chain and quadratic Majorana chains with boundary Lindblad driving, whose dynamics exhibit asymmetric (skewed) light cone behaviour. We also numerically treat the two-dimensional XY model and the XY spin chain with additional Dzyaloshinskii–Moriya interactions. The latter exhibits a new non-thermal phase transition which can be understood in terms of bifurcations of the quasi-particle dispersion relation. Finally, considering in some detail the periodic quadratic fermionic ring with dissipation at a single (arbitrary) site, we present analytical expressions for first order corrections (in the strength of dissipation) to the spectrum and non-equilibrium steady state (NESS) correlation functions.

arXiv:1708.06919•Journal of Statistical Mechanics Theory and Experiment•Dec 1, 2017•32 citations

Self-consistent microscopic derivation of Markovian master equations for open quadratic quantum systems

Antonio D’Abbruzzo, Davide Rossini

We provide a rigorous construction of Markovian master equations for a wide class of quantum systems that encompass quadratic models of finite size, linearly coupled to an environment modeled by a set of independent thermal baths. Our theory can be applied for both fermionic and bosonic models in any number of physical dimensions and does not require any particular spatial symmetry of the global system. We show that, for nondegenerate systems under a full secular approximation, the effective Lindblad operators are the normal modes of the system, with coupling constants that explicitly depend on the transformation matrices that diagonalize the Hamiltonian. Both the dynamics and the steady-state (guaranteed to be unique) properties can be obtained with a polynomial amount of resources in the system size. We also address the particle and energy current flowing through the system in a minimal two-bath scheme and find that they hold the structure of Landauer's formula, being thermodynamically consistent.

arXiv:2101.09303•Physical Review A•May 17, 2021•29 citations

Non-diagonal Lindblad master equations in quantum reservoir engineering

Diego N. Bernal-García, Lujun Huang, Andrey E. Miroshnichenko +1 more

Reservoir engineering has proven to be a practical approach to control open quantum systems, preserving quantum coherence by appropriately manipulating the reservoir and system-reservoir interactions. In this context, for systems comprised of different parts, it is common to describe the dynamics of a subsystem of interest by performing an adiabatic elimination of the remaining components of the system. This procedure often leads to an effective master equation for the subsystem that is not in the diagonal form of the Gorini-Kossakowski-Lindblad-Sudarshan master equation (here called diagonal Lindblad form). Instead, it has a more general structure (here called non-diagonal Lindblad form), which explicitly reveals the dissipative coupling between the various components of the subsystem. In this work, we present a set of dynamical equations for the first and second moments of the canonical variables for linear Gaussian systems, bosonic and fermionic, described by non-diagonal Lindblad master equations. Our method is efficient and allows one to obtain analytical solutions for the steady state. We supplement our findings with a review of covariance matrix methods, focusing on those related to the measurement of entanglement. Notably, our exploration yields a surprising byproduct: the Duan criterion, commonly applied to bosonic systems for verification of entanglement, is found to be equally valid for fermionic systems. We conclude with a practical example, where we revisit two-mode mechanical entanglement in an optomechanical setup. Our approach, which employs adiabatic elimination for systems governed by time-dependent Hamiltonians, opens the door to examine physical regimes that have not been explored before.

arXiv:2111.04041•Nov 7, 2021•0 citations

Third quantization of open quantum systems: Dissipative symmetries and connections to phase-space and Keldysh field-theory formulations

Alexander McDonald, Aashish A. Clerk

The connections between standard theoretical tools used to study open quantum systems can sometimes seem opaque. Whether it is a Lindblad master equation, the equation of motion for the Wigner function, or a dissipative Keldysh action, features evident in one formalism are often masked in another. Here, we reformulate the technique of third quantization in a way that explicitly connects all three methods. We first show that our formulation reveals a fundamental dissipative symmetry present in all quadratic bosonic or fermionic Lindbladians. This symmetry can then be used to easily diagonalize these models, and it provides an intuitive way to demonstrate the separation of dissipation and fluctuations in linear systems. For bosons, we then show that the Wigner function and the characteristic function can be thought of as ``wave functions'' of the density matrix in the eigenbasis of the third-quantized superoperators we introduce. The field-theory representation of the time-evolution operator in this basis is then the Keldysh path integral. To highlight the utility of our approach, we apply our version of third quantization to a dissipative nonlinear oscillator, and we use it to obtain new exact results.

arXiv:2302.14047•Physical Review Research•Aug 16, 2023•33 citations
Abstract The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k -point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.
We study the steady states of translation-invariant open quantum many-body systems governed by Lindblad master equations, where the Hamiltonian is quadratic in the ladder operators, and the Lindblad operators are either linear or quadratic and Hermitian. These systems are called quasifree and quadratic, respectively. We find that steady states of one-dimensional systems with finite-range interactions necessarily have exponentially decaying Green's functions. For the quasifree case without quadratic Lindblad operators, we show that fermionic systems with finite-range interactions are noncritical for any number of spatial dimensions and provide bounds on the correlation lengths. Quasifree bosonic systems can be critical in D>1 dimensions. Last, we address the question of phase transitions in quadratic systems and find that, without symmetry constraints beyond invariance under single-particle basis and particle-hole transformations, all gapped Liouvillians belong to the same phase.
Abstract Open quantum systems that interact with a Markovian environment can be described by a Lindblad master equation. The generator of time-translation is given by a Liouvillian superoperator acting on the density matrix of the system. As the Fock space for a single bosonic mode is already infinite-dimensional, the diagonalization of the Liouvillian has to be done on the creation- and annihilation-superoperators, a process called ‘third quantization’. We propose a method to solve the Liouvillian for quadratic systems using a single symplectic transformation. We show that the non-Hermitian effective Hamiltonian of the system, next to incorporating the dynamics of the system, is a tool to analyze its symmetries. As an example, we use the effective Hamiltonian to formulate a -‘symmetry’ of an open system. We describe how the inclusion of source terms allows us to obtain the cumulant generating function for observables such as the photon current.
Abstract We formulate and discuss the explicit computation of dynamic correlation functions in open quadratic fermionic systems which are driven and dissipated by Lindblad jump processes that are linear in canonical fermionic operators. Dynamic correlators are interpreted in terms of a local quantum quench where the pre-quench state is a non-equilibrium steady state, i.e. a fixed point of the Liouvillian. As an example, we study the XY spin 1/2 chain and quadratic Majorana chains with boundary Lindblad driving, whose dynamics exhibit asymmetric (skewed) light cone behaviour. We also numerically treat the two-dimensional XY model and the XY spin chain with additional Dzyaloshinskii–Moriya interactions. The latter exhibits a new non-thermal phase transition which can be understood in terms of bifurcations of the quasi-particle dispersion relation. Finally, considering in some detail the periodic quadratic fermionic ring with dissipation at a single (arbitrary) site, we present analytical expressions for first order corrections (in the strength of dissipation) to the spectrum and non-equilibrium steady state (NESS) correlation functions.
We provide a rigorous construction of Markovian master equations for a wide class of quantum systems that encompass quadratic models of finite size, linearly coupled to an environment modeled by a set of independent thermal baths. Our theory can be applied for both fermionic and bosonic models in any number of physical dimensions and does not require any particular spatial symmetry of the global system. We show that, for nondegenerate systems under a full secular approximation, the effective Lindblad operators are the normal modes of the system, with coupling constants that explicitly depend on the transformation matrices that diagonalize the Hamiltonian. Both the dynamics and the steady-state (guaranteed to be unique) properties can be obtained with a polynomial amount of resources in the system size. We also address the particle and energy current flowing through the system in a minimal two-bath scheme and find that they hold the structure of Landauer's formula, being thermodynamically consistent.
Reservoir engineering has proven to be a practical approach to control open quantum systems, preserving quantum coherence by appropriately manipulating the reservoir and system-reservoir interactions. In this context, for systems comprised of different parts, it is common to describe the dynamics of a subsystem of interest by performing an adiabatic elimination of the remaining components of the system. This procedure often leads to an effective master equation for the subsystem that is not in the diagonal form of the Gorini-Kossakowski-Lindblad-Sudarshan master equation (here called diagonal Lindblad form). Instead, it has a more general structure (here called non-diagonal Lindblad form), which explicitly reveals the dissipative coupling between the various components of the subsystem. In this work, we present a set of dynamical equations for the first and second moments of the canonical variables for linear Gaussian systems, bosonic and fermionic, described by non-diagonal Lindblad master equations. Our method is efficient and allows one to obtain analytical solutions for the steady state. We supplement our findings with a review of covariance matrix methods, focusing on those related to the measurement of entanglement. Notably, our exploration yields a surprising byproduct: the Duan criterion, commonly applied to bosonic systems for verification of entanglement, is found to be equally valid for fermionic systems. We conclude with a practical example, where we revisit two-mode mechanical entanglement in an optomechanical setup. Our approach, which employs adiabatic elimination for systems governed by time-dependent Hamiltonians, opens the door to examine physical regimes that have not been explored before.
The connections between standard theoretical tools used to study open quantum systems can sometimes seem opaque. Whether it is a Lindblad master equation, the equation of motion for the Wigner function, or a dissipative Keldysh action, features evident in one formalism are often masked in another. Here, we reformulate the technique of third quantization in a way that explicitly connects all three methods. We first show that our formulation reveals a fundamental dissipative symmetry present in all quadratic bosonic or fermionic Lindbladians. This symmetry can then be used to easily diagonalize these models, and it provides an intuitive way to demonstrate the separation of dissipation and fluctuations in linear systems. For bosons, we then show that the Wigner function and the characteristic function can be thought of as ``wave functions'' of the density matrix in the eigenbasis of the third-quantized superoperators we introduce. The field-theory representation of the time-evolution operator in this basis is then the Keldysh path integral. To highlight the utility of our approach, we apply our version of third quantization to a dissipative nonlinear oscillator, and we use it to obtain new exact results.