Hopping in a Supercooled Lennard-Jones Liquid: Metabasins, Waiting Time Distribution, and Diffusion

We investigate the jump motion among potential energy minima of a Lennard-Jones model glass former by extensive computer simulation. From the time series of minima energies, it becomes clear that the energy landscape is organized in superstructures, called metabasins. We show that diffusion can be pictured as a random walk among metabasins, and that the whole temperature dependence resides in the distribution of waiting times. The waiting time distribution exhibits algebraic decays: $\tau^{-1/2}$ for very short times and $\tau^{-\alpha}$ for longer times, where $\alpha\approx2$ near $T_c$. We demonstrate that solely the waiting times in the very stable basins account for the temperature dependence of the diffusion constant.Comment: to be published in Phys. Rev.

What does the potential energy landscape tell us about the dynamics of supercooled liquids and glasses?

For a model glass-former we demonstrate via computer simulations how macroscopic dynamic quantities can be inferred from a PEL analysis. The essential step is to consider whole superstructures of many PEL minima, called metabasins, rather than single minima. We show that two types of metabasins exist: some allowing for quasi-free motion on the PEL (liquid-like), the others acting as traps (solid-like). The activated, multi-step escapes from the latter metabasins are found to dictate the slowing down of dynamics upon cooling over a much broader temperature range than is currently assumed

Finite-Size Effects in a Supercooled Liquid

We study the influence of the system size on various static and dynamic properties of a supercooled binary Lennard-Jones liquid via computer simulations. In this way, we demonstrate that the treatment of systems as small as N=65 particles yields relevant results for the understanding of bulk properties. Especially, we find that a system of N=130 particles behaves basically as two non-interacting systems of half the size.Comment: Proceedings of the III Workshop on Non Equilibrium Phenomena in Supercooled Fluids, Glasses and Amorphous Materials, Sep 2002, Pis

Local Properties of the Potential Energy Landscape of a Model Glass: Understanding the Low Temperature Anomalies

Though the existence of two-level systems (TLS) is widely accepted to explain low temperature anomalies in the sound absorption, heat capacity, thermal conductivity and other quantities, an exact description of their microscopic nature is still lacking. We performed computer simulations for a binary Lennard-Jones system, using a newly developed algorithm to locate double-well potentials (DWP) and thus two-level systems on a systematic basis. We show that the intrinsic limitations of computer simulations like finite time and finite size problems do not hamper this analysis. We discuss how the DWP are embedded in the total potential energy landscape. It turns out that most DWP are connected to the dynamics of the smaller particles and that these DWP are rather localized. However, DWP related to the larger particles are more collective

Backward correlations and dynamic heterogeneities: a computer study of ion dynamics

We analyse the correlated back and forth dynamics and dynamic heterogeneities, i.e. the presence of fast and slow ions, for a lithium metasilicate system via computer simulations. For this purpose we define, in analogy to previous work in the field of glass transition, appropriate three-time correlation functions. They contain information about the dynamics during two successive time intervals. First we apply them to simple model systems in order to clarify their information content. Afterwards we use this formalism to analyse the lithium trajectories. A strong back-dragging effect is observed, which also fulfills the time-temperature superposition principle. Furthermore, it turns out that the back-dragging effect is long-ranged and exceeds the nearest neighbor position. In contrast, the strength of the dynamic heterogeneities does not fulfill the time-temperature superposition principle. The lower the temperature, the stronger the mobility difference between fast and slow ions. The results are then compared with the simple model systems considered here as well as with some lattice models of ion dynamics.Comment: 12 pages, 10 figure

Multiple-Point and Multiple-Time Correlations Functions in a Hard-Sphere Fluid

A recent mode coupling theory of higher-order correlation functions is tested on a simple hard-sphere fluid system at intermediate densities. Multi-point and multi-time correlation functions of the densities of conserved variables are calculated in the hydrodynamic limit and compared to results obtained from event-based molecular dynamics simulations. It is demonstrated that the mode coupling theory results are in excellent agreement with the simulation results provided that dissipative couplings are included in the vertices appearing in the theory. In contrast, simplified mode coupling theories in which the densities obey Gaussian statistics neglect important contributions to both the multi-point and multi-time correlation functions on all time scales.Comment: Second one in a sequence of two (in the first, the formalism was developed). 12 pages REVTeX. 5 figures (eps). Submitted to Phys.Rev.

Simple Lattice-Models of Ion Conduction: Counter Ion Model vs. Random Energy Model

The role of Coulomb interaction between the mobile particles in ionic conductors is still under debate. To clarify this aspect we perform Monte Carlo simulations on two simple lattice models (Counter Ion Model and Random Energy Model) which contain Coulomb interaction between the positively charged mobile particles, moving on a static disordered energy landscape. We find that the nature of static disorder plays an important role if one wishes to explore the impact of Coulomb interaction on the microscopic dynamics. This Coulomb type interaction impedes the dynamics in the Random Energy Model, but enhances dynamics in the Counter Ion Model in the relevant parameter range.Comment: To be published in Phys. Rev.

Complex lithium ion dynamics in simulated LiPO3 glass studied by means of multi-time correlation functions

Molecular dynamics simulations are performed to study the lithium jumps in LiPO3 glass. In particular, we calculate higher-order correlation functions that probe the positions of single lithium ions at several times. Three-time correlation functions show that the non-exponential relaxation of the lithium ions results from both correlated back-and-forth jumps and the existence of dynamical heterogeneities, i.e., the presence of a broad distribution of jump rates. A quantitative analysis yields that the contribution of the dynamical heterogeneities to the non-exponential depopulation of the lithium sites increases upon cooling. Further, correlated back-and-forth jumps between neighboring sites are observed for the fast ions of the distribution, but not for the slow ions and, hence, the back-jump probability depends on the dynamical state. Four-time correlation functions indicate that an exchange between fast and slow ions takes place on the timescale of the jumps themselves, i.e., the dynamical heterogeneities are short-lived. Hence, sites featuring fast and slow lithium dynamics, respectively, are intimately mixed. In addition, a backward correlation beyond the first neighbor shell for highly mobile ions and the presence of long-range dynamical heterogeneities suggest that fast ion migration occurs along preferential pathways in the glassy matrix. In the melt, we find no evidence for correlated back-and-forth motions and dynamical heterogeneities on the length scale of the next-neighbor distance.Comment: 12 pages, 13 figure

Mode-coupling theory for multiple-time correlation functions of tagged particle densities and dynamical filters designed for glassy systems

The theoretical framework for higher-order correlation functions involving multiple times and multiple points in a classical, many-body system developed by Van Zon and Schofield [Phys. Rev. E 65, 011106 (2002)] is extended here to include tagged particle densities. Such densities have found an intriguing application as proposed measures of dynamical heterogeneities in structural glasses. The theoretical formalism is based upon projection operator techniques which are used to isolate the slow time evolution of dynamical variables by expanding the slowly-evolving component of arbitrary variables in an infinite basis composed of the products of slow variables of the system. The resulting formally exact mode-coupling expressions for multiple-point and multiple-time correlation functions are made tractable by applying the so-called N-ordering method. This theory is used to derive for moderate densities the leading mode coupling expressions for indicators of relaxation type and domain relaxation, which use dynamical filters that lead to multiple-time correlations of a tagged particle density. The mode coupling expressions for higher order correlation functions are also succesfully tested against simulations of a hard sphere fluid at relatively low density.Comment: 15 pages, 2 figure