# Trapped Bose gases with large positive scattering length

###### Abstract

We calculate the energy and the condensate fraction of a system of trapped bosons interacting via a short-range two-body potential with positive scattering length. The potential is attractive and has a two-body bound state. When the scattering length is small compared to the trap length the system is model independent: all potential models – attractive, repulsive and zero-range – provide similar results. When the scattering length is large the attractive model differs qualitatively from the repulsive and zero-range models. In this regime the system with attractive potential becomes independent of the scattering length, with both the energy and the condensate fraction converging towards finite constants.

###### pacs:

03.75.Hh, 03.75.Nt, 21.45.+v

## 1 Introduction

The density of trapped cold gases is generally low under typical experimental conditions [1], such that the parameter is small, where is the particle density and is the range of the inter-particle potential. In other words the typical distance between particles is much larger than the range of the potential, and the typical relative momentum between particles is much smaller than the inverse range of the potential.

In this regime the system of particles exhibits universality (also called model independence or shape independence) [2]. The system is not sensitive to the details of the potential and the properties of the system are essentially determined by only few low-energy parameters of the potential. Very different interaction models then provide quantitatively similar results as soon as the low-energy parameters are the same.

In a two-body system the universality is manifested in the well known effective range expansion, where the low-energy -wave phase shift is determined by only two parameters, the scattering length and the effective range ,

(1) |

where is the relative momentum between particles. The effective range is typically of the order of the range of the potential while the scattering length can vary greatly.

In three-body systems the universality manifests itself in the Thomas [3] and Efimov [4] effects, the Phillips line [5], and other low-energy phenomena [6], also in two dimensions [7].

Cold gases also exhibit universality: in the dilute limit their properties, in particular the energy per particle, are independent on the shape of the inter-particle potential and are determined by the scattering length alone. This universality is customarily employed by using the zero-range (pseudo-)potential for theoretical descriptions of the Bose-Einstein condensates. The zero-range potential has only one parameter, the scattering length. Combined with the Hartree-Fock product wave-function the zero-range potential model is known as the Gross-Pitaevskii equation [8].

The scattering length in a trapped atomic gas can be tuned to an essentially arbitrary value using the technique of Feshbach resonances [1]. This gives a possibility to investigate trapped systems with very large scattering lengths. In this case the first term in the effective range expansion eq. (1) vanishes and the system becomes sensitive to the effective range of the potential.

The limits of the zero-range model have been tested by numerical calculations with finite-range potentials. In particular, for the large positive scattering length repulsive potentials have been employed within Monte-Carlo methods [9, 10, 11]. These investigations showed that as the scattering length is increased the energy of the Bose gas with a repulsive potential exceeds the zero-range predictions and the condensate fraction becomes considerably depleted.

However, repulsive potentials have a problem when modelling an increasingly large positive scattering length: the effective range of the potential has to be increased essentially linearly with the scattering length. This does not seem to match the experimental conditions where the scattering length is adjusted by tuning the atomic resonances in an external magnetic field. The range of the inter-atomic interaction is then left essentially unchanged.

Instead, an attractive finite-range potential might be a more realistic interaction model for descriptions of trapped Bose gases with Feshbach resonances. Indeed, with an attractive potential an arbitrary large positive scattering length can be achieved by fine-tuning the energy of the bound two-body state, while maintaining the given realistic effective range.

However, attractive potentials with bound states bring in a major complication for numerical calculations: a large number of many-body self-bound negative-energy states appears in the system and the condensate state in the trap becomes a highly excited state.

For a homogeneous Bose gas an approximate Jastrow-type wave-function was employed where the pair-correlation function was essentially a solution of the two-body equation [12, 10]. In contrast we propose a direct numerical diagonalisation of the many-body Hamiltonian where the condensate state of trapped bosons appears as a many-body excited state which is automatically orthogonal to all the self-bound negative-energy states.

The purpose of this paper is to investigate the energy and the condensate fraction of a system of trapped bosons with attractive potentials as function of the scattering length and the number of bosons, and compare the results with the zero-range and repulsive model.

## 2 System and numerical techniques

### 2.1 The system

We consider a system of identical bosons with mass and coordinates , , in a spherical harmonic trap with frequency . The Hamiltonian of the system is given by

(2) |

where the system parameters are taken from [9]: amu, Hz, the trap length au.

### 2.2 Two-body potentials

The zero-range potential model,

(3) |

has only one length parameter, . For dilute bosonic systems this parameter is customarily chosen to be equal the inter-atomic scattering length. The zero-range potential provides then the correct low-energy scattering amplitude in the first order Born approximation. The zero-range potential can only be used with an appropriate non-correlated functional space [13], like the Hartree-Fock product wave-functions. In the latter case it leads to the famous Gross-Pitaevskii equation [8].

For the finite-range potential model we use a Gaussian,

(4) |

with the range au and a varied negative strength . The variation of the strength is limited to the region where the potential provides exactly one two-body bound state and a positive scattering length.

### 2.3 Stochastic variational method

The wave-function of the system is represented as a linear combination of basis-functions taken in the form of symmetrised correlated Gaussian,

(5) |

where is the symmetrisation operator, and and are variational parameters. The linear parameters are determined by an ordinary diagonalisation of the Hamiltonian eq. (2) while the non-linear parameters are optimised stochastically by random sampling [14, 15] from a region that covers the distances from to . The center-of-mass motion is assumed to be in the oscillator’s ground state.

The zero-range potential eq. (3) requires an uncorrelated wave-function which we chose in the form of the linear combination of the hyper-radial basis-functions,

(6) |

where is the hyper-radius of the system. This function is totally symmetric and thus does not require the symmetrisation operator . The zero-range potential with the hyper-radial variational wave-function eq. (6) provides results similar to Gross-Pitaevskii equation [15].

The calculation of a highly excited state with the fully correlated basis eq. (5) is a difficult numerical task and is only possible for relatively small number of particles.

However for a typical system of trapped atoms even when the scattering length is large the density of the system remains small, , and one can assume that only binary collisions play a significant role in the system’s dynamics. In this approximation the variational wave-function can be simplified by only allowing two-body correlations in the basis-functions,

(7) | |||||

where and are the nonlinear parameters. The symmetrisation of this function can be done analytically [15] which greatly simplifies the numerical calculations.

During the calculation of a given system the number of Gaussian in the basis is increased and the stochastic optimisation is carried out until the number of negative energy states and the energy of the lowest state with positive energy is converged. The convergence within four digits typically requires about Gaussian and about random trials per nonlinear parameter. The stochastic optimisation algorithm is easily parallelisable with close to linear scalability.

### 2.4 Condensate fraction

We first calculate the one-body density matrix defined as

(8) | |||||

where are the coordinates of the atoms measured from the center of the trap. The density matrix is normalised to unity independent of the number of particles.

The density matrix is then diagonalised, meaning that its single-particle eigenfunctions and the corresponding eigenvalues are calculated,

(9) |

The condensate fraction, , is then defined as the largest eigenvalue.

### 2.5 Bose-Einstein condensate state

For repulsive potential models it is simply the ground state of the trapped many-boson system that is identified as the Bose-Einstein condensate state (BEC-state).

With deep attractive two-body potentials, however, the many-body system in a trap has a large number of self-bound negative-energy states and identification of the BEC-state is not obvious.

The typical spectrum of a trapped many-body system with attractive potentials is shown on fig. 1. The system has a number of deeply bound states with negative energies and then a positive quasi-continuum spectrum which starts at about and has the characteristic distance between levels of the order .

Apparently the BEC-state should then be the lowest state of the quasi-continuum spectrum or, equivalently, the lowest state with positive energy.

To verify this assumption we calculate the central density, , of the system for the negative- and positive-energy states around zero energy. The results are shown on fig. 2 in the form of the inverse central density (the volume per particle) . In the BEC-state the atoms should occupy the whole volume of the trap and thus the volume per particle should be close to one (in the correspondingly scaled oscillator units). The states with negative energy are self-bound states with much higher density and thus much smaller volume per particle. And indeed that is what the figure shows – a sharp increase in the volume per particle from a small value to about unity exactly at the lowest state with positive energy where the quasi-continuum starts.

Another test is the very condensate fraction, shown for several states around zero energy on fig. 3. The self-bound states with negative energy must have smaller condensate fraction compared to the BEC-state, and the excitations from the BEC-state must gradually deplete the condensate fraction. Apparently this is what is seen on the figure – a sharp increase of the condensate fraction to about 100% at the lowest state with positive energy with the subsequent gradual depletion.

We have thus verified that in the case of attractive potentials the BEC-state of a system of trapped bosons is the lowest state with positive energy.

### 2.6 Accuracy of the two-body correlated basis

GP | HS | ZR | A | |
---|---|---|---|---|

3 | 4.51032 | 4.51036(2) | 4.5103 | 4.510 |

5 | 7.53432 | 7.53443(4) | 7.5342 | 7.534 |

10 | 15.1534 | 15.1537(2) | 15.1533 | 15.154 |

20 | 30.638 | 30.640(1) | 30.6394 | 30.640 |

In ref. [9] the energies of several low-density systems of trapped bosons were calculated using the Gross-Pitaevskii and repulsive hard-sphere models with the same ”natural” scattering length of 100 au. In this regime the systems exhibit universality and the energies calculated in both models were very close.

To test the accuracy of our two-body correlated basis eq. (7) (which is expected to be a good approximation in the low-density regime) as well as the identification of the BEC-state for attractive potentials we consider the same systems with the same scattering length but with the attractive potential eq. (4) and calculate the energy of the BEC-state according to our prescription. The BEC-state is now an excited state and is identified in the calculations as the lowest state with positive energy. We also calculate the energies for the zero-range potential model eq. (3) with hyper-radial trial wave-function eq. (6).

The results are given in table 1. As expected, these low-density systems with relatively short scattering length exhibit universality as all potential models give essentially the same results. We conclude that we do a correct identification of the BEC-state and that that the two-body correlated basis has an adequate accuracy.

, au | , au | (2b) | (full) |
---|---|---|---|

-1.400e-7 | 119.4 | 6.025 | 6.025 |

-1.300e-7 | 327.0 | 6.067 | 6.067 |

-1.290e-7 | 402.4 | 6.083 | 6.083 |

-1.280e-7 | 525.4 | 6.108 | 6.108 |

-1.270e-7 | 761.0 | 6.155 | 6.156 |

-1.260e-7 | 1400 | 6.282 | 6.283 |

-1.255e-7 | 2430 | 6.478 | 6.481 |

-1.252e-7 | 4370 | 6.818 | 6.848 |

-1.251e-7 | 5962 | 7.059 | 7.112 |

To check the accuracy of the two-body correlated basis also for large scattering lengths we perform a test calculation for 4 particles with fully correlated and with two-body correlated basis for vastly different scattering lengths. The results are given in table 2.

Although the accuracy of the two-body correlated basis decreases somewhat with the increase of the scattering length, yet the relative accuracy is better than 1% even for exceedingly large scattering lengths.

## 3 Results

### 3.1 Energy

Fig. 4 shows the energy per particle for different and for two interaction models, zero-range potential and attractive potential. For small scattering lengths the different potential models give the same universal results – the system is model independent. For larger scattering lengths the energies from the attractive model are systematically below the zero-range model, quite unlike the repulsive model which goes above the zero-range model [9, 11]. For very large scattering lengths the attractive model, unlike the zero-range and repulsive models, becomes insensitive to the scattering length and the energies converge to a constant. This is consistent with the Jastrow-type approximation of ref. [12].

In the regime close to the two-body threshold, where the scattering length is large, an arbitrary large change in the scattering length needs only an infinitesimally small change of the depth of the attractive potential (see table 2). Therefore when the scattering length is larger than the trap length it ceases to be a physical length scale for the system which is then only subjected to a very small change of the depth of the two-body potential which rapidly converges towards the two-body threshold value. Since the external oscillator potential turns all continuum states into discrete states all singularities due to various thresholds are removed. Then clearly an infinitesimally small change in the potential, despite the large change in the scattering length, only leads to a linearly infinitesimal change in the system which then becomes independent on the scattering length.

On the contrary for the finite-range repulsive potential models an increase of the scattering length needs an almost proportional increase in the effective range of the potential. Thus the system never ceases to depend on the scattering length.

If we plot the data using a different parametrisation, namely as function of , the energy data points seem to follow a ”universal” curve as shown on fig. 5.

### 3.2 Condensate fraction

Our results for the condensate fraction of a system of trapped identical bosons in a BEC-state is shown on fig. 6. For small scattering lengths the system is 100% condensate. Then when is about 3 the condensate fraction rapidly drops a few percent before stabilising again. Although we could not reach further due to numerical convergence problems, we expect that the condensate fraction will not appreciably change with the further increase of the scattering length. This seems consistent with the energies being stagnated when approaching the two-body threshold.

This behavior is qualitatively different from what happens within the repulsive models. In the hard-spheres Monte-Carlo simulations [11] the condensate fraction starts to deviate from 100% much later, when is of the order of 0.1, and with further increase of the scattering length the condensate is completely quenched.

Again our results are different from the Jastrow-type approximation [12], where the condensate also becomes fully quenched. However the authors of ref. [12] believe that their estimates for the condensate fraction are ”rather crude”.

We have to note also that our wave-function includes only two-body correlations which may lead to an over-estimate of the condensate fraction.

## 4 Conclusions

We have calculated the energy and the condensate fraction of a system of bosons in a harmonic trap as function of the number of bosons and the scattering length . Specifically we considered the regime where the scattering length is positive and not small compared to the trap length. The positive scattering length is modeled using an attractive two-body potential with a bound two-body state. The many-body system then has a large number of negative-energy self-bound states and the condensate in the trap is identified as the lowest excited state with positive energy.

When the scattering length is small compared to the trap length the system shows model independence (universality) – the results from the attractive potential model are very close to those from the zero-range and the repulsive potential models.

In the limit of large scattering length the system becomes independent of the scattering length, contrary to the zero-range and the repulsive models.

The condensate fraction decreases with the scattering length and reaches a finite constant at large scattering lengths contrary to the repulsive models where it reaches zero in this limit.

For the attractive potentials the energy per particle of the system of trapped bosons follows a universal curve.

###### Acknowledgements.

Numerous fruitful discussions with H.H. Soerensen and T. Kjaergaard are acknowledged.## References

- [1] Ph. Courteille, R. S. Freeland, and D. J. Heinzen, Phys. Rev. Lett. 81 (1998) 69.
- [2] E. Braaten, and H.-W. Hammer, Physics Reports 428 (2006) 259.
- [3] L.H. Thomas, Phys. Rev. 47 (1935) 903.
- [4] V.N. Efimov, Phys. Lett. 33B (1970) 563.
- [5] C.R.C. Chen et al, Phys. Rev. C 44 (1991) 50; P.F. Bedaque, H.W. Hammer, U van Kolck, Nuclear Physics A 676 (2000) 357; D.V. Fedorov, and A.S. Jensen, Nucl.Phys. A 697 (2002) 783.
- [6] E. Nielsen, and J.H. Macek, Phys. Rev. Lett. 83 (1999) 1566.
- [7] E. Nielsen, D.V. Fedorov, A.S. Jensen, Phys. Rev. A 56 (1997) 3287.
- [8] F Dalfovo, S Giorgini, LP Pitaevskii, and S Stringari, Rev.Mod.Phys. 71 (1999) 463.
- [9] D. Blume, and C.H. Greene, Phys. Rev. A 63 (2001) 063601.
- [10] S. Giorgini, J. Boronat, and J. Casulleras, Phys. Rev. A 60 (1999) 5129.
- [11] J.L. DuBois and H.R. Glyde, Phys.Rev. A 63 (2001) 023602.
- [12] S. Cowell et.al., Phys.Rev.Lett. 88 (2002) 210403.
- [13] D.V. Fedorov and A.S. Jensen, Phys. Rev. A 63 (2001) 063608.
- [14] Y. Suzuki, and K. Varga, Sthochastic Variational Approach to Quantum-Mechanical Few-Body Problems, (Springer-Verlag, Berlin, Heidelberg, 1998).
- [15] H.H. Sørensen, D.V. Fedorov, and A.S. Jensen, AIP Conference Proceedings, 777 (2005) 12.