Elastically Collective Nonlinear Langevin Equation Theory of Glass-Forming Liquids: Transient Localization, Thermodynamic Mapping, and Cooperativity

Effect of the nature of the solid substrate on spatially heterogeneous activated dynamics in glass forming supported films

We extend the force-level elastically collective nonlinear Langevin equation theory to treat the spatial gradients of the alpha relaxation time and glass transition temperature, and the corresponding film-averaged quantities, to the geometrically asymmetric case of finite thickness supported films with variable fluid-substrate coupling. The latter typically nonuniversally slows down motion near the solid-liquid interface as modeled via modification of the surface dynamic free energy caging constraints that are spatially transferred into the film and which compete with the accelerated relaxation gradient induced by the vapor interface. Quantitative applications to the foundational hard sphere fluid and a polymer melt are presented. The strength of the effective fluid-substrate coupling has very large consequences for the dynamical gradients and film-averaged quantities in a film thickness and thermodynamic state dependent manner. The interference of the dynamical gradients of opposite nature emanating from the vapor and solid interfaces is determined, including the conditions for the disappearance of a bulk-like region in the film center. The relative importance of surface-induced modification of local caging vs the generic truncation of the long range collective elastic component of the activation barrier is studied. The conditions for the accuracy and failure of a simple superposition approximation for dynamical gradients in thin films are also determined. The emergence of near substrate dead layers, large gradient effects on film-averaged response functions, and a weak non-monotonic evolution of dynamic gradients in thick and cold films are briefly discussed. The connection of our theoretical results to simulations and experiments is briefly discussed, as is the extension to treat more complex glass-forming systems under nanoconfinement.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

Combined description of pressure-volume-temperature and dielectric relaxation of several polymeric and low-molecular-weight organic glass-formers using SL-TS2 approach

We apply our recently-developed mean-field "SL-TS2" (two-state Sanchez-Lacombe) model to simultaneously describe dielectric α-relaxation time, τ_α, and pressure-volume-temperature (PVT) data in four polymers (polystyrene, poly(methylmethacrylate), poly(vinyl acetate) and poly(cyclohexane methyl acrylate)) and four organic molecular glass formers (ortho-terphenyl, glycerol, PCB-62, and PDE). Previously, it has been shown that for all eight materials, the Casalini-Roland thermodynamical scaling, τ_α = f(Tvγsp) (where T is temperature and v_sp is specific volume) is satisfied (R. Casalini and C. M. Roland, Phys. Rev. E, 2004, 69(6), 62501). It has also been previously shown that the same scaling emerges naturally (for sufficiently low pressures) within the "SL-TS2" framework (V. V. Ginzburg, Soft Matter, 2021, 17, 9094-9106). Here, we fit the ambient pressure curves for the relaxation time and the specific volume as functions of temperature for the eight materials and observe a good agreement between theory and experiment. We then use the Casalini-Roland scaling to convert those results into "master curves", thus enabling predictions of relaxation times and specific volumes at elevated pressures. The proposed approach can be used to describe other glass-forming materials, both low-molecular-weight and polymeric.

Approaching Polymer Dynamics Combining Artificial Neural Networks and Elastically Collective Nonlinear Langevin Equation

The analysis of structural relaxation dynamics of polymers gives an insight into their mechanical properties, whose characterization is used to qualify a given material for its practical scope. The dynamics are usually expressed in terms of the temperature dependence of the relaxation time, which is only available through time-consuming experimental processes following polymer synthesis. However, it would be advantageous to estimate their dynamics before synthesizing them when designing new materials. In this work, we propose a combined approach of artificial neural networks and the elastically collective nonlinear Langevin equation (ECNLE) to estimate the temperature dependence of the main structural relaxation time of polymers based only on the knowledge of the chemical structure of the corresponding monomer.

Confinement Effects on the Spatially Inhomogeneous Dynamics in Metallic Glass Films

This work develops the elastically collective nonlinear Langevin equation theory to investigate, for the first time, the glassy dynamics in capped metallic glass thin films. Finite-size effects on the spatial gradient of structural relaxation time and glass transition temperature (T_g) are calculated at different temperatures and vitrification criteria. Molecular dynamics is significantly slowed down near rough solid surfaces, and the dynamics at location far from the interfaces is sped up. In thick films, the mobility gradient normalized by the bulk value obeys the double-exponential form since interference effects between two surfaces are weak. Reducing the film thickness induces a strong dynamic coupling between two surfaces and flattens the relaxation gradient. The normalized gradient of the glass transition temperature is independent of vitrification time scale criterion and can be fitted by a superposition function as the films are not ultrathin. The local fragility is found to remain unchanged with location. This finding suggests that one can use Angell plots of bulk relaxation time and the T_g spatial gradient to characterize glassy dynamics in metallic glass films. Our computational results agree well with experimental data and simulation.

Theory of Pressure-Induced Rejuvenation and Strain Hardening in Metallic Glasses

We theoretically investigate high-pressure effects on the atomic dynamics of metallic glasses. The theory predicts compression-induced rejuvenation and the resulting strain hardening that have been recently observed in metallic glasses. Structural relaxation under pressure is mainly governed by local cage dynamics. The external pressure restricts the dynamical constraints and slows down the atomic mobility. In addition, the compression induces a rejuvenated metastable state (local minimum) at a higher energy in the free-energy landscape. Thus, compressed metallic glasses can rejuvenate and the corresponding relaxation is reversible. This behavior leads to strain hardening in mechanical deformation experiments. Theoretical predictions agree well with experiments.

Determination of Young's Modulus of Active Pharmaceutical Ingredients by Relaxation Dynamics at Elevated Pressures

A new approach is theoretically proposed to study the glass transition of active pharmaceutical ingredients and a glass-forming anisotropic molecular liquid at high pressures. We describe amorphous materials as a fluid of hard spheres. Effects of nearest neighbor interactions and cooperative motions of particles on glassy dynamics are quantified through a local and collective elastic barrier calculated using the elastically collective nonlinear Langevin equation theory. Inserting two barriers into Kramer's theory gives the structural relaxation time. Then, we formulate a new mapping based on the thermal expansion process under pressure to intercorrelate particle density, temperature, and pressure. This analysis allows us to determine the pressure and temperature dependence of α relaxation. From this, we estimate the effective elastic modulus of amorphous materials and capture the effects of conformation on the relaxation process. Remarkably, our theoretical results agree well with experiments.

Theoretical and Experimental Study of Compression Effects on Structural Relaxation of Glass-Forming Liquids

We develop the elastically collective nonlinear Langevin equation theory of bulk relaxation of glass-forming liquids to investigate molecular mobility under compression conditions. The applied pressure restricts more molecular motion and therefore significantly slows down the molecular dynamics when increasing the pressure. We quantitatively determine the temperature and pressure dependence of the structural relaxation time. To validate our model, dielectric spectroscopy experiments for three rigid and nonpolymeric supramolecules are carried out at ambient and elevated pressures. The numerical results quantitatively agree with experimental data.

Theory of Structural and Secondary Relaxation in Amorphous Drugs under Compression

Compression effects on alpha and beta relaxation process of amorphous drugs are theoretically investigated by developing the elastically collective nonlinear Langevin equation theory. We describe the structural relaxation as a coupling between local and nonlocal activated process. Meanwhile, the secondary beta process is mainly governed by the nearest-neighbor interactions of a molecule. This assumption implies the beta relaxation acts as a precursor of the alpha relaxation. When external pressure is applied, a small displacement of a molecule is additionally exerted by a pressure-induced mechanical work in the dynamic free energy, which quantifies interactions between a molecule with its nearest neighbors. The local dynamics has more restriction and it induces stronger effects of collective motions on single-molecule dynamics. Thus, the alpha and beta relaxation times are significantly slowed down with increasing compression. We apply this approach to determine the temperature and pressure dependence of the alpha and beta relaxation time for curcumin, glibenclamide, and indomethacin, and compare numerical results with prior experimental studies. Both qualitative and quantitative agreement between theoretical calculations and experiments validate our assumptions and reveal their limitations. Our approach would pave the way for the development of the drug formulation process.

Chemical Langevin Equation for Complex Reactions

The stochastic dynamical behaviors of an elementary reaction system can be investigated by the chemical Langevin equation (CLE). However, most of the reactions in engineering belong to complex reactions. It is not appropriate to describe the random evolution process of a complex chemical reaction system by directly using CLE because the foundation of the deviation of CLE is the equilibrium equation of the number of molecules in elementary reaction systems. In the study, the chemical Langevin equation for complex reactions (CLE-CR) is proposed based on the random process theory by introducing the extent of reactions to express the reaction rates of complex reactions. The reaction rates of complex reactions are regarded as some random variables following Poisson distribution. To illustrate the essential consistency of CLE-CR and CLE, the physical meaning of the propensity function in CLE is comprehensively discussed. A numerical example from chemical engineering is employed to demonstrate the effectiveness of CLE-CR and the solving procedure. The results show that CLE-CR can be conveniently applied into engineering to investigate the stochastic dynamical behaviors of complex reaction systems, giving the probabilistic information of the concentration evolution of chemical constituents.

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.

Do String-like Cooperative Motions Predict Relaxation Times in Glass-Forming Liquids?

The Adam-Gibbs theory of glass formation posits that the growth in the activation barrier of fragile liquids on cooling emerges from a loss of configurational entropy and concomitant growth in "cooperatively rearranging regions" (CRRs). A body of literature over 2 decades has suggested that "string-like" cooperatively rearranging clusters observed in molecular simulations may be these CRRs-a scenario that would have profound implications for the understanding of the glass transition. The central element of this postulate is the report of an apparent zero-parameter relationship between the mass of string-like CRRs and the relaxation time. Here, we show, based on molecular dynamics simulations of multiple glass-forming liquids, that this finding is the result of an implicit adjustable parameter-a "replacement distance". This parameter is equivalent to an adjustable exponent within a generalized Adam-Gibbs relation, such that it tunes the entire functional form of the relation. Moreover, we are unable to find any objective criterion, based on the radial distribution function or the cluster fractal dimension, for selecting this replacement distance across multiple systems. We conclude that the present data do not establish that string-like cooperative rearrangements, as presently defined, are predictive of segmental relaxation via an Adam-Gibbs-like physical model.

Coupling between structural relaxation and diffusion in glass-forming liquids under pressure variation

We theoretically investigate structural relaxation and activated diffusion of glass-forming liquids at different pressures using both Elastically Collective Nonlinear Langevin Equation (ECNLE) theory and molecular dynamics (MD) simulations.

Theoretical Model for the Structural Relaxation Time in Coamorphous Drugs

We propose a simple approach to investigate the structural relaxation time and glass transition of amorphous drugs. Amorphous materials are modeled as a set of equal sized hard spheres. The structural relaxation time over many decades in hard-sphere fluids is theoretically calculated using the elastically collective nonlinear Langevin equation theory associated with Kramer's theory. Then, new thermal mapping from a real material to an effective hard-sphere fluid provides temperature-dependent relaxation time, which can be compared to experiments. Numerical results quantitatively agree with previous experiments for pharmaceutical binary mixtures having different weight ratios. We carry out experiments to test our calculations for an ezetimibe-simvastatin-Kollidon VA64 mixture. Our approach would provide a simple but comprehensive description of glassy dynamics in amorphous composites.

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

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

Effects of cooling rate on structural relaxation in amorphous drugs: elastically collective nonlinear langevin equation theory and machine learning study

Theoretical approaches are formulated to investigate the molecular mobility under various cooling rates of amorphous drugs. We describe the structural relaxation of a tagged molecule as a coupled process of cage-scale dynamics and collective molecular rearrangement beyond the first coordination shell. The coupling between local and non-local dynamics behaves distinctly in different substances. Theoretical calculations for the structural relaxation time, glass transition temperature, and dynamic fragility are carried out over twenty-two amorphous drugs and polymers. Numerical results have a quantitatively good accordance with experimental data and the extracted physical quantities using the Vogel–Fulcher–Tammann fit function and machine learning. The machine learning method reveals the linear relation between the glass transition temperature and the melting point, which is a key factor for pharmaceutical solubility. Our predictive approaches are reliable tools for developing drug formulations.

Theoretical approaches are formulated to investigate the molecular mobility under various cooling rates of amorphous drugs.