Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ab initio calculations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ab initio calculations for void swelling bias in α - and δ -plutonium

Void swelling can develop in materials under persistent irradiation when nonequilibrium vacancy and self-interstitial populations migrate under sufficiently asymmetric interaction biases. In conventional metals, the propensity is determined to the first approximation by comparing point-defect relaxation strains. We thus present DFT-based calculations of structures and formation energies and volumes of point defects in the α and the δ phases of plutonium. Here we discuss the pros and cons of various levels of electronic structure theory: spin polarization, spin-orbit coupling, and orbital polarization. Our results show that lattice defects in δ-Pu, in contrast to most fcc metals, have surprisingly small formation volumes. Equally unexpected are the large defect formation volumes found in the low-symmetry α-Pu phase. Both these unusual properties can be satisfactorily explained from defect-induced spin/orbital moment formation and destruction in the Pu phases. Surprisingly, the point defects in α-Pu are found to induce far larger transformation of the local electronic structure than in δ-Pu. When we use the calculated defect properties to estimate the classic void swelling bias in each of the phases, we find it to be unusually small in δ-Pu but likely much larger in α-Pu. Hence, swelling rates and mechanisms can diverge dramatically between the different phases of Pu. Especially in the transient regime before the formation of large defect clusters, the swelling rate of α-Pu can reliably be expected to be much larger than δ-Pu. However, accurate forecasts over longer times will require the conventional void-swelling theory to be modified to handle the complexities presented by the different Pu phases. As a case in point, we show the possible anomalous temperature dependence of vacancy properties in δ-Pu, caused by entropic contributions from defect-induced spin-lattice fluctuations. Such complications may affect defect-defect interactions and thus alter the void swelling bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A ray-tracing algorithm for ab-initio calculation of thermal load in undulator-based synchrotron beamlines

The OASYS suite and its powerful integration features are used to implement a ray-tracing algorithm to accurately calculate the thermal load in any component of an undulator-based synchrotron beamline. This is achieved by sampling and converting the SRW source of a given energy into a Shadow source and using the latter code to ray trace the full beamline. The accuracy of the algorithm is proved by reconstructing the full undulator radiation distribution through an aperture and comparing the result with direct calculaton of the total power using SRW. The algorithm is particularly suited to analyze cases with complex beamline layouts and optical elements, such as crystals, multilayers, and compound refractive lenses. Examples of its use to calculate the power load on elements of two of the feature beamlines at the Advanced Photon Source Upgrade Project and a comparison of the results with analytical calculations are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Harnessing the power of ab initio calculations, distributed computing and machine learning to efficiently locate extreme molecules for use in carbon-based solar cells (Final Technical Report)

The use of high-throughput virtual screening (HTVS) tools is a powerful tool to expedite the materials discovery of commercially relevant materials. In previous years, our group has developed a molecular discovery platform to generate libraries in order to obtain suitable candidates for different applications, starting from the Harvard Clean Energy Project. This platform is suitable to test in-silico on traditional supercomputing clusters and shared resources, for example, in the IBM World Community Grid. In this project, we used the molecular discovery platform to create and screen a library of candidates of organic photovoltaic (OPVs) molecules. Based on a set of candidates created with combinations of molecular moieties, we were able to filter, by conformation stability, the energy of electronic orbitals and approximated power conversion efficiencies (PCE). To improve the predictions of orbital energies calculated and the PCEs, we used Gaussian Process regression and two sets of molecules. These sets correspond to electronic structure calculations of a higher level of theory and experimental PCE values, respectively. Finally, we selected a subset of the best candidates (molecules with a PCE higher than 10%) to understand its absorbance properties with TD-DFT. This project has demonstrated the capabilities of our molecular discovery platform for HTVS. Finally, machine learning can help us to introduce more complex effects included in bulk conditions and computational intensive calculations on models.

14 SOLAR ENERGY↗
DOE OSTI2020

Comparing event generator predictions and ab-initio calculations of $\nu$-$^{12}$C neutral current quasi-elastic scattering at 1 GeV

The measurement of neutrino oscillations and exotic physics searches are important parts of the physics program in the near future, with new state-of-the-art experiments planned within the next decade. Future and modern experiments in these fields will make use of nuclear targets. Event Generators (EGs) are software used in the analysis of neutrino oscillation experiments. EGs use to predict kinematic observables for a range of neutrino energies. These simulations may lack physics captured by more rigorous theoretical calculations. This work compares EG performance to nuclear theory calculations by comparing observables generated in the two frameworks. We provide a common set of definitions between theory and experiment and assess the physics contained in EG simulations. Neutral current quasi-elastic (NCQE) scattering events for neutrinos and anti-neutrinos on a $^{12}$C target are simulated with a specific EG, NEUT, used by the T2K experiment for its analysis. The simulated cross sections are compared to analytic calculations from nuclear theory within the factorization scheme. We compare the NEUT implementation of two different models on nuclear spectral functions: the Relativistic Fermi Gas (RFG) and the correlated basis spectral function (CBF) to analytic calculations of the same models in the factorization scheme. For both nuclear physics models, we compare the appearance of features in the distributions relevant to experimental analyses. Qualitatively, the shape of the simulated distribution is similar to the one obtained through theory calculations; however, there are some discrepancies between the theory calculations and the NEUT simulation. While the EG simulations and analytic calculations with the same model of nuclear dynamics show similar overall features, there are still differences between the two.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A = 3 ( e , e ' ) x B ≥ 1 cross-section ratios and the isospin structure of short-range correlations

We study the relation between measured high-x B , high-Q 2 , helium-3 to tritium, (e,e') inclusive-scattering cross-section ratios and the relative abundance of high-momentum neutron-proton (np) and proton-proton (pp) short-range correlated (SRC) nucleon pairs in three-body (A=3) nuclei. In this study, analysis of this data using a simple pair-counting cross-section model suggested a much smaller np/pp ratio than previously measured in heavier nuclei, questioning our understanding of A=3 nuclei and, by extension, all other nuclei. Here we examine this finding using spectral-function-based cross-section calculations, with both an ab initio A=3 spectral function and effective Generalized Contact Formalism (GCF) spectral functions using different nucleon-nucleon interaction models. The ab initio calculation agrees with the data, showing good understanding of the structure of A=3 nuclei. An 8% uncertainty on the simple pair-counting model, as implied by the difference between it and the ab initio calculation, gives a factor of 5 uncertainty in the extracted np/pp ratio. Thus we see no evidence for the claimed "unexpected structure in the high-momentum wavefunction for hydrogen-3 and helium-3."

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The effect of rare-earth substitution on the Debye temperature of inorganic phosphors

In the quest to predict photoluminescent efficiency in rare-earth substituted inorganic phosphors, research has shown that materials consisting of ordered, rigid crystal structures tend to possess the highest photoluminescent quantum yields. A compound's Debye temperature (Θ D ), which can be calculated using ab initio calculations, is an ideal proxy for quantitatively comparing structural rigidity among different inorganic compounds, allowing potentially efficient phosphors to be selected from large crystal structure databases. However, the high computational cost of these calculations limits estimating Θ D for unsubstituted host crystal structures only. It is assumed that the low substitution concentration of the rare-earth luminescent center does not significantly influence a material's Debye temperature. This work evaluates the validity of this approximation by examining the effect of luminescent center substitution on a host structure's Θ D . Two well-known phosphors, (Y 1– x Ce x ) 3 Al 5 O 12 ( x = 0 – 0.05) and Ba 1-δ Eu δ MgAl 10 O 17 (δ = 0 – 0.15), were synthesized with varying rare-earth concentrations, while Θ D was computationally estimated and then determined by ultrasonic pulse-echo speed-of-sound and low-temperature heat capacity measurements. Overall, the ensuing results provide key implications for using Θ D as a proxy for structural rigidity in substituted inorganic compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Superconductivity in Y 4 RuGe 8 with a Vacancy-Ordered CeNiSi 2 -Type Superstructure

In this work, we report a new compound, Y 4 RuGe 8 , with a transition metal vacancy-ordered CeNiSi 2 -type superstructure, which has a superconducting transition at 1.3 K. Y 4 RuGe 8 crystals were grown by indium flux at relatively low temperatures (below 1273 K), which makes it possible to stabilize such a vacancy-ordered phase. The crystal structure of Y 4 RuGe 8 was solved by single-crystal X-ray diffraction and confirmed by transmission electron microscopy. The as-grown Y 4 RuGe 8 crystals are always twinned, crystallizing in the space group $P\bar{1}$(no. 2) with the lattice parameters a = 5.7680(1) Å, b = 8.2042(2) Å, c = 11.5093(3) Å, α = 79.696(1)degrees, β = 88.491(1)degrees, and γ = 79.637(2)degrees; this structure is a superstructure deriving from the higher symmetry CeNiSi 2 -type structure (Cmcm, no. 63) due to the ordering of Ru vacancies. The ordering of Ru sites breaks slightly distorted Ge planes in the CeNiSi 2 prototype into infinite cis-trans Ge chains in Y 4 RuGe 8 . The presence of bulk superconductivity in Y 4 RuGe 8 is well supported by zero resistance and a jump in specific heat at the critical transition temperature. The Sommerfeld coefficient (19 mJ K -2 mol -1 ) of the specific heat is greater than that (11 mJ K -2 mol -1 ) estimated using the bare density of states (4.7 states/eV/f.u.) from first-principles calculations. The ab initio calculations indicate that 4d electrons of both Y and Ru and 4p electrons of Ge are the main contributors to the total density of states at the Fermi level in Y 4 RuGe 8 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

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↗
DOE OSTI2021

Giant orbital polarization of Ni 2+ in a square planar environment

Finding large orbital polarization in Ni-based oxides has become a topic of paramount interest in recent years due to the prospect of finding superconductivity. In this study, we investigate the electronic structure of single-crystalline samples of Sr 2 CuO 3 and Ni-doped Sr 2 CuO 3 containing Cu/NiO 4 square planar units. Our polarization-dependent x-ray absorption spectroscopy experiments reveal extremely large orbital polarization (~63%) for Ni 2+ . Cluster calculations and ab initio calculations find that the giant orbital polarization arises due to the low-spin (S=0) configuration with two holes in the Ni 3d x 2 -y 2 orbital, contrary to the expected high-spin (S=1) state from Hund's first rule.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Directed gas phase preparation of ethynylallene (H 2 CCCHCCH; X 1 A') via the crossed molecular beam reaction of the methylidyne radical (CH; X 2 Π) with vinylacetylene (H 2 CCHCCH; X 1 A')

The gas-phase bimolecular reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A') was conducted at a collision energy of 20.3 kJ mol –1 under single collision conditions exploiting the crossed molecular beam experimental results merged with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations. The laboratory data reveal that the bimolecular reaction proceeds barrierlessly via indirect scattering dynamics through long-lived C 5 H 5 reaction intermediate(s) ultimately dissociating to C 5 H 4 isomers along with atomic hydrogen with the latter predominantly originating from the vinylacetylene reactant as confirmed by the isotopic substitution experiments in the D1-methylidyne–vinylacetylene reaction. Combined with ab initio calculations of the potential energy surface (PES) and statistical Rice–Ramsperger–Kassel–Marcus (RRKM) calculations, the experimental determined reaction energy of –146 ± 26 kJ mol –1 along with the distribution minimum of T(θ) at 90° and isotopic substitution experiments suggest ethynylallene (p1; ΔrG = –230 ± 4 kJ mol –1 ) as the dominant product. The ethynylallene (p1) may be formed with extensive rovibrational excitation, which would result in a lower maximum translational energy. Further, AIMD simulations reveal that the reaction dynamics leads to p1 (ethynylallene, 75%) plus atomic hydrogen with the dominant initial complex being i1 formed by methylidyne radical addition to the double C[double bond, length as m-dash]C bond in vinylacetylene. Altogether, combining the crossed molecular beam experimental results with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations, ethynylallene (p1) is expected to represent the dominant product in the reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A').

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Disentangling the confounding spectroscopy of C 1 molecules: Without symmetry as a guide, everything is allowed

The spectra of C 1 molecules are confounding in that each of the fundamental vibrational modes transform as the same irreducible representation (A) and hence each band consists of a seemingly random distribution of a-, b-, and c-type transitions. This is in contrast to higher symmetry molecules for which band types are readily deduced by simple symmetry rules. Herein, we present a method to simulate the convoluted rotational contours in the gas-phase spectra of C 1 molecules by combining existing ab initio calculations with Colin Western’s PGOPHER rotational contour program. Specifically, ab initio calculations in the NWC HEM suite of programs were employed to predict the components of the dipole moment derivatives along the principal axes of the moments of inertia. This information was then input into PGOPHER to model the fundamental band contours as a sum of a-, b-, and c-type transitions. This method was applied to simulate the rotational contour spectra of a series of representative C 1 molecules which were then compared against both ab initio stick spectra and experimentally measured broadband IR spectra from the Pacific Northwest National Laboratory infrared gas-phase database. In addition to providing further insight beyond what is revealed in a typical stick spectrum, the simulated contours showed good agreement with the measured spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Analysis of the rotationally-resolved infrared spectrum of CD 3 79 Br and CD 3 81 Br between 900 and 1400 cm -1

We report the infrared spectrum of CD 3 79 Br and CD 3 81 Br between 900 cm -1 and 1400 cm -1 has been analyzed at high resolution. In this region, two fundamental vibrational states, v 2 = 1 and v 5 = 1, an overtone state, v 3 = 2, and a combination state, v 3 = v 6 = 1, have been analyzed for both isotopologues. As found in other halomethane molecules, strong Coriolis resonances couple the v 2 = 1 state to the v 5 = 1 state, and similar resonances were also used to describe the coupling found between the v 5 = 1 state to the v 3 = 2 state. A new determination of the K-dependent constant A 0 was performed for each isotopologue through the use of perturbation-allowed transitions. The values for A 0 determined from the analysis are 2.6001898(26) and 2.6001905(27) for CD 3 79 Br and CD 3 81 Br, respectively. Ab initio calculations were performed to examine the efficacy of high-level theoretical calculations to accurately predict spectroscopic constants. The ab initio calculations had variable success with the largest errors associated with the prediction of the ΔB rotational constants of the interacting states

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Persistent exchange splitting in the chiral helimagnet $\mathrm{Cr}_{1/3}\mathrm{NbS}_2$

Here, using high-resolution angle-resolved photoemission spectroscopy and ab initio calculation, we systematically investigate the electronic structure of the chiral helimagnet Cr 1/3 NbS 2 and its temperature evolution. The comparison with NbS 2 suggests that the electronic structure of Cr 1/3 NbS 2 is strongly modified by the intercalation of Cr atoms. Our ab initio calculation, consistent with experimental results, suggests strong hybridization between Nb- and Cr-derived states near the Fermi level. In the chiral helimagnetic state (below the Curie temperature, T c ), we observe exchange splitting of the energy bands crossing the Fermi level, which follows the temperature evolution of the magnetic moment, suggesting a strong interaction between the conduction electrons and Cr spin moments. Interestingly, the exchange splitting persists far above T c with weak temperature dependence, in drastic contrast to the itinerant ferromagnetism described by the Stoner model, indicating the existence of short-range magnetic order. Our results provide important insights into the interplay between the electronic structure and magnetism in Cr 1/3 NbS 2 , which is helpful for understanding the microscopic mechanism of chiral helimagnetic ordering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Impact of Clustering on the Li 8 β Decay and Recoil Form Factors

We place unprecedented constraints on recoil corrections in the β decay of 8 Li, by identifying a strong correlation between them and the 8 Li ground state quadrupole moment in large-scale ab initio calculations. The results are essential for improving the sensitivity of high-precision experiments that probe the weak interaction theory and test physics beyond the standard model. In addition, our calculations predict a 2 + state of the α + α system that is energetically accessible to β decay but has not been observed in the experimental 8 Be energy spectrum, and has an important effect on the recoil corrections and β decay for the A = 8 systems. Furthermore, this state and an associated 0 + state are notoriously difficult to model due to their cluster structure and collective correlations, but become feasible for calculations in the ab initio symmetry-adapted no-core shell-model framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Local position-space two-nucleon potentials from leading to fourth order of chiral effective field theory

We present local, position-space chiral NN potentials through four orders of chiral effective field theory ranging from leading order (LO) to next-to-next-to-next-to-leading order (N 3 LO, fourth order) of the Δ-less version of the theory. The long-range parts of these potentials are fixed by the very accurate πN LECs as determined in the Roy-Steiner equations analysis. At the highest order (N 3 LO), the NN data below 190 MeV laboratory energy are reproduced with the respectable χ 2 /datum of 1.45. A comparison of the N 3 LO potential with the phenomenological Argonne v 18 (AV18) potential reveals substantial agreement between the two potentials in the intermediate range ruled by chiral symmetry, thus, providing a chiral underpinning for the phenomenological AV18 potential. Our chiral NN potentials may serve as a solid basis for systematic ab initio calculations of nuclear structure and reactions that allow for a comprehensive error analysis. In particular, the order by order development of the potentials will make possible a reliable determination of the truncation error at each order. Our new family of local position-space potentials differs from existing potentials of this kind by a weaker tensor force as reflected in relatively low D-state probabilities of the deuteron (P D ≲ 4.0 % for our N 3 LO potentials) and predictions for the triton binding energy above 8.00 MeV (from two-body forces alone). As a consequence, our potentials may lead to different predictions when applied to light and intermediate-mass nuclei in ab initio calculations and, potentially, help solve some of the outstanding problems in microscopic nuclear structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗