Biomedical subjects
Publications and source records attributed to Grzegorz Szamel.
We study the four-point dynamic susceptibility chi4(t) obtained from Brownian dynamics computer simulations of the Kob-Andersen Lennard-Jones mixture. We compare the results of the simulations with qualitative predictions of the mode-coupling theory. In addition, we test an estimate of the four-point susceptibility recently proposed by Berthier [Science 310, 1797 (2005)].
Molecular dynamics simulations of high density hard sphere fluids clearly show a breakdown of the Stokes-Einstein equation (SE). This result has been conjectured to be due to the presence of mobile particles, i.e., ones which have the propensity to "hop" distances that are integer multiples of the interparticle distance. We conclusively show that the sedentary particles, i.e., ones complementary to the "hoppers," obey the SE relationship to a good approximation, even though the fluid as a whole violates the SE equation at high densities. These results support the notion that the unusual diffusive behavior of supercooled liquids is dominated by the hopping particles.
We propose a quantitative measure of a time scale on which Fickian diffusion sets in for supercooled liquids, and we use Brownian dynamics computer simulations to determine the temperature dependence of this onset time in a Lennard-Jones binary mixture. The time for the onset of Fickian diffusion ranges between 6.5 and 31 times the relaxation time (the relaxation time is the characteristic relaxation time of the incoherent intermediate scattering function). The onset time increases faster with decreasing temperature than the relaxation time. Mean-squared displacement at the onset time increases with decreasing temperature.
We solved the mode-coupling equations for the Kob-Andersen binary mixture using structure factors calculated from Brownian dynamics simulations of the same system. We found, as was previously observed, that the mode-coupling temperature T(c) inferred from simulations is about two times greater than that predicted by the theory. However, we find that many time-dependent quantities agree reasonably well with the predictions of the mode-coupling theory if they are compared at the same reduced temperature epsilon = (T - T(c))/T(c), and if epsilon is not too small. Specifically, the simulation results for the incoherent intermediate scattering function, the mean square displacement, the relaxation time, and the self-diffusion coefficient agree reasonably well with the predictions of the mode-coupling theory. We find that there are substantial differences for the non-Gaussian parameter. At small reduced temperatures the probabilities of the logarithm of single particle displacements demonstrate that there is hopping-like motion present in the simulations, and this motion is not predicted by the mode-coupling theory. The wave-vector-dependent relaxation time is shown to be qualitatively different from the predictions of the mode-coupling theory for temperatures where hopping-like motion is present.
In a recent publication C. C. Yu and H. M. Carruzzo [Phys. Rev. E 69, 051201 (2004)] determined that the minimum sampling time to calculate the specific heat in a supercooled liquid using energy fluctuations is on the order of 10(3) alpha relaxation times, which is much longer than the sampling time used in most simulations. We demonstrate how the specific heat can be calculated from simulations run to 15 alpha relaxation times, and use statistical arguments to explain the systematic deviation of the specific heat calculated from a simulation of finite length from the true expectation value of the specific heat.
We examine the relaxation of the Kob-Andersen Lennard-Jones binary mixture using Brownian dynamics computer simulations. We find that in accordance with mode-coupling theory the self-diffusion coefficient and the relaxation time show power-law dependence on temperature. However, different mode-coupling temperatures and power laws can be obtained from the simulation data depending on the range of temperatures chosen for the power-law fits. The temperature that is commonly reported as this system's mode-coupling transition temperature, in addition to being obtained from a power law fit, is a crossover temperature at which there is a change in the dynamics from the high-temperature homogeneous, diffusive relaxation to a heterogeneous, hopping-like motion. The hopping-like motion is evident in the probability distributions of the logarithm of single-particle displacements: approaching the commonly reported mode-coupling temperature these distributions start exhibiting two peaks. Notably, the temperature at which the hopping-like motion appears for the smaller particles is slightly higher than that at which the hopping-like motion appears for the larger ones. We define and calculate a new non-Gaussian parameter whose maximum occurs approximately at the time at which the two peaks in the probability distribution of the logarithm of displacements are most evident.
A four-time correlation function was calculated using a computer simulation of a binary Lennard-Jones mixture. The information content of the four-time correlation function is similar to that of four-time correlation functions measured in NMR experiments. The correlation function selects a subensemble and analyzes its dynamics after some waiting time. The lifetime of the subensemble selected by the four-time correlation function is calculated, and compared to the lifetimes of slow subensembles selected using two different definitions of mobility, and to the alpha relaxation time.
Using memory-function formalism we show that in sheared colloidal suspensions the fluctuation-dissipation theorem for self-diffusion, i.e., Einstein's relation between self-diffusion and mobility tensors, is violated and propose a new way to measure this violation in Brownian dynamics simulations. We derive mode-coupling expressions for the tagged particle friction tensor and for an effective, shear-rate dependent temperature.
We study the simplest model of dynamic heterogeneities in glass forming liquids: one-spin facilitated kinetic Ising model introduced by Fredrickson and Andersen [G. H. Fredrickson and H. C. Andersen, Phys. Rev. Lett. 53, 1244 (1984); J. Chem. Phys. 83, 5822 (1985)]. We show that the low-temperature, long-time behavior of the density autocorrelation function predicted by a scaling approach can be obtained from a self-consistent mode-couplinglike approximation.
A new theory for the dynamics of concentrated colloidal suspensions and the colloidal glass transition is proposed. The starting point is the memory function representation of the density correlation function. The memory function can be expressed in terms of a time-dependent pair-density correlation function. An exact, formal equation of motion for this function is derived and a factorization approximation is applied to its evolution operator. In this way a closed set of equations for the density correlation function and the memory function is obtained. The theory predicts an ergodicity breaking transition similar to that predicted by mode-coupling theory, but at a higher density.