Structure, dynamics, and thermodynamics of a family of potentials with tunable softness

The role of the diffusion in the predictions of the classical nucleation theory for quasi-real systems differ in dipole moment value

In this paper, we examine the crystallization tendency for two quasi-real systems, which differ exclusively in the dipole moment's value. The main advantage of the studied system is the fact that despite that their structures are entirely identical, they exhibit different physical properties. Hence, the results obtained for one of the proposed model systems cannot be scaled to reproduce the results for another corresponding system, as it can be done for simple model systems, where structural differences are modeled by the different parameters of the intermolecular interactions. Our results show that both examined systems exhibit similar stability behavior below the melting temperature. This finding is contrary to the predictions of the classical nucleation theory, which suggests a significantly higher crystallization tendency for a more polar system. Our studies indicate that the noted discrepancies are caused by the kinetic aspect of the classical nucleation theory, which overestimates the role of diffusion in the nucleation process.

Fragile-to-strong crossover, growing length scales, and dynamic heterogeneity in Wigner glasses

Colloidal particles, which are ubiquitous, have become ideal testing grounds for the structural glass transition theories. In these systems glassy behavior arises as the density of the particles is increased. Thus, soft colloidal particles with varying degree of softness capture diverse glass-forming properties, observed normally in molecular glasses. Brownian dynamics simulations for a binary mixture of micron-sized charged colloidal suspensions show that tuning the softness of the interaction potential, achievable by changing the monovalent salt concentration results in a continuous transition from fragile to strong behavior. Remarkably, this is found in a system where the well characterized interaction potential between the colloidal particles is isotropic. We also show that the predictions of the random first-order transition (RFOT) theory quantitatively describes the universal features such as the growing correlation length, ξ∼(ϕ_{K}/ϕ-1)^{-ν} with ν=2/3 where ϕ_{K}, the analog of the Kauzmann temperature, depends on the salt concentration. As anticipated by the RFOT predictions, we establish a causal relationship between the growing correlation length and a steep increase in the relaxation time and dynamic heterogeneity as the system is compressed. The broad range of fragility observed in Wigner glasses is used to draw analogies with molecular and polymer glasses. The large variations in the fragility are normally found only when the temperature dependence of the viscosity is examined for a large class of diverse glass-forming materials. In sharp contrast, this is vividly illustrated in a single system that can be experimentally probed. Our work also shows that the RFOT predictions are accurate in describing the dynamics over the entire density range, regardless of the fragility of the glasses.

The role of the dipole moment orientations in the crystallization tendency of the van der Waals liquids - molecular dynamics simulations

Computer simulations of model systems play a remarkable role in the contemporary studies of structural, dynamic and thermodynamic properties of supercooled liquids. However, the commonly employed model systems, i.e., simple-liquids, do not reflect the internal features of the real molecules, e.g., structural anisotropy and spatial distribution of charges, which might be crucial for the behavior of real materials. In this paper, we use the new model molecules of simple but anisotropic structure, to studies the effect of dipole moment orientation on the crystallization tendency. Our results indicate that proper orientation of the dipole moment could totally change the stability behavior of the system. Consequently, the exchange of a single atom within the molecule causing the change of dipole moment orientation might be crucial for controlling the crystallization tendency. Moreover, employing the classical nucleation theory, we explain the reason for this behavior.

Universal localization transition accompanying glass formation: insights from efficient molecular dynamics simulations of diverse supercooled liquids

The origin of the precipitous dynamic arrest known as the glass transition is a grand open question of soft condensed matter physics. It has long been suspected that this transition is driven by an onset of particle localization and associated emergence of a glassy modulus. However, progress towards an accepted understanding of glass formation has been impeded by an inability to obtain data sufficient in chemical diversity, relaxation timescales, and spatial and temporal resolution to validate or falsify proposed theories for its physics. Here we first describe a strategy enabling facile high-throughput simulation of glass-forming liquids to nearly unprecedented relaxation times. We then perform simulations of 51 glass-forming liquids, spanning polymers, small organic molecules, inorganics, and metallic glass-formers, with longest relaxation times exceeding one microsecond. Results identify a universal particle-localization transition accompanying glass formation across all classes of glass-forming liquid. The onset temperature of non-Arrhenius dynamics is found to serve as a normalizing condition leading to a master collapse of localization data. This transition exhibits a non-universal relationship with dynamic arrest, suggesting that the nonuniversality of supercooled liquid dynamics enters via the dependence of relaxation times on local cage scale. These results suggest that a universal particle-localization transition may underpin the glass transition, and they emphasize the potential for recent theoretical developments connecting relaxation to localization and emergent elasticity to finally explain the origin of this phenomenon. More broadly, the capacity for high-throughput prediction of glass formation behavior may open the door to computational inverse design of glass-forming materials.

Liquid state isomorphism, Rosenfeld-Tarazona temperature scaling, and Riemannian thermodynamic geometry

Aspects of isomorph theory, Rosenfeld-Tarazona temperature scaling, and thermodynamic geometry are comparatively discussed on the basis of the Lennard-Jones potential. The first two approaches approximate the high-density fluid state well when the repulsive interparticle interactions become dominant, which is typically the case close to the freezing line. However, previous studies of Rosenfeld-Tarazona scaling for the isochoric heat capacity and its relation to isomorph theory reveal deviations for the temperature dependence. It turns out that a definition of a state region in which repulsive interactions dominate is required for achieving consistent results. The Riemannian thermodynamic scalar curvature R allows for such a classification, indicating predominantly repulsive interactions by R>0. An analysis of the isomorphic character of the freezing line and the validity of Rosenfeld-Tarazona temperature scaling show that these approaches are consistent only in a small state region.

Effects of the bond polarity on the structural and dynamical properties of silica-like liquids

Silica is a network-forming liquid that shares many properties with water due to its tetrahedral structure. It undergoes a transition from a fragile to a strong liquid as the temperature is decreased, which is accompanied by a structural change to lower density and higher tetrahedral order. In order to disentangle the effects of Coulomb and van der Waals interactions on the structure and dynamics of liquid silica, we modify the bond polarity by changing the partial charges assigned to each atom. Using molecular dynamics simulations, we show that density, tetrahedral order, and structural relaxation times decrease when reducing bond polarity. Moreover, we find that the density maximum and the fragile-to-strong transition move to lower temperatures until they eventually vanish when the partial charges are decreased below approximately 75% of their regular value. Irrespective of whether strong or fragile behavior exists, structural relaxation is governed by hopping motion at sufficiently low temperatures. As long as there is a strong regime, the energy barrier associated with strong dynamics decreases with decreasing partial charges, but the dependence on the bond polarity differs from that of the activation energy in the Arrhenius regime at high temperatures. We show that the fragile-to-strong transition is associated with structural changes occurring between the first and second coordination shells that lead to a decrease in density and an increase in tetrahedral order. In particular, independent of the value of the partial charges, the distribution of the local structures is the same at this dynamic crossover, but we find no evidence that the effect occurs upon crossing the Widom line. In the fragile regime at intermediate temperatures, the relaxation times are well described by a previously proposed model which decomposes the apparent activation energy into a constant single-particle contribution and a temperature-dependent collective contribution. However, our results for silica-like melts do not obey several common relations of the model parameters reported for molecular glass formers.

Distributions of pore sizes and atomic densities in binary mixtures revealed by molecular dynamics simulations

We report on the results of a molecular dynamics simulation study of porous glassy media, formed in the process of isochoric rapid quenching from a high-temperature liquid state. The transition to a porous solid occurs due to the concurrent processes of phase separation and material solidification. The study is focused on topographies of the model porous structures and their dependence on temperature and average density. To quantify the pore-size distributions, we put forth a scaling relation that provides a satisfactory data collapse in systems with high porosity. We also find that the local density of the solid domains in the porous structures is broadly distributed, and, with increasing average density, a distinct peak in the local density distribution is displaced toward higher values.

Strength of the repulsive part of the interatomic potential determines fragility in metallic liquids

The dynamical behaviour of liquids is frequently characterized by the fragility, which can be defined from the temperature dependence of the shear viscosity, η (ref. ). For a strong liquid, the activation energy for η changes little with cooling towards the glass transition temperature, T_g. The change is much greater in fragile liquids, with the activation energy becoming very large near T_g. While fragility is widely recognized as an important concept-believed, for example, to play an important role in glass formation-the microscopic origin of fragility is poorly understood. Here, we present new experimental evidence showing that fragility reflects the strength of the repulsive part of the interatomic potential, which can be determined from the steepness of the pair distribution function near the hard-sphere cutoff. On the basis of an analysis of scattering data from ten different metallic alloy liquids, we show that stronger liquids have steeper repulsive potentials.

Kinetic and structural fragility-a correlation between structures and dynamics in metallic liquids and glasses

The liquid phase remains poorly understood. In many cases, the densities of liquids and their crystallized solid phases are similar, but since they are amorphous they lack the spatial order of the solid. Their dynamical properties change remarkably over a very small temperature range. At high temperatures, near their melting temperature, liquids flow easily under shear. However, only a few hundred degrees lower flow effectively ceases, as the liquid transforms into a solid-like glass. This temperature-dependent dynamical behavior is frequently characterized by the concept of kinetic fragility (or, generally, simply fragility). Fragility is believed to be an important quantity in glass formation, making it of significant practical interest. The microscopic origin of fragility remains unclear, however, making it also of fundamental interest. It is widely (although not uniformly) believed that the dynamical behavior is linked to the atomic structure of the liquid, yet experimental studies show that although the viscosity changes by orders of magnitude with temperature, the structural change is barely perceptible. In this article the concept of fragility is discussed, building to a discussion of recent results in metallic glass-forming liquids that demonstrate the presumed connection between structural and dynamical changes. In particular, it becomes possible to define a structural fragility parameter that can be linked with the kinetic fragility.

Glass-Forming Tendency of Molecular Liquids and the Strength of the Intermolecular Attractions

When we cool down a liquid below the melting temperature, it can either crystallize or become supercooled, and then form a disordered solid called glass. Understanding what makes a liquid to crystallize readily in one case and form a stable glass in another is a fundamental problem in science and technology. Here we show that the crystallization/glass-forming tendencies of the molecular liquids might be correlated with the strength of the intermolecular attractions, as determined from the combined experimental and computer simulation studies. We use van der Waals bonded propylene carbonate and its less polar structural analog 3-methyl-cyclopentanone to show that the enhancement of the dipole-dipole forces brings about the better glass-forming ability of the sample when cooling from the melt. Our finding was rationalized by the mismatch between the optimal temperature range for the nucleation and crystal growth, as obtained for a modeled Lennard-Jones system with explicitly enhanced or weakened attractive part of the intermolecular 6–12 potential.

Bending and Gaussian rigidities of confined soft spheres from second-order virial series

We use virial series to study the equilibrium properties of confined soft-spheres fluids interacting through the inverse-power potentials. The confinement is induced by hard walls with planar, spherical, and cylindrical shapes. We evaluate analytically the coefficients of order two in density of the wall-fluid surface tension γ and analyze the curvature contributions to the free energy. Emphasis is in bending and Gaussian rigidities, which are found analytically at order two in density. Their contribution to γ(R) and the accuracy of different truncation procedures to the low curvature expansion are discussed. Finally, several universal relations that apply to low-density fluids are analyzed.

Hidden scale invariance in condensed matter

Recent developments show that many liquids and solids have an approximate "hidden" scale invariance that implies the existence of lines in the thermodynamic phase diagram, so-called isomorphs, along which structure and dynamics in properly reduced units are invariant to a good approximation. This means that the phase diagram becomes effectively one-dimensional with regard to several physical properties. Liquids and solids with isomorphs include most or all van der Waals bonded systems and metals, as well as weakly ionic or dipolar systems. On the other hand, systems with directional bonding (hydrogen bonds or covalent bonds) or strong Coulomb forces generally do not exhibit hidden scale invariance. The article reviews the theory behind this picture of condensed matter and the evidence for it coming from computer simulations and experiments.

Relaxation processes in liquids: variations on a theme by Stokes and Einstein

We investigate numerically the temperature and density dependence of the Stokes-Einstein ratio, Dη∕T, and of two commonly-used variants thereof, Dτ and Dτ∕T, where D is a diffusivity, η the shear viscosity, and τ a structural relaxation time. We consider a family of atomic binary mixtures with systematically-softened repulsive interactions, and the Lewis-Wahnström model of ortho-terphenyl (OTP). The three quantities grow significantly as the temperature decreases in the supercooled regime, a well-known phenomenon. At higher temperatures, Dτ exhibits negative violations of Stokes-Einstein behavior, i.e., decrease upon cooling, for the atomic systems, though not for OTP. We consider two choices for the relaxation time, one based on the decay of the self-intermediate scattering function, and the other on the integral of the stress autocorrelation function. The instantaneous shear modulus exhibits appreciable temperature dependence for the two classes of systems investigated here. Our results suggest that commonly-invoked assumptions, such as τ ∼ η and τ ∼ η∕T, should be critically evaluated across a wide spectrum of systems and thermodynamic conditions. We find the Stokes-Einstein ratio, Dη∕T, to be constant across a broad range of temperatures and densities for the two classes of systems investigated here.

The fragility of liquids and colloids and its relation to the softness of the potential

A parameter that is often used to characterize the dynamics of supercooled liquids is the dynamic fragility, however it is still debated how the fragility is related to other physical properties. Recent experimental data on colloidal systems have found that fragility decreases with increasing softness of the intermolecular potential. This result is in apparent disagreement with recent molecular dynamics simulations reporting the opposite behavior. Herein, using the thermodynamical scaling exponent γ as a measure of the steepness of the potential we show how these different results can be reconciled and also agree with previous results obtained for the dynamics of supercooled liquids at high pressures.

Superfragile glassy dynamics of a one-component system with isotropic potential: competition of diffusion and frustration

We investigate glassy dynamical properties of one-component three-dimensional system of particles interacting via pair repulsive potential by the molecular dynamic simulation in the wide region of densities. The glass state is superfragile and it has high glass-forming ability. The glass transition temperature T(g) has a pronounced minimum at densities where the frustration is maximal.

Dependence of the fragility of a glass former on the softness of interparticle interactions

We study the influence of the softness of the interparticle interactions on the fragility of a glass former by considering three model binary mixture glass formers. The interaction potential between particles is a modified Lennard-Jones type potential, with the repulsive part of the potential varying with an inverse power q of the interparticle distance, and the attractive part varying with an inverse power p. We consider the combinations (12,11) (model I), (12,6) (model II), and (8,5) (model III) for (q,p) such that the interaction potential becomes softer from model I to III. We evaluate the kinetic fragilities from the temperature variation of diffusion coefficients and relaxation times, and a thermodynamic fragility from the temperature variation of the configurational entropy. We find that the kinetic fragility increases with increasing softness of the potential, consistent with previous results for these model systems, but at variance with the thermodynamic fragility, which decreases with increasing softness of the interactions, as well as expectations from earlier results. We rationalize our results by considering the full form of the Adam-Gibbs relation, which requires, in addition to the temperature dependence of the configurational entropy, knowledge of the high temperature activation energies in order to determine fragility. We show that consideration of the scaling of the high temperature activation energy with the liquid density, analyzed in recent studies, provides a partial rationalization of the observed behavior.

Two liquid states of matter: a dynamic line on a phase diagram

It is generally agreed that the supercritical region of a liquid consists of one single state (supercritical fluid). On the other hand, we show here that liquids in this region exist in two qualitatively different states: "rigid" and "nonrigid" liquids. Rigid to nonrigid transition corresponds to the condition τ≈τ(0), where τ is the liquid relaxation time and τ(0) is the minimal period of transverse quasiharmonic waves. This condition defines a new dynamic crossover line on the phase diagram and corresponds to the loss of shear stiffness of a liquid at all available frequencies and, consequently, to the qualitative change in many important liquid properties. We analyze this line theoretically as well as in real and model fluids and show that the transition corresponds to the disappearance of high-frequency sound, to the disappearance of roton minima, qualitative changes in the temperature dependencies of sound velocity, diffusion, viscous flow, and thermal conductivity, an increase in particle thermal speed to half the speed of sound, and a reduction in the constant volume specific heat to 2k(B) per particle. In contrast to the Widom line that exists near the critical point only, the new dynamic line is universal: It separates two liquid states at arbitrarily high pressure and temperature and exists in systems where liquid-gas transition and the critical point are absent altogether. We propose to call the new dynamic line on the phase diagram "Frenkel line".