Inherent Structure Entropy of Supercooled Liquids

We present a quantitative description of the thermodynamics in a supercooled binary Lennard Jones liquid via the evaluation of the degeneracy of the inherent structures, i.e. of the number of potential energy basins in configuration space. We find that for supercooled states, the contribution of the inherent structures to the free energy of the liquid almost completely decouples from the vibrational contribution. An important byproduct of the presented analysis is the determination of the Kauzmann temperature for the studied system. The resulting quantitative picture of the thermodynamics of the inherent structures offers new suggestions for the description of equilibrium and out-of-equilibrium slow-dynamics in liquids below the Mode-Coupling temperature.Comment: 11 pages of Latex, 3 figure

The relationship between fragility, configurational entropy and the potential energy landscape of glass forming liquids

Glass is a microscopically disordered, solid form of matter that results when a fluid is cooled or compressed in such a fashion that it does not crystallise. Almost all types of materials are capable of glass formation -- polymers, metal alloys, and molten salts, to name a few. Given such diversity, organising principles which systematise data concerning glass formation are invaluable. One such principle is the classification of glass formers according to their fragility\cite{fragility}. Fragility measures the rapidity with which a liquid's properties such as viscosity change as the glassy state is approached. Although the relationship between features of the energy landscape of a glass former, its configurational entropy and fragility have been analysed previously (e. g.,\cite{speedyfr}), an understanding of the origins of fragility in these features is far from being well established. Results for a model liquid, whose fragility depends on its bulk density, are presented in this letter. Analysis of the relationship between fragility and quantitative measures of the energy landscape (the complicated dependence of energy on configuration) reveal that the fragility depends on changes in the vibrational properties of individual energy basins, in addition to the total number of such basins present, and their spread in energy. A thermodynamic expression for fragility is derived, which is in quantitative agreement with {\it kinetic} fragilities obtained from the liquid's diffusivity.Comment: 8 pages, 3 figure

Liquid Limits: The Glass Transition and Liquid-Gas Spinodal Boundaries of Metastable Liquids

The liquid-gas spinodal and the glass transition define ultimate boundaries beyond which substances cannot exist as (stable or metastable) liquids. The relation between these limits is analyzed {\it via} computer simulations of a model liquid. The results obtained indicate that the liquid - gas spinodal and the glass transition lines intersect at a finite temperature, implying a glass - gas mechanical instability locus at low temperatures. The glass transition lines obtained by thermodynamic and dynamic criteria agree very well with each other.Comment: 5 pages, 4 figures, to appear in Phys. Rev. Let

Maximum Valency Lattice Gas Models

We study lattice gas models with the imposition of a constraint on the maximum number of bonds (nearest neighbor interactions) that particles can participate in. The critical parameters, as well as the coexistence region are studied using the mean field approximation and the Bethe-Peierls approximation. We find that the reduction of the number of interactions suppresses the temperature-density region where the liquid and gas phases coexist. We confirm these results from simulations using the histogram reweighting method employing grand Canonical Monte Carlo simulations

Energy landscapes, ideal glasses, and their equation of state

Using the inherent structure formalism originally proposed by Stillinger and Weber [Phys. Rev. A 25, 978 (1982)], we generalize the thermodynamics of an energy landscape that has an ideal glass transition and derive the consequences for its equation of state. In doing so, we identify a separation of configurational and vibrational contributions to the pressure that corresponds with simulation studies performed in the inherent structure formalism. We develop an elementary model of landscapes appropriate to simple liquids which is based on the scaling properties of the soft-sphere potential complemented with a mean-field attraction. The resulting equation of state provides an accurate representation of simulation data for the Lennard-Jones fluid, suggesting the usefulness of a landscape-based formulation of supercooled liquid thermodynamics. Finally, we consider the implications of both the general theory and the model with respect to the so-called Sastry density and the ideal glass transition. Our analysis shows that a quantitative connection can be made between properties of the landscape and a simulation-determined Sastry density, and it emphasizes the distinction between an ideal glass transition and a Kauzmann equal-entropy condition.Comment: 11 pages, 3 figure

Supercooled Water and the Kinetic Glass Transition II: Collective Dynamics

In this article we study in detail the Q-vector dependence of the collective dynamics in simulated deeply supercooled SPC/E water. The evolution of the system has been followed for 250 ns at low T, allowing a clear identification of a two step relaxation process. We present evidence in favor of the use of the mode coupling theory for supercooled liquid as framework for the description of the slow alpha-relaxation dynamics in SPC/E water, notwithstanding the fact that the cage formation in this system is controlled by the formation of an open network of hydrogen bonds as opposed to packing constraints, as in the case of simple liquids.Comment: rev-tex + 9 figure

Freezing by Monte Carlo Phase-Switch

We describe a Monte Carlo procedure which allows sampling of the disjoint configuration spaces associated with crystalline and fluid phases, within a single simulation. The method utilises biased sampling techniques to enhance the probabilities of gateway states (in each phase) which are such that a global switch (to the other phase) can be implemented. Equilibrium freezing-point parameters can be determined directly; statistical uncertainties prescribed transparently; and finite-size effects quantified systematically. The method is potentially quite general; we apply it to the freezing of hard spheres.Comment: 5 pages, 2 figure

Mean field theory of hard sphere glasses and jamming

Hard spheres are ubiquitous in condensed matter: they have been used as models for liquids, crystals, colloidal systems, granular systems, and powders. Packings of hard spheres are of even wider interest, as they are related to important problems in information theory, such as digitalization of signals, error correcting codes, and optimization problems. In three dimensions the densest packing of identical hard spheres has been proven to be the FCC lattice, and it is conjectured that the closest packing is ordered (a regular lattice, e.g, a crystal) in low enough dimension. Still, amorphous packings have attracted a lot of interest, because for polydisperse colloids and granular materials the crystalline state is not obtained in experiments for kinetic reasons. We review here a theory of amorphous packings, and more generally glassy states, of hard spheres that is based on the replica method: this theory gives predictions on the structure and thermodynamics of these states. In dimensions between two and six these predictions can be successfully compared with numerical simulations. We will also discuss the limit of large dimension where an exact solution is possible. Some of the results we present here have been already published, but others are original: in particular we improved the discussion of the large dimension limit and we obtained new results on the correlation function and the contact force distribution in three dimensions. We also try here to clarify the main assumptions that are beyond our theory and in particular the relation between our static computation and the dynamical procedures used to construct amorphous packings.Comment: 59 pages, 25 figures. Final version published on Rev.Mod.Phy

Liquid-Solid Phase Transition of the System with Two particles in a Rectangular Box

We study the statistical properties of two hard spheres in a two dimensional rectangular box. In this system, the relation like Van der Waals equation loop is obtained between the width of the box and the pressure working on side walls. The auto-correlation function of each particle's position is calculated numerically. By this calculation near the critical width, the time at which the correlation become zero gets longer according to the increase of the height of the box. Moreover, fast and slow relaxation processes like $\alpha$ and $\beta$ relaxations observed in supper cooled liquid are observed when the height of the box is sufficiently large. These relaxation processes are discussed with the probability distribution of relative position of two particles.Comment: 6 figure