An efficient theoretical technique for calculating the one-body density matrix and momentum distributions is developed for general one-dimensional spinor quantum gases in the strong-interaction regime by exploring the connection between the density matrix and one-dimensional spinless hard-core anyon gases. The technique is expected to provide significant new insights into the properties of the one-dimensional spinor gases.
An efficient theoretical technique for calculating the one-body density matrix and momentum distributions is developed for general one-dimensional spinor quantum gases in the strong-interaction regime by exploring the connection between the density matrix and one-dimensional spinless hard-core anyon gases. The technique is expected to provide significant new insights into the properties of the one-dimensional spinor gases.
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.
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.
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.
Abstract Quantum many-body systems in one dimension (1D) exhibit some peculiar properties. In this article, we review some of our work on strongly interacting 1D spinor quantum gas. First, we discuss a generalized Bose–Fermi mapping that maps the charge degrees of freedom to a spinless Fermi gas and the spin degrees of freedom to a spin chain model. This also maps the strongly interacting system into a weakly interacting one, which is amenable for perturbative calculations. Next, based on this mapping, we construct an ansatz wavefunction for the strongly interacting system, using which many physical quantities can be conveniently calculated. We showcase the usage of this ansatz wavefunction by considering the collective excitations and quench dynamics of a harmonically trapped system.
Ultra-cold atomic Fermi gases at nano-Kelvin temperatures provide a new paradigm for the foundations of quantum many-body theory. Here, we summarize our recent theoretical efforts to understand the intriguing properties of these strongly correlated Fermi gases. A number of powerful theoretical techniques have been developed, including (i) a quantitative diagrammatic theory taken to infinite order; (ii) a virial expansion at high temperature; (iii) exact solutions for one-dimensional many-body systems; and (iv) new quantum simulations using the Gaussian Fermi representation. We have employed these techniques to predict universal thermodynamics in the strongly interacting limit and exotic quantum superfluid phases with imbalanced spin populations. More insights into strongly interacting Fermi gases are anticipated, by extending the Gaussian Fermi simulations to new regimes.