Source author record

Grzegorz Szamel

Grzegorz Szamel 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

33works
4topics
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

33 published item(s)

preprint2022arXiv

An alternative, dynamic density functional-like theory for time-dependent density fluctuations in glass-forming fluids

We propose an alternative theory for the relaxation of density fluctuations in glass-forming fluids. We derive an equation of motion for the density correlation function which is local in time and is similar in spirit to the equation of motion for the average non-uniform density profile derived within the dynamic density functional theory. We identify the Franz-Parisi free energy functional as the non-equilibrium free energy for the evolution of the density correlation function. An appearance of a local minimum of this functional leads to a dynamic arrest. Thus, the ergodicity breaking transition predicted by our theory coincides with the dynamic transition of the static approach based on the same non-equilibrium free energy functional.

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.

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

A theory for the dynamics of dense systems of athermal self-propelled particles

We present a derivation of a recently proposed theory for the time dependence of density fluctuations in stationary states of strongly interacting, athermal, self-propelled particles. The derivation consists of two steps. First, we start from the equation of motion for the joint distribution of particles' positions and self-propulsions and we integrate out the self-propulsions. In this way we derive an approximate, many-particle evolution equation for the probability distribution of the particles' positions. Second, we use this evolution equation to describe the time dependence of steady-state density correlations. We derive a memory function representation of the density correlation function and then we use a factorization approximation to obtain an approximate expression for the memory function. In the final equation of motion for the density correlation function the non-equilibrium character of the active system manifests itself through the presence of a new steady-state correlation function that quantifies spatial correlations of the velocities of the particles. This correlation function enters into the frequency term, and thus it describes the dependence of the short-time dynamics on the properties of the self-propulsions. More importantly, the correlation function of particles' velocities enters into the vertex of the memory function and through the vertex it modifies the long-time glassy dynamics.

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

On the coarse-grained density and compressibility of a non-ideal crystal

The isothermal compressibility of a general crystal is derived and discussed within (classical) density functional theory. Starting from the microscopic particle density, we carefully coarse grain to obtain the thermodynamic compressibility and the long wavelength limits of the correlation functions of elasticity theory. We explicitly show that the isothermal compressibility is not the wavevector to zero limit of the (total) density correlation function. The latter also cannot be obtained from the static structure factor measured in a scattering experiment.

preprint2015arXiv

On the coarse-grained density and compressibility of non-ideal crystals: general theory and an application to cluster crystals

The isothermal compressibility of a general crystal is analyzed within classical density functional theory. Our approach can be used for homogeneous and unstrained crystals containing an arbitrarily high density of local defects. We start by coarse-graining the microscopic particle density and then obtain the long wavelength limits of the correlation functions of elasticity theory and the thermodynamic derivatives. We explicitly show that the long wavelength limit of the microscopic density correlation function differs from the isothermal compressibility. It also cannot be obtained from the static structure factor measured in a scattering experiment. We apply our theory to crystals consisting of soft particles which can multiply occupy lattice sites ('cluster crystals'). The multiple occupancy results in a strong local disorder over an extended range of temperatures. We determine the cluster crystals' isothermal compressibility, the fluctuations of the lattice occupation numbers and their correlation functions, and the dispersion relations. We also discuss their low-temperature phase diagram.

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

A self-propelled particle in an external potential: is there an effective temperature?

We study a stationary state of a single self-propelled, athermal particle in linear and quadratic external potentials. The self-propulsion is modeled as a fluctuating force evolving according to the Ornstein-Uhlenbeck process, independently of the state of the particle. Without an external potential, in the long time limit, the self-propelled particle moving in a viscous medium performs diffusive motion, which allows one to identify an effective temperature. We show that in the presence of a linear external potential the stationary state distribution has an exponential form with the sedimentation length determined by the effective temperature of the free self-propelled particle. In the presence of a quadratic external potential the stationary state distribution has a Gaussian form. However, in general, this distribution is not determined by the effective temperature of the free self-propelled particle.

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.

preprint2012arXiv

Mode-coupling theory and beyond: a diagrammatic approach

For almost thirty years, mode-coupling theory has been the most widely discussed and used but also the most controversial theory of the glass transition. In this paper we briefly review the reasons for both its popularity and its controversy. We emphasize the need for the development of approaches that would be able to evaluate corrections to and extensions of the existing (standard) mode-coupling theory. Next, we review our diagrammatic formulation of the dynamics of interacting Brownian particles. We show that within this approach the standard mode-coupling theory can be derived in a very simple way. Finally, we use our diagrammatic approach to calculate two corrections to the mode-coupling theory's expression for the so-called irreducible memory function. These corrections involve re-summations of well defined classes of non-mode-coupling diagrams.

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

Effective temperatures of a driven, strongly anisotropic Brownian system

We use Brownian Dynamics computer simulations of a moderately dense colloidal system undergoing steady shear flow to investigate the uniqueness of the so-called effective temperature. We compare effective temperatures calculated from the fluctuation-dissipation ratios and from the linear response to a static, long wavelength, external perturbation along two directions: the shear gradient direction and the vorticity direction. At high shear rates, when the system is strongly anisotropic, the fluctuation-dissipation ratio derived effective temperatures are approximately wave-vector independent, but the temperatures along the gradient direction are somewhat higher than those along the vorticity direction. The temperatures derived from the static linear response show the same dependence on the direction as those derived from the fluctuation-dissipation ratio. However, the former and the latter temperatures are different. Our results suggest that the presently used formulae for effective temperatures may not be applicable for strongly anisotropic, driven systems.

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

Tagged particle in a sheared suspension: effective temperature determines density distribution in a slowly varying external potential beyond linear response

We consider a sheared colloidal suspension under the influence of an external potential that varies slowly in space in the plane perpendicular to the flow and acts on one selected (tagged) particle of the suspension. Using a Chapman-Enskog type expansion we derive a steady state equation for the tagged particle density distribution. We show that for potentials varying along one direction only, the tagged particle distribution is the same as the equilibrium distribution with the temperature equal to the effective temperature obtained from the violation of the Einstein relation between the self-diffusion and tagged particle mobility coefficients. We thus prove the usefulness of this effective temperature for the description of the tagged particle behavior beyond the realm of linear response. We illustrate our theoretical predictions with Brownian dynamics computer simulations.

preprint2010arXiv

Dynamic glass transition: bridging the gap between mode-coupling theory and the replica approach

We clarify the relation between the ergodicity breaking transition predicted by mode-coupling theory and the so-called dynamic transition predicted by the static replica approach. Following Franz and Parisi [Phys. Rev. Lett. 79, 2486 (1997)], we consider a system of particles in a metastable state characterized by non-trivial correlations with a quenched configuration. We show that the assumption that in a metastable state particle currents vanish leads to an expression for the replica off-diagonal direct correlation function in terms of a replica off-diagonal static four-point correlation function. A factorization approximation for this function results in an approximate closure for the replica off-diagonal direct correlation function. The replica off-diagonal Ornstein-Zernicke equation combined with this closure coincides with the equation for the non-ergodicity parameter derived using the mode-coupling theory.

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.

preprint2010arXiv

Increasing the density melts ultrasoft colloidal glasses

We use theory and simulations to investigate the existence of amorphous glassy states in ultrasoft colloids. We combine the hyper-netted chain approximation with mode-coupling theory to study the dynamic phase diagram of soft repulsive spheres interacting with a Hertzian potential, focusing on low temperatures and large densities. At constant temperature, we find that an amorphous glassy state is entered upon compression, as in colloidal hard spheres, but the glass unexpectedly melts when density increases further. We attribute this re-entrant fluid-glass transition to particle softness, and correlate this behaviour to previously reported anomalies in soft systems, thus emphasizing its generality. The predicted fluid-glass-fluid sequence is confirmed numerically.

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.