Perspective: Supercooled liquids and glasses

Decoupling of component diffusion in a glass-forming Zr(46.75)Ti(8.25)Cu(7.5)Ni(10)Be(27.5) melt far above the liquidus temperature

We report (95)Zr and (57)Co radiotracer diffusivities and viscosity data in the equilibrium liquid state of a bulk metallic glass forming Zr(46.75)Ti(8.25)Cu(7.5)Ni(10)Be(27.5) melt (Vitreloy 4) far above the liquidus temperature T(l) that are not affected by convection, as evidenced via quasielastic neutron scattering. Zr diffusion is strongly decoupled from diffusion of the smaller components by more than a factor of 4 at T(l), although it obeys the Stokes-Einstein equation. The results suggest that, in the present Zr-based metallic glass forming systems, diffusion and viscous flow start to develop solidlike, i.e., energy-landscape-controlled, features already in the stable liquid state more than 300 K above the mode coupling temperature T(c).

Non-Gaussian rotational diffusion in heterogeneous media

We employ a simple model for rotational diffusivity DR of dumbbells in porous media in order to study spatially heterogeneous and non-Gaussian dynamics at Fickian time scales. We obtain the distribution P(DR) of DR's of single dumbbells for both ergodic and nonergodic systems. When a pore percolating network disappears beyond the pore percolation transition and the rotational dynamics becomes nonergodic, each single dumbbell undergoes Gaussian rotational dynamics but with different DR, which depends solely on the local pore structure. We also construct a map of heterogeneous dynamic regions and illustrate that such seemingly Fickian but non-Gaussian dynamics could be understood as the linear combination of the Gaussian rotational displacement distribution functions of each dumbbell. With a percolating pore network, the rotational dynamics becomes ergodic, and P(DR) is a δ function at the average value of DR.

Dependence of the atomic level Green-Kubo stress correlation function on wavevector and frequency: molecular dynamics results from a model liquid

We report on a further investigation of a new method that can be used to address vibrational dynamics and propagation of stress waves in liquids. The method is based on the decomposition of the macroscopic Green-Kubo stress correlation function into the atomic level stress correlation functions. This decomposition, as was demonstrated previously for a model liquid studied in molecular dynamics simulations, reveals the presence of stress waves propagating over large distances and a structure that resembles the pair density function. In this paper, by performing the Fourier transforms of the atomic level stress correlation functions, we elucidate how the lifetimes of the stress waves and the ranges of their propagation depend on their frequency, wavevector, and temperature. These results relate frequency and wavevector dependence of the generalized viscosity to the character of propagation of the shear stress waves. In particular, the results suggest that an increase in the value of the frequency dependent viscosity at low frequencies with decrease of temperature is related to the increase in the ranges of propagation of the stress waves of the corresponding low frequencies. We found that the ranges of propagation of the shear stress waves of frequencies less than half of the Einstein frequency extend well beyond the nearest neighbor shell even above the melting temperature. The results also show that the crossover from quasilocalized to propagating behavior occurs at frequencies usually associated with the Boson peak.

Single molecule rotational probing of supercooled liquids

Much of the interesting behavior that has been observed in supercooled liquids appears to be related to dynamic heterogeneity, the presence of distinct dynamic environments - with no apparent underlying structural basis - in these systems. To most directly interrogate these environments, proposed to span regions just a few nanometers across, molecular length scale probes are required. Single molecule fluorescent microscopy was introduced to the field a decade ago and has provided strong evidence of dynamic heterogeneity in supercooled systems. However, only more recently has the full set of challenges associated with interpreting results of these experiments been described. With a fuller understanding of these challenges in hand, single molecule measurements can be employed to provide a more precise picture of dynamic heterogeneity in supercooled liquids and other complex systems. In this tutorial review, experimental and data analysis details are presented for the most commonly employed single molecule approach to studying supercooled liquids, the measurement of rotational dynamics of single molecule probes. Guidance is provided in experimental set-up and probe selection, with a focus on choices that affect data interpretation and probe sensitivity to dynamic heterogeneity.

The β-relaxation in metallic glasses

Focusing on metallic glasses as model systems, we review the features and mechanisms of the β-relaxations, which are intrinsic and universal to supercooled liquids and glasses, and demonstrate their importance in understanding many crucial unresolved issues in glassy physics and materials science, including glass transition phenomena, mechanical properties, shear-banding dynamics and deformation mechanisms, diffusion and the breakdown of the Stokes–Einstein relation as well as crystallization and stability of glasses. We illustrate that it is an attractive prospect to incorporate these insights into the design of new glassy materials with extraordinary properties. We also outline important questions regarding the nature of β-relaxations and highlight some emerging research directions in this still-evolving field.

Silicon-Chip-Based Dielectric Spectroscopy for Conductivity and Molecular Dynamics Studies of Organic Films

Interdigital electrodes fabricated by standard lithography on silicon chips are employed to probe the dipolar molecular dynamics and electric conduction properties of thin rhodamine films grown with two different methods. The conductivity is due to electronic charge carriers, and at around room-temperature, it is higher by 1 order of magnitude in solution-deposited films than in thermally evaporated ones. The organic material exhibits two intrinsic dynamic processes, of which the one at higher temperature is due to the orientational motion of the dipole moment of the rhodamine units, while the one at lower temperature is due to the motion of a local dipole associated with the chlorine counterions and is absent in thermally evaporated films. Our results show that thin-film dielectric spectroscopy is an easily implementable and versatile tool to extract valuable information on thin organic films.

Renewal events in glass-forming liquids

On cooling toward the glass transition temperature, glass-forming liquids display long periods of localized motion interrupted by fast "jumps" in the single-particle trajectories. Several theoretical models based on these single-particle jumps have been proposed, most prominently the continuous-time random walk (CTRW). The central assumption of the CTRW is that jumps are renewal events, i.e. that the internal clock of a particle can be reset upon a jump. In this paper, I present an easy-to-implement method to test whether jumps detected in a supercooled liquid or glass are renewal events or not. The test was applied to molecular dynamics simulations of a short-chain polymer melt, demonstrating that the jumps can in fact be treated as renewal events. The test further revealed that additional relaxation processes are present which are not accounted for in the CTRW picture, highlighting the limitations of this approach. The notion of renewal events in glass-forming systems could be a very important building block for the interpretation of aging and the glass transition. Furthermore, it could have practical implications for the study of non-equilibrium dynamics in glasses as well as mechanical rejuvenation.

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.

Physical aging, the local dynamics of glass-forming polymers under nanoscale confinement

The glass transition temperature marks a point below which a material's properties change significantly, and it is well-established that confinement to the nanoscale modifies the properties of glass-forming materials. We use molecular dynamics simulations to investigate the dynamics and aging behavior of model glass-forming polymers near and below the glass transition temperature of bulk and confined films. We show that both relaxation times and physical age rates vary similarly throughout a free-standing polymer film at temperatures close to the bulk glass transition temperature, where the surfaces have both lower relaxation times and physical age rates. Moreover, we provide evidence suggesting that string lengths in the bulk control dynamic length scales in the film. This realization, combined with the similarity between aging behavior and dynamic profiles, has implications for design rationale in the microelectronics industry.

Suppression of tunneling two-level systems in ultrastable glasses of indomethacin

Glasses and other noncrystalline solids exhibit thermal and acoustic properties at low temperatures anomalously different from those found in crystalline solids, and with a remarkable degree of universality. Below a few kelvin, these universal properties have been successfully interpreted using the tunneling model, which has enjoyed (almost) unanimous recognition for decades. Here we present low-temperature specific-heat measurements of ultrastable glasses of indomethacin that clearly show the disappearance of the ubiquitous linear contribution traditionally ascribed to the existence of tunneling two-level systems (TLS). When the ultrastable thin-film sample is thermally converted into a conventional glass, the material recovers a typical amount of TLS. This remarkable suppression of the TLS found in ultrastable glasses of indomethacin is argued to be due to their particular anisotropic and layered character, which strongly influences the dynamical network and may hinder isotropic interactions among low-energy defects, rather than to the thermodynamic stabilization itself. This explanation may lend support to the criticisms by Leggett and others [Yu CC, Leggett AJ (1988) Comments Condens Matter Phys 14(4):231-251; Leggett AJ, Vural DC (2013) J Phys Chem B 117(42):12966-12971] to the standard tunneling model, although more experiments in different kinds of ultrastable glasses are needed to ascertain this hypothesis.

Continuous-time random-walk approach to supercooled liquids. I. Different definitions of particle jumps and their consequences

Single-particle trajectories in supercooled liquids display long periods of localization interrupted by "fast moves." This observation suggests a modeling by a continuous-time random walk (CTRW). We perform molecular dynamics simulations of equilibrated short-chain polymer melts near the critical temperature of mode-coupling theory Tc and extract "moves" from the monomer trajectories. We show that not all moves comply with the conditions of a CTRW. Strong forward-backward correlations are found in the supercooled state. A refinement procedure is suggested to exclude these moves from the analysis. We discuss the repercussions of the refinement on the jump-length and waiting-time distributions as well as on characteristic time scales, such as the average waiting time ("exchange time") and the average time for the first move ("persistence time"). The refinement modifies the temperature (T) dependence of these time scales. For instance, the average waiting time changes from an Arrhenius-type to a Vogel-Fulcher-type T dependence. We discuss this observation in the context of the bifurcation of the α process and (Johari) β process found in many glass-forming materials to occur near Tc. Our analysis lays the foundation for a study of the jump-length and waiting-time distributions, their temperature and chain-length dependencies, and the modeling of the monomer dynamics by a CTRW approach in the companion paper [J. Helfferich et al., Phys. Rev. E 89, 042604 (2014)].

Probing the glass transition from structural and vibrational properties of zero-temperature glasses

We find that the density dependence of the glass transition temperature of Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) systems can be predicted from properties of the zero-temperature (T=0) glasses. Below a crossover density ρs, LJ and WCA glasses show different structures, leading to different vibrational properties and consequently making LJ glasses more stable with higher glass transition temperatures than WCA ones. Above ρs, structural and vibrational quantities of the T=0 glasses show scaling collapse. From scaling relations and dimensional analysis, we predict a density scaling of the glass transition temperature, in excellent agreement with simulation results. We also propose an empirical expression of the glass transition temperature using structural and vibrational properties of the T=0 glasses, which works well over a wide range of densities.

Structural relaxation of vapor-deposited molecular glasses and supercooled liquids

The properties of vapor-deposited molecular glasses largely depend on deposition conditions, and stable and/or dense glasses are formed with several compounds.

Note: Effects of adding a viscosity-increasing 2 nm-size molecule on dielectric relaxation features and the dynamic heterogeneity view

To investigate the extent to which dynamic heterogeneity determines the relaxation features, we studied the effect of adding a viscosity-increasing, ∼2 nm size polyhederal oligomeric silsesquioxane (POSS) molecule on dielectric relaxation. Addition of POSS increased the relaxation time of the α-process, τ(m,α), and decreased the non-exponential response parameter, β, which is compatible with the dynamic heterogeneity view. But β varied monotonically with the composition, and there was little variation of β with the temperature, which seems incompatible with that view. Effects of adding nm-size particles in simulations may reveal how the size of dynamically heterogeneous regions affects τ(m,α) and β.

Ultrastable metallic glass

A new metallic glass, which was created by vapour deposition at an appropriately high substrate temperature, shows exceptional thermal stability, and enhanced glass transition temperature and elastic modulus. Comparing this new material with other organic glasses prepared by similar routes and known as ultrastable glasses demonstrates the formation of ultrastable glassy materials correlates to the important concept of fragility.

Manifestation of random first-order transition theory in Wigner glasses

We use Brownian dynamics simulations of a binary mixture of highly charged spherical colloidal particles to test some of the predictions of the random first-order transition (RFOT) theory [Phys. Rev. Lett. 58, 2091 (1987); Phys. Rev. A 40, 1045 (1989)]. In accord with mode-coupling theory and RFOT, we find that as the volume fraction of the colloidal particles ϕ approaches the dynamical transition value ϕ(A), three measures of dynamics show an effective ergodic to nonergodic transition. First, there is a dramatic slowing down of diffusion, with the translational diffusion constant decaying as a power law as ϕ→ϕ(A)(-). Second, the energy metric, a measure of ergodicity breaking in classical many-body systems, shows that the system becomes effectively nonergodic as ϕ(A) is approached. Finally, the time t(*), at which the four-point dynamical susceptibility achieves a maximum, also increases as a power law near ϕ(A). Remarkably, the translational diffusion coefficients, ergodic diffusion coefficient, and (t(*))(-) all vanish as (ϕ(-1)-ϕ(A)(-1))(γ) with both ϕ(A)(≈0.1) and γ being the roughly the same for all three quantities. Above ϕ(A), transport involves crossing free energy barriers. In this regime, the density-density correlation function decays as a stretched exponential [exp-(t/τ(α))(β)] with β≈0.45. The ϕ dependence of the relaxation time τ(α) could be fit using the Vogel-Tamman-Fulcher law with the ideal glass transition at ϕ(K)≈0.47. By using a local entropy measure, we show that the law of large numbers is not obeyed above ϕ(A), and gives rise to subsample to subsample fluctuations in all physical observables. We propose that dynamical heterogeneity is a consequence of violation of law of large numbers.

Crossover to surface flow in supercooled unentangled polymer films

We study the driven flow of an unentangled glassy polymer film with a free upper surface and supported below by a substrate using nonequilibrium molecular dynamics simulations based on a bead-spring model. Above the glass transition temperature T(g), simple Poiseuille laminar flow is observed with the film mobility defined as the flow current density per unit pressure gradient scaling as h(3) with the film thickness h. Below T(g), the film mobility becomes independent of h, signifying surface transport. This is in full agreement with recent experiments on the time evolution of capillary waves in polystyrene films supported by silica. A mobile layer is found responsible for the surface transport, as previously conjectured. Our result also shows that it has a velocity profile decaying exponentially into the bulk.

Density dependence of dynamical heterogeneity in fluid methanol

Brillouin and Raman scattering experiments on methanol through its glass transition under pressure are reported. The Brillouin scattering data were analyzed using viscoelastic theory and a fit to the Vinet equation of state. The variation in the linewidth of the longitudinal acoustic mode with pressure shows a broad maximum centered around 3 GPa. The pressure evolution of the relaxation time in the GHz range is obtained, and the Raman data are analyzed in terms of the Boson peak and its associated relaxation time in the THz range. The pressure evolution of these two relaxation processes extends previous determinations of relaxations at lower frequency based on dielectric measurements in supercooled methanol. The relaxation processes in glass-forming methanol have now been investigated over a wide frequency range and their evolution followed over a large variation of density.

Heterogeneous relaxation dynamics in amorphous materials under cyclic loading

Molecular dynamics simulations are performed to investigate heterogeneous dynamics in amorphous glassy materials under oscillatory shear strain. We consider three-dimensional binary Lennard-Jones mixture well below the glass transition temperature. The structural relaxation and dynamical heterogeneity are quantified by means of the self-overlap order parameter and the dynamic susceptibility. We found that at sufficiently small strain amplitudes, the mean square displacement exhibits a broad subdiffusive plateau and the system undergoes nearly reversible deformation over about 10(4) cycles. Upon increasing strain amplitude, the transition to the diffusive regime occurs at shorter time intervals and the relaxation process involves intermittent bursts of large particle displacements. The detailed analysis of particle hopping dynamics and the dynamic susceptibility indicates that mobile particles aggregate into clusters whose sizes increase at larger strain amplitudes. Finally, the correlation between particle mobilities in consecutive time intervals demonstrates that dynamic facilitation becomes increasingly pronounced at larger strain amplitudes.

An explanation of the differences in diffusivity of the components of the metallic glass Pd43Cu27Ni10P20

Bartsch et al. [Phys. Rev. Lett. 104, 195901 (2010)] reported measurements of the diffusivities of different components of the multi-component bulk metallic glass Pd43Cu27Ni10P20. The diffusion of the largest Pd and the smallest P was found to be drastically different. The Stokes-Einstein relation breaks down when considering the P constituent atom, while the relation is obeyed by the Pd atom over 14 orders of magnitude of change in Pd diffusivity. This difference in behavior of Pd and P poses a problem challenging for explanation. With the assist of a recent finding in metallic glasses that the β-relaxation and the diffusion of the smallest component are closely related processes by Yu et al. [Phys. Rev. Lett. 109, 095508 (2012)], we use the Coupling Model to explain the observed difference between P and Pd quantitatively. The same model also explains the correlation between property of the β-relaxation with fragility found in the family of (CexLa1-x)68Al10Cu20Co2 with 0 ≤ x ≤ 1.

Anomalous diffusion of driven particles in supercooled liquids

We perform nonequilibrium dynamics simulations of a binary Lennard-Jones mixture in which an external force is applied on a single tagged particle. For the diffusive properties of this particle parallel to the force, superdiffusive behavior at intermediate times as well as giant long-time diffusivity is observed. A quantitative description of this nontrivial behavior is given by a continuous time random walk analysis of the system in configuration space. We further demonstrate that the same physical properties which are responsible for the superdiffusivity in nonequilibrium systems also determine the non-Gaussian parameter in equilibrium systems.