Ground-state properties of trapped Bose-Fermi mixtures: role of exchange-correlation

We introduce Density Functional Theory for inhomogeneous Bose-Fermi mixtures, derive the associated Kohn-Sham equations, and determine the exchange-correlation energy in local density approximation. We solve numerically the Kohn-Sham system and determine the boson and fermion density distributions and the ground-state energy of a trapped, dilute mixture beyond mean-field approximation. The importance of the corrections due to exchange--correlation is discussed by comparison with current experiments; in particular, we investigate the effect of of the repulsive potential energy contribution due to exchange--correlation on the stability of the mixture against collapse.Comment: 6 pages, 4 figures (final version as published in Physical Review

Mean-field analysis of the stability of a K-Rb Fermi-Bose mixture

We compare the experimental stability diagram of a Fermi-Bose mixture of K-40 and Rb-87 atoms with attractive interaction to the predictions of a mean-field theoretical model. We discuss how this comparison can be used to give a better estimate of the interspecies scattering length, which is currently known from collisional measurements with larger uncertainty.Comment: 5 pages, 4 figure

Finite temperature effects on the collapse of trapped Bose-Fermi mixtures

By using the self-consistent Hartree-Fock-Bogoliubov-Popov theory, we present a detailed study of the mean-field stability of spherically trapped Bose-Fermi mixtures at finite temperature. We find that, by increasing the temperature, the critical particle number of bosons (or fermions) and the critical attractive Bose-Fermi scattering length increase, leading to a significant stabilization of the mixture.Comment: 5 pages, 4 figures; minor changes, proof version, to appear in Phys. Rev. A (Nov. 1, 2003

Thermodynamics of a Trapped Bose-Fermi Mixture

By using the Hartree-Fock-Bogoliubov equations within the Popov approximation, we investigate the thermodynamic properties of a dilute binary Bose-Fermi mixture confined in an isotropic harmonic trap. For mixtures with an attractive Bose-Fermi interaction we find a sizable enhancement of the condensate fraction and of the critical temperature of Bose-Einstein condensation with respect to the predictions for a pure interacting Bose gas. Conversely, the influence of the repulsive Bose-Fermi interaction is less pronounced. The possible relevance of our results in current experiments on trapped $^{87}{\rm Rb}-^{40}${\rm K} mixtures is discussed.Comment: 5 pages + 4 figures; minor changes, final version to appear in Phys. Rev. A; the extension work on the finite-temperature low-lying excitations can be found in cond-mat/030763

Collective ferromagnetism in two-component Fermi-degenerate gas trapped in finite potential

Spin asymmetry of the ground states is studied for the trapped spin-degenerate (two-component) gases of the fermionic atoms with the repulsive interaction between different components, and, for large particle number, the asymmetric (collective ferromagnetic) states are shown to be stable because it can be energetically favorable to increase the fermi energy of one component rather than the increase of the interaction energy between up-down components. We formulate the Thomas-Fermi equations and show the algebraic methods to solve them. From the Thomas-Fermi solutions, we find three kinds of ground states in finite system: 1) paramagnetic (spin-symmetric), 2) ferromagnetic (equilibrium) and 3) ferromagnetic (nonequilibrium) states. We show the density profiles and the critical atom numbers for these states obtained analytically, and, in ferromagnetic states, the spin-asymmetries are shown to occur in the central regions of the trapped gas, and grows up with increasing particle number. Based on the obtained results, we discuss the experimental conditions and current difficulties to realize the ferromagnetic states of the trapped atom gas, which should be overcome.Comment: submit to PR

Hydrodynamic excitations of trapped dipolar fermions

A single-component Fermi gas of polarized dipolar particles in a harmonic trap can undergo a mechanical collapse due to the attractive part of the dipole-dipole interaction. This phenomenon can be conveniently manipulated by the shape of the external trapping potential. We investigate the signatures of the instability by studying the spectrum of low-lying collective excitations of the system in the hydrodynamic regime. To this end, we employ a time-dependent variational method as well as exact numerical solutions of the hydrodynamic equations of the system.Comment: 4 pages, 2 eps figures, final versio

Random-phase approximation study of collective excitations in the Bose-Fermi mixed condensate of alkali-metal gases

We perform Random Phase Approximation (RPA) study of collective excitations in the bose-fermi mixed degenerate gas of Alkali-metal atoms at T=0. The calculation is done by diagonalization in a model space composed of particle-hole type excitations from the ground state, the latter being obtained from the coupled Gross-Pitaevskii and Thomas-Fermi equations. We investigate strength distributions for different combinations of bose and fermi multipole ($L$) operators with $L=0,1,2,3$. Transition densities and dynamical structure factors are calculated for collective excitations. Comparison with the sum rule prediction for the collective frequency is given. Time dependent behavior of the system after an external impulse is studied.Comment: 28 pages, 13 figures, submitted to Phys. Rev.

Pairing in two-dimensional boson-fermion mixtures

The possibilities of pairing in two-dimensional boson-fermion mixtures are carefully analyzed. It is shown that the boson-induced attraction between two identical fermions dominates the p-wave pairing at low density. For a given fermion density, the pairing gap becomes maximal at a certain optimal boson concentration. The conditions for observing pairing in current experiments are discussedComment: 10 pages, 5 figs, revtex

Helioseismology of Sunspots: A Case Study of NOAA Region 9787

Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot.Comment: 28 pages, 18 figure

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.Comment: 73 pages, 19 figures, accepted by Solar Physic