Saddles in the energy landscape probed by supercooled liquids

We numerically investigate the supercooled dynamics of two simple model liquids exploiting the partition of the multi-dimension configuration space in basins of attraction of the stationary points (inherent saddles) of the potential energy surface. We find that the inherent saddles order and potential energy are well defined functions of the temperature T. Moreover, decreasing T, the saddle order vanishes at the same temperature (T_MCT) where the inverse diffusivity appears to diverge as a power law. This allows a topological interpretation of T_MCT: it marks the transition from a dynamics between basins of saddles (T>T_MCT) to a dynamics between basins of minima (T<T_MCT).Comment: 4 pages, 3 figures, to be published on PR

Instantaneous Normal Mode Analysis of Supercooled Water

We use the instantaneous normal mode approach to provide a description of the local curvature of the potential energy surface of a model for water. We focus on the region of the phase diagram in which the dynamics may be described by the mode-coupling theory. We find, surprisingly, that the diffusion constant depends mainly on the fraction of directions in configuration space connecting different local minima, supporting the conjecture that the dynamics are controlled by the geometric properties of configuration space. Furthermore, we find an unexpected relation between the number of basins accessed in equilibrium and the connectivity between them.Comment: 5 pages, 4 figure

Analytic computation of the Instantaneous Normal Modes spectrum in low density liquids

We analytically compute the spectrum of the Hessian of the Hamiltonian for a system of N particles interacting via a purely repulsive potential in one dimension. Our approach is valid in the low density regime, where we compute the exact spectrum also in the localized sector. We finally perform a numerical analysis of the localization properties of the eigenfunctions.Comment: 4 RevTeX pages, 4 EPS figures. Revised version to appear on Phys. Rev. Let

Energy landscape of a Lennard-Jones liquid: Statistics of stationary points

Molecular dynamics simulations are used to generate an ensemble of saddles of the potential energy of a Lennard-Jones liquid. Classifying all extrema by their potential energy u and number of unstable directions k, a well defined relation k(u) is revealed. The degree of instability of typical stationary points vanishes at a threshold potential energy, which lies above the energy of the lowest glassy minima of the system. The energies of the inherent states, as obtained by the Stillinger-Weber method, approach the threshold energy at a temperature close to the mode-coupling transition temperature Tc.Comment: 4 RevTeX pages, 6 eps figures. Revised versio

Molecular structural order and anomalies in liquid silica

The present investigation examines the relationship between structural order, diffusivity anomalies, and density anomalies in liquid silica by means of molecular dynamics simulations. We use previously defined orientational and translational order parameters to quantify local structural order in atomic configurations. Extensive simulations are performed at different state points to measure structural order, diffusivity, and thermodynamic properties. It is found that silica shares many trends recently reported for water [J. R. Errington and P. G. Debenedetti, Nature 409, 318 (2001)]. At intermediate densities, the distribution of local orientational order is bimodal. At fixed temperature, order parameter extrema occur upon compression: a maximum in orientational order followed by a minimum in translational order. Unlike water, however, silica's translational order parameter minimum is broad, and there is no range of thermodynamic conditions where both parameters are strictly coupled. Furthermore, the temperature-density regime where both structural order parameters decrease upon isothermal compression (the structurally anomalous regime) does not encompass the region of diffusivity anomalies, as was the case for water.Comment: 30 pages, 8 figure

Mean-atom-trajectory model for the velocity autocorrelation function of monatomic liquids

We present a model for the motion of an average atom in a liquid or supercooled liquid state and apply it to calculations of the velocity autocorrelation function $Z(t)$ and diffusion coefficient $D$. The model trajectory consists of oscillations at a distribution of frequencies characteristic of the normal modes of a single potential valley, interspersed with position- and velocity-conserving transits to similar adjacent valleys. The resulting predictions for $Z(t)$ and $D$ agree remarkably well with MD simulations of Na at up to almost three times its melting temperature. Two independent processes in the model relax velocity autocorrelations: (a) dephasing due to the presence of many frequency components, which operates at all temperatures but which produces no diffusion, and (b) the transit process, which increases with increasing temperature and which produces diffusion. Because the model provides a single-atom trajectory in real space and time, including transits, it may be used to calculate all single-atom correlation functions.Comment: LaTeX, 8 figs. This is an updated version of cond-mat/0002057 and cond-mat/0002058 combined Minor changes made to coincide with published versio

Harmonic Vibrational Excitations in Disordered Solids and the "Boson Peak"

We consider a system of coupled classical harmonic oscillators with spatially fluctuating nearest-neighbor force constants on a simple cubic lattice. The model is solved both by numerically diagonalizing the Hamiltonian and by applying the single-bond coherent potential approximation. The results for the density of states $g(\omega)$ are in excellent agreement with each other. As the degree of disorder is increased the system becomes unstable due to the presence of negative force constants. If the system is near the borderline of stability a low-frequency peak appears in the reduced density of states $g(\omega)/\omega^2$ as a precursor of the instability. We argue that this peak is the analogon of the "boson peak", observed in structural glasses. By means of the level distance statistics we show that the peak is not associated with localized states

String-like Clusters and Cooperative Motion in a Model Glass-Forming Liquid

A large-scale molecular dynamics simulation is performed on a glass-forming Lennard-Jones mixture to determine the nature of dynamical heterogeneities which arise in this model fragile liquid. We observe that the most mobile particles exhibit a cooperative motion in the form of string-like paths (``strings'') whose mean length and radius of gyration increase as the liquid is cooled. The length distribution of the strings is found to be similar to that expected for the equilibrium polymerization of linear polymer chains.Comment: 6 pages of RevTex, 6 postscript figures, uses epsf.st

Implementing biomarkers to predict motor recovery after stroke

BACKGROUND: There is growing interest in using biomarkers to predict motor recovery and outcomes after stroke. The PREP2 algorithm combines clinical assessment with biomarkers in an algorithm, to predict upper limb functional outcomes for individual patients. To date, PREP2 is the first algorithm to be tested in clinical practice, and other biomarker-based algorithms are likely to follow. PURPOSE: This review considers how algorithms to predict motor recovery and outcomes after stroke might be implemented in clinical practice. FINDINGS: There are two tasks: first the prediction information needs to be obtained, and then it needs to be used. The barriers and facilitators of implementation are likely to differ for these tasks. We identify specific elements of the Consolidated Framework for Implementation Research that are relevant to each of these two tasks, using the PREP2 algorithm as an example. These include the characteristics of the predictors and algorithm, the clinical setting and its staff, and the healthcare environment. CONCLUSIONS: Active, theoretically underpinned implementation strategies are needed to ensure that biomarkers are successfully used in clinical practice for predicting motor outcomes after stroke, and should be considered in parallel with biomarker developmen

LEM-3 is a midbody-tethered DNA nuclease that resolves chromatin bridges during late mitosis

Chromosome segregation and genome maintenance require the removal of DNA bridges that link chromosomes just before cells divide. Here the authors show that the LEM-3/Ankle1 nuclease processes DNA bridges before cells divide and define a previously undescribed genome integrity mechanism