Communication: Correlation of the instantaneous and the intermediate-time elasticity with the structural relaxation in glassforming systems

Testing Theories of the Glass Transition with the Same Liquid but Many Kinetic Rules

We study the glass transition by exploring a broad class of kinetic rules that can significantly modify the normal dynamics of supercooled liquids while maintaining thermal equilibrium. Beyond the usual dynamics of liquids, this class includes dynamics in which a fraction (1-f_{R}) of the particles can perform pairwise exchange or "swap moves," while a fraction f_{P} of the particles can move only along restricted directions. We find that (i) the location of the glass transition varies greatly but smoothly as f_{P} and f_{R} change and (ii) it is governed by a linear combination of f_{P} and f_{R}. (iii) Dynamical heterogeneities (DHs) are not governed by the static structure of the material; their magnitude correlates instead with the relaxation time. (iv) We show that a recent theory for temporal growth of DHs based on thermal avalanches holds quantitatively throughout the (f_{R},f_{P}) diagram. These observations are negative items for some existing theories of the glass transition, particularly those reliant on growing thermodynamic order or locally favored structure, and open new avenues to test other approaches, as we illustrate.

Confinement effect of inter-arm interactions on glass formation in star polymer melts

We utilized molecular dynamic simulation to investigate the glass formation of star polymer melts in which the topological complexity is varied by altering the number of star arms (f). Emphasis was placed on how the "confinement effect" of repulsive inter-arm interactions within star polymers influences the thermodynamics and dynamics of star polymer melts. All the characteristic temperatures of glass formation were found to progressively increase with increasing f, but unexpectedly the fragility parameter KVFT was found to decrease with increasing f. As previously observed, stars having more than 5 or 6 arms adopt an average particle-like structure that is more contracted relative to the linear polymer size having the same mass and exhibit a strong tendency for intermolecular and intramolecular segregation. We systematically analyzed how varying f alters collective particle motion, dynamic heterogeneity, the decoupling exponent ζ phenomenologically linking the slow β- and α-relaxation times, and the thermodynamic scaling index γt. Consistent with our hypothesis that the segmental dynamics of many-arm star melts and thin supported polymer films should exhibit similar trends arising from the common feature of high local segmental confinement, we found that ζ increases considerably with increasing f, as found in supported polymer films with decreasing thickness. Furthermore, increasing f led to greatly enhanced elastic heterogeneity, and this phenomenon correlates strongly with changes in ζ and γt. Our observations should be helpful in building a more rational theoretical framework for understanding how molecular topology and geometrical confinement influence the dynamics of glass-forming materials more broadly.

Temperature Dependence of Structural Relaxation in Glass-Forming Liquids and Polymers

Understanding the microscopic mechanism of the transition of glass remains one of the most challenging topics in Condensed Matter Physics. What controls the sharp slowing down of molecular motion upon approaching the glass transition temperature Tg, whether there is an underlying thermodynamic transition at some finite temperature below Tg, what the role of cooperativity and heterogeneity are, and many other questions continue to be topics of active discussions. This review focuses on the mechanisms that control the steepness of the temperature dependence of structural relaxation (fragility) in glass-forming liquids. We present a brief overview of the basic theoretical models and their experimental tests, analyzing their predictions for fragility and emphasizing the successes and failures of the models. Special attention is focused on the connection of fast dynamics on picosecond time scales to the behavior of structural relaxation on much longer time scales. A separate section discusses the specific case of polymeric glass-forming liquids, which usually have extremely high fragility. We emphasize the apparent difference between the glass transitions in polymers and small molecules. We also discuss the possible role of quantum effects in the glass transition of light molecules and highlight the recent discovery of the unusually low fragility of water. At the end, we formulate the major challenges and questions remaining in this field.

Fast evaluation technique for the shear viscosity and ionic conductivity of electrolyte solutions

With the growing need to obtain ideal materials for various applications, there is an increasing interest in computational methods to rapidly and accurately search for materials. Molecular dynamics simulation is one of the successful methods used to investigate liquid electrolytes with high transport properties applied in lithium-ion batteries. However, further reduction in computational cost is required to find a novel material with the desired properties from a large number of combinations. In this study, we demonstrate an effective fast evaluation technique for shear viscosity and ionic conductivity by molecular dynamics simulation for an exhaustive search of electrolyte materials with high transport properties. The proposed model was combined with a short-time correlation function of the stress tensor and empirical relationships to address the issues of inefficient and uncertain evaluation by conventional molecular dynamics methods. Because we focus on liquid electrolytes consisting of organic solvents and lithium salts, our model requires dissociation ratio and effective diffusion size of lithium salts. Our method is applied to search for the compositional combinations of electrolytes with superior transport properties even at low temperatures. These results correlate well with experimental results.

The initiation of shear band formation in deformed metallic glasses from soft localized domains

It has long been thought that shear band (SB) formation in amorphous solids initiates from relatively "soft" regions in the material in which large-scale non-affine deformations become localized. The test of this hypothesis requires an effective means of identifying "soft" regions and their evolution as the material is deformed to varying degrees, where the metric of "softness" must also account for the effect of temperature on local material stiffness. We show that the mean square atomic displacement on a caging timescale ⟨u²⟩, the "Debye-Waller factor," provides a useful method for estimating the shear modulus of the entire material and, by extension, the material stiffness at an atomic scale. Based on this "softness" metrology, we observe that SB formation indeed occurs through the strain-induced formation of localized soft regions in our deformed metallic glass free-standing films. Unexpectedly, the critical strain condition for SB formation occurs when the softness (⟨u²⟩) distribution within the emerging soft regions approaches that of the interfacial region in its undeformed state, initiating an instability with similarities to the transition to turbulence. Correspondingly, no SBs arise when the material is so thin that the entire material can be approximately described as being "interfacial" in nature. We also quantify relaxation in the glass and the nature and origin of highly non-Gaussian particle displacements in the dynamically heterogeneous SB regions at times longer than the caging time.

Stress hyperuniformity and transient oscillatory-exponential correlation decay as signatures of strength vs fragility in glasses

We examine and compare the local stress autocorrelation in the inherent states of a fragile and a strong glass: the Kob-Andersen (KA) binary mixture and the Beest-Kramer-Santen model of silica. For both systems, local (domain-averaged) stress fluctuations asymptotically reach the normal inverse-volume decay in the large domain limit; accordingly, the real-space stress autocorrelation presents long-range power law tails. However, in the case of silica, local stress fluctuations display a high degree of hyperuniformity, i.e., their asymptotic (normal) decay is disproportionately smaller than their bond level amplitude. This property causes the asymptotic power law tails of the real-space stress autocorrelation to be swamped, up to very large distances (several nanometers), by an intermediate oscillatory-exponential decay regime. Similar contributions exist in the KA stress autocorrelation, but they never can be considered as dominating the power law decay and fully disappear when stress is coarse-grained beyond one interatomic distance. Our observations document that the relevance of power-law stress correlation may constitute a key discriminating feature between strong and fragile glasses. Meanwhile, they highlight that the notion of local stress in atomistic systems involves by necessity a choice of observation (coarse-graining) scale, the relevant value of which depends, in principle, on both the model and the phenomenon studied.

Understanding the scaling of boson peak through insensitivity of elastic heterogeneity to bending rigidity in polymer glasses

Amorphous materials exhibit peculiar mechanical and vibrational properties, including non-affine elastic responses and excess vibrational states, i.e., the so-called boson peak (BP). For polymer glasses, these properties are considered to be affected by the bending rigidity of the constituent polymer chains. In our recent work [Tomoshige,et al2019,Sci. Rep.919514], we have revealed simple relationships between the variations of vibrational properties and the global elastic properties: the response of the BP scales only with that of the global shear modulus. This observation suggests that the spatial heterogeneity of the local shear modulus distribution is insensitive to changes in the bending rigidity. Here, we demonstrate the insensitivity of elastic heterogeneity by directly measuring the local shear modulus distribution. We also study transverse sound wave propagation, which is also shown to scale only with the global shear modulus. Through these analyses, we conclude that the bending rigidity does not alter the spatial heterogeneity of the local shear modulus distribution, which yields vibrational and acoustic properties that are controlled solely by the global shear modulus of a polymer glass.

Predictive relation for the α-relaxation time of a coarse-grained polymer melt under steady shear

We examine the influence of steady shear on structural relaxation in a simulated coarse-grained unentangled polymer melt over a wide range of temperature and shear rates. Shear is found to progressively suppress the α-relaxation process observed in the intermediate scattering function, leading ultimately to a purely inertially dominated β-relaxation at high shear rates, a trend similar to increasing temperature. On the basis of a scaling argument emphasizing dynamic heterogeneity in cooled liquids and its alteration under material deformation, we deduce and validate a parameter-free scaling relation for both the structural relaxation time τα from the intermediate scattering function and the "stretching exponent" β quantifying the extent of dynamic heterogeneity over the entire range of temperatures and shear rates that we can simulate.

Metallic glass-formers in 2D exhibit the same scaling as in 3D between vibrational dynamics and structural relaxation

Glass-forming systems approaching their glass transition exhibit universal correlations between picosecond vibrational dynamics and long-time structural relaxation, which can be described by the same master curve in the bulk or confined conditions. In this work, we study at a fundamental level the effects of the reduction of spatial dimensionality on this phenomenon. We perform molecular dynamics simulations of a metallic glass-formers in two dimensions (2D). We show that in the supercooled regime particle localization in the cage and structural relaxation are blurred by long-wavelength fluctuations specific to low-dimensional systems. Once these effects are properly removed, we demonstrate that the fast dynamics and slow relaxation comply, without any adjustment, with same scaling between the structural relaxation time and the Debye-Waller factor, originally observed in three-dimensions (3D).

Boson peak, elasticity, and glass transition temperature in polymer glasses: Effects of the rigidity of chain bending

The excess low-frequency vibrational spectrum, called boson peak, and non-affine elastic response are the most important particularities of glasses. Herein, the vibrational and mechanical properties of polymeric glasses are examined by using coarse-grained molecular dynamics simulations, with particular attention to the effects of the bending rigidity of the polymer chains. As the rigidity increases, the system undergoes a glass transition at a higher temperature (under a constant pressure), which decreases the density of the glass phase. The elastic moduli, which are controlled by the decrease of the density and the increase of the rigidity, show a non-monotonic dependence on the rigidity of the polymer chain that arises from the non-affine component. Moreover, a clear boson peak is observed in the vibrational density of states, which depends on the macroscopic shear modulus G. In particular, the boson peak frequency ω_BP is proportional to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{G}$$\end{document}G. These results provide a positive correlation between the boson peak, shear elasticity, and the glass transition temperature.

Vibrational scaling of the heterogeneous dynamics detected by mutual information

The correlations detected by the mutual information in the propensities of a molecular viscous liquid are studied by molecular-dynamics simulations. Dynamic heterogeneity is evidenced and two particle fractions with different mobility and relaxation identified. The two fractions exhibit the scaling of their relaxation in terms of the rattling amplitude of the particle trapped in the cage of the first neighbours 〈u²〉 . The scaling master curve does not differ from the one found for bulk systems, thus confirming identical results previously reported in other systems with strong dynamic heterogeneity as thin molecular films. The excitation of planar and globular structures at short and long times with respect to structural relaxation, respectively, is revealed. Some of the globular structures are different from the ones evidenced in atomic mixtures. States with equal 〈u²〉 are found to have identical time dependence of several quantities, referring to both bulk and the two fractions with heterogeneous dynamics, at least up to the structural relaxation time [Formula: see text].

Fast Vibrational Modes and Slow Heterogeneous Dynamics in Polymers and Viscous Liquids

Many systems, including polymers and molecular liquids, when adequately cooled and/or compressed, solidify into a disordered solid, i.e., a glass. The transition is not abrupt, featuring progressive decrease of the microscopic mobility and huge slowing down of the relaxation. A distinctive aspect of glass-forming materials is the microscopic dynamical heterogeneity (DH), i.e., the presence of regions with almost immobile particles coexisting with others where highly mobile ones are located. Following the first compelling evidence of a strong correlation between vibrational dynamics and ultraslow relaxation, we posed the question if the vibrational dynamics encodes predictive information on DH. Here, we review our results, drawn from molecular-dynamics numerical simulation of polymeric and molecular glass-formers, with a special focus on both the breakdown of the Stokes-Einstein relation between diffusion and viscosity, and the size of the regions with correlated displacements.

Mutual information does not detect growing correlations in the propensity of a model molecular liquid

The dynamical spatial correlations detected by the mutual information (MI) in the isoconfigurational particle displacements of a monodisperse molecular viscous liquid are studied via molecular-dynamics simulations by changing considerably both the molecular mobility and the degree of dynamical heterogeneity. Different from atomic liquids, the MI correlation length does not grow on approaching the glass transition by considering the liquid both in full detail as a collection of monomers and as a coarse-grained ensemble of molecular centers of mass. In the detailed picture, it is found that: (i) the MI correlations between monomers are largely due to inter-molecular correlations, (ii) the MI length scale is numerically identical, within the errors, to the correlation length scale of the displacement direction, as drawn by conventional correlation functions. The time evolution of the MI spatial correlations complies with the scaling between the fast vibrational dynamics and the long-time relaxation. Our findings suggest that the characteristics of the MI length scale are markedly system-dependent and not obviously related to dynamical heterogeneity.

Mechanical heterogeneity in ionic liquids

Molecular dynamics (MD) simulations of five ionic liquids based on 1-alkyl-3-methylimidazolium cations, [C _n C₁im]⁺, have been performed in order to calculate high-frequency elastic moduli and to evaluate heterogeneity of local elastic moduli. The MD simulations of [C _n C₁im][NO₃], n = 2, 4, 6, and 8, assessed the effect of domain segregation when the alkyl chain length increases, and [C₈C₁im][PF₆] assessed the effect of strength of anion-cation interaction. Dispersion curves of excitation energies of longitudinal and transverse acoustic, LA and TA, modes were obtained from time correlation functions of mass currents at different wavevectors. High-frequency sound velocity of LA modes depends on the alkyl chain length, but sound velocity for TA modes does not. High-frequency bulk and shear moduli, K _∞ and G _∞ , depend on the alkyl chain length because of a density effect. Both K _∞ and G _∞ are strongly dependent on the anion. The calculation of local bulk and shear moduli was accomplished by performing bulk and shear deformations of the systems cooled to 0 K. The simulations showed a clear connection between structural and elastic modulus heterogeneities. The development of nano-heterogeneous structure with increasing length of the alkyl chain in [C _n C₁im][NO₃] implies lower values for local bulk and shear moduli in the non-polar domains. The mean value and the standard deviations of distributions of local elastic moduli decrease when [NO₃]^- is replaced by the less coordinating [PF₆]^- anion.

A wrinkling-based method for investigating glassy polymer film relaxation as a function of film thickness and temperature

We investigate the relaxation dynamics of thin polymer films at temperatures below the bulk glass transition T_g by first compressing polystyrene films supported on a polydimethylsiloxane substrate to create wrinkling patterns and then observing the slow relaxation of the wrinkled films back to their final equilibrium flat state by small angle light scattering. As with recent relaxation measurements on thin glassy films reported by Fakhraai and co-workers, we find the relaxation time of our wrinkled films to be strongly dependent on film thickness below an onset thickness on the order of 100 nm. By varying the temperature between room temperature and T_g (≈100 °C), we find that the relaxation time follows an Arrhenius-type temperature dependence to a good approximation at all film thicknesses investigated, where both the activation energy and the relaxation time pre-factor depend appreciably on film thickness. The wrinkling relaxation curves tend to cross at a common temperature somewhat below T_g, indicating an entropy-enthalpy compensation relation between the activation free energy parameters. This compensation effect has also been observed recently in simulated supported polymer films in the high temperature Arrhenius relaxation regime rather than the glassy state. In addition, we find that the film stress relaxation function, as well as the height of the wrinkle ridges, follows a stretched exponential time dependence and the short-time effective Young's modulus derived from our modeling decreases sigmoidally with increasing temperature-both characteristic features of glassy materials. The relatively facile nature of the wrinkling-based measurements in comparison to other film relaxation measurements makes our method attractive for practical materials development, as well as fundamental studies of glass formation.

Energy Renormalization Method for the Coarse-Graining of Polymer Viscoelasticity

Developing temperature transferable coarse-grained (CG) models is essential for the computational prediction of polymeric glass-forming (GF) material behavior, but their dynamics are often greatly altered from those of all-atom (AA) models mainly because of the reduced fluid configurational entropy under coarse-graining. To address this issue, we have recently introduced an energy renormalization (ER) strategy that corrects the activation free energy of the CG polymer model by renormalizing the cohesive interaction strength ε as a function of temperature T, i.e., ε(T), thus semiempirically preserving the T-dependent dynamics of the AA model. Here we apply our ER method to consider-in addition to T-dependency-the frequency f-dependent polymer viscoelasticity. Through smallamplitude oscillatory shear molecular dynamics simulations, we show that changing the imposed oscillation f on the CG systems requires changes in ε values (i.e., ε(T, f)) to reproduce the AA viscoelasticity. By accounting for the dynamic fragility of polymers as a material parameter, we are able to predict ε(T, f) under coarse-graining in order to capture the AA viscoelasticity, and consequently the activation energy, across a wide range of T and f in the GF regime. Specifically, we showcase our achievements on two representative polymers of distinct fragilities, polybutadiene (PB) and polystyrene (PS), and show that our CG models are able to sample viscoelasticity up to the megahertz regime, which approaches state-of-the-art experimental resolutions, and capture results sampled via AA simulations and prior experiments.

Thermodynamic scaling of relaxation: insights from anharmonic elasticity

Using molecular dynamics simulations of a molecular liquid, we investigate the thermodynamic scaling (TS) of the structural relaxation time [Formula: see text] in terms of the quantity [Formula: see text], where T and ρ are the temperature and density, respectively. The liquid does not exhibit strong virial-energy correlations. We propose a method for evaluating both the characteristic exponent [Formula: see text] and the TS master curve that uses experimentally accessible quantities that characterise the anharmonic elasticity and does not use details about the microscopic interactions. In particular, we express the TS characteristic exponent [Formula: see text] in terms of the lattice Grüneisen parameter [Formula: see text] and the isochoric anharmonicity [Formula: see text]. An analytic expression of the TS master curve of [Formula: see text] with [Formula: see text] as the key adjustable parameter is found. The comparison with the experimental TS master curves and the isochoric fragilities of 34 glassformers is satisfying. In a few cases, where thermodynamic data are available, we test (i) the predicted characteristic exponent [Formula: see text] and (ii) the isochoric anharmonicity [Formula: see text], as drawn by the best fit of the TS of the structural relaxation, against the available thermodynamic data. A linear relation between the isochoric fragility and the isochoric anharmonicity [Formula: see text] is found and compared favourably with the results of experiments with no adjustable parameters. A relation between the increase of the isochoric vibrational heat capacity due to anharmonicity and the isochoric fragility is derived.

Direct links between dynamical, thermodynamic, and structural properties of liquids: Modeling results

We develop an approach to liquid thermodynamics based on collective modes. We perform extensive molecular-dynamics simulations of noble, molecular, and metallic liquids, and we provide direct evidence that liquid energy and specific heat are well-described by the temperature dependence of the Frenkel (hopping) frequency. The agreement between predicted and calculated thermodynamic properties is seen in the notably wide range of temperature spanning tens of thousands of Kelvin. The range includes both subcritical liquids and supercritical fluids. We discuss the structural crossover and interrelationships between the structure, dynamics, and thermodynamics of liquids and supercritical fluids.

Long-Lived Neighbors Determine the Rheological Response of Glasses

Glasses exhibit a liquidlike structure but a solidlike rheological response with plastic deformations only occurring beyond yielding. Thus, predicting the rheological behavior from the microscopic structure is difficult, but important for materials science. Here, we consider colloidal suspensions and propose to supplement the static structural information with the local dynamics, namely, the rearrangement and breaking of the cage of neighbors. This is quantified by the mean squared nonaffine displacement and the number of particles that remain nearest neighbors for a long time, i.e., long-lived neighbors, respectively. Both quantities are followed under shear using confocal microscopy and are the basis to calculate the affine and nonaffine contributions to the elastic stress, which is complemented by the viscoelastic stress to give the total stress. During start-up of shear, the model predicts three transient regimes that result from the interplay of affine, nonaffine, and viscoelastic contributions. Our prediction quantitatively agrees with rheological data and their dependencies on volume fraction and shear rate.

Energy-Renormalization for Achieving Temperature Transferable Coarse-Graining of Polymer Dynamics

The bottom-up prediction of the properties of polymeric materials based on molecular dynamics simulation is a major challenge in soft matter physics. Coarse-grained (CG) models are often employed to access greater spatiotemporal scales required for many applications, but these models normally experience significantly altered thermodynamics and highly accelerated dynamics due to the reduced number of degrees of freedom upon coarse-graining. While CG models can be calibrated to meet certain properties at particular state points, there is unfortunately no temperature transferable and chemically specific coarse-graining method that allows for modeling of polymer dynamics over a wide temperature range. Here, we pragmatically address this problem by "correcting" for deviations in activation free energies that occur upon coarse-graining the dynamics of a model polymeric material (polystyrene). In particular, we propose a new strategy based on concepts drawn from the Adam-Gibbs (AG) theory of glass formation. Namely we renormalize the cohesive interaction strength and effective interaction length-scale parameters to modify the activation free energy. We show that this energy-renormalization method for CG modeling allows accurate prediction of atomistic dynamics over the Arrhenius regime, the non-Arrhenius regime of glass formation, and even the non-equilibrium glassy regime, thus allowing for the predictive modeling of dynamic properties of polymer over the entire range of glass formation. Our work provides a practical scheme for establishing temperature transferable coarse-grained models for predicting and designing the properties of polymeric materials.

Nonaffine displacements and the nonlinear response of a strained amorphous solid

We demonstrate that irreversible structural reorganization is not necessary for the observation of yield behavior in an amorphous solid. While the majority of solids strained to their yield point do indeed undergo an irreversible reorganization, we find that a significant fraction of solids exhibits yield via a reversible strain. We also demonstrate that large instantaneous strains in excess of the yield stress can result in complete stress relaxation, a result of the large nonaffine motions driven by the applied strain. The empirical similarity of the dependence of the ratio of stress over strain on the nonaffine mean-square displacement to that for the shear modulus obtained from quiescent liquid at nonzero temperature supports the proposition that rigidity depends on the size of the sampled configurational space only and is insensitive to how this space is sampled.

Collective modes and thermodynamics of the liquid state

Strongly interacting, dynamically disordered and with no small parameter, liquids took a theoretical status between gases and solids with the historical tradition of hydrodynamic description as the starting point. We review different approaches to liquids as well as recent experimental and theoretical work, and propose that liquids do not need classifying in terms of their proximity to gases and solids or any categorizing for that matter. Instead, they are a unique system in their own class with a notably mixed dynamical state in contrast to pure dynamical states of solids and gases. We start with explaining how the first-principles approach to liquids is an intractable, exponentially complex problem of coupled non-linear oscillators with bifurcations. This is followed by a reduction of the problem based on liquid relaxation time τ representing non-perturbative treatment of strong interactions. On the basis of τ, solid-like high-frequency modes are predicted and we review related recent experiments. We demonstrate how the propagation of these modes can be derived by generalizing either hydrodynamic or elasticity equations. We comment on the historical trend to approach liquids using hydrodynamics and compare it to an alternative solid-like approach. We subsequently discuss how collective modes evolve with temperature and how this evolution affects liquid energy and heat capacity as well as other properties such as fast sound. Here, our emphasis is on understanding experimental data in real, rather than model, liquids. Highlighting the dominant role of solid-like high-frequency modes for liquid energy and heat capacity, we review a wide range of liquids: subcritical low-viscous liquids, supercritical state with two different dynamical and thermodynamic regimes separated by the Frenkel line, highly-viscous liquids in the glass transformation range and liquid-glass transition. We subsequently discuss the fairly recent area of liquid-liquid phase transitions, the area where the solid-like properties of liquids have become further apparent. We then discuss gas-like and solid-like approaches to quantum liquids and theoretical issues that are similar to the classical case. Finally, we summarize the emergent view of liquids as a unique system with a mixed dynamical state, and list several areas where interesting insights may appear and continue the extraordinary liquid story.

Slow Water Dynamics near a Glass Transition or a Solid Interface: A Common Rationale

Performing molecular dynamics simulations, we investigate the enormous slowdowns of water dynamics when approaching a glass transition or a solid interface. We show that both effects can be described on common grounds within a theoretical framework, which was recently proposed by Schweizer et al. and considers coupled local hopping and elastic distortion. For confined water, we correctly describe the variation of the α-relaxation time, τα, as a function of both temperature and position with respect to the interface. Exploiting our knowledge of a cooperative length scale ξ(T) from the confinement studies, we quantitatively rationalize the glassy slowdown, τα(T), and the Stokes-Einstein breakdown of bulk water. For both confined and bulk liquid, variations of the α-relaxation time are intimately related to changes of the cage-rattling amplitude.

The kinetic fragility of liquids as manifestation of the elastic softening

We show that the fragility m , the steepness of the viscosity and relaxation time close to the vitrification, increases with the degree of elastic softening, i.e. the decrease of the elastic modulus with increasing temperature, in a universal way. This provides a novel connection between the thermodynamics, via the modulus, and the kinetics. The finding is evidenced by numerical simulations and comparison with the experimental data of glassformers with widely different fragilities (33 ≤ m ≤ 115), leading to a fragility-independent elastic master curve extending over eighteen decades in viscosity and relaxation time. The master curve is accounted for by a cavity model pointing out the roles of both the available free volume and the cage softness. A major implication of our findings is that ultraslow relaxations, hardly characterised experimentally, become predictable by linear elasticity. As an example, the viscosity of supercooled silica is derived over about fifteen decades with no adjustable parameters.

Weak links between fast mobility and local structure in molecular and atomic liquids

We investigate by molecular-dynamics simulations, the fast mobility-the rattling amplitude of the particles temporarily trapped by the cage of the neighbors-in mildly supercooled states of dense molecular (linear trimers) and atomic (binary mixtures) liquids. The mixture particles interact by the Lennard-Jones potential. The non-bonded particles of the molecular system are coupled by the more general Mie potential with variable repulsive and attractive exponents in a range which is a characteristic of small n-alkanes and n-alcohols. Possible links between the fast mobility and the geometry of the cage (size and shape) are searched. The correlations on a per-particle basis are rather weak. Instead, if one groups either the particles in fast-mobility subsets or the cages in geometric subsets, the increase of the fast mobility with both the size and the asphericity of the cage is revealed. The observed correlations are weak and differ in states with equal relaxation time. Local forces between a tagged particle and the first-neighbour shell do not correlate with the fast mobility in the molecular liquid. It is concluded that the cage geometry alone is unable to provide a microscopic interpretation of the known, universal link between the fast mobility and the slow structural relaxation. We suggest that the particle fast dynamics is affected by regions beyond the first neighbours, thus supporting the presence of collective, extended fast modes.

Fluctuation-dissipation theorem and the dielectric response in supercooled liquids

We consider the correlation between static conductivity and dynamic dielectric relaxation in a number of polar organic liquids. Experimental evidence suggests that in the simple cases the linear dependence between characteristic frequency of relaxation process and the value of static susceptibility is observed. However, this proportionality can be broken due to the appearance of additional relaxation processes (secondary or high-frequency ones) so it can be confused with the "fractional" variant of Debye-Stokes-Einstein relation.

Structural relaxation is a scale-free process

We show that in deeply supercooled liquids, structural relaxation proceeds via the accumulation of Eshelby events, i.e. local rearrangements that create long-ranged and anisotropic stresses in the surrounding medium. Such events must be characterized using tensorial observables and we construct an analytical framework to probe their correlations using local stress data. By analyzing numerical simulations, we then demonstrate that events are power-law correlated in space, with a time-dependent amplitude which peaks at the alpha relaxation time τ(α). This effect, which becomes stronger near the glass transition, results from the increasingly important role of local stress fluctuations in facilitating relaxation events. Our finding precludes the existence of any length scale beyond which the relaxation process decorrelates.

Elastic signature of flow events in supercooled liquids under shear

Using numerical simulation of a 2D Lennard-Jones system, we study the crossover from shear thinning to Newtonian flow. We find that the short-time elastic response of our system essentially does not change through this crossover, and show that, in the Newtonian regime, thermal activation triggers shear transformations, i.e., local irreversible shear events that produce Eshelby (long-ranged, anisotropic) deformation fields as previously seen in low-T glasses. Quite surprisingly, these Eshelby fields are found to persist much beyond the α-relaxation time, and shear thinning to coincide with the emergence of correlations between shear relaxation centers.

Disorder-assisted melting and the glass transition in amorphous solids

The mechanical response of solids depends on temperature, because the way atoms and molecules respond collectively to deformation is affected at various levels by thermal motion. This is a fundamental problem of solid state science and plays a crucial role in materials science. In glasses, the vanishing of shear rigidity upon increasing temperature is the reverse process of the glass transition. It remains poorly understood due to the disorder leading to nontrivial (nonaffine) components in the atomic displacements. Our theory explains the basic mechanism of the melting transition of amorphous (disordered) solids in terms of the lattice energy lost to this nonaffine motion, compared to which thermal vibrations turn out to play only a negligible role. The theory is in good agreement with classic data on melting of amorphous polymers (for which no alternative theory can be found in the literature) and offers new opportunities in materials science.

Entropic vs. elastic models of fragility of glass-forming liquids: two sides of the same coin?

The two most influential atomistic models that have been proposed in the literature to explain the temperature dependent activation energy of viscous flow of a glass-forming liquid, i.e., its fragility, are the configurational entropy model of Adam and Gibbs [J. Chem. Phys. 43, 139 (1965)] and the elastic "shoving" model of Dyre et al. [J. Non-Cryst. Solids 352, 4635 (2006)]. Here we demonstrate a qualitative equivalence between these two models starting from the well-established general relationships between the interatomic potentials, elastic constants, structural rearrangement, and entropy in amorphous materials. The unification of these two models provides important predictions that are consistent with experimental observations and shed new light into the problem of glass transition.

Isoviscosity lines and the liquid-glass transition in simple liquids

This article presents the study of the generic behavior of viscosity of liquids based on some simple theoretical models, the soft-spheres and Lennard-Jones systems. The use of these simple models allows us to investigate in detail the viscosity behavior in a wide range of temperatures and pressures including the high-temperature-high-pressure limits. Based on the simulation results, we discuss the shape of isoviscosity lines and analyze the glass transition at high temperatures and high pressures. Despite the fact that the viscosity drastically increases in the limit of high temperatures and high pressures along the melting line, the relaxation time rapidly decreases in this region, and the system becomes further from the glass transition.