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

One-particle density matrix and momentum distribution of the out-of-equilibrium one-dimensional Tonks-Girardeau gas: Analytical results at large N

Stefano Scopa, Paola Ruggiero, Pasquale Calabrese, Jérôme Dubail

Published Jul 31, 2023• 10 citations
arXivHTMLPDF

Abstract

In one-dimensional (1D) quantum gases, the momentum distribution (MD) of the atoms is a standard experimental observable, routinely measured in various experimental setups. The MD is sensitive to correlations, and it is notoriously hard to compute theoretically for large numbers of atoms $N$, which often prevents direct comparison with experimental data. Here we report significant progress on this problem for the 1D Tonks-Girardeau (TG) gas in the asymptotic limit of large $N$, at zero temperature and driven out of equilibrium by a quench of the confining potential. We find an exact analytical formula for the one-particle density matrix $\ensuremath{\langle}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\mathrm{\ensuremath{\Psi}}}}^{\ifmmode\dagger\else\textdagger\fi{}}(x)\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\mathrm{\ensuremath{\Psi}}}({x}^{\ensuremath{'}})\ensuremath{\rangle}$ of the out-of-equilibrium TG gas in the $N\ensuremath{\rightarrow}\ensuremath{\infty}$ limit, valid on distances $|x\ensuremath{-}{x}^{\ensuremath{'}}|$ much larger than the interparticle distance. By comparing with time-dependent Bose-Fermi mapping numerics, we demonstrate that our analytical formula can be used to compute the out-of-equilibrium MD with great accuracy for a wide range of momenta (except in the tails of the distribution at very large momenta). For a quench from a double-well potential to a single harmonic well, which mimics a ``quantum Newton cradle'' setup, our method predicts the periodic formation of peculiar, multiply peaked, momentum distributions.

DOI
10.1103/physreva.108.013324

Related papers

Momentum Distribution Dynamics of a Tonks-Girardeau Gas: Bragg Reflections of a Quantum Many-Body Wave Packet

Robert Pezer, Hrvoje Buljan

The momentum distribution (MD) dynamics of a Tonks-Girardeau (TG) gas is studied in the context of Bragg reflections of a many-body wave packet. We find strong suppression of a Bragg reflection peak for a large and dense TG wave packet; our observation illustrates the dependence of the MD on the interactions and wave function symmetry. The MD is calculated from the reduced single-particle density matrix (RSPDM). We develop a method for calculating the RSPDM of a TG gas, which is operative for a large number of particles, and does not depend on the external potential and the state of the system. The method is based on a formula expressing the RSPDM via a dynamically evolving single-particle basis.

arXiv:cond-mat/0611584•Physical Review Letters•Jun 14, 2007•81 citations

Fermionization of a strongly interacting Bose-Fermi mixture in a one-dimensional harmonic trap

Bess Fang, Patrizia Vignolo, Christian Miniatura +1 more

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.

arXiv:0809.4419•Physical Review A•Feb 26, 2009•20 citations

Exact nonequilibrium dynamics of finite-temperature Tonks-Girardeau gases

Y. Y. Atas, Dimitri M. Gangardt, Isabelle Bouchoule +1 more

Describing finite-temperature nonequilibrium dynamics of interacting many-particle systems is a notoriously challenging problem in quantum many-body physics. Here we provide an exact solution to this problem for a system of strongly interacting bosons in one dimension in the Tonks-Girardeau regime of infinitely strong repulsive interactions. Using the Fredholm determinant approach and the Bose-Fermi mapping, we show how the problem can be reduced to a single-particle basis, wherein the finite-temperature effects enter the solution via an effective ``dressing'' of the single-particle wave functions by the Fermi-Dirac occupation factors. We demonstrate the utility of our approach and its computational efficiency in two nontrivial out-of-equilibrium scenarios: collective breathing-mode oscillations in a harmonic trap and collisional dynamics in the Newton's cradle setting involving real-time evolution in a periodic Bragg potential.

arXiv:1608.08720•Physical Review A•Apr 18, 2017•46 citations

Universal off-diagonal long-range-order behavior for a trapped Tonks-Girardeau gas

A. Colcelli, Jacopo Viti, Giuseppe Mussardo +1 more

The scaling of the largest eigenvalue ${\ensuremath{\lambda}}_{0}$ of the one-body density matrix of a system with respect to its particle number $N$ defines an exponent $\mathcal{C}$ and a coefficient $\mathcal{B}$ via the asymptotic relation ${\ensuremath{\lambda}}_{0}\ensuremath{\sim}\mathcal{B}\phantom{\rule{0.16em}{0ex}}{N}^{\mathcal{C}}$. The case $\mathcal{C}=1$ corresponds to off-diagonal long-range order. For a one-dimensional homogeneous Tonks-Girardeau gas, a well-known result also confirmed by bosonization gives instead $\mathcal{C}=1/2$. Here we investigate the inhomogeneous case, initially addressing the behavior of $\mathcal{C}$ in the presence of a general external trapping potential $V$. We argue that the value $\mathcal{C}=1/2$ characterizes the hard-core system independently of the nature of the potential $V$. We then define the exponents $\ensuremath{\gamma}$ and $\ensuremath{\beta}$, which describe the scaling of the peak of the momentum distribution with $N$ and the natural orbital corresponding to ${\ensuremath{\lambda}}_{0}$, respectively, and we derive the scaling relation $\ensuremath{\gamma}+2\ensuremath{\beta}=\mathcal{C}$. Taking as a specific case the power-law potential $V(x)\ensuremath{\propto}{x}^{2n}$, we give analytical formulas for $\ensuremath{\gamma}$ and $\ensuremath{\beta}$ as functions of $n$. Analytical predictions for the coefficient $\mathcal{B}$ are also obtained. These formulas are derived by exploiting a recent field theoretical formulation and checked against numerical results. The agreement is excellent.

arXiv:1809.01592•Physical Review A•Dec 26, 2018•17 citations

One-dimensional magnetized Tonks–Girardeau gas properties in an external magnetic field

Zhao Liang Wang, An Min Wang

With Girardeau's Fermi–Bose mapping, we have constructed the eigenstates of the Tonks–Girardeau (TG) gas in an external magnetic field. When the number of bosons N is commensurate with the number of potential cycles M , the probability of this TG gas in the ground state is bigger than the TG gas raised by Girardeau in 1960. Through the comparison of properties between this TG gas and Fermi gas, we find that the following issues are always the same: their average value of particle coordinates and potential energies, system total momenta, single-particle densities and pair distribution functions. But the reduced single-particle density matrices and their momentum distributions between them are different.

arXiv:1001.0116•Journal of Physics A Mathematical and Theoretical•Mar 29, 2010•2 citations

Harmonically confined Tonks-Girardeau gas: A simulation study based on Nelson's stochastic mechanics

Wolfgang Paul

We present a simulation study of the properties of the Tonks-Girardeau gas, a system of $N$ hard-core bosons confined in a one-dimensional harmonic trap. For this system the ground-state wave function is exactly known based on a Bose-Fermi mapping theorem. We employ Nelsons' interpretation of quantum mechanics in which this $N$-body wave function gives rise to a set of stochastic differential equations for the positions of the $N$ particles, which can be simulated by standard methods. In particular, real space densities and momentum distributions can simply be determined by time averages along particle trajectories. Employing this approach we are able to significantly extend the range of particle numbers $N$ treated numerically compared to earlier approaches, while reproducing all exactly known results for this model. We also show that for the bosons in a harmonic trap, contrary to what has been assumed so far, the momentum distribution reflects the system size scaling of the occupation numbers of the natural orbitals; that is, it can be used to decide on the presence of a Bose-Einstein condensate.

Physical Review A•Jul 6, 2012•5 citations
In one-dimensional (1D) quantum gases, the momentum distribution (MD) of the atoms is a standard experimental observable, routinely measured in various experimental setups. The MD is sensitive to correlations, and it is notoriously hard to compute theoretically for large numbers of atoms NNN, which often prevents direct comparison with experimental data. Here we report significant progress on this problem for the 1D Tonks-Girardeau (TG) gas in the asymptotic limit of large NNN, at zero temperature and driven out of equilibrium by a quench of the confining potential. We find an exact analytical formula for the one-particle density matrix \ensuremath⟨\ensuremathΨ\ifmmode^\elseˆ\fi\ifmmode†\else\textdagger\fi(x)\ensuremathΨ\ifmmode^\elseˆ\fi(x\ensuremath′)\ensuremath⟩\ensuremath{\langle}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\mathrm{\ensuremath{\Psi}}}}^{\ifmmode\dagger\else\textdagger\fi{}}(x)\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\mathrm{\ensuremath{\Psi}}}({x}^{\ensuremath{'}})\ensuremath{\rangle}\ensuremath⟨\ensuremathΨ\ifmmode^\elseˆ\fi​\ifmmode†\else\textdagger\fi(x)\ensuremathΨ\ifmmode^\elseˆ\fi​(x\ensuremath′)\ensuremath⟩ of the out-of-equilibrium TG gas in the N\ensuremath→\ensuremath∞N\ensuremath{\rightarrow}\ensuremath{\infty}N\ensuremath→\ensuremath∞ limit, valid on distances ∣x\ensuremath−x\ensuremath′∣|x\ensuremath{-}{x}^{\ensuremath{'}}|∣x\ensuremath−x\ensuremath′∣ much larger than the interparticle distance. By comparing with time-dependent Bose-Fermi mapping numerics, we demonstrate that our analytical formula can be used to compute the out-of-equilibrium MD with great accuracy for a wide range of momenta (except in the tails of the distribution at very large momenta). For a quench from a double-well potential to a single harmonic well, which mimics a ``quantum Newton cradle'' setup, our method predicts the periodic formation of peculiar, multiply peaked, momentum distributions.
The momentum distribution (MD) dynamics of a Tonks-Girardeau (TG) gas is studied in the context of Bragg reflections of a many-body wave packet. We find strong suppression of a Bragg reflection peak for a large and dense TG wave packet; our observation illustrates the dependence of the MD on the interactions and wave function symmetry. The MD is calculated from the reduced single-particle density matrix (RSPDM). We develop a method for calculating the RSPDM of a TG gas, which is operative for a large number of particles, and does not depend on the external potential and the state of the system. The method is based on a formula expressing the RSPDM via a dynamically evolving single-particle basis.
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 CCC 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.
Describing finite-temperature nonequilibrium dynamics of interacting many-particle systems is a notoriously challenging problem in quantum many-body physics. Here we provide an exact solution to this problem for a system of strongly interacting bosons in one dimension in the Tonks-Girardeau regime of infinitely strong repulsive interactions. Using the Fredholm determinant approach and the Bose-Fermi mapping, we show how the problem can be reduced to a single-particle basis, wherein the finite-temperature effects enter the solution via an effective ``dressing'' of the single-particle wave functions by the Fermi-Dirac occupation factors. We demonstrate the utility of our approach and its computational efficiency in two nontrivial out-of-equilibrium scenarios: collective breathing-mode oscillations in a harmonic trap and collisional dynamics in the Newton's cradle setting involving real-time evolution in a periodic Bragg potential.
The scaling of the largest eigenvalue ${\ensuremath{\lambda}}_{0}$ of the one-body density matrix of a system with respect to its particle number NNN defines an exponent C\mathcal{C}C and a coefficient B\mathcal{B}B via the asymptotic relation ${\ensuremath{\lambda}}_{0}\ensuremath{\sim}\mathcal{B}\phantom{\rule{0.16em}{0ex}}{N}^{\mathcal{C}}$. The case C=1\mathcal{C}=1C=1 corresponds to off-diagonal long-range order. For a one-dimensional homogeneous Tonks-Girardeau gas, a well-known result also confirmed by bosonization gives instead C=1/2\mathcal{C}=1/2C=1/2. Here we investigate the inhomogeneous case, initially addressing the behavior of C\mathcal{C}C in the presence of a general external trapping potential VVV. We argue that the value C=1/2\mathcal{C}=1/2C=1/2 characterizes the hard-core system independently of the nature of the potential VVV. We then define the exponents $\ensuremath{\gamma}$ and $\ensuremath{\beta}$, which describe the scaling of the peak of the momentum distribution with NNN and the natural orbital corresponding to ${\ensuremath{\lambda}}_{0}$, respectively, and we derive the scaling relation $\ensuremath{\gamma}+2\ensuremath{\beta}=\mathcal{C}$. Taking as a specific case the power-law potential $V(x)\ensuremath{\propto}{x}^{2n}$, we give analytical formulas for $\ensuremath{\gamma}$ and $\ensuremath{\beta}$ as functions of nnn. Analytical predictions for the coefficient B\mathcal{B}B are also obtained. These formulas are derived by exploiting a recent field theoretical formulation and checked against numerical results. The agreement is excellent.
With Girardeau's Fermi–Bose mapping, we have constructed the eigenstates of the Tonks–Girardeau (TG) gas in an external magnetic field. When the number of bosons N is commensurate with the number of potential cycles M , the probability of this TG gas in the ground state is bigger than the TG gas raised by Girardeau in 1960. Through the comparison of properties between this TG gas and Fermi gas, we find that the following issues are always the same: their average value of particle coordinates and potential energies, system total momenta, single-particle densities and pair distribution functions. But the reduced single-particle density matrices and their momentum distributions between them are different.
We present a simulation study of the properties of the Tonks-Girardeau gas, a system of NNN hard-core bosons confined in a one-dimensional harmonic trap. For this system the ground-state wave function is exactly known based on a Bose-Fermi mapping theorem. We employ Nelsons' interpretation of quantum mechanics in which this NNN-body wave function gives rise to a set of stochastic differential equations for the positions of the NNN particles, which can be simulated by standard methods. In particular, real space densities and momentum distributions can simply be determined by time averages along particle trajectories. Employing this approach we are able to significantly extend the range of particle numbers NNN treated numerically compared to earlier approaches, while reproducing all exactly known results for this model. We also show that for the bosons in a harmonic trap, contrary to what has been assumed so far, the momentum distribution reflects the system size scaling of the occupation numbers of the natural orbitals; that is, it can be used to decide on the presence of a Bose-Einstein condensate.