Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 685 records · Page 38

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↗

Properties of infrared extrapolations in a harmonic oscillator basis

Here, the success and utility of effective field theory (EFT) in explaining the structure and reactions of few-nucleon systems has prompted the initiation of EFT-inspired extrapolations to larger model spaces in ab initio methods such as the no-core shell model (NCSM). In this contribution, we review and continue our studies of infrared (ir) and ultraviolet (uv) regulators of NCSM calculations in which the input is phenomenological NN and NNN interactions fitted to data. We extend our previous findings that an extrapolation in the ir cutoff with the uv cutoff above the intrinsic uv scale of the interaction is quite successful, not only for the eigenstates of the Hamiltonian but also for expectation values of operators, such as r 2 , considered long range. The latter results are obtained with Hamiltonians transformed by the similarity renormalization group (SRG) evolution. On the other hand, a possible extrapolation of ground state energies in the uv cutoff when the ir cutoff is below the intrinsic ir scale is not robust and does not agree with the ir extrapolation of the same data or with independent calculations using other methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Vibrational Frequencies and Spectroscopic Constants for 1(sup 3)A' HNC and 1(sup 3)A' HOC+ from High-Accuracy Quartic Force Fields

The spectroscopic constants and vibrational frequencies for the 1(sup 3)A' states of HNC, DNC, HOC+, and DOC+ are computed and discussed in this work. The reliable CcCR quartic force field based on high-level coupled cluster ab initio quantum chemical computations is exclusively utilized to provide the anharmonic potential. Then, second order vibrational perturbation theory and vibrational configuration interaction methods are employed to treat the nuclear Schroedinger equation. Second-order perturbation theory is also employed to provide spectroscopic data for all molecules examined. The relationship between these molecules and the corresponding 1(sup 3)A' HCN and HCO+ isomers is further developed here. These data are applicable to laboratory studies involving formation of HNC and HOC+ as well as astronomical observations of chemically active astrophysical environments.

Astrochemistry↗

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗

Quasielastic neutron scattering measurements and ab initio MD-simulations on single ion motions in molten NaF

The ionic stochastic motions in the molten alkali halide NaF are investigated by quasielastic neutron scattering and first principles molecular dynamics simulation. Quasielastic neutron scattering was employed to extract the diffusion behavior of the sodium ions in the melt. An extensive first principles based simulation on a box of up to 512 particles has been performed to complement the experimental data. From that large box, a smaller 64-particle box has then been simulated over a runtime of 60 ps. A good agreement between calculated and neutron data on the level of spectral shape has been obtained. The obtained sodium diffusion coefficients agree very well. The simulation predicts a fluorine diffusion coefficient similar to the sodium one. Applying the Nernst-Einstein equation, a remarkable large cross correlation between both ions can be deduced. The velocity cross correlations demonstrate a positive correlation between the ions over a period of 0.1 ps. That strong correlation is evidence that the unlike ions do not move completely statistically independent and have a strong association over a short period of time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗