Continuous demixing transition of binary liquids: finite-size scaling from the analysis of sub-systems

Abstract

A binary liquid near its consolute point exhibits critical fluctuations of the local composition; the diverging correlation length has always challenged simulations. The method of choice for the calculation of critical points in the phase diagram is a scaling analysis of finite-size corrections, based on a sequence of widely different system sizes. Here, we discuss an alternative using cubic sub-systems of one large simulation as facilitated by modern, massively parallel hardware. We exemplify the method for a symmetric binary liquid at critical composition and compare different routes to the critical temperature: (1) fitting the critical divergences of the correlation length and the susceptibility encoded in the composition structure factor of the whole system, (2) testing data collapse and scaling of moments of the composition fluctuations in sub-volumes, and (3) applying the cumulant intersection criterion to the sub-systems. For the last route, two difficulties arise: sub-volumes are open systems with free boundary conditions, for which no precise estimate of the critical Binder cumulant $U_c$ is available. Second, the periodic boundaries of the simulation box interfere with the sub-volumes, which we resolve by a two-parameter finite-size scaling. The implied modification to the data analysis restores the common intersection point, and we estimate $U_c=0.201\pm 0.001$, universal for cubic Ising-like systems with free boundaries. Confluent corrections to scaling, which arise for small sub-system sizes, are quantified at leading order and our data for the critical susceptibility are compatible with the universal correction exponent $\omega \approx 0.83$.

add discussion and data on confluent corrections to scaling

References in corpus (16)

• Strong scaling of general-purpose molecular dynamics simulations on GPUs
• Solving the 3d Ising Model with the Conformal Bootstrap II. c-Minimization and Precise Critical Exponents
• Static and Dynamic Critical Behavior of a Symmetrical Binary Fluid: A Computer Simulation
• Critical Dynamics in a Binary Fluid: Simulations and Finite-size Scaling
• Enhanced wavelength-dependent surface tension of liquid-vapour interfaces
• Emergence of molecular friction in liquids: bridging between the atomistic and hydrodynamic pictures
• Rotational hot Brownian motion
• Statics and Dynamics of Colloid-Polymer Mixtures Near Their Critical Point of Phase Separation: A Computer Simulation Study of a Continuous AO Model
• Cluster-resolved dynamic scaling theory and universal corrections for transport on percolating systems
• Correction-to-scaling exponent for two-dimensional percolation
• Molecular dynamics of open systems: construction of a mean-field particle reservoir
• Open Systems out of Equilibrium: Theory and Simulation
• Simulation of phase transitions in highly asymmetric fluid mixtures
• Estimation of the critical behavior in an active colloidal system with Vicsek-like interactions
• Finite-size corrections for the static structure factor of a liquid slab with open boundaries
• Anisotropy and universality: Critical Binder cumulant of the two-dimensional Ising model

Cited by in corpus (6)

• Finite-size corrections for the static structure factor of a liquid slab with open boundaries
• Learning the bulk and interfacial physics of liquid-liquid phase separation with neural density functionals
• Structure of liquid--vapor interfaces: perspectives from liquid state theory, large-scale simulations, and potential grazing-incidence X-ray diffraction
• Nonequilibrium induced by reservoirs: Physico-mathematical models and numerical tests
• Gauge invariance and hyperforce correlation theory for equilibrium fluid mixtures
• Critical surface adsorption of confined binary liquids with locally conserved mass and composition