No Arabic abstract
Population annealing is a sequential Monte Carlo scheme well-suited to simulating equilibrium states of systems with rough free energy landscapes. Here we use population annealing to study a binary mixture of hard spheres. Population annealing is a parallel version of simulated annealing with an extra resampling step that ensures that a population of replicas of the system represents the equilibrium ensemble at every packing fraction in an annealing schedule. The algorithm and its equilibration properties are described and results are presented for a glass-forming fluid composed of a 50/50 mixture of hard spheres with diameter ratio of 1.4:1. For this system, we obtain precise results for the equation of state in the glassy regime up to packing fractions $varphi approx 0.60$ and study deviations from the BMCSL equation of state. For higher packing fractions, the algorithm falls out of equilibrium and a free volume fit predicts jamming at packing fraction $varphi approx 0.667$. We conclude that population annealing is an effective tool for studying equilibrium glassy fluids and the jamming transition.
We simulate the motion of spherical particles in a phase-separating binary mixture. By combining cell dynamical equations with Langevin dynamics for particles, we show that the addition of hard particles significantly changes both the speed and the morphology of the phase separation. At the late stage of the spinodal decomposition process, particles significantly slow down the domain growth, in qualitative agreement with earlier experimental data.
For binary fluid mixtures of spherical particles in which the two species are sufficiently different in size, the dominant wavelength of oscillations of the pair correlation functions is predicted to change from roughly the diameter of the large species to that of the small species along a sharp crossover line in the phase diagram [C. Grodon, M. Dijkstra, R. Evans & R. Roth, J.Chem.Phys. 121, 7869 (2004)]. Using particle-resolved colloid experiments in 3d we demonstrate that crossover exists and that its location in the phase diagram is in quantitative agreement with the results of both theory and our Monte-Carlo simulations. In contrast with previous work [J. Baumgartl, R. Dullens, M. Dijkstra, R. Roth & C. Bechinger, Phys.Rev.Lett. 98, 198303 (2007)], where a correspondence was drawn between crossover and percolation of both species, in our 3d study we find that structural crossover is unrelated to percolation.
Population annealing is a recent addition to the arsenal of the practitioner in computer simulations in statistical physics and beyond that is found to deal well with systems with complex free-energy landscapes. Above all else, it promises to deliver unrivaled parallel scaling qualities, being suitable for parallel machines of the biggest calibre. Here we study population annealing using as the main example the two-dimensional Ising model which allows for particularly clean comparisons due to the available exact results and the wealth of published simulational studies employing other approaches. We analyze in depth the accuracy and precision of the method, highlighting its relation to older techniques such as simulated annealing and thermodynamic integration. We introduce intrinsic approaches for the analysis of statistical and systematic errors, and provide a detailed picture of the dependence of such errors on the simulation parameters. The results are benchmarked against canonical and parallel tempering simulations.
Following the work of Leutheusser [Physica A 127, 667 (1984)], the solution to the Percus-Yevick equation for a seven-dimensional hard-sphere fluid is explicitly found. This allows the derivation of the equation of state for the fluid taking both the virial and the compressibility routes. An analysis of the virial coefficients and the determination of the radius of convergence of the virial series are carried out. Molecular dynamics simulations of the same system are also performed and a comparison between the simulation results for the compressibility factor and theoretical expressions for the same quantity is presented.
A cylindrical nanopore immersed in a non-additive hard sphere binary fluid is studied by means of integral equation theories and Monte Carlo simulations. It is found that at low and intermediate values of the bulk total number density the more concentrated bulk species is preferentially absorbed by the pore, as expected. However, further increments of the bulk number density lead to an abrupt population inversion in the confined fluid and an entropy driven prewetting transition at the outside wall of the pore. These phenomena are a function of the pore size, the non-additivity parameter, the bulk number density, and particles relative number fraction. We discuss our results in relation to the phase separation in the bulk.