BCS-BEC crossover in bilayers of cold fermionic polar molecules

We investigate the quantum and thermal phase diagram of fermionic polar molecules loaded in a bilayer trapping potential with perpendicular dipole moment. We use both a BCS-theory approach that is most reliable at weak coupling and a strong-coupling approach that considers the two-body bound dimer states with one molecule in each layer as the relevant degree of freedom. The system ground state is a Bose-Einstein condensate (BEC) of dimer bound states in the low-density limit and a paired superfluid (BCS) state in the high-density limit. At zero temperature, the intralayer repulsion is found to broaden the regime of BCS-BEC crossover and can potentially induce system collapse through the softening of roton excitations. The BCS theory and the strongly coupled dimer picture yield similar predictions for the parameters of the crossover regime. The Berezinskii-Kosterlitz-Thouless transition temperature of the dimer superfluid is also calculated. The crossover can be driven by many-body effects and is strongly affected by the intralayer interaction which was ignored in previous studies

Few-body bound state stability of dipolar molecules in two dimensions

Bound structures among dipolar molecules in multilayers are a topic of great interest in the light of recent experiments that have demonstrated the feasibility of the setup. While it is known that two molecules in two adjacent layers will always bind, larger complexes have only been scarcely addressed thus far. Here we prove rigorously that three- and four-body states will never be bound when the dipoles are oriented perpendicular to the layers. The technique employed is general and can be used for more molecules/layers and other geometries. Our analytical findings are supported by numerical calculations for both fermionic and bosonic molecules. Furthermore, we calculate the reduction in intralayer repulsion necessary to bind large complexes and estimate the influence of bound complexes in systems with many layers.Comment: 5 pages, 4 figures, final versio

Borromean ground state of fermions in two dimensions

The study of quantum mechanical bound states is as old as quantum theory itself. Yet, it took many years to realize that three-body borromean systems that are bound when any two-body subsystem is unbound are abundant in nature. Here we demonstrate the existence of borromean systems of spin-polarized (spinless) identical fermions in two spatial dimensions. The ground state with zero orbital (planar) angular momentum exists in a borromean window between critical two- and three-body strengths. The doubly degenerate first excited states of angular momentum one appears only very close to the two-body threshold. They are the lowest in a possible sequence of so-called super-Efimov states. While the observation of the super-Efimov scaling could be very difficult, the borromean ground state should be observable in cold atomic gases and could be the basis for producing a quantum gas of three-body states in two dimensions.Comment: 9 pages, 3 figures, published versio

Three-body recombination at finite energy within an optical model

We investigate three-boson recombination of equal mass systems as function of (negative) scattering length, mass, finite energy, and finite temperature. An optical model with an imaginary potential at short distance reproduces experimental recombination data and allows us to provide a simple parametrization of the recombination rate as function of scattering length and energy. Using the two-body van der Waals length as unit we find that the imaginary potential range and also the potential depth agree to within thirty percent for Lithium and Cesium atoms. As opposed to recent studies suggesting universality of the threshold for bound state formation, our results suggest that the recombination process itself could have universal features.Comment: 5 pages, 5 figure

Shell-Model Monte Carlo Simulations of BCS-BEC Crossover in Few-Fermion Systems

We study a trapped system of fermions with a zero-range two-body interaction using the shell-model Monte Carlo method, providing {\em ab initio} results for the low particle number limit where mean-field theory is not applicable. We present results for the $N$-body energies as function of interaction strength, particle number, and temperature. The subtle question of renormalization in a finite model space is addressed and the convergence of our method and its applicability across the BCS-BEC crossover is discussed. Our findings indicate that very good quantitative results can be obtained on the BCS side, whereas at unitarity and in the BEC regime the convergence is less clear. Comparison to N=2 analytics at zero and finite temperature, and to other calculations in the literature for $N>2$ show very good agreement.Comment: 6 pages, 5 figures, Revtex4, final versio

Analytic Harmonic Approach to the N-body problem

We consider an analytic way to make the interacting N-body problem tractable by using harmonic oscillators in place of the relevant two-body interactions. The two body terms of the N-body Hamiltonian are approximated by considering the energy spectrum and radius of the relevant two-body problem which gives frequency, center position, and zero point energy of the corresponding harmonic oscillator. Adding external harmonic one-body terms, we proceed to solve the full quantum mechanical N-body problem analytically for arbitrary masses. Energy eigenvalues, eigenmodes, and correlation functions like density matrices can then be computed analytically. As a first application of our formalism, we consider the N-boson problem in two- and three dimensions where we fit the two-body interactions to agree with the well-known zero-range model for two particles in a harmonic trap. Subsequently, condensate fractions, spectra, radii, and eigenmodes are discussed as function of dimension, boson number N, and scattering length obtained in the zero-range model. We find that energies, radii, and condensate fraction increase with scattering length as well as boson number, while radii decrease with increasing boson number. Our formalism is completely general and can also be applied to fermions, Bose-Fermi mixtures, and to more exotic geometries.Comment: 30 pages, 12 figures, updated reference