Fluid coexistence close to criticality: scaling algorithms for precise simulation

A novel algorithm is presented that yields precise estimates of coexisting liquid and gas densities, ρ(T), from grand canonical Monte Carlo simulations of model fluids near criticality. The algorithm utilizes data for the isothermal minima of the moment ratio Q(T;<)≡<m>L2/< in L×⋯×L boxes, where m=ρ-<. When L→∞ the minima, Qm±(T;L), tend to zero while their locations, ρm±(T;L), approach ρ(T) and ρ(T). Finite-size scaling relates the ratio Y=(ρm+-ρm-)/Δρ(T)universally to 1/2 >(Qm++Qm-), where Δρ=ρ(T)-ρ(T) is the desired width of the coexistence curve. Utilizing the exact limiting (L→∞) form, the corresponding scaling function can be generated in recursive steps by fitting overlapping data for three or more box sizes, L,L,…,L. Starting at a T sufficiently far below T and suitably choosing intervals ΔT=T-T>0 yields Δρ(T) and precisely locates T. The algorithm has been applied to simulation data for a hard-core square-well fluid and the restricted primitive model electrolyte for sizes up to L/a=8-12 (where a is the hard-core diameter): the coexistence curves can be computed to a precision of ±1-2% of ρ up to |T-T|/T=10 and 10 ^-3, respectively. Universality of the scaling functions and the exponent β is verified and the (T,ρ) estimates confirm previous values based on data from above T. The algorithm extends directly to calculating the diameter, ρ(T)≡1/2 >(ρ+ρ), and can lead to estimates of the Yang-Yang ratio. Furthermore, a new, explicit approximant for the basic scaling function Y permits straightforward estimates of Δρ(T) from limited Q-data when Ising-type criticality may be assumed.