Influence of string-like cooperative atomic motion on surface diffusion in the (110) interfacial region of crystalline Ni

Collective motion and its connection to the energy landscape in 2D soft crystals

We examine the collective motion in computational models of a two-dimensional dusty plasma crystal and a charged colloidal suspension as they approach their respective melting transitions. To unambiguously identify rearrangement events in the crystal, we map the trajectory of configurations from an equilibrium molecular dynamics simulation to the corresponding sequence of configurations of local potential energy minima ("inherent structures"). This inherent structure (IS) trajectory eliminates the ambiguity that arises from localized vibrational motion. We find that the evolution of the IS trajectory in the crystal can be split into comparatively longer-lived ground states and shorter-lived discrete excited states. These discrete excited energy levels are a consequence of discrete numbers of defect clusters in the crystal. We find that the collective rearrangement occurs through different mechanisms: (i) small closed-loop motion in the ground states without the facilitation of defects, and (ii) much larger and complex open-ended particle motions in excited states that are facilitated by clusters of defects. In both cases, clusters of displacing particles can be separated into much smaller groups of replacing particles with a loop-like structure. In contrast to glass-forming liquids, the mass of the rearranging groups grows on heating towards the melting temperature rather than cooling. We find that crystal melting in these systems can be anticipated by the merging of the average time the crystal spends in the ground state with the average time in the excited states.

A dynamical system approach to relaxation in glass-forming liquids

The "classical" thermodynamic and statistical mechanical theories of Gibbs and Boltzmann are both predicated on axiomatic assumptions whose applicability is hard to ascertain. Theoretical objections and an increasing number of observed deviations from these theories have led to sustained efforts to develop an improved mathematical and physical foundation for them, and the search for appropriate extensions that are generally applicable to condensed materials at low temperatures (T) and high material densities where the assumptions of these theories start to become particularly questionable. These theoretical efforts have largely focused on minimal models of condensed material systems, such as the Fermi-Ulam-Pasta-Tsingou model, and other simplified models of condensed materials that are amenable to numerical and analytic treatments and that can serve to illuminate essential features of relaxation processes in condensed materials under conditions approaching integrable dynamics where clear departures from classical thermodynamics and dynamics can be generally expected. These studies indicate an apparently general multi-step relaxation process, corresponding to an initial "fast" relaxation process (termed the fast β-relaxation in the context of cooled liquids), followed by a longer "equipartition time", namely, the α-relaxation time τα in the context of cooled liquids. This relaxation timescale can be enormously longer than the fast β-relaxation time τβ so that τα is the primary parameter governing the rate at which the material comes into equilibrium, and thus is a natural focus of theoretical attention. Since the dynamics of these simplified dynamical systems, originally intended as simplified models of real crystalline materials exhibiting anharmonic interactions, greatly resemble the observed relaxation dynamics of both heated crystals and cooled liquids, we adapt this dynamical system approach to the practical matter of estimating relaxation times in both cooled liquids and crystals at elevated temperatures, which we identify as weakly non-integrable dynamical systems.

Understanding Relaxation in the Kob-Andersen Liquid Based on Entropy, String, Shoving, Localization, and Parabolic Models

We assess the validity of a range of models of glass formation based on molecular dynamics simulation results of the Kob-Andersen (KA) model system under a wide range of constant volume and constant pressure conditions. These models include the Adam-Gibbs model emphasizing configurational entropy, the string model emphasizing collective particle exchange motion, the shoving model emphasizing material elasticity, the localization model emphasizing dynamical free volume, and parabolic models based on the ideas of dynamic facilitation and, alternatively, the hypothesis that glass formation involves an avoided critical point. We demonstrate that these seemingly disparate models all provide a reasonable description of structural relaxation and diffusion data for the KA model system under all simulation conditions considered. Hence, the present study points to some unity in our understanding of the relationship between leading models of glass formation, supporting inferences drawn from previous studies of polymeric glass-forming liquids.

Influence of Density and Pressure on Glass Formation in the Kob-Andersen Model

We systematically study glass formation in the well-known Kob-Andersen model over a wide range of densities and pressures as a basis for judging the "universality" of glass formation through a comparison to a recent systematic study of model polymeric glass-forming liquids. Our purpose is to establish general characteristics of glass formation, to identify new relations, and to discern which properties of glass-forming liquids are material-specific and which are "universal." To this end, we analyze a number of characteristic properties of glass formation, such as the structural relaxation time, self-diffusion coefficient, viscosity, characteristic temperatures, and fragility. We also consider a suite of properties presumably related to dynamic heterogeneity in an attempt to better understand its relation to structural relaxation. We demonstrate that the glassy dynamics in the Kob-Andersen model exhibit many of the essential trends observed in polymeric glass-forming liquids, pointing to a remarkable "universality" of many aspects of glass formation.

A closer examination of the nature of atomic motion in the interfacial region of crystals upon approaching melting

Although crystalline materials are often conceptualized as involving a static lattice configuration of particles, it has recently become appreciated that string-like collective particle exchange motion is a ubiquitous and physically important phenomenon in both the melting and interfacial dynamics of crystals. This type of collective motion has been evidenced in melting since early simulations of hard disc melting by Alder et al. [Phys. Rev. Lett. 11(6), 241-243 (1963)], but a general understanding of its origin, along with its impact on melting and the dynamics of crystalline materials, has been rather slow to develop. We explore this phenomenon further by focusing on the interfacial dynamics of a model crystalline Cu material using molecular dynamics simulations where we emphasize the geometrical nature and spatial extent of the atomic trajectories over the timescale that they are caged, and we also quantify string-like collective motion on the timescale of the fast β-relaxation time, τf, i.e., "stringlets." Direct visualization of the atomic trajectories in their cages over the timescale over which the cage persists indicates that they become progressively more anisotropic upon approaching the melting temperature Tm. The stringlets, dominating the large amplitude atomic motion in the fast dynamics regime, are largely localized to the crystal interfacial region and correspond to "excess" modes in the density of states that give rise to a "boson peak." Moreover, interstitial point defects occur in direct association with the stringlets, demonstrating a link between classical defect models of melting and more recent studies of melting emphasizing the role of this kind of collective motion.

Molecular Dynamics Simulation of the Influence of Temperature and Salt on the Dynamic Hydration Layer in a Model Polyzwitterionic Polymer PAEDAPS

We investigate the hydration of poly(3-[2-(acrylamido) ethyldimethylammonio] propanesulfonate) over a range of temperatures in pure water and with the inclusion of 0.1 mol/L NaCl using atomistic molecular dynamics simulation. Drawing on concepts drawn from the field of glass-forming liquids, we use the Debye-Waller parameter () for describing the water mobility gradient around the polybetaine backbone extending to an overall distance ≈18 Å. The water mobility in this layer is defined through the mean-square water molecule displacement at a time on the order of water's β-relaxation time. The brushlike topology of polybetaines leads to two regions in the dynamic hydration layer. The inner region of ≈10.5 Å is explored by pendant group conformational motions, and the outer region of ≈7.5 Å represents an extended layer of reduced water mobility relative to bulk water. The dynamic hydration layer extends far beyond the static hydration layer, adjacent to the polymer.

Localization model description of the interfacial dynamics of crystalline Cu and Cu 64 Zr 36 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< u 2 > on a ps 'caging' timescale relative to the square of the average interparticle distanceσ 2 , 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τ α , 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τ α to D. In the current study, we employed classical molecular dynamics simulation to investigate the structural relaxation and diffusion of modelCu 64 Zr 36 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 (Cu 64 Zr 36 ) NPs to a good approximation, further confirming the generality of the model.

Factors correlating to enhanced surface diffusion in metallic glasses

The enhancement of surface diffusion (D_S) over the bulk (D_V) in metallic glasses (MGs) is well documented and likely to strongly influence the properties of glasses grown by vapor deposition. Here, we use classical molecular dynamics (MD) simulations to identify different factors influencing the enhancement of surface diffusion in MGs. MGs have a simple atomic structure and belong to the category of moderately fragile glasses that undergo pronounced slowdown of bulk dynamics with cooling close to the glass transition temperature (T_g). We observe that D_S exhibits a much more moderate slowdown compared to D_V when approaching T_g, and D_S/D_V at T_g varies by two orders of magnitude among the MGs investigated. We demonstrate that both the surface energy and the fraction of missing bonds for surface atoms show good correlation to D_S/D_V, implying that the loss of nearest neighbors at the surface directly translates into higher mobility, unlike the behavior of network-bonded and hydrogen-bonded organic glasses. Fragility, a measure of the slowdown of bulk dynamics close to T_g, also correlates to D_S/D_V, with more fragile systems having larger surface enhancement of mobility. The deviations observed in the fragility-D_S/D_V relationship are shown to be correlated to the extent of segregation or depletion of the mobile element at the surface. Finally, we explore the relationship between the diffusion pre-exponential factor (D₀) and the activation energy (Q) and compare it to a ln(D₀)-Q correlation previously established for bulk glasses, demonstrating similar correlations from MD as in the experiments and that the surface and bulk have very similar ln(D₀)-Q correlations.

Localization model description of the interfacial dynamics of crystalline Cu and Cu64Zr36 metallic glass films

Recent studies of structural relaxation in Cu-Zr metallic glass materials having a range of compositions and over a wide range of temperatures and in crystalline UO₂ under superionic conditions have indicated that the localization model (LM) can predict the structural relaxation time τ_α of these materials from the intermediate scattering function without any free parameters from the particle mean square displacement ⟨r²⟩ at a caging time on the order of ps, i.e., the "Debye-Waller factor" (DWF). In the present work, we test whether this remarkable relation between the "fast" picosecond dynamics and the rate of structural relaxation τ_α in these model amorphous and crystalline materials can be extended to the prediction of the local interfacial dynamics of model amorphous and crystalline films. Specifically, we simulate the free-standing amorphous Cu₆₄Zr₃₆ and crystalline Cu films and find that the LM provides an excellent parameter-free prediction for τ_α of the interfacial region. We also show that the Tammann temperature, defining the initial formation of a mobile interfacial layer, can be estimated precisely for both crystalline and glass-forming solid materials from the condition that the DWFs of the interfacial region and the material interior coincide.

Universal nature of dynamic heterogeneity in glass-forming liquids: A comparative study of metallic and polymeric glass-forming liquids

Glass-formation is a ubiquitous phenomenon that is often observed in a broad class of materials ranging from biological matter to commonly encountered synthetic polymer, as well as metallic and inorganic glass-forming (GF) materials. Despite the many regularities in the dynamical properties of GF materials, the structural origin of the universal dynamical properties of these materials has not yet been identified. Recent simulations of coarse-grained polymeric GF liquids have indicated the coexistence of clusters of mobile and immobile particles that appear to be directly linked, respectively, to the rate of molecular diffusion and structural relaxation. The present work examines the extent to which these distinct types of "dynamic heterogeneity" (DH) arise in metallic GF liquids (Cu-Zr, Ni-Nb, and Pd-Si alloys) having a vastly different molecular structure and chemistry. We first identified mobile and immobile particles and their transient clusters and found the DH in the metallic alloys to be remarkably similar in form to polymeric GF liquids, confirming the "universality" of the DH phenomenon. Furthermore, the lifetime of the mobile particle clusters was found to be directly related to the rate of diffusion in these materials, while the lifetime of immobile particles was found to be proportional to the structural relaxation time, providing some insight into the origin of decoupling in GF liquids. An examination of particles having a locally preferred atomic packing, and clusters of such particles, suggests that there is no one-to-one relation between these populations of particles so that an understanding of the origin of DH in terms of static fluid structure remains elusive.

Collective Motion in the Interfacial and Interior Regions of Supported Polymer Films and Its Relation to Relaxation

To understand the role of collective motion in the often large changes in interfacial molecular mobility observed in polymer films, we investigate the extent of collective motion in the interfacial regions of a thin supported polymer film and within the film interior by molecular dynamics simulation. Contrary to commonly stated expectations, we find that the extent of collective motion, as quantified by string-like molecular exchange motion, is similar in magnitude in the polymer-air interfacial layer as the film interior and distinct from the bulk material. This finding is consistent with Adam-Gibbs description of the segmental dynamics within mesoscopic film regions, where the extent of collective motion is related to the configurational entropy of the film as a whole rather than a locally defined extent of collective motion or configurational entropy.

Comparative Study of the Collective Dynamics of Proteins and Inorganic Nanoparticles

Molecular dynamics simulations of ubiquitin in water/glycerol solutions are used to test the suggestion by Karplus and coworkers that proteins in their biologically active state should exhibit a dynamics similar to 'surface-melted' inorganic nanoparticles (NPs). Motivated by recent studies indicating that surface-melted inorganic NPs are in a 'glassy' state that is an intermediate dynamical state between a solid and liquid, we probe the validity and significance of this proposed analogy. In particular, atomistic simulations of ubiquitin in solution based on CHARMM36 force field and pre-melted Ni NPs (Voter-Chen Embedded Atom Method potential) indicate a common dynamic heterogeneity, along with other features of glass-forming (GF) liquids such as collective atomic motion in the form of string-like atomic displacements, potential energy fluctuations and particle displacements with long range correlations ('colored' or 'pink' noise), and particle displacement events having a power law scaling in magnitude, as found in earthquakes. On the other hand, we find the dynamics of ubiquitin to be even more like a polycrystalline material in which the α-helix and β-sheet regions of the protein are similar to crystal grains so that the string-like collective atomic motion is concentrated in regions between the α-helix and β-sheet domains.