Perspective: Supercooled liquids and glasses

In-layer inhomogeneity of molecular dynamics in quasi-liquid layers of ice

Quasi-liquid layers (QLLs) are present on the surface of ice and play a significant role in its distinctive chemical and physical properties. These layers exhibit considerable heterogeneity across different scales ranging from nanometers to millimeters. Although the formation of partially ice-like structures has been proposed, the molecular-level understanding of this heterogeneity remains unclear. Here, we examined the heterogeneity of molecular dynamics on QLLs based on molecular dynamics simulations and machine learning analysis of the simulation data. We demonstrated that the molecular dynamics of QLLs do not comprise a mixture of solid- and liquid water molecules. Rather, molecules having similar behaviors form dynamical domains that are associated with the dynamical heterogeneity of supercooled water. Nonetheless, molecules in the domains frequently switch their dynamical state. Furthermore, while there is no observable characteristic domain size, the long-range ordering strongly depends on the temperature and crystal face. Instead of a mixture of static solid- and liquid-like regions, our results indicate the presence of heterogeneous molecular dynamics in QLLs, which offers molecular-level insights into the surface properties of ice.

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

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

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

A method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob-Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

Cooperative activated hopping dynamics in binary glass-forming liquids: effects of the size ratio, composition, and interparticle interactions

Slow dynamics in supercooled and glassy liquids is an important research topic in soft matter physics. Compared to the traditionally focused one-component systems, glassy dynamics in mixture systems adds in a rich set of new complexities, which are fundamentally interesting and also relevant for many technological applications. In this paper, we apply the recently developed self-consistent cooperative hopping theory (SCCHT) to systematically investigate the effects of the size ratio, composition and interparticle interactions on the cooperative activated hopping dynamics of matrix (in larger size) and penetrant (in smaller size) particles in varied binary sphere mixture model systems, with a specific focus on ultrahigh mixture packing fractions that mimic the deeply supercooled glass transition conditions for molecular/polymeric mixture materials. Analysis shows that in these high activation barrier cases, the long-range elastic distortion associated with a matrix particle hopping over its cage confinement always generates an elastic barrier of a nonnegligible magnitude, although the ratio between the elastic barrier and local barrier contribution is sensitively dependent on all three mixture-specific system factors considered in this work. SCCHT predicts two general scenarios of penetrant-matrix cooperative activated hopping dynamics: matrix/penetrant co-hopping (regime 1) or the penetrant mean barrier hopping time shorter than that of the matrix (regime 2). Increasing the penetrant-to-matrix size ratio or the penetrant-matrix cross-attraction strength is found to universally enlarge the composition window of regime 1. Diverse dynamical properties characterising different aspects of the cooperative activated hopping process, including the penetrant and matrix transient localization lengths, penetrant and matrix hopping jump distances, different types of local and elastic activated barriers, and matrix long-time diffusivity, relaxation time and dynamic fragility are quantitatively studied against a wide range of variations over the three system factors. Of particular interest is the universal "anti-plasticization" phenomenon achievable for sufficiently strong cross-attractive interactions. The prospects this work opens for the exploration of a wide variety of polymer-based mixture materials are briefly discussed at the end.

Importance of Many Particle Correlations to the Collective Debye-Waller Factor in a Single-Particle Activated Dynamic Theory of the Glass Transition

We theoretically study the importance of many body correlations on the collective Debye-Waller (DW) factor in the context of the Nonlinear Langevin Equation (NLE) single-particle activated dynamics theory of glass transition and its extension to include collective elasticity (ECNLE theory). This microscopic force-based approach envisions structural alpha relaxation as a coupled local-nonlocal process involving correlated local cage and longer range collective barriers. The crucial question addressed here is the importance of the deGennes narrowing contribution versus a literal Vineyard approximation for the collective DW factor that enters the construction of the dynamic free energy in NLE theory. While the Vineyard-deGennes approach-based NLE theory and its ECNLE theory extension yields predictions that agree well with experimental and simulation results, use of a literal Vineyard approximation for the collective DW factor massively overpredicts the activated relaxation time. The current study suggests many particle correlations are crucial for a reliable description of activated dynamics theory of model hard sphere fluids.

Non-simple flow behavior in a polar van der Waals liquid: Structural relaxation under scrutiny

The non-exponential character of the structural relaxation is considered one of the hallmarks of the glassy dynamics, and in this context, the relatively narrow shape observed by dielectric techniques for polar glass formers has attracted the attention of the community for long time. This work addresses the phenomenology and role of specific non-covalent interactions in the structural relaxation of glass-forming liquids by the study of polar tributyl phosphate. We show that dipole interactions can couple to shear stress and modify the flow behavior, preventing the occurrence of the simple liquid behavior. We discuss our findings in the general framework of glassy dynamics and the role of intermolecular interactions.

Probing excitations and cooperatively rearranging regions in deeply supercooled liquids

Upon approaching the glass transition, the relaxation of supercooled liquids is controlled by activated processes, which become dominant at temperatures below the so-called dynamical crossover predicted by Mode Coupling theory (MCT). Two of the main frameworks rationalising this behaviour are dynamic facilitation theory (DF) and the thermodynamic scenario which give equally good descriptions of the available data. Only particle-resolved data from liquids supercooled below the MCT crossover can reveal the microscopic mechanism of relaxation. By employing state-of-the-art GPU simulations and nano-particle resolved colloidal experiments, we identify the elementary units of relaxation in deeply supercooled liquids. Focusing on the excitations of DF and cooperatively rearranging regions (CRRs) implied by the thermodynamic scenario, we find that several predictions of both hold well below the MCT crossover: for the elementary excitations, their density follows a Boltzmann law, and their timescales converge at low temperatures. For CRRs, the decrease in bulk configurational entropy is accompanied by the increase of their fractal dimension. While the timescale of excitations remains microscopic, that of CRRs tracks a timescale associated with dynamic heterogeneity, [Formula: see text]. This timescale separation of excitations and CRRs opens the possibility of accumulation of excitations giving rise to cooperative behaviour leading to CRRs.

Dynamic Heterogeneity and Kovacs' Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. We used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water's structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

Molecular Mechanism of the Debye Relaxation in Monohydroxy Alcohols Revealed from Rheo-Dielectric Spectroscopy

Rheo-dielectric spectroscopy is employed to investigate the effect of external shear on Debye-like relaxation of a model monohydroxy alcohol, i.e., the 2-ethyl-1-hexanol (2E1H). Shear deformation leads to strong acceleration in the structural relaxation, the Debye relaxation, and the terminal relaxation of 2E1H. Moreover, the shear-induced reduction in structural relaxation time, τ_{α}, scales quadratically with that of Debye time, τ_{D}, and the terminal flow time, τ_{f}, suggesting a relationship of τ_{D}^{2}∼τ_{α}. Further analyses reveal τ_{D}^{2}/τ_{α} of 2E1H follows Arrhenius temperature dependence that applies remarkably well to many other monohydroxy alcohols with different molecular sizes, architectures, and alcohol types. These results cannot be understood by the prevailing transient chain model, and suggest a H-bonding breakage facilitated sub-supramolecular reorientation as the origin of Debye relaxation of monohydroxy alcohols, akin to the molecular mechanism for the terminal relaxation of unentangled "living" polymers.

On the bridge hypothesis in the glass transition of freestanding polymer films

Freestanding thin polymer films with high molecular weights exhibit an anomalous decrease in the glass-transition temperature with film thickness. Specifically, in such materials, the measured glass-transition temperature evolves in an affine way with the film thickness, with a slope that weakly depends on the molecular weight. De Gennes proposed a sliding mechanism as the hypothetical dominant relaxation process in these systems, where stress kinks could propagate in a reptation-like fashion through so-called bridges, i.e. from one free interface to the other along the backbones of polymer macromolecules. Here, by considering the exact statistics of finite-sized random walks within a confined box, we investigate in details the bridge hypothesis. We show that the sliding mechanism cannot reproduce the basic features appearing in the experiments, and we exhibit the fundamental reasons behind such a fact.

A fractal structural feature related to dynamic crossover in metallic glass-forming liquids

The dynamic crossover in supercooled liquids initially predicted by model coupling theory has been widely accepted, but its underlying structural origin is still an open issue for glass-forming liquids. By molecular dynamics simulations of binary CuZr liquids, the present work verifies that high pressure could enhance this crossover, facilitating the studies on the structural features at the crossover temperature T_c. We discover that the topological connectivity of icosahedral clusters is responsible for this dynamic crossover, rather than all clusters. T_c is the temperature at which the connectivity degree between these clusters reaches a maximum and the dynamic heterogeneity begins to keep stable. Below T_c, the fractal topological structures appear in the medium-range order scale. The icosahedral clusters with a certain connectivity pattern can be regarded as a fractal structural unit. By employing the established fractal analysis method, the fractal dimension D of the icosahedral network is calculated. Our results indicate that the D value increases monotonically with increasing pressure and the fractal behavior of the icosahedral network is an inherent feature of metallic glasses. We also find similar fractal behavior in clusters with high local five-fold symmetry. Our findings shed light on the origin of a dynamic crossover in the deep supercooled region of metallic glasses and also demonstrate the important role of icosahedral clusters in uncovering the fractal behavior of metallic glass.

Absorption of pressurized methane in normal and supercooled p-xylene revealed via high-resolution neutron imaging

Supercooling of liquids leads to peculiarities which are scarcely studied under high-pressure conditions. Here, we report the surface tension, solubility, diffusivity, and partial molar volume for normal and supercooled liquid solutions of methane with p-xylene. Liquid bodies of perdeuterated p-xylene (p-C₈D₁₀), and, for comparison, o-xylene (o-C₈D₁₀), were exposed to pressurized methane (CH₄, up to 101 bar) at temperatures ranging 7.0-30.0 °C and observed at high spatial resolution (pixel size 20.3 μm) using a non-tactile neutron imaging method. Supercooling led to the increase of diffusivity and partial molar volume of methane. Solubility and surface tension were insensitive to supercooling, the latter substantially depended on methane pressure. Overall, neutron imaging enabled to reveal and quantify multiple phenomena occurring in supercooled liquid p-xylene solutions of methane under pressures relevant to the freeze-out in the production of liquefied natural gas.

Developments on Constitutive Material Model for Architectural Soda-Lime Silicate (SLS) Glass and Evaluation of Key Modelling Parameters

Architectural soda-lime silicate glass (SLS) is increasingly taking on complex shapes that require more detailed numerical analysis. Glass modeling is a thoroughly described topic with validated constitutive models. However, these models require a number of precise material parameters for SLS glass, and these are very sensitive to changes in glass composition. The currently available information is based on SLS glass tested in the late 1990s. As a result, most current publications are based on the above data. The object of this work was to analyze the available sources and update the information on selected key parameters for modeling. Using the currently utilized SLS glass in construction, the coefficient of thermal expansion (CTE), glass transition temperature, and the Young's modulus have been experimentally investigated. The updated material parameters will allow for more accurate modeling of the SLS glass currently used in construction, and in consequence will make the prototyping process for glass with complex geometries possible to be transferred from the production stage to the design stage, resulting in shorter production times.

Effect of the polar group content on the glass transition temperature of ROMP copolymers

Polar groups have long been recognized to greatly influence the glass transition temperature (T_g) of polymers, but understanding the underlying physical mechanism remains a challenge. Here, we study the glass formation of ring-opening metathesis polymerization (ROMP) copolymers containing polar groups by employing all-atom molecular dynamics simulations. We show that although the number of hydrogen bonds (N_HB) and the cohesive energy density increase linearly as the content of polar groups (f_pol) increases, the T_g of ROMP copolymers increases with the increase of f_pol in a nonlinear fashion, and tends to plateau for sufficiently high f_pol. Importantly, we find that the increase rate of Gibbs free energy for HB breaking gradually slows down with the increase of f_pol, indicating that the HB is gradually stabilized. Therefore, T_g is jointly determined by N_HB and the strength of HBs in the system, while the latter dominates. Although N_HB increases linearly with increasing f_pol, the HB strength increases slowly with increasing f_pol, which leads to a decreasing rate of increase in T_g.

Designing Glass and Crystalline Phases of Metal-Bis(acetamide) Networks to Promote High Optical Contrast

Owing to their high tunability and predictable structures, metal-organic materials offer a powerful platform to study glass formation and crystallization processes and to design glasses with unique properties. Here, we report a novel series of glass-forming metal-ethylenebis(acetamide) networks that undergo reversible glass and crystallization transitions below 200 °C. The glass-transition temperatures, crystallization kinetics, and glass stability of these materials are readily tunable, either by synthetic modification or by liquid-phase blending, to form binary glasses. Pair distribution function (PDF) analysis reveals extended structural correlations in both single and binary metal-bis(acetamide) glasses and highlights the important role of metal-metal correlations during structural evolution across glass-crystal transitions. Notably, the glass and crystalline phases of a Co-ethylenebis(acetamide) binary network feature a large reflectivity contrast ratio of 4.8 that results from changes in the local coordination environment around Co centers. These results provide new insights into glass-crystal transitions in metal-organic materials and have exciting implications for optical switching, rewritable data storage, and functional glass ceramics.

Two-temperature molecular dynamics simulations of crystal growth in a tungsten supercooled melt

In this work we use the two-temperature model (TTM) coupled to molecular dynamics (MD) with sinks at the boundaries of the electronic subsystem to study crystal-growth rate in a quasi-one-dimensional tungsten crystal into a supercooled melt. The possibility of varying the extension of the electronic grid along with the sinks allows a more realistic description of the electronic thermal transport away from the system, providing a considerable heat dissipation from the crystallization front. Based on this approach, our results regarding crystal-growth rates are not affected even if the size of the system is changed. Moreover, comparisons are established with respect to MD and standard TTM simulations. For these comparisons between models, something remarkable is found, and it is that the temperature and the value of the maximum growth rate are the same. In contrast, the inclusion of sinks has a great impact with respect to the standard approaches specially reflected at low temperatures, where a frustration of the liquid-crystal interface dynamics is seen until a state of zero crystal growth is reached, which is not possible to characterize quantitatively since a kind of stochastic behavior is present.

Photochemistry of pyruvic acid is governed by photo-induced intermolecular electron transfer through hydrogen bonds

Despite more than 85 years of research, the mechanism behind the photodecarboxylation of pyruvic acid remains elusive. Most studies focused on the gas and liquid phase of diluted solutions of pyruvic acid to understand the impact of sun light on the degradation of this molecule in the atmosphere. By analyzing concentrated supercooled solutions at 77 K, we demonstrate that instead of decarboxylating, the pyruvic acid molecule plays the role of electron donor and transfers an electron to an acceptor molecule that subsequently degrades to form CO₂. We show that this electron transfer occurs via hydrogen bonding and that in aqueous solutions of pyruvic acid, the hydrated form is the electron acceptor. These findings demonstrate that photo-induced electron transfer via hydrogen bonding can occur between two simple carboxylic acids and that this mechanism governs the photochemistry of pyruvic acid, providing unexplored alternative pathways for the decarboxylation of photo-inactive molecules.

When supercooled pyruvic acid is photo-irradiated, a radical detectable by ESR forms following the transfer of an electron from a molecule in its keto form to a molecule in its hydrated form. The latter subsequently degrades to CO₂ and acetic acid.

Structural Changes in Metallic Glass-Forming Liquids on Cooling and Subsequent Vitrification in Relationship with Their Properties

The present review is related to the studies of structural changes observed in metallic glass-forming liquids on cooling and subsequent vitrification in terms of radial distribution function and its analogues. These structural changes are discussed in relationship with liquid's properties, especially the relaxation time and viscosity. These changes are found to be directly responsible for liquid fragility: deviation of the temperature dependence of viscosity of a supercooled liquid from the Arrhenius equation through modification of the activation energy for viscous flow. Further studies of this phenomenon are necessary to provide direct mathematical correlation between the atomic structure and properties.

Annealing glasses by cyclic shear deformation

A major challenge in simulating glassy systems is the ability to generate configurations that may be found in equilibrium at sufficiently low temperatures, in order to probe static and dynamic behavior close to the glass transition. A variety of approaches have recently explored ways of surmounting this obstacle. Here, we explore the possibility of employing mechanical agitation, in the form of cyclic shear deformation, to generate low energy configurations in a model glass former. We perform shear deformation simulations over a range of temperatures, shear rates, and strain amplitudes. We find that shear deformation induces faster relaxation toward low energy configurations, or overaging, in simulations at sufficiently low temperatures, consistently with previous results for athermal shear. However, for temperatures at which simulations can be run until a steady state is reached with or without shear deformation, we find that the inclusion of shear deformation does not result in any speed up of the relaxation toward low energy configurations. Although we find the configurations from shear simulations to have properties indistinguishable from an equilibrium ensemble, the cyclic shear procedure does not guarantee that we generate an equilibrium ensemble at a desired temperature. In order to ensure equilibrium sampling, we develop a hybrid Monte Carlo algorithm that employs cyclic shear as a trial generation step and has acceptance probabilities that depend not only on the change in internal energy but also on the heat dissipated (equivalently, work done). We show that such an algorithm, indeed, generates an equilibrium ensemble.

Evidence suggesting kinetic unfreezing of water mobility in two distinct processes in pressure-amorphized clathrate hydrates

Type II clathrate hydrates (CHs) with tetrahydrofuran (THF), cyclobutanone (CB) or 1,3-dioxolane (DXL) guest molecules collapse to an amorphous state near 1 GPa on pressurization below 140 K. On subsequent heating in the 0.2-0.7 GPa range, thermal conductivity and heat capacity results of the homogeneous amorphous solid show two glass transitions, first a thermally weak glass transition, GT1, near 130 K; thereafter a thermally strong glass transition, GT2, which implies a transformation to an ultraviscous liquid on heating. Here we compare the GTs of normal and deuterated samples and samples with different guest molecules. The results show that GT1 and GT2 are unaffected by deuteration of the THF guest and exchange of THF with CB or DXL, whereas the glass transition temperatures (T_gs) shift to higher temperatures on deuteration of water; T_g of GT2 increases by 2.5 K. These results imply that both GTs are associated with the water network. This is corroborated by the fact that GT2 is detected only in the state which is the amorphized CH's counterpart of expanded high density amorphous ice. The results suggest a rare transition sequence of an orientational glass transition followed by a glass to liquid transition, i.e., kinetic unfreezing of H₂O reorientational and translational mobility in two distinct processes.

A density scaling conjecture for aging glasses

The aging rate of glasses has traditionally been modeled as a function of temperature, T , andfictive temperature, while density, ρ, is not explicitly included as a parameter. However, this de-scription does not naturally connect to the modern understanding of what governs the relaxationrate in equilibrium. In equilibrium it is well known that the relaxation rate, γeq , depends on tem-perature and density. In addition a large class of systems obey density scaling which means therate specifically depends on the scaling parameter, Γ = e(ρ)/T , where e(ρ) is a system specificfunction. This paper present a generalization of the fictive temperature concept in terms of a fic-tive scaling paramter, Γfic , and a density scaling conjecture for aging glasses in which the agingrate depends on Γ and Γfic .

Vapor-Deposited Thin Films: Studying Crystallization and α-relaxation Dynamics of the Molecular Drug Celecoxib

Crystallization is one of the major challenges in using glassy solids for technological applications. Considering pharmaceutical drugs, maintaining a stable amorphous form is highly desirable for improved solubility. Glasses prepared by the physical vapor deposition technique got attention because they possess very high stability, taking thousands of years for an ordinary glass to achieve. In this work, we have investigated the effect of reducing film thickness on the α-relaxation dynamics and crystallization tendency of vapor-deposited films of celecoxib (CXB), a pharmaceutical substance. We have scrutinized its crystallization behavior above and below the glass-transition temperature (T_g). Even though vapor deposition of CXB cannot inhibit crystallization completely, we found a significant decrease in the crystallization rate with decreasing film thickness. Finally, we have observed striking differences in relaxation dynamics of vapor-deposited thin films above the T_g compared to spin-coated counterparts of the same thickness.

Local-Average Free Volume Correlates with Dynamics in Glass Formers

Glass formers exhibit a pronounced slowdown in dynamics, accompanied by progressive heterogeneity as they approach the glass transition. There is intense debate over whether the dramatic slowdown is caused by dynamical heterogeneity and whether the enhanced dynamical heterogeneity originates from structural causes. However, the connection between dynamical heterogeneity and the spatial distribution of the single-particle free volume (a purely static structural quantity) was found to be rather weak, which raises the question of whether dynamic heterogeneity has a purely structural origin. Here, by introducing the concept of local-average free volume, we present numerical evidence that long-time dynamic heterogeneity shows significantly enhanced correlation with the average local free volume over a length scale of a few neighboring shells. Our results resolve the long-standing controversy about whether free volume plays an important role in particle rearrangements associated with the activated hopping relaxation. The concept of "local average" can be applied to other local structural descriptors to better correlate with dynamic heterogeneity in glass-forming liquids.

Universal scaling law of glass rheology

The similarity in atomic/molecular structure between liquids and glasses has stimulated a long-standing hypothesis that the nature of glasses may be more fluid-like, rather than the apparent solid. In principle, the nature of glasses can be characterized by the dynamic response of their rheology in a wide rate range, but this has not been realized experimentally, to the best of our knowledge. Here we report the dynamic response of shear stress to the shear strain rate of metallic glasses over a timescale of nine orders of magnitude, equivalent to hundreds of years, by broadband stress relaxation experiments. The dynamic response of the metallic glasses, together with other 'glasses', follows a universal scaling law within the framework of fluid dynamics. The universal scaling law provides comprehensive validation of the conjecture on the jamming (dynamic) phase diagram by which the dynamic behaviours of a wide variety of 'glasses' can be unified under one rubric parameterized by the thermodynamic variables of temperature, volume and stress in the trajectory space.

Intrinsic viscuit probability distribution functions for transport coefficients of liquids and solids

A reformulation of the Green–Kubo expressions for the transport coefficients of liquids in terms of a probability distribution function (PDF) of short trajectory contributions, which were named “viscuits,” has been explored in a number of recent publications. The viscuit PDF, P, is asymmetric on the two sides of the distribution. It is shown here using equilibrium 3D and 2D molecular dynamics simulations that the viscuit PDF of a range of simple molecular single component and mixture liquid and solid systems can be expressed in terms of the same intrinsic PDF ( P0), which is derived from P with the viscuit normalized by the standard deviation separately on each side of the distribution. P0 is symmetric between the two sides and can be represented for not very small viscuit values by the same gamma distribution formulated in terms of a single disposable parameter. P0 tends to an exponential in the large viscuit wings. Scattergrams of the viscuits and their associated single trajectory correlation functions are shown to distinguish effectively between liquids, solids, and glassy systems. The so-called viscuit square root method for obtaining the transport coefficients is shown to be a useful probe of small and statistically zero self-diffusion coefficients of molecules in the liquid and solid states, respectively. The results of this work suggest that the transport coefficients have a common underlying physical origin, reflecting at a coarse-grained level the traversal statistics of the system through its high-dimensioned potential energy landscape.

How do interfaces alter the dynamics of supercooled water?

The structure of liquid water in the proximity of an interface can deviate significantly from that of bulk water, with surface-induced structural perturbations typically converging to bulk values at about ∼1 nm from the interface. While these structural changes are well established it is, in contrast, less clear how an interface perturbs the dynamics of water molecules within the liquid. Here, through an extensive set of molecular dynamics simulations of supercooled bulk and interfacial water films and nano-droplets, we observe the formation of persistent, spatially extended dynamical domains in which the average mobility varies as a function of the distance from the interface. This is in stark contrast with the dynamical heterogeneity observed in bulk water, where these domains average out spatially over time. We also find that the dynamical response of water to an interface depends critically on the nature of the interface and on the choice of interface definition. Overall these results reveal a richness in the dynamics of interfacial water that opens up the prospect of tuning the dynamical response of water through specific modifications of the interface structure or confining material.

Understanding, quantifying, and controlling the molecular ordering of semiconducting polymers: from novices to experts and amorphous to perfect crystals

Molecular packing and texture of semiconducting polymers are often critical to the performance of devices using these materials. Although frameworks exist to quantify the ordering, interpretations are often just qualitative, resulting in imprecise use of terminology. Here, we reemphasize the significance of quantifying molecular ordering in terms of degree of crystallinity (volume fractions that are ordered) and quality of ordering and their relation to the size scale of an ordered region. We are motivated in part by our own imprecise and inconsistent use of terminology in the past, as well as the need to have a primer or tutorial reference to teach new group members. We strive to develop and use consistent terminology with regards to crystallinity, semicrystallinity, paracrystallinity, and related characteristics. To account for vastly different quality of ordering along different directions, we classify paracrystals into 2D and 3D paracrystals and use paracrystallite to describe the spatial extent of molecular ordering in 1-10 nm. We show that a deeper understanding of molecular ordering can be achieved by combining grazing-incidence wide-angle X-ray scattering and differential scanning calorimetry, even though not all aspects of these measurements are consistent, and some classification appears to be method dependent. We classify a broad range of representative polymers under common processing conditions into five categories based on the quantitative analysis of the paracrystalline disorder parameter (g) and thermal transitions. A small database is presented for 13 representative conjugated and insulating polymers ranging from amorphous to semi-paracrystalline. Finally, we outline the challenges to rationally design more perfect polymer crystals and propose a new molecular design approach that envisions conceptual molecular grafting that is akin to strained and unstrained hetero-epitaxy in classic (compound) semiconductors thin film growth.

The energy landscapes of bidisperse particle assemblies on a sphere

The interplay between crystalline ordering, curvature, and size dispersity make the packing of bidisperse mixtures of particles on a sphere a varied and complex phenomenon. These structures have functional significance in a broad range of systems, such as cellular organisation in spherical epithelia, catalytic activity in binary colloidosomes, and chemical activity in heterofullerenes. In this contribution, we elucidate the potential energy landscapes for systems of repulsive, bidisperse particles confined to the surface of a sphere. It is commonly asserted that particle size dispersity destroys ordered arrangements, leading to glassy landscapes. Surprisingly, across a range of compositions, we find highly ordered global minima. Moreover, a minority of small particles is able to passivate defects, stabilising bidisperse global minima relative to monodisperse systems. However, our landscape analysis also reveals that bidispersity introduces numerous defective, low-lying states that are expected to cause broken ergodicity in corresponding experimental and computational systems. Probing the global minimum structures further, particle segregation is energetically preferred at intermediate compositions, contrasting with the approximate icosahedral global packing at either end of the composition range. Finally, changing the composition has a dramatic effect on the heat capacity: systems with low-symmetry global minima have melting temperatures an order of magnitude lower than monodisperse or high-symmetry systems. This observation may provide a further example of the principle of maximum symmetry: higher symmetry global minima exhibit a larger energy separation from the minima that define the high-entropy phase-like region of configuration space, raising the transition temperature.

Impact of Low Concentration of Strongly Hydrogen-Bonded Water Molecules on the Dynamics of Amorphous Terfenadine: Insights from Molecular Dynamics Simulations and Dielectric Relaxation Spectroscopy

The impact of low water concentration of strongly hydrogen-bonded water molecules on the dynamical properties of amorphous terfenadine (TFD) is investigated through complementary molecular dynamics (MD) simulations and dielectric relaxation spectroscopy (DRS) experiments. In this article, we especially highlight the important role played by some residual water molecules in the concentration of 1-2% (w/w) trapped in the TFD glassy matrix, which are particularly difficult to remove experimentally without a specific heating/drying process. From MD computations and analyses of the hydrogen bonding (HB) interactions, different categories of water molecules are revealed and particularly the presence of strongly HB water molecules. These latter localize themselves in small pockets in empty spaces existing in between the TFD molecules due to the poor packing of the glassy state and preferentially interact with the polar groups close to the flexible central part of the TFD molecules. We present a simple model which rationalizes at the molecular scale the effect of these strongly HB water molecules on dynamics and how they give rise to a supplementary relaxation process (namely process S) which is detected for the first time in the glassy state of TFD annealed at room temperature while this process is completely absent in a non-annealed glass. It also explains how this supplementary relaxation is coupled with the intramolecular motion (namely process γ) of the very flexible central part of the TFD molecule. The present findings help to understand more generally the microscopic origin of the secondary relaxations often detected by DRS in the glassy states of molecular compounds for which the exact nature is still debated.

Non-hyperuniform metastable states around a disordered hyperuniform state of densely packed spheres: stochastic density functional theory at strong coupling

The disordered and hyperuniform structures of densely packed spheres near and at jamming are characterized by vanishing of long-wavelength density fluctuations, or equivalently by long-range power-law decay of the direct correlation function (DCF). We focus on previous simulation results that exhibit the degradation of hyperuniformity in jammed structures while maintaining the long-range nature of the DCF to a certain length scale. Here we demonstrate that the field-theoretic formulation of stochastic density functional theory is relevant to explore the degradation mechanism. The strong-coupling expansion method of stochastic density functional theory is developed to obtain the metastable chemical potential considering the intermittent fluctuations in dense packings. The metastable chemical potential yields the analytical form of the metastable DCF that has a short-range cutoff inside the sphere while retaining the long-range power-law behavior. It is confirmed that the metastable DCF provides the zero-wavevector limit of the structure factor in quantitative agreement with the previous simulation results of degraded hyperuniformity. We can also predict the emergence of soft modes localized at the particle scale by plugging this metastable DCF into the linearized Dean-Kawasaki equation, a stochastic density functional equation.

Theory of structural relaxation in glass from the thermodynamics of irreversible processes

This work proposes a fundamental thermodynamic description of structural relaxation in glasses by establishing a link between the Prony series solution to volume relaxation derived from the principles of irreversible thermodynamics and asymmetric Lévy stable distribution of relaxation rates. Additionally, it is shown that the bulk viscosity of glass, and not the shear viscosity, is the transport coefficient governing structural relaxation. We also report the distribution of relaxation times and energy barrier heights underpinning stretched exponential relaxation. It is proposed that this framework may be used for qualitative and quantitative descriptions of the relaxation kinetics in glass.

Reorientational dynamics in highly asymmetric binary low-molecular mixtures-A quantitative comparison of dielectric and NMR spectroscopy results

Previously, we scrutinized the dielectric spectra of a binary glass former made by a low-molecular high-T_g component 2-(m-tertbutylphenyl)-2'-tertbutyl-9,9'-spirobi[9H]fluorene (m-TPTS; T_g = 350 K) and low-T_g tripropyl phosphate (TPP; T_g = 134 K) [Körber et al., Phys. Chem. Chem. Phys. 23, 7200 (2021)]. Here, we analyze nuclear magnetic resonance (NMR) spectra and stimulated echo decays of deuterated m-TPTS-d₄ (²H) and TPP (³¹P) and attempt to understand the dielectric spectra in terms of component specific dynamics. The high-T_g component (α₁) shows relaxation similar to that of neat systems, yet with some broadening upon mixing. This correlates with high-frequency broadening of the dielectric spectra. The low-T_g component (α₂) exhibits highly stretched relaxations and strong dynamic heterogeneities indicated by "two-phase" spectra, reflecting varying fractions of fast and slow liquid-like reorienting molecules. Missing for the high-T_g component, such two-phase spectra are identified down to w_TPP = 0.04, indicating that isotropic reorientation prevails in the rigid high-T_g matrix stretching from close to T_g ^TPP to T_g1 w_TPP. This correlates with low-frequency broadening of the dielectric spectra. Two T_g values are defined: T_g1 (w_TPP) displays a plasticizer effect, whereas T_g2 (w_TPP) passes through a maximum, signaling extreme separation of the component dynamics at low w_TPP. We suggest understanding the latter counter-intuitive feature by referring to a crossover from "single glass" to "double glass" scenario revealed by recent MD simulations. Analyses reveal that a second population of TPP molecules exists, which is associated with the dynamics of the high-T_g component. However, the fractions are lower than suggested by the dielectric spectra. We discuss this discrepancy considering the role of collective dynamics probed by dielectric but not by NMR spectroscopy.

Decoupling between calorimetric and dynamical glass transitions in high-entropy metallic glasses

Glass transition is one of the unresolved critical issues in solid-state physics and materials science, during which a viscous liquid is frozen into a solid or structurally arrested state. On account of the uniform arrested mechanism, the calorimetric glass transition temperature (Tg) always follows the same trend as the dynamical glass transition (or α-relaxation) temperature (Tα) determined by dynamic mechanical analysis (DMA). Here, we explored the correlations between the calorimetric and dynamical glass transitions of three prototypical high-entropy metallic glasses (HEMGs) systems. We found that the HEMGs present a depressed dynamical glass transition phenomenon, i.e., HEMGs with moderate calorimetric Tg represent the highest Tα and the maximum activation energy of α-relaxation. These decoupled glass transitions from thermal and mechanical measurements reveal the effect of high configurational entropy on the structure and dynamics of supercooled liquids and metallic glasses, which are associated with sluggish diffusion and decreased dynamic and spatial heterogeneities from high mixing entropy. The results have important implications in understanding the entropy effect on the structure and properties of metallic glasses for designing new materials with plenteous physical and mechanical performances.

From Lab Formulation Development to CTM Manufacturing of KO-947 Injectable Drug Products: a Case Study and Lessons Learned

Formulation development of KO-947-K mesylate injectable drug products was described. Solution formulations were initially attempted, and key parameters such as drug concentration, buffer, pH, complexing agent, and tonicity modifying agent were carefully evaluated in the lab setting, mainly focusing on solubility and chemical stability. A lead solution formulation was advanced to a scaleup campaign. An unexpected stability issue was encountered, and the root cause was attributed to the heterogeneous liquid freezing process of the formulated solution at -20°C, which had not been captured in the lab setting. A lyophilized product was then designed to overcome the issue and supplied to the phase I clinical trial.

The effects of vacancies and their mobility on the dynamic heterogeneity in 1,3-dimethylimidazolium hexafluorophosphate organic ionic plastic crystals

Organic ionic plastic crystals (OIPCs) are the crystals of electrolytes with a long-range translational order. The rotational modes of ions in OIPCs are, however, activated even in solid phases such that the diffusion of dopants such as lithium ions may be facilitated. OIPCs have been, therefore, considered as good candidates for solid electrolytes. Recent experiments and theoretical studies have suggested that both the translational and the rotational diffusion of ions are quite heterogeneous: the diffusion of some ions are quite fast while other ions of the same kind hardly diffuse, either rotationally or translationally. Such dynamic heterogeneity would be a key to the transport mechanism of dopants in solid state electrolytes. In this work, we investigate the effects of defects on the dynamic heterogeneity of OIPCs. We perform atomistic molecular dynamics simulation of 1,3-dimethylimidazolium hexafluorophosphate ([MMIM][PF6]) with a pair of cation and anion vacancies. At low temperature, vacancies undergo hopping motions toward each other and form a charge-neutral cluster. At high temperature, two vacancies act like a loosely bonded molecule and diffuse together via hopping motions. We find that the translational diffusion of ions is correlated strongly with the vacancy diffusion and becomes heterogeneous when the vacancies hop. The rotation of ions also becomes active when the ions are close to vacancies such that the rotational dynamic heterogeneity strengthens.

Introducing Memory in Coarse-Grained Molecular Simulations

Preserving the correct dynamics at the coarse-grained (CG) level is a pressing problem in the development of systematic CG models in soft matter simulation. Starting from the seminal idea of simple time-scale mapping, there have been many efforts over the years toward establishing a meticulous connection between the CG and fine-grained (FG) dynamics based on fundamental statistical mechanics approaches. One of the most successful attempts in this context has been the development of CG models based on the Mori-Zwanzig (MZ) theory, where the resulting equation of motion has the form of a generalized Langevin equation (GLE) and closely preserves the underlying FG dynamics. In this Review, we describe some of the recent studies in this regard. We focus on the construction and simulation of dynamically consistent systematic CG models based on the GLE, both in the simple Markovian limit and the non-Markovian case. Some recent studies of physical effects of memory are also discussed. The Review is aimed at summarizing recent developments in the field while highlighting the major challenges and possible future directions.

Experimental test of a predicted dynamics-structure-thermodynamics connection in molecularly complex glass-forming liquids

Understanding in a unified manner the generic and chemically specific aspects of activated dynamics in diverse glass-forming liquids over 14 or more decades in time is a grand challenge in condensed matter physics, physical chemistry, and materials science and engineering. Large families of conceptually distinct models have postulated a causal connection with qualitatively different "order parameters" including various measures of structure, free volume, thermodynamic properties, short or intermediate time dynamics, and mechanical properties. Construction of a predictive theory that covers both the noncooperative and cooperative activated relaxation regimes remains elusive. Here, we test using solely experimental data a recent microscopic dynamical theory prediction that although activated relaxation is a spatially coupled local-nonlocal event with barriers quantified by local pair structure, it can also be understood based on the dimensionless compressibility via an equilibrium statistical mechanics connection between thermodynamics and structure. This prediction is found to be consistent with observations on diverse fragile molecular liquids under isobaric and isochoric conditions and provides a different conceptual view of the global relaxation map. As a corollary, a theoretical basis is established for the structural relaxation time scale growing exponentially with inverse temperature to a high power, consistent with experiments in the deeply supercooled regime. A criterion for the irrelevance of collective elasticity effects is deduced and shown to be consistent with viscous flow in low-fragility inorganic network-forming melts. Finally, implications for relaxation in the equilibrated deep glass state are briefly considered.

Dimensionality dependence of the Kauzmann temperature: A case study using bulk and confined water

The Kauzmann temperature (T_K) of a supercooled liquid is defined as the temperature at which the liquid entropy becomes equal to that of the crystal. The excess entropy, the difference between liquid and crystal entropies, is routinely used as a measure of the configurational entropy, whose vanishing signals the thermodynamic glass transition. The existence of the thermodynamic glass transition is a widely studied subject, and of particular recent interest is the role of dimensionality in determining the presence of a glass transition at a finite temperature. The glass transition in water has been investigated intensely and is challenging as the experimental glass transition appears to occur at a temperature where the metastable liquid is strongly prone to crystallization and is not stable. To understand the dimensionality dependence of the Kauzmann temperature in water, we study computationally bulk water (three-dimensions), water confined in the slit pore of the graphene sheet (two-dimensions), and water confined in the pore of the carbon nanotube of chirality (11,11) having a diameter of 14.9 Å (one-dimension), which is the lowest diameter where amorphous water does not always crystallize into nanotube ice in the supercooled region. Using molecular dynamics simulations, we compute the entropy of water in bulk and under reduced dimensional nanoscale confinement to investigate the variation of the Kauzmann temperature with dimension. We obtain a value of T_K (133 K) for bulk water in good agreement with experiments [136 K (C. A. Angell, Science 319, 582-587 (2008) and K. Amann-Winkel et al., Proc. Natl. Acad. Sci. U. S. A. 110, 17720-17725 (2013)]. However, for confined water, in two-dimensions and one-dimension, we find that there is no finite temperature Kauzmann point (in other words, the Kauzmann temperature is 0 K). Analysis of the fluidicity factor, a measure of anharmonicity in the oscillation of normal modes, reveals that the Kauzmann temperature can also be computed from the difference in the fluidicity factor between amorphous and ice phases.

Structural relaxation and crystallization in supercooled water from 170 to 260 K

The origin of water's anomalous properties has been debated for decades. Resolution of the problem is hindered by a lack of experimental data in a crucial region of temperatures, T, and pressures where supercooled water rapidly crystallizes-a region often referred to as "no man's land." A recently developed technique where water is heated and cooled at rates greater than 10⁹ K/s now enables experiments in this region. Here, it is used to investigate the structural relaxation and crystallization of deeply supercooled water for 170 K < T < 260 K. Water's relaxation toward a new equilibrium structure depends on its initial structure with hyperquenched glassy water (HQW) typically relaxing more quickly than low-density amorphous solid water (LDA). For HQW and T > 230 K, simple exponential relaxation kinetics is observed. For HQW at lower temperatures, increasingly nonexponential relaxation is observed, which is consistent with the dynamics expected on a rough potential energy landscape. For LDA, approximately exponential relaxation is observed for T > 230 K and T < 200 K, with nonexponential relaxation only at intermediate temperatures. At all temperatures, water's structure can be reproduced by a linear combination of two, local structural motifs, and we show that a simple model accounts for the complex kinetics within this context. The relaxation time, τ _rel , is always shorter than the crystallization time, τ _xtal For HQW, the ratio, τ _xtal /τ _rel , goes through a minimum at ∼198 K where the ratio is about 60.

Molecular heterogeneity in aqueous cosolvent systems

Aqueous cosolvent systems (ACoSs) are mixtures of small polar molecules such as amides, alcohols, dimethyl sulfoxide, or ions in water. These liquids have been the focus of fundamental studies due to their complex intermolecular interactions as well as their broad applications in chemistry, medicine, and materials science. ACoSs are fully miscible at the macroscopic level but exhibit nanometer-scale spatial heterogeneity. ACoSs have recently received renewed attention within the chemical physics community as model systems to explore the relationship between intermolecular interactions and microscopic liquid-liquid phase separation. In this perspective, we provide an overview of ACoS spatial segregation, dynamic heterogeneity, and multiscale relaxation dynamics. We describe emerging approaches to characterize liquid microstructure, H-bond networks, and dynamics using modern experimental tools combined with molecular dynamics simulations and network-based analysis techniques.

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.