Source author record

Walter Kob

Walter Kob 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

40works
9topics
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

40 published item(s)

preprint2023arXiv

From creep to flow: Granular materials under cyclic shear

Granular materials such as sand, powders, and grains are omnipresent in daily life, industrial applications, and earth-science [1]. When unperturbed, they form stable structures that resemble the ones of other amorphous solids like metallic and colloidal glasses [2]. It is commonly conjectured that all these amorphous materials show a universal mechanical response when sheared slowly, i.e., to have an elastic regime, followed by yielding [3]. Here we use X-ray tomography to determine the microscopic dynamics of a cyclically sheared granular system in three dimensions. Independent of the shear amplitude $Γ$, the sample shows a cross-over from creep to diffusive dynamics, indicating that granular materials have no elastic response and always yield, in stark contrast to other glasses. The overlap function [4] reveals that at large $Γ$ yielding is a simple cross-over phenomenon, while for small $Γ$ it shows features of a first order transition with a critical point at $Γ\approx 0.1$ at which one finds a pronounced slowing down and dynamical heterogeneity. Our findings are directly related to the surface roughness of granular particles which induces a micro-corrugation to the potential energy landscape, thus creating relaxation channels that are absent in simple glasses. These processes must be understood for reaching an understanding of the complex relaxation dynamics of granular systems.

preprint2022arXiv

Creating bulk ultrastable glasses by random particle bonding

A recent breakthrough in glass science has been the synthesis of ultrastable glasses via physical vapor deposition techniques. These samples display enhanced thermodynamic, kinetic and mechanical stability, with important implications for fundamental science and technological applications. However, the vapor deposition technique is limited to atomic, polymer and organic glass-formers and is only able to produce thin film samples. Here, we propose a novel approach to generate ultrastable glassy configurations in the bulk, via random particle bonding, and using computer simulations we show that this method does indeed allow for the production of ultrastable glasses. Our technique is in principle applicable to any molecular or soft matter system, such as colloidal particles with tunable bonding interactions, thus opening the way to the design of a large class of ultrastable glasses.

preprint2021arXiv

The critical role of the interaction potential and simulation protocol for the structural and mechanical properties of sodosilicate glasses

We compare the ability of various interaction potentials to predict the structural and mechanical properties of silica and sodium silicate glasses. While most structural quantities show a relatively mild dependence on the potential used, the mechanical properties such as the failure stress and strain as well as the elastic moduli depend very strongly on the potential, once finite size effects have been taken into account. We find that to avoid such finite size effects, samples of at least 75,000 atoms are needed. Finally we probe how the simulation ensemble influences the fracture properties of the glasses and conclude that fracture simulations should be carried out in the constant pressure ensemble.

preprint2020arXiv

Friction-controlled entropy-stability competition in granular systems

Using cyclic shear to drive a two dimensional granular system, we determine the structural characteristics for different inter-particle friction coefficients. These characteristics are the result of a competition between mechanical stability and entropy, with the latter's effect increasing with friction. We show that a parameter-free maximum-entropy argument alone predicts an exponential cell order distribution, with excellent agreement with the experimental observation. We show that friction only tunes the mean cell order and, consequently, the exponential decay rate and the packing fraction. We further show that cells, which can be very large in such systems, are short-lived, implying that our systems are liquid-like rather than glassy.

preprint2020arXiv

Revealing the three-dimensional structure of liquids using four-point correlation functions

Disordered systems like liquids, gels, glasses, or granular materials are not only ubiquitous in daily life and in industrial applications but they are also crucial for the mechanical stability of cells or the transport of chemical and biological agents in living organisms. Despite the importance of these systems, their microscopic structure is understood only on a rudimentary level, thus in stark contrast to the case of gases and crystals. Since scattering experiments and analytical calculations usually give only structural information that is spherically averaged, the three dimensional (3D) structure of disordered systems is basically unknown. Here we introduce a simple method that allows to probe the 3D structure of such systems. Using computer simulations we find that hard-sphere-like liquids have on intermediate and large scales an intricate structural order given by alternating layers with icosahedral and dodecahedral symmetries, while open network liquids like silica have a structural order with tetrahedral symmetry. These results show that liquids have a highly non-trivial 3D structure and that this structural information is encoded in non-standard correlation functions.

preprint2020arXiv

Roughness and scaling properties of oxide glass surfaces at the nanoscale

Using atomistic computer simulations we determine the roughness and topographical features of melt-formed (MS) and fracture surfaces (FS) of oxide glasses. We find that the topography of the MS is described well by the frozen capillary wave theory. The FS are significant rougher than the MS and depend strongly on glass composition. The height-height correlation function for the FS shows an unexpected logarithmic dependence on distance, in contrast to the power-law found in experiments. We thus conclude that on length scales less than 10 nm FS are not self-affine fractals.

preprint2019arXiv

Addition to Structure and dynamics of a polymer-nanoparticle composite: Effect of nanoparticle size and volume fraction

In our previous publication (Ref. 1) we have shown that the data for the normalized diffusion coefficient of the polymers, $D_p/D_{p0}$, falls on a master curve when plotted as a function of $h/λ_d$, where $h$ is the mean interparticle distance and $λ_d$ is a dynamic length scale. In the present note we show that also the normalized diffusion coefficient of the nanoparticles, $D_N/D_{N0}$, collapses on a master curve when plotted as a function of $h/R_h$, where $R_h$ is the hydrodynamic radius of the nanoparticles.

preprint2019arXiv

Characterizing the mesh size of polymer solutions via the pore size distribution

In order to characterize the geometrical mesh size $ξ$, we simulate a solution of coarse-grained polymers with densities ranging from the dilute to the concentrated regime and for different chain lengths. Conventional ways to estimate $ξ$ rely either on scaling assumptions which give $ξ$ only up to an unknown multiplicative factor, or on measurements of the monomer density fluctuation correlation length $ξ_c$. We determine $ξ_c$ from the monomer structure factor and from the radial distribution function, and find that the identification $ξ=ξ_c$ is not justified outside of the semidilute regime. In order to better characterize $ξ$, we compute the pore size distribution (PSD) following two different definitions, one by Torquato et al. (Ref.1) and one by Gubbins et al. (Ref.2). We show that the mean values of the two distributions, $\langle r \rangle_T$ and $\langle r \rangle_G$, both display the behavior predicted for $ξ$ by scaling theory, and argue that $ξ$ can be identified with either one of these quantities. This identification allows to interpret the PSD as the distribution of mesh sizes, a quantity which conventional methods cannot access. Finally, we show that it is possible to map a polymer solution on a system of hard or overlapping spheres, for which Torquato's PSD can be computed analytically and reproduces accurately the PSD of the solution. We give an expression that allows $\langle r \rangle_T$ to be estimated with great accuracy in the semidilute regime by knowing only the radius of gyration and the density of the polymers.

preprint2019arXiv

New interaction potentials for borate glasses with mixed network formers

We adapt and apply a recently developed optimization scheme used to obtain effective potentials for aluminosilicate glasses to include the network former boron into the interaction parameter set. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by ab initio molecular dynamics simulations, and density, coordination and elastic modulus of glass at room temperature from experiments. The new interaction potentials are shown to reproduce reliably the structure, coordination and mechanical properties over a wide range of compositions for binary alkali borates. Furthermore, the transferability of these new interaction parameters allows mixing to reliably reproduce properties of various boroaluminate and borosilicate glasses.

preprint2016arXiv

First-principles simulations of glass-formers

In this article we review results of computer simulation of glasses carried out using first principles approaches, notably density functional theory. We start with a brief introduction to this method and compare the pros and cons of this approach with the ones of simulations with classical potentials. This is followed by a discussion of simulation results of various glass-forming systems that have been obtained via ab initio simulations and that demonstrate the usefulness of this approach to understand the properties of glasses on the microscopic level.

preprint2015arXiv

Realistic Tunneling States for the Magnetic Effects in Non-Metallic Real Glasses

The discovery of magnetic and compositional effects in the low temperature properties of multi-component glasses has prompted the need to extend the standard two-level systems (2LSs) tunneling model. A possible extension \cite{Jug2004} assumes that a subset of tunneling quasi-particles is moving in a three-welled potential (TWP) associated with the ubiquitous inhomogeneities of the disordered atomic structure of the glass. We show that within an alternative, cellular description of the intermediate-range atomic structure of glasses the tunneling TWP can be fully justified. We then review how the experimentally discovered magnetic effects can be explained within the approach where only localized atomistic tunneling 2LSs and quasi-particles tunneling in TWPs are allowed. We discuss the origin of the magnetic effects in the heat capacity, dielectric constant (real and imaginary parts), polarization echo and SQUID magnetization in several glassy systems. We conclude by commenting on a strategy to reveal the mentioned tunneling states (2LSs and TWPs) by means of atomistic computer simulations and discuss the microscopic nature of the tunneling states in the context of the potential energy landscape of glass-forming systems.

preprint2015arXiv

Relaxation dynamics in a transient network fluid with competing gel and glass phases

We use computer simulations to study the relaxation dynamics of a model for oil-in-water microemulsion droplets linked with telechelic polymers. This system exhibits both gel and glass phases and we show that the competition between these two arrest mechanisms can result in a complex, three-step decay of the time correlation functions, controlled by two different localization lengthscales. For certain combinations of the parameters, this competition gives rise to an anomalous logarithmic decay of the correlation functions and a subdiffusive particle motion, which can be understood as a simple crossover effect between the two relaxation processes. We establish a simple criterion for this logarithmic decay to be observed. We also find a further logarithmically slow relaxation related to the relaxation of floppy clusters of particles in a crowded environment, in agreement with recent findings in other models for dense chemical gels. Finally, we characterize how the competition of gel and glass arrest mechanisms affects the dynamical heterogeneities and show that for certain combination of parameters these heterogeneities can be unusually large. By measuring the four-point dynamical susceptibility, we probe the cooperativity of the motion and find that with increasing coupling this cooperativity shows a maximum before it decreases again, indicating the change in the nature of the relaxation dynamics. Our results suggest that compressing gels to large densities produces novel arrested phases that have a new and complex dynamics.

preprint2014arXiv

Crossovers in the dynamics of supercooled liquids probed by an amorphous wall

We study the relaxation dynamics of a binary Lennard-Jones liquid in the presence of an amorphous wall generated from equilibrium particle configurations. In qualitative agreement with the results presented in Nature Phys. {\bf 8}, 164 (2012) for a liquid of harmonic spheres, we find that our binary mixture shows a saturation of the dynamical length scale close to the mode-coupling temperature $T_c$. Furthermore we show that, due to the broken symmetry imposed by the wall, signatures of an additional change in dynamics become apparent at a temperature well above $T_c$. We provide evidence that this modification in the relaxation dynamics occurs at a recently proposed dynamical crossover temperature $T_s > T_c$, which is related to the breakdown of the Stokes-Einstein relation. We find that this dynamical crossover at $T_s$ is also observed for a system of harmonic spheres as well as a WCA liquid, showing that it may be a general feature of glass-forming systems.

preprint2014arXiv

Non-linear dynamic response of glass-forming liquids to random pinning

We use large scale computer simulations of a glass-forming liquid in which a fraction c of the particles has been permanently pinned. We find that the relaxation dynamics shows an exponential dependence on c. This result can be rationalized by means of a simple theoretical Ansatz and we discuss its implication for thermodynamic theories for the glass-transition. For intermediate and low temperatures we find that the slowing down of the dynamics due to the pinning saturates and that the cooperativity decreases with increasing c, results which indicate that in glass-forming liquids there is a dynamic crossover at which the shape of the relaxing entities changes.

preprint2013arXiv

Gelling by Heating

We introduce a simple model, a binary mixture of patchy particles, which has been designed to form a gel upon heating. Due to the specific nature of the particle interactions, notably the number and geometry of the patches as well as their interaction energies, the system is a fluid both at high and at low temperatures, whereas at intermediate temperatures the system forms a solid-like disordered open network structure, i.e. a gel. Using molecular dynamics we investigate the static and dynamic properties of this system.

preprint2013arXiv

Intermittent dynamics and logarithmic domain growth during the spinodal decomposition of a glass-forming liquid

We use large-scale molecular dynamics simulations of a simple glass-forming system to investigate how its liquid-gas phase separation kinetics depends on temperature. A shallow quench leads to a fully demixed liquid-gas system whereas a deep quench makes the dense phase undergo a glass transition and become an amorphous solid. This glass has a gel-like bicontinuous structure that evolves very slowly with time and becomes fully arrested in the limit where thermal fluctuations become negligible. We show that the phase separation kinetics changes qualitatively with temperature, the microscopic dynamics evolving from a surface tension-driven diffusive motion at high temperature to a strongly intermittent, heterogeneous and thermally activated dynamics at low temperature, with a logarithmically slow growth of the typical domain size. These results shed light on recent experimental observations of various porous materials produced by arrested spinodal decomposition, such as nonequilibrium colloidal gels and bicontinuous polymeric structures, and they elucidate the microscopic mechanisms underlying a specific class of viscoelastic phase separation.

preprint2013arXiv

Phase Diagram of a Reentrant Gel of Patchy Particles

We study the phase diagram of a binary mixture of patchy particles which has been designed to form a reversible gel. For this we perform Monte Carlo and molecular dynamics simulations to investigate the thermodynamics of such a system and compare our numerical results with predictions based on the analytical parameter-free Wertheim theory. We explore a wide range of the temperature-density-composition space that defines the three-dimensional phase diagram of the system. As a result, we delimit the region of thermodynamic stability of the fluid. We find that for a large region of the phase diagram the Wertheim theory is able to give a quantitative description of the system. For higher densities our simulations show that the system is crystallizing into a BCC structure. Finally we study the relaxation dynamics of the system by means of the density and temperature dependences of the diffusion coefficient. We show that there exists a density range where the system passes reversibly from a gel to a fluid upon both heating and cooling, encountering neither demixing nor phase separation.

preprint2013arXiv

Probing a liquid to glass transition in equilibrium

We use computer simulations to investigate the static properties of a simple glass-forming fluid in which the positions of a finite fraction of the particles has been frozen in. By probing the equilibrium distribution of the overlap between independent configurations of the liquid, we find strong evidence that this random pinning induces a glass transition. At low temperatures, our numerical findings are consistent with the existence of a random first order phase transition rounded by finite size effects.

preprint2012arXiv

Finite size effects in the dynamics of glass-forming liquids

We present a comprehensive theoretical study of finite size effects in the relaxation dynamics of glass-forming liquids. Our analysis is motivated by recent theoretical progress regarding the understanding of relevant correlation length scales in liquids approaching the glass transition. We obtain predictions both from general theoretical arguments and from a variety of specific perspectives: mode-coupling theory, kinetically constrained and defect models, and random first order transition theory. In the latter approach, we predict in particular a non-monotonic evolution of finite size effects across the mode-coupling crossover due to the competition between mode-coupling and activated relaxation. We study the role of competing relaxation mechanisms in giving rise to non-monotonic finite size effects by devising a kinetically constrained model where the proximity to the mode-coupling singularity can be continuously tuned by changing the lattice topology. We use our theoretical findings to interpret the results of extensive molecular dynamics studies of four model liquids with distinct structures and kinetic fragilities. While the less fragile model only displays modest finite size effects, we find a more significant size dependence evolving with temperature for more fragile models, such as Lennard-Jones particles and soft spheres. Finally, for a binary mixture of harmonic spheres we observe the predicted non-monotonic temperature evolution of finite size effects near the fitted mode-coupling singularity, suggesting that the crossover from mode-coupling to activated dynamics is more pronounced for this model. Finally, we discuss the close connection between our results and the recent report of a non-monotonic temperature evolution of a dynamic length scale near the mode-coupling crossover in harmonic spheres.

preprint2012arXiv

Spatial correlations in glass-forming liquids across the mode-coupling crossover

We discuss a novel approach, the point-to-set correlation functions, that allows to determine relevant static and dynamic length scales in glass-forming liquids. We find that static length scales increase monotonically when the temperature is lowered, whereas the measured dynamic length scale shows a maximum around the critical temperature of mode-coupling theory. We show that a similar non-monotonicity is found in the temperature evolution of certain finite size effects in the relaxation dynamics. These two independent sets of results demonstrate the existence of a change in the transport mechanism when the glass-former is cooled from moderately to deeply supercooled states across the mode-coupling crossover and clarify the status of the theoretical calculations done at the mean field level.

preprint2011arXiv

Computer Simulation Study of the Phase Behavior and Structural Relaxation in a Gel-Former Modeled by Three Body Interactions

We report a computer simulation study of a model gel-former obtained by modifying the three-body interactions of the Stillinger-Weber potential for silicon. This modification reduces the average coordination number and consequently shifts the liquid-gas phase coexistence curve to low densities, thus facilitating the formation of gels without phase separation. At low temperatures and densities, the structure of the system is characterized by the presence of long linear chains interconnected by a small number of three coordinated junctions at random locations. At small wave-vectors the static structure factor shows a non-monotonic dependence on temperature, a behavior which is due to the competition between the percolation transition of the particles and the stiffening of the formed chains. We compare in detail the relaxation dynamics of the system as obtained from molecular dynamics with the one obtained from Monte Carlo dynamics. We find that the bond correlation function displays stretched exponential behavior at moderately low temperatures and densities, but exponential relaxation at low temperatures. The bond lifetime shows an Arrhenius behavior, independent of the microscopic dynamics. For the molecular dynamics at low temperatures, the mean squared displacement and the (coherent and incoherent) intermediate scattering function display at intermediate times a dynamics with ballistic character and we show that this leads to compressed exponential relaxation. For the Monte Carlo dynamics we find always an exponential or stretched exponential relaxation. Thus we conclude that the compressed exponential relaxation observed in experiments is due to the out-of-equilibrium dynamics.

preprint2011arXiv

Influence of the glass transition on the liquid-gas spinodal decomposition

We use large-scale molecular dynamics simulations to study the kinetics of the liquid-gas phase separation if the temperature is lowered across the glass transition of the dense phase. We observe a gradual change from phase separated systems at high temperatures to nonequilibrium, gel-like structures that evolve very slowly at low temperatures. The microscopic mechanisms responsible for the coarsening strongly depend on temperature, and change from diffusive motion at high temperature to a strongly intermittent, heterogeneous and thermally activated dynamics at low temperature, leading to logarithmically slow growth of the typical domain size.

preprint2011arXiv

Non-monotonic temperature evolution of dynamic correlations in glass-forming liquids

The viscosity of glass-forming liquids increases by many orders of magnitude if their temperature is lowered by a mere factor of 2-3 [1,2]. Recent studies suggest that this widespread phenomenon is accompanied by spatially heterogeneous dynamics [3,4], and a growing dynamic correlation length quantifying the extent of correlated particle motion [5-7]. Here we use a novel numerical method to detect and quantify spatial correlations which reveal a surprising non-monotonic temperature evolution of spatial dynamical correlations, accompanied by a second length scale that grows monotonically and has a very different nature. Our results directly unveil a dramatic qualitative change in atomic motions near the mode-coupling crossover temperature [8] which involves no fitting or indirect theoretical interpretation. Our results impose severe new constraints on the theoretical description of the glass transition, and open several research perspectives, in particular for experiments, to confirm and quantify our observations in real materials.

preprint2011arXiv

Predicting the Structure of Alloys using Genetic Algorithms

We discuss a novel genetic algorithm that can be used to find global minima on the potential energy surface of disordered ceramics and alloys using a real-space symmetry adapted crossover. Due to a high number of symmetrically equivalent solutions of many alloys a conventional genetic algorithms using reasonable population sizes are unable to locate the global minima for even the smallest systems. We demonstrate the superior performance of the use of symmetry adapted crossover by the comparison of that of a conventional GA for finding global minima of two binary Ising-type alloys that either order or phase separate at low temperature. Comparison of different representations and crossover operations show that the use of real-space crossover outperforms crossover operators working on binary representations by several orders of magnitude.

preprint2011arXiv

Static point-to-set correlations in glass-forming liquids

We analyze static point-to-set correlations in glass-forming liquids. The generic idea is to freeze the position of a set of particles in an equilibrium configuration and to perform sampling in the presence of this additional constraint. Qualitatively different geometries for the confining set of particles are considered and a detailed comparison of resulting static and dynamic correlation functions is performed. Our results reveal the existence of static spatial correlations not detected by conventional two-body correlators, which appear to be decoupled from, and shorter-ranged than, dynamical length scales characterizing dynamic heterogeneity. We find that the dynamics slows down dramatically under confinement, which suggests new ways to investigate the glass transition. Our results indicate that the geometry in which particles are randomly pinned is the best candidate to study static correlations.

preprint2011arXiv

The vibrational density of states of a disordered gel model

We investigate the vibrational density of states (vDOS) in harmonic approximation of a binary mixture of colloidal patchy particles with two and three patches for different relative compositions $x_2$. At low temperature, this system forms a thermo-reversible gel, i.e. a fully bonded network of chains of two-patches particles, in which the branching points are provided by three-patches particles. For all the compositions we find in the vDOS a pronounced peak at low frequency whose height grows on increasing the fraction of two-functional particles, or equivalently with the average length of the chains. To identify the various spectral features, we compare the vDOS of the whole system with the one of small representative structures of the network and with the vDOS of a long linear chain of two-patches particles and we find that these structures are indeed able to rationalize the various peaks in the vDOS of the full system. At large $x_2$ the vDOS of the gel and of the long chain show remarkable similarities. Analyzing the dispersion relations and the spectrum of the linear chain we show that the excess of low frequency modes, the analog of the boson peak in glassy disordered systems, arises from the strong coupling between rotations and translations.

preprint2010arXiv

A microscopic model for colloidal gels with directional effective interactions: Network induced glassy dynamics

By means of molecular dynamics, we study the structure and the dynamics of a microscopic model for colloidal gels at low volume fractions. The presence of directional interactions leads to the formation of a persistent interconnected network at temperatures where phase separation does not occur. We find that large scale spatial correlations strongly depend on the volume fraction and characterize the formation of the persistent network. We observe a pre-peak in the static structure factor and relate it to the network structure. The slow dynamics at gelation is characterized by the coexistence of fast collective motion of the mobile parts of the network structure (chains) with large scale rearrangements producing stretched exponential relaxations. We show that, once the network is sufficiently persistent, it induces slow, cooperative processes related to the network nodes. We suggest that this peculiar glassy dynamics is a hallmark of the physics of colloidal gels at low volume fractions.

preprint2009arXiv

Structural Relaxation of a Gel Modeled by Three Body Interactions

We report a molecular dynamics simulation study of a model gel whose interaction potential is obtained by modifying the three body Stillinger-Weber model potential for silicon. The modification reduces the average coordination number, and suppresses the liquid-gas phase coexistence curve. The low density, low temperature equilibrium gel that can thus form exhibits interesting dynamical behavior, including compressed exponential relaxation of density correlations. We show that motion responsible for such relaxation has ballistic character, and arises from the motion of chain segments in the gel without the restructuring of the gel network.

preprint2009arXiv

When gel and glass meet: A mechanism for multistep relaxation

We use computer simulations to study the dynamics of a physical gel at high densities where gelation and the glass transition interfere. We report and provide detailed physical understanding of complex relaxation patterns for time correlation functions which generically decay in a three-step process. For certain combinations of parameters we find logarithmic decays of the correlators and subdiffusive particle motion.

preprint2008arXiv

On the relaxation dynamics of glass-forming systems: Insights from computer simulations

We discuss the relaxation dynamics of a simple lattice gas model for glass-forming systems and show that with increasing density of particles this dynamics slows down very quickly. By monitoring the trajectory of tagged particles we find that their motion is very heterogeneous in space and time, leading to regions in space in which there is a fast dynamics and others in which it is slow. We determine how the geometric properties of these quickly relaxing regions depend on density and time. Motivated by this heterogeneous hopping dynamics, we use a simple model, a variant of a continuous time random walk, to characterize the relaxation dynamics. In particular we find from this model that for large displacements the self part of the van Hove function shows an exponential tail, in agreement with recent findings from experiments and simulations of glass-forming systems.

preprint2007arXiv

Heterogeneous diffusion in a reversible gel

We introduce a microscopically realistic model of a physical gel and use computer simulations to study its static and dynamic properties at thermal equilibrium. The phase diagram comprises a sol phase, a coexistence region ending at a critical point, a gelation line determined by geometric percolation, and an equilibrium gel phase unrelated to phase separation. The global structure of the gel is homogeneous, but the stress is supported by a fractal network. The gel dynamics is highly heterogeneous and we propose a theoretical model to quantitatively describe dynamic heterogeneity in gels. We elucidate several differences between the dynamics of gels and that of glass-formers.

preprint2007arXiv

Structure and dynamics of l-Ge: Neutron scattering experiments and ab initio molecular dynamics simulations

We report the first measurements of the dynamics of liquid germanium (l-Ge) by quasi-elastic neutron scattering on time-of-flight and triple-axis spectrometers. These results are compared with simulation data of the structure and dynamics of l-Ge which have been obtained with ab initio density functional theory methods. The simulations accurately reproduce previous results from elastic and inelastic scattering experiments, as well as the q-dependence of the width of the quasi-elastic signal of the new experimental data. In order to understand some special features of the structure of the liquid we have also simulated amorphous Ge. Overall we find that the atomistic model represents accurately the average structure of real l-Ge as well as the time dependent structural fluctuations. The new quasi-elastic neutron scattering data allows us to investigate to what extent simple theoretical models can be used to describe diffusion in l-Ge.

preprint2006arXiv

Structure and Relaxation Dynamics of a Colloidal Gel

Using molecular dynamics computer simulations we investigate the structural and dynamical properties of a simple model for a colloidal gel at low volume fraction. We find that at low T the system is forming an open percolating cluster, without any sign of a phase separation. The nature of the relaxation dynamics depends strongly on the length scale/wave-vector considered and can be directly related to the geometrical properties of the spanning cluster.

preprint2006arXiv

The Monte-Carlo dynamics of a binary Lennard-Jones glass-forming mixture

We use a standard Monte-Carlo algorithm to study the slow dynamics of a binary Lennard-Jones glass-forming mixture at low temperature. We find that Monte-Carlo is by far the most efficient way to simulate a stochastic dynamics since relaxation is about 10 times faster than in Brownian Dynamics and about 30 times faster than in Stochastic Dynamics. Moreover, the average dynamical behaviour of the system is in quantitative agreement with the one obtained using Newtonian dynamics, apart at very short times where thermal vibrations are suppressed. We show, however, that dynamic fluctuations quantified by four-point dynamic susceptibilities do retain a dependence on the microscopic dynamics, as recently predicted theoretically.

preprint1997arXiv

Molecular Dynamics Computer Simulation of the Dynamics of Supercooled Silica

We present the results of a large scale computer simulation of supercooled silica. We find that at high temperatures the diffusion constants show a non-Arrhenius temperature dependence whereas at low temperature this dependence is also compatible with an Arrhenius law. We demonstrate that at low temperatures the intermediate scattering function shows a two-step relaxation behavior and that it obeys the time temperature superposition principle. We also discuss the wave-vector dependence of the nonergodicity parameter and the time and temperature dependence of the non-Gaussian parameter.

preprint1996arXiv

A quantitative test of the mode-coupling theory of the ideal glass transition for a binary Lennard-Jones system

Using a molecular dynamics computer simulation we determine the temperature dependence of the partial structure factors for a binary Lennard-Jones system. These structure factors are used as input data to solve numerically the wave-vector dependent mode-coupling equations in the long time limit. Using the so determined solutions, we compare the predictions of mode-coupling theory (MCT) with the results of a previously done molecular dynamics computer simulation [Phys. Rev. E 51, 4626 (1995), ibid. 52, 4134 (1995)]. From this comparison we conclude that MCT gives a fair estimate of the critical coupling constant, a good estimate of the exponent parameter, predicts the wave-vector dependence of the various nonergodicity parameters very well, except for very large wave-vectors, and gives also a very good description of the space dependence of the various critical amplitudes. In an attempt to correct for some of the remaining discrepancies between the theory and the results of the simulation, we investigate two small (ad hoc) modifications of the theory. We find that one modification gives a worse agreement between theory and simulation, whereas the second one leads to an improved agreement.

preprint1995arXiv

Testing Mode-Coupling Theory for a Supercooled Binary Lennard-Jones Mixture I: The van Hove Correlation Function

We report the results of a large scale computer simulation of a binary supercooled Lennard-Jones liquid. We find that at low temperatures the curves for the mean squared displacement of a tagged particle for different temperatures fall onto a master curve when they are plotted versus rescaled time $tD(T)$, where $D(T)$ is the diffusion constant. The time range for which these curves follow the master curve is identified with the $α$-relaxation regime of mode-coupling theory (MCT). This master curve is fitted well by a functional form suggested by MCT. In accordance with idealized MCT, $D(T)$ shows a power-law behavior at low temperatures. The critical temperature of this power-law is the same for both types of particle and also the critical exponents are very similar. However, contrary to a prediction of MCT, these exponents are not equal to the ones determined previously for the divergence of the relaxation times of the intermediate scattering function [Phys. Rev. Lett. {\bf 73}, 1376 (1994)]. At low temperatures the van Hove correlation function (self as well as distinct part) shows hardly any sign of relaxation in a time interval that extends over about three decades in time. This time interval can be interpreted as the $β$-relaxation regime of MCT. From the investigation of these correlation functions we conclude that hopping processes are not important on the time scale of the $β$-relaxation for this system and for the temperature range investigated. We test whether the factorization property predicted by MCT holds and find that this is indeed the case for all correlation functions investigated. The distance dependence of the critical amplitudes are in qualitative accordance with the ones predicted by MCT for some other mixtures. The non-gaussian parameter for the self part of the van