Search

We undertake a theoretical study of a finite-time quantum Otto engine cycle driven by inter-particle interactions in a weakly interacting one-dimensional Bose gas in the quasicondensate regime. Utilizing a $c$-field approach, we simulate the entire Otto cycle, i.e. the two work strokes and the two equilibration strokes. More specifically, the interaction-induced work strokes are modelled by treating the working fluid as an isolated quantum many-body system undergoing unitary evolution. The equilibration strokes, on the other hand, are modelled by treating the working fluid as an open quantum system tunnel-coupled to another quasicondensate which acts as either the hot or cold reservoir, albeit of finite size. We find that, unlike a uniform 1D Bose gas, a harmonically trapped quasicondensate cannot operate purely as a \emph{heat} engine; instead, the engine operation is enabled by additional \emph{chemical} work performed on the working fluid, facilitated by the inflow of particles from the hot reservoir. The microscopic treatment of dynamics during equilibration strokes enables us to evaluate the characteristic operational time scales of this Otto chemical engine, crucial for characterizing its power output, without any \emph{ad hoc} assumptions about typical thermalization timescales. We analyse the performance and quantify the figures of merit of the proposed Otto chemical engine, finding that it offers a favourable trade-off between efficiency and power output, particularly when the interaction-induced work strokes are implemented via a sudden quench. We further demonstrate that in the sudden quench regime, the engine operates with an efficiency close to the near-adiabatic (near maximum efficiency) limit, while concurrently achieving maximum power output., Comment: 24 pages, 7 figures

We present a comprehensive review on the state-of-the-art of the approximate analytic approaches describing the finite-temperature thermodynamic quantities of the Lieb-Liniger model of the one-dimensional (1D) Bose gas with contact repulsive interactions. This paradigmatic model of quantum many-body-theory plays an important role in many areas of physics -- thanks to its integrability and possible experimental realization using, e.g., ensembles of ultracold bosonic atoms confined to quasi-1D geometries. The thermodynamics of the uniform Lieb-Liniger gas can be obtained numerically using the exact thermal Bethe ansatz (TBA) method, first derived in 1969 by Yang and Yang. However, the TBA numerical calculations do not allow for the in-depth understanding of the underlying physical mechanisms that govern the thermodynamic behavior of the Lieb-Liniger gas at finite temperature. Our work is then motivated by the insights that emerge naturally from the transparency of closed-form analytic results, which are derived here in six different regimes of the gas and which exhibit an excellent agreement with the TBA numerics. Our findings can be further adopted for characterising the equilibrium properties of inhomogeneous (e.g., harmonically trapped) 1D Bose gases within the local density approximation and for the development of improved hydrodynamic theories, allowing for the calculation of breathing mode frequencies which depend on the underlying thermodynamic equation of state. Our analytic approaches can be applied to other systems including impurities in a quantum bath, liquid helium-4, and ultracold Bose gas mixtures., Comment: 55 pages, 9 figures; minor revisions in v2 in response to referee's comments

Particle-particle correlations, characterized by Glauber's second-order correlation function,play an important role in the understanding of various phenomena in radio and optical astronomy, quantum and atom optics, particle physics, condensed matter physics, and quantum many-body theory. However, the relevance of such correlations to quantum thermodynamics has so far remained illusive. Here, we propose and investigate a class of quantum many-body thermal machines whose operation is directly enabled by second-order atom-atom correlations in an ultracold atomic gas. More specifically, we study quantum thermal machines that operate in a sudden interaction-quench Otto cycle and utilize a one-dimensional Lieb-Liniger gas of repulsively interacting bosons as the working fluid. The atom-atom correlations in such a gas are different to those of a classical ideal gas, and are a result of the interplay between interparticle interactions, quantum statistics, and thermal fluctuations. We show that operating these thermal machines in the intended regimes, such as a heat engine, refrigerator, thermal accelerator, or heater, would be impossible without such atom-atom correlations. Our results constitute a step forward in the design of conceptually new quantum thermodynamic devices which take advantage of uniquely quantum resources such as quantum coherence, correlations, and entanglement., Comment: 20 pages, 6 figures

We continue the study of dispersive quantum shock waves in a one-dimensional Bose gas beyond the mean-field approximation. In a recent work by Simmons et al. [Phys. Rev. Let. 125, 180401 (2020)], the oscillatory shock wave train developing in this system from an initial localized density bump on a uniform background was interpreted as a result of quantum mechanical self-interference, wherein the interference contrast would diminish with the loss of matter-wave phase coherence. Such loss of coherence, relative to the mean-field Gross-Pitaevskii description, occurs due to either quantum or thermal fluctuations, as well as in the strongly interacting regime. In this work, we extend the analysis of dispersive quantum shock waves in this context to other dynamical scenarios. More specifically, the scenarios studied include evolution of a sufficiently high density bump, known to lead to the so-called ``vacuum point'' in the mean-field description, and evolution of an initial density dip, known to shed a train of grey solitons in the same mean-field approximation. We study the fate of these nonlinear wave structures in the presence of quantum and thermal fluctuations, as well as at intermediate and strong interactions, and show that both the vacuum point and grey solitons cease to manifest themselves beyond the mean-field approach. On the other hand, we find that a vacuum point can occur in an ideal (noninteracting) Bose gas evolving from a ground state of a localized dimple potential. Due to the ubiquity of dispersive shock waves in nature, our results should provide useful insights and perspectives for a variety of other physical systems known to display nonlinear wave phenomena., Comment: 20 pages, 10 figures

This article reviews the recent developments in the theory of generalised hydrodynamics (GHD) with emphasis on the repulsive one-dimensional Bose gas. We discuss the implications of GHD on the mechanisms of thermalisation in integrable quantum many-body systems as well as its ability to describe far-from-equilibrium behaviour of integrable and near integrable systems in a variety of quantum quench scenarios. We outline the experimental tests of GHD in cold-atom gases and its benchmarks with other microscopic theoretical approaches. Finally, we offer some perspectives on the future direction of the development of GHD., Comment: 16 pages, 11 figures

Generalized hydrodynamics (GHD) is a recent theoretical approach that is becoming a go-to tool for characterizing out-of-equilibrium phenomena in integrable and near-integrable quantum many-body systems. Here, we benchmark its performance against an array of alternative theoretical methods, for an interacting one-dimensional Bose gas described by the Lieb-Liniger model. In particular, we study the evolution of both a localized density bump and dip, along with a quantum Newton's cradle setup, for various interaction strengths and initial equilibrium temperatures. We find that GHD generally performs very well at sufficiently high temperatures or strong interactions. For low temperatures and weak interactions, we highlight situations where GHD, while not capturing interference phenomena on short lengthscales, can describe a coarse-grained behaviour based on convolution averaging that mimics finite imaging resolution in ultracold atom experiments. In a quantum Newton's cradle setup based on a double-well to single-well trap quench, we find that GHD with diffusive corrections demonstrates excellent agreement with the predictions of a classical field approach., Comment: 8 pages, 4 figures, plus 6 pages of Supplemental Material

We study the breathing (monopole) oscillations and their damping in a harmonically trapped one-dimensional (1D) Bose gas in the quasicondensate regime using a finite-temperature classical field approach. By characterising the oscillations via the dynamics of the density profile's rms width over long time, we find that the rms width displays beating of two distinct frequencies. This means that 1D Bose gas oscillates not at a single breathing mode frequency, as found in previous studies, but as a superposition of two distinct breathing modes, one oscillating at frequency close to $\simeq\!\sqrt{3}\omega$ and the other at $\simeq\!2\omega$, where $\omega$ is the trap frequency. The breathing mode at $\sim\!\sqrt{3}\omega$ dominates the beating at lower temperatures, deep in the quasicondensate regime, and can be attributed to the oscillations of the bulk of the density distribution comprised of particles populating low-energy, highly-occupied states. The breathing mode at $\simeq\!2\omega$, on the other hand, dominates the beating at higher temperatures, close to the nearly ideal, degenerate Bose gas regime, and is attributed to the oscillations of the tails of the density distribution comprised of thermal particles in higher energy states. The two breathing modes have distinct damping rates, with the damping rate of the bulk component being approximately four times larger than that of the tails component., Comment: 16 pages, 12 figures

Hydrodynamic theories offer successful approaches that are capable of simulating the otherwise difficult-to-compute dynamics of quantum many-body systems. In this work we derive, within the positive-P phase-space formalism, a new stochastic hydrodynamic method for the description of interacting Bose gases. It goes beyond existing hydrodynamic approaches, such as superfluid hydrodynamics or generalized hydrodynamics, in its capacity to simulate the full quantum dynamics of these systems: it possesses the ability to compute non-equilibrium quantum correlations, even for short-wavelength phenomena. Using this description, we derive a linearized stochastic hydrodynamic scheme which is able to simulate such non-equilibrium situations for longer times than the full positive-P approach, at the expense of approximating the treatment of quantum fluctuations, and show that this linearized scheme can be directly connected with existing Bogoliubov approaches. Furthermore, we go on to demonstrate the usefulness and advantages of this formalism by exploring the correlations that arise in a quantum shock wave scenario and comparing its predictions to other established quantum many-body approaches., Comment: 23 pages, 16 figures

We study the non-equilibrium dynamics of two tunnel-coupled one-dimensional quasicondensates following a quench of the coupling strength from zero to a fixed finite value. More specifically, starting from two independent quasicondensates in thermal equilibrium, with initial temperature and chemical potential imbalance, we suddenly switch on the tunnel-coupling and analyze the post-quench equilibration in terms of particle number and energy imbalances. We find that, in certain parameter regimes, the net energy can flow from the colder quasicondensate to the hotter one and is governed by the surplus of low energy particles flowing from the cold to the hot system relative to the high-energy particles flowing in the reverse direction. In all cases, the approach to the new thermal equilibrium occurs through transient, damped oscillations. We also find that for a balanced initial state the coupled quasicondensates can relax into a final thermal equilibrium state in which they display a thermal phase coherence length that is larger than their initial phase coherence length, even though the new equilibrium temperature is higher. The increase in the phase coherence length occurs due to phase locking which manifests itself via an increased degree of correlation between the local relative phases of the quasicondensates at two arbitrary points., Comment: 14 pages, 6 figures; final published version

Shock waves are examples of the far-from-equilibrium behaviour of matter; they are ubiquitous in nature, yet the underlying microscopic mechanisms behind their formation are not well understood. Here, we study the dynamics of dispersive quantum shock waves in a one-dimensional Bose gas, and show that the oscillatory train forming from a local density bump expanding into a uniform background is a result of quantum mechanical self-interference. The amplitude of oscillations, i.e., the interference contrast, decreases with the increase of both the temperature of the gas and the interaction strength due to the reduced phase coherence length. Furthermore, we show that vacuum and thermal fluctuations can significantly wash out the interference contrast, seen in the mean-field approaches, due to shot-to-shot fluctuations in the position of interference fringes around the mean., Comment: Final published version, 6 pages, 3 figures, plus Supplementary Material

We develop a general approach for calculating the characteristic function of the work distribution of quantum many-body systems in a time-varying potential, whose many-body wave function can be cast in the Slater determinant form. Our results are applicable to a wide range of systems including an ideal gas of spinless fermions in one dimension (1D), the Tonks-Girardeau (TG) gas of hard-core bosons, as well as a 1D gas of hard-core anyons. In order to illustrate the utility of our approach, we focus on the TG gas confined to an arbitrary time-dependent trapping potential. In particular, we use the determinant representation of the many-body wave function to characterize the nonequilibrium thermodynamics of the TG gas and obtain exact and computationally tractable expressions---in terms of Fredholm determinants---for the mean work, the work probability distribution function, the nonadiabaticity parameter, and the Loschmidt amplitude. When applied to a harmonically trapped TG gas, our results for the mean work and the nonadiabaticity parameter reduce to those derived previously using an alternative approach. We next propose to use periodic modulation of the trap frequency in order to drive the system to highly non-equilibrium states by taking advantage of the phenomenon of parametric resonance. Under such driving protocol, the nonadiabaticity parameter may reach large values, which indicates a large amount of irreversible work being done on the system as compared to sudden quench protocols considered previously. This scenario is realizable in ultracold atom experiments, aiding fundamental understanding of all thermodynamic properties of the system., Comment: Modified title and abstract, final accepted version (to appear in PRA), 19 pages, 6 figures

We investigate the dynamics of atomic twin beams produced from a phase-fluctuating source, specifically a 1D Bose gas in the quasi-condensate regime, motivated by the experiment reported in Nature Physics 7, 608 (2011). A short-time analytic model is constructed, which is a modified version of the undepleted pump approximation widely used in quantum and atom optics, except that here we take into account the initial phase fluctuations of the pump source as opposed to assuming long-range phase coherence. We use this model to make quantitative and qualitative predictions of how phase-fluctuations of the source impact the two-particle correlations of scattered atom-pairs. The model is benchmarked against detailed numerical simulations using stochastic phase-space methods, and is shown to validate the intuitive notion that the broadening of momentum-space correlation functions between atoms scattered from a quasi-condensate is driven by the broadened momentum width of the source compared to a true phase coherent condensate. Finally, we combine these theoretical tools and results to investigate the effect phase fluctuations of the twin-beam source can have on a proposed demonstration of a violation of a Bell inequality, which intrinsically relies on phase-sensitive pair correlations., Comment: 14 pages, 4 figures

We study the out-of-equilibrium dynamics of a finite-temperature harmonically trapped Tonks-Girardeau gas induced by periodic modulation of the trap frequency. We give explicit exact solutions for the real-space density and momentum distributions of this interacting many-body system and characterize the stability diagram of the dynamics by mapping the many-body solution to the solution and stability diagram of Mathieu's differential equation. The mapping allows one to deduce the exact structure of parametric resonances in the parameter space characterized by the driving amplitude and frequency of the modulation. Furthermore, we analyze the same problem within the finite-temperature hydrodynamic approach and show that the respective solutions to the hydrodynamic equations can be mapped to the same Mathieu equation. Accordingly, the stability diagram and the structure of resonances following from the hydrodynamic approach is exactly the same as those obtained from the exact many-body solution., Comment: 15 pages, 8 figures

We analyse the breathing-mode oscillations of a harmonically quenched Tonks-Giradeau (TG) gas using an exact finite-temperature dynamical theory. We predict a striking collective manifestation of impenetrability---a collective many-body bounce effect. The effect, while being invisible in the evolution of the in-situ density profile of the gas, can be revealed through a nontrivial periodic narrowing of its momentum distribution, taking place at twice the rate of the fundamental breathing-mode frequency. We identify physical regimes for observing the many-body bounce and construct the respective nonequilibrium phase diagram as a function of the quench strength and the initial temperature of the gas. We also develop a finite-temperature hydrodynamic theory of the TG gas, wherein the many-body bounce is explained by an increased thermodynamic pressure of the gas during the isentropic compression, which acts as a potential barrier at the inner turning points of the breathing cycle., Comment: 5 pages, 4 figures, and Supplemental Material. arXiv admin note: substantial text overlap with arXiv:1608.08720

Describing finite-temperature nonequilibrium dynamics of interacting many-particle systems is a notoriously challenging problem in quantum many-body physics. Here we provide an exact solution to this problem for a system of strongly interacting bosons in one dimension in the Tonks-Girardeau regime of infinitely strong repulsive interactions. Using the Fredholm determinant approach and the Bose-Fermi mapping we show how the problem can be reduced to a single-particle basis, wherein the finite-temperature effects enter the solution via an effective "dressing" of the single-particle wavefunctions by the Fermi-Dirac occupation factors. We demonstrate the utility of our approach and its computational efficiency in two nontrivial out-of-equilibrium scenarios: collective breathing mode oscillations in a harmonic trap and collisional dynamics in the Newton's cradle setting involving real-time evolution in a periodic Bragg potential., Comment: Final published version in PRA style; moved Supplemental Material into main text; 6 pages, 3 figures

We consider the Wigner quasi-probability distribution function of a single mode of an electromagnetic or matter-wave field to address the question of whether a direct stochastic sampling and binning of the absolute square of the complex field amplitude can yield a distribution function $\tilde{P}_n$ that closely approximates the true particle number probability distribution $P_n$ of the underlying quantum state. By providing an operational definition of the binned distribution $\tilde{P}_n$ in terms of the Wigner function, we explicitly calculate the overlap between $\tilde{P}_n$ and ${P}_n$ and hence quantify the statistical distance between the two distributions. We find that there is indeed a close quantitative correspondence between $\tilde{P}_n$ and $P_n$ for a wide range of quantum states that have smooth and broad Wigner function relative to the scale of oscillations of the Wigner function for the relevant Fock state. However, we also find counterexamples, including states with high mode occupation, for which $\tilde{P}_n$ does not closely approximate $P_n$., Comment: 9 pages, 7 figures. arXiv admin note: text overlap with arXiv:1503.05647. Accepted for publication in Phys. Rev. A

We develop a finite-temperature hydrodynamic approach for a harmonically trapped one-dimensional quasicondensate and apply it to describe the phenomenon of frequency doubling in the breathing-mode oscillations of its momentum distribution. The doubling here refers to the oscillation frequency relative to the oscillations of the real-space density distribution, invoked by a sudden confinement quench. We find that the frequency doubling is governed by the quench strength and the initial temperature, rather than by the crossover from the ideal Bose gas to the quasicondensate regime. The hydrodynamic predictions are supported by the results of numerical simulations based on a finite-temperature c-field approach, and extend the utility of the hydrodynamic theory for low-dimensional quantum gases to the description of finite-temperature systems and their dynamics in momentum space., Comment: Final published version; 5 pages and 1 figure in the main text, plus Supplemental Material

We use the coordinate Bethe ansatz to exactly calculate matrix elements between eigenstates of the Lieb-Liniger model of one-dimensional bosons interacting via a two-body delta-potential. We investigate the static correlation functions of the zero-temperature ground state and their dependence on interaction strength, and analyze the effects of system size in the crossover from few-body to mesoscopic regimes for up to seven particles. We also obtain time-dependent nonequilibrium correlation functions for five particles following quenches of the interaction strength from two distinct initial states. One quench is from the non-interacting ground state and the other from a correlated ground state near the strongly interacting Tonks-Girardeau regime. The final interaction strength and conserved energy are chosen to be the same for both quenches. The integrability of the model highly constrains its dynamics, and we demonstrate that the time-averaged correlation functions following quenches from these two distinct initial conditions are both nonthermal and moreover distinct from one another., Comment: 17 pages, 8 figures. v2: Final version. Revised in response to referee's comments

We examine the interpretation of individual phase-space trajectories of the Wigner function as corresponding to possible outcomes of single experimental trials. To this end, we investigate the relation between the true (measured) particle number distribution $P_n$ for a single-mode state and that obtained by discretely binning the individual stochastic realisations of squared mode amplitudes $|\alpha|^2$ of the sampled Wigner distribution $W(\alpha)$, which we denote via $\tilde{P}_n$. We provide an operational definition of $\tilde{P}_n$ in terms of the underlying Wigner function, which allows us to explicitly calculate the overlap between the two number distributions and hence quantify the statistical distance between them. We find that there is indeed a close quantitative correspondence between $P_n$ and $\tilde{P}_n$ for a wide range of states, justifying the broadly accepted view that, for highly occupied modes, individual stochastic realisations of Wigner trajectories should approximately correspond to outcomes of single experiments. However, we also find counterexamples for which high mode occupation may not be sufficient for such an interpretation, we find instead that a more relevant and sufficient requirement is the smoothness and broadness of the Wigner function $W(\alpha)$ for the state of interest relative to the scale of oscillations of the Wigner functions for the relevant Fock states., Comment: Our interpretation of results in this article was not adequately articulated. Confusion over terms and phrasing used could lead the reader to tangential issues which were not the original intention or motivation of this paper. A new version can be found as arxiv:1605.07276

We analyze free expansion of a trapped one-dimensional Bose gas after a sudden release from the confining trap potential. By using the stationary phase and local density approximations, we show that the long-time asymptotic density profile and the momentum distribution of the gas are determined by the initial distribution of Bethe rapidities (quasimomenta) and hence can be obtained from the solutions to the Lieb-Liniger equations in the thermodynamic limit. For expansion from a harmonic trap, and in the limits of very weak and very strong interactions, we recover the self-similar scaling solutions known from the hydrodynamic approach. For all other power-law traps and arbitrary interaction strengths, the expansion is not self-similar and shows strong dependence of the density profile evolution on the trap anharmonicity. We also characterize dynamical fermionization of the expanding cloud in terms of correlation functions describing phase and density fluctuations., Comment: Final published version with modified title and a couple of other minor changes. 5 pages, 2 figures, and Supplemental Material

We propose and theoretically simulate an experiment for demonstrating a motional-state Bell inequality violation for pairs of momentum-entangled atoms produced in Bose-Einstein condensate collisions. The proposal is based on realizing an atom-optics analog of the Rarity-Tapster optical scheme: it uses laser-induced Bragg pulses to implement two-particle interferometry on the underlying Bell-state for two pairs of atomic scattering modes with equal but opposite momenta. The collision dynamics and the sequence of Bragg pulses are simulated using the stochastic Bogoliubov approach in the positive-P representation. We predict values of the Clauser-Horne-Shimony-Holt (CHSH) parameter up to S~2.5 for experimentally realistic parameter regimes, showing a strong violation of the CSHS-Bell inequality bounded classically by S<2., Comment: Final published version; 11 pages, 6 figures

In 1963 Glauber introduced the modern theory of quantum coherence, which extended the concept of first-order (one-body) correlations, describing phase coherence of classical waves, to include higher-order (n-body) quantum correlations characterizing the interference of multiple particles. Whereas the quantum coherence of photons is a mature cornerstone of quantum optics, the quantum coherence properties of massive particles remain largely unexplored. To investigate these properties, here we use a uniquely correlated source of atoms that allows us to observe n-body correlations up to the sixth-order at the ideal theoretical limit (n!). Our measurements constitute a direct demonstration of the validity of one of the most widely used theorems in quantum many-body theory--Wisck's theorem--for a thermal ensemble of massive particles. Measurements involving n-body correlations may play an important role in the understanding of thermalization of isolated quantum systems and the thermodynamics of exotic many-body systems, such as Efimov trimers., Comment: 10 pages, 3 figures, plus Supplementary Information

The Hong-Ou-Mandel (HOM) effect is a striking demonstration of destructive quantum interference between pairs of indistinguishable bosons, realised so far only with massless photons. Here we propose an experiment which can realise this effect in the matter-wave regime using pair-correlated atoms produced via a collision of two Bose-Einstein condensates and subjected to two laser induced Bragg pulses. We formulate a novel measurement protocol appropriate for the multimode matter-wave field, which---unlike the typical two-mode optical case---bypasses the need for repeated measurements under different displacement settings of the beam-splitter, thus dramatically reducing the number of experimental runs required to map out the interference visibility. The protocol can be utilised in related matter-wave schemes; here we focus on condensate collisions and by simulating the entire experiment we predict a HOM-dip visibility of ~69%. By being larger than 50%, such a visibility highlights strong quantum correlations between the atoms and paves the way for a possible demonstration of a Bell inequality violation with massive particles in a related Rarity-Tapster setup., Comment: Essentially the same version as v2, except in Nature Communications style; for Supplementary Information see the source file

A fundamental property of a three-dimensional Bose-Einstein condensate (BEC) is long-range coherence, however, in systems of lower dimensionality, not only is the long range coherence destroyed, but additional states of matter are predicted to exist. One such state is a `transverse condensate', first predicted by van Druten and Ketterle [Phys. Rev. Lett. 79, 549 (1997)], in which the gas condenses in the transverse dimensions of a highly anisotropic trap while remaining thermal in the longitudinal dimension. Here we detect the transition from a three-dimensional thermal gas to a gas undergoing transverse condensation by probing Hanbury Brown--Twiss correlations., Comment: 5 pages, 3 figures

We examine the prospect of demonstrating Einstein-Podolsky-Rosen (EPR) entanglement for massive particles using spin-changing collisions in a spinor Bose-Einstein condensate. Such a demonstration has recently been attempted by Gross et al. [Nature 480, 219 (2011)] using a condensate of Rb-87 atoms trapped in an optical lattice potential. For the condensate initially prepared in the (F,m_{F})=(2,0) hyperfine state, with no population in the m_{F}=+-1 states, we predict a significant suppression of the product of inferred quadrature variances below the Heisenberg uncertainty limit, implying strong EPR entanglement. However, such EPR entanglement is lost when the collisions are initiated in the presence of a small (currently undetectable) thermal population n_{th} in the m_{F}=+-1 states. For condensates containing 150 to 200 atoms, we predict an upper bound of n_{th}~1 that can be tolerated in this experiment before EPR entanglement is lost., Comment: 7 pages, 4 figures. Modifications to Figs. 2, 3 and 4 and text

We analyze the two-body momentum correlation function for a uniform weakly interacting one-dimensional Bose gas. We show that the strong positive correlation between opposite momenta, expected in a Bose-Einstein condensate with a true long-range order, almost vanishes in a phase-fluctuating quasicondensate where the long-range order is destroyed. Using the Luttinger liquid approach, we derive an analytic expression for the momentum correlation function in the quasicondensate regime, showing (i) the reduction and broadening of the opposite-momentum correlations (compared to the singular behavior in a true condensate) and (ii) an emergence of anticorrelations at small momenta. We also numerically investigate the momentum correlations in the crossover between the quasicondensate and the ideal Bose-gas regimes using a classical field approach and show how the anticorrelations gradually disappear in the ideal-gas limit., Comment: Final published version

We predict the existence of a dip below unity in the second-order coherence function of a partially condensed ideal Bose gas in harmonic confinement, signaling the anticorrelation of density fluctuations in the sample. The dip in the second-order coherence function is revealed in a canonical-ensemble calculation, corresponding to a system with fixed total number of particles. In a grand-canonical ensemble description, this dip is obscured by the occupation-number fluctuation catastrophe of the ideal Bose gas. The anticorrelation is most pronounced in highly anisotropic trap geometries containing small particle numbers. We explain the fundamental physical mechanism which underlies this phenomenon, and its relevance to experiments on interacting Bose gases., Comment: 10 pages, 5 figures. v2: Minor changes and corrections to figures and text. To appear in PRA

The Cauchy-Schwarz (CS) inequality -- one of the most widely used and important inequalities in mathematics -- can be formulated as an upper bound to the strength of correlations between classically fluctuating quantities. Quantum mechanical correlations can, however, exceed classical bounds.Here we realize four-wave mixing of atomic matter waves using colliding Bose-Einstein condensates, and demonstrate the violation of a multimode CS inequality for atom number correlations in opposite zones of the collision halo. The correlated atoms have large spatial separations and therefore open new opportunities for extending fundamental quantum-nonlocality tests to ensembles of massive particles., Comment: Final published version (with minor changes). 5 pages, 3 figures, plus Supplementary Material

We describe the use of the exact Yang-Yang solutions for the one-dimensional Bose gas to enable accurate kinetic-energy thermometry based on the root-mean-square width of an experimentally measured momentum distribution. Furthermore, we use the stochastic projected Gross-Pitaevskii theory to provide the first quantitative description of the full momentum distribution measurements of Van Amerongen et al., Phys. Rev. Lett. 100, 090402 (2008). We find the fitted temperatures from the stochastic projected Gross-Pitaevskii approach are in excellent agreement with those determined by Yang-Yang kinetic-energy thermometry., Comment: 5 pages, 3 figures. v2: Updated to published version

A Gaussian operator basis provides a means to formulate phase-space simulations of the real- and imaginary-time evolution of quantum systems. Such simulations are guaranteed to be exact while the underlying distribution remains well-bounded, which defines a useful simulation time. We analyse the application of the Gaussian phase-space representation to the dynamics of the dissociation of an ultra-cold molecular gas. We show how the choice of mapping to stochastic differential equations can be used to tailor the stochastic behaviour, and thus the useful simulation time. In the phase-space approach, it is only averages of stochastic trajectories that have a direct physical meaning. Whether particular constants of the motion are satisfied by individual trajectories depends on the choice of mapping, as we show in examples., Comment: Final published version; Comput. Phys. Commun. (2011)

We demonstrate that the quantum dynamics of a many-body Fermi-Bose system can be simulated using a Gaussian phase-space representation method. In particular, we consider the application of the mixed fermion-boson model to ultracold quantum gases and simulate the dynamics of dissociation of a Bose-Einstein condensate of bosonic dimers into pairs of fermionic atoms. We quantify deviations of atom-atom pair correlations from Wick's factorization scheme, and show that atom-molecule and molecule-molecule correlations grow with time, in clear departures from pairing mean-field theories. As a first-principles approach, the method provides benchmarking of approximate approaches and can be used to validate dynamical probes for characterizing strongly correlated phases of fermionic systems., Comment: Final published version

We theoretically analyze dissociation of a harmonically trapped Bose-Einstein condensate of molecular dimers and examine how the spatial inhomogeneity of the molecular condensate affects the conversion dynamics and the atom-atom pair correlations in the short-time limit. Both fermionic and bosonic statistics of the constituent atoms are considered. Using the undepleted molecular-field approximation, we obtain explicit analytic results for the asymptotic behavior of the second-order correlation functions and for the relative number squeezing between the dissociated atoms in one, two and three spatial dimensions. Comparison with the numerical results shows that the analytic approach employed here captures the main underlying physics and provides useful insights into the dynamics of dissociation for conversion efficiencies up to 10%. The results show explicitly how the strength of atom-atom correlations and relative number squeezing degrade with the reduction of the size of the molecular condensate., Comment: Final published version; 23 pages, 13 figures

The behavior of the spatial two-particle correlation function is surveyed in detail for a uniform 1D Bose gas with repulsive contact interactions at finite temperatures. Both long-, medium-, and short-range effects are investigated. The results span the entire range of physical regimes, from ideal gas, to strongly interacting, and from zero temperature to high temperature. We present perturbative analytic methods, available at strong and weak coupling, and first-principle numerical results using imaginary time simulations with the gauge-P representation in regimes where perturbative methods are invalid. Nontrivial effects are observed from the interplay of thermally induced bunching behavior versus interaction induced antibunching., Comment: 22 pages, 13 figures

We describe a pairing mean-field theory related to the Hartree-Fock-Bogoliubov approach, and apply it to the dynamics of dissociation of a molecular Bose-Einstein condensate (BEC) into correlated bosonic atom pairs. We also perform the same simulation using two stochastic phase-space techniques for quantum dynamics -- the positive P-representation method and the truncated Wigner method. By comparing the results of our calculations we are able to assess the relative strength of these theoretical techniques in describing molecular dissociation in one spatial dimension. An important aspect of our analysis is the inclusion of atom-atom interactions which can be problematic for the positive-P method. We find that the truncated Wigner method mostly agrees with the positive-P simulations, but can be simulated for significantly longer times. The pairing mean-field theory results diverge from the quantum dynamical methods after relatively short times., Comment: 11 pages, 7 figures, corrected typos, minor content changes

Ultra-cold clouds of dimeric molecules can dissociate into quantum mechanically correlated constituent atoms that are either bosons or fermions. We theoretically model the dissociation of cigar shaped molecular condensates, for which this difference manifests as complementary geometric structures of the dissociated atoms. For atomic bosons beams form along the long axis of the molecular condensate. For atomic fermions beams form along the short axis. This directional beaming simplifies the measurement of correlations between the atoms through relative number squeezing., Comment: Fixed scrambled reference list

We analyze atom-atom correlations in the s-wave scattering halo of two colliding condensates. By developing a simple perturbative approach, we obtain explicit analytic results for the collinear (CL) and back-to-back (BB) correlations corresponding to realistic density profiles of the colliding condensates with interactions. The results in the short time limit are in agreement with the first-principles simulations using the positive-$P$ representation and provide analytic insights into the experimental observations of Perrin et al. [Phys. Rev. Lett. 99, 150405 (2007)]. For long collision durations, we predict that the BB correlation becomes broader than the CL correlation., Comment: Minor corrections to the text; final published version

We study atom-atom correlations and relative number squeezing in the dissociation of a Bose-Einstein condensate (BEC) of molecular dimers made of either bosonic or fermionic atom pairs. Our treatment addresses the role of the spatial inhomogeneity of the molecular BEC on the strength of correlations in the short time limit. We obtain explicit analytic results for the density-density correlation functions in momentum space, and show that the correlation widths and the degree of relative number squeezing are determined merely by the shape of the molecular condensate., Comment: Minor corrections, final published version

We analytically calculate the spatial nonlocal pair correlation function for an interacting uniform 1D Bose gas at finite temperature and propose an experimental method to measure nonlocal correlations. Our results span six different physical realms, including the weakly and strongly interacting regimes. We show explicitly that the characteristic correlation lengths are given by one of four length scales: the thermal de Broglie wavelength, the mean interparticle separation, the healing length, or the phase coherence length. In all regimes, we identify the profound role of interactions and find that under certain conditions the pair correlation may develop a global maximum at a finite interparticle separation due to the competition between repulsive interactions and thermal effects., Comment: Final published version, modified title

We develop a pairing mean-field theory to describe the quantum dynamics of the dissociation of molecular Bose-Einstein condensates into their constituent bosonic or fermionic atoms. We apply the theory to one, two, and three-dimensional geometries and analyze the role of dimensionality on the atom production rate as a function of the dissociation energy. As well as determining the populations and coherences of the atoms, we calculate the correlations that exist between atoms of opposite momenta, including the column density correlations in 3D systems. We compare the results with those of the undepleted molecular field approximation and argue that the latter is most reliable in fermionic systems and in lower dimensions. In the bosonic case we compare the pairing mean-field results with exact calculations using the positive-$P$ stochastic method and estimate the range of validity of the pairing mean-field theory. Comparisons with similar first-principle simulations in the fermionic case are currently not available, however, we argue that the range of validity of the present approach should be broader for fermions than for bosons in the regime where Pauli blocking prevents complete depletion of the molecular condensate., Comment: 16 pages, 10 figures

We investigate the behavior of a weakly interacting nearly one-dimensional (1D) trapped Bose gas at finite temperature. We perform in situ measurements of spatial density profiles and show that they are very well described by a model based on exact solutions obtained using the Yang-Yang thermodynamic formalism, in a regime where other, approximate theoretical approaches fail. We use Bose-gas focusing [Shvarchuck etal., Phys. Rev. Lett. 89, 270404 (2002)] to probe the axial momentum distribution of the gas, and find good agreement with the in situ results., Comment: extended introduction and conclusions, and minor changes throughout; accepted for publication in Phys. Rev. Lett

We perform first-principles quantum simulations of dissociation of trapped, spatially inhomogeneous Bose-Einstein condensates of molecular dimers. Specifically, we study spatial pair correlations of atoms produced in dissociation after time of flight. We find that the observable correlations may significantly degrade in systems with spatial inhomogeneity compared to the predictions of idealized uniform models. We show how binning of the signal can enhance the detectable correlations and lead to the violation of the classical Cauchy-Schwartz inequality and relative number squeezing., Comment: Final published version

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 non-opposite 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., Comment: Final published version, with updated references and minor modifications

We study a fermionic atom optics counterpart of parametric down-conversion with photons. This can be realized through dissociation of a Bose-Einstein condensate of molecular dimers consisting of fermionic atoms. We present a theoretical model describing the quantum dynamics of dissociation and find analytic solutions for mode occupancies and atomic pair correlations, valid in the short time limit. The solutions are used to identify upper bounds for the correlation functions, which are applicable to any fermionic system and correspond to ideal particle number-difference squeezing, Comment: Changes in response to referees' comments, updated references

We calculate the density profiles and density correlation functions of the one-dimensional Bose gas in a harmonic trap, using the exact finite-temperature solutions for the uniform case, and applying a local density approximation. The results are valid for a trapping potential which is slowly varying relative to a correlation length. They allow a direct experimental test of the transition from the weak coupling Gross-Pitaevskii regime to the strong coupling, 'fermionic' Tonks-Girardeau regime. We also calculate the average two-particle correlation which characterizes the bulk properties of the sample, and find that it can be well approximated by the value of the local pair correlation in the trap center., Comment: Final published version; updated references; 19 pages, 12 figures

In the Comment by M. Mackie \textit{et al.} [arXiv: physics/0212111 v.4], the authors suggest that the molecular conversion efficiency in atom-molecule STIRAP can be improved by lowering the initial atomic density, which in turn requires longer pulse durations to maintain adiabaticity. Apart from the fact that the mean-field approximation becomes questionable at low densities, we point out that a low-density strategy with longer pulses has several problems. It generally requires higher pulse energies, and increases radiative losses. We also show that even within the approximations used in the Comment, their example leads to no efficiency improvement compared to our high-density case. In a more careful analysis including radiative losses neglected in the Comment, the proposed strategy gives almost no conversion owing to the longer pulse durations required., Comment: Accepted for publication in Phys. Rev. A

We present a theoretical analysis of three-dimensional (3D) matter-wave solitons and their stability properties in coupled atomic and molecular Bose-Einstein condensates (BEC). The soliton solutions to the mean-field equations are obtained in an approximate analytical form by means of a variational approach. We investigate soliton stability within the parameter space described by the atom-molecule conversion coupling, atom-atom s-wave scattering, and the bare formation energy of the molecular species. In terms of ordinary optics, this is analogous to the process of sub/second-harmonic generation in a quadratic non-linear medium modified by a cubic nonlinearity, together with a phase mismatch term between the fields. While the possibility of formation of multidimensional spatio-temporal solitons in pure quadratic media has been theoretically demonstrated previously, here we extend this prediction to matter-wave interactions in BEC systems where higher-order non-linear processes due to interparticle collisions are unavoidable and may not be neglected. The stability of the solitons predicted for repulsive atom-atom interactions is investigated by direct numerical simulations of the equations of motion in a full 3D lattice. Our analysis also leads to a possible technique for demonstrating the ground state of the Schroedinger-Newton and related equations that describe Bose-Einstein condensates with non-local inter-particle forces., Comment: Final published version (minor modifications to the text)

Recent experimental measurements of atomic intensity correlations through atom shot noise suggest that atomic quadrature phase correlations may soon be measured with a similar precision. We propose a test of local realism with mesoscopic numbers of massive particles based on such measurements. Using dissociation of a Bose-Einstein condensate of diatomic molecules into bosonic atoms, we demonstrate that strongly entangled atomic beams may be produced which possess Einstein-Podolsky-Rosen (EPR) correlations in field quadratures, in direct analogy to the position and momentum correlations originally considered by EPR., Comment: Final published version (corrections in Ref. [32], updated references)

We analyze molecular bound states of atomic quantum gases near a Feshbach resonance. A simple, renormalizable field theoretic model is shown to have exact solutions in the two-body sector, whose binding energy agrees well with observed experimental results in both Bosonic and Fermionic cases. These solutions, which interpolate between BEC and BCS theories, also provide a more general variational ansatz for resonant superfluidity and related problems., Comment: Minor changes -- to match the final published version

We review progress towards direct simulation of quantum dynamics in many-body systems, using recently developed stochastic gauge techniques. We consider master equations, canonical ensemble calculations and reversible quantum dynamics are compared, as well the general question of strategies for choosing the gauge., Comment: 11 pages, 2 figures, to be published in Proceedings of the 16th International Conference on Laser Spectroscopy (ICOLS), Palm Cove, Australia (2003)