No Arabic abstract
We present in detail two variants of the lattice Monte Carlo method aimed at tackling systems in external trapping potentials: a uniform-lattice approach with hard-wall boundary conditions, and a non-uniform Gauss-Hermite lattice approach. Using those two methods, we compute the ground-state energy and spatial density profile for systems of N=4 - 8 harmonically trapped fermions in one dimension. From the favorable comparison of both energies and density profiles (particularly in regions of low density), we conclude that the trapping potential is properly resolved by the hard-wall basis. Our work paves the way to higher dimensions and finite temperature analyses, as calculations with the hard-wall basis can be accelerated via fast Fourier transforms, the cost of unaccelerated methods is otherwise prohibitive due to the unfavorable scaling with system size.
The formation of bosonic bound states underlies the formation of a superfluid ground state in the many-body phase diagram of ultracold Fermi gases. We study bound-state formation in a spin- and mass-imbalanced ultracold Fermi gas confined in a box with hard-wall boundary conditions. Because of the presence of finite Fermi spheres, the center-of-mass momentum of the potentially formed bound states can be finite, depending on the parameters controlling mass and spin imbalance as well as the coupling strength. We exploit this observation to estimate the potential location of inhomogeneous phases in the many-body phase diagram as a function of spin- and mass imbalance as well as the box size. Our results suggest that a hard-wall box does not alter substantially the many-body phase diagram calculated in the thermodynamic limit. Therefore, such a box may serve as an ideal trap potential to bring experiment and theory closely together and facilitate the search for exotic inhomogeneous ground states.
We investigate the strongly interacting hard-core anyon gases in a one dimensional harmonic potential at finite temperature by extending thermal Bose-Fermi mapping method to thermal anyon-ferimon mapping method. With thermal anyon-fermion mapping method we obtain the reduced one-body density matrix and therefore the momentum distribution for different statistical parameters and temperatures. At low temperature hard-core anyon gases exhibit the similar properties as those of ground state, which interpolate between Bose-like and Fermi-like continuously with the evolution of statistical properties. At high temperature hard-core anyon gases of different statistical properties display the same reduced one-body density matrix and momentum distribution as those of spin-polarized fermions. The Tans contact of hard-core anyon gas at finite temperature is also evaluated, which take the simple relation with that of Tonks-Girardeau gas $C_b$ as $C=frac12(1-coschipi)C_b$.
In this paper we study a mixed system of bosons and fermions with up to six particles in total. All particles are assumed to have the same mass. The two-body interactions are repulsive and are assumed to have equal strength in both the Bose-Bose and the Fermi-Boson channels. The particles are confined externally by a harmonic oscillator one-body potential. For the case of four particles, two identical fermions and two identical bosons, we focus on the strongly interacting regime and analyze the system using both an analytical approach and DMRG calculations using a discrete version of the underlying continuum Hamiltonian. This provides us with insight into both the ground state and the manifold of excited states that are almost degenerate for large interaction strength. Our results show great variation in the density profiles for bosons and fermions in different states for strongly interacting mixtures. By moving to slightly larger systems, we find that the ground state of balanced mixtures of four to six particles tends to separate bosons and fermions for strong (repulsive) interactions. On the other hand, in imbalanced Bose-Fermi mixtures we find pronounced odd-even effects in systems of five particles. These few-body results suggest that question of phase separation in one-dimensional confined mixtures are very sensitive to system composition, both for the ground state and the excited states.
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 kF)^(-1), where kF is the Fermi momentum, while repulsive interactions tend to generate Wigner-crystal oscillations of period (4 kF)^(-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 Tans contact density and, upon integration, Tans 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.
In the dilute limit, the properties of fermionic lattice models with short-range attractive interactions converge to those of a dilute Fermi gas in continuum space. We investigate this connection using mean-field and we show that the existence of a finite lattice spacing has consequences down to very small densities. In particular we show that the reduced translational invariance associated to the lattice periodicity has a pivotal role in the finite-density corrections to the universal zero-density limit. For a parabolic dispersion with a sharp cut-off, we provide an analytical expression for the leading-order corrections in the whole BCS-BEC crossover. These corrections, which stem only from the unavoidable cut-off, contribute to the leading-order corrections to the relevant observables. In a generic lattice we find a universal power-law behavior $n^{1/3}$ which leads to significant corrections already for small densities. Our results pose strong constraints on lattice extrapolations of dilute Fermi gas properties.