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

Numerical methods and analytic results for one-dimensional strongly interacting spinor gases

Ovidiu I. Pâţu

Published Nov 5, 2024• 1 citation
arXivHTMLPDF

Abstract

One of quantum physics' fundamental, but largely unsolved, problems is the computation of the correlation functions in many-body systems. We address this problem in the case of one-dimensional spinor gases with repulsive contact interactions in the presence of a confining potential. We take advantage of the fact that in the strong-coupling limit, the wave function factorizes with the charge degrees of freedom expressed as a Slater determinant of spinless fermions and the spin sector described by a spin chain of Sutherland type with exchange coefficients that depend only on the trapping potential. This factorization is also present in the expressions for the correlation functions. Still, analytical and numerical investigations were hindered by the fact that the local exchange coefficients and the charge component of the correlators are expressed as $N\ensuremath{-}1$ multidimensional integrals, with $N$ the number of particles, which are notoriously hard to compute using conventional methods. We introduce a different approach to calculating these integrals that is extremely simple, scales polynomially with the number of particles, and is several orders of magnitude faster than the previous methods reported in the literature. This allows us to investigate the static and dynamic properties, temperature dependence, and nonequilibrium dynamics for systems with a larger number of particles than previously considered and discover new phenomena. We show that contrary to natural expectations, the momentum distribution of strongly interacting trapped spinor gases becomes narrower as we increase the temperature and derive simple determinant representations for the correlators in the spin incoherent regime that is valid for both equilibrium and nonequilibrium situations. For homogeneous systems of impenetrable particles at zero temperature, we analytically compute the large distance asymptotics of the correlators, finding that the leading term of the asymptotic expansion is proportional to ${e}^{\ensuremath{-}x{k}_{F}ln\ensuremath{\kappa}/\ensuremath{\pi}}$, with $\ensuremath{\kappa}$ the number of components of the system, which precludes any singularity of the momentum distribution.

DOI
10.1103/physreva.110.053303

Related papers

Correlation functions of one-dimensional strongly interacting two-component gases

Ovidiu I. Pâţu

We address the problem of calculating the correlation functions of one-dimensional two-component gases with strong repulsive contact interactions. The model considered in this paper describes particles with fractional statistics and in appropriate limits reduces to the Gaudin-Yang model or the spinor Bose gas. In the case of impenetrable particles we derive a Fredholm determinant representation for the temperature-, time-, and space-dependent correlation functions which is very easy to implement numerically and constitutes the starting point for the analytical investigation of the asymptotics. Making use of this determinant representation and the solution of an associated Riemann-Hilbert problem we derive the low-energy asymptotics of the correlators in the spin-incoherent regime characterized by near ground-state charge degrees of freedom but a highly thermally disordered spin sector. The asymptotics present features reminiscent of spin-charge separation with the spin part exponentially decaying in space separation and oscillating with a period proportional to the statistics parameter, while the charge part presents scaling with anomalous exponents which cannot be described by any unitary conformal field theory. The momentum distribution and the Fourier transform of the dynamical Green's function are asymmetrical for arbitrary statistics, a direct consequence of the broken space-reversal symmetry. Due to the exponential decay the momentum distribution $n(k)$ at zero temperature does not present algebraic singularities but the tails obey the universal decay ${lim}_{k\ensuremath{\rightarrow}\ifmmode\pm\else\textpm\fi{}\ensuremath{\infty}}n(k)\ensuremath{\sim}C/{k}^{4}$ with the amplitude $C$ given by Tan's contact. As a function of the statistics parameter the contact is a monotonic function reaching its minimum for the fermionic system and the maximum for the bosonic system.

arXiv:1910.10042•Physical Review A•Dec 23, 2019•22 citations

Momentum distributions and numerical methods for strongly interacting one-dimensional spinor gases

F. Deuretzbacher, D. Becker, Luis Santos

One-dimensional spinor gases with strong $\ensuremath{\delta}$ interaction fermionize and form a spin chain. The spatial degrees of freedom of this atom chain can be described by a mapping to spinless noninteracting fermions and the spin degrees of freedom are described by a spin-chain model with nearest-neighbor interactions. Here, we compute momentum and occupation-number distributions of up to 16 strongly interacting spinor fermions and bosons as a function of their spin imbalance, the strength of an externally applied magnetic field gradient, the length of their spin, and for different excited states of the multiplet. We show that the ground-state momentum distributions resemble those of the corresponding noninteracting systems, apart from flat background distributions, which extend to high momenta. Moreover, we show that the spin order of the spin chain---in particular antiferromagnetic spin order---may be deduced from the momentum and occupation-number distributions of the system. Finally, we present efficient numerical methods for the calculation of the single-particle densities and one-body density matrix elements and of the local exchange coefficients of the spin chain for large systems containing more than 20 strongly interacting particles in arbitrary confining potentials.

arXiv:1602.06816•Physical Review A•Aug 4, 2016•42 citations

Strongly Interacting One-Dimensional Spinor Quantum Gases

Yang Li

Quantum many-body physics has been studied for many decades. A lot of intriguing phenomena have been observed and theories have been developed. Yet many problems remain unsolved, which is largely due to a lack of general efficient classical computation method. One-dimensional (1D) systems have also drawn much attention over the past few decades. One reason is that many unique strongly correlated quantum phenomena only appear in low dimensions. Another reason is that many exact results can be obtained in 1D for cross benchmarking, such as Bethe Ansatz, Bosonization, Bose-Fermi mapping, etc. In addition, many numerical methods, such as Matrix-Product-State based and Monte Carlo methods, work most efficiently in 1D. In recent years, 1D quantum gases have been realized in many cold atom labs, providing experimental motivation for their studies. One of the most mysterious assumptions about a system of identical quantum particles is that the wavefunction must be symmetric (for bosons) or anti-symmetric (for fermions). This means that no two fermions can occupy the same state, while multiple occupancy is allowed for bosons. Bosons and fermions are therefore generally very different. But in 1D this distinction could become rather subtle, provided that the bosonic multiple occupancy is suppressed, which can happen if strong repulsion exists between them. However, when considering spin degrees of freedom, things become more complicated. In this work, we develop a generalized Bose-Fermi mapping theory, under which, the charge degrees of freedom is mapped to a spinless fermions, while the spin degrees of freedom to an effective spin chain. This mapping works for arbitrary spin, arbitrary trapping potentials, arbitrary spin-dependent interaction potential, and for either spinor bosons or fermions. In the strong interaction limit, the wavefunction of the system is represented by the strong coupling ansatz wavefunction (SCAW), which takes the form of a direct product of spinless fermions wavefunction and spin wavefunction. Using this mapping technique, we study the dynamics and collective modes of the system. Furthermore, we develop a very efficient method to calculate the one-body density matrix, from which we can calculate the momentum distribution of the system.

•Feb 5, 2019•0 citations

Dynamical Fermionization in One-Dimensional Spinor Gases at Finite Temperature

Ovidiu I. Pâţu

Following the removal of axial confinement the momentum distribution of a Tonks-Girardeau gas approaches that of a system of noninteracting spinless fermions in the initial harmonic trap. This phenomenon, called dynamical fermionization, has been experimentally confirmed in the case of the Lieb-Liniger model and theoretically predicted in the case of multicomponent systems at zero temperature. We prove analytically that for all spinor gases with strong repulsive contact interactions at finite temperature the momentum distribution after release from the trap asymptotically approaches that of a system of spinless fermions at the same temperature but with a renormalized chemical potential which depends on the number of components of the spinor system. In the case of the Gaudin-Yang model we check numerically our analytical predictions using the results obtained from a nonequilibrium generalization of Lenard's formula describing the time evolution of the field-field correlators.

arXiv:2211.01975•Physical Review Letters•Apr 21, 2023•5 citations

Expansion of one-dimensional spinor gases from power-law traps

Ovidiu I. Pâţu, Gianni Aupetit-Diallo

Free expansion following the removal of axial confinement represents a fundamental nonequilibrium scenario in the study of many-body ultracold gases. Using the stationary phase approximation, we analytically demonstrate that for all one-dimensional spinor gases with repulsive contact interactions, whether bosonic or fermionic, the asymptotic density and momentum distribution can be directly determined from the quasimomentum distribution (Bethe rapidities) of the trapped gas. We efficiently obtain the quasimomentum distribution numerically by solving the integral equations that characterize the ground state of the integrable system within the local density approximation. Additionally, we derive analytical solutions for both weakly and strongly interacting regimes. Unlike in bosonic gases, where rapidity distributions and density profiles vary significantly across interaction regimes, fermionic gases maintain similar profiles in both weakly and strongly interacting limits. Notably, the gas expands self-similarly only when released from a harmonic trap. For other power-law trapping potentials, the asymptotic density profile is strongly influenced by the initial confinement geometry. Our results extend readily to Bose-Fermi mixtures and finite temperatures.

arXiv:2505.05186•Physical Review A•Jul 17, 2025•0 citations

One-body correlations and momentum distributions of trapped one-dimensional Bose gases at finite temperature

Attila Takács, Yicheng Zhang, Pasquale Calabrese +3 more

We introduce a general approximate method for calculating the one-body correlations and the momentum distributions of one-dimensional Bose gases at finite interaction strengths and temperatures trapped in smooth confining potentials. Our method combines asymptotic techniques for the long-distance behavior of the gas (similar to Luttinger liquid theory) with known short-distance expansions. We derive analytical results for the limiting cases of strong and weak interactions and provide a general procedure for calculating one-body correlations at any interaction strength. A step-by-step explanation of the numerical method used to compute Green's functions (needed as input to our theory) is included. We benchmark our method against exact numerical calculations and compare its predictions to recent experimental results.

arXiv:2409.16929•Physical Review A•Mar 20, 2025•1 citation
One of quantum physics' fundamental, but largely unsolved, problems is the computation of the correlation functions in many-body systems. We address this problem in the case of one-dimensional spinor gases with repulsive contact interactions in the presence of a confining potential. We take advantage of the fact that in the strong-coupling limit, the wave function factorizes with the charge degrees of freedom expressed as a Slater determinant of spinless fermions and the spin sector described by a spin chain of Sutherland type with exchange coefficients that depend only on the trapping potential. This factorization is also present in the expressions for the correlation functions. Still, analytical and numerical investigations were hindered by the fact that the local exchange coefficients and the charge component of the correlators are expressed as N\ensuremath−1N\ensuremath{-}1N\ensuremath−1 multidimensional integrals, with NNN the number of particles, which are notoriously hard to compute using conventional methods. We introduce a different approach to calculating these integrals that is extremely simple, scales polynomially with the number of particles, and is several orders of magnitude faster than the previous methods reported in the literature. This allows us to investigate the static and dynamic properties, temperature dependence, and nonequilibrium dynamics for systems with a larger number of particles than previously considered and discover new phenomena. We show that contrary to natural expectations, the momentum distribution of strongly interacting trapped spinor gases becomes narrower as we increase the temperature and derive simple determinant representations for the correlators in the spin incoherent regime that is valid for both equilibrium and nonequilibrium situations. For homogeneous systems of impenetrable particles at zero temperature, we analytically compute the large distance asymptotics of the correlators, finding that the leading term of the asymptotic expansion is proportional to e\ensuremath−xkFln\ensuremathκ/\ensuremathπ{e}^{\ensuremath{-}x{k}_{F}ln\ensuremath{\kappa}/\ensuremath{\pi}}e\ensuremath−xkF​ln\ensuremathκ/\ensuremathπ, with \ensuremathκ\ensuremath{\kappa}\ensuremathκ the number of components of the system, which precludes any singularity of the momentum distribution.
We address the problem of calculating the correlation functions of one-dimensional two-component gases with strong repulsive contact interactions. The model considered in this paper describes particles with fractional statistics and in appropriate limits reduces to the Gaudin-Yang model or the spinor Bose gas. In the case of impenetrable particles we derive a Fredholm determinant representation for the temperature-, time-, and space-dependent correlation functions which is very easy to implement numerically and constitutes the starting point for the analytical investigation of the asymptotics. Making use of this determinant representation and the solution of an associated Riemann-Hilbert problem we derive the low-energy asymptotics of the correlators in the spin-incoherent regime characterized by near ground-state charge degrees of freedom but a highly thermally disordered spin sector. The asymptotics present features reminiscent of spin-charge separation with the spin part exponentially decaying in space separation and oscillating with a period proportional to the statistics parameter, while the charge part presents scaling with anomalous exponents which cannot be described by any unitary conformal field theory. The momentum distribution and the Fourier transform of the dynamical Green's function are asymmetrical for arbitrary statistics, a direct consequence of the broken space-reversal symmetry. Due to the exponential decay the momentum distribution n(k)n(k)n(k) at zero temperature does not present algebraic singularities but the tails obey the universal decay ${lim}_{k\ensuremath{\rightarrow}\ifmmode\pm\else\textpm\fi{}\ensuremath{\infty}}n(k)\ensuremath{\sim}C/{k}^{4}$ with the amplitude CCC given by Tan's contact. As a function of the statistics parameter the contact is a monotonic function reaching its minimum for the fermionic system and the maximum for the bosonic system.
One-dimensional spinor gases with strong $\ensuremath{\delta}$ interaction fermionize and form a spin chain. The spatial degrees of freedom of this atom chain can be described by a mapping to spinless noninteracting fermions and the spin degrees of freedom are described by a spin-chain model with nearest-neighbor interactions. Here, we compute momentum and occupation-number distributions of up to 16 strongly interacting spinor fermions and bosons as a function of their spin imbalance, the strength of an externally applied magnetic field gradient, the length of their spin, and for different excited states of the multiplet. We show that the ground-state momentum distributions resemble those of the corresponding noninteracting systems, apart from flat background distributions, which extend to high momenta. Moreover, we show that the spin order of the spin chain---in particular antiferromagnetic spin order---may be deduced from the momentum and occupation-number distributions of the system. Finally, we present efficient numerical methods for the calculation of the single-particle densities and one-body density matrix elements and of the local exchange coefficients of the spin chain for large systems containing more than 20 strongly interacting particles in arbitrary confining potentials.
Quantum many-body physics has been studied for many decades. A lot of intriguing phenomena have been observed and theories have been developed. Yet many problems remain unsolved, which is largely due to a lack of general efficient classical computation method. One-dimensional (1D) systems have also drawn much attention over the past few decades. One reason is that many unique strongly correlated quantum phenomena only appear in low dimensions. Another reason is that many exact results can be obtained in 1D for cross benchmarking, such as Bethe Ansatz, Bosonization, Bose-Fermi mapping, etc. In addition, many numerical methods, such as Matrix-Product-State based and Monte Carlo methods, work most efficiently in 1D. In recent years, 1D quantum gases have been realized in many cold atom labs, providing experimental motivation for their studies. One of the most mysterious assumptions about a system of identical quantum particles is that the wavefunction must be symmetric (for bosons) or anti-symmetric (for fermions). This means that no two fermions can occupy the same state, while multiple occupancy is allowed for bosons. Bosons and fermions are therefore generally very different. But in 1D this distinction could become rather subtle, provided that the bosonic multiple occupancy is suppressed, which can happen if strong repulsion exists between them. However, when considering spin degrees of freedom, things become more complicated. In this work, we develop a generalized Bose-Fermi mapping theory, under which, the charge degrees of freedom is mapped to a spinless fermions, while the spin degrees of freedom to an effective spin chain. This mapping works for arbitrary spin, arbitrary trapping potentials, arbitrary spin-dependent interaction potential, and for either spinor bosons or fermions. In the strong interaction limit, the wavefunction of the system is represented by the strong coupling ansatz wavefunction (SCAW), which takes the form of a direct product of spinless fermions wavefunction and spin wavefunction. Using this mapping technique, we study the dynamics and collective modes of the system. Furthermore, we develop a very efficient method to calculate the one-body density matrix, from which we can calculate the momentum distribution of the system.
Following the removal of axial confinement the momentum distribution of a Tonks-Girardeau gas approaches that of a system of noninteracting spinless fermions in the initial harmonic trap. This phenomenon, called dynamical fermionization, has been experimentally confirmed in the case of the Lieb-Liniger model and theoretically predicted in the case of multicomponent systems at zero temperature. We prove analytically that for all spinor gases with strong repulsive contact interactions at finite temperature the momentum distribution after release from the trap asymptotically approaches that of a system of spinless fermions at the same temperature but with a renormalized chemical potential which depends on the number of components of the spinor system. In the case of the Gaudin-Yang model we check numerically our analytical predictions using the results obtained from a nonequilibrium generalization of Lenard's formula describing the time evolution of the field-field correlators.
Free expansion following the removal of axial confinement represents a fundamental nonequilibrium scenario in the study of many-body ultracold gases. Using the stationary phase approximation, we analytically demonstrate that for all one-dimensional spinor gases with repulsive contact interactions, whether bosonic or fermionic, the asymptotic density and momentum distribution can be directly determined from the quasimomentum distribution (Bethe rapidities) of the trapped gas. We efficiently obtain the quasimomentum distribution numerically by solving the integral equations that characterize the ground state of the integrable system within the local density approximation. Additionally, we derive analytical solutions for both weakly and strongly interacting regimes. Unlike in bosonic gases, where rapidity distributions and density profiles vary significantly across interaction regimes, fermionic gases maintain similar profiles in both weakly and strongly interacting limits. Notably, the gas expands self-similarly only when released from a harmonic trap. For other power-law trapping potentials, the asymptotic density profile is strongly influenced by the initial confinement geometry. Our results extend readily to Bose-Fermi mixtures and finite temperatures.
We introduce a general approximate method for calculating the one-body correlations and the momentum distributions of one-dimensional Bose gases at finite interaction strengths and temperatures trapped in smooth confining potentials. Our method combines asymptotic techniques for the long-distance behavior of the gas (similar to Luttinger liquid theory) with known short-distance expansions. We derive analytical results for the limiting cases of strong and weak interactions and provide a general procedure for calculating one-body correlations at any interaction strength. A step-by-step explanation of the numerical method used to compute Green's functions (needed as input to our theory) is included. We benchmark our method against exact numerical calculations and compare its predictions to recent experimental results.