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.

Degenerate mixtures of rubidium and ytterbium for engineering open quantum systems

In the last two decades, experimental progress in controlling cold atoms and ions now allows us to manipulate fragile quantum systems with an unprecedented degree of precision. This has been made possible by the ability to isolate small ensembles of atoms and ions from noisy environments, creating truly closed quantum systems which decouple from dissipative channels. However in recent years, several proposals have considered the possibility of harnessing dissipation in open systems, not only to cool degenerate gases to currently unattainable temperatures, but also to engineer a variety of interesting many-body states. This thesis will describe progress made towards building a degenerate gas apparatus that will soon be capable of realizing these proposals. An ultracold gas of ytterbium atoms, trapped by a species-selective lattice will be immersed into a Bose-Einstein condensate (BEC) of rubidium atoms which will act as a bath. Here we describe the challenges encountered in making a degenerate mixture of rubidium and ytterbium atoms and present two experiments performed on the path to creating a controllable open quantum system. The first experiment will describe the measurement of a tune-out wavelength where the light shift of $\Rb{87}$ vanishes. This wavelength was used to create a species-selective trap for ytterbium atoms. Furthermore, the measurement of this wavelength allowed us to extract the dipole matrix element of the $5s \rightarrow 6p$ transition in $\Rb{87}$ with an extraordinary degree of precision. Our method to extract matrix elements has found use in atomic clocks where precise knowledge of transition strengths is necessary to account for minute blackbody radiation shifts. The second experiment will present the first realization of a degenerate Bose-Fermi mixture of rubidium and ytterbium atoms. Using a three-color optical dipole trap (ODT), we were able to create a highly-tunable, species-selective potential for rubidium and ytterbium atoms which allowed us to use $\Rb{87}$ to sympathetically cool $\Yb{171}$ to degeneracy with minimal loss. This mixture is the first milestone creating the lattice-bath system and will soon be used to implement novel cooling schemes and explore the rich physics of dissipation.

Interaction-free measurements by quantum Zeno stabilization of ultracold atoms

Quantum mechanics predicts that our physical reality is influenced by events that can potentially happen but factually do not occur. Interaction-free measurements (IFMs) exploit this counterintuitive influence to detect the presence of an object without requiring any interaction with it. Here we propose and realize an IFM concept based on an unstable many-particle system. In our experiments, we employ an ultracold gas in an unstable spin configuration, which can undergo a rapid decay. The object-realized by a laser beam-prevents this decay because of the indirect quantum Zeno effect and thus, its presence can be detected without interacting with a single atom. Contrary to existing proposals, our IFM does not require single-particle sources and is only weakly affected by losses and decoherence. We demonstrate confidence levels of 90%, well beyond previous optical experiments.

Quantum reflection: The invisible quantum barrier

We construct an invisible quantum barrier which represents the phenomenon of quantum reflection using available data on atom-wall and Bose-Einstein-condensate-wall reflection. We use the Abel equation to invert the data. The resulting invisible quantum barrier is double valued in both axes. We study this invisible barrier in the case of atom and Bose-Einstein condensate (BEC) reflection from a solid silicon surface. A time-dependent, one-spatial-dimension Gross-Pitaevskii equation is solved for the BEC case. We found that the BEC behaves very similarly to the single atom except for size effects, which manifest themselves in a maximum in the reflectivity at small distances from the wall. The effect of the atom-atom interaction on the BEC reflection and correspondingly on the invisible barrier is found to be appreciable at low velocities and comparable to the finite-size effect. The trapping of an ultracold atom or BEC between two walls is discussed.

Inelastic scattering of atoms in a double well

We study a mixture of two light spin-1/2 fermionic atoms and two heavy atoms in a double-well potential. Inelastic scattering processes between both atomic species excite the heavy atoms and renormalize the tunneling rate and the interaction of the light atoms (polaron effect). The effective interaction of the light atoms changes its sign and becomes attractive for strong inelastic scattering. This is accompanied by a crossing of the energy levels from singly occupied sites at weak inelastic scattering to a doubly occupied and an empty site for stronger inelastic scattering. We are able to identify the polaron effect and the level crossing in the quantum dynamics.

Spin Hall effect due to intersubband-induced spin-orbit interaction in symmetric quantum wells

We investigate the intrinsic spin Hall effect in two-dimensional electron gases in quantum wells with two subbands, where a new intersubband-induced spin-orbit coupling is operative. The bulk spin Hall conductivity sigma(z)(xy) is calculated in the ballistic limit within the standard Kubo formalism in the presence of a magnetic field B and is found to remain finite in the B=0 limit, as long as only the lowest subband is occupied. Our calculated sigma(z)(xy) exhibits a nonmonotonic behavior and can change its sign as the Fermi energy (the carrier areal density n(2D)) is varied between the subband edges. We determine the magnitude of sigma(z)(xy) for realistic InSb quantum wells by performing a self-consistent calculation of the intersubband-induced spin-orbit coupling.

Free-fall in a uniform gravitational field in noncommutative quantum mechanics

We study the free-fall of a quantum particle in the context of noncommutative quantum mechanics (NCQM). Assuming noncommutativity of the canonical type between the coordinates of a two-dimensional configuration space, we consider a neutral particle trapped in a gravitational well and exactly solve the energy eigenvalue problem. By resorting to experimental data from the GRANIT experiment, in which the first energy levels of freely falling quantum ultracold neutrons were determined, we impose an upper-bound on the noncommutativity parameter. We also investigate the time of flight of a quantum particle moving in a uniform gravitational field in NCQM. This is related to the weak equivalence principle. As we consider stationary, energy eigenstates, i.e., delocalized states, the time of flight must be measured by a quantum clock, suitably coupled to the particle. By considering the clock as a small perturbation, we solve the (stationary) scattering problem associated and show that the time of flight is equal to the classical result, when the measurement is made far from the turning point. This result is interpreted as an extension of the equivalence principle to the realm of NCQM. (C) 2010 American Institute of Physics. [doi:10.1063/1.3466812]

Time-reversal test for stochastic quantum dynamics

The calculation of quantum dynamics is currently a central issue in theoretical physics, with diverse applications ranging from ultracold atomic Bose-Einstein condensates to condensed matter, biology, and even astrophysics. Here we demonstrate a conceptually simple method of determining the regime of validity of stochastic simulations of unitary quantum dynamics by employing a time-reversal test. We apply this test to a simulation of the evolution of a quantum anharmonic oscillator with up to 6.022×1023 (Avogadro's number) of particles. This system is realizable as a Bose-Einstein condensate in an optical lattice, for which the time-reversal procedure could be implemented experimentally.

Algebraic Bethe ansatz method for the exact calculation of energy spectra and form factors: applications to models of Bose-Einstein condensates and metallic nanograins

In this review we demonstrate how the algebraic Bethe ansatz is used for the calculation of the-energy spectra and form factors (operator matrix elements in the basis of Hamiltonian eigenstates) in exactly solvable quantum systems. As examples we apply the theory to several models of current interest in the study of Bose-Einstein condensates, which have been successfully created using ultracold dilute atomic gases. The first model we introduce describes Josephson tunnelling between two coupled Bose-Einstein condensates. It can be used not only for the study of tunnelling between condensates of atomic gases, but for solid state Josephson junctions and coupled Cooper pair boxes. The theory is also applicable to models of atomic-molecular Bose-Einstein condensates, with two examples given and analysed. Additionally, these same two models are relevant to studies in quantum optics; Finally, we discuss the model of Bardeen, Cooper and Schrieffer in this framework, which is appropriate for systems of ultracold fermionic atomic gases, as well as being applicable for the description of superconducting correlations in metallic grains with nanoscale dimensions.; In applying all the above models to. physical situations, the need for an exact analysis of small-scale systems is established due to large quantum fluctuations which render mean-field approaches inaccurate.

Three-dimensional quantum solitons with parametric coupling

We consider the quantum field theory of two bosonic fields interacting via both parametric (cubic) and quartic couplings. In the case of photonic fields in a nonlinear optical medium, this corresponds to the process of second-harmonic generation (via chi((2)) nonlinearity) modified by the chi((3)) nonlinearity. The quantum solitons or energy eigenstates (bound-state solutions) are obtained exactly in the simplest case of two-particle binding, in one, two, and three space dimensions. We also investigate three-particle binding in one space dimension. The results indicate that the exact quantum solitons of this field theory have a singular, pointlike structure in two and three dimensions-even though the corresponding classical theory is nonsingular. To estimate the physically accessible radii and binding energies of the bound states, we impose a momentum cutoff on the nonlinear couplings. In the case of nonlinear optical interactions, the resulting radii and binding energies of these photonic particlelike excitations in highly nonlinear parametric media appear to be close to physically observable values.

Signature of mott-insulator transition with ultracold fermions in a one-dimensional optical lattice

Using the exact Bethe ansatz solution of the Hubbard model and Luttinger liquid theory, we investigate the density profiles and collective modes of one-dimensional ultracold fermions confined in an optical lattice with a harmonic trapping potential. We determine a generic phase diagram in terms of a characteristic filling factor and a dimensionless coupling constant. The collective oscillations of the atomic mass density, a technique that is commonly used in experiments, provide a signature of the quantum phase transition from the metallic phase to the Mott-insulator phase. A detailed experimental implementation is proposed.

Dynamical tunnelling of ultracold atoms

The divergence of quantum and classical descriptions of particle motion is clearly apparent in quantum tunnelling(1,2) between two regions of classically stable motion. An archetype of such nonclassical motion is tunnelling through an energy barrier. In the 1980s, a new process, 'dynamical' tunnelling(1-3), was predicted, involving no potential energy barrier; however, a constant of the motion (other than energy) still forbids classically the quantum-allowed motion. This process should occur, for example, in periodically driven, nonlinear hamiltonian systems with one degree of freedom(4-6). Such systems may be chaotic, consisting of regions in phase space of stable, regular motion embedded in a sea of chaos. Previous studies predicted(4) dynamical tunnelling between these stable regions. Here we observe dynamical tunnelling of ultracold atoms from a Bose-Einstein condensate in an amplitude-modulated optical standing wave. Atoms coherently tunnel back and forth between their initial state of oscillatory motion (corresponding to an island of regular motion) and the state oscillating 180 degrees out of phase with the initial state.

Twist mode in atomic Fermi gases

The quantum-kinetic energy of a finite number of trapped fermionic atoms provides a restoring force for shear motion due to a distortion of the momentum distribution. In analogy to the twist mode of nuclear physics, it is proposed that counter rotating the upper and lower hemisphere of a spherical atomic cloud yields a finite-frequency mode closely related to transverse zero sound waves in bulk Fermi liquids.

Energy and structure of dilute hard- and soft-sphere gases

The energy and structure of dilute hard- and soft-sphere Bose gases are systematically studied in the framework of several many-body approaches, such as the variational correlated theory, the Bogoliubov model, and the uniform limit approximation, valid in the weak-interaction regime. When possible, the results are compared with the exact diffusion Monte Carlo ones. Jastrow-type correlation provides a good description of the systems, both hard- and soft-spheres, if the hypernetted chain energy functional is freely minimized and the resulting Euler equation is solved. The study of the soft-sphere potentials confirms the appearance of a dependence of the energy on the shape of the potential at gas paremeter values of x~0.001. For quantities other than the energy, such as the radial distribution functions and the momentum distributions, the dependence appears at any value of x. The occurrence of a maximum in the radial distribution function, in the momentum distribution, and in the excitation spectrum is a natural effect of the correlations when x increases. The asymptotic behaviors of the functions characterizing the structure of the systems are also investigated. The uniform limit approach is very easy to implement and provides a good description of the soft-sphere gas. Its reliability improves when the interaction weakens.