Ground-state properties of a dilute homogeneous Bose gas of hard disks in two dimensions

The energy and structure of a dilute hard-disks Bose gas are studied in the framework of a variational many-body approach based on a Jastrow correlated ground-state wave function. The asymptotic behaviors of the radial distribution function and the one-body density matrix are analyzed after solving the Euler equation obtained by a free minimization of the hypernetted chain energy functional. Our results show important deviations from those of the available low density expansions, already at gas parameter values x~0.001 . The condensate fraction in 2D is also computed and found generally lower than the 3D one at the same x.

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.

Simulations of Bose fields at finite temperature

We introduce a time-dependent projected Gross-Pitaevskii equation to describe a partially condensed homogeneous Bose gas, and find that this equation will evolve randomized initial wave functions to equilibrium. We compare our numerical data to the predictions of a gapless, second order theory of Bose-Einstein condensation [S. A. Morgan, J. Phys. B 33, 3847 (2000)], and find that we can determine a temperature when the theory is valid. As the Gross-Pitaevskii equation is nonperturbative, we expect that it can describe the correct thermal behavior of a Bose gas as long as all relevant modes are highly occupied. Our method could be applied to other boson fields.

Edge excitations and topological order in rotating Bose gas

The edge excitations and related topological orders of correlated states of a fast rotating Bose gas are studied. Using exact diagonalization of small systems, we compute the energies and number of edge excitations, as well as the boson occupancy near the edge for various states. The chiral Luttinger-liquid theory of Wen is found to be a good description of the edges of the bosonic Laughlin and other states identified as members of the principal Jain sequence for bosons. However, we find that in a harmonic trap the edge of the state identified as the Moore-Read (Pfaffian) state shows a number of anomalies. An experimental way of detecting these correlated states is also discussed.

Ordered structures in rotating ultracold Bose gases

Two-dimentional systems of trapped samples of few cold bosonic atoms submitted to strong rotation around the perpendicular axis may be realized in optical lattices and microtraps. We investigate theoretically the evolution of ground state structures of such systems as the rotational frequency Omega increases. Various kinds of ordered structures are observed. In some cases, hidden interference patterns exhibit themselves only in the pair correlation function; in some other cases explicit broken-symmetry structures appear that modulate the density. For N < 10 atoms, the standard scenario, valid for large sytems is absent, and is only gradually recovered as N increases. On the one hand, the Laughlin state in the strong rotational regime contains ordered structures much more similar to a Wigner molecule than to a fermionic quantum liquid. On the other hand, in the weak rotational regime, the possibility to obtain equilibrium states, whose density reveals an array of vortices, is restricted to the vicinity of some critical values of the rotational frequency Omega.

Edge excitations and topological order in rotating Bose gas

The edge excitations and related topological orders of correlated states of a fast rotating Bose gas are studied. Using exact diagonalization of small systems, we compute the energies and number of edge excitations, as well as the boson occupancy near the edge for various states. The chiral Luttinger-liquid theory of Wen is found to be a good description of the edges of the bosonic Laughlin and other states identified as members of the principal Jain sequence for bosons. However, we find that in a harmonic trap the edge of the state identified as the Moore-Read (Pfaffian) state shows a number of anomalies. An experimental way of detecting these correlated states is also discussed.

The BCS-Bose crossover theory

We contrast four distinct versions of the BCS-Bose statistical crossover theory according to the form assumed for the electron-number equation that accompanies the BCS gap equation. The four versions correspond to explicitly accounting for two-hole-(2h) as well as two-electron-(2e) Cooper pairs (CPs), or both in equal proportions, or only either kind. This follows from a recent generalization of the Bose-Einstein condensation (GBEC) statistical theory that includes not boson-boson interactions but rather 2e- and also (without loss of generality) 2h-CPs interacting with unpaired electrons and holes in a single-band model that is easily converted into a two-band model. The GBEC theory is essentially an extension of the Friedberg-Lee 1989 BEC theory of superconductors that excludes 2h-CPs. It can thus recover, when the numbers of 2h- and 2e-CPs in both BE-condensed and non-condensed states are separately equal, the BCS gap equation for all temperatures and couplings as well as the zero-temperature BCS (rigorous-upper-bound) condensation energy for all couplings. But ignoring either 2h- or 2e-CPs it can do neither. In particular, only half the BCS condensation energy is obtained in the two crossover versions ignoring either kind of CPs. We show how critical temperatures T-c from the original BCS-Bose crossover theory in 2D require unphysically large couplings for the Cooper/BCS model interaction to differ significantly from the T(c)s of ordinary BCS theory (where the number equation is substituted by the assumption that the chemical potential equals the Fermi energy). (c) 2007 Published by Elsevier B.V.

Chaos in collapsing Bose-condensed gas

We reinvestigate the dynamics of the grow and collapse of Bose-Einstein condensates in a system of trapped ultracold atoms with negative scattering lengths, and found a new behavior in the long time scale evolution: the number of atoms can go far beyond the static stability limit. The condensed state is described by the solution of the time-dependent nonlinear Schrödinger equation, in a model that includes atomic feeding and three-body dissipation. Our results for the model show that, by changing the feeding parameter and when a substantial depletion of the ground-state exists, a chaotic behavior is found. We consider a criterion proposed by Deissler and Kaneko [Phys. Lett. A 119, 397 (1987)] to diagnose spatiotemporal chaos. ©2000 The American Physical Society.

Bright vortex solitons in bose condensates

We suggest the possibility of observing and studying bright vortex solitons in attractive Bose-Einstein condensates in three dimensions with a radial trap. Such systems lie on the verge of critical stability and we discuss the conditions of their stability. We study the interaction between two such solitons. Unlike the text-book solitons in one dimension, the interaction between two radially trapped and axially free three-dimensional solitons is inelastic in nature and involves exchange of particles and deformation in shape. The interaction remains repulsive for all phase δ between them except for δ ≈ 0.

Quantum critical scaling at a Bose-glass/superfluid transition: Theory and experiment for a model quantum magnet

In this paper we investigate the quantum phase transition from magnetic Bose Glass to magnetic Bose-Einstein condensation induced by amagnetic field in NiCl2 center dot 4SC(NH2)(2) (dichloro-tetrakis-thiourea-nickel, or DTN), doped with Br (Br-DTN) or site diluted. Quantum Monte Carlo simulations for the quantum phase transition of the model Hamiltonian for Br-DTN, as well as for site-diluted DTN, are consistent with conventional scaling at the quantum critical point and with a critical exponent z verifying the prediction z = d; moreover the correlation length exponent is found to be nu = 0.75(10), and the order parameter exponent to be beta = 0.95(10). We investigate the low-temperature thermodynamics at the quantum critical field of Br-DTN both numerically and experimentally, and extract the power-law behavior of the magnetization and of the specific heat. Our results for the exponents of the power laws, as well as previous results for the scaling of the critical temperature to magnetic ordering with the applied field, are incompatible with the conventional crossover-scaling Ansatz proposed by Fisher et al. [Phys. Rev. B 40, 546 (1989)]. However they can all be reconciled within a phenomenological Ansatz in the presence of a dangerously irrelevant operator.

The stochastic Gross-Pitaevskii equation: II

We provide a derivation of a more accurate version of the stochastic Gross-Pitaevskii equation, as introduced by Gardiner et al (2002 J. Phys. B: At. Mol. Opt. Phys. 35 1555). This derivation does not rely on the concept of local energy and momentum conservation and is based on a quasiclassical Wigner function representation of a 'high temperature' master equation for a Bose gas, which includes only modes below an energy cut-off ER that are sufficiently highly occupied (the condensate band). The modes above this cutoff (the non-condensate band) are treated as being essentially thermalized. The interaction between these two bands, known as growth and scattering processes, provides noise and damping terms in the equation of motion for the condensate band, which we call the stochastic Gross-Pitaevskii equation. This approach is distinguished by the control of the approximations made in its derivation and by the feasibility of its numerical implementation.

Properties of the stochastic Gross-Pitaevskii equation: finite temperature Ehrenfest relations and the optimal plane wave representation

We present Ehrenfest relations for the high temperature stochastic Gross-Pitaevskii equation description of a trapped Bose gas, including the effect of growth noise and the energy cutoff. A condition for neglecting the cutoff terms in the Ehrenfest relations is found which is more stringent than the usual validity condition of the truncated Wigner or classical field method-that all modes are highly occupied. The condition requires a small overlap of the nonlinear interaction term with the lowest energy single particle state of the noncondensate band, and gives a means to constrain dynamical artefacts arising from the energy cutoff in numerical simulations. We apply the formalism to two simple test problems: (i) simulation of the Kohn mode oscillation for a trapped Bose gas at zero temperature, and (ii) computing the equilibrium properties of a finite temperature Bose gas within the classical field method. The examples indicate ways to control the effects of the cutoff, and that there is an optimal choice of plane wave basis for a given cutoff energy. This basis gives the best reproduction of the single particle spectrum, the condensate fraction and the position and momentum densities.

MULTI-BAND BOSE HUBBARD MODELS AND EFFECTIVE THREE-BODY INTERACTIONS

Experiments with ultracold atoms in optical lattice have become a versatile testing ground to study diverse quantum many-body Hamiltonians. A single-band Bose-Hubbard (BH) Hamiltonian was first proposed to describe these systems in 1998 and its associated quantum phase-transition was subsequently observed in 2002. Over the years, there has been a rapid progress in experimental realizations of more complex lattice geometries, leading to more exotic BH Hamiltonians with contributions from excited bands, and modified tunneling and interaction energies. There has also been interesting theoretical insights and experimental studies on “un- conventional” Bose-Einstein condensates in optical lattices and predictions of rich orbital physics in higher bands. In this thesis, I present our results on several multi- band BH models and emergent quantum phenomena. In particular, I study optical lattices with two local minima per unit cell and show that the low energy states of a multi-band BH Hamiltonian with only pairwise interactions is equivalent to an effec- tive single-band Hamiltonian with strong three-body interactions. I also propose a second method to create three-body interactions in ultracold gases of bosonic atoms in a optical lattice. In this case, this is achieved by a careful cancellation of two contributions in the pair-wise interaction between the atoms, one proportional to the zero-energy scattering length and a second proportional to the effective range. I subsequently study the physics of Bose-Einstein condensation in the second band of a double-well 2D lattice and show that the collision aided decay rate of the con- densate to the ground band is smaller than the tunneling rate between neighboring unit cells. Finally, I propose a numerical method using the discrete variable repre- sentation for constructing real-valued Wannier functions localized in a unit cell for optical lattices. The developed numerical method is general and can be applied to a wide array of optical lattice geometries in one, two or three dimensions.

Nonlinear diffraction of light beams propagating in photorefractive media with embedded reflecting wire

The theory of nonlinear diffraction of intensive light beams propagating through photorefractive media is developed. Diffraction occurs on a reflecting wire embedded in the nonlinear medium at a relatively small angle with respect to the direction of the beam propagation. It is shown that this process is analogous to the generation of waves by a flow of a superfluid past an obstacle. The ""equation of state"" of such a superfluid is determined by the nonlinear properties of the medium. On the basis of this hydrodynamic analogy, the notion of the ""Mach number"" is introduced where the transverse component of the wave vector plays the role of the fluid velocity. It is found that the Mach cone separates two regions of the diffraction pattern: inside the Mach cone oblique dark solitons are generated and outside the Mach cone the region of ""optical ship waves"" (the wave pattern formed by a two-dimensional packet of linear waves) is situated. Analytical theory of the ""optical ship waves"" is developed and two-dimensional dark soliton solutions of the generalized two-dimensional nonlinear Schrodinger equation describing the light beam propagation are found. Stability of dark solitons with respect to their decay into vortices is studied and it is shown that they are stable for large enough values of the Mach number.

Charged Kaon Interferometric Probes of Space-Time Evolution in Au plus Au Collisions at s(NN)=200 GeV

Bose-Einstein correlations of charged kaons are used to probe Au+Au collisions at s(NN)=200 GeV and are compared to charged pion probes, which have a larger hadronic scattering cross section. Three-dimensional Gaussian source radii are extracted, along with a one-dimensional kaon emission source function. The centrality dependences of the three Gaussian radii are well described by a single linear function of N(part)(1/3) with a zero intercept. Imaging analysis shows a deviation from a Gaussian tail at r greater than or similar to 10 fm, although the bulk emission at lower radius is well described by a Gaussian. The presence of a non-Gaussian tail in the kaon source reaffirms that the particle emission region in a heavy-ion collision is extended, and that similar measurements with pions are not solely due to the decay of long-lived resonances.