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:2001.04615 · Physical Review A

Simulating fermions in spin-dependent potentials with spin models on an energy lattice

Michael L. Wall

Published Aug 26, 2020• 6 citations
arXivHTMLPDF

Abstract

We study spin-1/2 fermions in spin-dependent potentials under the spin-model approximation, in which interatomic collisions that change the total occupation of single-particle modes are ignored. The spin-model approximation maps the interacting fermion problem to an ensemble of lattice spin models in energy space, where spin-spin interactions are long range and spin anisotropic. We consider a Ramsey-type sequence in which quantum degenerate, initially spin-polarized fermions are rotated to a spin superposition, and the differences in potential for the two spin states induce $s$-wave interactions. Applied to such a setup, we show that the spin-model approximation qualitatively captures the collective spin dynamics for weak interactions compared to the harmonic-oscillator frequency, and often captures quantitative features at timescales much longer than would be expected from perturbation theory. We explore corrections to the spin model, and the relative importance of corrections for both collective spin and fidelity metrics when realistic anharmonic potential corrections are taken into account, demonstrating that anharmonicity can significantly improve the fidelity of the spin-model approximation. Additionally, we present numerical techniques that are useful for analysis of spin models on an energy lattice, including enacting a change of single-particle basis on a many-body state as an effective time evolution, and fitting of spatially inhomogeneous long-range interactions with exponentials. This latter technique is useful for constructing matrix product operators for use in density-matrix renormalization group analyses, and may have broader applicability within the tensor network community.

DOI
10.1103/physreva.102.023329

Related papers

Dynamics of Interacting Fermions in Spin-Dependent Potentials

Andrew Koller, Michael L. Wall, Joshua Mundinger +1 more

Recent experiments with dilute trapped Fermi gases observed that weak interactions can drastically modify spin transport dynamics and give rise to robust collective effects including global demagnetization, macroscopic spin waves, spin segregation, and spin self-rephasing. In this Letter, we develop a framework for studying the dynamics of weakly interacting fermionic gases following a spin-dependent change of the trapping potential which illuminates the interplay between spin, motion, Fermi statistics, and interactions. The key idea is the projection of the state of the system onto a set of lattice spin models defined on the single-particle mode space. Collective phenomena, including the global spreading of quantum correlations in real space, arise as a consequence of the long-ranged character of the spin model couplings. This approach achieves good agreement with prior measurements and suggests a number of directions for future experiments.

arXiv:1601.01004•Physical Review Letters•Nov 4, 2016•27 citations

Auxiliary-fermion approach to critical fluctuations in the two-dimensional quantum antiferromagnetic Heisenberg model

Jan Brinckmann, P. Wölfle

The nearest-neighbor quantum antiferromagnetic (AF) Heisenberg model for spin-$1∕2$ on a two-dimensional square lattice is studied in the auxiliary-fermion representation. Expressing spin operators by canonical fermionic particles requires a constraint on the fermion charge ${Q}_{i}=1$ on each lattice site $i$, which is imposed approximately through the thermal average. The resulting interacting fermion system is first treated in mean-field theory (MFT), which yields an AF ordered ground state and spin waves in quantitative agreement with conventional spin-wave theory. At finite temperature a self-consistent approximation beyond mean field is required in order to fulfill the Mermin-Wagner theorem. We first discuss a fully self-consistent approximation, where fermions are renormalized due to fluctuations of their spin density, in close analogy to FLEX. While static properties like the correlation length, $\ensuremath{\xi}(T)\ensuremath{\propto}\mathrm{exp}(aJ∕T)$, come out correctly, the dynamical response lacks the magnon-like peaks which would reflect the appearance of short-range order at low $T$. This drawback, which is caused by overdamping, is overcome in a ``minimal self-consistent approximation'' (MSCA), which we derive from the equations of motion. The MSCA features dynamical scaling at small energy and temperature and is qualitatively correct both in the regime of order-parameter relaxation at long wavelengths $\ensuremath{\lambda}>\ensuremath{\xi}$ and in the short-range-order regime at $\ensuremath{\lambda}<\ensuremath{\xi}$. We also discuss the impact of vertex corrections and the problem of pseudo-gap formation in the single-particle density of states due to long-range fluctuations. Finally we show that the (short-range) magnetic order in MFT and MSCA helps to fulfill the constraint on the local fermion occupancy.

arXiv:cond-mat/0405438•Physical Review B•Nov 1, 2004•7 citations

Beyond the Spin Model Approximation for Ramsey Spectroscopy

Andrew Koller, Michael E. Beverland, Alexey Vyacheslavovich Gorshkov +1 more

Ramsey spectroscopy has become a powerful technique for probing nonequilibrium dynamics of internal (pseudospin) degrees of freedom of interacting systems. In many theoretical treatments, the key to understanding the dynamics has been to assume the external (motional) degrees of freedom are decoupled from the pseudospin degrees of freedom. Determining the validity of this approximation-known as the spin model approximation-has not been addressed in detail. Here we shed light in this direction by calculating Ramsey dynamics exactly for two interacting spin-1/2 particles in a harmonic trap. We focus on s-wave-interacting fermions in quasi one- and two-dimensional geometries. We find that in one dimension the spin model assumption works well over a wide range of experimentally relevant conditions, but can fail at time scales longer than those set by the mean interaction energy. Surprisingly, in two dimensions a modified version of the spin model is exact to first order in the interaction strength. This analysis is important for a correct interpretation of Ramsey spectroscopy and has broad applications ranging from precision measurements to quantum information and to fundamental probes of many-body systems.

arXiv:1312.0887•Physical Review Letters•Mar 26, 2014•6 citations

Interacting chiral fermions on the lattice with matrix product operator norms

Jutho Haegeman, Laurens Lootens, Quinten Mortier +3 more

We develop a Hamiltonian formalism for simulating interacting chiral fermions on the lattice while preserving unitarity and locality and without breaking the chiral symmetry. The fermion doubling problem is circumvented by constructing a Fock space endowed with a semi-definite norm. When projecting our theory on the the single-particle sector, we recover the framework of Stacey fermions, and we demonstrate that the scaling limit of the free model recovers the chiral fermion field. Technically, we make use of a matrix product operator norm to mimick the boundary of a higher dimensional topological theory. As a proof of principle, we consider a single Weyl fermion on a periodic ring with Hubbard-type nearest-neighbor interactions and construct a variational generalized DMRG code to demonstrate that the ground state for large system sizes can be determined efficiently. As our tensor network approach does not exhibit any sign problem, we can add a chemical potential and study real-time evolution.

arXiv:2405.10285•May 16, 2024•1 citation

Slow semiclassical dynamics of a two-dimensional Hubbard model in disorder-free potentials

Aleksander Kaczmarek, Adam S. Sajna

The quench dynamics of the Hubbard model in tilted and harmonic potentials is discussed within the semiclassical picture. Applying the fermionic truncated Wigner approximation (fTWA), the dynamics of imbalances for charge and spin degrees of freedom is analyzed and its time evolution is compared with the exact simulations in a one-dimensional lattice. Quench from charge or spin density waves is considered. We show that introduction of harmonic and spin-dependent linear potentials sufficiently validates fTWA for longer times. Such an improvement of fTWA is also obtained for the higher-order correlations in terms of quantum Fisher information for charge and spin channels. This allows us to discuss the dynamics of larger system sizes and connect our discussion to the recently introduced Stark many-body localization. In particular, we focus on a finite two-dimensional system and show that at intermediate linear potential strength, the addition of a harmonic potential and spin dependence of the tilt results in subdiffusive dynamics, similar to that of disordered systems. Moreover, for specific values of harmonic potential, we observed phase separation of ergodic and nonergodic regions in real space. The latter fact is especially important for ultracold-atom experiments in which harmonic confinement can be easily imposed, causing a significant change in relaxation times for different lattice locations.

arXiv:2210.01082•Physical review. B./Physical review. B•Oct 9, 2023•5 citations

A generalized phase space approach for solving quantum spin dynamics

Bihui Zhu, Ana Maria Rey, Johannes Schachenmayer

Abstract Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitrary S ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it for S > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in large S > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that large S systems can feature larger entanglement than corresponding S = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail.

arXiv:1905.08782•New Journal of Physics•Jul 25, 2019•59 citations
We study spin-1/2 fermions in spin-dependent potentials under the spin-model approximation, in which interatomic collisions that change the total occupation of single-particle modes are ignored. The spin-model approximation maps the interacting fermion problem to an ensemble of lattice spin models in energy space, where spin-spin interactions are long range and spin anisotropic. We consider a Ramsey-type sequence in which quantum degenerate, initially spin-polarized fermions are rotated to a spin superposition, and the differences in potential for the two spin states induce sss-wave interactions. Applied to such a setup, we show that the spin-model approximation qualitatively captures the collective spin dynamics for weak interactions compared to the harmonic-oscillator frequency, and often captures quantitative features at timescales much longer than would be expected from perturbation theory. We explore corrections to the spin model, and the relative importance of corrections for both collective spin and fidelity metrics when realistic anharmonic potential corrections are taken into account, demonstrating that anharmonicity can significantly improve the fidelity of the spin-model approximation. Additionally, we present numerical techniques that are useful for analysis of spin models on an energy lattice, including enacting a change of single-particle basis on a many-body state as an effective time evolution, and fitting of spatially inhomogeneous long-range interactions with exponentials. This latter technique is useful for constructing matrix product operators for use in density-matrix renormalization group analyses, and may have broader applicability within the tensor network community.
Recent experiments with dilute trapped Fermi gases observed that weak interactions can drastically modify spin transport dynamics and give rise to robust collective effects including global demagnetization, macroscopic spin waves, spin segregation, and spin self-rephasing. In this Letter, we develop a framework for studying the dynamics of weakly interacting fermionic gases following a spin-dependent change of the trapping potential which illuminates the interplay between spin, motion, Fermi statistics, and interactions. The key idea is the projection of the state of the system onto a set of lattice spin models defined on the single-particle mode space. Collective phenomena, including the global spreading of quantum correlations in real space, arise as a consequence of the long-ranged character of the spin model couplings. This approach achieves good agreement with prior measurements and suggests a number of directions for future experiments.
The nearest-neighbor quantum antiferromagnetic (AF) Heisenberg model for spin-1∕21∕21∕2 on a two-dimensional square lattice is studied in the auxiliary-fermion representation. Expressing spin operators by canonical fermionic particles requires a constraint on the fermion charge Qi=1{Q}_{i}=1Qi​=1 on each lattice site iii, which is imposed approximately through the thermal average. The resulting interacting fermion system is first treated in mean-field theory (MFT), which yields an AF ordered ground state and spin waves in quantitative agreement with conventional spin-wave theory. At finite temperature a self-consistent approximation beyond mean field is required in order to fulfill the Mermin-Wagner theorem. We first discuss a fully self-consistent approximation, where fermions are renormalized due to fluctuations of their spin density, in close analogy to FLEX. While static properties like the correlation length, $\ensuremath{\xi}(T)\ensuremath{\propto}\mathrm{exp}(aJ∕T)$, come out correctly, the dynamical response lacks the magnon-like peaks which would reflect the appearance of short-range order at low TTT. This drawback, which is caused by overdamping, is overcome in a ``minimal self-consistent approximation'' (MSCA), which we derive from the equations of motion. The MSCA features dynamical scaling at small energy and temperature and is qualitatively correct both in the regime of order-parameter relaxation at long wavelengths $\ensuremath{\lambda}>\ensuremath{\xi}$ and in the short-range-order regime at $\ensuremath{\lambda}<\ensuremath{\xi}$. We also discuss the impact of vertex corrections and the problem of pseudo-gap formation in the single-particle density of states due to long-range fluctuations. Finally we show that the (short-range) magnetic order in MFT and MSCA helps to fulfill the constraint on the local fermion occupancy.
Ramsey spectroscopy has become a powerful technique for probing nonequilibrium dynamics of internal (pseudospin) degrees of freedom of interacting systems. In many theoretical treatments, the key to understanding the dynamics has been to assume the external (motional) degrees of freedom are decoupled from the pseudospin degrees of freedom. Determining the validity of this approximation-known as the spin model approximation-has not been addressed in detail. Here we shed light in this direction by calculating Ramsey dynamics exactly for two interacting spin-1/2 particles in a harmonic trap. We focus on s-wave-interacting fermions in quasi one- and two-dimensional geometries. We find that in one dimension the spin model assumption works well over a wide range of experimentally relevant conditions, but can fail at time scales longer than those set by the mean interaction energy. Surprisingly, in two dimensions a modified version of the spin model is exact to first order in the interaction strength. This analysis is important for a correct interpretation of Ramsey spectroscopy and has broad applications ranging from precision measurements to quantum information and to fundamental probes of many-body systems.
We develop a Hamiltonian formalism for simulating interacting chiral fermions on the lattice while preserving unitarity and locality and without breaking the chiral symmetry. The fermion doubling problem is circumvented by constructing a Fock space endowed with a semi-definite norm. When projecting our theory on the the single-particle sector, we recover the framework of Stacey fermions, and we demonstrate that the scaling limit of the free model recovers the chiral fermion field. Technically, we make use of a matrix product operator norm to mimick the boundary of a higher dimensional topological theory. As a proof of principle, we consider a single Weyl fermion on a periodic ring with Hubbard-type nearest-neighbor interactions and construct a variational generalized DMRG code to demonstrate that the ground state for large system sizes can be determined efficiently. As our tensor network approach does not exhibit any sign problem, we can add a chemical potential and study real-time evolution.
The quench dynamics of the Hubbard model in tilted and harmonic potentials is discussed within the semiclassical picture. Applying the fermionic truncated Wigner approximation (fTWA), the dynamics of imbalances for charge and spin degrees of freedom is analyzed and its time evolution is compared with the exact simulations in a one-dimensional lattice. Quench from charge or spin density waves is considered. We show that introduction of harmonic and spin-dependent linear potentials sufficiently validates fTWA for longer times. Such an improvement of fTWA is also obtained for the higher-order correlations in terms of quantum Fisher information for charge and spin channels. This allows us to discuss the dynamics of larger system sizes and connect our discussion to the recently introduced Stark many-body localization. In particular, we focus on a finite two-dimensional system and show that at intermediate linear potential strength, the addition of a harmonic potential and spin dependence of the tilt results in subdiffusive dynamics, similar to that of disordered systems. Moreover, for specific values of harmonic potential, we observed phase separation of ergodic and nonergodic regions in real space. The latter fact is especially important for ultracold-atom experiments in which harmonic confinement can be easily imposed, causing a significant change in relaxation times for different lattice locations.
Abstract Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitrary S ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it for S > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in large S > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that large S systems can feature larger entanglement than corresponding S = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail.