Probing cooperative liquid dynamics with the mean square displacement

Superstructure and Correlated Na+ Hopping in a Layered Mg-Substituted Sodium Manganate Battery Cathode are Driven by Local Electroneutrality

In this work, we present a variable-temperature 23Na NMR and variable-temperature and variable-frequency electron paramagnetic resonance (EPR) analysis of the local structure of a layered P2 Na-ion battery cathode material, Na0.67[Mg0.28Mn0.72]O2 (NMMO). For the first time, we elucidate the superstructure in this material by using synchrotron X-ray diffraction and total neutron scattering and show that this superstructure is consistent with NMR and EPR spectra. To complement our experimental data, we carry out ab initio calculations of the quadrupolar and hyperfine 23Na NMR shifts, the Na+ ion hopping energy barriers, and the EPR g-tensors. We also describe an in-house simulation script for modeling the effects of ionic mobility on variable-temperature NMR spectra and use our simulations to interpret the experimental spectra, available upon request. We find long-zigzag-type Na ordering with two different types of Na sites, one with high mobility and the other with low mobility, and reconcile the tendency toward Na+/vacancy ordering to the preservation of local electroneutrality. The combined magnetic resonance methodology for studying local paramagnetic environments from the perspective of electron and nuclear spins will be useful for examining the local structures of materials for devices.

Configurational Entropy Effects on Glass Transition in Metallic Glasses

Configurational entropy (S_conf) is known to be a key thermodynamic factor governing a glass transition process. However, this significance remains speculative because S_conf is not directly measurable. In this work, we demonstrate the role of S_conf theoretically and experimentally by a comparative study of a Zr-Ti-Cu-Ni-Be high-entropy metallic glass (HE-MG) with one of its conventional MG counterparts. It is revealed that the higher S_conf leads to a glass that is energetically more stable and structurally more ordered. This is manifested by ab initio molecular dynamics simulations, showing that ∼60% fewer atoms are agitated above T_g, and experimental results of smaller heat capacity jump, inconspicuous stiffness loss, insignificant structural change during glass transition, and a more depressed boson peak in the HE-MG than its counterpart. We accordingly propose a model to explain that a higher S_conf promotes a faster degeneracy-dependent kinetics for exploration of the potential energy landscape upon glass transition.

Understanding the role of cross-link density in the segmental dynamics and elastic properties of cross-linked thermosets

Cross-linking is known to play a pivotal role in the relaxation dynamics and mechanical properties of thermoset polymers, which are commonly used in structural applications because of their light weight and inherently strong nature. Here, we employ a coarse-grained (CG) polymer model to systematically explore the effect of cross-link density on basic thermodynamic properties as well as corresponding changes in the segmental dynamics and elastic properties of these network materials upon approaching their glass transition temperatures (T_g). Increasing the cross-link density unsurprisingly leads to a significant slowing down of the segmental dynamics, and the fragility K of glass formation shifts in lockstep with T_g, as often found in linear polymer melts when the polymer mass is varied. As a consequence, the segmental relaxation time τ_α becomes almost a universal function of reduced temperature, (T - T_g)/T_g, a phenomenon that underlies the applicability of the "universal" Williams-Landel-Ferry (WLF) relation to many polymer materials. We also test a mathematical model of the temperature dependence of the linear elastic moduli based on a simple rigidity percolation theory and quantify the fluctuations in the local stiffness of the network material. The moduli and distribution of the local stiffness likewise exhibit a universal scaling behavior for materials having different cross-link densities but fixed (T - T_g)/T_g. Evidently, T_g dominates both τ_α and the mechanical properties of our model cross-linked polymer materials. Our work provides physical insights into how the cross-link density affects glass formation, aiding in the design of cross-linked thermosets and other structurally complex glass-forming materials.

A general structural order parameter for the amorphous solidification of a supercooled liquid

The persistent problem posed by the glass transition is to develop a general atomic level description of amorphous solidification. The answer proposed in this paper is to measure a configuration's capacity to restrain the motion of the constituent atoms. Here, we show that the instantaneous normal modes can be used to define a measure of atomic restraint that accounts for the difference between fragile and strong liquids and the collective length scale of the supercooled liquid. These results represent a significant simplification of the description of amorphous solidification and provide a powerful systematic treatment of the influence of microscopic factors on the formation of an amorphous solid.

Application of cell models to the melting and sublimation lines of the Lennard-Jones and related potential systems

Harmonic cell models (HCMs) are shown to predict the melting line of the Lennard-Jones (LJ) but not the sublimation line. In addition, even for the melting line, the HCMs are found to be physically unrealistic for inverse power potential systems near the hard-sphere limit, and for the Weeks-Chandler-Andersen system at extremely low temperatures. Despite this, the HCM accurately predicts the LJ mean-square displacement (MSD) from molecular-dynamics (MD) simulations along both lines after simple scaling corrections, to include the effects of anharmonicity and correlated dynamics of the atoms, are applied. Single caged atom molecular dynamics and Monte Carlo simulations provide further quantitative characterization of these additional effects, which go beyond harmonicity. The melting indicator and a modification of the cell model in a similar form are shown to be approximately constant along the melting line, which indicates an isomorph. The less well studied LJ sublimation line is shown not to be an isomorph, yet it still can be represented analytically very accurately by the relationship k_{B}T=aρ^{4}+bρ^{2}, where a and b are constants (k_{B} is Boltzmann's constant, T is the temperature, and ρ is the number density). This relationship has been found previously for the melting line, but the two constants have opposite signs for the sublimation and melting lines. This simple formula is also predicted using a nonharmonic static lattice expression for the pressure. The probability distribution function of the melting factor indicates departures from harmonic or Gaussian behavior in the lower wing. Nevertheless, the mean melting factor is shown to follow a simple MSD Debye-Waller factor dependence along both the melting and sublimation lines. This work combining HCM and MD simulations provides a comparison of the melting and sublimation lines of the LJ system, which could provide the foundations for a more unified statistical mechanical description of these two solid boundary lines.

Activation free energy gradient controls interfacial mobility gradient in thin polymer films

We examine the mobility gradient in the interfacial region of substrate-supported polymer films using molecular dynamics simulations and interpret these gradients within the string model of glass-formation. No large gradients in the extent of collective motion exist in these simulated films, and an analysis of the mobility gradient on a layer-by-layer basis indicates that the string model provides a quantitative description of the relaxation time gradient. Consequently, the string model indicates that the interfacial mobility gradient derives mainly from a gradient in the high-temperature activation enthalpy ΔH₀ and entropy ΔS₀ as a function of depth z, an effect that exists even in the high-temperature Arrhenius relaxation regime far above the glass transition temperature. To gain insight into the interfacial mobility gradient, we examined various material properties suggested previously to influence ΔH₀ in condensed materials, including density, potential and cohesive energy density, and a local measure of stiffness or u²(z)^-3/2, where u²(z) is the average mean squared particle displacement at a caging time (on the order of a ps). We find that changes in local stiffness best correlate with changes in ΔH₀(z) and that ΔS₀(z) also contributes significantly to the interfacial mobility gradient, so it must not be neglected.

From soft- to hard-sphere fluids: Crossover evidenced by high-frequency elastic moduli

The conventional (Zwanzig-Mountain) expressions for instantaneous elastic moduli of simple fluids predict their divergence as the limit of hard-sphere (HS) interaction is approached. However, elastic moduli of a true HS fluid are finite. Here we demonstrate that this paradox reveals the soft-to-hard-sphere crossover in fluid excitations and thermodynamics. With extensive in silico study of fluids with repulsive power-law interactions (∝r^{-n}), we locate the crossover at n≃10-20 and develop a simple and accurate model for the HS regime. The results open prospects to deal with the elasticity and related phenomena in various systems, from simple fluids to melts and glasses.

Sound Velocities of Generalized Lennard-Jones (n - 6) Fluids Near Freezing

In a recent paper [S. Khrapak, Molecules 25, 3498 (2020)], the longitudinal and transverse sound velocities of a conventional Lennard-Jones system at the liquid-solid coexistence were calculated. It was shown that the sound velocities remain almost invariant along the liquid-solid coexistence boundary lines and that their magnitudes are comparable with those of repulsive soft-sphere and hard-sphere models at the fluid-solid phase transition. This implies that attraction does not considerably affect the magnitude of the sound velocities at the fluid-solid phase transition. This paper provides further evidence to this by examining the generalized Lennard-Jones (n - 6) fluids with n ranging from 12 to 7 and demonstrating that the steepness of the repulsive term has only a minor effect on the magnitude of the sound velocities. Nevertheless, these minor trends are identified and discussed.

Localization model description of the interfacial dynamics of crystalline Cu and [Formula: see text] metallic glass nanoparticles

Many of the special properties of nanoparticles (NPs) and nanomaterials broadly derive from the significant fraction of particles (atoms, molecules or segments of polymeric molecules) in the NP interfacial region in which the interparticle interactions are characteristically highly anharmonic in comparison to the bulk material. This leads to relatively large mean square particle displacements relative to the material interior, often resulting in a strong increase interfacial mobility and reactivity in both crystalline and glass NPs. The 'Debye-Waller factor', or the mean square particle displacement [Formula: see text] on a ps 'caging' timescale relative to the square of the average interparticle distance [Formula: see text], provides an often experimentally accessible measure of the strength of this anharmonic interaction. The Localization Model (LM) of the dynamics of condensed materials relates this thermodynamic property to the structural relaxation time [Formula: see text], determined from the intermediate scattering function, without any free parameters. Moreover, the LM allows for the prediction of the diffusion coefficient D when combined with the 'decoupling' or Fractional Stokes-Einstein relation linking [Formula: see text] to D. In the current study, we employed classical molecular dynamics simulation to investigate the structural relaxation and diffusion of model [Formula: see text] metallic glass and Cu crystalline NPs with different sizes. As with previous studies validating the LM on model bulk and crystalline materials, and for the interfacial dynamics of thin crystalline and metallic glass films, we find the LM model also describes the interfacial dynamics of model crystalline metal (Cu) and metallic glass ([Formula: see text] NPs to a good approximation, further confirming the generality of the model.

Sound Velocities of Lennard-Jones Systems Near the Liquid-Solid Phase Transition

Longitudinal and transverse sound velocities of Lennard-Jones systems are calculated at the liquid-solid coexistence using the additivity principle. The results are shown to agree well with the "exact" values obtained from their relations to excess energy and pressure. Some consequences, in particular in the context of the Lindemann's melting rule and Stokes-Einstein relation between the self-diffusion and viscosity coefficients, are discussed. Comparison with available experimental data on the sound velocities of solid argon at melting conditions is provided.

An inelastic neutron scattering, Raman, far-infrared, and molecular dynamics study of the intermolecular dynamics of two ionic liquids

The intermolecular dynamics in the THz frequency range of the ionic liquids n-butyl-trimethylammonium bis(trifluoromethanesulfonyl)imide, [N1114][NTf2], and methyl-tributylammonium bis(trifluoromethanesulfonyl)imide, [N1444][NTf2], were investigated by a combined usage of inelastic neutron scattering (INS), Raman, and far-infrared (FIR) spectroscopies and the power spectrum calculated by molecular dynamics (MD) simulations. The collective dynamics of the simulated systems is also discussed by the calculation of time correlation functions of charge and mass currents that are projected onto acoustic- and optic-like motions. The INS and Raman measurements have been performed as a function of temperature in the glassy, crystalline, and liquid phases. The excess in the vibrational density of states over the expectation of the Debye theory, the so-called boson peak, is found in the INS and Raman spectra as a peak at ∼2 meV (∼16 cm-1) and also in the direct measurement of heat capacity at very low temperatures (4-20 K). This low-frequency vibration is incorporated into the curve fits of Raman, FIR, and MD data at room temperature. Fits of spectra from these different sources in the range below 100 cm-1 are consistently achieved with three components at ca. 25, 50, and 80 cm-1, but with distinct relative intensities among the different techniques. It is proposed as the collective nature of the lowest-frequency component and the anion-cation intermolecular vibration nature of the highest-frequency component. The MD results indicate that there is no clear distinction between acoustic and optic vibrations in the spectral range investigated in this work for the ionic liquids [N1114][NTf2] and [N1444][NTf2]. The analysis carried out here agrees in part, but not entirely, with other propositions in the literature, mainly from optical Kerr effect (OKE) and FIR spectroscopies, concerning the intermolecular dynamics of ionic liquids.

Elastic properties of dense hard-sphere fluids

A new analysis of elastic properties of dense hard-sphere (HS) fluids is presented, based on the expressions derived by Miller [J. Chem. Phys. 50, 2733 (1969)JCPSA60021-960610.1063/1.1671437]. Important consequences for HS fluids in terms of sound waves propagation, Poisson's ratio, Stokes-Einstein relation, and generalized Cauchy identity are explored. Conventional expressions for high-frequency elastic moduli for simple systems with continuous and differentiable interatomic interaction potentials are known to diverge when approaching the HS repulsive limit. The origin of this divergence is identified here. It is demonstrated that these conventional expressions are only applicable for sufficiently soft interactions and should not be applied to HS systems. The reported results can be of interest in the context of statistical physics, physics of fluids, soft condensed matter, and granular materials.

Low-frequency Raman spectra of a glass-forming ionic liquid at low temperature and high pressure

The frequency range below ∼100 cm^-1 of the Raman spectrum of a glass-forming liquid exhibits two features that characterize the short-time (THz) dynamics: the quasi-elastic scattering (QES) tail and the boson peak (BP). In this work, we follow temperature and pressure effects on the intermolecular dynamics of a typical ionic liquid, 1-butyl-1-methylpiperidinium bis(trifluoromethanesulfonyl)imide, [Pip₁₄][[NTf₂]. The glass transition temperature of [Pip₁₄][[NTf₂] at atmospheric pressure is T_g = 198 K, and the pressure of glass transition at room temperature is P_g = 1.1 GPa. Raman spectra obtained while cooling the liquid or heating the glass exhibit hysteresis in QES and BP intensities, I_QES and I_BP. The dependence of I_QES, I_BP, and the BP frequency, ω_BP, with pressure up to the glass transition is steeper than the temperature dependence due to the stronger pressure effect on density within the GPa range. The temperature and pressure behaviors of the parameters I_QES, I_BP, and ω_BP obtained here for [Pip₁₄][[NTf₂] are discussed in light of known results for other glass-formers.

Communication: Glass transition and melting lines of an ionic liquid

The phase diagram of the ionic liquid 1-butyl-1-methylpyrrolidinium bis(trifluoromethanesufonyl)imide, [Pyrr_1,4][NTf₂], was explored by synchroton X-ray diffraction and Raman scattering measurements as a function of temperature and pressure. Glass transition T_g(p) and melting T_m(p) temperatures were obtained from atmospheric pressure up to ca. 2.0 GPa. We found that both the T_g(p) and T_m(p) curves follow essentially the same pressure dependence. The similarity of pressure coefficients, dT_g/dp ≈ dT_m/dp, is explained within the non-equilibrium thermodynamics approach for the glass transition by assuming that one of the Ehrenfest equations is appropriated for T_g(p), whereas T_m(p) follows the Clausius-Clapeyron equation valid for the first-order transitions. The results highlight that ionic liquids are excellent model systems to address fundamental questions related to the glass transition.

Highly Surface-Active Ca(OH)2 Monolayer as a CO2 Capture Material

Greenhouse gas emissions originating from fossil fuel combustion contribute significantly to global warming, and therefore the design of novel materials that efficiently capture CO₂ can play a crucial role in solving this challenge. Here, we show that reducing the dimensionality of bulk crystalline portlandite results in a stable monolayer material, named portlandene, that is highly effective at capturing CO₂. On the basis of theoretical analysis comprised of ab initio quantum mechanical calculations and force-field molecular dynamics simulations, we show that this single-layer phase is robust and maintains its stability even at high temperatures. The chemical activity of portlandene is seen to further increase upon defect engineering of its surface using vacancy sites. Defect-containing portlandene is capable of separating CO and CO₂ from a syngas (CO/CO₂/H₂) stream, yet is inert to water vapor. This selective behavior and the associated mechanisms have been elucidated by examining the electronic structure, local charge distribution, and bonding orbitals of portlandene. Additionally, unlike conventional CO₂ capturing technologies, the regeneration process of portlandene does not require high temperature heat treatment because it can release the captured CO₂ by application of a mild external electric field, making portlandene an ideal CO₂ capturing material for both pre- and postcombustion processes.

Structural relaxation and highly viscous flow

The highly viscous flow is due to thermally activated Eshelby transitions which transform a region of the undercooled liquid to a different structure with a different elastic misfit to the viscoelastic surroundings. A self-consistent determination of the viscosity in this picture explains why the average structural relaxation time is a factor of eight longer than the Maxwell time. The physical reason for the short Maxwell time is the very large contribution of strongly strained inherent states to the fluidity (the inverse viscosity). At the Maxwell time, the viscous no-return processes coexist with the back-and-forth jumping retardation processes.

Structural Origin of Enhanced Dynamics at the Surface of a Glassy Alloy

The enhancement of mobility at the surface of an amorphous alloy is studied using a combination of molecular dynamic simulations and normal mode analysis of the nonuniform distribution of Debye-Waller factors. The increased mobility at the surface is found to be associated with the appearance of Arrhenius temperature dependence. We show that the transverse Debye-Waller factor exhibits a peak at the surface. Over the accessible temperature range, we find that the bulk and surface diffusion coefficients obey the same empirical relationship with the respective Debye-Waller factors. Extrapolating this relationship to lower T, we argue that the observed decrease in the constraint at the surface is sufficient to account for the experimentally observed surface enhancement of mobility.

Universal structural parameter to quantitatively predict metallic glass properties

Quantitatively correlating the amorphous structure in metallic glasses (MGs) with their physical properties has been a long-sought goal. Here we introduce ‘flexibility volume' as a universal indicator, to bridge the structural state the MG is in with its properties, on both atomic and macroscopic levels. The flexibility volume combines static atomic volume with dynamics information via atomic vibrations that probe local configurational space and interaction between neighbouring atoms. We demonstrate that flexibility volume is a physically appropriate parameter that can quantitatively predict the shear modulus, which is at the heart of many key properties of MGs. Moreover, the new parameter correlates strongly with atomic packing topology, and also with the activation energy for thermally activated relaxation and the propensity for stress-driven shear transformations. These correlations are expected to be robust across a very wide range of MG compositions, processing conditions and length scales.

Various known structural descriptors of metallic glasses have limitations in quantitatively predicting properties. Here authors define a physically-motivated measure and show it to correlate strongly with elastic properties, local structure and relaxation kinetics over a wide range of simulated compositions.