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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

HTESP (High-throughput electronic structure package): A package for high-throughput ab initio calculations

High-throughput ab initio calculations are the indispensable parts of data-driven discovery of new materials with desirable properties, as reflected in the establishment of several online material databases. The accumulation of extensive theoretical data through computations enables data-driven discovery by constructing machine learning and artificial intelligence models to predict novel compounds and forecast their properties. Efficient usage and extraction of data from these existing online material databases can accelerate the next stage materials discovery that targets different and more advanced properties, such as electron–phonon coupling for phonon-mediated superconductivity. However, extracting data from these databases, generating tailored input files for different ab initio calculations, performing such calculations, and analyzing new results can be demanding tasks. Here, in this work, we introduce a software package named “HTESP” (High-Throughput Electronic Structure Package) written in Python and Bash languages, which automates the entire workflow including data extraction, input file generation, calculation submission, result collection and plotting. Our HTESP will help speed up future computational materials discovery processes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ab initio calculations of scattering and fusion reactions

We review recent advancements achieved within the ab initio no-core shell model with continuum (NCSMC), a unified first-principle (or, ab initio) approach to nuclear bound and continuum states. We highlight its application to reactions of astrophysical interests, including the 7 Be(p, γ) 8 B solar fusion and α-α scattering, and to the structure of the exotic 12 Be nucleus. Furthermore, our findings provide a solid foundation for integrating ab initio calculations with experimental measurements, allowing for more precise evaluations of key thermonuclear reaction rates at astrophysical energies, and offer new insights into the factors that contribute to the emergence of clustering.

Quaglioni, Sofia [Lawrence Livermore National Labo↗
DOE OSTI2025

Ab initio calculations of the carbon and oxygen isotopes: Energies, correlations, and superfluid pairing

We perform ab initio nuclear lattice calculations of the neutron-rich carbon and oxygen isotopes using high-fidelity chiral interactions. We find good agreement with the observed binding energies and compute correlations associated with each two-nucleon interaction channel. For the isospin T = 1 channels, we show that the dependence on T z provides a measure of the correlations among the extra neutrons in the neutron-rich nuclei. For the spin-singlet S-wave channel, we observe that any paired neutron interacts with the nuclear core as well as its neutron pair partner, while any unpaired neutron interacts primarily with only the nuclear core. For the other partial waves, the correlations among the extra neutrons grow more slowly and smoothly with the number of neutrons. These general patterns are observed in both the carbon and oxygen isotopes and may be universal features that appear in many neutron-rich nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Ab initio calculations of overlap integrals for μ → e conversion in nuclei

The rate for μ → e conversion in nuclei is set to provide the most stringent test of lepton-flavor symmetry and a window into physics beyond the Standard Model. However, to disentangle new lepton-flavor-violating interactions, in combination with information from μ → ey and μ → 3e, it is critical that uncertainties at each step of the analysis be controlled and fully quantified. In this regard, nuclear response functions related to the coupling to neutrons are notoriously problematic, since they are not directly constrained by experiment. We address these shortcomings by combining ab initio calculations with a recently improved determination of charge distributions from electron scattering by exploiting strong correlations among charge, point-proton, and point-neutron radii and densities. We present overlap integrals for 27 Al, 48 Ca, and 48 Ti including full covariance matrices, allowing, for the first time, for a comprehensive consideration of nuclear structure uncertainties in the interpretation of μ → e experiments.

ab initio calculations↗
DOE OSTI2026

Ab initio calculation of atomic solid hydrogen phases based on Gutzwiller many-body wave functions

We apply two ab initio many-body methods based on Gutzwiller wave functions, i.e., correlation matrix renormalization theory (CMRT) and Gutzwiller conjugate gradient minimization (GCGM), to the study of crystalline phases of atomic hydrogen. Both methods avoid empirical Hubbard U parameters and are free from double-counting issues. CMRT employs a Gutzwiller-type approximation that enables efficient calculations, while GCGM goes beyond this approximation to achieve higher accuracy at higher computational cost. By benchmarking against available quantum Monte Carlo (QMC) results, we demonstrate that while both methods are more accurate than the widely used density-functional theory, GCGM systematically captures additional correlation energy missing in CMRT, leading to significantly improved total energy predictions. We also show that by including the correlation energy Ec from local density approximation in the CMRT calculation, CMRT + E c produces energy in better agreement with the QMC results in these hydrogen lattice systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Improved structure of calcium isotopes from ab initio calculations

The in-medium similarity renormalization group (IMSRG) is a powerful and flexible many-body method to compute the structure of nuclei starting from nuclear forces. Recent developments have extended the IMSRG from its standard truncation at the normal-ordered two-body level, the IMSRG(2), to a precision approximation including normal-ordered three-body operators, the IMSRG(3)-N 7 . This improvement provides a more precise solution to the many-body problem and makes it possible to quantify many-body uncertainties in IMSRG calculations. We explore the structure of 44,48,52 Ca using the IMSRG(3)-N 7 , focusing on understanding existing discrepancies of the IMSRG(2) to experimental results. We find a significantly better description of the first 2 + excitation energy of 48 Ca, improving the description of the shell closure at N=28. At the same time, we find that the IMSRG(3)-N 7 corrections to charge radii do not resolve the systematic underprediction of the puzzling large charge radius difference between 52 Ca and 48 Ca. We present estimates of many-body uncertainties of IMSRG(2) calculations applicable also to other systems based on the size extensivity of the method.

39 ≤ A ≤ 58↗
DOE OSTI2024

Direct ab initio calculation of the 4 He nuclear electric dipole polarizability

The calculation of nuclear electromagnetic sum rules by directly diagonalizing the nuclear Hamiltonian in a large basis is numerically challenging and has not been performed for A>2 nuclei. With the significant progress of high performance computing, we show that calculating sum rules using numerous discretized continuum states obtained by directly diagonalizing the ab initio no-core shell model Hamiltonian is achievable numerically. Specifically, we calculate the 4 He electric dipole (E1) polarizability, that is an inverse energy weighted sum rule, employing the Daejeon16 NN interaction. We demonstrate that the calculations are numerically tractable as the dimension of the basis increases and are convergent. Our results for the 4 He electric dipole polarizability are consistent with the most recent experimental data and are compared with those of other theoretical studies employing different techniques and various interactions.

Astronomy & Astrophysics↗
DOE OSTI2026

Ab initio calculations of monopole sum rules: From finite nuclei to infinite nuclear matter

We compute moments of the isoscalar monopole response of 𝑁 = 𝑍 closed-shell nuclei based on chiral nucleon-nucleon plus three-nucleon interactions. We employ the random-phase approximation (RPA) and two ab initio many-body approaches, the in-medium similarity renormalization group (IMSRG) and coupled-cluster theory (CC). In the IMSRG framework, the moments are obtained as ground-state expectation values, whereas in the CC approach, they are evaluated through excited-state calculations. We find good agreement between the IMSRG and CC results across all nuclei studied. RPA provides a reasonable approximation to the correlated methods if the interaction is soft. From the calculated moments, we extract average energies of the monopole response, compute finite-nucleus incompressibilities, and estimate the incompressibility of symmetric nuclear matter by a fit to a leptodermous expansion. Our extrapolated values are lower than those obtained in nuclear-matter calculations with the same interactions, but the values are consistent with phenomenological ranges.

Bonaiti, Francesca [Michigan State Univ., East Lan↗
DOE OSTI2024

A comprehensive study on three typical photoacid generators using photoelectron spectroscopy and ab initio calculations

Conducting a comprehensive molecular-level evaluation of a photoacid generator (PAG) and its subsequent impact on lithography performance can facilitate the rational design of a promising 193 nm photoresist tailored to specific requirements. In this study, we integrated spectroscopy and computational techniques to meticulously investigate the pivotal factors of three prototypical PAG anions, p-toluenesulfonate (pTS - ), 2-(trifluoromethyl)benzene-1-sulfonate (TFMBS - ), and triflate (TF - ), in the lithography process. Our findings reveal a significant redshift in the absorption spectra caused by specific PAG anions, attributed to their involvement in electronic transition processes, thereby enhancing the transparency of the standard PAG cation, triphenylsulfonium (TPS + ), particularly at ~193 nm. Furthermore, the electronic stability of PAG anions can be enhanced by solvent effects with varying degrees of strength. Here we observed the lowest vertical detachment energy of 6.6 eV of pTS - in PGMEA solution based on the polarizable continuum model, which prevents anion loss at 193 nm lithography. In addition, our findings indicate gas-phase proton affinity values of 316.4 kcal/mol for pTS - , 308.1 kcal/mol for TFMBS - , and 303.2 kcal/mol for TF - , which suggest the increasing acidity strength, yet even the weakest acid pTS - is still stronger than strong acid HBr. The photolysis of TPS + -based PAG, TPS + ·pTS - , generated an excited state leading to homolysis bond cleavage with the lowest reaction energy of 83 kcal/mol. Overall, the PAG anion pTS - displayed moderate acidity, possessed the lowest photolysis reaction energy, and demonstrated an appropriate redshift. These properties collectively render it a promising candidate for an effective acid producer.

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

Edge dislocation mediated anomalous charge transfer in face centered cubic high entropy alloys

The charge transfer in alloys dictates their structural and functional properties. The presence of the line defect can influence the charge transfer characteristics in alloys. The edge dislocation-volume misfit interaction in high entropy alloys is a critical area of interest concerning the development of accurate solid solution strengthening models for these materials. The evolution of the volume misfit of atoms is either studied empirically using Vegard's law or mechanistically using \emph{ab initio} calculations, without considering the effect of edge dislocation on charge transfer explicitly. In this work, we show the influence of edge dislocation on the charge transfer in CoCrFeMnNi, CoCrNi and CoNi using large-scale \emph{ab initio} calculations. The anomalous charge transfer characteristics, in terms of deviations from electronegativity trends, are demonstrated, which are dictated by electronegativity equalisation rather than pair-atom interactions. The deviation of atomic volume in the compressive and tensile stress regions of the edge dislocation is rationalised in terms of anomalous magneto-volume fluctuations in such alloys.

Ghosh, Swarnava [ORNL] (ORCID:0000000338005264)↗
DOE OSTI2024

Reversible hydrogen storage in multilayer graphane: Lattice dynamics, compressibility, and heat capacity studies

Multilayer graphane (hydride of graphite) is a crystalline hydrocarbon of composition CH, which can be synthesized from graphite and molecular hydrogen at pressures above 2GPa [V.E. Antonov et al. Carbon 100 (2016) 465]. Using X-ray diffraction, this compound was tentatively identified as the “graphane II” phase of 3D-graphane predicted by ab initio calculations [X.-D. Wen et al. PNAS 108 (2011) 6833] and consisting of layers of 2D-graphane in the “chair” conformation. When heated in a vacuum, the compound does not form any intermediate hydrocarbons and reversibly decomposes back into graphite and hydrogen at 770–920 K. In the present work, almost single-phase samples of graphite hydride and deuteride were synthesized at 7.4 GPa and 870 K. Their investigation by inelastic neutron scattering supplemented by ab initio calculations gave spectra g(E) of the phonon density of states with a gap of about 15 meV at approx. 100 meV, which is a unique identifier for the chair form of graphane. The equation of state V(P) of the hydride was studied at room temperature and hydrogen pressures up to 53 GPa by synchrotron X-ray diffraction in a diamond anvil cell. Further, the graphane II phase did not react with the surrounding hydrogen and did not undergo any phase transformations upon the compression and after heating to 1500 K at 53GPa. The high thermal and pressure stability of this exotic phase makes it an important part of the C–H system. The obtained g(E) spectra of graphite hydride and deuteride were used to calculate temperature dependences of their heat capacity. Measurements of the heat capacity at temperatures 120–673 K confirmed the good accuracy of these calculations.

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

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2025

Long-Range Interaction between Molecules with Electronic Degeneracy: Beyond the Born–Oppenheimer Approximation

Long-range intermolecular interaction plays an essential role in (ultra)cold collisions. For monomers with degenerate electronic states, an accurate description of the long-range interaction can be challenging due to the potential breakdown of the Born–Oppenheimer approximation. In this work, we developed a general method to construct long-range interaction potential energy surfaces (IPESs) between arbitrary molecules using the degenerate perturbation theory of Schrödinger. To further simplify the application of this method, a group-theoretical method is introduced to identify nonzero components of electrostatic properties of monomers, especially for the nonlinear molecules. To demonstrate this approach, we report the long-range IPESs of the O( 3 P)–OH( 2 Π) and OH( 2 Π)–OH( 2 Π) systems. The excellent agreement with ab initio calculations in large separations, even near electronic degeneracies, confirms the accuracy of our approach. The method is expected to lay the groundwork for the accurate investigation of (ultra)cold collisional dynamics involving open-shell species.

ab initio calculations↗