No Arabic abstract
We characterize the high-temperature thermodynamics of rotating bosons and fermions in two- (2D) and three-dimensional (3D) isotropic harmonic trapping potentials. We begin by calculating analytically the conventional virial coefficients $b_n$ for all $n$ in the noninteracting case, as functions of the trapping and rotational frequencies. We also report on the virial coefficients for the angular momentum and associated moment of inertia. Using the bn coefficients, we analyze the deconfined limit (in which the angular frequency matches the trapping frequency) and derive explicitly the limiting form of the partition function, showing from the thermodynamic standpoint how both the 2D and 3D cases become effectively homogeneous 2D systems. To tackle the virial coefficients in the presence of weak interactions, we implement a coarse temporal lattice approximation and obtain virial coefficients up to third order.
By generalizing our automated algebra approach from homogeneous space to harmonically trapped systems, we have calculated the fourth- and fifth-order virial coefficients of universal spin-1/2 fermions in the unitary limit, confined in an isotropic harmonic potential. We present results for said coefficients as a function of trapping frequency (or, equivalently, temperature), which compare favorably with previous Monte Carlo calculations (available only at fourth order) as well as with our previous estimates in the untrapped limit (high temperature, low frequency). We use our estimates of the virial expansion, together with resummation techniques, to calculate the compressibility and spin susceptibility.
The virial expansion characterizes the high-temperature approach to the quantum-classical crossover in any quantum many-body system. Here, we calculate the virial coefficients up to the fifth-order of Fermi gases in 1D, 2D, and 3D, with attractive contact interactions, as relevant for a variety of applications in atomic and nuclear physics. To that end, we discretize the imaginary-time direction and calculate the relevant canonical partition functions. In coarse discretizations, we obtain analytic results featuring relationships between the interaction-induced changes $Delta b_3$, $Delta b_4$, and $Delta b_5$ as functions of $Delta b_2$, the latter being exactly known in many cases by virtue of the Beth-Uhlenbeck formula. Using automated-algebra methods, we push our calculations to progressively finer discretizations and extrapolate to the continuous-time limit. We find excellent agreement for $Delta b_3$ with previous calculations in all dimensions and we formulate predictions for $Delta b_4$ and $Delta b_5$ in 1D and 2D. We also provide, for a range of couplings,the subspace contributions $Delta b_{31}$, $Delta b_{22}$, $Delta b_{41}$, and $Delta b_{32}$, which determine the equation of state and static response of polarized systems at high temperature. As a performance check, we compare the density equation of state and Tan contact with quantum Monte Carlo calculations, diagrammatic approaches, and experimental data where available. Finally, we apply Pade and Pade-Borel resummation methods to extend the usefulness of the virial coefficients to approach and in some cases go beyond the unit-fugacity point.
The relationship between 2D $SO(2,1)$ conformal anomalies in nonrelativistic systems and the virial expansion is explored using recently developed path-integral methods. In the process, the Beth-Uhlenbeck formula for the shift of the second virial coefficient $delta b_2$ is obtained, as well as a virial expansion for the Tan contact. A possible extension of these techniques for higher orders in the virial expansion is discussed.
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 overcome the sign problem, 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 study the non-equilibrium dynamics of cold atoms held in an optical lattice potential. The expansion of an initially confined atom cloud occurs in two phases: an initial quadratic expansion followed by a ballistic behaviour at long times. Accounting for this gives a good description of recent experimental results, and provides a robust method to extract the effective intersite tunneling from time-of-flight measurements.