In this article we present a Monte Carlo calculation of the critical temperature and other thermodynamic quantities for the unitary Fermi gas with a population imbalance (unequal number of fermions in the two spin components). We describe an improved worm-type algorithm that is less prone to autocorrelations than the previously available methods and show how this algorithm can be applied to simulate the unitary Fermi gas in presence of a small imbalance. Our data indicate that the critical temperature remains almost constant for small imbalances h = Δμ/ε F &z.lap; 0.2. We obtain the continuum result T c = 0.171(5)ε F in units of Fermi energy and derive a lower bound on the deviation of the critical temperature from the balanced limit, T c(h) - T c(0) > -0.5ε F h 2. Using an additional assumption a tighter lower bound can be obtained. We also calculate the energy per particle and the chemical potential in the balanced and imbalanced cases.
|Journal||Physical Review A - Atomic, Molecular, and Optical Physics|
|State||Published - 2010|