We consider a strongly interacting one-dimensional (1D) Bose-Fermi mixture confined in a harmonic trap. It consists of a Tonks-Girardeau (TG) gas (1D Bose gas with repulsive hard-core interactions) and of a noninteracting Fermi gas (1D spin-aligned Fermi gas), both species interacting through hard-core repulsive interactions. Using a generalized Bose-Fermi mapping, we determine the one-body density matrices, exact particle density profiles, momentum distributions, and behavior of the mixture under 1D expansion when opening the trap. In real space, bosons and fermions do not display any phase separation: the respective density profiles extend over the same region and they both present a number of peaks equal to the total number of particles in the trap. In momentum space the bosonic component has the typical narrow TG profile, while the fermionic component shows a broad distribution with fermionic oscillations at small momenta. Due to the large boson-fermion repulsive interactions, both the bosonic and the fermionic momentum distributions decay as $C{p}^{\ensuremath{-}4}$ at large momenta, like in the case of a pure bosonic TG gas. The coefficient
C is related to the two-body density matrix and to the bosonic concentration in the mixture. When opening the trap, both momentum distributions ``fermionize'' under expansion and turn into that of a Fermi gas with a particle number equal to the total number of particles in the mixture.
After release from the trap the momentum distribution of an impenetrable gas asymptotically approaches that of a spinless noninteracting Fermi gas in the initial trap. This phenomenon is called dynamical fermionization and, very recently, has been experimentally confirmed in the case of the Lieb-Liniger model in the Tonks-Girardeau regime. We prove analytically and confirm numerically that following the removal of axial confinement the strongly interacting Bose-Fermi mixture exhibits dynamical fermionization and the asymptotical momentum distribution of each component has the same shape as its density profile at
t=0. Under a sudden change of the trap frequency to a new nonzero value the dynamics of both fermionic and bosonic momentum distributions presents characteristics which are similar to the case of single component bosons experiencing a similar quench. Our results are derived using a product representation for the correlation functions which, in addition to analytical considerations, can be implemented numerically very easily with complexity which scales polynomially in the number of particles.
Dynamical fermionization refers to the phenomenon in Tonks-Girardeau gases where, upon release from harmonic confinement, the gases' momentum density profile evolves asymptotically to that of an ideal Fermi gas in the initial trap. This phenomenon has been demonstrated theoretically in hardcore and anyonic Tonks-Girardeau gases and was recently experimentally observed in a strongly interacting Bose gas. We extend this study to a one-dimensional spinor gas of arbitrary spin in the strongly interacting regime and analytically prove that the total momentum distribution after the harmonic trap is turned off approaches that of a spinless ideal Fermi gas, while the asymptotic momentum distribution of each spin component takes the same shape of the initial real space density profile of that spin component. Our work demonstrates the rich physics arising from the interplay between the spin and the charge degrees of freedom in a spinor system.
We study a one-dimensional system of interacting anyons with short-range interactions under external confinement. This model, which describes anyons with zero-range
s-wave and
p-wave scattering interactions in the strong interaction regime, continuously interpolates between spin-polarized fermions with
p-wave interactions and free bosons. At zero temperature, the correlation functions decay exponentially with distance, with oscillations governed by the statistics parameter. The decay rate is maximal for
p-wave fermions and decreases monotonically as the statistics parameter approaches the bosonic limit, where it vanishes. The momentum distribution is asymmetric, a hallmark of one-dimensional anyons, and takes the form of a shifted Lorentzian with universal power-law tails, ${lim}_{k\ensuremath{\rightarrow}\ifmmode\pm\else\textpm\fi{}\ensuremath{\infty}}n(k)\ensuremath{\sim}C/{k}^{2}$. We prove analytically that, following release from a harmonic trap, the asymptotic momentum distribution converges to that of free bosons in the same trap, a phenomenon known as dynamical bosonization. We also establish the universality of the ground-state
n-particle reduced density matrices: their natural occupations are independent of the confining potential, while the associated natural
n-point functions for different confinements are related through a simple analytical transformation. In particular, for the one-particle reduced density matrix, we derive exact expressions for both the natural occupations and the natural orbitals at arbitrary particle number. These results extend and unify earlier partial findings for
p-wave fermions, and they provide a clear conceptual explanation of the double degeneracy observed in their spectrum.
We show by means of an exact numerical approach that the momentum distribution of a free expanding gas of hard-core bosons on a one-dimensional lattice approaches that of noninteracting fermions, acquiring a Fermi edge. Yet there is a power-law decay of the one-particle density matrix ${\ensuremath{\rho}}_{x}\ensuremath{\sim}1/\sqrt{x}$, as usual for hard-core bosons in the ground state, which accounts for a large occupation of the lowest natural orbitals for all expansion times. The fermionization of the momentum distribution function, which is not observed in equilibrium, is analyzed in detail.
We study the nonequilibrium quench dynamics of a one-dimensional anyonic gas. We focus on the integrable anyonic Lieb-Liniger model and consider the quench from noninteracting to hard-core anyons. We study the dynamics of the local properties of the system. By means of integrability-based methods, we compute analytically the one-body density matrix and the density-density correlation function at all times after the quench and in particular at infinite time. Our results show that the system evolves from an initial state where the local momentum distribution function is nonsymmetric to a steady state where it becomes symmetric. Furthermore, while the initial momentum distribution functions (and the equilibrium ones) explicitly depend on the anyonic parameter, the final ones do not. This is reminiscent of the dynamical fermionization observed in the context of free expansions after release from a confining trap.