Potential Energy Landscape Formalism Research Articles

Is the glassy dynamics same in 2D as in 3D? The Adam Gibbs relation test.

It has been recognized of late that even amorphous, glass-forming materials in two dimensions (2D) are affected by Mermin-Wagner-type long wavelength thermal fluctuation, which is inconsequential in three dimensions (3D). We consider the question of whether the effect of spatial dimension on dynamics is only limited to such fluctuations or if the nature of glassy dynamics is intrinsically different in 2D. To address it, we study the relationship between dynamics and thermodynamics using the Adam-Gibbs (AG) relation and the random first order transition (RFOT) theory. Using two model glass-forming liquids, we find that even after removing the effect of long wavelength fluctuations, the AG relation breaks down in two dimensions. Next, we consider the effect of anharmonicity of vibrational entropy-a second factor that affects the thermodynamics but not dynamics. Using the potential energy landscape formalism, we explicitly compute the configurational entropy, both with and without the anharmonic correction. We show that even with both the corrections, the AG relation still breaks down in 2D. The extent of deviation from the AG relation crucially depends on the attractive vs repulsive nature of interparticle interactions, choice of representative timescale (diffusion coefficient vs α-relaxation time), and implies that the RFOT scaling exponents also depend on these factors. Thus, our results suggest that some differences in the nature of glassy dynamics between 2D and 3D remain that are not explained by long wavelength fluctuations.

Potential energy landscape of a flexible water model: Equation of state, configurational entropy, and Adam-Gibbs relationship.

The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid-liquid critical point (LLCP) in the supercooled regime, at (Pc = 150 MPa, Tc = 190 K, and ρc = 1.04 g/cm3) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water and the corresponding location of the LLCP are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian, which allows us to calculate the configurational entropy of the system, Sconf. The Sconf of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam-Gibbs relation, which relates the diffusion coefficient D with Sconf, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high frequencies lead to unphysical (negative) entropy for q-TIP4P/F water when using classical statistical mechanics (yet, the PEL formalism can be applied successfully).

A continuum of amorphous ices between low-density and high-density amorphous ice

Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15−1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated at P ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.

The Harmonic and Gaussian Approximations in the Potential Energy Landscape Formalism for Quantum Liquids.

The potential energy landscape (PEL) formalism has been used in the past to describe the behavior of classical low-temperature liquids and glasses. Here, we extend the PEL formalism to describe the behavior of liquids and glasses that obey quantum mechanics. In particular, we focus on the (i) harmonic and (ii) Gaussian approximations of the PEL, which have been commonly used to describe classical systems, and show how these approximations can be applied to quantum liquids/glasses. Contrary to the case of classical liquids/glasses, the PEL of quantum liquids is temperature-dependent, and hence, the main expressions resulting from approximations (i) and (ii) depend on the nature (classical vs quantum) of the system. The resulting theoretical expressions from the PEL formalism are compared with results from path-integral Monte Carlo (PIMC) simulations of a monatomic model liquid. In the PIMC simulations, every atom of the quantum liquid is represented by a ring-polymer. Our PIMC simulations show that at the local minima of the PEL (inherent structures, or IS), sampled over a wide range of temperatures and volumes, the ring-polymers are collapsed. This considerably facilitates the description of quantum liquids using the PEL formalism. Specifically, the normal modes of the ring-polymer system/quantum liquid at an IS can be calculated analytically if the normal modes of the classical liquid counterpart are known (as obtained, e.g., from classical MC or molecular dynamics simulations of the corresponding atomic liquid).

Potential energy landscape formalism for quantum liquids

The authors extend the potential energy landscape formalism to the case of quantum liquids and how they approach their glass state.

Basis glass states: New insights from the potential energy landscape

Abstract Using the potential energy landscape formalism we show that, in the temperature range in which the dynamics of a glass forming system is thermally activated, there exists a unique set of ‘basis glass states’ each of which is confined to a single metabasin of the energy landscape of a glass forming system. These basis glass states tile the entire configuration space of the system, exhibit only secondary relaxation and are solid-like. Any macroscopic state of the system (whether liquid or glass) can be represented as a superposition of basis glass states and can be described by a probability distribution over these states. During cooling of a liquid from a high temperature, the probability distribution freezes at sufficiently low temperatures describing the process of liquid to glass transition. The time evolution of the probability distribution towards the equilibrium distribution during subsequent aging describes the primary relaxation of a glass.

Glass polymorphism in TIP4P/2005 water: A description based on the potential energy landscape formalism.

The potential energy landscape (PEL) formalism is a statistical mechanical approach to describe supercooled liquids and glasses. Here, we use the PEL formalism to study the pressure-induced transformations between low-density amorphous ice (LDA) and high-density amorphous ice (HDA) using computer simulations of the TIP4P/2005 molecular model of water. We find that the properties of the PEL sampled by the system during the LDA-HDA transformation exhibit anomalous behavior. In particular, at conditions where the change in density during the LDA-HDA transformation is approximately discontinuous, reminiscent of a first-order phase transition, we find that (i) the inherent structure (IS) energy, eIS(V), is a concave function of the volume and (ii) the IS pressure, PIS(V), exhibits a van der Waals-like loop. In addition, the curvature of the PEL at the IS is anomalous, a nonmonotonic function of V. In agreement with previous studies, our work suggests that conditions (i) and (ii) are necessary (but not sufficient) signatures of the PEL for the LDA-HDA transformation to be reminiscent of a first-order phase transition. We also find that one can identify two different regions of the PEL, one associated with LDA and another with HDA. Our computer simulations are performed using a wide range of compression/decompression and cooling rates. In particular, our slowest cooling rate (0.01 K/ns) is within the experimental rates employed in hyperquenching experiments to produce LDA. Interestingly, the LDA-HDA transformation pressure that we obtain at T = 80 K and at different rates extrapolates remarkably well to the corresponding experimental pressure.

State variables for glasses: The case of amorphous ice.

Glasses are out-of-equilibrium systems whose state cannot be uniquely defined by the usual set of equilibrium state variables. Here, we seek to identify an expanded set of variables that uniquely define the state of a glass. The potential energy landscape (PEL) formalism is a useful approach within statistical mechanics to describe supercooled liquids and glasses. We use the PEL formalism and computer simulations to study the transformations between low-density amorphous ice (LDA) and high-density amorphous ice (HDA). We employ the ST2 water model, which exhibits an abrupt first-order-like phase transition from LDA to HDA, similar to that observed in experiments. We prepare a number of distinct samples of both LDA and HDA that have completely different preparation histories. We then study the evolution of these LDA and HDA samples during compression and decompression at temperatures sufficiently low that annealing is absent and also during heating. We find that the evolution of each glass sample, during compression/decompression or heating, is uniquely determined by six macroscopic properties of the initial glass sample. These six quantities consist of three conventional thermodynamic state variables, the number of molecules N, the system volume V, and the temperature T, as well as three properties of the PEL, the inherent structure (IS) energy EIS, the IS pressure PIS, and the average curvature of the PEL at the IS SIS. In other words, (N,V,T,EIS,PIS,SIS) are state variables that define the glass state in the case of amorphous ice. An interpretation of our results in terms of the PEL formalism is provided. Since the behavior of water in the glassy state is more complex than for most substances, our results suggest that these six state variables may be applicable to amorphous solids in general and that there may be situations in which fewer than six variables would be sufficient to define the state of a glass.

Potential energy landscape of the apparent first-order phase transition between low-density and high-density amorphous ice.

The potential energy landscape (PEL) formalism is a valuable approach within statistical mechanics to describe supercooled liquids and glasses. Here we use the PEL formalism and computer simulations to study the pressure-induced transformations between low-density amorphous ice (LDA) and high-density amorphous ice (HDA) at different temperatures. We employ the ST2 water model for which the LDA-HDA transformations are remarkably sharp, similar to what is observed in experiments, and reminiscent of a first-order phase transition. Our results are consistent with the view that LDA and HDA configurations are associated with two distinct regions (megabasins) of the PEL that are separated by a potential energy barrier. At higher temperature, we find that low-density liquid (LDL) configurations are located in the same megabasin as LDA, and that high-density liquid (HDL) configurations are located in the same megabasin as HDA. We show that the pressure-induced LDL-HDL and LDA-HDA transformations occur along paths that interconnect these two megabasins, but that the path followed by the liquid is different from the path followed by the amorphous solid. At higher pressure, we also study the liquid-to-ice-VII first-order phase transition, and find that the behavior of the PEL properties across this transition is qualitatively similar to the changes found during the LDA-HDA transformation. This similarity supports the interpretation that the LDA-HDA transformation is a first-order phase transition between out-of-equilibrium states. Finally, we compare the PEL properties explored during the LDA-HDA transformations in ST2 water with those reported previously for SPC/E water, for which the LDA-HDA transformations are rather smooth. This comparison illuminates the previous work showing that, at accessible computer times scales, a liquid-liquid phase transition occurs in the case of ST2 water, but not for SPC/E water.

Comparing landscape calculations with calorimetric data on ortho-terphenyl, and the question of the configurational fraction of the excess entropy

Mossa et al. [Phys. Rev. E 65, 041205 (2002)] have calculated the total and configurational entropies of supercooled ortho-terphenyl liquid using the potential-energy landscape formalism and a simplified model of the intermolecular potential. I show here that the agreement of their calculated configurational entropy with the experimental data depends on what is assumed about the configurational fraction of the excess entropy and its temperature dependence. In particular, if the configurational fraction is taken as 0.70 and independent of temperature the agreement is excellent; if a marked temperature dependence of that fraction inferred from calorimetric data is assumed the agreement is only fair at best. This marked temperature dependence of the configurational fraction also implies some implausible behavior of contributions to the excess entropy at the Kauzmann temperature, but no obvious reason for disregarding it presents itself.

Liquid stability in a model for ortho-terphenyl.

We report an extensive study of the phase diagram of a simple model for ortho-terphenyl, focusing on the limits of stability of the liquid state. Reported data extend previous studies of the same model to both lower and higher densities and to higher temperatures. We estimate the location of the homogeneous liquid-gas nucleation line and of the spinodal locus. Within the potential energy landscape formalism, we calculate the distributions of depth, number, and shape of the potential energy minima and show that the statistical properties of the landscape are consistent with a Gaussian distribution of minima over a wide range of volumes. We report the volume dependence of the parameters entering in the Gaussian distribution (amplitude, average energy, variance). We finally evaluate the locus where the configurational entropy vanishes, the so-called Kauzmann line, and discuss the relative location of the spinodal and Kauzmann loci.

Aging and energy landscapes: application to liquids and glasses

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.