BoseFermi Pair Correlations in Attractively Interacting BoseFermi Atomic Mixtures
Abstract
We study static properties of attractively interacting BoseFermi mixtures of uniform atomic gases at zero temperature. Using Green’s function formalism we calculate bosonfermion scattering amplitude and fermion selfenergy in the medium to lowest order of the hole line expansion. We study ground state energy and pressure as functions of the scattering length for a few values of the bosonfermion mass ratio and the number ratio . We find that the attractive contribution to energy is greatly enhanced for small values of the mass ratio. We study the role of the BoseFermi pair correlations in the mixture by calculating the pole of the bosonfermion scattering amplitude in the medium. The pole shows a standard quasiparticle dispersion for a BoseFermi pair, for . For small values of the mass ratio, on the other hand, a BoseFermi pair with a finite centerofmass momentum experiences a strong attraction, implying large medium effects. In addition, we also study the fermion dispersion relation. We find two dispersion branches with the possibility of the avoided crossings. This strongly depends on the number rario .
pacs:
PACS number: 03.75.Hh, 05.30.FkI Introduction
Recent developments in the field of cold atomic gases have proven that this system provides an ideal laboratory for the studies of quantum manybody systemspethick (). This is due to the experimental facilities which allow to control various parameters characterizing the manybody system, e.g., external potentials including optical lattices, choice of atoms obeying Bose or Fermi statistics and their mixtures, variable particle densities, etc. The use of Feshbach resonances, in particular, makes the atomic gases an extremely flexible system as it provides a means to control atomic interactionsfesh_th (); fesh_exp (). One thus was, for instance, able to study the BECBCS crossover process in the twocomponent Fermi system, which has been under intense investigation for decadescrossover (). By changing the resonance energies through the external magnetic field, one can in principle change the magnitude and the sign of the scattering length of the interacting particles, keeping track all the way from a resonating fermion pair to a bound composite particle, a bosonic molecule.
The aim of the present paper is the study of pair correlations in a different system, a BoseFermi(BF) mixture. Degenerate mixtures of bosons and fermions have been created since several years, and studies of static and dynamic properties have been performedBFexp (). Among those are the studies of attractively interacting BF systems, where one finds a sudden loss of fermions as the BF attractive interaction is effectively increasedBFcollapse (). Detailed studies of the dynamics of this system are still missing, however. Recently the finding of Feshbach resonances and formation of the bosonfermion molecules have been reportedBFres (); BFmol (). It is thus expected that by controling the BF interaction one may realize an analog of the process found in twocomponent Fermi systems. What should be expected if one replaced fermion pairs in the BECBCS crossover process by BF pairs? Such studies have indeed been performed theoreticallyschuck () (see also bf_pair ().) By adopting a Cooper type twoparticle problem on top of the bosonfermion degenerate system, it was shown that a stable correlated BF pair is created even before the threshold for the BF bound state. In contrast to the BCS case, however, the system allows only one correlated BF pair with a given centerofmass (CM) momentum because of the fermionic nature of the composite particle. It is then suggested that by increasing the BF attractive interaction, one may create BF pairs with different CM momentum stepwise, until finally a new Fermi sea of the BF pairs is completed.
In Ref.schuck () a separable BF interaction has been adopted to elucidate the mechanism of the creation of BF pairs. In the present paper we adopt a standard pseudopotential for the interaction, and calculate energy and pressure of the system for various values of input parameters. We use Green’s function formalism for this system and calculate perturbatively relevant diagrams to lowest order of the holeline expansion. Such formalism has been developed in albus () together with the calculation of the energies including BoseBose (BB) interaction. Our formulation is similar to albus (), but we use the renormalization procedure of randeria () in relating the pseudopotential strength to the wave scattering length. This allows us to formally take the limit , the unitarity limit unit_th (), which is necessary when one considers a (nearly) bound state of a pair of atoms. We then calculate the poles of the BF pair scattering amplitude in the BF medium, which may be compared with the results of schuck (). Studies of the behavior of the poles as a function of input parameters give us suggestions on the role of the BF pair correlations in the static properties of the system.
The content of the paper is as follows: In the next section we present our model based on the Hamiltonian without BoseBose interaction. We calculate the BF scattering amplitude in the BF mixture in ladder approximation, and give formulas for physical quantities in terms of the amplitude. In section 3 we show numerical results for the ground state energy and pressure for various choices of the boson/fermion masses and the values of the BoseFermi interaction. We then study BoseFermi pair correlation in Section 4 by focusing on the pole structure of the bosonfermion scattering amplitude in the mixture. We also calculate the pole of the single fermion Green’s function and study the role of the BoseFermi pair and its dispersion. We summarize our results in section 5 together with a comment on the effects of the BoseBose interaction. Detailed expressions for the scattering amplitude are given in the appendix.
Ii Formulation
We consider a uniform system of a polarized BoseFermi mixture of atomic gases with attractive bosonfermion interaction. The model Hamiltonian of the system is given by
(1) 
where and are the boson and fermion field operators, respectively, and denote bosonic and fermionic kinetic energies, while denotes bosonfermion interaction with strength of the bosonfermion pseudopotential. Effects of the bosonboson interaction will be mentioned later, while the fermionfermion interaction is omitted throughout as we consider onecomponent (polarized) fermions. We will adopt the Bogoliubov approximation in treating the BoseEinstein condensate (BEC), and therefore include in the bosonic chemical potential .
ii.1 Green’s function formalism in the BoseFermi mixture
To treat condensed bosons, we adopt the conventional Bogoliubov method by separating the zero momentum mode from the remainder :
(2) 
together with its conjugate. is the number density of bosons with momentum . As usual we omit the fluctuation of the boson number in the condensate. The boson number operator writes
(3) 
and the Hamiltonian takes the form
(4) 
where
(5) 
and
(6) 
Physical quantities can be expressed in terms of Green’s functions. We define the boson and fermion Green’s functions by
(7)  
(8) 
where are the field operators in the Heisenberg picture, and represents the interacting ground state.
The energy of the system can be expressed in terms of the Green’s functions. The fermion and boson kinetic energies are calculated according to the standard procedure fetter () as:
(9)  
(10) 
where ’s are the Fourier transform of the Green’s functions, is a positive infinitesimal, and we set . The different signs in the two expressions come from the ordering of the field operators.
To calculate the interaction energy, we first consider the Heisenberg equation of motion for the fermion field:
(11) 
Multiplying by and integrating over , we obtain
(12) 
We now introduce fermion and boson selfenergies and through
(13) 
and
(14) 
In the integrand of Eq.(12) one may use the relation from Eq.(13)
(15) 
and finds
(16) 
the first term in the r.h.s. in Eq.(15) giving a null contribution to the integral. The total energy of the system is finally obtained as
(17) 
The thermodynamic potential at zero temperature is given by
(18) 
where
(19) 
The system is characterized by the boson and fermion particle numbers, and . The number of bosons satisfies the thermodynamic relation
(20) 
The parameter should be chosen to minimize the thermodynamic potential
(21) 
which leads to an explicit expression for as shown below.
We also will calculate the pressure to discuss the stability of the system: As usual, it is obtained from the thermodynamic relation
(22) 
ii.2 Selfenergy in the ladder approximation
To obtain the total energy of the system we calculate the fermion selfenergy in ladder approximation. Here the selfenergy is expressed in terms of the twoparticle scattering amplitude, , in the medium of the BoseFermi mixture as shown in fig.1. The interaction energy is accordingly calculated up to the lowest twoparticle correlation diagram, fig.2, in the spirit of the holeline expansion fetter ().
The scattering amplitude in the present model obeys the integral equation albus (); schuck (); fetter (); galli (),
(23) 
where denotes a fourmomentum of the centerofmass motion of the interacting particles, while and are the relative threemomentum in the final and intial states. The boson and fermion free Green’s functions in medium are given by
(24)  
(25) 
where the Fermi momentum is fixed by the fermion density . After the integration over , Eq.(23) becomes
(26) 
with and . We dropped the hole propagation part in accordance with the present approximation. With respect to the Tmatrix equation in schuck (), we notice that there the phase space factor in (26) is replaced by . This is natural, because in schuck () the shift operation (2) for the bosons has not been performed and therefore the free boson occupancy appears additionally. The two formulations are, however, essentially equivalent. From the structure of Eq.(26), one easily finds that depends only on the variable , and we hereafter write simply . One also finds that the integral in Eq.(26) requires a momentum cutoff, which originates from the use of the zerorange interaction. We can remedy this shortcoming by employing the observable Swave scattering length , instead of the pseudopotential coupling constant . We perform this renormalization following the procedure adopted in randeria () (see also, ohashi ()), slightly different from the one in albus (). The Swave scattering length is related to the twoparticle scattering amplitude in vacuum by the relation
(27) 
where is the reduced mass, and obeys the equation similar to Eq.(23) with replaced with the free Green’s function in vacuum. By solving the equation for one obtains
(28) 
where the integral in the denominator involves again the implicit momentum cutoff. Now one may combine the above expression with Eq.(26), and eliminate in favor of the scattering length , and finally obtains
(29) 
with
(30) 
Since the integral in the denominator is convergent at large , we can let the momentum cutoff go to infinity. The expression (29) involves all orders in the scattering length and allows us to formally take the unitarity limit in the following section. This limit has been studied for twocomponent Fermi systems in relation with the BEC to BCS crossover phenomenon. If a similar phenomenon is expected or not for BoseFermi pairs will be studied in the next section.
Using above vertex function, we can calculate the proper selfenergies for the fermion and the boson as
(31)  
(32) 
Expression (26) implies that is analytic in the upper half plane, and Eq.(32) reduces to
(33) 
with . This shows that , and hence, also is analytic in the upper half plane. One then finds from Eq.(20) that
(34) 
Iii Results for Total Energy and Pressure
We calculate the energy of the system in the leading order of the holeline expansion, that is we replace Green’s functions in Eq.(17) with the free one in Eq.(25), and obtain
(35) 
where, and is the Fermi energy. Within the same approximations, the thermodynamic potential at zero temperature is given by
(36) 
Thus, the equilibrium condition (21) for leads to the integral equation for ,
(37) 
where depends also on . The total energy of the system is then finally given by
(38) 
Details of the calculation and the analytic expression for are given in the appendix.
We may rewrite Eq.(37) in a scaled form as
(39) 
where we introduced tilde (dimensionless) quantities through , , , . The expression shows that the scaled chemical potential depends only on the mass ratio and the dimensionless scattering length . We solved Eq.(39) for numerically as a function of the bosonfermion mass ratio for different values of the interaction strength represented by . In terms of the scaled quantities, the ground state energy per particle is expressed from Eq.(38) as
(40) 
where the dimensionless parameter is given by
(41) 
We first show the results for energy and pressure in the unitarity limit, . In this limit, assuming Swave scattering and neglecting effective range, we are left with only one length scale, , or in terms of the density unit_th (); univ_hyp () for a given mass ratio . Note that the chemical potential has no dependence in the unitarity limit, and the parameter depends only on the number ratio . Thus the ground state energy per particle, Eq.(40), is proportional to , and the dependence on the parameters are all absorbed in a simple multiplicative factor .
We show in fig.3 the chemical potential and in fig.4 the beta parameter as functions of the mass ratio , both in the unitarity limit. Note that the results are independent of the magnitude of the individual mass parameters as we show dimensionless quantities scaled with . Figure 3 shows that the boson chemical potential is always negative. This fact reflects the attractive bosonfermion interaction in the unitarity limit, in accordance with Eq.(28) which implies negative . The behavior of in fig.4 simply follows the one of the chemical potential. The results suggest that the attractive interaction becomes more effective for small values of the mass ratio , and the effect is greatly enhanced as the particle number ratio becomes larger, that is as the number of bosons increases with respect to the fermions. The dependence on may partly be understood by noting that the relative phase space available for the intermediate states in the twobody scattering in the mixture will be larger for small , i.e., for a relatively larger , because of the lower Fermi energy and higher level density.
We next consider the pressure to study the stability of the system. Since the total energy takes a universal form and the parameter has no volume dependence in the unitarity limit, the pressure is simply given by
(42) 
For the pressure becomes negative, and the system collapses. This happens especially for larger values of , where the pressure becomes always negative irrespective of the mass ratio . In actual experiments, e.g., for the KRb mixture which has , the number ratio is typically and the system would collapse in the unitarity limit. This is not in contradiction to recent experimental results BFcollapse ().
Another feature in the unitarity limit seen from fig.3 and fig.4 is that the boson chemical potential and the ground state energy saturate when the mass ratio becomes large. It is natural that the bosonic degree of freedom gets frozen and a universal fermionic description appears in this case, since the fermionic effects on the bosons would become negligible, and the bosons would act as a static external field for fermions. This behavior at large also has been discussed in albus ().
We now study the case with an arbitrary value of the scattering length. First, we show the chemical potential as a function of in fig.5, where the mass ratio and are set to 1. Then, we show the energy and pressure as a function of in figs.6, 7 and 8, for the mass ratio with different values of . For , the pressure is given by
(43) 
where the term dependent on reflects that the chemical potential, and hence the parameter, depends on . From the results on energies we see that the system becomes more attractive at smaller values of the mass ratio as in the unitarity limit, although the effect is not large in this parameter range. We find a strong increase of the attraction as the parameter passes through zero, the unitarity limit, from negative to positive. This is in accord with a naive picture where the positive values of the scattering length imply a newly formed bound state. One should however note that even in the present approximation the effects of the medium modify the twobody scattering amplitude, and a simple picture of independent bound pairs does not hold in general. Turning now to the pressure, a comparison of Eq.(43) with the corresponding expression (42) in the unitarity limit shows that the strong attraction for comes from the large negative values of as well as from the coherence of the two terms in the integrand.
Iv BoseFermi Pair Correlation
Results of the previous section indicate that the strong attraction in the mixture will show up especially for positive , which may eventually lead to a collapse of the system. We now consider another scenario for the attractively interacting mixture, the possibility of a BoseFermi pair formationschuck (); bf_pair (). For this purpose we study in this section the behavior of the pole of the BoseFermi scattering amplitude in the mixture.
In unpolarized (or twocomponent) Fermi systems, an infinitesimal attraction around the Fermi surface leads to formation of Cooper pairs with a centerofmass (CM) momentum , causing a transition to the BCS state. In a BoseFermi mixture, on the other hand, the difference in the momentum distribution of the two particles and the fermionic character of the BoseFermi pair, in particular, predict quite a different scenario for the formation of the pairs in the mixture. It requires consideration of the balance of the kinetic energies of different kinds of particles and the magnitude of the attractive interaction.
iv.1 Preliminary considerations
Let us give a brief picture on the formation of BoseFermi pairs in the mixture. Following the idea of schuck (), we may take a BoseFermi mixture with , focusing only one pair of a boson and a fermion with CM momentum , putting other particles as a free background. The Hamiltonian in this system may be written as
(44) 
where we explicitly write the kinetic energy and the interaction for the two particles, while the background Hamiltonian acts only to impose a Pauliprinciple constraint. We neglect here the effect of the boson chemical potential for simplicity. If the effect of on the pair were negligible, the energy of the BoseFermi pair with CM momentum would be simply
(45) 
which has been called free branch in schuck (), since the boson will remain at in the condensate. When the effect of the interaction becomes important, one may rewrite Eq.(44) as
(46) 
where the interaction is contained in the Hamiltonian of the relative motion. By replacing the latter with its eigenvalue (the effect of the medium may be included here), one obtains a different dispersion curve for the BoseFermi pair, the collective branch, which is given by
(47) 
The calculation for the model with a separable interaction in schuck () shows that the two dispersion curves, Eqs.(45) and (47), coexist except for the region of CM momentum given by a solution of the equation
(48) 
There is a mixture of the two branches around , and the dispersion curve of the pair deviates from and . For a sufficiently attractive interaction, the solution of Eq.(48) satisfies the condition , which occurs for
(49) 
One may thus expect that the system will lower the energy by converting the free BoseFermi particles in the range into collective pairs. (Note that we neglect the interaction of the pairs in this simple argument.) If is so strong as to allow for a bound state of the pair, i.e., , all the free bosons and fermions would be replaced with the bound BF pairs, and the new Fermi sea of the pairs will be formed.
This picture may be modified even in this simple model, however, for a system with many BoseFermi pairs. As the pairs can occupy lowmomentum states having lower energies without changing the total momentum of the system, the formation of the collective pairs may start if the condition
(50) 
is satisfied, even before the condition (49). Similarly, the formation of the Fermi sea of the pair will be completed only when
(51) 
is satisfied.
We note that the collective branch (47) may suffer a Laudau type damping for sufficiently large . This is because the free BoseFermi mixture with Fermi momentum has a continuum of the bosonfermion excitation with momentum which starts at the energy schuck ()
(52) 
This implies that the collective branch remains undamped only when
(53) 
a condition which is approximately realized for the pole of the bosonfermion scattering amplitude as shown below.
iv.2 Pole behavior of the twoparticle scattering amplitude
Now we study the behavior of the pole of . A first study concerns the pole condition. That is
(54) 
For , we show in fig.9 the right hand side of (54) as a function of for .
We see the development of a logarithmic divergency as aproaches . This stems from the fact that above dispersion integral has exactly the same structure as the one encountered in the problem of Cooper for a fermion pair in a Fermisea. We want, however, to point out that the pole corresponds to a composite fermion what has important consequences for the physics. Nevertheless the fact is there that a stable collective BF pair developes for any infinitesimal attraction, i.e, even in the limit quite in analogy to the original Cooper pole fetter ().
Let us now discuss the collective pole contained in . Be the pole of with CM momentum . represents the total energy of a bosonfermion pair corresponding to the collective branch aside from the chemical potential. One may then define a dependent binding energy (including total kinetic energy) measured from the last filled free BoseFermi pair by
(55) 
see fig.10. Positive value of would signal a formation of the BoseFermi pair.
We solved numerically. Here we include chemical potential , and show in units of for in figs.11, 12 and 13. The number ratio is fixed at . Left panel of each figure shows results for several values of . One finds that each line extends up to a maximum value of , the energy of which corresponds to in Eq.(52), i.e., the point where the pole hits the continuum and obtains a finite lifetime due to the nonvanishing imaginary part. In the right panels of these figures, we show the dispersion curves at the threshold values of , where turns from negative to positive for the first time.
We see from these figures that the effective binding increases as the interaction becomes more attractive, and eventually leads to a formation of the BoseFermi pair. The threshold value of is lower for a smaller value of in agreement with the result for the total energy. We find also that is always negative in the unitarity limit, suggesting that the BoseFermi pair may not be formed in this limit.
A peculiar feature seen from the figures is the nonmonotonic dependence of against CM momentum , which is in contrast to the dependence of Eq.(47) in the simple picture. This is particularly apparent for small values of . Missing in the simple model is the dependence of the pair binding energy which should be present due, e.g., to the phase space available for the twobody scattering in medium. The medium effect becomes stronger for larger (positive) values of , especially for small as seen from the figure. It may be mentioned in addition, that the results in the next subsection suggest a strong interplay of the collective and free branches for the Fermion dispersion relation.
iv.3 Fermion dispersion and level crossing
Let us now investigate the pole structure of the single particle Green’s function of (13). The pole condition reads,
(56) 
In order to be consistent with our publication in schuck (), we here neglect the chemical potential in of (29). This means that we replace in the Boson and Fermion propagators by totally free ones with the kinetic energies of boson and fermion and which is the uncorrelated value.
With the definition of in (13) with and the expression of (29), we easily can investigate the poles of for various system parameters. We first consider cases, like in schuck (), where .
We expect two branches: one which corresponds in week coupling to and one which corresponds to the collective branch, i.e. the pole contained in .






In fig.14, we show on the left panel the case of and on the right panel the case of . From top to botom, we have . Here mass ratio . The dotted area is the region where the imaginary part of is different from zero. As in our ealier work schuck (), we see an avoided crossing of the two branches at some finite value of . We also see that the interaction between the two branches becomes stronger as increases.
In fig.15, we show the case and for left panel and for right panel. We see that no crossing feature is visible any longer and the two branches probably become completely hybridised. In a future publication we intend to investigate in a systematic way the nature of the two branches as a function of the system parameters.
V Summary, Discussion and Conclusion
We studied in this paper the static properties of the BoseFermi mixture in the lowest order of the holeline expansion, i.e. in the Tmatrix approximation. The interaction parameter is expressed in terms of the scattering length up to infinite order using the renormalization procedure of ref.randeria (), so as to allow for the calculation around the unitarity limit. The Tmatrix approach is a common approximation often used in the past in Fermi and Bose systems fetter (). It has, however, not been investigated very much in BoseFermi mixtures and, therefore, one has not much experience about its quality in that situation. Recently appeared, however, a study for a one dimensional system pollet (), from where exact Quantum Monte Carlo results are available for comparison. In ref xavier () it is shown that the Tmatrix approach yields quite reasonable results also in the BoseFermi case, even in a one dimensional case which is probably the worst situation possible.
With this background in mind, we first studied energy and pressure as functions of the inverse scattering length for several choices of the mass ratio and the number ratio . The energy of the system becomes strongly attractive as the inverse scattering length changes sign from negative to positive, i.e. around the unitary limit. As one increases the number of bosons with respect to the fermions, arrives a point where the pressure becomes negative, i.e. the system becomes unstable (collapse). The effect is stronger for small values of . This is not in contradiction with experiments BFcollapse ().
Next we studied the possibility of stable BFpairs, as in schuck (). We, indeed, also found in the present model that even for infinitesimal BFattraction a stable BFmode appears, reminiscent of the Cooper pole in a two component Fermi gas, since in both cases its origin stems from the presence of a sharp Fermi surface. However, in the BF case the BFpair is a composite fermion, whereas in the original Cooper problem, one has a composite fermion pair, i.e. a bosonlike cluster. In the latter case, many pairs can well be treated by the usual BCS formalism. On the other hand , the case of many BFpairs still has to be worked out. This shall be done in future work. But we here studied the BF pair formation by observing its binding energy measured from the last filled free BF pair. For some finite values of the attractive interaction, there occurs the formation of stable BF pairs. For the pair shows a standard dispersion of a quasiparticle, while for the pair with finite centerofmass momenta feel stronger attraction. This effect is not clearly seen in the energy or pressure, where the singular effect may have been averaged out.
Up to this point we assumed an ideal case where the BoseFermi interaction is dominant, while the BoseBose interaction was neglected. To compare with a real atomic gas system, the BoseBose interaction cannot be discarded even when the BoseFermi interaction is enhanced, e.g., through Feshbach resonances. We checked the effect of the BoseBose interaction on the energy of the system up to the first order in . We took the BoseBose scattering length in the range , and the number ratio 2.0, 1.0 and 0.5, and repeated the calculation of the parameter and the pressure watanabe (). The calculation shows that the effect of the BoseBose interaction is small within the adopted parameter values: The value at the unitarity limit, for instance, deviates less than 10% for and , and does not change conclusions obtained at .
In summary, we suggest that the energy gain in the BoseFermi mixture at positive values of is related to the formation of the resonant BoseFermi pairs, and that the centerofmass momenta of the pairs are dependent on the ratio due to the statistics of the two kinds of the particles.
Acknowledgements.
We thank X. BarillierPertuisel, K. Suziki, T. Nishimura, T. Maruyama and H. Yabu for useful discussions. We also thank J. Dukelsky and S. Pittel for their general interest and contributions to the subject of BF correlations.*
Appendix A Expression for the Scattering Amplitude
Let us calculate the (29) of the scattering amplitude
(57) 
As we are interested in the real part of the pole of , we hereafter omit in the denominator. Each term in the integrand shows an ultraviolet divergence, and we formally introduce a cutoff which will be taken to infinity later.
(58) 
The divergent term at coming from the first term in the integrand is cancelled out by the second term, and we obtain the finite quantit
(59) 
where
(60) 
with
(61) 
The final expression for is given in terms of by