Three-dimensional correlated-fermion phase separation from analysis of the geometric mean of the individual susceptibilities
A quasi-Gaussian approximation scheme is formulated to study the strongly correlated imbalanced fermions thermodynamics, where the mean-field theory is not applicable. The non-Gaussian correlation effects are understood to be captured by the statistical geometric mean of the individual susceptibilities. In the three-dimensional unitary fermions ground state, an universal non-linear scaling transformation relates the physical chemical potentials with the individual Fermi kinetic energies. For the partial polarization phase separation to full polarization, the calculated critical polarization ratio is . The defines the ratio of the symmetric ground state energy density to that of the ideal fermion gas.
pacs:05.30.Fk, 03.75.Hh, 21.65.+f
Phase separation resulting from the microscopic dynamics involves a wide range of natural phenomena. Currently, it also serves as the pivotal topic in the strongly interacting fermions quantum many-body theory.
In the past few years, considerable efforts understanding the crossover physics from Bardeen-Cooper-Schrieffer to Bose-Einstein condensation(BCS-BEC crossover) with ultra-cold atomic fermi gases have been made. At the Feshbach resonance point, the divergent -wave scattering length with the existence of a zero-energy bound state can exhibit the thermodynamic universalityHo2004 (); Giorgini2007 ().
Another exciting issue is on the asymmetric fermions phase diagram with unequal populations, which persists as a fundamental theme for a long timeSarma (); Fulde (); Larkin (). With the Feshbach resonance techniques, tuning interaction strength and controlling the population or mass imbalance among the components offer a playground to test the non-perturbative many-body theory for the asymmetric fermions thermodynamicsPartridge2006 (); Zwierlein (). The experiments show that at zero temperature the system can undergo a first-order quantum phase transition from a fully-paired superfluidity to a partially polarized normal gasYong2007 (). Furthermore, when the imbalance population ratio between the two spin components reaches a critical value , the state can transfer from the partially polarized phase separation to the fully polarized normal one.
The three-dimensional strongly interacting asymmetric fermions ground state is one of the prominent unresolved problemsPethick2002 (). Meanwhile, the zero-temperature ground state energy of the partial polarization phase separation is of particular importance for explaining experimentsYong2007 (); Yong2008 (). For the equally populated spin states, the Monte Carlo methods are believed to be capable of doing the simulation calculations. For the asymmetric system, the numerical simulations will suffer from the serious fermion sign problemBraaten2008 (). Making a definitive theoretical conclusion or finding an unbiased analytical solution remains as an intriguing but herculean task. As far as we are aware, until recently, any analytical attempt was not known.
There can exist a relation between the and dimensionless constant at unitarity with Cohen2005 (); reddy (); Bulgac2007 (); Yong2008 (); He2008 (). The defines the ratio of the symmetric fermions energy density to that of the ideal noninteracting ones(). By exploring the complex non-Gaussian correlations, the purpose of this Letter is offering a novel tool to calibrate the ground state energy of strongly interacting asymmetric fermions and fix the relation between and . The method is based on extending the formalism developed inchen2007 (); chen2008 (); chen2007-1 (). The calculations are performed with .
To address the imbalanced fermions thermodynamics, we formulate the universal medium-scaling Hamiltonian
In Eq.(1), represent the hyperfine spin projection Ising-variables. The Hamiltonian is the same as the Bethe-Peierls contact interaction, except that the bare coupling constant is substituted by a medium regulated functional . The is the reduced mass. In the vacuum limit , the reduces to the original version possessing a global or gauge symmetry (assuming without a loss of generality).
Due to the medium dependence of , the non-trivial enforce the energy-momentum conservation lawBrown2002 (). Routinely, the strong correlation effects are further understood as a spontaneously generated one-body correlation potential from the density functional theory viewpointchen2007 (); chen2008 (), with which the general thermodynamic consistencies are strictly guaranteed within the non-perturbative procedures.
The medium-scaling formalism allows a natural implementation of the strong correlation physics. Initially, we need to constitute the functional. The main observation is that the environment frustrating characteristic can be realized by the twisted contact interaction
With Eq.(2), one can define the in-medium inverse scattering length notation according to .
With the “negative” sign in the denominator of Eq.(2), the interaction is alternatively renormalized by the strong correlation effects or the background fluctuations induced in the counterpart Fermi sea. To a great extent, the susceptibility is interpreted as an external field modulator simulating the mutual non-linear response. Economically, it characterizes the complications of the single particle spectrum properties and correlating couplings.
At unitarity, the scaling thermodynamic quantities depend solely on the susceptibilities. To quantify the complicated correlation effects beyond the Gaussian statistics, the global variable of interest is captured by the statistical geometric mean(instead of arithmetic mean) of the individual and for the and subsystems
The susceptibility itself can be separately calculated from the Lindhard correlation response function in terms of the random phase or spin “wave” approximation, which is also the final goalchen2007-1 (). Alternatively, it is more convenient to calculate the susceptibility with the generalized Ward-Identity
The is the Fermi integrals with and the thermodynamical de Broglie wavelength is . The collective dynamical variables are defined by the gap equations Eqs.(7). The effective fugacity is analogous to the fugacity chen2007 (); chen2008 ().
At , is the density of states(DOS) geometric mean
Here, is the familiar un-perturbated DOS near the Fermi surface for the one component fermionsPethick2002 ().
From the action Eq.(1), the grand thermodynamical potential or pressure can be presented as the coupled parametric equations with the instantaneous quasi-Gaussian approximation methodchen2007 (); chen2008 ()
The dynamical gap equations Eqs.(7) for the corresponding single particle Green function give the definition of the effective chemical potential , respectively. The number density is expressed in terms of the defined quasi-particle Fermi-Dirac distribution function
Sticking to the effective interaction Eq.(2) renormalized by the mixing susceptibility Eq.(3) with Eq.(4), the strengths of the correction terms to the Gaussian term are self-consistently derivedchen2007 (); chen2008 ()
A few remarks are re-emphasized on the grand thermodynamical potential because it plays an important role. The is not the naive polynomial expanded according to the bare vacuum interaction strength . The collective dynamical variable mixes the low and high order contributions. Hidden in the factors and as well as in , the collective correlations are combined with the infinite individual dynamical high order effects. The equations Eq.(6)-Eq.(7) are highly non-linear although they appear a set of coupled algebra ones. The dependence of can be numerically eliminated in favor of the physical chemical potential .
In the following, we will examine the zero-temperature asymmetric interacting system thermodynamics, which can be derived analytically. At , the pressure Eq.(6) and energy density Eq.(9) expressions read
where the “Fermi” kinetic energy , particle number density and rearrangement factor are
The reduced chemical potentials Eq.(7) and are
The is the corresponding Fermi momentum.
With unitary , the above expressions can be further reduced to an oversimplified compact formalism
The salient feature of the main results Eq.(13) is that they appear as the component contributions “plus” those of . The zero-temperature unitary fermions thermodynamics obeys the new form of universality, i.e., the Dalton partial pressure law of the non-interacting ideal gas. Implicitly, the interaction and collective correlation information is incorporated through the chemical potentials Eq.(14).
Let us further discuss the physical chemical potentials. They are related with the Fermi kinetic energies of the two subsystems through the non-linear transformations Eq.(14). With Eq.(14), the energy density and pressure can be expressed as the regular functions of the Fermi kinetic energies according to Eq.(13).
With the population imbalance polarization ratio (, assuming and )
the numerical solutions of the chemical potentials are presented in Fig.1. One can see the minority chemical potential() decreases very rapidly while crossing the transverse axis.
For -dimensional, the reciprocal relations of the transformations Eq.(14) are very involved; i.e., it is a hard task to express the Fermi kinetic energies in terms of the chemical potentials and . In principle, the Fermi kinetic energies and particle number densities in terms of the physical chemical potentials can be the functional formalisms such as . In addition to the mathematical complication, the inverse transformation can be singular and nontrivial when . The minority fermions chemical potential will change sign at this point indicated by . In other words, the minority component will quickly “collapse” due to the strongly attractive interaction and the system consists of a single Fermi surface; below it, there can be two distinct Fermi surfaces and the phase separation is favoredSheehy2007 ().
Namely, the vanishing determines the transferring criterion position from the phase separation-partial polarization to the full polarization state. The analytical relation between and can be derived from Eq.(14)
With , the calculated critical ratio is . It is in agreement with the Monte Carlo(MC) result Lobo2006 (); Pilati2008 (). The mean-field theory gives Sheehy2007 (), while the recent experimental value is Yong2007 ().
with the analytical scaling function
It is easy to verify that corresponds to the critical concentration defined in Refs.Yong2008 (); Lobo2006 () with , where the scaling function takes the minimum value as indicated in Fig. 3(b). The convex behavior of coincides with the experimental measurementYong2008 (). The numerical solution of the versus is displayed in Fig.2.b. These results contribute to understanding the realistic trapped system thermodynamics.
To conclude, the complicated non-Gaussian fluctuation and correlation effects beyond the canonical Gaussian techniques are encoded with the statistical geometric mean of the individual susceptibilities. The simple medium-renormalized effective action constitutes the bridge towards fixing the energy density functional of the strong interaction fermions system.
The non-linear scaling transformation identity Eq.(14) relates the physical chemical potentials with the Fermi kinetic energies. The reassuring association between the critical proportional ratio and universal coefficient agrees well with the Monte Carlo calculations and experimental measurements. This non-Gaussian correlation perspective will stimulate further insights in understanding the novel quantum phase separation dynamics.
Acknowledgements.J.-s Chen acknowledges the discussions with L.-y He, D.-f Hou, J.-r Li and X.-j Xia. Supported by the Fund of Central China Normal University, Natural Science Foundation of China under Grant No. 10875050, 10675052 and MOE of China under projects No.IRT0624.
- (1) T.-L. Ho, Phys. Rev. Lett. 92, 090402 (2004).
- (2) S. Giorgini, L.P. Pitaevskii and S. Stringari, Rev. Mod. Phys. 80, 1215 (2008) and references therein.
- (3) G. Sarma, J. Phys. Chem. Solids. 9, 266 (1962).
- (4) A.I. Larkin and Y.N. Ovchinnikovm, Phys. JETP 20, 762 (1965).
- (5) P. Fulde and R.A. Ferrell, Phys. Rev. 135, A550 (1964).
- (6) G.B. Partridge, W. Li, R.I. Kamar, Y.-A Liao, and R.G. Hulet, Science 311, 503 (2006).
- (7) M.W. Zwierlein, A. Schirotzek, C.H. Schunck and W. Ketterle, Science 311, 492 (2006).
- (8) Y.I. Shin, C.H. Schunck, A. Schirotzek and W. Ketterle, Nature 451, 689 (2008).
- (9) C.J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases (Cambridge University Press, Cambridge, England, 2002).
- (10) Y.I. Shin, Phys. Rev. A 77, 041603(R) (2008).
- (11) E. Braaten and L. Platter, Phys. Rev. Lett. 100, 205301 (2008).
- (12) T.D. Cohen, Phys. Rev. Lett. 95, 120403 (2005).
- (13) J. Carlson and S. Reddy, Phys. Rev. Lett. 95, 060401 (2005).
- (14) A. Bulgac and M.M. Forbes, Phys. Rev. A 75, 031605(R) (2007).
- (15) L.-Y He and P.-F Zhuang, e-print arXiv:0803.4052.
- (16) J.-S Chen, C.-M Cheng, J.-R Li, and Y.-P Wang, Phys. Rev. A 76, 033617 (2007); J.-S Chen, J.-R Li, Y.-P Wang and X.-J Xia, J. Stat. Mech. (2008) P12008.
- (17) J.-S Chen, e-print arXiv:0807.4781.
- (18) J.-S Chen, Chin. Phys. Lett. 24, 1825 (2007) and Commun. Theor. Phys. 48, 99 (2007).
- (19) G.E. Brown and M. Rho, Phys. Rep. 363, 85 (2002).
- (20) D.E. Sheehy and L. Radzihovsky, Ann. Phys. (N.Y.) 322, 1790 (2007).
- (21) C. Lobo, A. Recati, S. Giorgini, and S. Stringari, Phys. Rev. Lett. 97, 200403 (2006).
- (22) S. Pilati and S. Giorgini, Phys. Rev. Lett. 100, 030401 (2008).