Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

In-Silico Analysis of High Refractive Index Materials Through Principles of Materials Design

The intent of the paper is to use specific principles of Materials Design that were developed and applied in the electronics industry for enabling understanding and design of improved high refractive index materials. Further, by combining first-principle based ab-initio, semiempirical interatomic potential methods, and machine learning approaches in conjunction with experimental data, we identified specific determinants of high refractive index materials, which can be critically applied for informing materials design and accelerating discovery. Specifically, it was demonstrated that chalcogenides and perovskites as bulk materials can exhibit higher refractive indices with appropriate engineering of specific aspects of the materials.

36 MATERIALS SCIENCE↗

Computational Studies of Solubilities of LiO 2 and Li 2 O 2 in Aprotic Solvents

Knowledge of the solubilities of Li 2 O 2 and LiO 2 in aprotic solvents is important for insight into the discharge and charge processes of Li-O 2 batteries, but these quantities are not well known. In this contribution, the solvation free energies of molecular LiO 2 and Li 2 O 2 in various organic solvents were calculated using various explicit and implicit solvent models, as well as ab initio molecular dynamics (AIMD) methods. Best estimates for the solvation energies from these calculations along with calculated lattice energies of Li 2 O 2 and LiO 2 were used to determine the solubility of bulk LiO 2 and Li 2 O 2 . The computed solubility of LiO 2 (1.8 × 10 -2 M) is about 15 orders higher than that of Li 2 O 2 (2.0 × 10 -17 M) due to a much less negative lattice energy of bulk LiO 2 compared to that of Li 2 O 2 . The difference in solubilities between LiO 2 and Li 2 O 2 likely will affect the nucleation and growth mechanisms and resulting morphologies of the products formed during battery discharge, influencing the performance of the battery cell. In conclusion, the calculated LiO 2 and Li 2 O 2 solubilities provide important information for fundamental studies of discharge and charge chemistries in Li-O 2 batteries.

25 ENERGY STORAGE↗

Thermodynamics of Plutonium Monocarbide from Anharmonic and Relativistic Theory

Thermodynamics of plutonium monocarbide is studied from first-principles theory that includes relativistic electronic structure and anharmonic lattice vibrations. Density-functional theory (DFT) is expanded to include orbital-orbital coupling in addition to the relativistic spin-orbit interaction for the electronic structure and it is combined with anharmonic, temperature dependent, lattice dynamics derived from the self-consistent ab initio lattice dynamics (SCAILD) method. The obtained thermodynamics are compared to results from simpler quasi-harmonic theory and experimental data. Formation enthalpy, specific heat, and Gibbs energy calculated from the anharmonic model are validated by direct comparison with a calculation of phase diagram (CALPHAD) assessment of PuC and sub-stochiometric PuC0.896. Overall, the theory reproduces CALPHAD results and measured data for PuC rather well, but the comparison is hampered by the sub-stoichiometric nature of plutonium monocarbide. It was found that a bare theoretical approach that ignores spin-orbit and orbital-orbital coupling (orbital polarization) of the plutonium 5f electrons promotes too soft phonons and Gibbs energies that are incompatible with that of the CALPHAD assessment of the experimental data. The investigation of PuC suggests that the electronic structure is accurately described by plutonium 5f electrons as “band like” and delocalized, but correlate through spin polarization, orbital polarization, and spin-orbit coupling, in analogy to previous findings for plutonium metal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermodynamics Modeling for Actinide Monocarbides and Mononitrides from First Principles

The high-temperature thermodynamical properties for the actinide monocarbides and mononitrides ThC, ThN, UC, UN, PuC, and PuN are calculated from first-principles electronic-structure theory. The electronic structure is modeled with density-functional theory (DFT) and is fully relativistic, including the spin-orbit interaction. Furthermore, the DFT is extended to account for orbital–orbital interactions, by means of a parameter-free orbital-polarization (OP) technique, that has proven to be essential for the 5f electrons in plutonium. Strong anharmonicity and the temperature dependence of the lattice vibrations are captured with the self-consistent ab initio lattice dynamics (SCAILD) method. The calculated free energies and heat capacities are compared to published results from quasi-harmonic (QH) theory, and experiments, where available. For the uranium and plutonium compounds, we make use of CALPHAD assessments to help evaluate the theory. Generally, our anharmonic relativistic approach compares well with both CALPHAD and experiments. For the thorium compounds, our theory is in good accord with QH modeling of the free energy at lower temperatures but for the heat capacity the comparison is less favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ab Initio Simulation of Structure and Properties in Ni-Based Superalloys: Haynes282 and Inconel740

The electronic structure, interatomic bonding, and mechanical properties of two supercell models of Ni-based superalloys are calculated using ab initio density functional theory methods. The alloys, Haynes282 and Inconel740, are face-centered cubic lattices with 864 atoms and eleven elements. These multi-component alloys have very complex electronic structure, bonding and partial-charge distributions depending on the composition and strength of the local bonding environment. We employ the novel concept of total bond order density (TBOD) and its partial components (PBOD) to ascertain the internal cohesion that controls the intricate balance between the propensity of metallic bonding between Ni, Cr and Co, and the strong bonds with C and Al. We find Inconel740 has slightly stronger mechanical properties than Haynes282. Both Inconel740 and Haynes282 show ductile natures based on Poisson’s ratio. Poisson’s ratio shows marginal correlation with the TBOD. Comparison with more conventional high entropy alloys with equal components are discussed.

36 MATERIALS SCIENCE↗

On the contribution of intramolecular zero point energy to the equation of state of solid H2

Experimental evidence shows that the internal zero-point energy of the H2 molecule exhibits a relatively strong pressure dependence in the solid as well as changing considerably upon condensation. It is shown that these effects contribute about 6% to the total sublimation energy and to the pressure in the solid state. Methods to modify the ab initio isolated pair potential to account for these environmental effects are discussed.

Chandrasekharan, V.↗

Theoretical Study of Indium Compounds of Interest for Organometallic Chemical Vapor Deposition

The structural. electronic and therinochemical properties of indium compounds which are of interest in halide transport and organometallic chemical vapor deposition processes have been studied by ab initio and statistical mechanics methods. The compounds reported include: indium halides and hydrides (InF, InCl, InCl3, InH, InH2, InH3); indium clusters (In2, In3); methylindium, dimethylindium, and their hydrogen derivatives [In(CH3), In(CH3)H, In(CH3)H2, In(CH3)2, In(CH3)2H]; dimethyl-indium dimer [In2(CH3)4], trimethyl-indium [In(CH3)3]; dehydrogenated methyl, dimethyl and trimethylindium [In(CH3)2CH2, In(CH3)CH2, In(CH2)], trimethylindium adducts with ammonia, trimethylamine and hydrazine [(CH3)3In:NH3, (CH3)3In:N(CH3)3, (CH3)3In:N(H2)N(H2)]; dimethylamino-indium and methylimino-indium [In(CH3)2(NH2), In(CH3)(NH)]; indium nitride and indium nitride dimer (InN, In2N2), indium phosphide, arsenide and antimonide ([InP, InAs, InSb). The predicted electronic properties are based on density functional theory calculations; the calculated thermodynamic properties are reported following the format of the JANAF (Joint Army, Navy, NASA, Air Force) Tables. Equilibrium compositions at two temperatures (298 and 1000 K) have been analyzed for groups of competing simultaneous reactions.

Cardelino, B. H.↗

Ab initio generation of thermal neutron scattering law for uranium dioxide

Neutronic analysis of moderated nuclear fission systems, e.g. nuclear reactors, relies upon an accurate prediction of the thermal neutron spectrum. Thermal neutrons have wavelengths on the order of the interatomic spacing between atoms in matter. Consequently, the scattering cross-sections of these neutrons are affected by the scattering law of the material (i.e. S(α,β), where α and β are dimensionless momentum and energy transfer variables), which describes the atomic structure and thermodynamic motion of the atomic nuclei. Uranium-dioxide (UO{sub 2}) is the primary fuel source for nuclear power reactors moderated by light water. In such systems (e.g. nuclear power reactors, spent fuel pools, etc.) the nuclear fuel comprises a large fraction of the system volume, approximately 1/3 in light-water reactors. Therefore, the scattering law of the fuel material may impact the thermal neutron spectrum and criticality of these systems. In this work previously established ab initio lattice dynamics (AILD) methods were utilized to calculate the phonon density of states (DOS), which acts as a fundamental input to the scattering law. Subsequently, the S(α,β) of UO{sub 2} was generated for publication as an ENDF thermal scattering library (TSL), which will be available for use in standard neutronics computer applications. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

Generation of Enrichment-Dependent Thermal Neutron Scattering Data

This work details the generation of enrichment-dependent thermal neutron scattering cross sections for several crucial uranium fuel compounds. The evaluations of the thermal scattering law (TSL) and associated cross sections for uranium dioxide (UO 2 ), uranium carbide (UC), and uranium nitride (UN) were performed using standard ab initio lattice dynamics (AILD) methods. The data for uranium metal was produced using a novel hybrid approach of molecular dynamics combined with lattice dynamics methods. 235 U enrichments of 5%, 10% (LEU+), 19.75% (HALEU), 93% (HEU), and 100% were considered, in addition to natural uranium. The enrichment-dependent masses and free atom cross sections were used in the generation of elastic and inelastic thermal neutron scattering cross sections, while the calculation of the phonon density of states (DOS) and resulting TSL considered only the natural isotopic composition of uranium. The use of an identical DOS for all enrichments is expected to have minimal impact on the final data, as the small change in uranium mass should not significantly affect lattice vibrations. The cross sections are shown to exhibit significant dependence on 235 U enrichment. The submission of this data to the National Nuclear Data Center (NNDC) for release in the ENDF/B-VIII.1 database should support the design of advanced reactor concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ab initio self-energy embedding for the photoemission spectra of NiO and MnO

The accurate ab initio simulation of periodic solids with strong correlations is one of the grand challenges of condensed matter. While mature methods exist for weakly correlated solids, the ab initio description of strongly correlated systems is an active field of research. In this work, we show results for the single-particle spectral function of the two correlated $\textit{d}$-electron solids NiO and MnO from self-energy embedding theory. Unlike earlier work, the theory does not use any adjustable parameters and is fully ab initio, while being able to treat both the strong correlation and the nonlocal screening physics of these materials. Here, we derive the method, discuss aspects of the embedding and choices of physically important orbitals, and compare our results to x-ray and angle-resolved photoemission spectroscopy as well as bremsstrahlung-isochromat spectroscopy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enabling Nuclear and Solar Thermal Propulsion with the Computational Laboratory: Faster and Cheaper Materials Characterization

Nuclear and solar thermal propulsion offer enhanced efficiencies for in-space travel which may open up the outer Solar System to both crewed and automated craft. However, these advanced propulsion systems require cutting-edge materials which are often not easily tested in the laboratory. Ab initio thermodynamics and other computational methods offer a way to test materials at high pressures and temperatures more quickly, and with reduced cost. Here, several case studies are presented which highlight the utility of modern simulation techniques. First, an extensive thermodynamic study was performed at the ab initio level to elucidate the role of hot hydrogen flow on the erosion of four refractory carbides to examine their appropriateness for use as channel coatings. A detailed analysis of reaction products informs estimates of erosion which are validated in comparison with heritage NERVA data. Specific material suggestions are made for both nuclear and solar designs. Secondly, we examined the mechanical behavior and chemical stability of various components of an idealized nuclear cermet core design using techniques spanning from the atomic to microscale. Next, the mechanical and chemical behavior of tungsten (a potential matrix material) and the stability of uranium mononitride (a potential ceramic fuel element), are characterized in the presence of hydrogen and at elevated temperatures with ab initio thermodynamics. Finally, the mechanical response and failure of tungsten under tensile strain is simulated at the micron scale with a combined finite element and dislocation dynamics approach. The resultant stress-strain data agrees well with experiments and leads to a workflow which can incorporate ab initio thermodynamic data into micron-scale simulations and thus provide real-world relevant estimates of erosion and stress-strain behavior. The power of this framework and our plans to use it in the future will be discussed.

William C Tucker↗

Ab-initio calculation of total x-ray scattering from molecules

Here, we present a method to calculate total x-ray scattering cross sections directly from ab-initio electronic wavefunctions in atoms and molecules. The approach can be used in conjunction with multiconfigurational wavefunctions and exploits analytical integrals of Gaussian-type functions over the scattering operator, which leads to accurate and efficient calculations. The results are validated by comparison to experimental results and previous theory for the molecules H 2 and CO 2 . Importantly, we find that the inelastic component of the total scattering varies strongly with molecular geometry. The approach is appropriate for use in conjunction with quantum molecular dynamics simulations for the analysis of new ultrafast x-ray scattering experiments, and to interpret accurate gas-phase scattering experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Use of an Adapting Database of Ab Initio Calculations To Generate Accurate Newtonian Dynamics

We develop and demonstrate a method to efficiently use density functional calculations to drive classical dynamics of complex atomic and molecular systems. The method has the potential to scale to systems and time scales unreachable with current ab initio molecular dynamics schemes. It relies on an adapting dataset of independently computed Hellmann–Feynman forces for atomic configurations endowed with a distance metric. The metric on configurations enables fast database lookup and robust interpolation of the stored forces. Here, we discuss mechanisms for the database to adapt to the needs of the evolving dynamics, while maintaining accuracy, and other extensions of the basic algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab Initio Potential Energy Surfaces and the Calculation of Accurate Vibrational Frequencies

Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within plus or minus 8 cm(exp -1) on average, and molecular bond distances are accurate to within plus or minus 0.001-0.003 Angstroms, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as vibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy will be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.

Lee, Timothy J.↗

Z2Pack: Numerical implementation of hybrid Wannier centers for identifying topological materials

The intense theoretical and experimental interest in topological insulators and semimetals has established band structure topology as a fundamental material property. Consequently, identifying band topologies has become an important, but often challenging, problem, with no exhaustive solution at the present time. Here, we compile a series of techniques, some previously known, that allow for a solution to this problem for a large set of the possible band topologies. The method is based on tracking hybrid Wannier charge centers computed for relevant Bloch states, and it works at all levels of materials modeling: continuous k ∙ p models, tight-binding models, and ab initio calculations. We apply the method to compute and identify Chern, $\mathbb{Z}$ 2 , and crystalline topological insulators, as well as topological semimetal phases, using real material examples. Moreover, we provide a numerical implementation of this technique (the Z2Pack software package) that is ideally suited for high-throughput screening of materials databases for compounds with nontrivial topologies. We expect that our work will allow researchers to (a) identify topological materials optimal for experimental probes, (b) classify existing compounds, and (c) reveal materials that host novel, not yet described, topological states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling surface motion effects in N{sub 2} dissociation on W(110): Ab initio molecular dynamics calculations and generalized Langevin oscillator model

Accurately modeling surface temperature and surface motion effects is necessary to study molecule-surface reactions in which the energy dissipation to surface phonons can largely affect the observables of interest. We present here a critical comparison of two methods that allow to model such effects, namely, the ab initio molecular dynamics (AIMD) method and the generalized Langevin oscillator (GLO) model, using the dissociation of N{sub 2} on W(110) as a benchmark. AIMD is highly accurate as the surface atoms are explicitly part of the dynamics, but this advantage comes with a large computational cost. The GLO model is much more computationally convenient, but accounts for lattice motion effects in a very approximate way. Results show that, despite its simplicity, the GLO model is able to capture the physics of the system to a large extent, returning dissociation probabilities which are in better agreement with AIMD than static-surface results. Furthermore, the GLO model and the AIMD method predict very similar energy transfer to the lattice degrees of freedom in the non-reactive events, and similar dissociation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimates of crystalline LiF thermal conductivity at high temperature and pressure by a Green-Kubo method

Here, given the unique optical properties of LiF, it is often used as an observation window in high-temperature and -pressure experiments; hence, estimates of its transmission properties are necessary to interpret observations. Since direct measurements of the thermal conductivity of LiF at the appropriate conditions are difficult, we resort to molecular simulation methods. Using an empirical potential validated against ab initio phonon density of states, we estimate the thermal conductivity of LiF at high temperatures (1000–4000 K) and pressures (100–400 GPa) with the Green-Kubo method. We also compare these estimates to those derived directly from ab initio data. To ascertain the correct phase of LiF at these extreme conditions, we calculate the (relative) phase stability of the B1 and B2 structures using a quasiharmonic ab initio model of the free energy. We also estimate the thermal conductivity of LiF in an uniaxial loading state that emulates initial stages of compression in high-stress ramp loading experiments and show the degree of anisotropy induced in the conductivity due to deformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗