A Bethe ansatz solvable model for superpositions of Cooper pairs and condensed molecular bosons

We introduce a general Hamiltonian describing coherent superpositions of Cooper pairs and condensed molecular bosons. For particular choices of the coupling parameters, the model is integrable. One integrable manifold, as well as the Bethe ansatz solution, was found by Dukelsky et al. [J. Dukelsky, G.G. Dussel, C. Esebbag, S. Pittel, Phys. Rev. Lett. 93 (2004) 050403]. Here we show that there is a second integrable manifold, established using the boundary quantum inverse scattering method. In this manner we obtain the exact solution by means of the algebraic Bethe ansatz. In the case where the Cooper pair energies are degenerate we examine the relationship between the spectrum of these integrable Hamiltonians and the quasi-exactly solvable spectrum of particular Schrodinger operators. For the solution we derive here the potential of the Schrodinger operator is given in terms of hyperbolic functions. For the solution derived by Dukelsky et al., loc. cit. the potential is sextic and the wavefunctions obey PT-symmetric boundary conditions. This latter case provides a novel example of an integrable Hermitian Hamiltonian acting on a Fock space whose states map into a Hilbert space of PE-symmetric wavefunctions defined on a contour in the complex plane. (c) 2006 Elsevier B.V. All rights reserved.

Critical Temperature of a Trapped Bose Gas: Comparison of Theory and Experiment

We apply the projected Gross-Pitaevskii equation (PGPE) formalism to the experimental problem of the shift in critical temperature T-c of a harmonically confined Bose gas as reported in Gerbier , Phys. Rev. Lett. 92, 030405 (2004). The PGPE method includes critical fluctuations and we find the results differ from various mean-field theories, and are in best agreement with experimental data. To unequivocally observe beyond mean-field effects, however, the experimental precision must either improve by an order of magnitude, or consider more strongly interacting systems. This is the first application of a classical field method to make quantitative comparison with experiment.

First-principles quantum dynamics in interacting Bose gases II: stochastic gauges

First principles simulations of the quantum dynamics of interacting Bose gases using the stochastic gauge representation are analysed. In a companion paper, we showed how the positive-P representation can be applied to these problems using stochastic differential equations. That method, however, is limited by increased sampling error as time evolves. Here, we show how the sampling error can be greatly reduced and the simulation time significantly extended using stochastic gauges. In particular, local stochastic gauges (a subset) are investigated. Improvements are confirmed in numerical calculations of single-, double- and multi-mode systems in the weak-mode coupling regime. Convergence issues are investigated, including the recognition of two modes by which stochastic equations produced by phase-space methods in general can diverge: movable singularities and a noise-weight relationship. The example calculated here displays wave-like behaviour in spatial correlation functions propagating in a uniform 1D gas after a sudden change in the coupling constant. This could in principle be tested experimentally using Feshbach resonance methods.

Three-body recombination of ultracold Bose gases using the truncated Wigner method

We apply the truncated Wigner method to the process of three-body recombination in ultracold Bose gases. We find that within the validity regime of the Wigner truncation for two-body scattering, three-body recombination can be treated using a set of coupled stochastic differential equations that include diffusion terms, and can be simulated using known numerical methods. As an example we investigate the behavior of a simple homogeneous Bose gas, finding a very slight increase of the loss rate compared to that obtained by using the standard method.

First-principles quantum dynamics in interacting Bose gases: I. The positive P representation

The performance of the positive P phase-space representation for exact many- body quantum dynamics is investigated. Gases of interacting bosons are considered, where the full quantum equations to simulate are of a Gross-Pitaevskii form with added Gaussian noise. This method gives tractable simulations of many-body systems because the number of variables scales linearly with the spatial lattice size. An expression for the useful simulation time is obtained, and checked in numerical simulations. The dynamics of first-, second- and third-order spatial correlations are calculated for a uniform interacting 1D Bose gas subjected to a change in scattering length. Propagation of correlations is seen. A comparison is made with other recent methods. The positive P method is particularly well suited to open systems as no conservation laws are hard-wired into the calculation. It also differs from most other recent approaches in that there is no truncation of any kind.

One-dimensional interacting anyon gas: Low-energy properties and haldane exclusion statistics

The low-energy properties of the one-dimensional anyon gas with a delta-function interaction are discussed in the context of its Bethe ansatz solution. It is found that the anyonic statistical parameter and the dynamical coupling constant induce Haldane exclusion statistics interpolating between bosons and fermions. Moreover, the anyonic parameter may trigger statistics beyond Fermi statistics for which the exclusion parameter alpha is greater than one. The Tonks-Girardeau and the weak coupling limits are discussed in detail. The results support the universal role of alpha in the dispersion relations.

Equation of state of a superfluid Fermi gas in the BCS-BEC crossover

We present a theory for a superfluid Fermi gas near the BCS-BEC crossover, including pairing fluctuation contributions to the free energy similar to that considered by Nozieres and Schmitt-Rink for the normal phase. In the strong coupling limit, our theory is able to recover the Bogoliubov theory of a weakly interacting Bose gas with a molecular scattering length very close to the known exact result. We compare our results with recent Quantum Monte Carlo simulations both for the ground state and at finite temperature. Excellent agreement is found for all interaction strengths where simulation results are available.