An improved Monte Carlo method for direct calculation of the density of states

We present an efficient Monte Carlo algorithm for determining the density of states which is based on the statistics of transition probabilities between states. By measuring the infinite temperature transition probabilities--that is, the probabilities associated with move proposal only--we are able to extract excellent estimates of the density of states. When this estimator is used in conjunction with a Wang-Landau sampling scheme [F. Wang and D. P. Landau, Phys. Rev. Lett. 86, 2050 (2001)], we quickly achieve uniform sampling of macrostates (e.g., energies) and systematically refine the calculated density of states. This approach requires only potential energy evaluations, continues to improve the statistical quality of its results as the simulation time is extended, and is applicable to both lattice and continuum systems. We test the algorithm on the Lennard-Jones liquid and demonstrate good statistical convergence properties.Comment: 7 pages, 4 figures. to appear in Journal of Chemical Physic

Spinodal of supercooled polarizable water

We develop a series of molecular dynamics computer simulations of liquid water, performed with a polarizable potential model, to calculate the spinodal line and the curve of maximum density inside the metastable supercooled region. After analysing the structural properties,the liquid spinodal line is followed down to T=210 K. A monotonic decrease is found in the explored region. The curve of maximum density bends on approaching the spinodal line. These results, in agreement with similar studies on non polarizable models of water, are consistent with the existence of a second critical point for water.Comment: 8 pages, 5 figures, 2 tables. To be published in Phys. Re

Saddles in the energy landscape: extensivity and thermodynamic formalism

We formally extend the energy landscape approach for the thermodynamics of liquids to account for saddle points. By considering the extensive nature of macroscopic potential energies, we derive the scaling behavior of saddles with system size, as well as several approximations for the properties of low-order saddles (i.e., those with only a few unstable directions). We then cast the canonical partition function in a saddle-explicit form and develop, for the first time, a rigorous energy landscape approach capable of reproducing trends observed in simulations, in particular the temperature dependence of the energy and fractional order of sampled saddles.Comment: 4 pages, 1 figur

Phase diagram of a polydisperse soft-spheres model for liquids and colloids

The phase diagram of soft spheres with size dispersion has been studied by means of an optimized Monte Carlo algorithm which allows to equilibrate below the kinetic glass transition for all sizes distribution. The system ubiquitously undergoes a first order freezing transition. While for small size dispersion the frozen phase has a crystalline structure, large density inhomogeneities appear in the highly disperse systems. Studying the interplay between the equilibrium phase diagram and the kinetic glass transition, we argue that the experimentally found terminal polydispersity of colloids is a purely kinetic phenomenon.Comment: Version to be published in Physical Review Letter

Liquid-Liquid Phase Transitions for Soft-Core Attractive Potentials

Using event driven molecular dynamics simulations, we study a three dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point and a liquid-liquid critical point. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and give us insight on the mechanisms ruling the dependence of the two critical points on the potential's parameters. The soft core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.Comment: 29 pages, 15 figure

Cooperative Origin of Low-Density Domains in Liquid Water

We study the size of clusters formed by water molecules possessing large enough tetrahedrality with respect to their nearest neighbors. Using Monte Carlo simulation of the SPC/E model of water, together with a geometric analysis based on Voronoi tessellation, we find that regions of lower density than the bulk are formed by accretion of molecules into clusters exceeding a minimum size. Clusters are predominantly linear objects and become less compact as they grow until they reach a size beyond which further accretion is not accompanied by a density decrease. The results suggest that the formation of "ice-like" regions in liquid water is cooperative.Comment: 16 pages, 6 figure

Exact Solution of a Jamming Transition: Closed Equations for a Bootstrap Percolation Problem

Jamming, or dynamical arrest, is a transition at which many particles stop moving in a collective manner. In nature it is brought about by, for example, increasing the packing density, changing the interactions between particles, or otherwise restricting the local motion of the elements of the system. The onset of collectivity occurs because, when one particle is blocked, it may lead to the blocking of a neighbor. That particle may then block one of its neighbors, these effects propagating across some typical domain of size named the dynamical correlation length. When this length diverges, the system becomes immobile. Even where it is finite but large the dynamics is dramatically slowed. Such phenomena lead to glasses, gels, and other very long-lived nonequilibrium solids. The bootstrap percolation models are the simplest examples describing these spatio-temporal correlations. We have been able to solve one such model in two dimensions exactly, exhibiting the precise evolution of the jamming correlations on approach to arrest. We believe that the nature of these correlations and the method we devise to solve the problem are quite general. Both should be of considerable help in further developing this field.Comment: 17 pages, 4 figure

Ion-induced nucleation in polar one-component fluids

We present a Ginzburg-Landau theory of ion-induced nucleation in a gas phase of polar one-component fluids, where a liquid droplet grows with an ion at its center. By calculating the density profile around an ion, we show that the solvation free energy is larger in gas than in liquid at the same temperature on the coexistence curve. This difference much reduces the nucleation barrier in a metastable gas.Comment: 9 pagers, 9 figures, to be published in J. Chem. Phy

Gas-Liquid Nucleation in Two Dimensional System

We study the nucleation of the liquid phase from a supersaturated vapor in two dimensions (2D). Using different Monte Carlo simulation methods, we calculate the free energy barrier for nucleation, the line tension and also investigate the size and shape of the critical nucleus. The study is carried out at an intermediate level of supersaturation(away from the spinodal limit). In 2D, a large cut-off in the truncation of the Lennard-Jones (LJ) potential is required to obtain converged results, whereas low cut-off (say, $2.5\sigma$ is generally sufficient in three dimensional studies, where $\sigma$ is the LJ diameter) leads to a substantial error in the values of line tension, nucleation barrier and characteristics of the critical cluster. It is found that in 2D, the classical nucleation theory (CNT) fails to provide a reliable estimate of the free energy barrier. It underestimates the barrier by as much as 70% at the saturation-ratio S=1.1 (defined as S=P/PC, where PC is the coexistence pressure at reduced temperature $T^{\star}= 0.427$). Interestingly, CNT has been found to overestimate the nucleation free energy barrier in three dimensional (3D)systems near the triple point. In fact, the agreement with CNT is worse in 2D than in 3D. Moreover, the existing theoretical estimate of the line tension overestimates the value significantly.Comment: 24 pages, 8 figure