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:2603.15554

Interaction-Enabled Hartree Fixed Points in Fermionic Resetting Dynamics

Jishad Kumar, Achilleas Lazarides, Tapio Ala-Nissila

quant-phSubmitted Mar 16, 2026• Announced Mar 17, 2026• 0 citations
arXivHTMLPDF

Abstract

In resetting dynamics, a system is repeatedly coupled to and decoupled from ancillary degrees of freedom that are reinitialized between interactions. This provides a versatile route to engineer nonequilibrium steady states and constitutes a powerful and analytically transparent framework for studying nonequilibrium dynamics in quadratic fermionic models. The baseline noninteracting resetting scheme yields an affine evolution for the subsystem single-particle density matrix (SPDM), with a clear operational interpretation: a finite environment block E mediates the interaction between the subsystem S and an ideal external thermal reservoir. In this work, we develop a controlled extension of such a framework to weakly interacting systems. We introduce a Hartree mean-field treatment of density-density interactions that preserves closure of the SPDM dynamics while producing genuinely nonlinear behavior. We further construct a completely positive (CP-safe) Gaussian Lindblad embedding that reproduces the resetting dynamics in the noninteracting limit and yields a continuous-time representation of environmental thermalization when interactions are present. Our analytical results are complemented by numerical studies of a ring segmentation geometry and a minimal two-site model, revealing interaction-enabled steady states that cannot be obtained in any purely quadratic setting. Together, these results establish a general and physically consistent route for incorporating weak interactions ino resetting-based approaches to open quantum system.

Comments
14 pages, 7 figures

Related papers

Quantum dynamics with stochastic reset

Bhaskar Mukherjee, K. Sengupta, Satya N. Majumdar

We study nonequilibrium dynamics of integrable and nonintegrable closed quantum systems whose unitary evolution is interrupted with stochastic resets, characterized by a reset rate $r$, that project the system to its initial state. We show that the steady-state density matrix of a nonintegrable system, averaged over the reset distribution, retains its off-diagonal elements for any finite $r$. Consequently a generic observable $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{O}$, whose expectation value receives a contribution from these off-diagonal elements, never thermalizes under such dynamics for any finite $r$. We demonstrate this phenomenon by exact numerical studies of experimentally realizable models of ultracold bosonic atoms in a tilted optical lattice. For integrable Dirac-like fermionic models driven periodically between such resets, the reset-averaged steady state is found to be described by a family of generalized Gibbs ensembles characterized by $r$. We also study the spread of particle density of a noninteracting one-dimensional fermionic chain, starting from an initial state where all fermions occupy the left half of the sample, while the right half is empty. When driven by resetting dynamics, the density profile approaches at long times to a nonequilibrium stationary profile that we compute exactly. We suggest concrete experiments that can possibly test our theory.

arXiv:1806.00019•Physical review. B./Physical review. B•Sep 26, 2018•108 citations

Universal nonequilibrium dynamics of pure states and density-dependent thermalization in the Sachdev-Ye-Kitaev model

Rishik Perugu, Arijit Haldar, Sumilan Banerjee

Nonequilibrium dynamics of unentangled and entangled pure states in interacting quantum systems is crucial for harnessing quantum information and to understand quantum thermalization. We develop a general Schwinger-Keldysh (SK) field theory for nonequilibrium dynamics of pure states of fermions. We apply our formalism to study the time evolution of initial density inhomogeneity and multipoint correlations of pure states in the complex Sachdev-Ye-Kitaev (SYK) models. We demonstrate a remarkable universality in the dynamics of pure states in the SYK model. We show that dynamics of almost all pure states in a fixed particle number sector is solely determined by a set of universal large-$N$ Kadanoff-Baym equations. Moreover, irrespective of the initial state the site- and disorder-averaged Green's function thermalizes instantaneously, whereas local and nonlocal Green's functions have finite thermalization rate. We provide understanding of our numerical and analytical large-$N$ results through random-matrix theory (RMT) analysis. Furthermore, we show that the thermalization of an initial pure product state in the noninteracting ${\mathrm{SYK}}_{2}$ model is independent of fermion filling and an initial density inhomogeneity decays with weak but long lived oscillations due to dephasing. In contrast, the interacting ${\mathrm{SYK}}_{q\ensuremath{\ge}4}$ model thermalizes slower than the noninteracting model and exhibits filling-dependent monotonic relaxation of initial inhomogeneity. For evolution of entangled pure states, we show that the initial entanglement is encoded in the nonlocal and/or multipoint quantum correlations that relax as the system thermalizes.

arXiv:2504.13258•Physical review. B./Physical review. B•Oct 16, 2025•2 citations

Lindblad-driven discretized leads for nonequilibrium steady-state transport in quantum impurity models: Recovering the continuum limit

Frauke Schwarz, Moshe Goldstein, Antonius Dorda +3 more

The description of interacting quantum impurity models in steady-state nonequilibrium is an open challenge for computational many-particle methods: the numerical requirement of using a finite number of lead levels and the physical requirement of describing a truly open quantum system are seemingly incompatible. One possibility to bridge this gap is the use of Lindblad-driven discretized leads (LDDL): one couples auxiliary continuous reservoirs to the discretized lead levels and represents these additional reservoirs by Lindblad terms in the Liouville equation. For quadratic models governed by Lindbladian dynamics, we present an elementary approach for obtaining correlation functions analytically. In a second part, we use this approach to explicitly discuss the conditions under which the continuum limit of the LDDL approach recovers the correct representation of thermal reservoirs. As an analytically solvable example, the nonequilibrium resonant level model is studied in greater detail. Lastly, we present ideas towards a numerical evaluation of the suggested Lindblad equation for interacting impurities based on matrix product states. In particular, we present a reformulation of the Lindblad equation, which has the useful property that the leads can be mapped onto a chain where both the Hamiltonian dynamics and the Lindblad driving are local at the same time. Moreover, we discuss the possibility to combine the Lindblad approach with a logarithmic discretization needed for the exploration of exponentially small energy scales.

arXiv:1604.02050•Physical review. B./Physical review. B•Oct 27, 2016•72 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

Gaussian phase-space representation of fermion dynamics: Beyond the time-dependent Hartree-Fock approximation

Saar Rahav, Shaul Mukamel

A Gaussian operator representation for the many-body density matrix of fermionic systems, developed by Corney and Drummond [Phys. Rev. Lett. 93, 260401 (2004)], is used to derive approximate decoupling schemes for their dynamics. In this approach the reduced single electron density-matrix elements serve as stochastic variables which satisfy an exact Fokker-Planck equation. The number of variables scales as $\ensuremath{\sim}{N}^{2}$ rather than $\ensuremath{\sim}\text{exp}(N)$ with the basis set size, and the time-dependent Hartree-Fock approximation (TDHF) is recovered in the ``classical'' limit. An approximate closed set of equations of motion for the one- and two-particle reduced density matrices provides a direct generalization of the TDHF.

arXiv:0902.3227•Physical Review B•Apr 10, 2009•10 citations

Coherence-Controlled Quantum Zeno Dynamics from Exact Reset Maps

Jishad Kumar, Achilleas Lazarides, Tapio Ala-Nissila

We develop an exact framework for quantum Zeno and anti-Zeno dynamics in a broad class of open systems, whose microscopic Hamiltonians are quadratic in bosonic or fermionic operators. We treat the environment through an exact stroboscopic resetting scheme acting at the level of the single-particle density matrix (SPDM). Within this framework, we consider two cases: a repeated-interaction (RI) protocol, in which the environment block is rethermalized and all system-environment coherences are erased after each step, and an evolving-correlation (EC) protocol, in which only the environment block is reset while system-environment coherences are preserved. For RI, we derive a general short-time Zeno law for the survival probability of a single-particle excitation and show that the corresponding decay rate scales linearly with the reset interval, implying Zeno freezing in the limit of infinitely frequent resets. Beyond the short-time regime, we formulate the RI dynamics directly in terms of the exact one-cycle propagator, which allows us to analyze finite-$τ$ anti-Zeno windows without additional approximations. For EC, we obtain a continuous-reset description in which the kept single-particle correlators obey a finite-dimensional linear differential equation. In this case the drift in the system block remains finite in the frequent-reset limit, so strict freezing is absent. We illustrate these results for a single fermionic level coupled to a semi-infinite tight-binding chain acting as a structured bath. Our results identify coherence erasure versus coherence retention as the key factor controlling the reset-induced Zeno physics.

quant-pharXiv:2603.27619•Mar 31, 2026•0 citations
In resetting dynamics, a system is repeatedly coupled to and decoupled from ancillary degrees of freedom that are reinitialized between interactions. This provides a versatile route to engineer nonequilibrium steady states and constitutes a powerful and analytically transparent framework for studying nonequilibrium dynamics in quadratic fermionic models. The baseline noninteracting resetting scheme yields an affine evolution for the subsystem single-particle density matrix (SPDM), with a clear operational interpretation: a finite environment block E mediates the interaction between the subsystem S and an ideal external thermal reservoir. In this work, we develop a controlled extension of such a framework to weakly interacting systems. We introduce a Hartree mean-field treatment of density-density interactions that preserves closure of the SPDM dynamics while producing genuinely nonlinear behavior. We further construct a completely positive (CP-safe) Gaussian Lindblad embedding that reproduces the resetting dynamics in the noninteracting limit and yields a continuous-time representation of environmental thermalization when interactions are present. Our analytical results are complemented by numerical studies of a ring segmentation geometry and a minimal two-site model, revealing interaction-enabled steady states that cannot be obtained in any purely quadratic setting. Together, these results establish a general and physically consistent route for incorporating weak interactions ino resetting-based approaches to open quantum system.
We study nonequilibrium dynamics of integrable and nonintegrable closed quantum systems whose unitary evolution is interrupted with stochastic resets, characterized by a reset rate rrr, that project the system to its initial state. We show that the steady-state density matrix of a nonintegrable system, averaged over the reset distribution, retains its off-diagonal elements for any finite rrr. Consequently a generic observable $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{O}$, whose expectation value receives a contribution from these off-diagonal elements, never thermalizes under such dynamics for any finite rrr. We demonstrate this phenomenon by exact numerical studies of experimentally realizable models of ultracold bosonic atoms in a tilted optical lattice. For integrable Dirac-like fermionic models driven periodically between such resets, the reset-averaged steady state is found to be described by a family of generalized Gibbs ensembles characterized by rrr. We also study the spread of particle density of a noninteracting one-dimensional fermionic chain, starting from an initial state where all fermions occupy the left half of the sample, while the right half is empty. When driven by resetting dynamics, the density profile approaches at long times to a nonequilibrium stationary profile that we compute exactly. We suggest concrete experiments that can possibly test our theory.
Nonequilibrium dynamics of unentangled and entangled pure states in interacting quantum systems is crucial for harnessing quantum information and to understand quantum thermalization. We develop a general Schwinger-Keldysh (SK) field theory for nonequilibrium dynamics of pure states of fermions. We apply our formalism to study the time evolution of initial density inhomogeneity and multipoint correlations of pure states in the complex Sachdev-Ye-Kitaev (SYK) models. We demonstrate a remarkable universality in the dynamics of pure states in the SYK model. We show that dynamics of almost all pure states in a fixed particle number sector is solely determined by a set of universal large-NNN Kadanoff-Baym equations. Moreover, irrespective of the initial state the site- and disorder-averaged Green's function thermalizes instantaneously, whereas local and nonlocal Green's functions have finite thermalization rate. We provide understanding of our numerical and analytical large-NNN results through random-matrix theory (RMT) analysis. Furthermore, we show that the thermalization of an initial pure product state in the noninteracting SYK2{\mathrm{SYK}}_{2}SYK2​ model is independent of fermion filling and an initial density inhomogeneity decays with weak but long lived oscillations due to dephasing. In contrast, the interacting ${\mathrm{SYK}}_{q\ensuremath{\ge}4}$ model thermalizes slower than the noninteracting model and exhibits filling-dependent monotonic relaxation of initial inhomogeneity. For evolution of entangled pure states, we show that the initial entanglement is encoded in the nonlocal and/or multipoint quantum correlations that relax as the system thermalizes.
The description of interacting quantum impurity models in steady-state nonequilibrium is an open challenge for computational many-particle methods: the numerical requirement of using a finite number of lead levels and the physical requirement of describing a truly open quantum system are seemingly incompatible. One possibility to bridge this gap is the use of Lindblad-driven discretized leads (LDDL): one couples auxiliary continuous reservoirs to the discretized lead levels and represents these additional reservoirs by Lindblad terms in the Liouville equation. For quadratic models governed by Lindbladian dynamics, we present an elementary approach for obtaining correlation functions analytically. In a second part, we use this approach to explicitly discuss the conditions under which the continuum limit of the LDDL approach recovers the correct representation of thermal reservoirs. As an analytically solvable example, the nonequilibrium resonant level model is studied in greater detail. Lastly, we present ideas towards a numerical evaluation of the suggested Lindblad equation for interacting impurities based on matrix product states. In particular, we present a reformulation of the Lindblad equation, which has the useful property that the leads can be mapped onto a chain where both the Hamiltonian dynamics and the Lindblad driving are local at the same time. Moreover, we discuss the possibility to combine the Lindblad approach with a logarithmic discretization needed for the exploration of exponentially small energy scales.
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.
A Gaussian operator representation for the many-body density matrix of fermionic systems, developed by Corney and Drummond [Phys. Rev. Lett. 93, 260401 (2004)], is used to derive approximate decoupling schemes for their dynamics. In this approach the reduced single electron density-matrix elements serve as stochastic variables which satisfy an exact Fokker-Planck equation. The number of variables scales as $\ensuremath{\sim}{N}^{2}$ rather than $\ensuremath{\sim}\text{exp}(N)$ with the basis set size, and the time-dependent Hartree-Fock approximation (TDHF) is recovered in the ``classical'' limit. An approximate closed set of equations of motion for the one- and two-particle reduced density matrices provides a direct generalization of the TDHF.
We develop an exact framework for quantum Zeno and anti-Zeno dynamics in a broad class of open systems, whose microscopic Hamiltonians are quadratic in bosonic or fermionic operators. We treat the environment through an exact stroboscopic resetting scheme acting at the level of the single-particle density matrix (SPDM). Within this framework, we consider two cases: a repeated-interaction (RI) protocol, in which the environment block is rethermalized and all system-environment coherences are erased after each step, and an evolving-correlation (EC) protocol, in which only the environment block is reset while system-environment coherences are preserved. For RI, we derive a general short-time Zeno law for the survival probability of a single-particle excitation and show that the corresponding decay rate scales linearly with the reset interval, implying Zeno freezing in the limit of infinitely frequent resets. Beyond the short-time regime, we formulate the RI dynamics directly in terms of the exact one-cycle propagator, which allows us to analyze finite-τττ anti-Zeno windows without additional approximations. For EC, we obtain a continuous-reset description in which the kept single-particle correlators obey a finite-dimensional linear differential equation. In this case the drift in the system block remains finite in the frequent-reset limit, so strict freezing is absent. We illustrate these results for a single fermionic level coupled to a semi-infinite tight-binding chain acting as a structured bath. Our results identify coherence erasure versus coherence retention as the key factor controlling the reset-induced Zeno physics.