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:2311.06331 · Physical review. B./Physical review. B

Relaxation to persistent currents in a Hubbard trimer coupled to fermionic baths

Nikodem Szpak, Gernot Schaller, Ralf Schützhold, Jürgen König

Published Sep 17, 2024• 2 citations
arXivHTMLPDF

Abstract

We consider a ring of fermionic quantum sites, modeled by the Fermi-Hubbard Hamiltonian, in which electrons can move and interact strongly via the Coulomb repulsion. The system is coupled to fermionic cold baths, which by the exchange of particles and energy induce relaxation in the system. We eliminate the environment and describe the system effectively by Lindblad master equations in various versions valid for different coupling parameter regimes. The early relaxation phase proceeds in a universal way, irrespective of the relative couplings and approximations. The system settles down to its low-energy sector and is consecutively well approximated by the Heisenberg model. We compare different Lindblad approaches, which, in the late relaxation, push the system towards different final states with opposite, extreme spin orders, from ferromagenetic to antiferromagnetic. Due to spin frustration in the trimer (a three site ring), degenerate ground states are formed by spin waves (magnons). The system described by the global coherent version of the Lindblad operators relaxes towards the final states carrying directed persistent spin currents. We numerically confirm these predictions.

DOI
10.1103/physrevb.110.115131

Related papers

Relaxation dynamics in a Hubbard dimer coupled to fermionic baths: Phenomenological description and its microscopic foundation

Eric Kleinherbers, Nikodem Szpak, Jürgen König +1 more

We study relaxation dynamics in a strongly interacting two-site Fermi-Hubbard model that is induced by coupling each site to a local fermionic bath. To derive the proper form of the Lindblad operators that enter an effective description of the system-bath coupling in different temperature regimes, we employ a diagrammatic real-time technique for the time evolution of the reduced density matrix. In spite of a local coupling to the baths, the found Lindblad operators are nonlocal in space. We compare with the local approximation, where those nonlocal effects are neglected. Furthermore, we propose an improvement on the commonly used secular approximation (rotating-wave approximation), referred to as coherent approximation, which turns out superior in all studied parameter regimes (and equivalent otherwise). We look at the relaxation dynamics for several important observables and compare the methods for early and late times in various temperature regimes.

arXiv:1910.04130•Physical review. B./Physical review. B•Mar 27, 2020•37 citations

Environment-induced decay dynamics of antiferromagnetic order in Mott-Hubbard systems

Gernot Schaller, Friedemann Queißer, Nikodem Szpak +2 more

We study the dissipative Fermi-Hubbard model in the limit of weak tunneling and strong repulsive interactions, where each lattice site is tunnel-coupled to a Markovian fermionic bath. For cold baths at intermediate chemical potentials, the Mott insulator property remains stable and we find a fast relaxation of the particle number towards half filling. On longer time scales, we find that the antiferromagnetic order of the Mott-N\'eel ground state on bipartite lattices decays, even at zero temperature. For zero and nonzero temperatures, we quantify the different relaxation time scales by means of waiting time distributions, which can be derived from an effective (non-Hermitian) Hamiltonian and obtain fully analytic expressions for the Fermi-Hubbard model on a tetramer ring.

arXiv:2112.04956•Physical review. B./Physical review. B•Mar 29, 2022•7 citations

From Heisenberg to Hubbard: An initial state for the shallow quantum simulation of correlated electrons

Bruno Murta, Joaquín Fernández-Rossier

The widespread use of the noninteracting ground state as the initial state for the digital quantum simulation of the Fermi-Hubbard model is largely due to the scarcity of alternative easy-to-prepare approximations to the exact ground state in the literature. Exploiting the fact that the spin-$\frac{1}{2}$ Heisenberg model is the effective low-energy theory of the Fermi-Hubbard model at half-filling in the strongly interacting limit, here we propose a three-step deterministic quantum routine to prepare an educated guess of the ground state of the Fermi-Hubbard model through a shallow circuit suitable for near-term quantum hardware. First, the ground state of the Heisenberg model is initialized via a hybrid variational method using an ansatz that explores only the correct symmetry subspace. Second, a general method is devised to convert a multi-spin-$\frac{1}{2}$ wave function into its fermionic version. Third, taking inspiration from the Baeriswyl ansatz, a constant-depth single-parameter layer that adds doublon-holon pairs is applied to this fermionic state. Numerical simulations on chains and ladders with up to 12 sites confirm the improvement over the noninteracting ground state of the overlap with the exact ground state for the intermediate values of the interaction strength at which quantum simulation is found to be most relevant. More broadly, the general scheme to convert a multi-spin-$\frac{1}{2}$ state into a half-filled fermionic state may bridge the gap between quantum spin models and lattice models of correlated fermions in the realm of digital quantum simulation.

arXiv:2310.16775•Physical review. B./Physical review. B•Jan 16, 2024•4 citations

Persistent Currents in Interacting Systems: Role of the Spin

G. Bouzerar, Didier Poilblanc

Persistent currents flowing through disordered mesoscopic rings threaded by a magnetic flux are investigated. Models of fermions with on-site interactions (Hubbard model) or models of spinless fermions with nearest neighbor interactions are considered on 2D cylinders with twisted boundary conditions in one direction to account for a magnetic flux. Self-consistent Hartree-Fock methods are used to treat the electron-electron interaction beyond first order. We show that the second harmonic of the current (which is relevant in the diffusive regime) is strongly suppressed by the interaction in the case of spinless fermions while it is significantly enhanced in the Hubbard model. Our data also strongly suggest that the reduction (increase) of this harmonic is related to a strong increase (reduction) of the spatial fluctuations of the charge density.

Journal de Physique I•Jul 1, 1997•2 citations

Persistent Currents in Interacting Systems: Role of the Spin

G. Bouzerar, Didier Poilblanc

Persistent currents flowing through disordered mesoscopic rings threaded by a magnetic flux are investigated. Models of fermions with on-site interactions (Hubbard model) or models of spinless fermions with nearest neighbor interactions are considered on 2D cylinders with twisted boundary conditions in one direction to account for a magnetic flux. Self-consistent Hartree-Fock methods are used to treat the electron-electron interaction beyond first order. We show that the second harmonic of the current (which is relevant in the diffusive regime) is {\it strongly suppressed} by the interaction in the case of spinless fermions while it is {\it significantly enhanced} in the Hubbard model. Our data also strongly suggest that the reduction (increase) of this harmonic is related to a strong increase (reduction) of the spacial fluctuations of the charge density.

arXiv:cond-mat/9512153•Dec 22, 1995•0 citations

Persistent current in disordered Aharonov-Bohm rings with interacting electrons

Qiang-Hua Wang, Z. D. Wang, Jian‐Xin Zhu

The role of repulsive on-site and nearest-neighbor Coulomb interactions in disordered half-filled Aharanov-Bohm rings is studied by world-line quantum Monte Carlo simulations. The diverse dependence of the equilibrium persistent current on the couplings is found to relate systematically to the magnetic phase of the model: the maximum charge stiffness (or the persistent current) coexists with the phase-transition line between the dominant charge-density-wave state and the dominant spin-density-wave state. The stiffness vanishes with an increasing departure from the transition line. Thus in the disordered rings the Coulomb interactions can enhance the charge stiffness over the noninteracting limit in such a way as to drive the system toward the phase-transition regime. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Sep 15, 1996•5 citations
We consider a ring of fermionic quantum sites, modeled by the Fermi-Hubbard Hamiltonian, in which electrons can move and interact strongly via the Coulomb repulsion. The system is coupled to fermionic cold baths, which by the exchange of particles and energy induce relaxation in the system. We eliminate the environment and describe the system effectively by Lindblad master equations in various versions valid for different coupling parameter regimes. The early relaxation phase proceeds in a universal way, irrespective of the relative couplings and approximations. The system settles down to its low-energy sector and is consecutively well approximated by the Heisenberg model. We compare different Lindblad approaches, which, in the late relaxation, push the system towards different final states with opposite, extreme spin orders, from ferromagenetic to antiferromagnetic. Due to spin frustration in the trimer (a three site ring), degenerate ground states are formed by spin waves (magnons). The system described by the global coherent version of the Lindblad operators relaxes towards the final states carrying directed persistent spin currents. We numerically confirm these predictions.
We study relaxation dynamics in a strongly interacting two-site Fermi-Hubbard model that is induced by coupling each site to a local fermionic bath. To derive the proper form of the Lindblad operators that enter an effective description of the system-bath coupling in different temperature regimes, we employ a diagrammatic real-time technique for the time evolution of the reduced density matrix. In spite of a local coupling to the baths, the found Lindblad operators are nonlocal in space. We compare with the local approximation, where those nonlocal effects are neglected. Furthermore, we propose an improvement on the commonly used secular approximation (rotating-wave approximation), referred to as coherent approximation, which turns out superior in all studied parameter regimes (and equivalent otherwise). We look at the relaxation dynamics for several important observables and compare the methods for early and late times in various temperature regimes.
We study the dissipative Fermi-Hubbard model in the limit of weak tunneling and strong repulsive interactions, where each lattice site is tunnel-coupled to a Markovian fermionic bath. For cold baths at intermediate chemical potentials, the Mott insulator property remains stable and we find a fast relaxation of the particle number towards half filling. On longer time scales, we find that the antiferromagnetic order of the Mott-N\'eel ground state on bipartite lattices decays, even at zero temperature. For zero and nonzero temperatures, we quantify the different relaxation time scales by means of waiting time distributions, which can be derived from an effective (non-Hermitian) Hamiltonian and obtain fully analytic expressions for the Fermi-Hubbard model on a tetramer ring.
The widespread use of the noninteracting ground state as the initial state for the digital quantum simulation of the Fermi-Hubbard model is largely due to the scarcity of alternative easy-to-prepare approximations to the exact ground state in the literature. Exploiting the fact that the spin-12\frac{1}{2}21​ Heisenberg model is the effective low-energy theory of the Fermi-Hubbard model at half-filling in the strongly interacting limit, here we propose a three-step deterministic quantum routine to prepare an educated guess of the ground state of the Fermi-Hubbard model through a shallow circuit suitable for near-term quantum hardware. First, the ground state of the Heisenberg model is initialized via a hybrid variational method using an ansatz that explores only the correct symmetry subspace. Second, a general method is devised to convert a multi-spin-12\frac{1}{2}21​ wave function into its fermionic version. Third, taking inspiration from the Baeriswyl ansatz, a constant-depth single-parameter layer that adds doublon-holon pairs is applied to this fermionic state. Numerical simulations on chains and ladders with up to 12 sites confirm the improvement over the noninteracting ground state of the overlap with the exact ground state for the intermediate values of the interaction strength at which quantum simulation is found to be most relevant. More broadly, the general scheme to convert a multi-spin-12\frac{1}{2}21​ state into a half-filled fermionic state may bridge the gap between quantum spin models and lattice models of correlated fermions in the realm of digital quantum simulation.
Persistent currents flowing through disordered mesoscopic rings threaded by a magnetic flux are investigated. Models of fermions with on-site interactions (Hubbard model) or models of spinless fermions with nearest neighbor interactions are considered on 2D cylinders with twisted boundary conditions in one direction to account for a magnetic flux. Self-consistent Hartree-Fock methods are used to treat the electron-electron interaction beyond first order. We show that the second harmonic of the current (which is relevant in the diffusive regime) is strongly suppressed by the interaction in the case of spinless fermions while it is significantly enhanced in the Hubbard model. Our data also strongly suggest that the reduction (increase) of this harmonic is related to a strong increase (reduction) of the spatial fluctuations of the charge density.
Persistent currents flowing through disordered mesoscopic rings threaded by a magnetic flux are investigated. Models of fermions with on-site interactions (Hubbard model) or models of spinless fermions with nearest neighbor interactions are considered on 2D cylinders with twisted boundary conditions in one direction to account for a magnetic flux. Self-consistent Hartree-Fock methods are used to treat the electron-electron interaction beyond first order. We show that the second harmonic of the current (which is relevant in the diffusive regime) is {\it strongly suppressed} by the interaction in the case of spinless fermions while it is {\it significantly enhanced} in the Hubbard model. Our data also strongly suggest that the reduction (increase) of this harmonic is related to a strong increase (reduction) of the spacial fluctuations of the charge density.
The role of repulsive on-site and nearest-neighbor Coulomb interactions in disordered half-filled Aharanov-Bohm rings is studied by world-line quantum Monte Carlo simulations. The diverse dependence of the equilibrium persistent current on the couplings is found to relate systematically to the magnetic phase of the model: the maximum charge stiffness (or the persistent current) coexists with the phase-transition line between the dominant charge-density-wave state and the dominant spin-density-wave state. The stiffness vanishes with an increasing departure from the transition line. Thus in the disordered rings the Coulomb interactions can enhance the charge stiffness over the noninteracting limit in such a way as to drive the system toward the phase-transition regime. \textcopyright{} 1996 The American Physical Society.