Source author record

Elijah Flenner

Elijah Flenner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

Microscopic analysis of sound attenuation in low-temperature amorphous solids reveals quantitative importance of non-affine effects

Sound attenuation in low temperature amorphous solids originates from their disordered structure. However, its detailed mechanism is still being debated. Here we analyze sound attenuation starting directly from the microscopic equations of motion. We derive an exact expression for the zero-temperature sound damping coefficient. We verify that the sound damping coefficients calculated from our expression agree very well with results from independent simulations of sound attenuation. The small wavevector analysis of our expression shows that sound attenuation is primarily determined by the non-affine displacements' contribution to the sound wave propagation coefficient coming from the frequency shell of the sound wave. Our expression involves only quantities that pertain to solids' static configurations. It can be used to evaluate the low temperature sound damping coefficients without directly simulating sound attenuation.

preprint2021arXiv

Single particle fluctuations dominate the long-time dynamic susceptibility in glass-forming liquids

Liquids near the glass transition exhibit dynamical heterogeneity, i.e. correlated regions in the liquid relax at either a much faster rate or a much slower rate than the average. This collective phenomenon has been characterized by measurements of a dynamic susceptibility $χ_4(t)$, which are sometimes interpreted in terms of the size of those relaxing regions and the intensity of the fluctuations. We show that the results of those measurements can be affected not only by the collective fluctuations in the relaxation rate, but also by density fluctuations in the initial state and by single-particle fluctuations. We also show that at very long times the average overlap $C(t)$ probing the similarity between an initial and a final state separated by a time interval $t$ decays as a power law $C(t) \sim t^{-d/2}$. This is much slower than the stretched exponential behavior $C(t) \sim {\rm e}^{-(t/τ)^β}$ previously observed at times within one or two orders of magnitude of the $α$-relaxation time $τ_α$. We find that for times longer than $10-100 τ_α$, the dynamic susceptibility $χ_4(t)$ is dominated by single particle fluctuations, and that $χ_4(t) \approx C(t) \sim t^{-d/2}$. Finally, we introduce a method to extract the collective relaxation contribution to the dynamic susceptibility $χ_4(t)$ by subtracting the effects of single-particle fluctuations and initial state density fluctuations. We apply this method to numerical simulations of two glass forming models: a binary hard sphere system and a Kob-Andersen Lennard-Jones system. This allows us to extend the analysis of numerical data to timescales much longer than previously possible, and opens the door for further future progress in the study of dynamic heterogeneities, including the determination of the exchange time.

preprint2020arXiv

Active matter: quantifying the departure from equilibrium

Active matter systems are driven out of equilibrium at the level of individual constituents. One widely studied class are systems of athermal particles that move under the combined influence of interparticle interactions and self-propulsions, with the latter evolving according to the Ornstein-Uhlenbeck stochastic process. Intuitively, these so-called active Ornstein-Uhlenbeck particles (AOUPs) systems are farther from equilibrium for longer self-propulsion persistence times. Quantitatively, this is confirmed by the increasing equal-time velocity correlations (which are trivial in equilibrium) and by the increasing violation of the Einstein relation between the self-diffusion and mobility coefficients. In contrast, the entropy production rate, calculated from the ratio of the probabilities of the position space trajectory and its time-reversed counterpart, has a non-monotonic dependence on the persistence time. Thus, it does not properly quantify the departure of AOUPs systems from equilibrium.

preprint2020arXiv

Sound attenuation in finite-temperature stable glasses

The temperature dependence of the thermal conductivity of amorphous solids is markedly different from that of their crystalline counterparts, but exhibits universal behaviour. Sound attenuation is believed to be related to this universal behaviour. Recent computer simulations demonstrated that in the harmonic approximation sound attenuation $Γ$ obeys quartic, Rayleigh scattering scaling for small wavevectors $k$ and quadratic scaling for wavevectors above the Ioffe-Regel limit. However, simulations and experiments do not provide a clear picture of what to expect at finite temperatures where anharmonic effects become relevant. Here we study sound attenuation at finite temperatures for model glasses of various stability, from unstable glasses that exhibit rapid aging to glasses whose stability is equal to those created in laboratory experiments. We find several scaling laws depending on the temperature and stability of the glass. First, we find the large wavevector quadratic scaling to be unchanged at all temperatures. Second, we find that at small wavectors $Γ\sim k^{1.5}$ for an aging glass, but $Γ\sim k^2$ when the glass does not age on the timescale of the calculation. For our most stable glass, we find that $Γ\sim k^2$ at small wavevectors, then a crossover to Rayleigh scattering scaling $Γ\sim k^4$, followed by another crossover to the quadratic scaling at large wavevectors. Our computational observation of this quadratic behavior reconciles simulation, theory and experiment, and will advance the understanding of the temperature dependence of thermal conductivity of glasses.

preprint2016arXiv

Dynamic heterogeneity in two-dimensional supercooled liquids: comparison of bond-breaking and bond-orientational correlations

We compare the spatial correlations of bond-breaking events and bond-orientational relaxation in a model two-dimensional liquid undergoing Newtonian dynamics. We find that the relaxation time of the bond-breaking correlation function is much longer than the relaxation time of the bond-orientational correlation function and self-intermediate scattering function. However, the relaxation time of the bond-orientational correlation function increases faster with decreasing temperature than the relaxation time of the bond-breaking correlation function and the self-intermediate scattering function. Moreover, the dynamic correlation length that characterizes the size of correlated bond-orientational relaxation grows faster with decreasing temperature than the dynamic correlation length that characterizes the size of correlated bond-breaking events. We also examine the ensemble-dependent and ensemble-independent dynamic susceptibilities for both bond-breaking correlations and bond-orientational correlations. We find that for both correlations, the ensemble-dependent and ensemble-independent susceptibilities exhibit a maximum at nearly the same time, and this maximum occurs at a time slightly shorter than the peak position of the dynamic correlation length.

preprint2016arXiv

Kinetic stability and energetics of simulated glasses created by constant pressure cooling

We use computer simulations to study the cooling rate dependence of the stability and energetics of model glasses created at constant pressure conditions and compare the results with glasses formed at constant volume conditions. To examine the stability, we determine the time it takes for a glass cooled and reheated at constant pressure to transform back into a liquid, $t_{\mathrm{trans}}$, and calculate the stability ratio $S = t_{\mathrm{trans}}/τ_α$, where $τ_α$ is the equilibrium relaxation time of the liquid. We find that, for slow enough cooling rates, cooling and reheating at constant pressure results in a larger stability ratio $S$ than for cooling and reheating at constant volume. We also compare the energetics of glasses obtained by cooling while maintaining constant pressure with those of glasses created by cooling from the same state point while maintaining constant volume. We find that cooling at constant pressure results in glasses with lower average potential energy and average inherent structure energy. We note that in model simulations of the vapor deposition process glasses are created under constant pressure conditions, and thus they should be compared to glasses obtained by constant pressure cooling.

preprint2016arXiv

The nonequilibrium glassy dynamics of self-propelled particles

We study the glassy dynamics taking place in dense assemblies of athermal active particles that are driven solely by a nonequilibrium self-propulsion mechanism. Active forces are modeled as an Ornstein-Uhlenbeck stochastic process, characterized by a persistence time and an effective temperature, and particles interact via a Lennard-Jones potential that yields well-studied glassy behavior in the Brownian limit, obtained as the persistence time vanishes. By increasing the persistence time, the system departs more strongly from thermal equilibrium and we provide a comprehensive numerical analysis of the structure and dynamics of the resulting active fluid. Finite persistence times profoundly affect the static structure of the fluid and give rise to nonequilibrium velocity correlations that are absent in thermal systems. Despite these nonequilibrium features, for any value of the persistence time we observe a nonequilibrium glass transition as the effective temperature is decreased. Surprisingly, increasing departure from thermal equilibrium is found to promote (rather than suppress) the glassy dynamics. Overall, our results suggest that with increasing persistence time, microscopic properties of the active fluid change quantitatively, but the broad features of the nonequilibrium glassy dynamics observed with decreasing the effective temperature remain qualitatively similar to those of thermal glass-formers.

preprint2015arXiv

Cooling-rate dependence of kinetic and mechanical stability of simulated glasses

Recently, ultrastable glasses have been created through vapor deposition. Subsequently, computer simulation algorithms have been proposed that mimic the vapor deposition process and result in simulated glasses with increased stability. In addition, random pinning has been used to generate very stable glassy configurations without the need for lengthy annealing or special algorithms inspired by vapor deposition. Kinetic and mechanical stability of experimental ultrastable glasses is compared to those of experimental glasses formed by cooling. We provide the basis for a similar comparison for simulated stable glasses: we analyze the kinetic and mechanical stability of simulated glasses formed by cooling at a constant rate by examining the transformation time to a liquid upon rapid re-heating, the inherent structure energies, and the shear modulus. The kinetic and structural stability increases slowly with decreasing cooling rate. The methods outlined here can be used to assess kinetic and mechanical stability of simulated glasses generated by using specialized algorithms.

preprint2015arXiv

Glassy dynamics of athermal self-propelled particles: Computer simulations and a nonequilibrium microscopic theory

We combine computer simulations and analytical theory to investigate the glassy dynamics in dense assemblies of athermal particles evolving under the sole influence of self-propulsion. The simulations reveal that when the persistence time of the self-propelled particles is increased, the local structure becomes more pronounced whereas the long-time dynamics first accelerates and then slows down. These seemingly contradictory evolutions are explained by constructing a nonequilibrium mode-coupling-like theory for interacting self-propelled particles. To predict the collective dynamics the theory needs the steady state structure factor and the steady state correlations of the local velocities. It yields nontrivial predictions for the glassy dynamics of self-propelled particles in qualitative agreement with the simulations.

preprint2015arXiv

Reduced strength and extent of dynamic heterogeneity in a strong glass former as compared to fragile glass formers

We examined dynamic heterogeneity in a model tetrahedral network glass-forming liquid. We used four-point correlation functions to extract dynamic correlation lengths xi_4^a(t) and susceptibilities chi_4^a(t) corresponding to structural relaxation on two length scales a. One length scale corresponds to structural relaxation at nearest neighbor distances and the other corresponds to relaxation of the tetrahedral structure. We find that the dynamic correlation length xi_4^{a} grows much slower with increasing relaxation time than for model fragile glass formers. We also find that chi_4^a ~ (xi_4^a)^z for a range of temperatures, but z < 3 at the lowest temperatures examined in this study. However, we do find evidence that the temperature where Stokes-Einstein violation begins marks a temperature where there is a change in the character of dynamically heterogeneous regions. Throughout the paper, we contrast the structure and dynamics of a strong glass former with that of a representative fragile glass former.

preprint2014arXiv

Long-range spatial correlations of particle displacements and the emergence of elasticity

We examine correlations of transverse particle displacements and their relationship to the shear modulus of a glass and the viscosity of a fluid. To this end we use computer simulations to calculate a correlation function of the displacements, $S_4(q;t)$, which is similar to functions used to study heterogeneous dynamics in glass-forming fluids. We show that in the glass the shear modulus can be obtained from the long-time, small-q limit of $S_4(q;t)$. By using scaling arguments, we argue that a four-point correlation length $ξ_4(t)$ grows linearly in time in a glass and grows as $\sqrt{t}$ at long times in a fluid, and we verify these results by analyzing $S_4(q;t)$ obtained from simulations. For a viscoelastic fluid, the simulation results suggest that the crossover to the long-time $\sqrt{t}$ growth of $ξ_4(t)$ occurs at a characteristic decay time of the shear stress autocorrelation function. Using this observation, we show that the amplitude of the long-time $\sqrt{t}$ growth is proportional to $\sqrtη$ where $η$ is the viscosity of the fluid.

preprint2013arXiv

Breakdown of a renormalized perturbation expansion around mode-coupling theory of the glass transition

We analyze a renormalized perturbation expansion around the mode-coupling theory of the glass transition. We focus on the long-time limit of the irreducible memory function. We discuss a renormalized diagrammatic expansion for this function and re-sum two infinite classes of diagrams. We show that the resulting contributions to the irreducible memory function diverge at the mode-coupling transition. A further re-summation of ladder diagrams constructed by iterating these divergent contributions gives a finite result which cancels the mode-coupling theory's expression for the irreducible memory function.

preprint2013arXiv

Universal features of dynamic heterogeneity in supercooled liquids

A few years ago it was showed that some systems that have very similar local structure, as quantified by the pair correlation function, exhibit vastly different slowing down upon supercooling [L. Berthier and G. Tarjus, Phys. Rev. Lett. 103, 170601 (2009); U.R. Pedersen, T.B. Schroder and J.C. Dyre, Phys. Rev. Lett. 105, 157801 (2010)]. Recently, a more subtle structural quantity, the so-called "point-to-set" length, was found to reliably correlate with the average dynamics [G.M. Hocky, T.E. Markland and D.R. Reichman, Phys. Rev. Lett. 108, 225506 (2012)]. Here we use computer simulations to examine the behavior of fluctuations around the average dynamics, i.e., dynamic heterogeneity. We study five model glass-forming liquids: three model liquids used in previous works and two additional model liquids with finite range interactions. Some of these systems have very similar local structure but vastly different dynamics. We show that for all these systems the spatial extent and the anisotropy of dynamic heterogeneity correlate very well with the average dynamics.

preprint2012arXiv

Dynamic heterogeneities above and below the mode-coupling temperature. Evidence of a dynamic crossover

We examine dynamic heterogeneities in a model glass-forming fluid, a binary harmonic sphere mixture, above and below the mode-coupling temperature T_c. We calculate the ensemble independent susceptibility chi_4(tau_alpha) and the dynamic correlation length xi_4(tau_alpha) at the alpha-relaxation time tau_alpha. We also examine in detail the temperature dependence of tau_alpha and the diffusion coefficient D. For higher temperatures we find that the standard Stokes-Einstein relationship, D ~ tau_alpha^{-1}, holds, but at lower temperatures a fractional Stokes-Einstein relationship, D ~ tau_alpha^{-sigma} with sigma = 0.65, emerges. By examining the relationships between tau_alpha, D, chi_4(tau_alpha), and xi_4(tau_alpha) we determine that the emergence of the fractional Stokes-Einstein relationship is accompanied by a dynamic crossover from tau_alpha ~ e^{k_2 xi_4} at higher temperatures to tau_alpha ~ e^{k_1 xi_4^{3/2}} at lower temperatures.

preprint2012arXiv

Glassy dynamics of partially pinned fluids: an alternative mode-coupling approach

We use a simple mode-coupling approach to investigate glassy dynamics of partially pinned fluid systems. Our approach is different from the mode-coupling theory developed by Krakoviack [Phys. Rev. Lett. 94, 065703 (2005), Phys. Rev. E 84, 050501(R) (2011)]. In contrast to Krakoviack's theory, our approach predicts a random pinning glass transition scenario that is qualitatively the same as the scenario obtained using a mean-field analysis of the spherical p-spin model and a mean-field version of the random first-order transition theory. We use our approach to calculate quantities which are often considered to be indicators of growing dynamic correlations and static point-to-set correlations. We find that the so-called static overlap is dominated by the simple, low pinning fraction contribution. Thus, at least for randomly pinned fluid systems, only a careful quantitative analysis of simulation results can reveal genuine, many-body point-to-set correlations.

preprint2011arXiv

Analysis of a growing dynamic length scale in a glass-forming binary hard-sphere mixture

We examine a length scale that characterizes the spatial extent of heterogeneous dynamics in a glass-forming binary hard-sphere mixture up to the mode-coupling volume fraction phi_c. First, we characterize the system's dynamics. Then, we utilize a new method [Phys. Rev. Lett. 105, 217801 (2010)] to extract and analyze the ensemble independent dynamic susceptibility chi_4(t) and the dynamic correlation length xi(t) for a range of times between the beta and alpha relaxation times. We find that in this time range the dynamic correlation length follows a volume fraction independent curve xi(t) ~ ln(t). For longer times, xi(t) departs from this curve and remains constant up to the largest time at which we can determine the length accurately. In addition to the previously established correlation tau_alpha ~ exp[xi(tau_alpha)] between the alpha relaxation time, tau_alpha, and the dynamic correlation length at this time, xi(tau_alpha), we also find a similar correlation for the diffusion coefficient D ~ exp[xi(tau_alpha)^theta] with theta approximately 0.6. We discuss the relevance of these findings for different theories of the glass transition.

preprint2011arXiv

Emergence of Long-range Correlations and Rigidity at the Dynamic Glass Transition

At the microscopic level, equilibrium liquid's translational symmetry is spontaneously broken at the so-called dynamic glass transition predicted by the mean-field replica approach. We show that this fact implies the emergence of Goldstone modes and long-range density correlations. We derive and evaluate a new statistical mechanical expression for the glass shear modulus.

preprint2011arXiv

Kinetic Monte Carlo and Cellular Particle Dynamics Simulations of Multicellular Systems

Computer modeling of multicellular systems has been a valuable tool for interpreting and guiding in vitro experiments relevant to embryonic morphogenesis, tumor growth, angiogenesis and, lately, structure formation following the printing of cell aggregates as bioink particles. Computer simulations based on Metropolis Monte Carlo (MMC) algorithms were successful in explaining and predicting the resulting stationary structures (corresponding to the lowest adhesion energy state). Here we present two alternatives to the MMC approach for modeling cellular motion and self-assembly: (1) a kinetic Monte Carlo (KMC), and (2) a cellular particle dynamics (CPD) method. Unlike MMC, both KMC and CPD methods are capable of simulating the dynamics of the cellular system in real time. In the KMC approach a transition rate is associated with possible rearrangements of the cellular system, and the corresponding time evolution is expressed in terms of these rates. In the CPD approach cells are modeled as interacting cellular particles (CPs) and the time evolution of the multicellular system is determined by integrating the equations of motion of all CPs. The KMC and CPD methods are tested and compared by simulating two experimentally well known phenomena: (1) cell-sorting within an aggregate formed by two types of cells with different adhesivities, and (2) fusion of two spherical aggregates of living cells.

preprint2010arXiv

Dynamic heterogeneity in a glass forming fluid: susceptibility, structure factor and correlation length

We investigate the growth of dynamic heterogeneity in a glassy hard-sphere mixture for volume fractions up to and including the mode-coupling transition. We use an 80 000 particle system to test a new procedure to evaluate a dynamic correlation length xi(t): we determine the ensemble independent dynamic susceptibility chi_4(t) and use it to facilitate evaluation of xi(t) from the small wave vector behavior of the four-point structure factor. We analyze relations between the alpha relaxation time tau_alpha, chi_4(tau_alpha), and xi(tau_alpha). We find that mode-coupling like power laws provide a reasonable description of the data over a restricted range of volume fractions, but the power laws' exponents differ from those predicted by the inhomogeneous mode-coupling theory. We find xi(tau_alpha) ~ ln(tau_alpha) over the full range of volume fractions studied, which is consistent with Adams-Gibbs-type relation.

preprint2009arXiv

Anisotropic spatially heterogeneous dynamics on the $α$ and $β$ relaxation time scales studied via a four-point correlation function

We examine the anisotropy of a four-point correlation function $G_4(\vec{k},\vec{r};t)$ and it's associated structure factor $S_4(\vec{k},\vec{q};t)$ calculated using Brownian Dynamics computer simulations of a model glass forming system. These correlation functions measure the spatial correlations of the relaxation of different particles, and we examine the time and temperature dependence of the anisotropy. We find that the anisotropy is strongest at nearest neighbor distances at time scales corresponding to the peak of the non-Gaussian parameter $α_2(t) = 3 < δr^4(t) >/[ 5 < δr^2(t) >^2] - 1$, but is still pronounced around the $α$ relaxation time. We find that the structure factor $S_4(\vec{k},\vec{q};t)$ is anisotropic even for the smallest wave vector accessible in our simulation suggesting that our system (and other systems commonly used in computer simulations) may be too small to extract the $\vec{q} \to 0$ limit of the structure factor. We find that the determination of a dynamic correlation length from $S_4(\vec{k},\vec{q};t)$ is influenced by the anisotropy. We extract an effective anisotropic dynamic correlation length from the small $q$ behavior of $S_4(\vec{k},\vec{q};t)$.

preprint2009arXiv

Diverging length scale of the inhomogeneous mode-coupling theory: a numerical investigation

Biroli et al.'s extension of the standard mode-coupling theory to inhomogeneous equilibrium states [Phys. Rev. Lett. 97, 195701 (2006)] allowed them to identify a characteristic length scale that diverges upon approaching the mode-coupling transition. We present a numerical investigation of this length scale. To this end we derive and numerically solve equations of motion for coefficients in the small q expansion of the dynamic susceptibility $χ_{\mathbf{q}}(\mathbf{k};t)$ that describes the change of the system's dynamics due to an external inhomogeneous potential. We study the dependence of the characteristic length scale on time, wave-vector, and on the distance from the mode-coupling transition. We verify scaling predictions of Biroli et al. In addition, we find that the numerical value of the diverging length scale qualitatively agrees with lengths obtained from four-point correlation functions. We show that the diverging length scale has very weak k dependence, which contrasts with very strong $k$ dependence of the $q\to 0$ limit of the susceptibility, $χ_{\mathbf{q}=0}(\mathbf{k};t)$. Finally, we compare the diverging length obtained from the small q expansion to that resulting from an isotropic approximation applied to the equation of motion for the dynamic susceptibility $χ_{\mathbf{q}}(\mathbf{k};t)$.

preprint2009arXiv

Scaling of the glassy dynamics of soft repulsive particles: a mode-coupling approach

We combine the hyper-netted chain approximation of liquid state theory with the mode-coupling theory of the glass transition to analyze the structure and dynamics of soft spheres interacting via harmonic repulsion. We determine the locus of the fluid-glass dynamic transition in a temperature -- volume fraction phase diagram. The zero-temperature (hard sphere) glass transition influences the dynamics at finite temperatures in its vicinity. This directly implies a form of dynamic scaling for both the average relaxation time and dynamic susceptibilities quantifying dynamic heterogeneity. We discuss several qualitative disagreements between theory and existing simulations at equilibrium. Our theoretical results are, however, very similar to numerical results for the driven athermal dynamics of repulsive spheres, suggesting that `mean-field' mode-coupling approaches might be good starting points to describe these nonequilibrium dynamics.

preprint2006arXiv

Anisotropic spatially heterogeneous dynamics in a model glass-forming binary mixture

We calculated a four-point correlation function G_4(k,r;t) and the corresponding structure factor S_4(k,q;t) for a model glass-forming binary mixture. These functions measure the spatial correlations of the relaxation of different particles. We found that these four-point functions are anisotropic and depend on the angle between vectors k and r (or q). The anisotropy is the strongest for times somewhat longer than the beta relaxation time but it is quite pronounced even for times comparable to the alpha relaxation time, tau_alpha. At the lowest temperatures S_4(k,q;tau_alpha) is strongly anisotropic even for the smallest wavevector q accessible in our simulation.