E. La Nave

Model for reversible colloidal gelation

We report a numerical study, covering a wide range of packing fraction Phi and temperature T, for a system of particles interacting via a square well potential supplemented by an additional constraint on the maximum number n(max) of bonded interactions. We show that, when n(max)<6, the liquid-gas coexistence region shrinks, giving access to regions of low Phi where dynamics can be followed down to low T without an intervening phase separation. We characterize these arrested states at low densities (gel states) in terms of structure and dynamical slowing down, pointing out features which are very different from the standard glassy states observed at high Phi values.

Dynamics and configurational entropy in the Lewis-Wahnström model for supercooled orthoterphenyl.

We study thermodynamic and dynamic properties of a rigid model of the fragile glass-forming liquid orthoterphenyl. This model, introduced by Lewis and Wahnström in 1993, collapses each phenyl ring to a single interaction site; the intermolecular site-site interactions are described by the Lennard-Jones potential whose parameters have been selected to reproduce some bulk properties of the orthoterphenyl molecule. A system of N=343 molecules is considered in a wide range of densities and temperatures, reaching simulation times up to 1 micros. Such long trajectories allow us to equilibrate the system at temperatures below the mode coupling temperature T(c) at which the diffusion constant reaches values of order 10(-10) cm(2)/s and thereby to sample in a significant way the potential energy landscape in the entire temperature range. Working within the inherent structures thermodynamic formalism, we present results for the temperature and density dependence of the number, depth and shape of the basins of the potential energy surface. We evaluate the total entropy of the system by thermodynamic integration from the ideal-noninteracting-gas state and the vibrational entropy approximating the basin free energy with the free energy of 6N-3 harmonic oscillators. We evaluate the configurational part of the entropy as a difference between these two contributions. We study the connection between thermodynamical and dynamical properties of the system. We confirm that the temperature dependence of the configurational entropy and of the diffusion constant, as well as the inverse of the characteristic structural relaxation time, are strongly connected in supercooled states; we demonstrate that this connection is well represented by the Adam-Gibbs relation, stating a linear relation between logD and the quantity 1/TS(c). This relation is found to hold both above and below the critical temperature T(c)-as previously found in the case of silica-supporting the hypothesis that a connection exists between the number of basins and the connectivity properties of the potential energy surface.

Instantaneous Normal Mode Analysis of Supercooled Water

We use the instantaneous normal mode approach to provide a description of the local curvature of the potential energy surface of a model for water. We focus on the region of the phase diagram in which the dynamics may be described by mode-coupling theory. We find that the diffusion constant depends on the fraction of directions in configuration space connecting different local minima, supporting the hypothesis that the dynamics are controlled by the geometric properties of configuration space. Furthermore, we find a relation between the number of basins accessed in equilibrium and the connectivity between them.

Energy landscape of a simple model for strong liquids

We calculate the statistical properties of the energy landscape of a minimal model for strong network-forming liquids. Dynamic and thermodynamic properties of this model can be computed with arbitrary precision even at low temperatures. A degenerate disordered ground state and logarithmic statistics for the local minima energy distribution are the landscape signatures of strong liquid behavior. Differences from fragile liquid properties are attributed to the presence of a discrete energy scale, provided by the particle bonds, and to the intrinsic degeneracy of topologically disordered networks.

Application of statistical physics to understand static and dynamic anomalies in liquid water

We present an overview of recent research applying ideas of statistical mechanics to try to better understand the statics and especially the dynamic puzzles regarding liquid water. We discuss recent molecular dynamics simulations using the Mahoney–Jorgensen transferable intermolecular potential with five points (TIP5P), which is closer to real water than previously-proposed classical pairwise additive potentials. Simulations of the TIP5P model for a wide range of deeply supercooled states, including both positive and negative pressures, reveal (i) the existence of a non-monotonic temperature of maximum density line and a non-reentrant spinodal, (ii) the presence of a low-temperature phase transition. The take-home message for the static aspects is that what seems to “matter” more than previously appreciated is local tetrahedral order, so that liquid water has features in common with SiO2 and P, as well as perhaps Si and C. To better understand dynamic aspects of water, we focus on the role of the number of diffusive directions in the potential energy landscape. What seems to “matter” most is not values of thermodynamic parameters such as temperature T and pressure P, but only the value of a parameter characterizing the potential energy landscape—just as near a critical point what matters is not the values of T and P but rather the values of the correlation length.

Statistical physics and liquid water: "What matters"

We present an overview of recent research applying ideas of statistical mechanics to try to better understand the statics and especially the dynamic puzzles regarding liquid water. The take-home message for the static aspects is that what seems to “matter” more than previously appreciated is local tetrahedral order, so that liquid water has features in common with SiO2 and P, as well as perhaps Si and C. For the dynamic aspects, what may “matter” is the number of diffusive directions in the potential energy landscape.

Test of nonequilibrium thermodynamics in glassy systems: The soft-sphere case

The scaling properties of the soft-sphere potential allow the derivation of an exact expression for the pressure of a frozen liquid, i.e., the pressure corresponding to configurations which are local minima in its multidimensional potential energy landscape. The existence of such a relation offers the unique possibility for testing the recently proposed extension of the liquid free energy to glassy out-of-equilibrium conditions and the associated expression for the temperature of the configurational degrees of freedom. We demonstrate that the nonequilibrium free energy provides an exact description of the soft-sphere pressure in glass states.

Non-Gaussian energy landscape of a simple model for strong network-forming liquids: Accurate evaluation of the configurational entropy

We present a numerical study of the statistical properties of the potential energy landscape of a simple model for strong network-forming liquids. The model is a system of spherical particles interacting through a square-well potential, with an additional constraint that limits the maximum number of bonds Nmax per particle. Extensive simulations have been carried out as a function of temperature, packing fraction, and Nmax. The dynamics of this model are characterized by Arrhenius temperature dependence of the transport coefficients and by nearly exponential relaxation of dynamic correlators, i.e., features defining strong glass-forming liquids. This model has two important features: (i) Landscape basins can be associated with bonding patterns. (ii) The configurational volume of the basin can be evaluated in a formally exact way, and numerically with an arbitrary precision. These features allow us to evaluate the number of different topologies the bonding pattern can adopt. We find that the number of fully bonded configurations, i.e., configurations in which all particles are bonded to Nmax neighbors, is extensive, suggesting that the configurational entropy of the low temperature fluid is finite. We also evaluate the energy dependence of the configurational entropy close to the fully bonded state and show that it follows a logarithmic functional form, different from the quadratic dependence characterizing fragile liquids. We suggest that the presence of a discrete energy scale, provided by the particle bonds, and the intrinsic degeneracy of fully bonded disordered networks differentiates strong from fragile behavior.

Aging and energy landscapes: application to liquids and glasses

Abstract:The equation of state for a liquid in equilibrium, written in the potential energy landscape formalism, is generalized to describe out-of-equilibrium conditions. The hypothesis that during aging the system explores basins associated to equilibrium configurations is the key ingredient in the derivation. Theoretical predictions are successfully compared with data from molecular dynamics simulations of different aging processes, such as temperature and pressure jumps.

Numerical evaluation of the statistical properties of a potential energy landscape

The techniques which allow the numerical evaluation of the statistical properties of the potential energy landscape for models of simple liquids are reviewed and critically discussed. Expressions for the liquid free energy and its vibrational and configurational components are reported. Finally, a possible model for the statistical properties of the landscape, which appears to describe correctly fragile liquids in the region where equilibrium simulations are feasible, is discussed.