-
Jammed solids with pins: Thresholds, Force networks and Elasticity
Authors:
Andy L. Zhang,
Sean A. Ridout,
Celia Parts,
Aarushi Sachdeva,
Cacey S. Bester,
Katharina Vollmayr-Lee,
Brian C. Utter,
Ted Brzinski,
Amy L. Graves
Abstract:
The role of fixed degrees of freedom in soft/granular matter systems has broad applicability and theoretical interest. Here we address questions of the geometrical role that a scaffolding of fixed particles plays in tuning the threshold volume fraction and force network in the vicinity of jamming. Our 2d simulated system consists of soft particles and fixed "pins", both of which harmonically repel…
▽ More
The role of fixed degrees of freedom in soft/granular matter systems has broad applicability and theoretical interest. Here we address questions of the geometrical role that a scaffolding of fixed particles plays in tuning the threshold volume fraction and force network in the vicinity of jamming. Our 2d simulated system consists of soft particles and fixed "pins", both of which harmonically repel overlaps. On one hand, we find that many of the critical scalings associated with jamming in the absence of pins continue to hold in the presence of even dense pin latices. On the other hand, the presence of pins lowers the jamming threshold, in a universal way at low pin densities and a geometry-dependent manner at high pin densities, producing packings with lower densities and fewer contacts between particles. The onset of strong lattice dependence coincides with the development of bond-orientational order. Furthermore, the presence of pins dramatically modifies the network of forces, with both unusually weak and unusually strong forces becoming more abundant. The spatial organization of this force network depends on pin geometry and is described in detail. Using persistent homology we demonstrate that pins modify the topology of the network. Finally, we observe clear signatures of this developing bond-orientational order and broad force distribution in the elastic moduli which characterize the linear response of these packings to strain.
△ Less
Submitted 25 August, 2022; v1 submitted 29 May, 2022;
originally announced May 2022.
-
Introduction to molecular dynamics simulations
Authors:
K. Vollmayr-Lee
Abstract:
We provide an introduction to molecular dynamics simulations in the context of the Kob-Andersen model of a glass. We introduce a complete set of tools for doing and analyzing the results of simulations at fixed NVE and NVT. The modular format of the paper allows readers to select sections that meet their needs. We start with an introduction to molecular dynamics independent of the programming lang…
▽ More
We provide an introduction to molecular dynamics simulations in the context of the Kob-Andersen model of a glass. We introduce a complete set of tools for doing and analyzing the results of simulations at fixed NVE and NVT. The modular format of the paper allows readers to select sections that meet their needs. We start with an introduction to molecular dynamics independent of the programming language, followed by introductions to an implementation using Python and then the freely available open source software package LAMMPS. We also describe analysis tools for the quick testing of the program during its development and compute the radial distribution function and the mean square displacement using both Python and LAMMPS.
△ Less
Submitted 21 May, 2020; v1 submitted 20 January, 2020;
originally announced January 2020.
-
Slow and Long-ranged Dynamical Heterogeneities in Dissipative Fluids
Authors:
Karina E. Avila,
Horacio E. Castillo,
Katharina Vollmayr-Lee,
Annette Zippelius
Abstract:
A two-dimensional bidisperse granular fluid is shown to exhibit pronounced long-ranged dynamical heterogeneities as dynamical arrest is approached. Here we focus on the most direct approach to study these heterogeneities: we identify clusters of slow particles and determine their size, $N_c$, and their radius of gyration, $R_G$. We show that $N_c\propto R_G^{d_f}$, providing direct evidence that t…
▽ More
A two-dimensional bidisperse granular fluid is shown to exhibit pronounced long-ranged dynamical heterogeneities as dynamical arrest is approached. Here we focus on the most direct approach to study these heterogeneities: we identify clusters of slow particles and determine their size, $N_c$, and their radius of gyration, $R_G$. We show that $N_c\propto R_G^{d_f}$, providing direct evidence that the most immobile particles arrange in fractal objects with a fractal dimension, $d_f$, that is observed to increase with packing fraction $φ$. The cluster size distribution obeys scaling, approaching an algebraic decay in the limit of structural arrest, i.e., $φ\toφ_c$. Alternatively, dynamical heterogeneities are analyzed via the four-point structure factor $S_4(q,t)$ and the dynamical susceptibility $χ_4(t)$. $S_4(q,t)$ is shown to obey scaling in the full range of packing fractions, $0.6\leqφ\leq 0.805$, and to become increasingly long-ranged as $φ\toφ_c$. Finite size scaling of $χ_4(t)$ provides a consistency check for the previously analyzed divergences of $χ_4(t)\propto (φ-φ_c)^{-γ_χ}$ and the correlation length $ξ\propto (φ-φ_c)^{-γ_ξ}$. We check the robustness of our results with respect to our definition of mobility. The divergences and the scaling for $φ\toφ_c$ suggest a non-equilibrium glass transition which seems qualitatively independent of the coefficient of restitution.
△ Less
Submitted 25 October, 2017;
originally announced October 2017.
-
Universal Scaling in the Aging of the Strong Glass Former SiO$_2$
Authors:
Katharina Vollmayr-Lee,
Christopher H. Gorman,
Horacio E. Castillo
Abstract:
We show that the aging dynamics of a strong glass former displays a strikingly simple scaling behavior, connecting the average dynamics with its fluctuations, namely the dynamical heterogeneities. We perform molecular dynamics simulations of SiO$_2$ with BKS interactions, quenching the system from high to low temperature, and study the evolution of the system as a function of the waiting time…
▽ More
We show that the aging dynamics of a strong glass former displays a strikingly simple scaling behavior, connecting the average dynamics with its fluctuations, namely the dynamical heterogeneities. We perform molecular dynamics simulations of SiO$_2$ with BKS interactions, quenching the system from high to low temperature, and study the evolution of the system as a function of the waiting time $t_{\rm w}$ measured from the instant of the quench. We find that both the aging behavior of the dynamic susceptibility $χ_4$ and the aging behavior of the probability distribution $P(f_{{\rm s},{\mathbf r}})$ of the local incoherent intermediate scattering function $f_{{\rm s},{\mathbf r}}$ can be described by simple scaling forms in terms of the global incoherent intermediate scattering function $C$. The scaling forms are the same that have been found to describe the aging of several fragile glass formers and that, in the case of $P(f_{{\rm s},{\mathbf r}})$, have been also predicted theoretically. A thorough study of the length scales involved highlights the importance of intermediate length scales. We also analyze directly the scaling dependence on particle type and on wavevector $q$, and find that both the average and the fluctuations of the slow aging dynamics are controlled by a unique aging clock, which is not only independent of the wavevector $q$, but is the same for O and Si atoms.
△ Less
Submitted 20 March, 2016;
originally announced March 2016.
-
Strong Dynamical Heterogeneity and Universal Scaling in Driven Granular Fluids
Authors:
Karina E. Avila,
Horacio E. Castillo,
Andrea Fiege,
Katharina Vollmayr-Lee,
Annette Zippelius
Abstract:
Large scale simulations of two-dimensional bidisperse granular fluids allow us to determine spatial correlations of slow particles via the four-point structure factor $S_4(q,t)$. Both cases, elastic ($\varepsilon=1$) as well as inelastic ($\varepsilon < 1$) collisions, are studied. As the fluid approaches structural arrest, i.e. for packing fractions in the range $0.6 \le φ\le 0.805$, scaling is s…
▽ More
Large scale simulations of two-dimensional bidisperse granular fluids allow us to determine spatial correlations of slow particles via the four-point structure factor $S_4(q,t)$. Both cases, elastic ($\varepsilon=1$) as well as inelastic ($\varepsilon < 1$) collisions, are studied. As the fluid approaches structural arrest, i.e. for packing fractions in the range $0.6 \le φ\le 0.805$, scaling is shown to hold: $S_4(q,t)/χ_4(t)=s(qξ(t))$. Both the dynamic susceptibility, $χ_4(τ_α)$, as well as the dynamic correlation length, $ξ(τ_α)$, evaluated at the $α$ relaxation time, $τ_α$, can be fitted to a power law divergence at a critical packing fraction. The measured $ξ(τ_α)$ widely exceeds the largest one previously observed for hard sphere 3d fluids. The number of particles in a slow cluster and the correlation length are related by a robust power law, $χ_4(τ_α) \approxξ^{d-p}(τ_α)$, with an exponent $d-p\approx 1.6$. This scaling is remarkably independent of $\varepsilon$, even though the strength of the dynamical heterogeneity increases dramatically as $\varepsilon$ grows.
△ Less
Submitted 17 April, 2014; v1 submitted 12 December, 2013;
originally announced December 2013.
-
Temperature-Dependent Defect Dynamics in the Network Glass SiO2
Authors:
Katharina Vollmayr-Lee,
Annette Zippelius
Abstract:
We investigate the long time dynamics of a strong glass former, SiO2, below the glass transition temperature by averaging single particle trajectories over time windows which comprise roughly 100 particle oscillations. The structure on this coarse-grained time scale is very well defined in terms of coordination numbers, allowing us to identify ill-coordinated atoms, called defects in the following…
▽ More
We investigate the long time dynamics of a strong glass former, SiO2, below the glass transition temperature by averaging single particle trajectories over time windows which comprise roughly 100 particle oscillations. The structure on this coarse-grained time scale is very well defined in terms of coordination numbers, allowing us to identify ill-coordinated atoms, called defects in the following. The most numerous defects are OO neighbors, whose lifetimes are comparable to the equilibration time at low temperature. On the other hand SiO and OSi defects are very rare and short lived. The lifetime of defects is found to be strongly temperature dependent, consistent with activated processes. Single-particle jumps give rise to local structural rearrangements. We show that in SiO2 these structural rearrangements are coupled to the creation or annihilation of defects, giving rise to very strong correlations of jumping atoms and defects.
△ Less
Submitted 8 July, 2013;
originally announced July 2013.
-
Microscopic picture of aging in SiO2
Authors:
K. Vollmayr-Lee,
R. Bjorkquist,
L. M. Chambers
Abstract:
We investigate the aging dynamics of amorphous SiO2 via molecular dynamics simulations of a quench from a high temperature T_i to a lower temperature T_f. We obtain a microscopic picture of aging dynamics by analyzing single particle trajectories, identifying jump events when a particle escapes the cage formed by its neighbors, and by determining how these jumps depend on the waiting time t_w, the…
▽ More
We investigate the aging dynamics of amorphous SiO2 via molecular dynamics simulations of a quench from a high temperature T_i to a lower temperature T_f. We obtain a microscopic picture of aging dynamics by analyzing single particle trajectories, identifying jump events when a particle escapes the cage formed by its neighbors, and by determining how these jumps depend on the waiting time t_w, the time elapsed since the temperature quench to T_f. We find that the only t_w-dependent microscopic quantity is the number of jumping particles per unit time, which decreases with age. Similar to previous studies for fragile glass formers, we show here for the strong glass former SiO2 that neither the distribution of jump lengths nor the distribution of times spent in the cage are t_w-dependent. We conclude that the microscopic aging dynamics is surprisingly similar for fragile and strong glass formers.
△ Less
Submitted 18 September, 2012;
originally announced September 2012.
-
Hydrodynamic Correlation Functions of a Driven Granular Fluid in Steady State
Authors:
Katharina Vollmayr-Lee,
Timo Aspelmeier,
Annette Zippelius
Abstract:
We study a homogeneously driven granular fluid of hard spheres at intermediate volume fractions and focus on time-delayed correlation functions in the stationary state. Inelastic collisions are modeled by incomplete normal restitution, allowing for efficient simulations with an event-driven algorithm. The incoherent scattering function, F_incoh(q,t), is seen to follow time-density superposition wi…
▽ More
We study a homogeneously driven granular fluid of hard spheres at intermediate volume fractions and focus on time-delayed correlation functions in the stationary state. Inelastic collisions are modeled by incomplete normal restitution, allowing for efficient simulations with an event-driven algorithm. The incoherent scattering function, F_incoh(q,t), is seen to follow time-density superposition with a relaxation time that increases significantly as volume fraction increases. The statistics of particle displacements is approximately Gaussian. For the coherent scattering function S(q,omega) we compare our results to the predictions of generalized fluctuating hydrodynamics which takes into account that temperature fluctuations decay either diffusively or with a finite relaxation rate, depending on wave number and inelasticity. For sufficiently small wave number q we observe sound waves in the coherent scattering function S(q,omega) and the longitudinal current correlation function C_l(q,omega). We determine the speed of sound and the transport coefficients and compare them to the results of kinetic theory.
△ Less
Submitted 27 July, 2010;
originally announced July 2010.
-
Aging to Equilibrium Dynamics of SiO2
Authors:
K. Vollmayr-Lee,
J. A. Roman,
J. Horbach
Abstract:
Molecular dynamics computer simulations are used to study the aging dynamics of SiO2 (modeled by the BKS model). Starting from fully equilibrated configurations at high temperatures T_i =5000K/3760K the system is quenched to lower temperatures T_f=2500K, 2750K, 3000K, 3250K and observed after a waiting time t_w. Since the simulation runs are long enough to reach equilibrium at T_f, we are able to…
▽ More
Molecular dynamics computer simulations are used to study the aging dynamics of SiO2 (modeled by the BKS model). Starting from fully equilibrated configurations at high temperatures T_i =5000K/3760K the system is quenched to lower temperatures T_f=2500K, 2750K, 3000K, 3250K and observed after a waiting time t_w. Since the simulation runs are long enough to reach equilibrium at T_f, we are able to study the transition from out-of-equilibrium to equilibrium dynamics. We present results for the partial structure factors, for the generalized incoherent intermediate scattering function C_q(t_w, t_w+t), and for the mean square displacement msd(t_w,t_w+t). We conclude that there are three different t_w regions: (I) At very short waiting times, C_q(t_w, t_w+t) decays very fast without forming a plateau. Similarly msd(t_w,t_w+t) increases without forming a plateau. (II) With increasing t_w a plateau develops in C_q(t_w, t_w+t) and msd(t_w,t_w+t). For intermediate waiting times the plateau height is independent of t_w and T_i. Time superposition applies, i.e. C_q=C_q(t/t_r) where t_r=t_r(t_w) is a waiting time dependent decay time. Furthermore C_q=C(q,t_w,t_w+t) scales as C_q=C(q,z(t_w,t) where z is a function of t_w and t only, i.e. independent of q. (III) At large t_w the system reaches equilibrium, i.e. C_q(t_w,t_w+t) and msd(t_w,t_w+t) are independent of t_w and T_i. For C_q(t_w,t_w+t) we find that the time superposition of intermediate waiting times (II) includes the equilibrium curve (III).
△ Less
Submitted 27 July, 2010; v1 submitted 5 January, 2010;
originally announced January 2010.
-
Self-Organized Criticality Below The Glass Transition
Authors:
Katharina Vollmayr-Lee,
Elizabeth A. Baker
Abstract:
We obtain evidence that the dynamics of glassy systems below the glass transition is characterized by self-organized criticality. Using molecular dynamics simulations of a model glass-former we identify clusters of cooperatively jumping particles. We find string-like clusters whose size is power-law distributed not only close to T_c but for ALL temperatures below T_c, indicating self-organized c…
▽ More
We obtain evidence that the dynamics of glassy systems below the glass transition is characterized by self-organized criticality. Using molecular dynamics simulations of a model glass-former we identify clusters of cooperatively jumping particles. We find string-like clusters whose size is power-law distributed not only close to T_c but for ALL temperatures below T_c, indicating self-organized criticality which we interpret as a freezing in of critical behavior.
△ Less
Submitted 26 March, 2006; v1 submitted 28 February, 2006;
originally announced February 2006.
-
Heterogeneities in the Glassy State
Authors:
Katharina Vollmayr-Lee,
Annette Zippelius
Abstract:
We study heterogeneities in a binary Lennard-Jones system below the glass transition using molecular dynamics simulations. We identify mobile and immobile particles and measure their distribution of vibrational amplitudes. For temperatures near the glass transition the distribution of vibrational amplitudes obeys scaling and compares reasonably well with a mean field theory for the amorphous sol…
▽ More
We study heterogeneities in a binary Lennard-Jones system below the glass transition using molecular dynamics simulations. We identify mobile and immobile particles and measure their distribution of vibrational amplitudes. For temperatures near the glass transition the distribution of vibrational amplitudes obeys scaling and compares reasonably well with a mean field theory for the amorphous solid state. To investigate correlations among the immobile and mobile particles we identify clusters and analyze their size and shape. For a fixed number of immobile particles we observe that the immobile particles cluster more strongly together as the temperature is increased which allows the particles to block each other more effectively, and to therefore stay immobile. For the mobile particles, on the other hand, the clustering is most pronounced at small temperatures, indicating that mobility at low temperatures can only be sustained in cooperative motion.
△ Less
Submitted 5 July, 2005;
originally announced July 2005.
-
Single Particle Jumps in a Binary Lennard-Jones System Below The Glass Transition
Authors:
K. Vollmayr-Lee
Abstract:
We study a binary Lennard-Jones system below the glass transition with molecular dynamics simulations. To investigate the dynamics we focus on events ("jumps") where a particle escapes the cage formed by its neighbors. Using single particle trajectories we define a jump by comparing for each particle its fluctuations with its changes in average position. We find two kinds of jumps: "reversible j…
▽ More
We study a binary Lennard-Jones system below the glass transition with molecular dynamics simulations. To investigate the dynamics we focus on events ("jumps") where a particle escapes the cage formed by its neighbors. Using single particle trajectories we define a jump by comparing for each particle its fluctuations with its changes in average position. We find two kinds of jumps: "reversible jumps," where a particle jumps back and forth between two or more average positions, and "irreversible jumps," where a particle does not return to any of its former average positions. For all investigated temperatures both kinds of particles jump and both irreversible and reversible jumps occur. With increasing temperature relaxation is enhanced by an increasing number of jumps, and growing jump lengths in position and potential energy. However, the waiting time between two successive jumps is independent of temperature. This temperature independence might be due to aging, which is present in our system. The ratio of irreversible to reversible jumps is also increasing with increasing temperature, which we interpret as a consequence of the increased likelihood of changes in the cages, i.e. a blocking of the "entrance" back into the previous cage. A comparison of the fluctuations of jumping particles and non-jumping particles indicates that jumping particles are more mobile even when not jumping. The jumps in energy normalized by their fluctuations are decreasing with increasing temperature, which is consistent with relaxation being increasingly driven by thermal fluctuations. In accordance with subdiffusive behavior are the distributions of waiting times and jump lengths in position.
△ Less
Submitted 28 August, 2003;
originally announced August 2003.
-
Dynamical Heterogeneities Below the Glass Transition
Authors:
Katharina Vollmayr-Lee,
Walter Kob,
Kurt Binder,
Annette Zippelius
Abstract:
We present molecular dynamics simulations of a binary Lennard-Jones mixture at temperatures below the kinetic glass transition. The ``mobility'' of a particle is characterized by the amplitude of its fluctuation around its average position. The 5% particles with the largest/smallest mean amplitude are thus defined as the relatively most mobile/immobile particles. We investigate for these 5% part…
▽ More
We present molecular dynamics simulations of a binary Lennard-Jones mixture at temperatures below the kinetic glass transition. The ``mobility'' of a particle is characterized by the amplitude of its fluctuation around its average position. The 5% particles with the largest/smallest mean amplitude are thus defined as the relatively most mobile/immobile particles. We investigate for these 5% particles their spatial distribution and find them to be distributed very heterogeneously in that mobile as well as immobile particles form clusters. The reason for this dynamic heterogeneity is traced back to the fact that mobile/immobile particles are surrounded by fewer/more neighbors which form an effectively wider/narrower cage. The dependence of our results on the length of the simulation run indicates that individual particles have a characteristic mobility time scale, which can be approximated via the non-Gaussian parameter.
△ Less
Submitted 25 September, 2001;
originally announced September 2001.
-
Using mean field theory to determine the structure of uniform fluids
Authors:
Katharina Vollmayr-Lee,
Kirill Katsov,
John D. Weeks
Abstract:
The structure of a uniform simple liquid is related to that of a reference fluid with purely repulsive intermolecular forces in a self-consistently determined external reference field (ERF) phi_ R. The ERF can be separated into a harshly repulsive part phi_ R0 generated by the repulsive core of a reference particle fixed at the origin and a more slowly varying part phi_ R1 arising from a mean fi…
▽ More
The structure of a uniform simple liquid is related to that of a reference fluid with purely repulsive intermolecular forces in a self-consistently determined external reference field (ERF) phi_ R. The ERF can be separated into a harshly repulsive part phi_ R0 generated by the repulsive core of a reference particle fixed at the origin and a more slowly varying part phi_ R1 arising from a mean field treatment of the attractive forces. We use a generalized linear response method to calculate the reference fluid structure, first determining the response to the smoother part phi_ R1 of the ERF alone, followed by the response to the harshly repulsive part. Both steps can be carried out very accurately, as confirmed by MD simulations, and good agreement with the structure of the full LJ fluid is found.
△ Less
Submitted 24 August, 2000;
originally announced August 2000.
-
Cooling Rate Dependence and Dynamic Heterogeneity Below the Glass Transition in a Lennard-Jones Glass
Authors:
K. Vollmayr-Lee,
W. Kob,
K. Binder,
A. Zippelius
Abstract:
We investigate a binary Lennard-Jones mixture with molecular dynamics simulations. We consider first a system cooled linearly in time with the cooling rate gamma. By varying gamma over almost four decades we study the influence of the cooling rate on the glass transition and on the resulting glass. We find for all investigated quantities a cooling rate dependence; with decreasing cooling rate th…
▽ More
We investigate a binary Lennard-Jones mixture with molecular dynamics simulations. We consider first a system cooled linearly in time with the cooling rate gamma. By varying gamma over almost four decades we study the influence of the cooling rate on the glass transition and on the resulting glass. We find for all investigated quantities a cooling rate dependence; with decreasing cooling rate the system falls out of equilibrium at decreasing temperatures, reaches lower enthalpies and obtains increasing local order. Next we study the dynamics of the melting process by investigating the most immobile and most mobile particles in the glass. We find that their spatial distribution is heterogeneous and that the immobile/mobile particles are surrounded by denser/less dense cages than an average particle.
△ Less
Submitted 20 December, 1999;
originally announced December 1999.