Efimov-driven phase transitions of the unitary

Full text

Turn on search term navigation

ARTICLE

Received 4 Sep 2013 | Accepted 24 Feb 2014 | Published 20 Mar 2014

Initially predicted in nuclear physics, Emov trimers are bound congurations of three quantum particles that fall apart when any one of them is removed. They open a window into a rich quantum world that has become the focus of intense experimental and theoretical research, as the region of unitary interactions, where Emov trimers form, is now accessible in cold-atom experiments. Here we use a path-integral Monte Carlo algorithm backed up by theoretical arguments to show that unitary bosons undergo a rst-order phase transition from a normal gas to a superuid Emov liquid, bound by the same effects as Emov trimers. A triple point separates these two phases and another superuid phase, the conventional BoseEinstein condensate, whose coexistence line with the Emov liquid ends in a critical point. We discuss the prospects of observing the proposed phase transitions in cold-atom systems.

DOI: 10.1038/ncomms4503

Emov-driven phase transitions of the unitary Bose gas

Swann Piatecki1,* & Werner Krauth1,*

1 Laboratoire de Physique Statistique,cole Normale Suprieure, UPMC, Universit Paris Diderot, CNRS, 24 rue Lhomond, 75005 Paris, France. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to W.K. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:3503 | DOI: 10.1038/ncomms4503 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4503

Astriking analogue to Borremean rings, Emov trimers are bound congurations of three quantum particles that form near the point of unitary interactions1, where the

pair potential becomes too weak to bind any two of them. Emov trimers open a window into a largely uncharted world of quantum physics based on many-particle bound states. They were discovered in nuclear physics24, but have also been discussed in quantum magnets5, biophysics of DNA6 and, most importantly, in ultra-cold atomic gases. In these systems, it has become possible to ne-tune both the sign and the value of the pair interactions through the Feshbach mechanism7, from a scattering length a close to zero up to the unitary point |a| N. In current

experiments, Emov trimers have not been seen directly, but their presence has been traced through the variations of the rate at which the gas loses particles as the interactions are scanned through the unitary region811.

Theoretical research on Emov physics has unveiled its universal nature, with the three-particle bound states that form an innite sequence at unitarity, and that disappear with zero energy at large negative scattering lengths. Beyond the physics of three-bound particles, the ground-state properties of small unitary clusters were studied numerically12, and their spectrum computed with a variational ansatz13 in a trap. However, the macroscopic many-body properties of the unitary Bose gas have remained unknown. The understanding of its thermodynamic behaviour is of great importance, especially as the experimental stability of the unitary Bose gas of cold atoms has been reported for appreciable time scales11,14,15.

In this work, we apply a dedicated Path-Integral Monte Carlo algorithm to the unitary Bose gas. This allows us to address the thermodynamics of unitary bosons at nite temperature, both above and below BoseEinstein condensation. We obtain the phase diagram of a nite number of trapped bosons, and back up our numerical calculations by a general theoretical model that yields a homogeneous phase diagram. At high temperature, we nd that the unitary Bose gas is very well described by the available virial coefcients16. At lower temperature, the unitary Bose gas undergoes a rst-order phase transition to a new superuid phase, the Emov liquid, held together by the same effects as

Emov trimers, and whose physical quantities are given by three-body observables. From these two phases, transition lines to a third phase, the conventional BoseEinstein condensate (BEC), start at a triple point. The coexistence line between the Emov liquid and the conventional BEC ends in a critical point at high temperature. At a difference with the experimental systems, our model is thermodynamically stable. This is assured through a parameter, the three-body hard core R0, that bounds from below the energies of the Emov states. In cold-atom systems, this parameter is on the order of the van der Waals length lvdW1719. Experimental cold-atom systems are metastable, and particles disappear from the trap into deeply-bound states via the notorious three-body losses. This intricate quantum dynamics, and the description of the losses, are naturally beyond our exact Quantum Monte Carlo approach, and we discuss them on a phenomenological level in order to assess the prospects for experimental tests of our predictions.

ResultsModel Hamiltonian. We describe the system of N interacting bosons by a Hamiltonian

H X

2m X

ioj

Va2rij X

iojok

V3Rijk

p2i

mo2 2

Xix2i; 1

where pi, xi and m are the momentum, the position and the mass of particle i, respectively. The pair interaction Va2 has zero range and may be viewed, as illustrated in Fig. 1ac, as a squarewell interaction potential whose range r0 and depth V0 are simultaneously taken to 0 and N while keeping the scattering length a constant. The unitary point, where the only bound state disappears with zero energy and innite extension, corresponds to an innite scattering length a. The three-body interaction V3 implements a hard-core hyperradial cutoff condition, Rijk4R0,

where the hyperradius Rijk of particles i, j and k corresponds to their root mean square pair distance 3R2ijk r2ij r2ik r2jk

This three-body hard core prevents the so-called Thomas collapse20 into a many-body state with vanishing extension and innite negative energy by setting a fundamental trimer

V2(r )

p R(R/R 0)

[afii9845] r (r/R 0)

102

Imaginary time ([afii9848])

[afii9826]

6 102

V2(r )

2 [afii9845] r ( r / R 0 )

(r/R 0)

0r0 0r0

0r0

Unbound particles R0/a 1

Unitarity R0/a = 0

Bound pair+particle R0/a 1

0 0 20 40 Root mean square distance R/R0

0 0 20 40 Pair separation r/R0

Figure 1 | Path-integral representation of three bosons at different scattering lengths. The bosons are also at R0/lth 1.3 10 2 (for the graphical

projection of (x, y, z and t) to three dimensions, see the Methods section). (a) At R0/aB 1, bosons are unbound and uctuate on a scale lth.

(b) At unitarity (R0/a 0), pairs of bosons form and break up throughout the imaginary time (yellow highlights, arrows), forming a three-body state bound

by pair effects. (c) At R0/aB1, two bosons bind into a stable dimer (red, blue) and one boson is unbound (green). Insets in (ac) correspond to nite-range versions of the zero-range interaction used in each case. Cyan levels correspond to the dimer energy. (d) Sample-averaged hyperradial (root mean square pair distance) probability pR(R) at constant t (solid yellow) compared with its analytic zero-temperature value (dashed blue). (e) Sample-averaged pair distribution rr(r), diverging p1/r2 as r-0 and asymptotically constant r-shell probability r2rr(r), in agreement with the BethePeierls condition.

Data in this gure concern co-cyclic particles in a shallow trap of oB0.

2 NATURE COMMUNICATIONS | 5:3503 | DOI: 10.1038/ncomms4503 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4503 ARTICLE

energy Et / R 20. The nal term in equation (1) models an

isotropic harmonic trapping potential of length ao

h= mo

as it is realized in ultra-cold bosons experiments. The properties of the system described by equation (1) are universal when R0 is

much smaller than all other length scales.

Path-integral representation and Emov trimers. In thermal equilibrium, within the path-integral representation of quantum systems that we use for our computations, position variables xi(t)

carry an imaginary time index tA[0, b 1/kBT], where kB is the

Boltzmann constant and T is the temperature. The uctuations of xi(t) along t account for the quantum uncertainty. The bosonic nature of many-particle systems manifests itself through the periodic boundary conditions in t and, in particular, through the permutation structure of particles. The length of permutation cycles correlates with the degree of quantum coherence2123. Interactions set the statistical weights of congurations24,25. Ensemble averaging, performed by a dedicated Quantum Monte Carlo algorithm, yields the complete thermodynamics of the system (see the Methods section for computational details). The N-body simulation code is massively run on a cluster of independent processors. It succeeds in equilibrating samples with up to a few hundred bosons.

Figure 1 presents snapshots of three bosons in a shallow trap at oB0 that illustrate the quantum uctuations in xi(t), and

characterize the Emov trimer (see the Methods section for a full description of the three bosons simulation). Indeed, for two-body interactions without a bound state, virtually free particles

uctuate on the scale of the de Broglie thermal wavelength lth

2p h2b=mq that diverges at low temperature (see Fig. 1a). In

contrast, for positive a, a bound state with energy E

dimer

:2/(ma2) forms in the two-body interaction potential. Two particles bind into a dimer, and the third particle is free (see

Fig. 1c). At unitarity, the bound state of the pair potential is at resonance Edimer 0, and the scattering length a is innite (see

Fig. 1b). At this point, the two-body interaction is scale-free. While two isolated particles do not bind in the three-particle system, pairs of particles approach each other and then dissociate, so that, between t 0 and t b, the identity of close-by partners

changes several times. This coherent particle-pair scattering process, the hallmark of the Emov effect4, is highlighted in Fig. 1b. At small temperatures, the uctuations of this bound state remain on a scale proportional to R0 and do not diverge as lth.

While the two-particle properties are universal at unitarity, the three-boson fundamental trimer state generally depends on the details of the pair interaction. Excited trimers form a geometric sequence of asymptotically universal Emov trimers with an asymptotic ratio of energies En/En 1E515.0 when n-N, where

En is the energy of the n-th excited trimer1. Owing to this large ratio, thermal averages cannot identify individual trimer states other than the fundamental trimer at low temperature. For the model Hamiltonian of equation (1), the ground-state trimer is virtually identical to a universal Emov trimer12, and we obtain excellent agreement of the probability distribution of the hyperradius R with its analytically known distribution pR(R)

(see Fig. 1d)1. This effectively validates our algorithm. Furthermore, the observed quadratic divergence of the pair distance distribution rr(r), leading to an asymptotically constant r2rr(r) for r-0 (see Fig. 1e), checks with the BethePeierls condition for the zero-range unitary potential26.

Equation of state of the unitary Bose gas. In local-density approximation, particles experience an effective chemical

potential m(r) m0 mo2r2/2 that depends on the distance r

from the centre of the trap. This allows us to obtain the grand canonical equation of state (pressure P as a function of m) from the doubly integrated density prole27 obtained from a single simulation run at temperature T. We nd that the equation of state of the normal gas is described very accurately by the virial expansion up to the third order in the fugacity ebm (see Fig. 2 and the Methods section)16. The third-order term is crucial to the description of Emov physics as it is the rst term at which three-body effects appear28. It depends explicitly on T and R0.

Phase transitions in the trapped unitary Bose gas. In the harmonic trap, particles can be in different thermodynamic phases depending on the distance r from its centre. We monitor the correlation between the pair distances and the position in the trap, and are able to track the creation of a drop of high-density liquid at rB0 (see Fig. 3ac). This drop grows as the temperature decreases.

The observed behaviour corresponds to a rst-order normal-gasto-superuid-liquid transition, and is fundamentally different from the second-order free BoseEinstein condensation (see Supplementary Fig. 1). All particles in the drop are linked through coherent close-by particle switches as in Fig. 1b, showing that the drop is superuid. Deep inside the liquid phase, the Quantum Monte Carlo simulation drops out of equilibrium on the available simulation times. Nevertheless, at its onset and for all values of R0, the peak of the pair correlation function is located around 10R0, which indicates that the liquid phase is of constant density n1 / R 30 (see Fig. 3c). At larger values of R0, the density difference

between the trap centre and the outside vanishes continuously, and the phase transition is no longer seen (see Fig. 3df). Beyond this critical point, the peak of the pair correlation also stabilizes around 10R0, which indicates a cross-over to liquid behaviour (see Fig. 3f and the Methods section for additional details on the characterization of the rst-order phase transition).

Model for the normal-gas-to-Emov-liquid transition. Our numerical ndings suggest a theoretical model for the competition between the unitary gas in third-order virial expansion and an incompressible liquid of density nl / R 30 and constant energy

per particle EpEt, as suggested by ground-state computations

for small clusters12 (see Supplementary Fig. 2). For simplicity, we neglect the entropic contributions to the liquid-state free energy EcTS, so that F1E NE. The phase equilibrium is due to the

difference in free energy and in specic volume at the saturated

2.4

1.6

[afii9826][afii9838] thP ([afii9839])

0.8

0.0

3.0 1.5 0.0[afii9826][afii9839]

Figure 2 | Equation of state of unitary bosons in a harmonic trap. The trap is at R0/ao and T=T0BEC 1:05 and N 100. The numerical equation of state

of the normal gas for N 100 (solid cyan line) obtained by ensemble

averaging of congurations27 is compared with the theoretical virial expansion up to the rst (dash-dotted black line), second (dashed black line) and third (solid black line) virial coefcients. The vertical grey line indicates the most central region of the trap used to determine m0 by

comparison with an ideal gas.

NATURE COMMUNICATIONS | 5:3503 | DOI: 10.1038/ncomms4503 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4503

2.5

130

r sep/R 0

130

2.0

Normal gas

Superuid gas

T/T0 BEC

1.5

0 2 4

Superuid Emov liquid

0 2 4

1.0

0.5

0 0.0 1.5 3.0 0.0 1.5 3.0 0.0 1.5 3.0 r/a[afii9853] r/a[afii9853] r/a[afii9853]

Solid

BEC

Figure 3 | Two-dimensional histogram of pair distances and centre-of-mass positions for N 100 bosons. The pair distance is given by rsep and the centre-of-mass position is given by r. Upper panels: rst-order phase transition for R0/ao 0.07. (a) At T=T0BEC 1:7, the distribution is that of

the normal phase. (b) At a slightly lower temperature T=T0BEC 1:6, a second peak with smaller pair distances B10R0 appears in the trap centre.

distances B10R0. Lower panels: smooth dependence of pair distances and densities on temperature. At R0/ao 0.23, the pair distances decrease

smoothly between T=T0BEC 0:9 (d) and T=T0BEC 0:8 (e) and has

stabilized around 10R0 at T=T0BEC 0:7 (f).

0.0 0.1 0.2 0.3 R0/a[afii9853]

m3 5

h3 /

0.4

Superuid Emov liquid

0.2

Normal gas

0.0 0 1

// [afii9838]th R0

kBT /

Figure 4 | Phase diagram of unitary bosons in the trap centre and in homogeneous space. (a) In the harmonic trap, depending on the value of R0/ao, we observe a normal-gas-to-superuid-liquid rst-order phase transition or a conventional second-order BoseEinstein phase transition (solid black lines). The error bars indicate the last value of T at which the system was a pure normal gas, and the rst value of Tat which we observed either the superuid Emov liquid or the BEC. The normal-gas-tosuperuid-liquid phase transition corresponds well to our theoretical model (solid blue line) at high temperatures. Consistently with theoretical predictions, the numerical coexistence lines are qualitatively continued to a triple point and a critical point (dashed black lines). (b) In a homogeneous system, the normal-gas-to-superuid-liquid coexistence line (solid blue) and the conventional BoseEinstein condensation line (solid red) are universal. The predicted divergence of the normal-gas-to-superuid-liquid coexistence line (dashed red) and the phase transition to a solid phase (dashed purple) are non-universal physics specic to our interaction model.

vapour pressure (see the Methods section and Supplementary Fig. 3). We extend the third-order virial expansion to describe the gaseous phase in the region where quantum correlations become important. Because the conventional BoseEinstein condensation is continuous, it still conveys qualitative information about the transition into the superuid Emov liquid in that region. At small R0-0, the coexistence line approaches innite temperatures as the fundamental trimer energy Et / R 20 diverges. For larger values of

R0, the density of the liquid decreases and approaches the one of the gas. The liquidgas transition line ends in a critical point, where both densities coincide. As the liquid is bound by quantum coherence intrinsic to the Emov effect, this critical point must always be inside a superuid, that is, between the Emov liquid and the BEC, which become indistinguishable. The agreement between this approximate theory and numerical calculations for the trap centre phase diagram is remarkable (see Fig. 4a). Beyond the critical point, we no longer observe a steep drop in the density on decreasing the temperature. We also notice that quantum coherence builds up in the gaseous phase, so that only a conventional BoseEinstein condensation takes place. Our numerical results suggest it occurs at a temperature slightly lower than that for the ideal Bose gas29, kBT0BEC ho 0:94N1=3 0:69.

Homogeneous phase diagram of the unitary Bose gas. Our theoretical model also yields a phase diagram for a homogeneous system of unitary bosons (see Fig. 4b) where, in addition, the conventional BoseEinstein condensation is simply modelled by that of free bosons30. In absence of a harmonic trap, only two independent dimensionless numbers may be built, kBT/E, and

P h3=

m3E5

p . As a consequence, the phase diagram in these two dimensionless numbers is independent of the choice of E, that is, of R0. The scaling of the dimensionless pressure with R50 and of the dimensionless temperature with R20 explains that the triple point appears much farther from the critical point in the homogeneous phase diagram than in the trap. We expect

model-dependent non-universal effects in two regions. At high temperature lthooR0, only a classical gas should exist as the quantum uctuations are too small to build up quantum coherence and, in particular, permutations between particles. At high pressure PpT, we expect a classical solid phase driven by entropic effects, as for conventional hard-sphere melting.

DiscussionTo situate the theoretical results presented in this work in an experimental context, we present the data of Fig. 4b in terms of R0n1/3 and lthn1/3 (see Fig. 5). In this diagram, systems with identical values of T and R0 and different densities correspond to straight lines passing through the origin. In the unstable region, we expect a homogeneous system described by the Hamiltonian of equation (1) to phase-separate on the same such line into the superuid Emov liquid and the normal gas or the BEC.

In cold-atom systems, the atomic interactions have deeply-bound states not present in our model, and the strict hyperradial cutoff is absent. Nevertheless, an effective three-body barrier at a universal value B2lvdW18 induces a universal relation between the fundamental trimer energy Et and the van der Waals length ldvW17,19,31,32. In Fig. 5, a vertical line marks the experimentally

4 NATURE COMMUNICATIONS | 5:3503 | DOI: 10.1038/ncomms4503 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4503 ARTICLE

3.0

2.5

Fundamental trimer of 133Cs

Emov liquid

2.0

[afii9838] thn1/3

1.5

Phase separation

BEC

[afii9838]th R0

1.0

0.5

Normal gas

0.0

0.00 0.05 0.10 0.15 0.20

R0n1/3

Figure 5 | Homogeneous phase diagram of unitary bosons: comparison with experimental densities. Data of Fig. 4b expressed through the dimensionless numbers lthn1/3 and R0n1/3. In the coexistence region, the system phase separates into the superuid Emov liquid (thick blue line)

and, depending on temperature, the normal gas or the BEC. At high density, the BEC undergoes a cross-over to Emov liquid behaviour (thick grey line). The value of R0n1/3B0.023 corresponding to the fundamental trimer of 133Cs at density n 1013 cm 3 is indicated by the vertical

magenta line.

realistic value R0n1/3B0.023 obtained for 133Cs atoms at a density n 1013cm 3, using R0B2lvdW and lvdW 101a07, where a0 is

the Bohr radius. Other atomic species yield smaller values of R0n1/3, for which this discussion is also valid. An experimental system may be quenched along this line into the unstable region by suddenly increasing the scattering length up to unitarity via the Feshbach mechanism. If this quench starts from a weakly interacting BEC at very low temperature, the four-body recombination process (that saturates at a rate K4 hl7th=m)

will dominate the three-body recombination process (that saturates at a rate K3 hl4th=m)11,3335, a condition that may

be written as nK4\K3 (cf. refs 36,37). Whereas the latter is responsible for the three-body losses into deeply-bound dimer states, the former may represent one possible strategy for creating Emov trimers.

In the unstable region, the recombination process would create liquid droplets of ever increasing size. Deep inside the unstable region, far away from the coexistence line, we expect this growth to be a fast, barrier-free, runaway process involving spinodal decomposition38,39, rather than a slow activation process of nucleation over a free-energy barrier produced by the competition of bulk and surface energies. In a thermodynamically stable system, complete separation on macroscopic length scales proceeds through a coarsening process on length scales that slowly increase with time38,40. At the BoseEinstein condensation temperature of the non-interacting gas, this would lead to a large liquid fraction 1 ng/n, where ng is the density of

the gas at coexistence, and n the density of the system before phase separation. However, this will certainly not be observable in current cold-atom experiments, as sufciently large droplets of dense liquid are unstable towards decay into deeply-bound atomic states. Nevertheless, the instability of the gas and the creation of microscopic liquid droplets might be quite realistic. These could be observed as they would lead to an important increase of the three-body losses, proportional to K3n3l in the liquid phase.

The time scales for the creation of a minority liquid phase after a quench is a classic problem of non-equilibrium statistical mechanics. In the present context, it is rendered even richer by the fact that the initial low-temperature phase is a weakly interacting superuid, and that even the normal gas phase is thermodynamically unstable. These questions are closely related to the very existence of the unitary Bose gas on time scales larger

than its thermalization time. Very recent experimental works indicate that the ultra-cold unitary Bose gas can indeed be stabilized on appreciable time scales11,14,15.

Methods

Path-integral Monte Carlo algorithm. In our path-integral quantum Monte Carlo simulation, the contribution of the three-body hard-core interaction to the statistical weights of discretized path congurations are computed using Trotters approximation, which consists in simply rejecting congurations with RoR0. The

contribution of the zero-range unitary interaction is computed using the pair-product approximation, which estimates the weight of two nearby particles without taking other particles into consideration. Both approximations are valid when the discretization step is small23.

In the simulation, new congurations are built from existing congurations according to several possible update moves. A new one, the compression-dilation move, was introduced to specically address the divergence of the pair correlation function at small distances (see Fig. 1e; Supplementary Fig. 4). For each set of parameters, the simulation was run on up to 16 independent processors for up to 10 days to reduce the statistical error.

Simulation of three bosons. In our simulations of unitary bosons, the system is contained in a harmonic trap. This regulates the available conguration space. For the same purpose, in the three-body calculations at oB0 presented in Fig. 1, we impose that the three bosons are on a single permutation cycle: in Fig. 1ac the blue, red and green bosons are, respectively, exchanged with the red, green and blue bosons at imaginary time b. This condition does not modify the properties of the fundamental trimers at unitarity, as other permutations could be sampled at no cost at points of close encounter such as those highlighted in Fig. 1b.

Figure 1ac presents four-dimensional co-cyclic path-integral congurations in three-dimensional plots. In this graphic representation, the centre of mass, whose motion is decoupled from the effect of the interactions, is set to zero at all t. The three spatial dimensions are then reduced to two dimensions by rotating the triangle formed by the three particles at each imaginary time to the same plane in a way that does not favour any of the three spatial dimensions while conserving the permutation cycle structure and the pair distances.

High-temperature equation of state. Within the local-density approximation, the grand canonical equation of state P(m) may be obtained from the numerical doubly integrated density prole

nx27 as

Pmx

mo2

nx; 2

where m(x) m0 mo2x2/2 is the local chemical potential along direction x,

and m0 the chemical potential at the centre of the trap, measured from a t of the equation of state to that of an ideal gas in the outer region of the trap.

We compare this numerical equation of state to the cluster expansion that expresses the pressure in terms of the fugacity ebm (see Fig. 2):

kBT

l3th

l 1

blelbm: 3

The l-th cluster integral bl follows from the virial coefcients of smaller order. It represents l-body effects that cannot be reduced to smaller non-interacting groups of interacting particles28. We use the analytical expressions of b2 and b3 at unitarity that have become available16.

Monitoring the phase transitions. When the densities of the gas and the superuid Emov liquid approach each other, observing directly the two-dimensional histogram of pair distances and centre-of-mass positions does not allow to distinguish between a weakly rst-order phase transition and a cross-over (see Supplementary Fig. 5). In this regime, we monitor the normal-gas-to-superuid-liquid phase transition more accurately by following the evolution of the rst peak of the pair correlation function (obtained by ensemble averaging) with temperature (see Supplementary Fig. 6). In Fig. 4, BoseEinstein condensation is assumed when particles lie on a permutation cycle of length 410 with probability 0.05 (ref. 22).

Analytical model for the transition into the Emov liquid. First-order phase transitions take place when the free energy F of a homogeneous physical system is not a convex function of its volume V. Splitting the system into two phases is then favourable over keeping it homogeneous (see Supplementary Fig. 3). In this situation, at coexistence and in absence of interface energy, the chemical potentials and the pressures of both phases are equal.

The virial expansion is an excellent approximation to describe the normal gas phase of unitary bosons far from the superuid transition. Although this expansion becomes irrelevent in the superuid gas, its analytic continuation conveys important qualitative features because of the continuous nature of conventional BoseEinstein condensation, and is therefore a suitable approximation for the conventional BEC.

NATURE COMMUNICATIONS | 5:3503 | DOI: 10.1038/ncomms4503 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4503

The theoretical model for the superuid liquid is that of an incompressible liquid of specic volume n1 and negligible entropic contribution to the free energy

F1 NE. Simulations yield v 1=31 5R0, and the negligible contribution of the

entropy to the free energy is ensured for non-pathological systems at low temperature. The results in ref. 12 may be extrapolated to E 10.1Et (see

Supplementary Fig. 2), a value adjusted to E 8Et in Fig. 4a to account for nite-

size effects.

In practice, we draw the transition line into the superuid liquid by nding the smallest chemical potential at which the presssures of the incompressible liquid and of the normal gas coincide at a temperature T:

m E

kBT

l3th

ebm b2e2bm b3e3bm: 4

As n 1/n @mP, the crossing to the regime where this equation has no

solution corresponds to the critical point, where both densities are equal.To draw the coexistence line for the trap centre with N 100 particles in

Fig. 4a, its chemical potential m0 is found from integrating the density throughout the trap:

4p l3th

Z r2dr ebmr 2b2e2bmr 3b3e3bmr

N; 5

where the local chemical potential m(r) m0 mo2r2/2 is computed within the

local density approximation.

References

1. Braaten, E. & Hammer, H.-W. Universality in few-body systems with large scattering length. Phys. Rep. 428, 259390 (2006).

2. Emov, V. Energy levels arising from resonant two-body forces in a three-body system. Phys. Lett. B 33, 563564 (1970).

3. Emov, V. Weakly-bound states of three resonantly-interacting particles. Sov. J. Nucl. Phys. 12, 589595 (1971).

4. Emov, V. Low-energy property of three resonantly interacting particles. Sov. J. Nucl. Phys 29, 10581069 (1979).

5. Nishida, Y., Kato, Y. & Batista, C. D. Emov effect in quantum magnets. Nat. Phys. 9, 9397 (2013).

6. Maji, J., Bhattacharjee, S. M., Seno, F. & Trovato, A. When a DNA triple helix melts: an analogue of the Emov state. New J. Phys. 12, 083057 (2010).

7. Chin, C., Grimm, R., Julienne, P. & Tiesinga, E. Feshbach resonances in ultracold gases. Rev. Mod. Phys. 82, 12251286 (2010).

8. Gross, N., Shotan, Z., Kokkelmans, S. & Khaykovich, L. Observation of universality in ultracold 7Li three-body recombination. Phys. Rev. Lett. 103, 163202 (2009).

9. Pollack, S. E., Dries, D. & Hulet, R. G. Universality in three- and four-body bound states of ultracold atoms. Science 326, 16831685 (2009).

10. Zaccanti, M. et al. Observation of an Emov spectrum in an atomic system. Nat. Phys. 5, 586591 (2009).

11. Rem, B. S. et al. Lifetime of the Bose gas with resonant interactions. Phys. Rev. Lett. 110, 163202 (2013).

12. von Stecher, J. Weakly bound cluster states of Emov character. J. Phys. B 43, 101002 (2010).

13. Thgersen, M., Fedorov, D. V. & Jensen, A. S. N-body Emov states of trapped bosons. Europhys. Lett. 83, 30012 (2008).

14. Fletcher, R. J., Gaunt, A. L., Navon, N., Smith, R. P. & Hadzibabic, Z. Stability of a unitary Bose gas. Phys. Rev. Lett. 111, 125303 (2013).

15. Makotyn, P., Klauss, C. E., Goldberger, D. L., Cornell, E. A. & Jin, D. S. Universal dynamics of a degenerate unitary Bose gas. Nat. Phys. 10, 116119 (2014).

16. Castin, Y. & Werner, F. Le troisime coefcient du viriel du gaz de Bose unitaire. Can. J. Phys. 91, 382389 (2013).

17. Berninger, M. et al. Universality of the three-body parameter for Emov states in ultracold cesium. Phys. Rev. Lett. 107, 120401 (2011).

18. Wang, J., DIncao, J. P., Esry, B. D. & Greene, C. H. Origin of the three-body parameter universality in Emov physics. Phys. Rev. Lett. 108, 263001 (2012).

19. Wild, R. J., Makotyn, P., Pino, J. M., Cornell, E. A. & Jin, D. S. Measurements of Tans contact in an atomic Bose-Einstein condensate. Phys. Rev. Lett. 108, 145305 (2012).

20. Thomas, L. H. The interaction between a neutron and a proton and the structure of H3. Phys. Rev. 47, 903909 (1935).

21. Feynman, R. P. Statistical Mechanics: a Set of Lectures (Addison-Wesley, 1982).22. Krauth, W. Quantum Monte Carlo calculations for a large number of bosons in a harmonic trap. Phys. Rev. Lett. 77, 36953699 (1996).

23. Krauth, W. Statistical Mechanics: Algorithms and Computations (Oxford Univ. Press, 2006).

24. Pollock, E. L. & Ceperley, D. M. Simulation of quantum many-body systems by pathintegral methods. Phys. Rev. B 30, 25552568 (1984).

25. Ceperley, D. M. Path integrals in the theory of condensed helium. Rev. Mod. Phys. 67, 279355 (1995).

26. Bethe, H. A. & Peierls, R. The scattering of neutrons by protons. Proc. R. Soc. Lond. A 149, 176183 (1935).

27. Ho, T.-L. & Zhou, Q. Obtaining the phase diagram and thermodynamic quantities of bulk systems from the densities of trapped gases. Nat. Phys. 6, 131134 (2010).

28. Huang, K. Statistical Mechanics 2nd ed. (Wiley, 1987).29. Dalfovo, F., Giorgini, S., Pitaevskii, L. P. & Stringari, S. Theory of Bose-Einstein condensation in trapped gases. Rev. Mod. Phys. 71, 463512 (1999).

30. Landau, L. D. & Lifshitz, L. M. Statistical Physics 3rd ed. Course of theoretical physics No. 5 (Butterworth-Heinemann, 1980).

31. Schmidt, R., Rath, S. P. & Zwerger, W. Emov physics beyond universality. Eur. Phys. J. B 85, 386 (2012).

32. Roy, S. et al. Test of the universality of the three-body Emov parameter at narrow Feshbach resonances. Phys. Rev. Lett. 111, 053202 (2013).

33. DIncao, J. P., Suno, H. & Esry, B. D. Limits on universality in ultracold three-boson recombination. Phys. Rev. Lett. 93, 123201 (2004).

34. Mehta, N. P., Rittenhouse, Seth T., DIncao, J. P., von Stecher, J. & Greene, Chris H. General theoretical description of n-body recombination. Phys. Rev. Lett. 103, 153201 (2009).

35. Li, W. & Ho, T.-L. Bose gases near unitarity. Phys. Rev. Lett. 108, 195301 (2012).

36. Ferlaino, F. et al. Evidence for universal four-body states tied to an Emov trimer. Phys. Rev. Lett. 102, 140401 (2009).

37. Greene, C. H. Universal insights from few-body land. Phys. Today 63, 40 (2010).

38. Binder, K. Theory of rst-order phase transitions. Rep. Prog. Phys. 50, 783859 (1987).

39. Cahn, J. W. & Hilliard, J. E. Free energy of a nonuniform system. I. Interfacial free energy. J. Chem. Phys. 28, 258267 (2004).

40. Lifshitz, I. M. & Slyozov, V. V. The kinetics of precipitation from supersaturated solid solutions. J. Phys. Chem. Solids 19, 3550 (1961).

Acknowledgements

We thank S. Balibar, Y. Castin, F. Chevy, I. Ferrier-Barbut, A.T. Grier, B. Rem,C. Salomon and F. Werner for very helpful discussions. W.K. acknowledges the hospitality of the Aspen Center for Physics, which is supported by the National Science Foundation Grant No. PHY-1066293.

Author contributions

Both authors contributed equally to this work.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/nature

Web End =http://www.nature.com/nature communications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Piatecki, S. and Krauth, W. Emov-driven phase transitions of the unitary Bose gas. Nat. Commun. 5:3503 doi: 10.1038/ncomms4503 (2014).

6 NATURE COMMUNICATIONS | 5:3503 | DOI: 10.1038/ncomms4503 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

Word count: 5756

Show less

Abstract

Translate

Initially predicted in nuclear physics, Efimov trimers are bound configurations of three quantum particles that fall apart when any one of them is removed. They open a window into a rich quantum world that has become the focus of intense experimental and theoretical research, as the region of 'unitary' interactions, where Efimov trimers form, is now accessible in cold-atom experiments. Here we use a path-integral Monte Carlo algorithm backed up by theoretical arguments to show that unitary bosons undergo a first-order phase transition from a normal gas to a superfluid Efimov liquid, bound by the same effects as Efimov trimers. A triple point separates these two phases and another superfluid phase, the conventional Bose-Einstein condensate, whose coexistence line with the Efimov liquid ends in a critical point. We discuss the prospects of observing the proposed phase transitions in cold-atom systems.

Details

Title

Efimov-driven phase transitions of the unitary Bose gas

Author

Piatecki, Swann; Krauth, Werner

Pages

3503

Publication year

2014

Publication date

Mar 2014

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms4503

ProQuest document ID

1508690079

Efimov-driven phase transitions of the unitary Bose gas

Jump to:

Full text

Abstract

Details

Suggested sources