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.

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.

First-principles quantum simulations of dissociation of molecular condensates: Atom correlations in momentum space

We investigate the quantum many-body dynamics of dissociation of a Bose-Einstein condensate of molecular dimers into pairs of constituent bosonic atoms and analyze the resulting atom-atom correlations. The quantum fields of both the molecules and atoms are simulated from first principles in three dimensions using the positive-P representation method. This allows us to provide an exact treatment of the molecular field depletion and s-wave scattering interactions between the particles, as well as to extend the analysis to nonuniform systems. In the simplest uniform case, we find that the major source of atom-atom decorrelation is atom-atom recombination which produces molecules outside the initially occupied condensate mode. The unwanted molecules are formed from dissociated atom pairs with nonopposite momenta. The net effect of this process-which becomes increasingly significant for dissociation durations corresponding to more than about 40% conversion-is to reduce the atom-atom correlations. In addition, for nonuniform systems we find that mode mixing due to inhomogeneity can result in further degradation of the correlation signal. We characterize the correlation strength via the degree of squeezing of particle number-difference fluctuations in a certain momentum-space volume and show that the correlation strength can be increased if the signals are binned into larger counting volumes.

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.

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.

BCS-BEC crossover in an asymmetric two-component Fermi gas

We discuss the superfluid phase transition of a strongly interacting Fermi gas with unequal ( asymmetric) chemical potentials in two pairing hyperfine states, and map out its phase diagram near the BCS-BEC crossover. Our approach includes the fluctuation contributions of preformed Cooper pairs to the thermodynamic potential at finite temperature. We show that, below a critical difference in chemical potentials between species, a normal gas is unstable towards the formation of either a finite-momentum paired Fulde-Ferrell-Larkin-Ovchinnikov superconducting phase or a uniform superfluid, depending on the asymmetry and interaction strengths. We determine the value of critical chemical potential mismatch, and find that it is consistent with a recent measurement by Zwierlein et al. ( Science, 311 ( 2006) 492).

Mean-field phase diagrams of imbalanced Fermi gases near a Feshbach resonance

We propose phase diagrams for an imbalanced (unequal number of atoms or Fermi surface in two pairing hyperfine states) gas of atomic fermions near a broad Feshbach resonance using mean-field theory. Particularly, in the plane of interaction and polarization we determine the region for a mixed phase composed of normal and superfluid components. We compare our prediction of phase boundaries with the recent measurement and find a good qualitative agreement.

Temperature of a trapped unitary Fermi gas at finite entropy

We present theoretical predictions for the equation of state of a harmonically trapped Fermi gas in the unitary limit. Our calculations compare Monte Carlo results with the equation of state of a uniform gas using three distinct perturbation schemes. We show that in experiments the temperature can be usefully calibrated by making use of the entropy, which is invariant during an adiabatic conversion into the weakly interacting limit of molecular BEC. We predict the entropy dependence of the equation of state.

Effect of extrusion parameters on flavour retention, functional and physical properties of mixtures of starch and D-limonene encapsulated in milk protein

The purpose of this research was to investigate the retention of flavour volatiles encapsulated in water-insoluble systems during high temperature–short time extrusion process. A protein precipitation method was used to produce water-insoluble capsules encapsulating limonene, and the capsules were added to the extruder feed material (cornstarch). A twin-screw extruder was used to evaluate the effect of capsule level of addition (0–5%), barrel temperature (125–145 °C) and screw speed (145–175 r.p.m.) on extruder parameters (torque, die pressure, specific mechanical energy, residence time distribution) and extrudate properties [flavour retention, texture, colour, density, expansion, water absorption index, water solubility index (WSI)]. Capsule level had a significant effect on extrusion conditions, flavour retention and extrudate physical properties. Flavour retention increased with the increase in capsule level from 0% to 2.5%, reached a maximum value at capsule level of 2.5% and decreased when the capsule level increased from 2.5% to 5%. The die pressure, torque, expansion ratio, hardness and WSI exhibited the opposite effect with the presence of capsules.

Experiments with Bose-Einstein condensates on an atom chip

We describe the production of BECs on a new type of atom chip based on silver foil. Our atom chip is fabricated with thick wires capable of carrying currents of several amperes without overheating. The silver surface is highly reflective to light resonant with optical transitions used for Rb. The pattern on the chip consists of two parallel Z-trap wires, capable of producing two-wire guide, and two additional endcap wires for varying the axial confinement. Condensates are produced in magnetic microtraps formed within 1 mm of surface of the chip. We have observed the fragmentation of cold atom clouds when brought close to the chip surface. This results from a perturbed trapping potential caused by nanometer deviations of the current path through the wires on the chip. We present results of fragmentation of cold clouds at distances below 100 µm from the wires and investigate the origin of the deviating current. The fragmentation has different characteristics to those seen with copper conductors. The dynamics of atoms in these microtraps is also investigated. ©2005 COPYRIGHT SPIE--The International Society for Optical Engineering. Downloading of the abstract is permitted for personal use only.