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

Ovidiu I. Patu

2024-08-12

Abstract:One of quantum physics's fundamental, but largely unsolved, problems is the computation of the correlation functions in many-body systems. In this paper 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 wavefunction 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-1$ multidimensional integrals, with $N$ the number of particles, which are notoriously hard to compute using conventional methods. We introduce a new 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 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 valid for both equilibrium and nonequilibrium situations.

Quantum Gases,Statistical Mechanics

What problem does this paper attempt to address?

The paper attempts to address the following issues:

1. **Calculation of Correlation Functions in Strongly Interacting Spin Gases**: A fundamental yet unresolved problem in quantum physics is the calculation of correlation functions in many-body systems. This paper investigates this issue for one-dimensional spin gases under repulsive contact interactions.

2. **Factorized Wave Functions**: In the strong coupling limit, the wave function can be factorized into charge degrees of freedom (represented as a Slater determinant of spinless fermions) and a spin part (described as a Sutherland-type spin chain). This factorization also appears in the expressions for correlation functions.

3. **High-Dimensional Integral Calculations**: Although the factorization of the wave function simplifies the calculations, traditional methods struggle due to the need to compute \(N-1\) dimensional integrals (where \(N\) is the number of particles) for local exchange coefficients and the charge component of the correlation functions. This paper proposes a new computational method that efficiently calculates these integrals.

4. **Study of Static and Dynamic Properties**: The new method allows the study of static and dynamic properties, temperature dependence, and non-equilibrium dynamics of systems with a larger number of particles, revealing new phenomena.

5. **Momentum Distribution Characteristics**: The study shows that the momentum distribution of strongly interacting trapped spin gases narrows with increasing temperature. Additionally, a simple determinant representation for the spin-incoherent state is derived, applicable to both equilibrium and non-equilibrium situations.

6. **Zero-Temperature Correlation Functions and Low-Energy Phases**: At zero temperature, multi-component bosons have a polarized ground state, and the energy of the spin part is much lower than that of the charge part under strong interactions. When the temperature is much higher than the spin excitation energy but much lower than the charge part energy, the system can be described by a spin-incoherent Luttinger liquid (SILL).

7. **Analysis of Effective Hamiltonian**: By analyzing the effective Hamiltonian, the separation of energy scales between the charge and spin parts is studied, allowing the calculation of correlation functions on the ground state manifold within a specific temperature range.

8. **Development of Numerical Methods**: The paper develops a new numerical method for efficiently calculating multi-dimensional integrals such as local exchange coefficients, single-particle density matrix elements, and single-particle densities. This method is simple and several orders of magnitude faster than existing methods.