Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evaluation of phase shifts for nonrelativistic elastic scattering using quantum computers

Simulations of scattering processes are essential in understanding the physics of our universe. Computing relevant scattering quantities from ab initio methods is extremely difficult on classical devices because of the substantial computational resources needed. Here, this work reports the development of an algorithm that makes it possible to obtain phase shifts for generic nonrelativistic elastic scattering processes on a quantum computer. This algorithm is based on extracting phase shifts from the direct implementation of the real-time evolution. The algorithm is improved by a variational procedure making it more accurate and resistant to the quantum noise. The reliability of the algorithm is first demonstrated by means of classical numerical simulations for different potentials and later tested on existing quantum hardware, specifically on IBM quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Few-nucleon and many-nucleon systems with semilocal coordinate-space regularized chiral nucleon-nucleon forces

We employ a variety of ab initio methods, including Faddeev-Yakubovsky equations, no-core configuration interaction approach, coupled-cluster theory, and in-medium similarity renormalization group, to perform a comprehensive analysis of the nucleon-deuteron elastic and breakup reactions and selected properties of light and medium-mass nuclei up to 48 C a using the recently constructed semilocal coordinate-space regularized chiral nucleon-nucleon potentials. We compare the results with those based on selected phenomenological and chiral EFT two-nucleon potentials, discuss the convergence pattern of the chiral expansion, and estimate the achievable theoretical accuracy at various chiral orders using an approach to quantify truncation errors of the chiral expansion without relying on cutoff variation. We also address the robustness of this method and explore alternative ways to estimate the theoretical uncertainty from the truncation of the chiral expansion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Impact of Modular Total Absorption Spectrometer measurements of β decay of fission products on the decay heat and reactor ν ¯ e flux calculation

Here, we report the results of a β-decay study of fission products 86 Br, 89 Kr, 89 Rb, 90 gs Rb, 90m Rb, 90 Kr, 92 Rb, 139 Xe, and 142 Cs performed with the Modular Total Absorption Spectrometer (MTAS) and on-line mass-separated ion beams. These radioactivities were assessed by the Nuclear Energy Agency as having high priority for decay heat analysis during a nuclear fuel cycle. We observe a substantial increase in β feeding to high excited states in all daughter isotopes in comparison to earlier data. This increases the average γ-ray energy emitted by the decay of fission fragments during the first 10 000 s after fission of 235 U and 239 Pu by approximately 2% and 1%, respectively, improving agreement between results of calculations and direct observations. New MTAS results reduce the reference reactor ¯ν e flux used to analyze reactor ¯ν e interaction with detector matter. The reduction determined by the ab initio method for the four nuclear fuel components, 235 U, 238 U, 239 Pu, and 241 Pu, amounts to 0.976, 0.986, 0.983, and 0.984, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab Initio Simulation of Amorphous Graphite

An amorphous graphite material has been predicted from molecular dynamics simulation using ab initio methods. Carbon materials reveal a strong proclivity to convert into a sp 2 network and then layer at temperatures near 3000 K within a density range of ca. 2.2–2.8 g/cm 3 . Each layer of amorphous graphite is a monolayer of amorphous graphene including pentagons and heptagons in addition to hexagons, and the planes are separated by about 3.1 Å. The layering transition has been studied using various structural and dynamical analyses. The transition is unique as one of partial ordering (long range order of planes and galleries, but topological disorder in the planes). The planes are quite flat, even though monolayer amorphous graphene puckers near pentagonal sites. Interplane cohesion is due partly to non-Van der Waals interactions. The structural disorder has been studied closely, especially the consequences of disorder to electronic transport. Finally, it is expected that the transition elucidated here may be salient to other layered materials.

42 ENGINEERING↗
DOE OSTI2022

Learning Grain-Boundary Segregation: From First Principles to Polycrystals

The segregation of solute atoms at grain boundaries (GBs) can strongly impact the structural and functional properties of polycrystals. Yet, due to the limited availability of simulation tools to study polycrystals at the atomistic scale (i.e., interatomic potentials), there is a minimal understanding of the variation of solute segregation tendencies across the very complex space of GB microenvironments and the large range of alloys in which it can occur. Here, we develop an algorithmic framework that can directly learn the full spectrum of segregation energies for a metal solute atom in a metal polycrystal from ab initio methods, bypassing the need for alloy interatomic potentials. This framework offers a pathway to a comprehensive catalog of GB solute segregation with quantum accuracy, for the entire alloy space. As an initial demonstration in this pursuit, we build an extensive GB segregation database for aluminum-based alloys across the periodic table, including dozens of alloys for which there are substantially no prior data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Intermolecular coupling and superconductivity in PbMo 6 S 8 and other Chevrel phase compounds

To understand superconductivity in Chevrel phase compounds and guide the search for interesting properties in materials created with Chevrel phase molecules as building blocks, we use ab initio methods to study the properties of single Mo 6 X 8 molecules with X=S, Se, Te as well as the bulk solid PbMo 6 S 8 . In bulk PbMo 6 S 8 , the different energy scales from strong to weak are the band kinetic energy, the intramolecular Coulomb interaction, the on-molecule Jahn-Teller energy, and the Hund's exchange coupling. The metallic state is stable with respect to Mott and polaronic insulating states. The bulk compound is characterized by a strong electron-phonon interaction with the largest coupling involving phonon modes with energies in the range from 11 to 17 meV and with a strong intermolecule (Peierls) character. A two-band Eliashberg equation analysis shows that the superconductivity is strong coupling, with different gaps on the two Fermi surface sheets. A Bergman-Rainer analysis of the functional derivative of the transition temperature with respect to the electron-phonon coupling reveals that the Peierls modes provide the most important contribution to the superconductivity. This work illustrates the importance of intermolecular coupling for collective phenomena in molecular solids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory

We extend the spherical coupled-cluster ab initio method for open-shell nuclei where two nucleons are removed from a shell subclosure. Following the recent implementation of the two-particle-attached approach [Phys. Rev. C 110, 044306 (2024)], we focus on the two-particle-removed method. Using the equations-of-motion framework, we address both nuclear structure and dipole response functions by coupling coupled-cluster theory with the Lorentz integral transform technique. We perform calculations using chiral interactions, including three-nucleon forces, and estimate many-body uncertainties by comparing different coupled-cluster truncation schemes. Here, we validate our approach by studying ground-state energies, excited states, and electric dipole polarizabilities in the oxygen and calcium isotopic chains. For binding energies and selected low-lying excited states, we achieve an accuracy comparable to that of the established closed-shell coupled-cluster theory and generally agree with experiment. Finally, we underestimate experimental data for electric dipole polarizabilities, particularly in calcium isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Grain-Boundary Shear-Migration Coupling in Al Bicrystals. Atomistic Modeling

The energy of grain boundary shears is calculated for symmetric grain boundaries (GBs) using ab initio methods and molecular-dynamic modeling in order to elucidate mechanisms that control GB shear-migration coupling in typical symmetric GBs, such as Σ3 (111), Σ5 (012), Σ5 (013) and Σ11 (113) tilt GBs, in Al bicrystal. The energy of generalized grain-boundary stacking faults (GB–SF) is determined, and the preferred directions and the energy barrier are established for grain-boundary slippage. It is shown that the relative slippage of neighboring grains at certain directions of particle shears is accompanied by conservative migration of GB in the direction perpendicular to its plain. The modeling data are comparative to known grain-boundary shear-migration coupling mechanisms in Al.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Mechanical and Dynamic Stability of Complete and Nonstoichiometric 3C-Si{sub x}C{sub y} from Ab Initio Calculations

The energies of formation of vacancies in the carbon and silicon sublattices, the independent elastic constants, the all-round compression, shear and Young’s moduli, and the anisotropy coefficients are determined for the complete and nonstoichiometric cubic phases of 3C-Si{sub x}C{sub y} (x, y = 1.0–0.75) by ab initio methods of the band theory. In the formalism of the density functional perturbation theory (DFPT), the phonon dispersion dependences are obtained for these phases (the comparison with the experiment is given for the complete phase). It is shown that the mechanical characteristics of the phases become strongly anisotropic upon the transition from 3C-SiC{sub 0.875} to 3C-SiC{sub 0.75}. It is established from the analysis of the phonon dispersion curves that the 3C-SiC{sub 0.875} and 3C-SiC{sub 0.75} phases, in contrast to the complete 3C-SiC phase, are dynamically unstable at T = 0 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Microscopic Description of the Mechanism of Transition between the 2H and 4H Polytypes of Silicon Carbide

Abstract—: The mechanism of displacement of one close-packed SiC layer from one minimum position to another on the example of SiC polytype transition 2H → 4H has been studied by ab initio methods. It has been shown that the intermediate state with monoclinic symmetry Cm greatly facilitates this displacement breaking it into two stages. Initially, the Si atom chiefly moves, only then—mainly the C atom. In this case, the Si–C bond is significantly tilted in comparison with the initial position, which allows the reducing of the compression of the SiC bonds in the (11 2-bar 0) plane. Two transition states of this process, which also possess the Cm symmetry, have been computed. It has been found that the height of the activation barrier of the process of moving the close-packed layer of SiC from one position to another is equal to 1.8 eV. The energy profile of this movement has been calculated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Techniques for Polytypic Transformations in Silicon Carbide

Abstract—: Two main polytype transformations in silicon carbide, namely, 2H → 6H and 3C → 6H, have been studied by ab initio methods. It has been shown that the intermediate phases with trigonal symmetry P3m1 and monoclinic symmetry Cm make it much easier to move the close-packed layers in such transitions by breaking them up into separate stages. It has been found that these two polytype transformations proceed in completely different ways. The links being relocated noticeably tilted compared to their initial position at the transition 2H → 6H, which allows the compression of the SiC links in the plane (11 2-bar 0). The transition 3C → 6H is carried out through the formation of Si–Si and C–C auxiliary links, living for a short time and helping densely packed layers to swap places. As a result, the activation barrier of the transformation 2H → 6H (1.7 eV/atom) is significantly less than the activation barrier of the transformation 3C → 6H (3.6 eV/atom), which means that the second transition should occur at the temperatures by 750–800°C higher than the first one. The energy profiles of this polytypic transformations, as well as the geometry of all intermediate and transition phases have been calculated. It has been shown that all the transition states have monoclinic symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

The Effect of Hydrogen on Fluctuation Embrittlement of Aluminum

In this paper, we describe the main processes occurring during the generation of vacancies in aluminum in the presence of hydrogen by ab initio methods using the SCAN functional. Hydrogen is shown to reduce the generation energy of vacancies from 2.8 to 0.8 eV. In this case, the eight hydrogen atoms located in the tetrahedral lattice voids around single aluminum atom greatly facilitate its movement to the interstitial site. The dependence of the activation energy of hydrogen embrittlement of aluminum on the hydrogen concentration in aluminum and temperature is calculated based on the kinetic strength concept. Hydrogen is shown to reduce the time of aluminum destruction only if its mole fraction in aluminum is greater than the critical level (~ 3 × 10{sup –4} at T = 293 K).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A Theoretical Analysis on the Dynamic Behavior and Electromechanical Response of PZT

Piezoelectric materials are used as a power source in stockpile components for safety and reliability assurances. The objective of this project is to gain insights into the fundamental pressure-induced behavior and electromechanical response of lead zirconate titanate (PZT). Specifically, to establish the basis for an accurate, physics-based material model for robust model-based design to rapidly optimize PZT-based materials and components for performance studies. Established ab-initio methods are used to interrogate and understand the dynamic behavior of PZT as a function of composition (50/50, 65/35, 80/20, 95/5) and dopants (La, Nb) to overcome the costly and time-consuming experimental methodologies. New cold curves for pure and doped single-crystal PZT are obtained as the reference equation-of-state (EOS). A negligible change in the pressure responses was observed for the systems and strains studied. Dielectric and piezoelectric responses of pure and doped single-crystal PZT were also calculated as a function of pressure. For undoped PZT, there is a clear and orderly pressure and composition dependence. For doped PZT, there is a significant increase in the responses, but the behavior is very disordered and inconclusive.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Studying electroweak few-body observables in chiral effective field theory

The use of nuclei to study electroweak probes is becoming increasingly relevant experimentally. The success of dark matter and neutrino experiments strongly depends on the ability to control nuclear effects in order to extract the fundamental parameters associated with external probes. Therefore, reliable theoretical calculations of nuclear structure and reactions, with well-controlled errors, are crucial for the success of experimental efforts. Currently, chiral effective field theory ($\chi$EFT) coupled with {\it ab-initio} methods represents one of the best approaches that fulfills these requirements. To use this approach as a tool for studying fundamental physics, it is essential to validate it against experimental data for which the calculations are well under control, such as the elastic scattering of electrons on nuclei. In this proceeding, I will present recent developments in the fitting of electromagnetic currents derived using $\chi$EFT and the calculation of electromagnetic form factors of light nuclei. The results of these calculations demonstrate the strength of the theory in describing the interaction of nuclei with electromagnetic probes over a broad range of momentum transfers and highlight the robustness of $\chi$EFT for analyzing future experimental data aimed at extracting fundamental parameters.

Gnech, Alex [Old Dominion Univ., Norfolk, VA (Unit↗
DOE OSTI2021

Quantum Calculations of VX Ammonolysis and Hydrolysis Pathways via Hydrated Lithium Nitride

Recently, lithium nitride (Li 3 N) has been proposed as a chemical warfare agent (CWA) neutralization reagent for its ability to produce nucleophilic ammonia molecules and hydroxide ions in aqueous solution. Quantum chemical calculations can provide insight into the Li 3 N neutralization process that has been studied experimentally. Here, we calculate reaction-free energies associated with the Li3N-based neutralization of the CWA VX using quantum chemical density functional theory and ab initio methods. We find that alkaline hydrolysis is more favorable to either ammonolysis or neutral hydrolysis for initial P-S and P-O bond cleavages. Reaction-free energies of subsequent reactions are calculated to determine the full reaction pathway. Notably, products predicted from favorable reactions have been identified in previous experiments.

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

How Accurate Can Crystal Structure Predictions Be for High-Energy Molecular Crystals?

Molecular crystals have shallow potential energy landscapes, with multiple local minima separated by very small differences in total energy. Predicting molecular packing and molecular conformation in the crystal generally requires ab initio methods of high accuracy, especially when polymorphs are involved. We used dispersion-corrected density functional theory (DFT-D) to assess the capabilities of an evolutionary algorithm (EA) for the crystal structure prediction (CSP) of well-known but challenging high-energy molecular crystals (HMX, RDX, CL-20, and FOX-7). While providing the EA with the experimental conformation of the molecule quickly re-discovers the experimental packing, it is more realistic to start instead from a naïve, flat, or neutral initial conformation, which reflects the limited experimental knowledge we generally have in the computational design of molecular crystals. By doing so, and using fully flexible molecules in fully variable unit cells, we show that the experimental structures can be predicted in fewer than 20 generations. Nonetheless, one must be aware that some molecular crystals have naturally hindered evolutions, requiring as many attempts as there are space groups of interest to predict their structures, and some may require the accuracy of all-electron calculations to discriminate between closely ranked structures. To save resources in this computationally demanding process, we showed that a hybrid xTB/DFT-D approach could be considered in a subsequent study to push the limits of CSP beyond 200+ atoms and for cocrystals.

42 ENGINEERING↗
DOE OSTI2017

Pressure Distortion of the H{sub 2}–He Collision-induced Absorption at the Photosphere of Cool White Dwarf Stars

Collision-induced absorption (CIA) from molecular hydrogen is a dominant opacity source in the atmosphere of cool white dwarfs. It results in a significant flux depletion in the near-IR and IR parts of their spectra. Because of the extreme conditions of helium-rich atmospheres (where the density can be as high as a few g cm{sup −3}), this opacity source is expected to undergo strong pressure distortion and the currently used opacities have not been validated at such extreme conditions. To check the distortion of the CIA opacity, we applied state-of-the-art ab initio methods of computational quantum chemistry to simulate the CIA opacity at high densities. The results show that the CIA profiles are significantly distorted above densities of 0.1 g cm{sup −3} in a way that is not captured by the existing models. The roto-translational band is enhanced and shifted to higher frequencies as an effect of the decrease of the interatomic separation of the H{sub 2} molecule. The vibrational band is blueward shifted and split into Q {sub R} and Q {sub P} branches, separated by a pronounced interference dip. Its intensity is also substantially reduced. The distortions result in a shift of the maximum of the absorption from 2.3 μm to 3–7 μm, which could potentially explain the spectra of some very cool, helium-rich white dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

Radiative lifetimes of the second negative system of O2(+)

The dipole moment of the A2Pi(u)-X2Pi(g) transition of O2(+) is calculated as a function of internuclear distance using ab initio methods. The band absorption oscillator strengths and band transition probabilities of the second negative system are derived and the resulting lifetimes are compared with experimental data. The high-lying v double prime levels of the ground state may decay into low-lying v prime levels of the excited state. The corresponding radiative lifetimes are calculated.

Wetmore, R. W.↗