Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio modeling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

First β -delayed neutron spectroscopy of 24 O

The β decay of 24 O was investigated at NSCL/MSU using a combination of neutron and γ-ray spectroscopy. For the first time, the β-delayed neutron energy spectrum of 24 O was measured, revealing three intensely populated, isolated neutron-unbound states in 24 F. This allowed for the extraction of the decay strength in 24 F up to 6.2 MeV. A comprehensive comparison of the experimental results with various nuclear theories, ranging from the empirical shell model to the most advanced ab initio calculations, was conducted. While most theoretical predictions align with the experimental data for low-lying states, discrepancies arise at higher excitation energies. Finally, in the transition from 24 O to 24 F, shell model calculations using the empirical USDB interaction predicted the structure of both nuclei without invoking the need for a stronger proton-neutron tensor force, which was postulated for the neighboring isotone 25 F.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Toward Realistic Models of the Electrocatalytic Oxygen Evolution Reaction

The electrocatalytic oxygen evolution reaction (OER) supplies the protons and electrons needed to transform renewable electricity into chemicals and fuels. However, the OER is kinetically sluggish; it operates at significant rates only when the applied potential far exceeds the reversible voltage. The origin of this overpotential is hidden in a complex mechanism involving multiple electron transfers and chemical bond making/breaking steps. Our desire to improve catalytic performance has then made mechanistic studies of the OER an area of major scientific inquiry, though the complexity of the reaction has made understanding difficult. While historically, mechanistic studies have relied solely on experiment and phenomenological models, over the past twenty years ab initio simulation has been playing an increasingly important role in developing our understanding of the electrocatalytic OER and its reaction mechanisms. In this Review we cover advances in our mechanistic understanding of the OER, organized by increasing complexity in the way through which the OER is modeled. We begin with phenomenological models built using experimental data before reviewing early efforts to incorporate ab initio methods into mechanistic studies. We go on to cover how the assumptions in these early ab initio simulations-no electric field, electrolyte, or explicit kinetics-have been relaxed. Through comparison with experimental literature, we explore the veracity of these different assumptions. We summarize by discussing the most critical open challenges in developing models to understand the mechanisms of the OER.

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

A density-functional theory study of the Al/AlOx/Al tunnel junction

The aluminum oxide tunnel junction is a key component of the majority of superconducting quantum devices. For high-quality, reproducible, and scalably manufacturable qubits, the ability to fabricate Josephson junctions (JJs) with a targeted critical current and high uniformity is essential. In this study, we use first-principles modeling to assess fundamental aspects of the atomic structure of both amorphous and crystalline aluminum oxide tunnel junctions and relate the structure to predicted performance metrics. We use modified ab initio molecular dynamics to develop realistic models of the tunnel junction, from which interface roughness and local thickness fluctuations are analyzed in an unbiased manner by training a neural network to identify the boundary between metal and oxide. We show that the effective thickness of the insulating part of the junction can be different from the apparent physical thickness. We calculate the rate of Cooper pair tunneling for the atomically resolved electrostatic potential using direct numerical solution in 3D, which shows a channeling effect that impacts the junction critical current. The predicted critical current is a useful JJ design parameter that can be accessed from the ab initio calculations without fitting parameters. To assess the limits of uniformity and fabrication choices (e.g., oxidation vs epitaxy), we compare the amorphous junctions to crystalline models, which show order of magnitude more efficient tunneling compared to the amorphous case, underlining the connection between atomistic structure and Cooper pair tunneling efficiency. Further, this work provides a foundation for ab initio materials design and evaluation to help accelerate future development of improved tunnel junctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Observing and Modeling the Sequential Pairwise Reactions that Drive Solid‐State Ceramic Synthesis

Abstract Solid‐state synthesis from powder precursors is the primary processing route to advanced multicomponent ceramic materials. Designing reaction conditions and precursors for ceramic synthesis can be a laborious, trial‐and‐error process, as heterogeneous mixtures of precursors often evolve through a complicated series of reaction intermediates. Here, ab initio thermodynamics is used to model which pair of precursors has the most reactive interface, enabling the understanding and anticipation of which non‐equilibrium intermediates form in the early stages of a solid‐state reaction. In situ X‐ray diffraction and in situ electron microscopy are then used to observe how these initial intermediates influence phase evolution in the synthesis of the classic high‐temperature superconductor YBa 2 Cu 3 O 6+ x (YBCO). The model developed herein rationalizes how the replacement of the traditional BaCO 3 precursor with BaO 2 redirects phase evolution through a low‐temperature eutectic melt, facilitating the formation of YBCO in 30 min instead of 12+ h. Precursor selection plays an important role in tuning the thermodynamics of interfacial reactions and emerges as an important design parameter in planning kinetically favorable synthesis pathways to complex ceramic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

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

Mirror symmetry breaking and lateral stacking shifts in twisted trilayer graphene

Here, we construct a continuum model of twisted trilayer graphene using ab initio density-functional-theory calculations, and we apply it to address twisted trilayer electronic structure. Our model accounts for moiré variation in site energies, hopping between outside layers, and hopping within layers. We focus on the role of a mirror symmetry present in ABA graphene trilayers with a middle layer twist. The mirror symmetry is lost intentionally when a displacement field is applied between layers, and unintentionally when the top layer is shifted laterally relative to the bottom layer. We calculate two band-structure characteristics that are directly relevant to transport measurements, namely the Drude weight and the weak-field Hall conductivity. We relate them via the Hall density to assess the influence of the accidental lateral stacking shifts currently present in all experimental devices on electronic properties. Finally, we comment on the role of the possible importance of accidental lateral stacking shifts for superconductivity in twisted trilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Imaging moiré flat bands in three-dimensional reconstructed WSe 2 /WS 2 superlattices

Moiré superlattices in transition metal dichalcogenide (TMD) heterostructures can host novel correlated quantum phenomena due to the interplay of narrow moiré flat bands and strong, long-range Coulomb interactions. However, microscopic knowledge of the atomically reconstructed moiré superlattice and resulting flat bands is still lacking, which is critical for fundamental understanding and control of the correlated moiré phenomena. Here we quantitatively study the moiré flat bands in three-dimensional (3D) reconstructed WSe 2 /WS 2 moiré superlattices by comparing scanning tunnelling spectroscopy (STS) of high-quality exfoliated TMD heterostructure devices with ab initio simulations of TMD moiré superlattices. A strong 3D buckling reconstruction accompanied by large in-plane strain redistribution is identified in our WSe 2 /WS 2 moiré heterostructures. STS imaging demonstrates that this results in a remarkably narrow and highly localized K-point moiré flat band at the valence band edge of the heterostructure. A series of moiré flat bands are observed at different energies that exhibit varying degrees of localization. Our observations contradict previous simplified theoretical models but agree quantitatively with ab initio simulations that fully capture the 3D structural reconstruction. Furthermore, our results reveal that the strain redistribution and 3D buckling in TMD heterostructures dominate the effective moiré potential and the corresponding moiré flat bands at the Brillouin zone K points.

36 MATERIALS SCIENCE↗
DOE OSTI2022

SS-HORSE extension of the no-core shell model: Application to resonances in He 7

Theoretical ab initio studies of resonances in the unbound 7 He nucleus are presented. We perform no-core shell-model calculations with NN interactions Daejeon16 and JISP16 and utilize the SS-HORSE method to calculate the S matrix for two-body channels n + 6 He and n + 6 He* with 6He respectively in the ground and excited 2 + states as well as for the four-body democratic decay channel 4 He + n + n + n. The resonant energies and widths are obtained by numerical location of the S-matrix poles. Here, we describe all experimentally known 7 He resonances and suggest an interpretation of an observed wide resonance of unknown spin-parity.

6 ≤ A ≤ 19↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2022

General, Rigorous Approach for the Treatment of Interfragment Covalent Bonds

Here, a generalized, projection-based transformation of the method-agnostic Fock operator in various ab initio fragment-based quantum chemistry methods has been developed for the treatment of interfragment covalent bonds. This transformation freezes the relevant localized molecular orbital associated with each interfragment bond, thereby restricting the variational subspace of the fragment wave functions, in order to maintain the proper physical characteristics of the involved covalent bonds. In addition, sets of orbitals that would lead to multiple occupancy of certain orbitals are explicitly removed from the variational space. The transformation is developed for the specific case of mutually orthonormal frozen and unfrozen orbitals within each fragment. The newly developed approach is then used to study model systems with two popular ab initio fragment-based methods, and the results of these calculations are compared to those obtained by existing methodologies. Analysis is focused on both quantitative and qualitative accuracy as well as computational scalability and stability. Other methods for which the developed formalisms are appropriate are outlined, and future extensions of the methods are discussed.

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

Efficacy of the symmetry-adapted basis for ab initio nucleon-nucleus interactions for light- and intermediate-mass nuclei

We study the efficacy of a new ab initio framework that combines the symmetry-adapted (SA) no-core shell-model approach with the resonating group method (RGM) for unified descriptions of nuclear structure and reactions. We obtain ab initio neutron-nucleus interactions for 4 He, 16 O, and 20 Ne targets, starting with realistic nucleon-nucleon potentials. We discuss the effect of increasing model space sizes and symmetry-based selections on the SA-RGM norm and direct potential kernels, as well as on phase shifts, which are the input to calculations of cross sections. We demonstrate the efficacy of the SA basis and its scalability with particle numbers and model space dimensions, with a view toward ab initio descriptions of nucleon scattering and capture reactions up through the medium-mass region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Simple and efficient algorithms for training machine learning potentials to force data

Machine learning models, trained on data from ab initio quantum simulations, are yielding molecular dynamics potentials with unprecedented accuracy. One limiting factor is the quantity of available training data, which can be expensive to obtain. A quantum simulation often provides all atomic forces, in addition to the total energy of the system. These forces provide much more information than the energy alone. It may appear that training a model to this large quantity of force data would introduce significant computational costs. Actually, training to all available force data should only be a few times more expensive than training to energies alone. Here, we present a new algorithm for efficient force training, and benchmark its accuracy by training to forces from real-world datasets for organic chemistry and bulk aluminum.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

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

Unraveling Local Structure of Molten Salts via X-ray Scattering, Raman Spectroscopy, and Ab Initio Molecular Dynamics

In this work, we resolve a long-standing issue concerning the local structure of molten MgCl 2 by employing a multimodal approach, including X-ray scattering and Raman spectroscopy, along with the theoretical modeling of the experimental spectra based on ab initio molecular dynamics (AIMD) simulations utilizing several density functional theory (DFT) methods. We demonstrate the reliability of AIMD simulations in achieving excellent agreement between the experimental and simulated spectra for MgCl 2 and 50 mol % MgCl 2 + 50 mol % KCl, and ZnCl 2 , thus allowing structural insights not directly available from experiment alone. A thorough computational analysis using five DFT methods provides a convergent view that octahedrally coordinated magnesium in pure MgCl 2 upon melting preferentially coordinates with five chloride anions to form distorted square pyramidal polyhedra that are connected via corners and to a lesser degree via edges. This is contrasted with the results for ZnCl 2 , which does not change its tetrahedral coordination on melting. Although the five-coordinate MgCl 5 3– complex was not considered in the early literature, together with an increasing tendency to form a tetrahedrally coordinated complex with decreasing the MgCl 2 content in the mixture with alkali metal chloride systems, current work reconciles the results of most previous seemingly contradictory experimental studies.

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

Derivation of low-energy Hamiltonians for heavy-fermion materials

Here, by utilizing a multiorbital periodic Anderson model with parameters obtained from ab initio band structure calculations, combined with degenerate perturbation theory, we derive effective Kondo-Heisenberg and spin Hamiltonians that capture the interaction among the effective magnetic moments. This derivation encompasses fluctuations via both nonmagnetic 4⁢𝑓 0 and magnetic 4⁢𝑓 2 virtual states, and its accuracy is confirmed through comparison with experimental data obtained from CeIn 3 . The significant agreement observed between experimental results and theoretical predictions underscores the potential of deriving minimal models from first-principles calculations for achieving a quantitative description of 4⁢𝑓 materials. Moreover, our microscopic derivation unveils the underlying origin of anisotropy in the exchange interaction between Kramers doublets, shedding light on the conditions under which this anisotropy may be weak compared to the isotropic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Partial-wave decomposition of the Keldysh ionization amplitude

We present an alternative way of calculating the Keldysh amplitude, i.e., the length-gauge form of the ionization amplitude in the strong-field approximation. The amplitude is evaluated exactly by expanding it in Fourier components and partial waves. Comparisons of the semianalytic model predictions with results of ab initio numerical simulations of the time-dependent Schrödinger equation for the interaction of electrons in short-range potentials with intense laser light yield excellent agreement, for wavelengths from the single photon to the multiphoton to the tunneling regime. Specifically, for ionization from initial states with higher angular momentum quantum number, e.g., p states, a significant improvement over predictions based on the popular saddle-point approximation is found. Furthermore, the current model rate allows for interpretation of the strong-field ionization process in terms of multiphoton absorption pathways and angular momentum selection rules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Neutrino Flavor Model Building and the Origins of Flavor and CP Violation

The neutrino sector offers one of the most sensitive probes of new physics beyond the Standard Model of Particle Physics (SM). The mechanism of neutrino mass generation is still unknown. The observed suppression of neutrino masses hints at a large scale, conceivably of the order of the scale of a rand unified theory (GUT), which is a unique feature of neutrinos that is not shared by the charged fermions. The origin of neutrino masses and mixing is part of the outstanding puzzle of fermion masses and mixings, which is not explained ab initio in the SM. Flavor model building for both quark and lepton sectors is important in order to gain a better understanding of the origin of the structure of mass hierarchy and flavor mixing, which constitute the dominant fraction of the SM parameters. Recent activities in neutrino flavor model building based on non-Abelian discrete flavor symmetries and modular flavor symmetries have been shown to be a promising direction to explore. The emerging models provide a framework that has a significantly reduced number of undetermined parameters in the flavor sector. In addition, such a framework affords a novel origin of CP violation from group theory due to the intimate connection between physical CP transformation and group theoretical properties of non-Abelian discrete groups. Model building based on non-Abelian discrete flavor symmetries and their modular variants enables the particle physics community to interpret the current and anticipated upcoming data from neutrino experiments. Non-Abelian discrete flavor symmetries and their modular variants can result from compactification of a higher-dimensional theory. Pursuit of flavor model building based on such frameworks thus also provides the connection to possible UV completions: in particular, to string theory. We emphasize the importance of constructing models in which the uncertainties of theoretical predictions are smaller than, or at most compatible with, the error bars of measurements in neutrino experiments. While there exist proof-of-principle versions of bottom-up models in which the theoretical uncertainties are under control, it is remarkable that the key ingredients of such constructions were discovered first in top-down model building. We outline how a successful unification of bottom-up and top-down ideas and techniques may guide us towards a new era of precision flavor model building in which future experimental results can give us crucial insights into the UV completion of the SM.

Almumin, Yahya (ORCID:0000000196396772)↗