## Gas-liquid coexistence for the bosons square-well fluid and the 4 He binodal anomaly

Riccardo Fantoni 1, ∗

1 Dipartimento di Scienze Molecolari e Nanosistemi, Universit` a Ca' Foscari Venezia, Calle Larga S. Marta DD2137, I-30123 Venezia, Italy

(Dated: August 4, 2014)

The binodal of a boson square-well fluid is determined as a function of the particle mass through the newly devised quantum Gibbs ensemble Monte Carlo algorithm [R. Fantoni and S. Moroni, to be published ]. In the infinite mass limit we recover the classical result. As the particle mass decreases the gas-liquid critical point moves at lower temperatures. We explicitely study the case of a quantum delocalization de Boer parameter close to the one of 4 He. For comparison we also determine the gas-liquid coexistence curve of 4 He for which we are able to observe the binodal anomaly below the λ -transition temperature.

PACS numbers: 05.30.Rt,64.60.-i,64.70.F-,67.10.Fj

Keywords: Quantum statistical physics, path integral Monte Carlo, Quantum Gibbs ensemble Monte Carlo, vapor-liquid phase transition, square well bosons, Helium-4, quantum fluids

Soon after Feynman rewriting of quantum mechanics and quantum statistical physics in terms of the path integral [1, 2] it was realized that the new mathematical object could be used as a powerful numerical instrument.

The statistical physics community soon realized that a path integral could be calculated using the Monte Carlo method [3].

Consider a fluid of N bosons at a given absolute temperature T = 1 /k B β with k B Boltzmann constant. Let the system of particles have a Hamiltonian ˆ H = -λ ∑ N i =1 ∇ 2 i + ∑ i&lt;j φ ( | r i -r j | ) symmetric under particle exchange, with λ = /planckover2pi1 2 / 2 m m , the mass of the particles, and φ ( | r i -r j | ) the pair-potential of interaction between particle i at r i and particle j at r j . The many-particles system will have spatial configurations { R } , with R ≡ ( r 1 , . . . , r N ) the coordinates of the N particles. The partition function of the fluid can be calculated [3] as a sum over the N ! possible particles permutations, P , of a path integral over closed many-particles paths X ≡ ( R , . . . , R 0 P ) in the imaginary time interval τ ∈ [0 , β = P/epsilon1 ], discretized into P intervals of equal length /epsilon1 , the time-step, with R P = P R 0 the β -periodic boundary condition.

extensions of Gibbs ensemble Monte Carlo to include quantum effects (some [8] only consider fluids with internal quantum states; others [9] successfully exploit the path integral Monte Carlo isomorphism between quantum particles and classical ring polymers, but lack the structure of particle exchanges which underlies Bose or Fermi statistics), the QGEMC scheme is viable even for systems with strong quantum delocalization in the degenerate regime of temperature. Details of the QGEMC algorithm will be presented elsewhere [5].

More recently a grand canonical ensemble algorithm has been devised by Massimo Boninsegni et al. [4] for the path integral Monte Carlo method. This paved the way to the development of a quantum Gibbs ensemble Monte Carlo algorithm (QGEMC) to study the gas-liquid coexistence of a generic boson fluid [5]. This algorithm is the quantum analogue of Athanassios Panagiotopoulos [6] method which has now been successfully used for several decades to study first order phase transitions in classical fluids [7]. However, like simulations in the grandcanonical ensemble, the method does rely on a reasonable number of successful particle insertions to achieve compositional equilibrium. As a consequence, the Gibbs ensemble Monte Carlo method cannot be used to study equilibria involving very dense phases. Unlike previous

In this communication we will apply the QGEMC method to the fluid of square well bosons in three spatial dimensions as an extension of the work of Vega et al. [10] on the classical fluid. The de Boer quantum delocalization parameter Λ = /planckover2pi1 /σ m ( E ) 1 / 2 , with E and σ measures of the energy and length scale of the potential energy, can be used to estimate the quantum mechanical effects on the thermodynamic properties of nearly classical liquids [11]. We will consider square well fluids with two values of the particle mass m : Λ = 1 / √ 50, close but different from zero, and Λ = 1 / √ 5. In the first case we compare our result with the one of Vega and in the second case with the one of 4 He which we consider in our second application. When studying the binodal of 4 He in three spatial dimensions we are able to reproduce the binodal anomaly appearing below the λ -point where the liquid branch of the coexistence curve shows a re-entrant behavior.

In our implementation of the QGEMC [5] algorithm we choose the primitive approximation to the path integral action discussed in Ref. [3]. The simulation is performed in two boxes (representing the two coexisting phases) of varying volumes V 1 and V 2 = V -V 1 and numbers of particles N 1 = V ρ 1 1 and N 2 = V ρ 2 2 = N -N 1 with V and N = V ρ constants. The Gibbs equilibrium conditions of pressures and chemical potentials equality between the two boxes is enforced by allowing changes in the volumes of the two boxes (the volume move , q = 5)

and by allowing exchanges of particles between the two boxes (the open-insert move , q = 1, plus the complementary close-remove move , q = 2, plus the advance-recede move , q = 3) while at the same time sampling the closed paths configuration space (the swap move , q = 4, plus the displace move , q = 6, plus the wiggle move , q = 7). We thus have a menu of seven, q = 1 2 , , . . . , 7, different Monte Carlo moves where a single random attempt of any one of them with probability G q = g / q ∑ 7 q =1 g q constitutes a Monte Carlo step.

/negationslash

Letting the system evolve at a given absolute temperature T from a given initial state (for example we shall take ρ 1 = ρ 2 = ρ ) we measure the densities of the two coexisting phases, ρ 1 &lt; ρ and ρ 2 &gt; ρ , which soon approach the coexistence equilibrium values.

We denote with V the maximum displacement of ln( V /V 1 2 ) in the volume move, with L ( p ) the maximum particle displacement in box p = 1 2 in the displacement , move, and with M q &lt; P the maximum number of time slices involved in the q = 5 6 move. , In order to fulfill detailed balance we must choose M 1 = M 2 .

First we study a system of bosons in three dimensions interacting with a square well pair-potential

$$\lambda \quad \phi ( r ) = \begin{cases} + \infty \ r < \sigma & \text{$\alpha$} \\ - \mathcal { A } \ \sigma \leq r < \sigma ( 1 + \Delta ) \quad & ( 1 ) \quad & \text{extr} \\ 0 \quad \sigma ( 1 + \Delta ) \leq r \quad & \text{in $F$} \\ \text{ich, for example, can be used as an effective poten-} \quad \text{me} \\ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \\ \end{cases}$$

which, for example, can be used as an effective potential for cold atoms [12] with a scattering length a = σ (1 + ∆)[1 -tan( σ ∆ √ A / 2 λ /σ ) (1 + ∆) √ A / 2 λ ]. We choose A &gt; 0 as the unit of energies and σ as the unit of lengths. We then introduce a reduced temperature T ∗ = k B T/ A and a reduced density ρ ∗ = ρσ 3 . When the mass of the boson is very big, i.e. λ ∗ = λ/ ( A σ 2 ) /lessmuch 1 we are in the classical limit. The classical fluid has been studied originally by Vega et al. [10] who found that the critical point of the gas-liquid coexistence moves at lower temperatures and higher densities as ∆ gets smaller. The quantum mechanical effects on the thermodynamic properties of nearly classical liquids can be estimated by the de Boer quantum delocalization parameter Λ = √ 2 λ ∗ .

During the subcritical temperature runs we register the densities of the gas, ρ g , and of the liquid, ρ l ( &gt; ρ g ), phase (box). When the densities of the two boxes are too close one another we may observe curves crossing which implies that the two boxes exchange identity. It is then necessary the computation of a density probability distribution function, created using the densities of both boxes. When we are at temperatures sufficiently below the critical point, this distribution appears to be bimodal, i.e. it has two peaks approximated by Gaussians. In some representative cases we checked that the peaks of the bimodal so calculated occur at the same densities as the peaks of the bimodal obtained from the single density distribution of the worm algorithm after a careful tuning of the chemical potential [13].

We study the model with ∆ = 0 5 near their classical . limit λ ∗ = 1 100 (Λ / ≈ 0 14 . , a ∗ = a/σ ≈ 1 44) and at an . intermediate case λ ∗ = 1 10 (Λ / ≈ 0 45 . , a ∗ ≈ 0 58). . We choose N = 50, ρ ∗ = 0 3, . L ( p ) = V 1 / 3 p / 10, V = 1 10, we / take all M q equal, adjusted so as to have the acceptance ratios of the wiggle move close to 50%, g 1 = g 2 = g 3 = g 4 = g 7 = 1, g 5 = 0 0001, and . g 6 = 0 1. . Moreover we choose the relative weight of the Z and G sectors of our extended worm algorithm, C [4], so as to have the Zsector acceptance ratios close to 50%. We started from an initial configuration where we have an equal number of particles in boxes of equal volumes at a total density ρ ∗ = 0 3. .

All our runs were made of 10 5 blocks of 10 5 MC steps with properties measurements every 10 2 steps [14]. The time needed to reach the equilibrium coexistence increases with P and in general with a lowering of the temperature.

If we choose λ ∗ = 1 100 and / P = 2, M q = 1 (in this case the advance-recede move cannot occur) we find that our algorithm gives results close to the ones of Vega [10] obtained with the classical statistical mechanics ( λ ∗ = 0) algorithm of Panagiotopoulos [6] [15]. As we diminish the time-step /epsilon1 ∗ = 1 /PT ∗ at a given temperature we can extrapolate to the zero time-step limit P →∞ as shown in Fig. 2. We thus obtain the fully quantum statistical mechanics result for the binodal shown in Fig. 1 which turns out to exist for T ∗ /lessorsimilar 1. This shows that the critical point due to the effect of the quantum statistics moves at lower temperatures. For the studied temperatures the superfluid fraction [16] of the system was always negligible as in the systems studied in Ref. [9] like Neon (Λ ≈ 0 095) and molecular Hydrogen (Λ . ≈ 0 276). .

/negationslash

/negationslash

In order to extrapolate the binodal to the critical point we used the law of 'rectilinear diameters', ρ l + ρ g = 2 ρ c + a T | -T c | , and the Fisher expansion [17], ρ l -ρ g = b T | -T c | β 1 ( | T -T c | + ) c β 0 -β 1 , with β 1 = 1 2 and / β 0 = 0 3265, and . a, b, c fitting parameters with c = 0 for λ = 0 and c = 0 for λ = 0.

Upon increasing λ ∗ to 1 / 10 the binodal now appears at T ∗ /lessorsimilar 0 008 where we had a non negligible superfluid . fraction [16] ( ρ /ρ s ≈ 0 32(2) at . T ∗ = 0 006 on the liq-. uid branch). As a consequence it proves necessary to use bigger P in the extrapolation to the zero time-step limit. Notice also that at lower temperature it is necessary to run longer simulations due to the longer paths and equilibration times. We generally expect that increasing λ ∗ the gas-liquid critical temperature decreases and the normalsuper fluid critical temperature increases. So the window of temperature for the normal liquid tends to close.

Our second study is on 4 He, for which λ ∗ = 6 0596. . We now take 1 ˚ as unit of lengths and A k B K as unit of energies. In this case σ ≈ 2 5 ˚ , . A E ≈ 10 9K, and Λ . ≈ 0 42. . A situation comparable to the square well case with λ ∗ = 1 10. We use / N = 128 and the Aziz HFDHE2

FIG. 1. (color online) Binodal for the square well fluid in three dimensions. Shown are the classical results of Vega et al. [10] at λ ∗ = 0 and our results in the P → ∞ limit for λ ∗ = 1 100 1 10. In the simulations we used / , / N = 50 and for the extrapolation to the zero time-step limit up to P = 20 for λ ∗ = 1 100 and / P = 500 for λ ∗ = 1 10. / The curves extrapolating to the critical point are obtained as described in the text. The filled triangles are the expected critical points.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0065v1-figure-0989413b-4f80-46c3-baad-322bc85e2b3d.jpg)

FIG. 2. (color online) Linear fit to the zero time-step limit P →∞ for T ∗ = 1 and λ ∗ = 1 100. /

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0065v1-figure-70b33c3e-00aa-4616-b2d4-d664bb6e422f.jpg)

pair-potential [18]

$$\phi ( r ) = \begin{cases} \, \epsilon \phi ^ { * } ( x ) \ r < r _ { \text{cut} } \\ 0 \, & r \geq r _ { \text{cut} } \\ & / C _ { \epsilon } \, & C _ { \epsilon } \, \cdot \, & C _ { \epsilon } \, \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \cdot \, & \dots \end{cases}$$

$$\phi ^ { * } ( x ) = A \exp ( - \alpha x ) - \left ( \frac { C _ { 6 } } { x ^ { 6 } } + \frac { C _ { 8 } } { x ^ { 8 } } + \frac { C _ { 1 0 } } { x ^ { 1 0 } } \right ) F ( x ), \, ( 3 ) \quad \text{a squ} \\ \quad \text{cle m}$$

$$F ( x ) = \begin{cases} \exp [ - ( D / x - 1 ) ^ { 2 } ] & x < D \\ 1 & x \geq D \end{cases}, \quad \begin{pmatrix} \Gamma & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phphi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \PhPh & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Php & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phph & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phip & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Ph i & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phigma & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phia & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Ph  & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phis & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phish & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phih & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phrib & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phig & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phj & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phie & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phsi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phs & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi & \Phi &$$

where x = r/r m , r m = 2 9673, . /epsilon1/k B = 10 8, . A = 0 5448504, . α = 13 353384, . C 6 = 1 3732412, . C 8 = 0 4253785, . C 10 = 0 178100, . D = 1 241314, and . r cut = 6 ˚ A (here we explicitly checked that during the simulation the conditions V 1 / 3 p &gt; 2 r cut for p = 1 2 are always satisfied). ,

In this case it proves convenient to choose ρ ∗ = 0 01, . L ( p ) = V 1 / 3 p / 10, V = 1 10, / g 1 = g 2 = g 3 = g 4 = g 7 = 1, g 5 = 0 0001, and . g 6 = 0 1. As for the SW case we observe . a decrease of the width of the coexistence curve ρ l -ρ g as the number of time slices increases. We thus work at a small (fixed) time-step /epsilon1 ∗ = 0 002 about 1 . / 1000 of the superfluid transition temperature as advised in Ref. [3] to be necessary when studying Helium with the primitive approximation for the action.

The results for the binodal are shown in Fig. 3. The experimental critical point is at T c = 5 25K and . ρ c = 17 3mol/l [19]. . Factors explaining the discrepancy with experiment could be the size error or the choice of the pair-potential. Choosing bigger sizes N it is possible to increase r cut and this shifts the simulated critical temperature to higher values. For the three dimensional 4 He we expect to have the superfluid below a λ -temperature T ∗ λ = 2 193(6) [4], . so our results again show that our method works well even in the presence of a non negligible superfluid fraction. Moreover as shown by the points at the two lowest temperatures we are observing the expected [20] binodal anomaly below the λ -point.

FIG. 3. (color online) Binodal for the 4 He of Aziz [18] in three dimensions. In our simulations we used N = 128, r ∗ cut = 6, and a time-step /epsilon1 ∗ = 0 002. . The continuous (red) curve extrapolating to the critical point are obtained as described in the text. The filled triangle is the estimated critical point. The experimental results from Ref. [19] are also shown as a dashed curve.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0065v1-figure-b31835ce-fe66-45e7-a67b-0f5dfdc9296e.jpg)

In conclusion we determined the gas-liquid binodal of a square well fluid of bosons as a function of the particle mass and of 4 He, in three spatial dimensions, from first principles. The critical point of the square well fluid moves to lower temperatures as the mass of the particles decreases, or as the de Boer parameter increases, while the critical density stays approximately constant.

Our results for 4 He compare well with the experimental critical density even if a lower critical temperature is observed in the simulation. We expect this to be due mainly to a finite size effect unavoidable in the simula-

tion. Nonetheless we are able to determine the binodal anomaly [20] occurring below the λ -transition temperature. The anomaly that we observe in the simulation appears to be more accentuated than in the experiment and the liquid branch of the binodal falls at slightly lower densities.

Even if our QGEMC method is more efficient at high temperatures it is able to detect the liquid phase at low temperatures even below the superfluid transition temperature. The new numerical method is extremely simple to use and unlike current methods does not need the matching of free energies calculated separately for each phase or the simulation of large systems containing both phases and their interface.

R.F. would like to acknowledge the use of the PLX computational facility of CINECA through the ISCRA grant. We are grateful to Michael Ellis Fisher for correspondence and helpful comments.

- ∗ rfantoni@ts.infn.it
- [1] R. P. Feynman, Rev. Mod. Phys. 20 , 367 (1948).
- [2] R. P. Feynman, Statistical Mechanics: A Set of Lectures , Frontiers in Physics, Vol. 36 (W. A. Benjamin, Inc., 1972) notes taken by R. Kikuchi and H. A. Feiveson, edited by Jacob Shaham.
- [3] D. M. Ceperley, Rev. Mod. Phys. 67 , 279 (1995).
- [4] M. Boninsegni, N. Prokof'ev, and B. Svistunov, Phys. Rev. Lett. 96 , 070601 (2006); M. Boninsegni, N. V. Prokof'ev, and B. V. Svistunov, Phys. Rev. E 74 , 036701 (2006).
- [5] R. Fantoni and S. Moroni, (2014), to be published.
- [6] A. Z. Panagiotopoulos, Mol. Phys. 61 , 813 (1987); A. Z. Panagiotopoulos, N. Quirke, M. Stapleton, and D. J. Tildesley, ibid . 63 , 527 (1988); B. Smit, Ph. De Smedt, and D. Frenkel, ibid . 68 , 931 (1989); B. Smit and D.
- Frenkel, ibid . 68 , 951 (1989); D. Frenkel and B. Smit, Understanding Molecular Simulation (Academic Press, San Diego, 1996).
- [7] A. Z. Panagiotopoulos, Mol. Sim. 9 , 1 (1992); F. Sciortino, A. Giacometti, and G. Pastore, Phys. Rev. Lett. 103 , 237801 (2009); R. Fantoni and G. Pastore, Phys. Rev. E 87 , 052303 (2013).
- [8] F. Schneider, D. Marx, and P. Nielaba, Phys. Rev. E 51 , 5162 (1995); P. Nielaba, Int. J. of Thermophys. 17 , 157 (1996).
- [9] Q. Wang and J. K. Johnson, Fluid Phase Equilibria 132 , 93 (1997); I. Georgescu, S. E. Brown, and V. A. Mandelshtam, J. Chem. Phys. 138 , 134502 (2013); P. Kowalczyk, P. A. Gauden, A. P. Terzyk, E. Pantatosaki, and G. K. Papadopoulos, J. Chem. Theory Comput. 9 , 2922 (2013).
- [10] L. Vega, E. de Miguel, L. F. Rull, G. Jackson, and I. A. McLure, J. Chem. Phys. 96 , 2296 (1992); H. Liu, S. Garde, and S. Kumar, ibid . 123 , 174505 (2005).
- [11] R. A. Young, Phys. Rev. Lett. 45 , 638 (1980).
- [12] C. J. Pethik and H. Smith, Bose-Einstein condensation in dilute gases (Cambridge University Press, Cambridge, 2002) chapter 5.
- [13] N. B. Wilding, Phys. Rev. E 52 , 602 (1995).
- [14] Our QGEMC code took ≈ 90 seconds of CPU time for one million steps of a system of size N = 50 , P = 10 , M q = 5 calculating properties every 100 steps, on an IBM iDataPlex DX360M3 Cluster (2.40GHz). The algorithm scales as N 2 , due to the potential energy calculation, and as P , due to the volume move.
- [15] Note that there is no difference between our algorithm in the limit P = 2 , M q = 1, and λ ∗ → 0 and the one of Panagiotopoulos [6].
- [16] E. L. Pollock and D. M. Ceperley, Phys. Rev. B 36 , 8343 (1987).
- [17] M. E. Fisher, Phys. Rev. Lett. 16 , 11 (1966).
- [18] R. A. Aziz, V. P. S. Nain, J. S. Carley, W. L. Taylor, and G. T. McConville, J. Chem. Phys. 70 , 4330 (1979).
- [19] R. D. McCarty, J. Phys. Chem. Ref. Data 2 (1973); NBS TN , 1024 (1980).
- [20] H. Stein, C. Porthun, and G. R¨pke, Eur. Phys. J. B o 2 , 393 (1998).
Riccardo Fantoni
Gas-liquid coexistence for the bosons square-well fluid and the $\mbox{}^4$He binodal anomaly
The binodal of a boson square-well fluid is determined as a function of the
particle mass through the newly devised quantum Gibbs ensemble Monte Carlo
algorithm [R. Fantoni and S. Moroni, {\sl to be published}]. In the infinite
mass limit we recover the classical result. As the particle mass decreases the
gas-liquid critical point moves at lower temperatures. We explicitely study the
case of a quantum delocalization de Boer parameter close to the one of
$\mbox{}^4$He. For comparison we also determine the gas-liquid coexistence
curve of $\mbox{}^4$He for which we are able to observe the binodal anomaly
below the $\lambda$-transition temperature.