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

Quantum Monte Carlo study of one-dimensional trapped fermions with attractive contact interactions

Michele Casula, David M. Ceperley, Erich J. Mueller

Published Sep 11, 2008• 116 citations
arXivHTMLPDF

Abstract

Using exact continuous quantum Monte Carlo techniques, we study the zero- and finite-temperature properties of a system of harmonically trapped one-dimensional spin-$\frac{1}{2}$ fermions with short-range interactions. Motivated by experimental searches for modulated Fulde-Ferrel-Larkin-Ovchinikov states, we systematically examine the impact of a spin imbalance on the density profiles. We quantify the accuracy of the Thomas-Fermi approximation, finding that for sufficiently large particle numbers $(N\ensuremath{\gtrsim}100)$ it quantitatively reproduces most features of the exact density profile. The Thomas-Fermi approximation fails to capture small Friedel-like spin and density oscillations and overestimates the size of the fully paired region in the outer shell of the trap. Based on our results, we suggest a range of experimentally tunable parameters to maximize the visibility of the double-shell structure of the system and the Fulde-Ferrel-Larkin-Ovchinikov state in the one-dimensional harmonic trap. Furthermore, we analyze the fingerprints of the attractive contact interactions in the features of the momentum and pair momentum distributions.

DOI
10.1103/physreva.78.033607

Related papers

Ground-state energy, density profiles, and momentum distribution of attractively interacting 1D Fermi gases with hard-wall boundaries: a Monte Carlo study

Shill,, Drut,, Porter Porter +1 more

Motivated by the realization of hard-wall boundary conditions in experiments with ultracold atoms, we investigate the ground-state properties of spin-1/2 fermions with attractive interactions in a one-dimensional box.We use lattice Monte Carlo methods to determine essential quantities like the energy, which we compute as a function of coupling strength and particle number in the regime from few to many particles.Many-fermion systems bound by hard walls display non-trivial density profiles characterized by so-called Friedel oscillations (which are similar to those observed in harmonic traps).In non-interacting systems, the characteristic length scale of the oscillations is set by (2k F ) -1 , where k F is the Fermi momentum, while repulsive interactions tend to generate Wignercrystal oscillations of period (4k F ) -1 .Based on the non-interacting result, we find a remarkably simple parametrization of the density profiles of the attractively interacting case, which we generalize to the one-body density matrix.While the total momentum is not a conserved quantity in the presence of hard walls, the magnitude of the momentum does provide a good quantum number.We are therefore able to provide a detailed characterization of the (quasi-)momentum distribution, which displays rather robust discontinuity at the Fermi surface.In addition, we determine the spatially varying on-site density-density correlation, which in turn yields Tan's contact density and, upon integration, Tan's contact.As is well known, the latter fully determines the short-range correlations and plays a crucial role in a multitude of equilibrium and non-equilibrium sum rules.

•Apr 18, 2020•0 citations

Ground-state energy, density profiles, and momentum distribution of attractively interacting 1D Fermi gases with hard-wall boundaries: a Monte Carlo study

J. R. McKenney, C R Shill, W J Porter +1 more

Abstract Motivated by the realization of hard-wall boundary conditions in experiments with ultracold atoms, we investigate the ground-state properties of spin-1/2 fermions with attractive interactions in a one-dimensional box. We use lattice Monte Carlo methods to determine essential quantities like the energy, which we compute as a function of coupling strength and particle number in the regime from few to many particles. Many-fermion systems bound by hard walls display non-trivial density profiles characterized by so-called Friedel oscillations (which are similar to those observed in harmonic traps). In non-interacting systems, the characteristic length scale of the oscillations is set by ( 2 k F ) − 1 , where k F is the Fermi momentum, while repulsive interactions tend to generate Wigner-crystal oscillations of period ( 4 k F ) − 1 . Based on the non-interacting result, we find a remarkably simple parametrization of the density profiles of the attractively interacting case, which we generalize to the one-body density matrix. While the total momentum is not a conserved quantity in the presence of hard walls, the magnitude of the momentum does provide a good quantum number. We are therefore able to provide a detailed characterization of the (quasi-)momentum distribution, which displays rather robust discontinuity at the Fermi surface. In addition, we determine the spatially varying on-site density–density correlation, which in turn yields Tan’s contact density and, upon integration, Tan’s contact. As is well known, the latter fully determines the short-range correlations and plays a crucial role in a multitude of equilibrium and non-equilibrium sum rules.

arXiv:1510.08503•Journal of Physics B Atomic Molecular and Optical Physics•Oct 24, 2016•11 citations

Harmonically trapped fermions in one dimension: A finite-temperature lattice Monte Carlo study

Felipe Attanasio, Marc Bauer, Renzo Kapust +1 more

We study a one-dimensional two-component Fermi gas in a harmonic trapping potential using finite-temperature lattice quantum Monte Carlo methods. We are able to compute observables in the canonical ensemble via an efficient projective approach. Results for density profiles, correlations, as well as energy-related observables are presented for systems with up to 80 particles and various temperatures. Our simulations reproduce known numerical results and compare well against available experimental data close to the ground state, while at higher temperature they are benchmarked against the exact solution of the two-particle system. This provides an indication that a standard lattice discretization is sufficient to capture the physics of the trapped system. In the special case of a spin-imbalanced gas, we find no sign problem in the studied parameter ranges, allowing access without the need of specialized methods. This includes simulations close to the ground state and at large population imbalance, where we present results for density correlations, indicating pairing at finite total momentum.

arXiv:2309.06766•Physical Review A•Mar 5, 2024•2 citations

Spin-density-functional theory for imbalanced interacting Fermi gases in highly elongated harmonic traps

Gao Xianlong, Reza Asgari

We numerically study imbalanced two component Fermi gases with attractive interactions in highly elongated harmonic traps. An accurate parametrization formula for the ground state energy is presented for a spin-polarized attractive Gaudin-Yang model. Our studies are based on an accurate microscopic spin-density-functional theory through the Kohn-Sham scheme which employs the one-dimensional homogeneous Gaudin-Yang model with a Luther-Emery-liquid ground-state correlation as a reference system. A Thomas-Fermi approximation is examined incorporating the exchange-correlation interaction. By studying the charge and spin density profiles of the system based on these methods, we gain a quantitative understanding of the role of attractive interactions and polarization on the formation of the two-shell structure, with the coexisted Fulde-Ferrell-Larkin-Ovchinnikov-type phase in the center of the trap and either the BCS superfluid phase or the normal phase at the edges of the trap. Our results are in good agreement with the recent theoretical consequences.

arXiv:0801.3060•Physical Review A•Mar 5, 2008•32 citations

Trapped imbalanced fermionic superfluids in one dimension: A variational approach

Kelly R. Patton, Dominique Gautreau, Stephen Kudla +1 more

We propose and analyze a variational wave function for a population-imbalanced one-dimensional Fermi gas that allows for Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) type pairing correlations among the two fermion species, while also accounting for the harmonic confining potential. In the strongly interacting regime, we find large spatial oscillations of the order parameter, indicative of an FFLO state. The obtained density profiles versus imbalance are consistent with recent experimental results as well as with theoretical calculations based on combining Bethe ansatz with the local density approximation. Our variational wave function displays no signature of the FFLO state in the densities of the two fermion species. Nonetheless, the oscillations of the order parameter appear in density-density correlations, both in situ and after free expansion. Furthermore, above a critical polarization, the value of which depends on the interaction, we find the unpaired Fermi-gas state to be energetically more favorable.

arXiv:1702.02933•Physical Review A•Jun 28, 2017•3 citations

Unconventional pairing in few-fermion systems at finite temperature

Daniel Pęcak, Tomasz Sowiński

Attractively interacting two-component mixtures of fermionic particles confined in a one-dimensional harmonic trap are investigated. Properties of balanced and imbalanced systems are systematically explored with the exact diagonalization approach, focusing on the finite-temperature effects. Using single- and two-particle density distributions, specific non-classical pairing correlations are analyzed in terms of the noise correlations-quantity directly accessible in state-of-the-art experiments with ultra-cold atoms. It is shown that along with increasing temperature, any imbalanced system hosting Fulde-Ferrel-Larkin-Ovchinnikov pairs crossovers to a standard Bardeen-Cooper-Schrieffer one characterized by zero net momentum of resulting pairs. By performing calculations for systems with different imbalances, the approximate boundary between the two phases on a phase diagram is determined.

arXiv:2202.07639•Scientific Reports•Oct 19, 2022•6 citations
Using exact continuous quantum Monte Carlo techniques, we study the zero- and finite-temperature properties of a system of harmonically trapped one-dimensional spin-12\frac{1}{2}21​ fermions with short-range interactions. Motivated by experimental searches for modulated Fulde-Ferrel-Larkin-Ovchinikov states, we systematically examine the impact of a spin imbalance on the density profiles. We quantify the accuracy of the Thomas-Fermi approximation, finding that for sufficiently large particle numbers (N\ensuremath≳100)(N\ensuremath{\gtrsim}100)(N\ensuremath≳100) it quantitatively reproduces most features of the exact density profile. The Thomas-Fermi approximation fails to capture small Friedel-like spin and density oscillations and overestimates the size of the fully paired region in the outer shell of the trap. Based on our results, we suggest a range of experimentally tunable parameters to maximize the visibility of the double-shell structure of the system and the Fulde-Ferrel-Larkin-Ovchinikov state in the one-dimensional harmonic trap. Furthermore, we analyze the fingerprints of the attractive contact interactions in the features of the momentum and pair momentum distributions.
Motivated by the realization of hard-wall boundary conditions in experiments with ultracold atoms, we investigate the ground-state properties of spin-1/2 fermions with attractive interactions in a one-dimensional box.We use lattice Monte Carlo methods to determine essential quantities like the energy, which we compute as a function of coupling strength and particle number in the regime from few to many particles.Many-fermion systems bound by hard walls display non-trivial density profiles characterized by so-called Friedel oscillations (which are similar to those observed in harmonic traps).In non-interacting systems, the characteristic length scale of the oscillations is set by (2k F ) -1 , where k F is the Fermi momentum, while repulsive interactions tend to generate Wignercrystal oscillations of period (4k F ) -1 .Based on the non-interacting result, we find a remarkably simple parametrization of the density profiles of the attractively interacting case, which we generalize to the one-body density matrix.While the total momentum is not a conserved quantity in the presence of hard walls, the magnitude of the momentum does provide a good quantum number.We are therefore able to provide a detailed characterization of the (quasi-)momentum distribution, which displays rather robust discontinuity at the Fermi surface.In addition, we determine the spatially varying on-site density-density correlation, which in turn yields Tan's contact density and, upon integration, Tan's contact.As is well known, the latter fully determines the short-range correlations and plays a crucial role in a multitude of equilibrium and non-equilibrium sum rules.
Abstract Motivated by the realization of hard-wall boundary conditions in experiments with ultracold atoms, we investigate the ground-state properties of spin-1/2 fermions with attractive interactions in a one-dimensional box. We use lattice Monte Carlo methods to determine essential quantities like the energy, which we compute as a function of coupling strength and particle number in the regime from few to many particles. Many-fermion systems bound by hard walls display non-trivial density profiles characterized by so-called Friedel oscillations (which are similar to those observed in harmonic traps). In non-interacting systems, the characteristic length scale of the oscillations is set by ( 2 k F ) − 1 , where k F is the Fermi momentum, while repulsive interactions tend to generate Wigner-crystal oscillations of period ( 4 k F ) − 1 . Based on the non-interacting result, we find a remarkably simple parametrization of the density profiles of the attractively interacting case, which we generalize to the one-body density matrix. While the total momentum is not a conserved quantity in the presence of hard walls, the magnitude of the momentum does provide a good quantum number. We are therefore able to provide a detailed characterization of the (quasi-)momentum distribution, which displays rather robust discontinuity at the Fermi surface. In addition, we determine the spatially varying on-site density–density correlation, which in turn yields Tan’s contact density and, upon integration, Tan’s contact. As is well known, the latter fully determines the short-range correlations and plays a crucial role in a multitude of equilibrium and non-equilibrium sum rules.
We study a one-dimensional two-component Fermi gas in a harmonic trapping potential using finite-temperature lattice quantum Monte Carlo methods. We are able to compute observables in the canonical ensemble via an efficient projective approach. Results for density profiles, correlations, as well as energy-related observables are presented for systems with up to 80 particles and various temperatures. Our simulations reproduce known numerical results and compare well against available experimental data close to the ground state, while at higher temperature they are benchmarked against the exact solution of the two-particle system. This provides an indication that a standard lattice discretization is sufficient to capture the physics of the trapped system. In the special case of a spin-imbalanced gas, we find no sign problem in the studied parameter ranges, allowing access without the need of specialized methods. This includes simulations close to the ground state and at large population imbalance, where we present results for density correlations, indicating pairing at finite total momentum.
We numerically study imbalanced two component Fermi gases with attractive interactions in highly elongated harmonic traps. An accurate parametrization formula for the ground state energy is presented for a spin-polarized attractive Gaudin-Yang model. Our studies are based on an accurate microscopic spin-density-functional theory through the Kohn-Sham scheme which employs the one-dimensional homogeneous Gaudin-Yang model with a Luther-Emery-liquid ground-state correlation as a reference system. A Thomas-Fermi approximation is examined incorporating the exchange-correlation interaction. By studying the charge and spin density profiles of the system based on these methods, we gain a quantitative understanding of the role of attractive interactions and polarization on the formation of the two-shell structure, with the coexisted Fulde-Ferrell-Larkin-Ovchinnikov-type phase in the center of the trap and either the BCS superfluid phase or the normal phase at the edges of the trap. Our results are in good agreement with the recent theoretical consequences.
We propose and analyze a variational wave function for a population-imbalanced one-dimensional Fermi gas that allows for Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) type pairing correlations among the two fermion species, while also accounting for the harmonic confining potential. In the strongly interacting regime, we find large spatial oscillations of the order parameter, indicative of an FFLO state. The obtained density profiles versus imbalance are consistent with recent experimental results as well as with theoretical calculations based on combining Bethe ansatz with the local density approximation. Our variational wave function displays no signature of the FFLO state in the densities of the two fermion species. Nonetheless, the oscillations of the order parameter appear in density-density correlations, both in situ and after free expansion. Furthermore, above a critical polarization, the value of which depends on the interaction, we find the unpaired Fermi-gas state to be energetically more favorable.
Attractively interacting two-component mixtures of fermionic particles confined in a one-dimensional harmonic trap are investigated. Properties of balanced and imbalanced systems are systematically explored with the exact diagonalization approach, focusing on the finite-temperature effects. Using single- and two-particle density distributions, specific non-classical pairing correlations are analyzed in terms of the noise correlations-quantity directly accessible in state-of-the-art experiments with ultra-cold atoms. It is shown that along with increasing temperature, any imbalanced system hosting Fulde-Ferrel-Larkin-Ovchinnikov pairs crossovers to a standard Bardeen-Cooper-Schrieffer one characterized by zero net momentum of resulting pairs. By performing calculations for systems with different imbalances, the approximate boundary between the two phases on a phase diagram is determined.