DYNASOR -- A tool for extracting dynamical structure factors and current correlation functions from molecular dynamics simulations

Perturbative treatments of the lattice dynamics are widely successful for many crystalline materials, their applicability is, however, limited for strongly anharmonic systems, metastable crystal structures and liquids. The full dynamics of these systems can, however, be accessed via molecular dynamics (MD) simulations using correlation functions, which includes dynamical structure factors providing a direct bridge to experiment. To simplify the analysis of correlation functions, here the dynasor package is presented as a flexible and efficient tool that enables the calculation of static and dynamical structure factors, current correlation functions as well as their partial counterparts from MD trajectories. The dynasor code can handle input from several major open source MD packages and thanks to its C/Python structure can be readily extended to support additional codes. The utility of dynasor is demonstrated via examples for both solid and liquid single and multi-component systems. In particular, the possibility to extract the full temperature dependence of phonon frequencies and lifetimes is emphasized

Dynasor - A Tool for Extracting Dynamical Structure Factors and Current Correlation Functions from Molecular Dynamics Simulations

Perturbative treatments of the lattice dynamics are widely successful for many crystalline materials; however, their applicability is limited for strongly anharmonic systems, metastable crystal structures and liquids. The full dynamics of these systems can, however, be accessed via molecular dynamics (MD) simulations using correlation functions, which includes dynamical structure factors providing a direct bridge to experiment. To simplify the analysis of correlation functions, here the dynasor package is presented as a flexible and efficient tool that enables the calculation of static and dynamical structure factors, current correlation functions as well as their partial counterparts from MD trajectories. The dynasor code can handle input from several major open source MD packages and thanks to its C/Python structure can be readily extended to support additional codes. The utility of dynasor is demonstrated via examples for both solid and liquid single and multi-component systems. In particular, the possibility to extract the full temperature dependence of phonon frequencies and lifetimes is\ua0emphasized

Ab initio simulations of clustering and precipitation in Al-Mg-Si alloys

A class of proposed coherent precipitate structures (Guinier-Preston zones) in the Al-Mg-Si alloy are investigated using first-prin- ciples density functional theory methods. The cluster expansion method is used to extract effective interaction parameters, providing the means for large scale energy calculations of alloy structures. The Mg1Si1 L1(0) structure and structures related to the Mg5Si6 beta \u27\u27 phase are studied in more detail, and e.g., precipitate/matrix interface energies are presented. Using direct first-principles calculations we show that the former phase is dynamically unstable and thus must be stabilized by the surrounding Al matrix. Monte Carlo simulations and free-energy techniques are used to study the Al rich side of the phase diagram with the current CE parameters, and kinetic Monte Carlo simulations are used to study clustering in the disordered phase. The implications of our findings are discussed in the framework of classical nucleation theories, and we outline possible nucleation mechanisms. (C) 2007 Elsevier B.V. All rights reserved

Early stages of phase separation using three-dimensional atom probe and atomistic modelling

The early stages of clustering in Al-14.8 Mg alloy have been studied using atomistic modelling and 3DAP analysis. Atomistic modelling was based on first-principles quantum mechanical calculations and Monte Carlo (MC) based techniques. A good agreement between the experimental results and simulated data was obtained showing the appropriateness of the theoretical approach used. It has been shown that after 720 h at 35 \ub0C (308 K) clustering of Mg was clearly detectable in the experimental alloy. The simulated data after 200 MC steps corresponded well to the experimental results obtained after ageing for 720 h