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:2607.06293 · Quantum Science and Technology

Calculating strongly correlated ground states from the non-Markovian dissipative dynamics of Gaussian fermions

Giuliano Chiriacò

quant-phcond-mat.str-elSubmitted Jul 7, 2026• Announced Jul 8, 2026• 0 citations• Revised Oct 15, 2026 (v2)
v1 · Jul 7, 2026v2 · Oct 15, 2026
arXivHTMLPDF

Abstract

We introduce a mapping between the ground state of interacting fermionic Hamiltonians and the non-equilibrium steady state of a purely dissipative open quantum system. Within the framework of third quantization, we map the Fermi-Hubbard Hamiltonian onto Lindblad jump operators acting on Majorana fermions. Remarkably, both hopping and interaction terms map onto jump operators that preserve the Gaussianity of Majorana fermions along individual quantum trajectories. As a result, the dynamics can be unravelled and each trajectory can be simulated efficiently using only two-point correlation functions, with a computational cost that scales polynomially with system size. We further show that finite particle number requires negative dissipative rates, leading to an intrinsically non-Markovian dynamics. The corresponding trajectory unravelling involves both positive and negative stochastic weights and exhibits a sign problem and large fluctuations, so that convergence requires an exponentially large number of trajectories. The overall computational cost remains exponential in system size despite the efficient Gaussian representation of individual trajectories, but is crucially dependent on the computational complexity of the non-Markovian unravelling, motivating further studies on the efficiency of such unravellings.

DOI
10.1088/2058-9565/aeb256

Related papers

Efficient Simulation of Nontrivial Dissipative Spin Chains via Stochastic Unraveling

Andrew Pocklington, Aashish A. Clerk

We present a new technique for efficiently simulating (in polynomial time) a class of one-dimensional (1D) dissipative spin chains that, when mapped to fermions, have quadratic Hamiltonians, with the only nonlinearity coming from Jordan-Wigner strings appearing in the jump operators, despite the fact that these models cannot be mapped to quadratic fermionic master equations. We show that many such Lindblad master equations admit an exact stochastic unraveling, with individual trajectories evolving as Gaussian fermionic states, even though the full master equation describes a system inequivalent to free fermions. This allows one to calculate arbitrary observables efficiently without sign problems, and with bounded sampling complexity. We utilize this new technique to study three paradigmatic dissipative effects: the melting of antiferromagnetic order in the presence of local loss, many-body subradiant phenomenon in systems with correlated loss, and nonequilibrium steady states of a 1D dissipative transverse-field Ising model. Beyond simply providing a powerful numerical technique, our method can also be used to gain both qualitative and quantitative insights into the role of interactions in these models.

arXiv:2503.23469•PRX Quantum•Aug 6, 2025•1 citation

Computable measures of non-Markovianity for Gaussian free fermion systems

Giuliano Chiriacò

Abstract We investigate measures of non-Markovianity in open quantum systems governed by Gaussian free fermionic dynamics. Standard indicators of non-Markovian behavior, such as the BLP and LFS measures, are revisited in this context. We show that for Gaussian states, trace-based distances—specifically the Hilbert–Schmidt norm—and second-order Rényi mutual information can be efficiently expressed in terms of two-point correlation functions, enabling practical computation even in systems where the full-density matrix is intractable. Crucially, this framework remains valid even when the density matrix of the system is an average over stochastic Gaussian trajectories, yielding a non-Gaussian state. We present efficient numerical protocols based on this structure and demonstrate their feasibility through a small-scale simulation. Our approach opens a scalable path to quantifying non-Markovianity in interacting or measured fermionic systems, with applications in quantum information and non-equilibrium quantum dynamics. Graphic Abstract

arXiv:2509.25953•The European Physical Journal B•Sep 1, 2025•1 citation

Open quantum dynamics of strongly coupled oscillators with multi-configuration time-dependent Hartree propagation and Markovian quantum jumps

Johan F. Triana, Felipe Herrera

Modeling the non-equilibrium dissipative dynamics of strongly interacting quantized degrees of freedom is a fundamental problem in several branches of physics and chemistry. We implement a quantum state trajectory scheme for solving Lindblad quantum master equations that describe coherent and dissipative processes for a set of strongly coupled quantized oscillators. The scheme involves a sequence of stochastic quantum jumps with transition probabilities determined by the system state and the system-reservoir dynamics. Between consecutive jumps, the wave function is propagated in a coordinate space using the multi-configuration time-dependent Hartree method. We compare this hybrid propagation methodology with exact Liouville space solutions for physical systems of interest in cavity quantum electrodynamics, demonstrating accurate results for experimentally relevant observables using a tractable number of quantum trajectories. We show the potential for solving the dissipative dynamics of finite size arrays of strongly interacting quantized oscillators with high excitation densities, a scenario that is challenging for conventional density matrix propagators due to the large dimensionality of the underlying Hilbert space.

arXiv:2208.01217•The Journal of Chemical Physics•Nov 1, 2022•7 citations

Quantum Trajectories for Time-Local Non-Lindblad Master Equations

Tobias Becker, Ché Netzer, André Eckardt

For the efficient simulation of open quantum systems, we often use quantum jump trajectories given by pure states that evolve stochastically to unravel the dynamics of the underlying master equation. In the Markovian regime, when the dynamics is described by a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation, this procedure is known as Monte Carlo wave function approach. However, beyond ultraweak system-bath coupling, the dynamics of the system is not described by an equation of GKSL type, but rather by the Redfield equation, which can be brought into pseudo-Lindblad form. Here, negative dissipation strengths prohibit the conventional approach. To overcome this problem, we propose a pseudo-Lindblad quantum trajectory (PLQT) unraveling. It does not require an effective extension of the state space, like other approaches, except for the addition of a single classical bit. We test the PLQT for the eternal non-Markovian master equation for a single qubit and an interacting Fermi-Hubbard chain coupled to a thermal bath and discuss its computational effort compared to solving the full master equation.

arXiv:2306.14876•Physical Review Letters•Oct 17, 2023•11 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

Dissipative ground state preparation in ab initio electronic structure theory

Hao-En Li, Yongtao Zhan, Lin Lin

Abstract Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.

arXiv:2411.01470•npj Quantum Information•Nov 20, 2025•8 citations
We introduce a mapping between the ground state of interacting fermionic Hamiltonians and the non-equilibrium steady state of a purely dissipative open quantum system. Within the framework of third quantization, we map the Fermi-Hubbard Hamiltonian onto Lindblad jump operators acting on Majorana fermions. Remarkably, both hopping and interaction terms map onto jump operators that preserve the Gaussianity of Majorana fermions along individual quantum trajectories. As a result, the dynamics can be unravelled and each trajectory can be simulated efficiently using only two-point correlation functions, with a computational cost that scales polynomially with system size. We further show that finite particle number requires negative dissipative rates, leading to an intrinsically non-Markovian dynamics. The corresponding trajectory unravelling involves both positive and negative stochastic weights and exhibits a sign problem and large fluctuations, so that convergence requires an exponentially large number of trajectories. The overall computational cost remains exponential in system size despite the efficient Gaussian representation of individual trajectories, but is crucially dependent on the computational complexity of the non-Markovian unravelling, motivating further studies on the efficiency of such unravellings.
We present a new technique for efficiently simulating (in polynomial time) a class of one-dimensional (1D) dissipative spin chains that, when mapped to fermions, have quadratic Hamiltonians, with the only nonlinearity coming from Jordan-Wigner strings appearing in the jump operators, despite the fact that these models cannot be mapped to quadratic fermionic master equations. We show that many such Lindblad master equations admit an exact stochastic unraveling, with individual trajectories evolving as Gaussian fermionic states, even though the full master equation describes a system inequivalent to free fermions. This allows one to calculate arbitrary observables efficiently without sign problems, and with bounded sampling complexity. We utilize this new technique to study three paradigmatic dissipative effects: the melting of antiferromagnetic order in the presence of local loss, many-body subradiant phenomenon in systems with correlated loss, and nonequilibrium steady states of a 1D dissipative transverse-field Ising model. Beyond simply providing a powerful numerical technique, our method can also be used to gain both qualitative and quantitative insights into the role of interactions in these models.
Abstract We investigate measures of non-Markovianity in open quantum systems governed by Gaussian free fermionic dynamics. Standard indicators of non-Markovian behavior, such as the BLP and LFS measures, are revisited in this context. We show that for Gaussian states, trace-based distances—specifically the Hilbert–Schmidt norm—and second-order Rényi mutual information can be efficiently expressed in terms of two-point correlation functions, enabling practical computation even in systems where the full-density matrix is intractable. Crucially, this framework remains valid even when the density matrix of the system is an average over stochastic Gaussian trajectories, yielding a non-Gaussian state. We present efficient numerical protocols based on this structure and demonstrate their feasibility through a small-scale simulation. Our approach opens a scalable path to quantifying non-Markovianity in interacting or measured fermionic systems, with applications in quantum information and non-equilibrium quantum dynamics. Graphic Abstract
Modeling the non-equilibrium dissipative dynamics of strongly interacting quantized degrees of freedom is a fundamental problem in several branches of physics and chemistry. We implement a quantum state trajectory scheme for solving Lindblad quantum master equations that describe coherent and dissipative processes for a set of strongly coupled quantized oscillators. The scheme involves a sequence of stochastic quantum jumps with transition probabilities determined by the system state and the system-reservoir dynamics. Between consecutive jumps, the wave function is propagated in a coordinate space using the multi-configuration time-dependent Hartree method. We compare this hybrid propagation methodology with exact Liouville space solutions for physical systems of interest in cavity quantum electrodynamics, demonstrating accurate results for experimentally relevant observables using a tractable number of quantum trajectories. We show the potential for solving the dissipative dynamics of finite size arrays of strongly interacting quantized oscillators with high excitation densities, a scenario that is challenging for conventional density matrix propagators due to the large dimensionality of the underlying Hilbert space.
For the efficient simulation of open quantum systems, we often use quantum jump trajectories given by pure states that evolve stochastically to unravel the dynamics of the underlying master equation. In the Markovian regime, when the dynamics is described by a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation, this procedure is known as Monte Carlo wave function approach. However, beyond ultraweak system-bath coupling, the dynamics of the system is not described by an equation of GKSL type, but rather by the Redfield equation, which can be brought into pseudo-Lindblad form. Here, negative dissipation strengths prohibit the conventional approach. To overcome this problem, we propose a pseudo-Lindblad quantum trajectory (PLQT) unraveling. It does not require an effective extension of the state space, like other approaches, except for the addition of a single classical bit. We test the PLQT for the eternal non-Markovian master equation for a single qubit and an interacting Fermi-Hubbard chain coupled to a thermal bath and discuss its computational effort compared to solving the full master equation.
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.
Abstract Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.