Perspective: Supercooled liquids and glasses

Low-Temperature Heat Capacity Anomalies in Ordered and Disordered Phases of Normal and Deuterated Thiophene

We measured the specific heat C_p of normal (C₄H₄S) and deuterated (C₄D₄S) thiophene in the temperature interval of 1 ≤ T, K ≤ 25. C₄H₄S exhibits a metastable phase II₂ and a stable phase V, both with frozen orientational disorder (OD), whereas C₄D₄S exhibits a metastable phase II₂, which is analogous to the OD phase II₂ of C₄H₄S and a fully ordered stable phase V. Our measurements demonstrate the existence of a large bump in the heat capacity of both stable and metastable C₄D₄S and C₄H₄S phases at temperatures of ∼10 K, which significantly departs from the expected Debye temperature behavior of C_p ≈ T³. This case study demonstrates that the identified low-temperature C_p anomaly, typically referred to as a "Boson-peak" in the context of glassy crystals, is not exclusive of disordered materials.

Emergence of cooperatively reorganizing cluster and super-Arrhenius dynamics of fragile supercooled liquids

In this paper, we develop a theory to calculate the structural relaxation time τ_{α} of fragile supercooled liquids. Using the information of the configurational entropy and structure, we calculate the number of dynamically free, metastable, and stable neighbors around a central particle. In supercooled liquids, the cooperatively reorganizing clusters (CRCs) in which the stable neighbors form "stable" nonchemical bonds with the central particle emerge. For an event of relaxation to take place, these bonds have to reorganize irreversibly; the energy involved in the processes is the effective activation energy of relaxation. The theory brings forth a temperature T_{a} and a temperature-dependent parameter ψ(T) which characterize slowing down of dynamics on cooling. It is shown that the value of ψ(T) is equal to 1 for T>T_{a}, indicating that the underlying microscopic mechanism of relaxation is dominated by the entropy-driven processes, while for T<T_{a}, ψ(T) decreases on cooling, indicating the emergence of the energy-driven processes. This crossover of ψ(T) from high to low temperatures explains the crossover seen in τ_{α}. The dynamics of systems that may have similar static structure but very different dynamics can be understood in terms of ψ(T). We present results for the Kob-Anderson model for three densities and show that the calculated values of τ_{α} are in excellent agreement with simulation values for all densities. We also show that when ψ(T), τ_{α}, and other quantities are plotted as a function of T/T_{a} (or T_{a}/T), the data collapse on master curves.

On the Mutual Relationships between Molecular Probe Mobility and Free Volume and Polymer Dynamics in Organic Glass Formers: cis-1,4-poly(isoprene)

We report on the reorientation dynamics of small spin probe 2,2,6,6-tetramethylpiperidinyl-1-oxyl (TEMPO) in cis-1,4-poly(isoprene) (cis-1,4-PIP10k) from electron spin resonance (ESR) and the free volume of cis-1,4-PIP10k from positron annihilation lifetime spectroscopy (PALS) in relation to the high-frequency relaxations of cis-1,4-PIP10k using light scattering (LS) as well as to the slow and fast processes from broadband dielectric spectroscopy (BDS) and neutron scattering (NS). The hyperfine coupling constant, 2A_zz '(T), and the correlation times, τ _c(T), of cis-1,4-PIP10k/TEMPO system as a function of temperature exhibit several regions of the distinct spin probe TEMPO dynamics over a wide temperature range from 100 K up to 350 K. The characteristic ESR temperatures of changes in the spin probe dynamics in cis-1,4-PIP10k/TEMPO system are closely related to the characteristic PALS ones reflecting changes in the free volume expansion from PALS measurement. Finally, the time scales of the slow and fast dynamics of TEMPO in cis-1,4-PIP10k are compared with all of the six known slow and fast relaxation modes from BDS, LS and NS techniques with the aim to discuss the controlling factors of the spin probe reorientation mobility in polymer, oligomer and small molecular organic glass-formers.

Sign flipping of spontaneous polarization in vapour-deposited films of small polar organic molecules

Films of polar molecules vapour-deposited on sufficiently cold substrates are not only amorphous, but also exhibit charge polarization across their thickness. This is an effect known for 50 years, but it is very poorly understood and no mechanism exists in the literature that can explain and predict it. We investigated this bulk effect for 18 small organic molecules as a function of substrate temperature (30-130 K). We found that, as a rule, alcohol films have the negative end on the vacuum side at all temperatures. Alkyl acetates and toluene showed positive voltages which reached a maximum around the middle of the temperature range investigated. Tetrahydrofuran showed positive voltages which dropped with increasing deposition temperature. Diethyl ether, acetone, propanal, and butanal showed positive film voltages at low temperatures, negative at intermediate temperatures and again positive voltages at higher temperatures. In all cases, film voltages were monitored during heating leading to film evaporation. Film voltages were irreversibly eliminated before film elimination, but voltage profiles during temperature ramps differed vastly depending on compound and deposition temperature. In general, there was a gradual voltage reduction, but propanal, butanal, and diethyl ether showed a change in voltage sign during temperature ramp in films deposited at low temperatures. All these data expand substantially the experimental information regarding spontaneous polarization in vapour-deposited films, but still require complementary measurements as well as numerical simulations for a detailed explanation of the phenomenon.

Visualizing free-energy landscapes for four hard disks

We present a simple model system with four hard disks moving in a circular region for which free-energy landscapes can be directly calculated and visualized in two and three dimensions. We construct several energy landscapes for our system, and we explore the strengths and limitations of each in terms of understanding system dynamics, in particular the relationship between state transitions and free-energy barriers. We also demonstrate the importance of distinguishing between system dynamics in real space and those in landscape coordinates, and we show that care must be taken to appropriately combine dynamics with barrier properties to understand the transition rates. This simple model provides an intuitive way to understand free-energy landscapes, and it illustrates the benefits that free-energy landscapes can have over potential energy landscapes.

Why Is the Range of Timescale So Wide in Glass-Forming Liquid?

The viscosity and the relaxation time of a glass-forming liquid vary over 15 orders of magnitude before the liquid freezes into a glass. The rate of the change with temperature is characterized by liquid fragility. The mechanism of such a spectacular behavior and the origin of fragility have long been discussed, but it remains unresolved because of the difficulty of carrying out experiments and constructing theories that bridge over a wide timescale from atomic (ps) to bulk (minutes). Through the x-ray diffraction measurement and molecular dynamics simulation for metallic liquids we suggest that large changes in viscosity can be caused by relatively small changes in the structural coherence which characterizes the medium-range order. Here the structural coherence does not imply that of atomic-scale structure, but it relates to the coarse-grained density fluctuations represented by the peaks in the pair-distribution function (PDF) beyond the nearest neighbors. The coherence length is related to fragility and increases with decreasing temperature, and it diverges only at a negative temperature. This analysis is compared with several current theories which predict a phase transition near the glass transition temperature.

Morphological aspect of crystal nucleation in wall-confined supercooled metallic film

In this paper, we simulate the nucleation and growth of crystalline nuclei in a molybdenum film cooled at different rates confined between two amorphous walls. We also compare the results for the wall-confined and wall-free systems. We apply the same methodology as in the work (Kirova and Pisarev 2019J. Cryst. Growth528125266) which is based on reconstructing the probability density function for the largest crystalline nucleus in the system. The size of the nucleus and the asphericity parameter are considered as the reaction coordinates. We demonstrate that in both the free and confined systems there are two mechanisms of crystal growth: the attachment of atoms to the biggest crystal from the amorphous phase and the merging of the biggest crystal cluster with small ones (coalescence). We show that the attachment mechanism is dominant in the melt cooled down at a slower rate, and the mechanism gradually shifts to coalescence as cooling rate increases. We also observe the formation of long-lived crystal clusters and demonstrate that amorphous walls do not affect their geometric characteristics. However, system confined between walls demonstrates higher glass-forming ability.

Dynamic phase transitions in freestanding polymer thin films

After more than two decades of study, many fundamental questions remain unanswered about the dynamics of glass-forming materials confined to thin films. Experiments and simulations indicate that free interfaces enhance dynamics over length scales larger than molecular sizes, and this effect strengthens at lower temperatures. The nature of the influence of interfaces, however, remains a point of significant debate. In this work, we explore the properties of the nonequilibrium phase transition in dynamics that occurs in trajectory space between high- and low-mobility basins in a set of model polymer freestanding films. In thick films, the film-averaged mobility transition is broader than the bulk mobility transition, while in thin films it is a variant of the bulk result shifted toward a higher bias. Plotting this transition's local coexistence points against the distance from the films' surface shows thick films have surface and film-center transitions, while thin films practically have a single transition throughout the film. These observations are reminiscent of thermodynamic capillary condensation of a vapor-liquid phase between parallel plates, suggesting they constitute a demonstration of such an effect in a trajectory phase transition in the dynamics of a structural glass former. Moreover, this transition bears similarities to several experiments exhibiting anomalous behavior in the glass transition upon reducing film thickness below a material-dependent onset, including the broadening of the glass transition and the homogenization of surface and bulk glass transition temperatures.

Local Structure of Glassy Lithium Phosphorus Oxynitride Thin Films: A Combined Experimental and Ab Initio Approach

Lithium phosphorus oxynitride (LiPON) is an amorphous solid-state lithium ion conductor displaying exemplary cyclability against lithium metal anodes. There is no definitive explanation for this stability due to the limited understanding of the structure of LiPON. Herein, we provide a structural model of RF-sputtered LiPON. Information about the short-range structure results from 1D and 2D solid-state NMR experiments. These results are compared with first principles chemical shielding calculations of Li-P-O/N crystals and ab initio molecular dynamics-generated amorphous LiPON models to unequivocally identify the glassy structure as primarily isolated phosphate monomers with N incorporated in both apical and as bridging sites in phosphate dimers. Structural results suggest LiPON's stability is a result of its glassy character. Free-standing LiPON films are produced that exhibit a high degree of flexibility, highlighting the unique mechanical properties of glassy materials.

Hierarchical glass transition of hard hemidisks with local assemblies

Using computer simulation, we investigate the glass transition of a two-dimensional hard-hemidisk system. Upon increasing the packing fraction of the system, we find that the system vitrifies into a glass with local assembled discal "dimers", which are free to rotate in a collective way. The rotational mean square displacement does not exhibit the typical plateau (slowdown) like what occurs in the translational mean square displacement. This effect induces a pronounced violation of the rotational Stokes-Einstein relationship compared with the translational degree of freedom at the supercooled region. However, the obtained glass transition points in these two freedom degrees are found to be the same within the numerical accuracy, which is due to the strong positive spatial and dynamic correlation between translational and rotational slow-moving particles. Moreover, we find that the locally assembled dimers can serve as fast rotating gears facilitating the orientational relaxation in the system, and this suggests that the locally favored finite structures play an important role in the hierarchical glass transition of anisotropic colloids.

Dynamical phase transitions and their relation to structural and thermodynamic aspects of glass physics

We review recent developments in structural-dynamical phase transitions in trajectory space based on dynamic facilitation theory. An open question is how the dynamic facilitation perspective on the glass transition may be reconciled with thermodynamic theories that posit collective reorganization accompanied by a growing static length scale and, eventually, a vanishing configurational entropy. In contrast, dynamic facilitation theory invokes a dynamical phase transition between an active phase (close to the normal liquid) and an inactive phase, which is glassy and whose order parameter is either a time-averaged dynamic or structural quantity. In particular, the dynamical phase transition in systems with non-trivial thermodynamics manifests signatures of a lower critical point that lies between the mode-coupling crossover and the putative Kauzmann temperature, at which a thermodynamic phase transition to an ideal glass state would occur. We review these findings and discuss such criticality in the context of the low-temperature decrease in configurational entropy predicted by thermodynamic theories of the glass transition.

The polarizability response of a glass-forming liquid reveals intrabasin motion and interbasin transitions on a potential energy landscape

Potential energy landscape (PEL) concepts have been useful in conceptualizing the effects of intermolecular interactions on dynamic and thermodynamic properties of liquids and glasses. "Basins", or regions of reduced potential energy associated with locally preferred molecular packing are important PEL features. The molecular configurations at the bottom of these basins are referred to as inherent structures (ISs). Experimental methods for directly characterizing PEL features such as these are rare, largely relegating PEL concepts to theory and simulation studies, and impeding their exploration in real systems. Recently, we showed that quasielastic neutron scattering (QENS) data from propylene carbonate (PC) exhibit signatures of picosecond timescale motion that are consistent with intrabasin motion and interbasin transitions [Cicerone et al., J. Chem. Phys., 2017, 146, 054502]. Here we present optically-heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy studies on PC. The data exhibit signatures of motion within and transitions between basins that agree quantitatively with and extend the QENS results. We show that the librational component of the OKE response corresponds to intrabasin dynamics, and the enigmatic intermediate OKE response corresponds to interbasin transition events. The OKE data extend the measurement range of these parameters and reveal their utility in characterizing PEL features of real systems.

Over What Length Scale Does an Inorganic Substrate Perturb the Structure of a Glassy Organic Semiconductor?

While the bulk structure of vapor-deposited glasses has been extensively studied, structure at buried interfaces has received little attention, despite being important for organic electronic applications. To learn about glass structure at buried interfaces, we study the structure of vapor-deposited glasses of the organic semiconductor DSA-Ph (1,4-di-[4-(N,N-diphenyl)amino]styrylbenzene) as a function of film thickness; the structure is probed with grazing incidence X-ray scattering. We deposit on silicon and gold substrates and span a film thickness range of 10-600 nm. Our experiments demonstrate that interfacial molecular packing in vapor-deposited glasses of DSA-Ph is more disordered compared to the bulk. At a deposition temperature near room temperature, we estimate ∼8 nm near the substrate can have modified molecular packing. Molecular dynamics simulations of a coarse-grained representation of DSA-Ph reveal a similar length scale. In both the simulations and the experiments, deposition temperature controls glass structure beyond this interfacial layer of a few nanometers.

Single trajectory transport coefficients and the energy landscape by molecular dynamics simulations

The Green-Kubo (GK) method is widely used to calculate the transport coefficients of model liquids by Molecular Dynamics (MD) simulation. A reformulation of GK was proposed by Heyes et al. [J. Chem. Phys. 150, 174504 (2019)], which expressed the shear viscosity in terms of a probability distribution function (PDF) of "single trajectory (ST) viscosities," called "viscuits." This approach is extended here to the bulk viscosity, thermal conductivity, and diffusion coefficient. The PDFs of the four STs expressed in terms of their standard deviations (calculated separately for the positive and negative sides) are shown by MD to be statistically the same for the Lennard-Jones fluid. This PDF can be represented well by a sum of exponentials and is independent of system size and state point in the equilibrium fluid regime. The PDF is not well reproduced by a stochastic model. The PDF is statistically the same as that derived from the potential energy, u, and other thermodynamic quantities, indicating that the transport coefficients are determined quantitatively by and follow closely the time evolution of the underlying energy landscape. The PDFs of out-of-equilibrium supercooled high density states are quite different from those of the equilibrium states.

Understanding and Designing Thermally Activated Delayed Fluorescence Emitters: Beyond the Energy Gap Approximation

In this article recent progress in the development of molecules exhibiting Thermally Activated Delayed Fluorescence (TADF) is discussed with a particular focus upon their application as emitters in highly efficient organic light emitting diodes (OLEDs). The key aspects controlling the desirable functional properties, e. g. fast intersystem crossing, high radiative rate and unity quantum yield, are introduced with a particular focus upon the competition between the key requirements needed to achieve high performance OLEDs. The design rules required for organic and metal organic materials are discussed, and the correlation between them outlined. Recent progress towards understanding the influence of the interaction between a molecule and its environment are explained as is the role of the mechanism for excited state formation in OLEDs. Finally, all of these aspects are combined to discuss the ability to implement high level design rules for achieving higher quality materials for commercial applications. This article highlights the significant progress that has been made in recent years, but also outlines the significant challenges which persist to achieve a full understanding of the TADF mechanism and improve the stability and performance of these materials.

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.

Complex molecular dynamics of a symmetric model discotic liquid crystal revealed by broadband dielectric, thermal and neutron spectroscopy

The molecular dynamics of the triphenylene-based discotic liquid crystal HAT6 is investigated by broadband dielectric spectroscopy, advanced dynamical calorimetry and neutron scattering. Differential scanning calorimetry in combination with X-ray scattering reveals that HAT6 has a plastic crystalline phase at low temperatures, a hexagonally ordered liquid crystalline phase at higher temperatures and undergoes a clearing transition at even higher temperatures. The dielectric spectra show several relaxation processes: a localized γ-relaxation at lower temperatures and a so called α2-relaxation at higher temperatures. The relaxation rates of the α2-relaxation have a complex temperature dependence and bear similarities to a dynamic glass transition. The relaxation rates estimated by Hyper DSC, Fast Scanning calorimetry and AC Chip calorimetry have a different temperature dependence than the dielectric α2-relaxation and follow the VFT-behavior characteristic for glassy dynamics. Therefore, this process is called α1-relaxation. Its relaxation rates show a similarity with that of polyethylene. For this reason, the α1-relaxation is assigned to the dynamic glass transition of the alkyl chains in the intercolumnar space. Moreover, this process is not observed by dielectric spectroscopy, which supports its assignment. The α2-relaxation is assigned to small scale translatorial and/or small angle fluctuations of the cores. The neutron scattering data reveal two relaxation processes. The process observed at shorter relaxation times is assigned to the methyl group rotation. The second relaxation process at longer time scales agree in the temperature dependence of its relaxation rates with that of the dielectric γ-relaxation.

A simple mean-field model of glassy dynamics and glass transition

We propose a phenomenological model to describe the equilibrium dynamic behavior of amorphous glassy materials. It is assumed that a material can be represented by a lattice of cooperatively re-arranging regions (CRRs), with each CRR having two states, the low-temperature "solid" and the high-temperature "liquid". At low temperatures, the material exhibits two characteristic relaxation times, corresponding to the slow large-scale motion between the "solid" CRRs (α-relaxation) and the faster local motion within individual CRRs (β-relaxation). At high temperatures, the α- and β-relaxation times merge, as observed experimentally and suggested by the "Coupling Model" framework. Our new approach is labeled "Two-state, two (time)scale model" or TS2. It is shown that the TS2 treatment can successfully describe the "two-Arrhenius" relaxation time behavior described in several recent experiments. We also apply TS2 to describe the pressure- and molecular-weight dependence of the glass transition temperature in bulk polymers, as well as its dependence on film thickness in thin films.

Density scaling of structure and dynamics of an ionic liquid

The lines in the pressure–temperature phase diagram with constant conductivity are found to be lines where other dynamic variables as well as the molecular structure factor peak are constant, while charge ordering changes.

Tuning the melting point of selected ionic liquids through adjustment of the cation's dipole moment

In previous work with thermally robust salts [Cassity et al., Phys. Chem. Chem. Phys., 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. In this work we further test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation-anion interactions, which are recovered upon melting. If, however, the dipole-dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

Progress towards a phenomenological picture and theoretical understanding of glassy dynamics and vitrification near interfaces and under nanoconfinement

The nature of alterations to dynamics and vitrification in the nanoscale vicinity of interfaces-commonly referred to as "nanoconfinement" effects on the glass transition-has been an open question for a quarter century. We first analyze experimental and simulation results over the last decade to construct an overall phenomenological picture. Key features include the following: after a metrology- and chemistry-dependent onset, near-interface relaxation times obey a fractional power law decoupling relation with bulk relaxation; relaxation times vary in a double-exponential manner with distance from the interface, with an intrinsic dynamical length scale appearing to saturate at low temperatures; the activation barrier and vitrification temperature T_g approach bulk behavior in a spatially exponential manner; and all these behaviors depend quantitatively on the nature of the interface. We demonstrate that the thickness dependence of film-averaged T_g for individual systems provides a poor basis for discrimination between different theories, and thus we assess their merits based on the above dynamical gradient properties. Entropy-based theories appear to exhibit significant inconsistencies with the phenomenology. Diverse free-volume-motivated theories vary in their agreement with observations, with approaches invoking cooperative motion exhibiting the most promise. The elastically cooperative nonlinear Langevin equation theory appears to capture the largest portion of the phenomenology, although important aspects remain to be addressed. A full theoretical understanding requires improved confrontation with simulations and experiments that probe spatially heterogeneous dynamics within the accessible 1-ps to 1-year time window, minimal use of adjustable parameters, and recognition of the rich quantitative dependence on chemistry and interface.

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].

The Proton Density of States in Confined Water (H2O)

The hydrogen density of states (DOS) in confined water has been probed by inelastic neutron scattering spectra in a wide range of its P–T phase diagram. The liquid–liquid transition and the dynamical crossover from the fragile (super-Arrhenius) to strong (Arrhenius) glass forming behavior have been studied, by taking into account the system polymorphism in both the liquid and amorphous solid phases. The interest is focused in the low energy region of the DOS (E<10 meV) and the data are discussed in terms of the energy landscape (local minima of the potential energy) approach. In this latest research, we consider a unit scale energy (EC) linked to the water local order governed by the hydrogen bonding (HB). All the measured spectra, scaled according to such energy, evidence a universal power law behavior with different exponents (γ) in the strong and fragile glass forming regions, respectively. In the first case, the DOS data obey the Debye squared-frequency law, whereas, in the second one, we obtain a value predicted in terms of the mode-coupling theory (MCT) (γ≃1.6).

Glass polymorphism and liquid-liquid phase transition in aqueous solutions: experiments and computer simulations

One of the most intriguing anomalies of water is its ability to exist as distinct amorphous ice forms (glass polymorphism or polyamorphism). This resonates well with the possible first-order liquid-liquid phase transition (LLPT) in the supercooled state, where ice is the stable phase. In this Perspective, we review experiments and computer simulations that search for LLPT and polyamorphism in aqueous solutions containing salts and alcohols. Most studies on ionic solutes are devoted to NaCl and LiCl; studies on alcohols have mainly focused on glycerol. Less attention has been paid to protein solutions and hydrophobic solutes, even though they reveal promising avenues. While all solutions show polyamorphism and an LLPT only in dilute, sub-eutectic mixtures, there are differences regarding the nature of the transition. Isocompositional transitions for varying mole fractions are observed in alcohol but not in ionic solutions. This is because water can surround alcohol molecules either in a low- or high-density configuration whereas for ionic solutes, the water ion hydration shell is forced into high-density structures. Consequently, the polyamorphic transition and the LLPT are prevented near the ions, but take place in patches of water within the solutions. We highlight discrepancies and different interpretations within the experimental community as well as the key challenges that need consideration when comparing experiments and simulations. We point out where reinterpretation of past studies helps to draw a unified, consistent picture. In addition to the literature review, we provide original experimental results. A list of eleven open questions that need further consideration is identified.

How Tethered Probes Report the Dynamics of a Polymer near the Glass Transition

How tethered probes report dynamics of host polymers near the glass transition was investigated by changing the length of the flexible linkers and the number of tethering points via imaging rotational fluorescence correlation microscopy and compared with free probes of different sizes. The results show that tethering did not alter the temperature-dependence of polymer dynamics and the shape of the correlation decay reported by the probe; however, the rotation slowed down up to ≈1 decade when both ends of the probe were restricted with short alkyl chain linkers. Upon comparison with the bigger free probe, the mechanism of the slowdown was attributed to the restricted motion upon tethering for tethered probes compared to averaging over different regions of the dynamic heterogeneity for the bigger probe. If the size of the probe was comparable to that of the dynamic heterogeneity of the system, tethered probes accurately report dynamics relevant to glass transition, regardless of tethering conditions.

Fluid Phase Coexistence in Biological Membrane: Insights from Local Nonaffine Deformation of Lipids

Lateral heterogeneities in biomembranes play a crucial role in various physiological functions of the cell. Such heterogeneities lead to demixing of lipid constituents and formation of distinct liquid domains in the membrane. We study lateral heterogeneities in terms of topological rearrangements of lipids to identify the liquid-liquid phase coexistence in model membranes. Using ideas from the physics of amorphous systems and glasses, we calculate the degree of nonaffine deformation associated with individual lipids to characterize the liquid-ordered (L_o) and liquid-disordered (L_d) regions in model lipid bilayers. We explore the usage of this method on all-atom and coarse-grained lipid bilayer trajectories. This method is helpful in defining the instantaneous L_o-L_d domain boundaries in complex multicomponent bilayer systems. The characterization is also used to highlight the effect of line-active molecules on the phase boundaries and domain mixing. Overall, we propose a framework to explore the molecular origin of spatial and dynamical heterogeneity in biomembrane systems, which can be exploited not only in computer simulations but also in experiments.

A comparative study on the dynamic heterogeneity of supercooled polymers under nanoconfinement

Dynamic heterogeneity (DH) is a universal property of glass transition phenomena. In this work, we perform a comparative analysis of DH for pure polymer and polymer/nanoparticle composite systems in both film and bulk states via molecular dynamics simulations. We find that the dynamic gradient and the faster average dynamics due to the presence of a free surface are two leading factors, resulting from a nanoconfinement effect, which influence different parts of DH in a film system. The dynamic gradient results from differences in dynamics at different distances from the mobile surface, which induces a large deviation from the Gaussian distribution for the displacement distribution in the film. At the same time, the maximum string size which describes the region size for cooperative motion (dynamic correlation) can also be influenced by the dynamic gradient, although this influence is much weaker than that on the displacement distribution. On the other hand, reflecting temporal fluctuations of dynamics or temporal parts of DH, characteristic peak times of the non-Gaussian parameter and string size, and the ratio between persistent times and exchange times which describe the dynamic exchange properties, are mainly influenced by the faster dynamics on average. Our results demonstrate that measuring different properties (dynamic distribution, dynamic correlation or dynamic exchange) place an emphasis on distinct temporal and spatial parts of DH. It is necessary to use combinational measurements of these properties to give a complete picture of DH in nanoconfinement environments.

Density fluctuations and random walks in an overdamped and supercooled simple liquid

In this work, the short-time dynamics of simple liquid is explored both analytically and numerically with the focus on the interplay between the density fluctuations in a volume surrounding a chosen particle and its random walk motion. The particles interact via the Lennard-Jones potential with parameters corresponding to liquid argon. For large times, analytical calculations based on the fluctuation theory provides an explicit expression reproducing isothermal change of the self-diffusion coefficient in liquid argon corresponding to the experimental data. These results lead to the conclusion that such behavior is based on the reduced mobility of particles reflected in their density fluctuations that can be equivalently achieved in the cases of either low temperatures and pressures (supercooling) or moderate temperatures and high pressures (overdamping).

Liquid-like behaviours of metallic glassy nanoparticles at room temperature

Direct atomic-scale observations and measurements on dynamics of amorphous metallic nanoparticles (a-NPs) are challenging owing to the insufficient consciousness to their striking characterizations and the difficulties in technological approaches. In this study, we observe coalescence process of the a-NPs at atomic scale. We measure the viscosity of the a-NPs through the particles coalescence by in situ method. We find that the a-NPs have fast dynamics, and the viscosity of the a-NPs exhibits a power law relationship with size of the a-NPs. The a-NPs with sizes smaller than 3 nm are in a supercooled liquid state and exhibit liquid-like behaviours with a decreased viscosity by four orders of magnitude lower than that of bulk glasses. These results reveal the intrinsic flow characteristics of glasses in low demension, and pave a way to understand the liquid-like behaviours of low dimension glass, and are also of key interest to develop size-controlled nanodevices.

Nanoscale materials often exhibit size-dependent behaviour. Here, the authors use electron microscopy to quantitatively study the size-related dynamics of amorphous metallic nanoparticles, finding that particles below a critical size are in a supercooled liquid state at room temperature, with a viscosity much lower than that of bulk glasses.

Configurational entropy of glass-forming liquids

The configurational entropy is one of the most important thermodynamic quantities characterizing supercooled liquids approaching the glass transition. Despite decades of experimental, theoretical, and computational investigation, a widely accepted definition of the configurational entropy is missing, its quantitative characterization remains fraught with difficulties, misconceptions, and paradoxes, and its physical relevance is vividly debated. Motivated by recent computational progress, we offer a pedagogical perspective on the configurational entropy in glass-forming liquids. We first explain why the configurational entropy has become a key quantity to describe glassy materials, from early empirical observations to modern theoretical treatments. We explain why practical measurements necessarily require approximations that make its physical interpretation delicate. We then demonstrate that computer simulations have become an invaluable tool to obtain precise, nonambiguous, and experimentally relevant measurements of the configurational entropy. We describe a panel of available computational tools, offering for each method a critical discussion. This perspective should be useful to both experimentalists and theoreticians interested in glassy materials and complex systems.

Coupling and decoupling between translational and rotational dynamics in supercooled monodisperse soft Janus particles

We perform dynamics simulations to investigate the translational and rotational glassy dynamics in a glass-forming liquid of monodisperse soft Janus particles. We find that, with decreasing temperature, the mean-square angular displacement shows no clear plateau in the caging region, in contrast with the apparent caging behavior of translational motion. By defining a reorientational mean-square angular displacement, the caging behavior of rotational motion can be recognized. On approaching the glass transition (decreasing temperature), the coupling between translational and rotational relaxation increases, while the coupling between translational and rotational diffusion decreases, whereas the coupling between translational and reorientational diffusion increases. The strong decoupling between translational and rotational diffusion is due to the suppressed translational mobility but promoted rotational mobility of soft Janus particles. We think that the low-T SE and SED decoupling is mainly attributed to hopping motion of soft Janus particles, whereas the high-T SE and SED decoupling is mainly attributed to collective cage motion of soft Janus particles. Our results demonstrate that interaction anisotropy has a critical effect on the translational and rotational dynamics of soft Janus particles.

Thermodynamic and dynamical properties of the hard sphere system revisited by molecular dynamics simulation

Revised thermodynamic and dynamical properties of the hard sphere (HS) system are obtained from extensive molecular dynamics calculations carried out with large system sizes (number of particles, N) and long times. Accurate formulas for the compressibility factor of the HS solid and fluid branches are proposed, which represent the metastable region and take into account its divergence at close packing. Some basic second-order thermodynamic properties are obtained and a maximum in some of their derivatives in the metastable fluid region is found. The thermodynamic parameters associated with the melting-freezing transition have been determined to four digit accuracy, which generates accurate new values for the coexistence properties of the HS system. For the self-diffusion coefficient, D, it is shown that relatively large systems (N > 104) are required to achieve an accurate linear extrapolation of D to the infinite size limit with a D vs. N-1/3 plot. Moreover, it is found that there is a density dependence of the value of the slope in the linear regime. The density dependent correction becomes practically insignificant at higher densities and the hydrodynamic formula found in the literature is still accurate. However, with decreasing density the density dependence of the size correction cannot be neglected, which indicates that other sources of N-dependence, apart from those derived on purely hydrodynamic grounds, may also be important (and as yet unaccounted for). A detailed analytic representation of the density dependence of the HS self-diffusion coefficient and the HS viscosity, η, is given. It is shown that the HS viscosity near freezing and in the metastable region can be described well by the Krieger-Dougherty equation. Both D and η start to scale at high densities and in the metastable region in such a way that Dηp = const, where p ≃ 0.97, and D → 0 and η → ∞ at a packing fraction of 0.58, which coincides with some previous predictions of the HS glass transition density.

Temperature Induced Dynamical Transition of Biomolecules in Polarizable and Nonpolarizable TIP3P Water

Temperature induced dynamical transition (DT), associated with a sharp rise in molecular flexibility, is well-known to be exhibited between 270 and 280 K in glycerol to 200-230 K in hydrated biomolecules and is controlled by diffusivity (viscosity) of the solvation layer. In the molecular dynamics (MD) community, especially for water as a solvent, this has been an intense area of research despite decades of investigations. However, in general, water in these studies is described by empirical nonpolarizable force fields in which electronic polarizability is treated implicitly with effective charges and related parameters. This might have led to the present trait of discovery that DTs of biomolecules, irrespective of the potential functions for water models used, occur within a narrow band of temperature variation (30-40 K). Whereas a water molecule in a biomolecular surface and one in bulk are polarized differently, therefore explicit treatment of water polarizability would be a powerful approach toward the treatment of hydration water, believed to cause the DT manifestation. Using MD simulations, we investigated the effects of polarizable water on the DT of biomolecules and the dynamic properties of hydration water. We chose two types of solutes: globular protein (lysozyme) and more open and flexible RNAs (a hairpin and a riboswitch) with different natures of hydrophilic sites than proteins in general. We found that the characteristic temperature of DT ( T_DT) for the solutes in polarizable water is always higher than that in its nonpolarizable counterpart. In particular, for RNAs, the variations are found to be ∼45 K between the two water models, whereas for the more compact lysozyme, it is only ∼4 K. The results are discussed in light of the enormous increase in relaxation times of a liquid upon cooling in the paradigm of dynamic switchover in hydration water with liquid-liquid phase transition, derived from the existence of the second critical point. Our result supports the idea that structures of biomolecules and their interactions with the hydration water determines T_DT and provides evidence for the decisive role of polarizable water on the onset of DT, which has been hitherto ignored.

Statistics of small length scale density fluctuations in supercooled viscous liquids

Many successful theories of liquids near the melting temperature assume that small length scale density fluctuations follow Gaussian statistics. This paper presents a numerical investigation of density fluctuations in the supercooled viscous regime using an enhanced sampling method. Five model systems are investigated: the single component Lennard-Jones liquid, the Kob-Andersen binary mixture, the Wahnström binary mixture, the Lewis-Wahnström model of ortho-terphenyl, and the TIP4P/Ice model of water. The results show that the Gaussian approximation persists to a good degree into the supercooled viscous regime; however, it is less accurate at low temperatures. The analysis suggests that non-Gaussian fluctuations are related to crystalline configurations. Implications for theories of the glass transition are discussed.

Super-Arrhenius behavior of molecular glass formers

A theory is developed to calculate values of the potential-energy barriers to structural relaxation in molecular glass formers from the data of static pair-correlation function. The barrier height is shown to increase due to an increase in the number of "stable bonds" a particle forms with its neighbors and the energy of each bond as liquids move deeper into the supercooled (supercompressed) region. We present results for a system of hard spheres and compare calculated values of the structural relaxation time with experimental and simulation results.