Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio data”

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 703 records · Page 39

Morse-Smale Analysis of Ion Diffusion in Ab Initio Battery Materials Simulations

Ab initio molecular dynamics (AIMD) simulations are increasingly useful in modeling, optimizing and synthesizing materials in energy sciences. In solving Schrödinger’s equation, they generate the electronic structure of the simulated atoms as a scalar field. However, methods for analyzing these volume data are not yet common in molecular visualization. The Morse-Smale complex is a proven, versatile tool for topological analysis of scalar fields. In this paper, we apply the discrete Morse-Smale complex to analysis of first-principles battery materials simulations. We consider a carbon nanosphere structure used in battery materials research, and employ Morse-Smale decomposition to determine the possible lithium ion diffusion paths within that structure. Our approach is novel in that it uses the wavefunction itself as opposed distance fields, and that we analyze the 1-skeleton of the Morse-Smale complex to reconstruct our diffusion paths. Furthermore, it is the first application where specific motifs in the graph structure of the complete 1-skeleton define features, namely carbon rings with specific valence. We compare our analysis of DFT data with that of a distance field approximation, and discuss implications on larger classical molecular dynamics simulations.

25 ENERGY STORAGE↗

Electronic Properties and Magnetocaloric Effect of GdNi{sub 4}Si: Ab Initio Calculations, Mean Field Approximation, and Monte Carlo Simulation

We investigated the magnetic properties and magnetocaloric effect of GdNi{sub 4}Si by using a combination of first-principle density functional theory (DFT) calculations and Monte Carlo Simulation. Ab initio calculations show the stability of the ferromagnetic phase in this material. The Monte Carlo simulations combined with ab initio exchange coupling constants as input parameters lead to discuss the transition temperature, critical exponents, and isothermal entropy changes. The measured critical temperature and the isothermal magnetic entropy changes for a field change of 14T in good agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Energy and mechanical properties predictions in Fe-Ni binary system by ab initio calculations

Here, the fcc random solid solutions in Fe-Ni binary system were investigated using ab initio calculations. Density functional theory (DFT) was applied along with supercells generated by the Special Quasi-Random Structure (SQS) method to predict the enthalpy of mixing and magnetic and elastic properties. The results were compared with relevant experimental and computational data produced by other groups. The effect of SQS supercell size and relaxation procedure on the predicted properties of the alloys were also studied. Our calculations show that 16-atoms SQS model provides adequate approach to the mixing energy of the alloy. However, calculation of mechanical properties requires larger supercell. Excellent agreement between experimental and calculated elastic properties was found in fcc region of Fe-Ni system.

36 MATERIALS SCIENCE↗

X-ray Spectroscopy Fingerprints of Pristine and Functionalized Graphene

In this work, we demonstrate how to identify and characterize the atomic structure of pristine and functionalized graphene materials from a combination of computational simulation of X-ray spectra, on the one hand, and computer-aided interpretation of experimental spectra, on the other. Despite the enormous scientific and industrial interest, the precise structure of these 2D materials remains under debate. As we show in this study, a wide range of model structures from pristine to heavily oxidized graphene can be studied and understood with the same approach. We move systematically from pristine to highly oxidized and defective computational models, and we compare the simulation results with experimental data. Comparison with experiments is valuable also the other way around; this method allows us to verify that the simulated models are close to the real samples, which in turn makes simulated structures amenable to several computational experiments. Our results provide ab initio semiquantitative information and a new platform for extended insight into the structure and chemical composition of graphene-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Text-mined dataset of inorganic materials synthesis recipes

Abstract Materials discovery has become significantly facilitated and accelerated by high-throughput ab-initio computations. This ability to rapidly design interesting novel compounds has displaced the materials innovation bottleneck to the development of synthesis routes for the desired material. As there is no a fundamental theory for materials synthesis, one might attempt a data-driven approach for predicting inorganic materials synthesis, but this is impeded by the lack of a comprehensive database containing synthesis processes. To overcome this limitation, we have generated a dataset of “codified recipes” for solid-state synthesis automatically extracted from scientific publications. The dataset consists of 19,488 synthesis entries retrieved from 53,538 solid-state synthesis paragraphs by using text mining and natural language processing approaches. Every entry contains information about target material, starting compounds, operations used and their conditions, as well as the balanced chemical equation of the synthesis reaction. The dataset is publicly available and can be used for data mining of various aspects of inorganic materials synthesis.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

First-principles predicted low-energy structures of NaSc(BH{sub 4}){sub 4}

According to previous interpretations of experimental data, sodium-scandium double-cation borohydride NaSc(BH{sub 4}){sub 4} crystallizes in the crystallographic space group Cmcm where each sodium (scandium) atom is surrounded by six scandium (sodium) atoms. A careful investigation of this phase based on ab initio calculations indicates that the structure is dynamically unstable and gives rise to an energetically and dynamically more favorable phase with C222{sub 1} symmetry and nearly identical x-ray diffraction pattern. By additionally performing extensive structural searches with the minima-hopping method we discover a class of new low-energy structures exhibiting a novel structural motif in which each sodium (scandium) atom is surrounded by four scandium (sodium) atoms arranged at the corners of either a rectangle with nearly equal sides or a tetrahedron. These new phases are all predicted to be insulators with band gaps of 7.9–8.2 eV. Finally, we estimate the influence of these structures on the hydrogen-storage performance of NaSc(BH{sub 4}){sub 4}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The sea of quarks and antiquarks in the nucleon

The quark and gluon structure of the proton has been under intense experimental and theoretical investigation for five decades. Even for the distributions of the well-studied valence quarks, challenges such as the value of the down quark to up quark ratio at high fractional momenta remain. Much of the sea of quark-antiquark pairs emerges from the splitting of gluons and is well described by perturbative evolution in quantum chromodynamics. However, experiments confirm that there is a non-perturbative component to the sea that is not well understood and hitherto has been difficult to calculate with ab initio non-perturbative methods. This non-perturbative structure shows up, perhaps most directly, in the flavor dependence of the sea antiquark distributions. Even though some of the general trends can be reproduced by models, there are features of the data that do not seem to be well described. Conclusively, this article discusses the experimental situation, the status of calculations and models, and the directions where these studies will progress in the near future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Renormalized approach to neutrinoless double-$β$ decay

The process at the heart of neutrinoless double-$β$ decay, $nn →\ pp\ e⁻ e⁻$ induced by a light Majorana neutrino, is investigated in pionless and chiral effective field theory. We show in various regularization schemes the need to introduce a short-range lepton-number-violating operator at leading order, confirming earlier findings. We demonstrate that such a short-range operator is only needed in spin-singlet $S$-wave transitions, while leading-order transitions involving higher partial waves depend solely on long-range currents. Calculations are extended to include next-to-leading-order corrections in perturbation theory, where to this order no additional undetermined parameters appear. We establish a connection based on chiral symmetry between neutrinoless double-$β$ decay and nuclear charge-independence breaking induced by electromagnetism. Data on the latter confirm the need for a leading-order short-range operator but do not allow for a full determination of the corresponding lepton-number-violating coupling. Using a crude estimate of this coupling, we perform ab initio calculations of the matrix elements for neutrinoless double-$β$ decay for 6 He and 12 Be. We speculate on the phenomenological impact of the leading short-range operator on the basis of these results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗

Dicarbon Formation in Collisions of Two Carbon Atoms

Abstract Radiative association cross sections and rates are computed, using a quantum approach, for the formation of C 2 molecules (dicarbon) during the collision of two ground-state C( 3 P ) atoms. We find that transitions originating in the C 1 Π g , d 3 Π g , and 1 5 Π u states are the main contributors to the process. The results are compared and contrasted with previous results obtained from a semiclassical approximation. New ab initio potential curves and transition dipole moment functions have been obtained for the present work using the multi-reference configuration interaction approach with the Davidson correction (MRCI+Q) and aug-cc-pCV5Z basis sets, substantially increasing the available molecular data on dicarbon. Applications of the current computations to various astrophysical environments and laboratory studies are briefly discussed, focusing on these rates.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization of silicon-carbon clusters by infrared laser spectroscopy. The nu 1 band of SiC4

The nu 1 fundamental vibration of linear SiC4 has been observed by infrared diode laser spectroscopy of a supersonic cluster beam. Twenty-four rovibrational transitions were measured in the spectral region of 2094.6 to 2097.1 cm-1, the rotational temperature was 10 K. A combined least-squares fit of these transitions with previously reported microwave data yielded the following molecular constants: nu 1 = 2095.45806(37) cm-1, B" = 0.051161131(52) cm-1, and B' = 0.0509157(96) cm-1. These results are compared to vibrational spectroscopy measurements of SiC4 trapped in a solid Ar matrix and to ab initio calculations.

NASA Discipline Exobiology↗

Ab initio calculation of thermodynamic potentials and entropies for superionic water

We construct thermodynamic potentials for two superionic phases of water [with body-centered cubic (bcc) and face-centered cubic (fcc) oxygen lattice] using a combination of density functional theory (DFT) and molecular dynamics simulations (MD). For this purpose, a generic expression for the free energy of warm dense matter is developed and parametrized with equation of state data from the DFT-MD simulations. A second central aspect is the accurate determination of the entropy, which is done using an approximate two-phase method based on the frequency spectra of the nuclear motion. The boundary between the bcc superionic phase and the ices VII and X calculated with thermodynamic potentials from DFT-MD is consistent with that directly derived from the simulations. As a result, differences in the physical properties of the bcc and fcc superionic phases and their impact on interior modeling of water-rich giant planets are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Understanding the presence of vacancy clusters in ZnO from a kinetic perspective

Vacancy clusters have been observed in ZnO, but their physical origin remains elusive. This is because under thermal equilibrium clusters typically have higher formation energies. Combining first-principles calculation and kinetic Monte Carlo (kMC) simulation, we determine the roles of non-equilibrium kinetics on the vacancy clustering. We find that clustering starts with the formation of Zn and O vacancy pairs (VZn-Vo), which further grow by attracting additional mono-vacancies. At this stage, vacancy diffusivity becomes crucial: due to the larger diffusivity of VZn compared to VO, more VZn-abundant clusters are formed than VO-abundant clusters. The large dissociation energies, e.g., 3 eV for (VZn-Vo), suggest that, once formed, it is difficult for the clusters to dissociate. As such, thermal annealing will increase the size of the clusters by promoting mono-vacancy diffusion, in agreement with experiment.

Ab initio calculations↗

OpenMP in VASP: Threading and SIMD

Abstract The Vienna Ab initio Simulation Package (VASP) is a widely used electronic structure code that originally exploits process‐level parallelism through the Message Passing Interface (MPI) for work distribution within and across nodes. Architectural changes of modern parallel processors urge programmers to address thread‐ and data‐level parallelism as well to benefit most from the available compute resources within a node. We describe for VASP how to approach for an MPI + OpenMP parallelization including data‐level parallelism through OpenMP SIMD constructs together with a generic high‐level vector coding scheme. We can demonstrate an improved scalability of VASP and more than 20% gain over the MPI‐only version as well as a 2× increased performance of collective operations using the multiple‐endpoint MPI feature. The high‐level vector coding scheme applied to VASP's general gradient approximation routine gives up 9× performance gain on AVX512 platforms with the Intel compiler.

Wende, Florian↗

Identification of high-dielectric constant compounds from statistical design

Abstract The discovery of high-dielectric materials is crucial to increasing the efficiency of electronic devices and batteries. Here, we report three previously unexplored materials with very high dielectric constants (69 < ϵ < 101) and large band gaps (2.9 < E g (eV) < 5.5) obtained by screening materials databases using statistical optimization algorithms aided by artificial neural networks (ANN). Two of these new dielectrics are mixed-anion compounds (Eu 5 SiCl 6 O 4 and HoClO) and are shown to be thermodynamically stable against common semiconductors via phase diagram analysis. We also uncovered four other materials with relatively large dielectric constants (20 < ϵ < 40) and band gaps (2.3 < E g (eV) < 2.7). While the ANN training-data are obtained from the Materials Project, the search-space consists of materials from the Open Quantum Materials Database (OQMD)—demonstrating a successful implementation of cross-database materials design. Overall, we report the dielectric properties of 17 materials calculated using ab initio calculations, that were selected in our design workflow. The dielectric materials with high-dielectric properties predicted in this work open up further experimental research opportunities.

36 MATERIALS SCIENCE↗

New insights into the electronic structure of α-U and δ-Pu

Here, this work presents the results of a theoretical study of the electronic structure of two actinide metals, α -U and δ -Pu. We compare our ab-initio results obtained with the recently developed self-consistent Vertex corrected GW approach with previously published experimental measurements such as photo-electron spectroscopy, for the occupied density of states, and bremsstralung isochromat spectroscopy (BIS) and inverse photo-electron spectroscopy (IPES), for the unoccupied density of states. Our ab-initio approach includes all important relativistic effects (it is based on Dirac’s equation) and it represents the first application of the Vertex corrected GW approach in the physics of actinides. Overall, our theoretical results are in good agreement with the experimental data, which supports the level of approximations which our theoretical method is based upon. By comparing our vertex corrected GW results with our results obtained with less sophisticated approaches (local density approximation and self-consistent GW) we differentiate the strength of correlation effects in Uranium and Plutonium. Also, our theoretical results allow us to elucidate the subtle differences between the previously published experimental BIS and IPES data on the unoccupied density of states in α -U.

36 MATERIALS SCIENCE↗

Precision mass measurements of magnesium isotopes and implications for the validity of the isobaric mass multiplet equation

If the mass excess of neutron-deficient nuclei and their neutron-rich mirror partners are both known, it can be shown that deviations of the isobaric mass multiplet equation (IMME) in the form of a cubic term can be probed. Such a cubic term was probed by using the atomic mass of neutron-rich magnesium isotopes measured using the TITAN Penning trap and the recently measured proton-separation energies of 29 Cl and 30 Ar. The atomic mass of 27 Mg was found to be within 1.6σ of the value stated in the Atomic Mass Evaluation. The atomic masses of 28,29 Mg were measured to be both within 1σ, while being 7 and 33 times more precise, respectively. Using the 29 Mg mass excess and previous measurements of 29 Cl, we uncovered a cubic coefficient of d = 28(7)keV, which is the largest known cubic coefficient of the IMME. This departure, however, could also be caused by experimental data with unknown systematic errors. Hence there is a need to confirm the mass excess of 28 S and the one-neutron separation energy of 29 Cl, which have both come from a single measurement. Lastly, our results were compared with ab initio calculations from the valence-space in-medium similarity renormalization group, resulting in a good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗