We present a nonperturbative computation of the equation of state of polarized, attractively interacting, nonrelativistic fermions in one spatial dimension at finite temperature. We show results for the density, spin magnetization, magnetic susceptibility, and Tan's contact. We compare with the second-order virial expansion, a next-to-leading-order lattice perturbation theory calculation, and interpret our results in terms of pairing correlations. Our lattice Monte Carlo calculations implement an imaginary chemical potential difference to avoid the sign problem. The thermodynamic results on the imaginary side are analytically continued to obtain results on the real axis. We focus on an intermediate- to strong-coupling regime, and cover a wide range of temperatures and spin imbalances.
We present a nonperturbative computation of the equation of state of polarized, attractively interacting, nonrelativistic fermions in one spatial dimension at finite temperature. We show results for the density, spin magnetization, magnetic susceptibility, and Tan's contact. We compare with the second-order virial expansion, a next-to-leading-order lattice perturbation theory calculation, and interpret our results in terms of pairing correlations. Our lattice Monte Carlo calculations implement an imaginary chemical potential difference to avoid the sign problem. The thermodynamic results on the imaginary side are analytically continued to obtain results on the real axis. We focus on an intermediate- to strong-coupling regime, and cover a wide range of temperatures and spin imbalances.
We calculate the finite-temperature density and polarization equations of state of one-dimensional fermions with a zero-range interaction, considering both attractive and repulsive regimes. In the path-integral formulation of the grand-canonical ensemble, a finite chemical potential asymmetry makes these systems intractable for standard Monte Carlo approaches due to the sign problem. Although the latter can be removed in one spatial dimension, we consider the one-dimensional situation in the present work to provide an efficient test for studies of the higher-dimensional counterparts. To that end, we use the complex Langevin approach, which we compare here with other approaches: imaginary-polarization studies, third-order perturbation theory, and the third-order virial expansion. We find very good qualitative and quantitative agreement across all methods in the regimes studied, which supports their validity.
We analyze the pressure and density equations of state of unpolarized nonrelativistic fermions at finite temperature in one spatial dimension with contact interactions. For attractively interacting regimes, we perform a third-order lattice perturbation theory calculation, assess its convergence properties by comparing with hybrid Monte Carlo results (there is no sign problem in this regime), and demonstrate agreement with real Langevin calculations. For repulsive interactions, we present lattice perturbation theory results as well as complex Langevin calculations, with a modified action to prevent uncontrolled excursions in the complex plane. Although perturbation theory is a common tool, our implementation of it is unconventional; we use a Hubbard-Stratonovich transformation to decouple the system and automate the application of Wick's theorem, thus generating the diagrammatic expansion, including symmetry factors, at any desired order. We also present an efficient technique to tackle nested Matsubara frequency sums without relying on contour integration, which is independent of dimension and applies to both relativistic and nonrelativistic systems, as well as all energy-independent interactions. We find exceptional agreement between perturbative and nonperturbative results at weak couplings, and furnish predictions based on complex Langevin at strong couplings. We additionally present perturbative calculations of up to the fifth-order virial coefficient for repulsive and attractive couplings. Both the lattice perturbation theory and complex Langevin formalisms can easily be extended to a variety of situations including polarized systems, bosons, and higher dimension.
We calculate the pressure and density of polarized non-relativistic systems of two-component fermions coupled via a contact interaction at finite temperature. For the unpolarized one-dimensional system with an attractive interaction, we perform a thirdorder lattice perturbation theory calculation and assess its convergence by comparing with hybrid Monte Carlo. In that regime, we also demonstrate agreement with real Langevin. For the repulsive unpolarized one-dimensional system, where there is a so-called complex phase problem, we present lattice perturbation theory as well as complex Langevin calculations. For our studies, we employ a Hubbard-Stratonovich transformation to decouple the interaction and automate the application of Wick's theorem for perturbative calculations, which generates the diagrammatic expansion at any order. We find excellent agreement between the results from our perturbative calculations and stochastic studies in the weakly interacting regime. In addition, we show predictions for the strong coupling regime as well as for the polarized one-dimensional system. Finally, we show a first estimate for the equation of state in three dimensions where we focus on the polarized unitary Fermi gas.
We calculate the pressure and density of polarized non-relativistic systems of two-component fermions coupled via a contact interaction at finite temperature. For the unpolarized one-dimensional system with an attractive interaction, we perform a thirdorder lattice perturbation theory calculation and assess its convergence by comparing with hybrid Monte Carlo. In that regime, we also demonstrate agreement with real Langevin. For the repulsive unpolarized one-dimensional system, where there is a so-called complex phase problem, we present lattice perturbation theory as well as complex Langevin calculations. For our studies, we employ a Hubbard-Stratonovich transformation to decouple the interaction and automate the application of Wick’s theorem for perturbative calculations, which generates the diagrammatic expansion at any order. We find excellent agreement between the results from our perturbative calculations and stochastic studies in the weakly interacting regime. In addition, we show predictions for the strong coupling regime as well as for the polarized one-dimensional system. Finally, we show a first estimate for the equation of state in three dimensions where we focus on the polarized unitary Fermi gas.
We present finite-temperature, lattice Monte Carlo calculations of the particle number density, compressibility, pressure, and Tan's contact of an unpolarized system of short-range, attractively interacting spin-1/2 fermions in one spatial dimension, i.e., the Gaudin-Yang model. In addition, we compute the second-order virial coefficients for the pressure and the contact, both of which are in excellent agreement with the lattice results in the low-fugacity regime. Our calculations yield universal predictions for ultracold atomic systems with broad resonances in highly constrained traps. We cover a wide range of couplings and temperatures and find results that support the existence of a strong-coupling regime in which the thermodynamics of the system is markedly different from the noninteracting case. We compare and contrast our results with identical systems in higher dimensions.
Using ab initio lattice methods, we calculate the finite temperature thermodynamics of homogeneous two-dimensional spin-1/2 fermions with attractive short-range interactions. We present results for the density, pressure, compressibility, and quantum anomaly (i.e., Tan's contact) for a wide range of temperatures (mostly above the superfluid phase, including the pseudogap regime) and coupling strengths, focusing on the unpolarized case. Within our statistical and systematic uncertainties, our prediction for the density equation of state differs quantitatively from the prediction by Luttinger-Ward theory in the strongly coupled region of parameter space, but otherwise agrees well with it. We also compare our calculations with the second- and third-order virial expansion, with which they are in excellent agreement in the low-fugacity regime.