Elementary excitations and crossover phenomenon in liquids

Emergence of Debye Scaling in the Density of States of Liquids under Nanoconfinement

In the realm of nanoscience, the dynamic behaviors of liquids at scales beyond the conventional structural relaxation time, τ, unfold a fascinating blend of solid-like characteristics, including the propagation of collective shear waves and the emergence of elasticity. However, in classical bulk liquids, where τ is typically of the order of 1 ps or less, this solid-like behavior remains elusive in the low-frequency region of the density of states (DOS). Here, we provide evidence for the emergent solid-like nature of liquids at short distances through inelastic neutron scattering measurements of the low-frequency DOS in liquid water and glycerol confined within graphene oxide membranes. In particular, upon increasing the strength of confinement, we observe a transition from a liquid-like DOS (linear in the frequency ω) to a solid-like behavior (Debye law, ∼ω2) in the range of 1-4 meV. Molecular dynamics simulations confirm these findings and reveal additional solid-like features, including propagating collective shear waves and a reduction in the self-diffusion constant. Finally, we show that the onset of solid-like dynamics is pushed toward low frequency along with the slowing-down of the relaxation processes upon confinement. This nanoconfinement-induced transition, aligning with k-gap theory, underscores the potential of leveraging liquid nanoconfinement in advancing nanoscale science and technology, building more connections between fluid dynamics and materials engineering.

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl2 from Tm + 85 K to Tm + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl2. We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. This result may contribute to new understandings of the structural relaxation of other network liquids.

Understanding dynamics in coarse-grained models. IV. Connection of fine-grained and coarse-grained dynamics with the Stokes-Einstein and Stokes-Einstein-Debye relations

Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Furthermore, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the "entropy-free" CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level.

Do Specific Ion Effects on Collective Relaxation Arise from Perturbation of Hydrogen-Bonding Network Structure?

The change in the transport properties (i.e., water diffusivity, shear viscosity, etc.) when adding salts to water has been used to classify ions as either being chaotropic or kosmotropic, a terminology based on the presumption that this phenomenon arises from respective breakdown or enhancement of the hydrogen-bonding network structure. Recent quasi-elastic neutron scattering measurements of the collective structural relaxation time, τC, in aqueous salt solutions were interpreted as confirming this proposed origin of ion effects on the dynamics of water. However, we find similar changes in τC in the same salt solutions based on molecular dynamics (MD) simulations using a coarse-grained water model in which no hydrogen bonding exists, challenging this conventional interpretation of mobility change resulting from the addition of salts to water. A thorough understanding of specific ion effects should be useful in diverse material manufacturing and biomedical applications, where these effects are prevalent, but poorly understood.

Normal-to-Supercooled Liquid Transition in Molecular Glass-Formers: A Hidden Structural Transformation Fuelled by Conformational Interconversion

Molecular dynamics and transport coefficients change significantly around the so-called Arrhenius crossover in glass-forming systems. In this article, we revisit the dynamic processes occurring in a glass-forming macrocyclic crown thiaether MeBzS2O above its glass transition, revealing two crossover temperatures: TB at 309 and TA at 333 K. We identify the second one as the Arrhenius crossover that is closely related to the normal-to-supercooled liquid transition in this compound. We show that the transformation occurring at this point goes far beyond molecular dynamics (where the temperature dependence of structural relaxation times changes its character from activation-like to super-Arrhenius), being reflected also in the internal structure and diffraction pattern. In this respect, we found a twofold local organization of the nearest-neighbor molecules via weak van der Waals forces, without the formation of any medium-range order or mesophases. The nearest surrounding of each molecule evolves structurally in time due to the ongoing fast conformational changes. We identify several conformers of MeBzS2O, demonstrating that its lowest-energy conformation is preferred mainly at lower temperatures, i.e., in the supercooled liquid state. Its increased prevalence modifies locally the short-range intermolecular order and promotes vitrification. Consequently, we indicate that the Arrhenius transition is fuelled rather by conformational changes in this glass-forming macrocyclic crown thiaether, which is a different scenario from the so-far existing concepts. Our studies combine broadband dielectric spectroscopy (BDS), X-ray diffraction, Fourier transform infrared (FTIR) spectroscopy, molecular dynamics (MD) simulations, and density functional theory (DFT) calculations.

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.

Universal scaling of transverse sound speed and its isomorphic property in Yukawa fluids

Molecular dynamical simulations are performed to investigate the scaling of the transverse sound speed in two-dimensional (2D) and 3D Yukawa fluids. From the calculated diagnostics of the radial distribution function, the mean-squared displacement, and the Pearson correlation coefficient, the approximate isomorphic curves for 2D and 3D liquidlike Yukawa systems are obtained. It is found that the structure and dynamics of 2D and 3D liquidlike Yukawa systems exhibit the isomorphic property under the conditions of the same relative coupling parameter Γ/Γ_{m}=const. It is demonstrated that the reduced transverse sound speed, i.e., the ratio of the transverse sound speed to the thermal speed, is an isomorph invariant, which is a quasiuniversal function of Γ/Γ_{m}. The obtained isomorph invariant of the reduced transverse sound speed can be useful to estimate the transverse sound speed, or determine the coupling strength, with applications to dusty (complex) plasma or colloidal systems.

Ripples in the bottom of the potential energy landscape of metallic glass

In the absence of periodicity, the structure of glass is ill-defined, and a large number of structural states are found at similar energy levels. However, little is known about how these states are connected to each other in the potential energy landscape. We simulate mechanical relaxation by molecular dynamics for a prototypical [Formula: see text] metallic glass and follow the mechanical energy loss of each atom to track the change in the state. We find that the energy barriers separating these states are remarkably low, only of the order of 1 meV, implying that even quantum fluctuations can overcome these potential energy barriers. Our observation of numerous small ripples in the bottom of the potential energy landscape puts many assumptions regarding the thermodynamic states of metallic glasses into question and suggests that metallic glasses are not totally frozen at the local atomic level.

Mesoscopic two-point collective dynamics of glass-forming liquids

The collective density-density and hydrostatic pressure-pressure correlations of glass-forming liquids are spatiotemporally mapped out using molecular dynamics simulations. It is shown that the sharp rise of structural relaxation time below the Arrhenius temperature coincides with the emergence of slow, nonhydrodynamic collective dynamics on mesoscopic scales. The observed long-range, nonhydrodynamic mode is independent of wave numbers and closely coupled to the local structural dynamics. Below the Arrhenius temperature, it dominates the slow collective dynamics on length scales immediately beyond the first structural peak in contrast to the well-known behavior at high temperatures. These results highlight a key connection between the qualitative change in mesoscopic two-point collective dynamics and the dynamic crossover phenomenon.

Fast dynamics and high effective dimensionality of liquid fluidity

Fluidity, the ability of liquids to flow, is the key property distinguishing liquids from solids. This fluidity is set by the mobile transit atoms moving from one quasi-equilibrium point to the next. The nature of this transit motion is unknown. Here, we show that flow-enabling transits form a dynamically distinct sub-ensemble where atoms move on average faster than the overall system, with a manifestly non-Maxwellian velocity distribution. This is in contrast to solids and gases where no distinction of different ensembles can be made and where the distribution is always Maxwellian. The non-Maxwellian distribution is described by an exponent [Formula: see text] corresponding to high dimensionality of space. This is generally similar to extra synthetic dimensions in topological quantum matter, albeit higher dimensionality in liquids is not integer but is fractional. The dimensionality is close to 4 at melting and exceeds 4 at high temperature. [Formula: see text] has a maximum as a function of temperature and pressure in liquid and supercritical states, returning to its Maxwell value in the solid and gas states.

Electromagnetic levitation containerless processing of metallic materials in microgravity: rapid solidification

Space levitation processing allows researchers to conduct benchmark tests in an effort to understand the physical phenomena involved in rapid solidification processing, including alloy thermodynamics, nucleation and growth, heat and mass transfer, solid/liquid interface dynamics, macro- and microstructural evolution, and defect formation. Supported by ground-based investigations, a major thrust is to develop and refine robust computational tools based on theoretical and applied approaches. This work is accomplished in conjunction with experiments designed for precise model validation with application to a broad range of industrial processes.

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.

Arrhenius Crossover Temperature of Glass-Forming Liquids Predicted by an Artificial Neural Network

The Arrhenius crossover temperature, TA, corresponds to a thermodynamic state wherein the atomistic dynamics of a liquid becomes heterogeneous and cooperative; and the activation barrier of diffusion dynamics becomes temperature-dependent at temperatures below TA. The theoretical estimation of this temperature is difficult for some types of materials, especially silicates and borates. In these materials, self-diffusion as a function of the temperature T is reproduced by the Arrhenius law, where the activation barrier practically independent on the temperature T. The purpose of the present work was to establish the relationship between the Arrhenius crossover temperature TA and the physical properties of liquids directly related to their glass-forming ability. Using a machine learning model, the crossover temperature TA was calculated for silicates, borates, organic compounds and metal melts of various compositions. The empirical values of the glass transition temperature Tg, the melting temperature Tm, the ratio of these temperatures Tg/Tm and the fragility index m were applied as input parameters. It has been established that the temperatures Tg and Tm are significant parameters, whereas their ratio Tg/Tm and the fragility index m do not correlate much with the temperature TA. An important result of the present work is the analytical equation relating the temperatures Tg, Tm and TA, and that, from the algebraic point of view, is the equation for a second-order curved surface. It was shown that this equation allows one to correctly estimate the temperature TA for a large class of materials, regardless of their compositions and glass-forming abilities.

Modeling the viscoelastic relaxation dynamics of soft particles via molecular dynamics simulation-informed multi-dimensional transition-state theory

Viscoelastic soft colloidal particles have been widely explored in mechanical, chemical, pharmaceutical and other engineering applications due to their unique combination of viscosity and elasticity. The characteristic viscoelastic relaxation time shows an Arrhenius-type (or super-Arrhenius due to temperature-dependent transition attempts) thermally-activated behavior, but a holistic explanation from the relevant transition-state theory remains elusive. In this paper, the viscoelastic relaxation times of Lennard-Jones soft colloidal particle systems, including a single particle type system and a binary particle mixture based on the Kob-Andersen model, are determined using molecular dynamics (MD) simulations as the benchmark. First, the particle systems show a non-Maxwellian behavior after comparing the MD-predicted viscoelastic relaxation time and dynamic moduli (storage and loss modulus) to the classic Maxwell viscoelastic model and the recent particle local connectivity theory. Surprisingly, neither the Maxwell relaxation time τ_Maxwell (obtained from the static shear viscosity η and the high-frequency shear modulus G_∞) nor the particle local connectivity lifetime τ_LC can capture the super-Arrhenius temperature-dependent behavior in the MD-predicted relaxation time τ_MD. Then, the particle dissociation and association transition kinetics, fractal dimensions of the particle systems, and neighbor particle structure (obtained from the radial distribution functions) are shown to collectively determine the viscoelastic relaxation time. These factors are embedded into a new multi-dimensional transition kinetics model to directly estimate the viscoelastic relaxation time τ_Model, which is found to agree with the MD-predicted τ_MD remarkably well. This work highlights the microscopic origin of viscoelastic relaxation dynamics of soft colloidal particles, and theoretically connects rheological dynamics and transition kinetics in soft matters.

Structural evolution in Au- and Pd-based metallic glass forming liquids and the case for improved molecular dynamics force fields

The results of a combined experimental and computational investigation of the structural evolution of Au81Si19, Pd82Si18, and Pd77Cu6Si17 metallic glass forming liquids are presented. Electrostatically levitated metallic liquids are prepared, and synchrotron x-ray scattering studies are combined with embedded atom method molecular dynamics simulations to probe the distribution of relevant structural units. Metal–metalloid based metallic glass forming systems are an extremely important class of materials with varied glass forming ability and mechanical processibility. High quality experimental x-ray scattering data are in poor agreement with the data from the molecular dynamics simulations, demonstrating the need for improved interatomic potentials. The first peak in the x-ray static structure factor in Pd77Cu6Si17 displays evidence for a Curie–Weiss type behavior but also a peak in the effective Curie temperature. A proposed order parameter distinguishing glass forming ability, [Formula: see text], shows a peak in the effective Curie temperature near a crossover temperature established by the behavior of the viscosity, T A.

Quantum mechanical interpretation of the minimum viscosity of metallic liquids

Possible fundamental quantum bounds for viscosity and many other physical properties have drawn serious considerations recently from diverse communities encompassing those studying quantum gravity, high-energy physics, condensed matter physics, strongly correlated electron systems, and "strange metals," to name a few. However, little attention has been paid by materials scientists and the fluid dynamics community, perhaps because of the general belief that quantum mechanics is of little consequence for classical fluid dynamics. Here, considering the extrapolated high-temperature viscosity of 32 metallic alloy liquids as representative of minimum viscosity, experimental results are presented and evaluated in terms of a number of quantum- and statistical-mechanics-based theories. The surprising result is that the experimental data are within one order of magnitude of estimates from those theories. That quantum mechanics could be of importance at such high temperatures in conventional classical fluids is quite interesting. Another surprise is that the minimum viscosities of metallic alloy liquids are not too different from an archetypal quantum liquid, such as He.

Liquid-like atoms in dense-packed solid glasses

Revealing the microscopic structural and dynamic pictures of glasses is a long-standing challenge for scientists^1,2. Extensive studies on the structure and relaxation dynamics of glasses have constructed the current classical picture^3-5: glasses consist of some 'soft zones' of loosely bound atoms embedded in a tightly bound atomic matrix. Recent experiments have found an additional fast process in the relaxation spectra^6-9, but the underlying physics of this process remains unclear. Here, combining extensive dynamic experiments and computer simulations, we reveal that this fast relaxation is associated with string-like diffusion of liquid-like atoms, which are inherited from the high-temperature liquids. Even at room temperature, some atoms in dense-packed metallic glasses can diffuse just as easily as they would in liquid states, with an experimentally determined viscosity as low as 10⁷ Pa·s. This finding extends our current microscopic picture of glass solids and might help establish the dynamics-property relationship of glasses⁴.

Quantum chemical investigation of the effect of alkali metal ions on the dynamic structure of water in aqueous solutions

We report quantum chemical molecular dynamics (MD) simulations based on the density-functional tight-binding (DFTB) method to investigate the effect of K+, Na+, and Mg2+ ions in aqueous solutions on the static and dynamic structure of bulk water at room temperature and with various concentrations. The DFTB/MD simulations were validated for the description of ion solvation in aqueous ionic solutions by comparing static pair distribution functions (PDFs) as well as the cation solvation shell between experimental and available ab initio DFT data. The effect of the cations on the water structure, as well as relative differences between K+, Na+, and Mg2+ cations, were analyzed in terms of atomically resolved PDFs as well as time-dependent Van Hove correlation functions (VHFs). The investigation of the VHFs reveals that salt ions generally slow down the dynamic decay of the pair correlations in the water solvation sphere, irrespective of the cation size or charge. The analysis of partial metal-oxygen VHFs indicates that there are long-lived correlations between water and Na+ over long distances, in contrast to K+ and Mg2+.

Double universality of the transition in the supercritical state

Universality aids consistent understanding of physical properties and states of matter where a theory predicts how a property of a phase (solid, liquid, and gas) changes with temperature or pressure. Here, we show that the matter above the critical point has a remarkable double universality. The first universality is the transition between the liquid-like and gas-like states seen in the crossover of the specific heat on the dynamical length with a fixed inversion point. The second universality is the operation of this effect in many supercritical fluids, including N₂, CO₂, Pb, H₂O, and Ar. Despite different structure and chemical bonding, the transition has the same fixed inversion point deep in the supercritical state. This advances our understanding of the supercritical state previously considered to be a featureless area on the phase diagram and a theoretical guide for improved deployment of supercritical fluids in green and environmental applications.

Real-Space Local Dynamics of Molten Inorganic Salts Using Van Hove Correlation Function

Molten inorganic salts are attracting resurgent attention because of their unique physicochemical properties, making them promising media for next-generation concentrating solar power systems and molten salt reactors. The dynamics of these highly disordered ionic media is largely studied by theoretical simulations, while the robust experimental techniques capable of observing local dynamics are not well-developed. To provide fundamental insights into the atomic-scale transport properties of molten salts, we report the real-space dynamics of molten magnesium chloride at high temperatures employing the Van Hove correlation function obtained by inelastic neutron scattering. Our results directly depict the distance-dependent dynamics of a molten salt on the picosecond time scale. This study demonstrates the capability of the developed approach to describe the locally correlated- and self-dynamics in molten salts, significantly improving our understanding of the interplay between microscopic structural parameters and their dynamics that ultimately control physical properties of condensed matter in extreme environments.

Nonaffine Strains Control Ductility of Metallic Glasses

The origin of limited plasticity in metallic glasses is elusive, with no apparent link to their atomic structure. We propose that the response of the glassy structure to applied stress, not the original structure itself, provides a gauge to predict the degree of plasticity. We carried out high-energy x-ray diffraction on various bulk metallic glasses (BMGs) under uniaxial compression within the elastic limit and evaluated the anisotropic pair distribution function. We show that the extent of local deviation from the affine (uniform) deformation in the elastic regime is strongly correlated with the plastic behavior of BMGs beyond yield, across chemical compositions and sample history. The results suggest that the propensity for collective local atomic rearrangements under stress promotes plasticity.

Medium-range atomic correlation in simple liquids. II. Theory of temperature dependence

The spatial atomic correlations in liquids and glasses extend often significantly beyond the nearest neighbors. Such correlations, called the medium-range order (MRO), affect many physical properties, but their nature is not well understood. In this article the variation of the MRO with temperature is calculated based upon the concept of the atomic-level pressure, focusing on simple liquids, such as metallic liquids. It is shown that the structural coherence length that characterizes MRO follows the Curie-Weiss law with a negative Curie temperature as observed by experiment and simulation. It is also shown that the glass transition is induced by freezing of the MRO, rather than the freezing of the nearest-neighbor shell. The implications of these results are discussed.

Breakdown of the Stokes-Einstein relationship and rapid structural ordering in CuZrAl metallic glass-forming liquids

The results of a combined structural and dynamical study of Cu-Zr-Al metallic glass forming liquids are presented. Containerless high-energy x-ray scattering experiments made using electrostatic levitation are combined with molecular dynamics simulations to probe the onset of rapid structural ordering as well as the temperature-dependent diffusivity and viscosity in three liquids: Cu₄₉Zr₄₅Al₆, Cu₄₇Zr₄₅Al₈, and Cu₄₃Zr₄₅Al₁₂. These compositions were chosen because they are reported to have dramatically different glass forming-ability. Experimental data show that the first peak in the x-ray static structure factor displays evidence for a Curie-Weiss type behavior, but also a peak in the effective Curie temperature. The evidence provided here for the onset of cooperativity, marked by a crossover temperature, T_A (which is usually above the liquidus temperature), is accompanied by the onset of development of more spatially extended structural order in the liquids. Based on the molecular dynamics simulations, each of the liquids exhibits a clear breakdown of the Stokes-Einstein relation at a temperature near, but below, the crossover temperature, T_A. The breakdown is manifest as a rapid reduction in the relative diffusion coefficients between Cu, Zr, and Al.

Universal interrelation between dynamics and thermodynamics and a dynamically driven "c" transition in fluids

Our very wide survey of the supercritical phase diagram and its key properties reveals a universal interrelation between dynamics and thermodynamics and an unambiguous transition between liquidlike and gaslike states. This is seen in the master plot showing a collapse of the data representing the dependence of specific heat on key dynamical parameters in the system for many different paths on the phase diagram. As a result, the observed transition is path independent. We call it a "c" transition due to the c-shaped curve parametrizing the dependence of the specific heat on key dynamical parameters. The c transition has a fixed inversion point and provides a new structure to the phase diagram, operating deep in the supercritical state (up to, at least, 2000 times the critical pressure and 50 times the critical temperature). The data collapse and path independence as well as the existence of a special inversion point on the phase diagram are indicative of either of a sharp crossover or a new phase transition in the deeply supercritical state.

Microscopic Structural Evolution during Ultrastable Metallic Glass Formation

By decreasing the rate of physical vapor deposition, ZrCuAl metallic glasses with improved stability and mechanical performances can be formed, while the microscopic structural mechanisms remain unclear. Here, with scanning transmission electron microscopy and high-energy synchrotron X-ray diffraction, we found that the metallic glass deposited at a higher rate exhibits a heterogeneous structure with compositional fluctuations at a distance of a few nanometers, which gradually disappear on decreasing the deposition rate; eventually, a homogeneous structure is developed approaching ultrastability. This microscopic structural evolution suggests the existence of the following two dynamical processes during ultrastable metallic glass formation: a faster diffusion process driven by the kinetic energy of the depositing atoms, which results in nanoscale compositional fluctuations, and a slower collective relaxation process that eliminates the compositional and structural heterogeneity, equilibrates the deposited atoms, and strengthens the local atomic connectivity.

The intimate relationship between structural relaxation and the energy landscape of monatomic liquid metals

The characteristic property of a liquid, discriminating it from a solid, is its fluidity, which can be expressed by a velocity field. The reaction of the velocity field on forces is enshrined in the transport parameter viscosity. In contrast, a solid reacts to forces elastically through a displacement field, the particles are trapped in their potential minimum. The flow in a liquid needs enough thermal energy to overcome the changing potential barriers, which is supported through a continuous rearrangement of surrounding particles. Cooling a liquid will decrease the fluidity of a particle and the mobility of the neighbouring particles, resulting in an increase of the viscosity until the system comes to an arrest. This process with a concomitant slowing down of collective particle rearrangements might already start deep inside the liquid state. The idea of the potential energy landscape provides an attractive picture for these dramatic changes. However, despite the appealing idea there is a scarcity of quantitative assessments, in particular, when it comes to experimental studies. Here we present results on a monatomic liquid metal through a combination of ab initio molecular dynamics, neutron spectroscopy and inelastic x-ray scattering. We investigated the collective dynamics of liquid aluminium to reveal the changes in dynamics when the high temperature liquid is cooled towards solidification. The results demonstrate the main signatures of the energy landscape picture, a reduction in the internal atomic structural energy, a transition to a stretched relaxation process and a deviation from the high-temperature Arrhenius behavior of the relaxation time. All changes occur in the same temperature range at about [Formula: see text], which can be regarded as the temperature when the liquid aluminium enters the landscape influenced phase and enters a more viscous liquid state towards solidification. The similarity in dynamics with other monatomic liquid metals suggests a universal dynamic crossover above the melting point.

Building and Breaking Bonds by Homogenous Nucleation in Glass-Forming Melts Leading to Transitions in Three Liquid States

The thermal history of melts leads to three liquid states above the melting temperatures T_m containing clusters-bound colloids with two opposite values of enthalpy +Δε_lg × ΔH_m and -Δε_lg × ΔH_m and zero. All colloid bonds disconnect at T_n+ > T_m and give rise in congruent materials, through a first-order transition at T_LL = T_n+, forming a homogeneous liquid, containing tiny superatoms, built by short-range order. In non-congruent materials, (T_n+) and (T_LL) are separated, T_n+ being the temperature of a second order and T_LL the temperature of a first-order phase transition. (T_n+) and (T_LL) are predicted from the knowledge of solidus and liquidus temperatures using non-classical homogenous nucleation. The first-order transition at T_LL gives rise by cooling to a new liquid state containing colloids. Each colloid is a superatom, melted by homogeneous disintegration of nuclei instead of surface melting, and with a Gibbs free energy equal to that of a liquid droplet containing the same magic atom number. Internal and external bond number of colloids increases at T_n+ or from T_n+ to T_g. These liquid enthalpies reveal the natural presence of colloid-colloid bonding and antibonding in glass-forming melts. The Mpemba effect and its inverse exist in all melts and is due to the presence of these three liquid states.

Relevance of hydrogen bonded associates to the transport properties and nanoscale dynamics of liquid and supercooled 2-propanol

2-Propanol was investigated, in both the liquid and supercooled states, as a model system to study how hydrogen bonds affect the structural relaxation and the dynamics of mesoscale structures, of approximately several Ångstroms, employing static and quasi-elastic neutron scattering and molecular dynamics simulation. Dynamic neutron scattering measurements were performed over an exchanged wave-vector range encompassing the pre-peak, indicative of the presence of H-bonding associates, and the main peak. The dynamics observed at the pre-peak is associated with the formation and disaggregation of the H-bonded associates and is measured to be at least one order of magnitude slower than the dynamics at the main peak, which is identified as the structural relaxation. The measurements indicate that the macroscopic shear viscosity has a similar temperature dependence as the dynamics of the H-bonded associates, which highlights the important role played by these structures, together with the structural relaxation, in defining the macroscopic rheological properties of the system. Importantly, the characteristic relaxation time at the pre-peak follows an Arrhenius temperature dependence whereas at the main peak it exhibits a non-Arrhenius behavior on approaching the supercooled state. The origin of this differing behavior is attributed to an increased structuring of the hydrophobic domains of 2-propanol accommodating a more and more encompassing H-bond network, and a consequent set in of dynamic cooperativity.

Observation of General Entropy-Enthalpy Compensation Effect in the Relaxation of Wrinkled Polymer Nanocomposite Films

Surface-textured polymer nanocomposite (PNC) films are utilized in many device applications, and therefore understanding the relaxation behavior of such films is important. By extending an in situ wrinkle relaxation method, we observed that the thermal stability of wrinkled PNC films, both above and below the glass transition temperature (T_g), is proportional to a film's nanoparticle (polymer grafted and bare) concentration, with a slope that changes sign at a compensation temperature (T_comp) that is determined to be in the vicinity of the film's T_g. This provides unambiguous confirmation of entropy-enthalpy compensation (EEC) as a general feature of PNC films, implying that the stability of PNC films changes from being enhanced to becoming diminished by simply passing through this characteristic temperature, a phenomenon having evident practical ramifications. We suggest EEC will also arise in films where residual stresses are associated with the film fabrication process, which is relevant to nanotech device applications.

Crossover of dynamical instability and chaos in the supercritical state

We calculate the maximal Lyapunov exponent for a bulk system of 256 Lennard-Jones particles in constant energy molecular dynamics simulations deep into the supercritical state. We find that the maximal Lyapunov exponent undergoes a crossover and that this crossover coincides with the dynamical crossover at the Frenkel line from liquid physics. We explain this crossover in terms of two different contributions to dynamical instability: diffusion in the liquidlike state below the Frenkel line and collisions in the gaslike state above. These results provide insight into the phase-space dynamics far from the melting line and densities where rare-gas approximation are inapplicable.

Origin of liquid fragility

Liquid fragility characterizes how steeply the viscosity of a glass-forming liquid decreases with increasing temperature above the glass transition. It is one of the most fundamental properties of a liquid, with high importance for science and application. Yet, its origin is unclear. Here we show that it is directly related to the structural coherence of the medium-range order (MRO) in liquid defined by the decay of the pair-distribution function with distance. The MRO can also be evaluated from the first peak of the structure function determined by x-ray or neutron diffraction, and it is a measure of the cooperativity of atomic motion in a diffusive event in supercooled liquids. These findings shed light on the mechanism of atomic transport in supercooled liquids.

Correlated atomic dynamics in liquid seen in real space and time

In liquids, the timescales for structure, diffusion, and phonon are all similar, of the order of a pico-second. This not only makes characterization of liquid dynamics difficult but also renders it highly questionable to describe liquids in these terms. In particular, the current definition of the structure of liquids by the instantaneous structure may need to be expanded because the liquid structure is inherently dynamic. Here, we advocate describing the liquid structure through the distinct-part of the Van Hove function, which can be determined by inelastic neutron and x-ray scattering measurements as well as by simulation. It depicts the dynamic correlation between atoms in space and time, starting with the instantaneous correlation function at t = 0. The observed Van Hove functions show that the atomic dynamics is strongly correlated in some liquids, such as water. The effect of atomic correlation on various transport properties of fluid, including viscosity and diffusivity, is discussed.

Influence of Kosmotrope and Chaotrope Salts on Water Structural Relaxation

The structural relaxation in water solutions of kosmotrope (structure maker) and chaotrope (structure breaker) salts, namely sodium chloride, potassium chloride, and cesium chloride, were studied through quasielastic neutron scattering measurements. We found that the collective dynamics relaxation time at the structure factor peak obtained using heavy water solutions shows a distinctively different behavior in the kosmotrope as opposed to the chaotrope solutions, increasing with the salt concentration in the former and decreasing in the latter. In both cases the trends are proportional to the concentration dependence of the relative viscosity of the solutions. These results indicate that kosmotropes and chaotropes influence the solution's viscosity by impacting in opposite ways the hydrogen bond network of water, strengthening it in one case and softening it in the other.

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.

Local self-motion of water through the Van Hove function

We show that the self-part of the Van Hove function-the correlation function describing the dynamics of a single molecule-of water can be determined through a high-resolution inelastic x-ray scattering experiment. The measurement of inelastic x-ray scattering up to 10Å^{-1} makes it possible to convert the inelastic x-ray scattering spectra into the Van Hove function, and its self-part is extracted from the short-range correlations. The diffusivity estimated from the short-range dynamics of water molecules is different from the long-range diffusivity measured by other methods. This approach using the experimentally determined self-part of the Van Hove function will be useful to the study of the local dynamics of atoms and molecules in liquids.

On Viscous Flow in Glass-Forming Organic Liquids

The two-exponential Sheffield equation of viscosity η(T) = A1·T·[1 + A2·exp(Hm/RT)]·[1 + C·exp(Hd/RT)], where A1, A2, Hm, C, and Hm are material-specific constants, is used to analyze the viscous flows of two glass-forming organic materials-salol and α-phenyl-o-cresol. It is demonstrated that the viscosity equation can be simplified to a four-parameter version: η(T) = A·T·exp(Hm/RT)]·[1 + C·exp(Hd/RT)]. The Sheffield model gives a correct description of viscosity, with two exact Arrhenius-type asymptotes below and above the glass transition temperature, whereas near the Tg it gives practically the same results as well-known and widely used viscosity equations. It is revealed that the constants of the Sheffield equation are not universal for all temperature ranges and may need to be updated for very high temperatures, where changes occur in melt properties leading to modifications of A and Hm for both salol and α-phenyl-o-cresol.

Ideality of liquid structure: A case study for metallic alloy liquids

It is difficult to characterize by experiment the structural features of liquids and glasses which lack long-range translational periodicity in the structure. Here, we suggest that the height and shape of the first peak of the structure function S(Q) carry significant information about the nature of the medium-range order and the coherence of density correlations. It is further proposed that they indicate how ideal the liquid structure is. Here, the ideal state is defined by long-range density correlations, not by structural coherence at the atomic level. The analysis is applied to the S(Q) of metallic alloy liquids determined by x-ray diffraction and simulation. The ideality index defined here may provide a common parameter to characterize structural coherence among various disparate groups of liquids and glasses.

Experimental determination of the temperature-dependent Van Hove function in a Zr80Pt20 liquid

Even though the viscosity is one of the most fundamental properties of liquids, the connection with the atomic structure of the liquid has proven elusive. By combining inelastic neutron scattering with the electrostatic levitation technique, the time-dependent pair-distribution function (i.e., the Van Hove function) has been determined for liquid Zr₈₀Pt₂₀. We show that the decay time of the first peak of the Van Hove function is directly related to the Maxwell relaxation time of the liquid, which is proportional to the shear viscosity. This result demonstrates that the local dynamics for increasing or decreasing the coordination number of local clusters by one determines the viscosity at high temperature, supporting earlier predictions from molecular dynamics simulations.

Curie-Weiss behavior of liquid structure and ideal glass state

We present the results of a structural study of metallic alloy liquids from high temperature through the glass transition. We use high energy X-ray scattering and electro-static levitation in combination with molecular dynamics simulation and show that the height of the first peak of the structure function, S(Q) - 1, follows the Curie-Weiss law. The structural coherence length is proportional to the height of the first peak, and we suggest that its increase with cooling may be related to the rapid increase in viscosity. The Curie temperature is negative, implying an analogy with spin-glass. The Curie-Weiss behavior provides a pathway to an ideal glass state, a state with long-range correlation without lattice periodicity, which is characterized by highly diverse local structures, reminiscent of spin-glass.

Resistivity Saturation in Metallic Liquids Above a Dynamical Crossover Temperature Observed in Measurements Aboard the International Space Station

Although a resistivity saturation (minimum conductivity) is often observed in disordered metallic solids, such phenomena in the corresponding liquids are not known. Here we report a saturation of the electrical resistivity in Zr_{64}Ni_{36} and Cu_{50}Zr_{50} liquids above a dynamical crossover temperature for the viscosity (T_{A}). The measurements were made for the levitated liquids under the microgravity conditions of the International Space Station. Based on recent molecular dynamics simulations, the saturation is likely due to the ineffectiveness of electron-phonon scattering above T_{A} when the phonon lifetime becomes too short compared to the electron relaxation time. This is different from the conventional resistivity saturation mechanisms in solids.

Identifying Water-Anion Correlated Motion in Aqueous Solutions through Van Hove Functions

Electrolyte solutions are ubiquitous in materials in daily use and in biological systems. However, the understanding of their molecular and ionic dynamics, particularly those of their correlated motions, are elusive despite extensive experimental, theoretical, and numerical studies. Here we report the real-space observations of the molecular/ionic-correlated dynamics of aqueous salt (NaCl, NaBr, and NaI) solutions using the Van Hove functions obtained by high-resolution inelastic X-ray scattering measurement and molecular dynamics simulation. Our results directly depict the distance-dependent dynamics of aqueous salt solutions on the picosecond time scale and identify the changes in the anion-water correlations. This study demonstrates the capability of the real-space Van Hove function analysis to describe the local correlated dynamics in aqueous salt solutions.

Heat capacity of simple liquids in light of hydrodynamics as U(1) gauge theory

We investigate the heat capacity of simple liquids through a theoretical approach based on a quasiparticle description. By interpreting the microscopic dynamics of particles in liquids in terms of quasiparticles, we suggest a simplified understanding of the number of degrees of freedom in liquids. A equivalence between hydrodynamics and U(1) gauge theory, which is proposed in the present paper, develops the quasiparticle description to construct a new Lagrangian which correctly reproduces the number of modes at the melting points and at the critical points. The heat capacity evaluated from this Lagrangian naturally interpolates between these two points, and agrees with the phonon theory of liquids [Sci. Rep. 2, 421 (2012)2045-232210.1038/srep00421].

Temperature Dependences of Peak Positions in Pair Distribution Function of Metallic Liquids

The temperature dependences of the peak positions in pair distribution functions G(r) of pure metallic zinc (Zn) and indium (In) liquids have been studied using high-energy X-ray diffraction together with ab initio molecular dynamic simulations. It has been demonstrated that the first peak positions in G(r) of both Zn and In move to small r, whereas the second peak positions exhibit opposite movements with increasing temperature, originating from different thermal responses of polyhedron connections. However, the third, above peaks in G(r) in both liquids shift to large r with the expansion coefficients smaller than the values of bulk liquids.

Anticrossing of Longitudinal and Transverse Modes in Simple Fluids

If interacting modes of the same symmetry cross, they repel from each other and become hybridized. This phenomenon is called anticrossing and is well-known for mechanical oscillations, electromagnetic circuits, waveguides, metamaterials, polaritons, and phonons in crystals, but it still remains poorly understood in simple fluids. Here, we show that structural disorder and anharmonicity, governing properties of fluids, lead to the anticrossing of longitudinal and transverse modes, which is accompanied by their hybridization and strong redistribution of excitation spectra. We combined theory and simulations for noble gases to prove the reliability of mode anticrossing in simple fluids, studied here for the first time. Our results open novel prospects in understanding collective dynamics, thermodynamics, and transport phenomena in various fluids, spanning from noble gas fluids and metallic melts to strongly coupled plasmas and molecular and complex fluids.

Absorption and secondary scattering of X-rays with an off-axis small beam for a cylindrical sample geometry

Expressions for X-ray absorption and secondary scattering are developed for cylindrical sample geometries. The incident-beam size is assumed to be smaller than the sample and in general directed off-axis onto the cylindrical sample. It is shown that an offset beam has a non-negligible effect on both the absorption and multiple scattering terms, resulting in an asymmetric correction that must be applied to the measured scattering intensities. The integral forms of the corrections are first presented. A small-beam limit is then developed for easier computation.

Reply to "Comment on 'Lagrangian formulation and symmetrical description of liquid dynamics' "

We make three remarks in reply to the comment by Bryk, Duviryak, and Mryglod (BDM) [Bryk, Duviryak, and Mryglod, Phys. Rev. E 99, 036102 (2019)10.1103/PhysRevE.99.036102]: (a) the discussion of shear liquid dynamics cannot be incorrect for the reason that this discussion does not include other effects such as longitudinal fluctuations; (b) the same point of relaxation time has been already discussed and published by Bryk et al. in their earlier comment [Bryk, Mryglod, Ruocco, and Scopigno, Phys. Rev. Lett. 120, 219601 (2018)PRLTAO0031-900710.1103/PhysRevLett.120.219601] and in our related reply [Yang, Dove, Brazhkin, and Trachenko, Phys. Rev. Lett. 120, 219602 (2018)PRLTAO0031-900710.1103/PhysRevLett.120.219602]; and (c) the field transformation for the complex scalar field theory used by BDM is incorrect.