Iron under Earth's core conditions: Liquid-state thermodynamics and high-pressure melting curve

{\em Ab initio} techniques based on density functional theory in the projector-augmented-wave implementation are used to calculate the free energy and a range of other thermodynamic properties of liquid iron at high pressures and temperatures relevant to the Earth's core. The {\em ab initio} free energy is obtained by using thermodynamic integration to calculate the change of free energy on going from a simple reference system to the {\em ab initio} system, with thermal averages computed by {\em ab initio} molecular dynamics simulation. The reference system consists of the inverse-power pair-potential model used in previous work. The liquid-state free energy is combined with the free energy of hexagonal close packed Fe calculated earlier using identical {\em ab initio} techniques to obtain the melting curve and volume and entropy of melting. Comparisons of the calculated melting properties with experimental measurement and with other recent {\em ab initio} predictions are presented. Experiment-theory comparisons are also presented for the pressures at which the solid and liquid Hugoniot curves cross the melting line, and the sound speed and Gr\"{u}neisen parameter along the Hugoniot. Additional comparisons are made with a commonly used equation of state for high-pressure/high-temperature Fe based on experimental data.Comment: 16 pages including 6 figures and 5 table

New Algebraic Formulation of Density Functional Calculation

This article addresses a fundamental problem faced by the ab initio community: the lack of an effective formalism for the rapid exploration and exchange of new methods. To rectify this, we introduce a novel, basis-set independent, matrix-based formulation of generalized density functional theories which reduces the development, implementation, and dissemination of new ab initio techniques to the derivation and transcription of a few lines of algebra. This new framework enables us to concisely demystify the inner workings of fully functional, highly efficient modern ab initio codes and to give complete instructions for the construction of such for calculations employing arbitrary basis sets. Within this framework, we also discuss in full detail a variety of leading-edge ab initio techniques, minimization algorithms, and highly efficient computational kernels for use with scalar as well as shared and distributed-memory supercomputer architectures

First-principles results for electromagnetic properties of sdsd shell nuclei

In this work we present $ab~initio$ shell-model calculations for electric quadrupole moments and magnetic dipole moments of $sd$ shell nuclei using valence-space Hamiltonians derived with two $ab~initio$ approaches: the in-medium similarity renormalization group (IM-SRG) and the coupled-cluster effective interaction (CCEI). Results are in a reasonable agreement with the available experimental data as well as with the results from the phenomenological USDB effective interaction. This work will add more information to the available $ab~initio$ results for the spectroscopy of $sd$ shell nuclei.Comment: 13 pages, 3 figures, 3 table

Raman spectra of BN-nanotubes: Ab-initio and bond-polarizability model calculations

We present it ab-initio calculations of the non-resonant Raman spectra of zigzag and armchair BN nanotubes. In comparison, we implement a generalized bond-polarizability model where the parameters are extracted from first-principles calculations of the polarizability tensor of a BN sheet. For light-polarization along the tube-axis, the agreement between model and it ab-initio spectra is almost perfect. For perpendicular polarization, depolarization effects have to be included in the model in order to reproduce the it ab-initio Raman intensities.Comment: 4 pages, submitted to Phys. Rev. B rapid com

Long-range three-body atom-diatom potential for doublet Li3{}_3

An accurate long-range {\em ab initio} potential energy surface has been calculated for the ground state ${}^2A'$ lithium trimer in the frozen diatom approximation using all electron RCCSD(T). The {\em ab initio} energies are corrected for basis set superposition error and extrapolated to the complete basis limit. Molecular van der Waals dispersion coefficients and three-body dispersion damping terms for the atom-diatomic dissociation limit are presented from a linear least squares fit and shown to be an essentially exact representation of the {\em ab initio} surface at large range

Modeling the sorption dynamics of NaH using a reactive force field

We have parametrized a reactive force field for NaH, ReaxFFNaH, against a training set of ab initio derived data. To ascertain that ReaxFFNaH is properly parametrized, a comparison between ab initio heats of formation of small representative NaH clusters with ReaxFFNaH was done. The results and trend of ReaxFFNaH are found to be consistent with ab initio values. Further validation includes comparing the equations of state of condensed phases of Na and NaH as calculated from ab initio and ReaxFFNaH. There is a good match between the two results, showing that ReaxFFNaH is correctly parametrized by the ab initio training set. ReaxFFNaH has been used to study the dynamics of hydrogen desorption in NaH particles. We find that ReaxFFNaH properly describes the surface molecular hydrogen charge transfer during the abstraction process. Results on heat of desorption versus cluster size shows that there is a strong dependence on the heat of desorption on the particle size, which implies that nanostructuring enhances desorption process. To gain more insight into the structural transformations of NaH during thermal decomposition, we performed a heating run in a molecular dynamics simulation. These runs exhibit a series of drops in potential energy, associated with cluster fragmentation and desorption of molecular hydrogen. This is consistent with experimental evidence that NaH dissociates at its melting point into smaller fragments

Ab initio mechanical response: internal friction and structure of divacancies in silicon

This letter introduces ab initio study of the full activation-volume tensor of crystalline defects as a means to make contact with mechanical response experiments. We present a theoretical framework for prediction of the internal friction associated with divacancy defects and give the first ab initio value for this quantity in silicon. Finally, making connection with defect alignment studies, we give the first unambiguous resolution of the debate surrounding ab initio verification of the ground-state structure of the defect.Comment: 5 pages, 2 figures, submitted to PR

An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics

We present a new method which combines Car-Parrinello and Born-Oppenheimer molecular dynamics in order to accelerate density functional theory based ab-initio simulations. Depending on the system a gain in efficiency of one to two orders of magnitude has been observed, which allows ab-initio molecular dynamics of much larger time and length scales than previously thought feasible. It will be demonstrated that the dynamics is correctly reproduced and that high accuracy can be maintained throughout for systems ranging from insulators to semiconductors and even to metals in condensed phases. This development considerably extends the scope of ab-initio simulations.Comment: 4 pages, 3 figures; Accepted by Phys. Rev. Lett. for publicatio