The Potential Energy Landscape and Mechanisms of Diffusion in Liquids

The mechanism of diffusion in supercooled liquids is investigated from the potential energy landscape point of view, with emphasis on the crossover from high- to low-T dynamics. Molecular dynamics simulations with a time dependent mapping to the associated local mininum or inherent structure (IS) are performed on unit-density Lennard-Jones (LJ). New dynamical quantities introduced include r2_{is}(t), the mean-square displacement (MSD) within a basin of attraction of an IS, R2(t), the MSD of the IS itself, and g_{loc}(t) the mean waiting time in a cooperative region. At intermediate T, r2_{is}(t) posesses an interval of linear t-dependence allowing calculation of an intrabasin diffusion constant D_{is}. Near T_{c} diffusion is intrabasin dominated with D = D_{is}. Below T_{c} the local waiting time tau_{loc} exceeds the time, tau_{pl}, needed for the system to explore the basin, indicating the action of barriers. The distinction between motion among the IS below T_{c} and saddle, or border dynamics above T_{c} is discussed.Comment: submitted to pr

The dynamical structure factor in topologically disordered systems

A computation of the dynamical structure factor of topologically disordered systems, where the disorder can be described in terms of euclidean random matrices, is presented. Among others, structural glasses and supercooled liquids belong to that class of systems. The computation describes their relevant spectral features in the region of the high frequency sound. The analytical results are tested with numerical simulations and are found to be in very good agreement with them. Our results may explain the findings of inelastic X-ray scattering experiments in various glassy systems.Comment: Version to be published in J. Chem. Phy

Potential energy landscape-based extended van der Waals equation

The inherent structures ({\it IS}) are the local minima of the potential energy surface or landscape, $U({\bf r})$, of an {\it N} atom system. Stillinger has given an exact {\it IS} formulation of thermodynamics. Here the implications for the equation of state are investigated. It is shown that the van der Waals ({\it vdW}) equation, with density-dependent $a$ and $b$ coefficients, holds on the high-temperature plateau of the averaged {\it IS} energy. However, an additional ``landscape'' contribution to the pressure is found at lower $T$. The resulting extended {\it vdW} equation, unlike the original, is capable of yielding a water-like density anomaly, flat isotherms in the coexistence region {\it vs} {\it vdW} loops, and several other desirable features. The plateau energy, the width of the distribution of {\it IS}, and the ``top of the landscape'' temperature are simulated over a broad reduced density range, $2.0 \ge \rho \ge 0.20$, in the Lennard-Jones fluid. Fits to the data yield an explicit equation of state, which is argued to be useful at high density; it nevertheless reproduces the known values of $a$ and $b$ at the critical point

Risk of Suicide Attempt in Adopted and Nonadopted Offspring

OBJECTIVE: We asked whether adoption status represented a risk of suicide attempt for adopted and nonadopted offspring living in the United States. We also examined whether factors known to be associated with suicidal behavior would mediate the relationship between adoption status and suicide attempt. METHODS: Participants were drawn from the Sibling Interaction and Behavior Study, which included 692 adopted and 540 nonadopted offspring and was conducted at the University of Minnesota from 1998 to 2008. Adoptees were systematically ascertained from records of 3 large Minnesota adoption agencies; nonadoptees were ascertained from Minnesota birth records. Outcome measures were attempted suicide, reported by parent or offspring, and factors known to be associated with suicidal behavior including psychiatric disorder symptoms, personality traits, family environment, and academic disengagement. RESULTS: The odds of a reported suicide attempt were ∼4 times greater in adoptees compared with nonadoptees (odds ratio: 4.23). After adjustment for factors associated with suicidal behavior, the odds of reporting a suicide attempt were reduced but remained significantly elevated (odds ratio: 3.70). CONCLUSIONS: The odds for reported suicide attempt are elevated in individuals who are adopted relative to those who are not adopted. The relationship between adoption status and suicide attempt is partially mediated by factors known to be associated with suicidal behavior. Continued study of the risk of suicide attempt in adopted offspring may inform the larger investigation of suicidality in all adolescents and young adults

Different Scenarios for Critical Glassy Dynamics

We study the role of different terms in the $N$-body potential of glass forming systems on the critical dynamics near the glass transition. Using a simplified spin model with quenched disorder, where the different terms of the real $N$-body potential are mapped into multi-spin interactions, we identified three possible scenarios. For each scenario we introduce a ``minimal'' model representative of the critical glassy dynamics near, both above and below, the critical transition lin e. For each ``minimal'' model we discuss the low temperature equilibrium dynamics.Comment: Completely revised version, 8 pages, 5 figures, typeset using EURO-LaTeX, Europhysics Letters (in press

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

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

Heterogeneous Dynamics, Marginal Stability and Soft Modes in Hard Sphere Glasses

In a recent publication we established an analogy between the free energy of a hard sphere system and the energy of an elastic network [1]. This result enables one to study the free energy landscape of hard spheres, in particular to define normal modes. In this Letter we use these tools to analyze the activated transitions between meta-bassins, both in the aging regime deep in the glass phase and near the glass transition. We observe numerically that structural relaxation occurs mostly along a very small number of nearly-unstable extended modes. This number decays for denser packing and is significantly lowered as the system undergoes the glass transition. This observation supports that structural relaxation and marginal modes share common properties. In particular theoretical results [2, 3] show that these modes extend at least on some length scale $l^*\sim (\phi_c-\phi)^{-1/2}$ where $\phi_c$ corresponds to the maximum packing fraction, i.e. the jamming transition. This prediction is consistent with very recent numerical observations of sheared systems near the jamming threshold [4], where a similar exponent is found, and with the commonly observed growth of the rearranging regions with compression near the glass transition.Comment: 6 pages, improved versio

NASA Langley Research Center Systems Analysis & Concepts Directorate Participation in the Exploration Systems Architecture Study

The NASA Langley Research Center (LaRC) Systems Analysis & Concepts Directorate (SACD) began studying human exploration missions beyond low Earth orbit (LEO) in the year 1999. This included participation in NASA s Decadal Planning Team (DPT), the NASA Exploration Team (NExT), Space Architect studies and Revolutionary Aerospace Systems Concepts (RASC) architecture studies that were used in formulating the new Vision for Space Exploration. In May of 2005, NASA initiated the Exploration Systems Architecture Study (ESAS). The primary outputs of the ESAS activity were concepts and functional requirements for the Crewed Exploration Vehicle (CEV), its supporting launch vehicle infrastructure and identification of supporting technology requirements and investments. An exploration systems analysis capability has evolved to support these functions in the past and continues to evolve to support anticipated future needs. SACD had significant roles in supporting the ESAS study team. SACD personnel performed the liaison function between the ESAS team and the Shuttle/Station Configuration Options Team (S/SCOT), an agency-wide team charged with using the Space Shuttle to complete the International Space Station (ISS) by the end of Fiscal Year (FY) 2010. The most significant of the identified issues involved the ability of the Space Shuttle system to achieve the desired number of flights in the proposed time frame. SACD with support from the Kennedy Space Center performed analysis showing that, without significant investments in improving the shuttle processing flow, that there was almost no possibility of completing the 28-flight sequence by the end of 2010. SACD performed numerous Lunar Surface Access Module (LSAM) trades to define top level element requirements and establish architecture propellant needs. Configuration trades were conducted to determine the impact of varying degrees of segmentation of the living capabilities of the combined descent stage, ascent stage, and other elements. The technology assessment process was developed and implemented by SACD as the ESAS architecture was refined. SACD implemented a rigorous and objective process which included (a) establishing architectural functional needs, (b) collection, synthesis and mapping of technology data, and (c) performing an objective decision analysis resulting in technology development investment recommendations. The investment recommendation provided budget, schedule, and center/program allocations to develop required technologies for the exploration architecture, as well as the identification of other investment opportunities to maximize performance and flexibility while minimizing cost and risk. A summary of the trades performed and methods utilized by SACD for the Exploration Systems Mission Directorate (ESAS) activity is presented along with how SACD is currently supporting the implementation of the Vision for Space Exploration