Publications

Results 51–75 of 100
Skip to search filters

New density functional theory approaches for enabling prediction of chemical and physical properties of plutonium and other actinides

Wills, Ann E.

Density Functional Theory (DFT) based Equation of State (EOS) construction is a prominent part of Sandia's capabilities to support engineering sciences. This capability is based on amending experimental data with information gained from computational investigations, in parts of the phase space where experimental data is hard, dangerous, or expensive to obtain. A prominent materials area where such computational investigations are hard to perform today because of limited accuracy is actinide and lanthanide materials. The Science of Extreme Environment Lab Directed Research and Development project described in this Report has had the aim to cure this accuracy problem. We have focused on the two major factors which would allow for accurate computational investigations of actinide and lanthanide materials: (1) The fully relativistic treatment needed for materials containing heavy atoms, and (2) the needed improved performance of DFT exchange-correlation functionals. We have implemented a fully relativistic treatment based on the Dirac Equation into the LANL code RSPt and we have shown that such a treatment is imperative when calculating properties of materials containing actinides and/or lanthanides. The present standard treatment that only includes some of the relativistic terms is not accurate enough and can even give misleading results. Compared to calculations previously considered state of the art, the Dirac treatment gives a substantial change in equilibrium volume predictions for materials with large spin-orbit coupling. For actinide and lanthanide materials, a Dirac treatment is thus a fundamental requirement in any computational investigation, including those for DFT-based EOS construction. For a full capability, a DFT functional capable of describing strongly correlated systems such as actinide materials need to be developed. Using the previously successful subsystem functional scheme developed by Mattsson et.al., we have created such a functional. In this functional the Harmonic Oscillator Gas is providing the necessary reference system for the strong correlation and localization occurring in actinides. Preliminary testing shows that the new Hao-Armiento-Mattsson (HAM) functional gives a trend towards improved results for the crystalline copper oxide test system we have chosen. This test system exhibits the same exchange-correlation physics as the actinide systems do, but without the relativistic effects, giving access to a pure testing ground for functionals. During the work important insights have been gained. An example is that currently available functionals, contrary to common belief, make large errors in so called hybridization regions where electrons from different ions interact and form new states. Together with the new understanding of functional issues, the Dirac implementation into the RSPt code will permit us to gain more fundamental understanding, both quantitatively and qualitatively, of materials of importance for Sandia and the rest of the Nuclear Weapons complex.

More Details
TYPE SAND Report YEAR 2012
OSTIDOI

Calculating Hugoniots for molecular crystals from first principles

Wills, Ann E.; Wixom, Ryan R.; Mattsson, Thomas M.

Density Functional Theory (DFT) has over the last few years emerged as an indispensable tool for understanding the behavior of matter under extreme conditions. DFT based molecular dynamics simulations (MD) have for example confirmed experimental findings for shocked deuterium, enabled the first experimental evidence for a triple point in carbon above 850 GPa, and amended experimental data for constructing a global equation of state (EOS) for water, carrying implications for planetary physics. The ability to perform high-fidelity calculations is even more important for cases where experiments are impossible to perform, dangerous, and/or prohibitively expensive. For solid explosives, and other molecular crystals, similar success has been severely hampered by an inability of describing the materials at equilibrium. The binding mechanism of molecular crystals (van der Waals forces) is not well described within traditional DFT. Among widely used exchange-correlation functionals, neither LDA nor PBE balances the strong intra-molecular chemical bonding and the weak inter-molecular attraction, resulting in incorrect equilibrium density, negatively affecting the construction of EOS for undetonated high explosives. We are exploring a way of bypassing this problem by using the new Armiento-Mattsson 2005 (AM05) exchange-correlation functional. The AM05 functional is highly accurate for a wide range of solids, in particular in compression. In addition, AM05 does not include any van der Waals attraction, which can be advantageous compared to other functionals: Correcting for a fictitious van der Waals like attraction with unknown origin can be harder than correcting for a complete absence of all types of van der Waals attraction. We will show examples from other materials systems where van der Waals attraction plays a key role, where this scheme has worked well, and discuss preliminary results for molecular crystals and explosives.

More Details
TYPE Conference YEAR 2010
OSTI

On calculating the equilibrium structure of molecular crystals

Wills, Ann E.; Wixom, Ryan R.; Mattsson, Thomas M.

The difficulty of calculating the ambient properties of molecular crystals, such as the explosive PETN, has long hampered much needed computational investigations of these materials. One reason for the shortcomings is that the exchange-correlation functionals available for Density Functional Theory (DFT) based calculations do not correctly describe the weak intermolecular van der Waals' forces present in molecular crystals. However, this weak interaction also poses other challenges for the computational schemes used. We will discuss these issues in the context of calculations of lattice constants and structure of PETN with a number of different functionals, and also discuss if these limitations can be circumvented for studies at non-ambient conditions.

More Details
TYPE Conference YEAR 2010
OSTI
Results 51–75 of 100
Results 51–75 of 100