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

High–Capacity Splitting of Mono– and Dibranched Hexane Isomers by a Robust Zinc–Based Metal–Organic Framework

High-efficiency separation of C 6 alkanes, particularly the mono- and dibranched isomers by using porous solids, is of paramount significance in the petrochemical industry and, remains a daunting challenge. In this work, we report the complete separation of linear/monobranched hexanes from their dibranched isomers through selective size-exclusion by a microporous MOF, Zn-tcpt (H3tcpt=2,4,6-tris(4-carboxyphenoxy)-1,3,5-triazine), with a two-fold interpenetrated structure of hms nets. Importantly, its adsorption capacity and selectivity are notably higher than those of the previously reported adsorbents that can split mono- and dibranched alkane isomers. Dynamic breakthrough measurements verify the excellent separation of C 6 alkane isomers by Zn-tcpt, and the size-exclusion based separation mechanism has been confirmed by ab initio materials modeling. Furthermore, the high-efficiency separation of alkane isomers by Zn-tcpt can be attributed to its optimal pore dimensions as well as high porosity.

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

Structure/property relationship of semicrystalline polymers during tensile deformation: a molecular dynamics approach

Molecular dynamics (MD) simulation with an ab initio numerical model is conducted to study the structure–property relationship of semicrystalline polymers during uniaxial deformation. The effects of chain length, temperature, and strain rate on mechanical properties are discussed. The influences of microstructural evolution such as bond length, entanglement density, and chain orientation are also studied quantitatively. The temperature will greatly affect the chain conformation in amorphous domains, and the results revealed that the interaction of amorphous and crystalline domains played a crucial role during stretching. Accordingly, the yielding of semicrystalline polymers follows different mechanisms at temperatures above and below T g . A melt-recrystallization scheme is observed during yielding at higher temperatures, while destruction of crystal structures is observed at lower temperatures at the yield point. The correlated effects of different temperatures and strain rates on mechanical properties are examined. This research is part of our efforts to develop a digital twin of a real experiment for efficient optimization of polymer material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

ElasTool v3.0: Efficient computational and visualization toolkit for elastic and mechanical properties of materials

Efficient computation and visualization of elastic and mechanical properties are crucial in the selection of materials and the design of new materials. Here, the ElasTool v3.0 toolkit marks a significant advancement in the computational analysis and visualization of elastic and mechanical properties of materials, essential in material selection and design. This enhanced version extends beyond standard calculations like elastic tensor, Young's modulus, bulk modulus, and Poisson's ratio. It introduces capabilities for computing minimum thermal conductivity, linear compressibility, rendering the Christoffel equation, and elastic energy density. Notably, it integrates advanced visualization tools, including compatibility with Plotly and Elate web platforms for interactive web-based property exploration. A key feature of ElasTool v3.0 is the implementation of second-order elastic constants (SOECs) for tubular 2D-based nanostructures and nanotubes. Leveraging high-efficiency strain-matrix sets (OHESS), the toolkit now facilitates efficient computation of elastic constants and mechanical properties at both zero and finite temperatures for 1D, 2D, and 3D dimensions. ElasTool is openly accessible on GitHub: https://github.com/gmp007/elastool.

1D, 2D, 3D, and tubular 2D nanostructure and nanot↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗
DOE OSTI2022

In Situ Determination of Speciation and Local Structure of NaCl–SrCl 2 and LiF–ZrF 4 Molten Salts

Understanding the local environment of the metal atoms in salt melts is important for modeling the properties of melts and predicting their behavior and thus helping enable the development of technologies such as molten salt reactors and solar-thermal power systems and new approaches to recycling rare-earth metals. Toward that end, we have developed an in situ approach for measuring the coordination of metals in molten salt coupling X-ray absorption spectroscopy (XAS) and Raman spectroscopy. Our approach was demonstrated for two salt mixtures (1.9 and 5 mol % SrCl 2 in NaCl, 0.8 and 5 mol % ZrF 4 in LiF) at up to 1100 °C. Near-edge (X-ray absorption near-edge structure, XANES) and extended X-ray absorption fine structure (EXAFS) spectra were measured. The EXAFS response was modeled using ab initio FEFF calculations. Strontium’s first shell is observed to be coordinated with chlorine (Sr 2+ –Cl – ) and zirconium’s first shell is coordinated by fluorine (Zr 4+ –F – ), both having coordination numbers that decrease with increasing temperature. Multiple zirconium complexes are believed to be present in the melt, which may interfere and distort the EXAFS spectra and result in an anomalously low zirconium first shell coordination number. Finally, the use of boron nitride (BN) powder as a salt diluent for XAFS measurements was found to not interfere with measurements and thus can be used for investigations of such systems.

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

X-ray Transient Absorption Studies of Exciton Self-Trapping

Here, we present optical pump/X-ray absorption probe studies to probe changes in electronic distribution and local structure upon formation of a self-trapped exciton. The experiments were carried out on the mixed-valence halide-bridged transition metal linear chain charge density wave material [Pt(en 2 )][Pt(en 2 )Cl 2 ]·(ClO 4 ) 4 . Photoinduced changes in the Pt L III XANES spectrum were measured on a picosecond time scale following generation of excitons via optical excitation of the intervalence charge transfer transition. Ab initio FEFF9 modeling of the spectra and the associated angular momentum projected density of states was used to interpret the response. We find photoinduced changes reflecting electronic redistribution and local lattice distortions that correspond to strong localization of the self-trapped state.

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

Surface Phase Stability of Fe 2 O 3 (001) in Hydrogen Reducing Environments: A DFT and XPS Analysis

Here, this study combines density functional theory (DFT) and ab initio thermodynamics calculations with X-ray photoelectron spectroscopy (XPS) investigations to identify the reduction properties of the Fe 2 O 3 (001) surface with implications for corrosion resistance, hydrogen transport, and energy safety. Ab initio thermodynamics modeling predicts fully hydroxylated surface stability across a broad range of pressures (1 × 10 –23 to 1 × 10 5 mbar) and temperatures below 700 K, consistent with previous experimental studies. Above 800 K, exposures to 1 × 10 –4 mbar H 2 , 1 × 10 –4 mbar O 2 , or 1 × 10 –4 mbar H 2 + 1 × 10 –4 mbar O 2 each yield unique XPS signals indicating a loss of −OH coverage, aligning with DFT predictions. Insight into the mechanism of reduction as a function of H 2 exposure is provided, as well as conditions that promote further reduction toward Fe 3 O 4 . Theoretical and experimental investigations indicate the ability to maintain the Fe 2 O 3 protective layer of iron oxides that have been exposed to H 2 environments by including trace amounts of aqueous O2.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Hot-Hole versus Hot-Electron Transport at Cu/GaN Heterojunction Interfaces

Among all plasmonic metals, copper (Cu) has the greatest potential for realizing optoelectronic and photochemical hot-carrier devices, thanks to its CMOS compatibility and outstanding catalytic properties. Yet, relative to gold (Au) or silver (Ag), Cu has rarely been studied and the fundamental properties of its photoexcited hot carriers are not well understood. Here, we demonstrate that Cu nanoantennas on p-type gallium nitride (p-GaN) enable hot-hole-driven photodetection across the visible spectrum. Importantly, we combine experimental measurements of the internal quantum efficiency (IQE) with ab initio theoretical modeling to clarify the competing roles of hot-carrier energy and mean-free path on the performance of hot-hole devices above and below the interband threshold of the metal. We also examine Cu-based plasmonic photodetectors on corresponding n-type GaN substrates that operate via the collection of hot electrons. By comparing hot hole and hot electron photodetectors that employ the same metal/semiconductor interface (Cu/GaN), we further elucidate the relative advantages and limitations of these complementary plasmonic systems. In particular, we find that harnessing hot holes with p-type semiconductors is a promising strategy for plasmon-driven photodetection across the visible and ultraviolet regimes. Given the technological relevance of Cu and the fundamental insights provided by our combined experimental and theoretical approach, we anticipate that our studies will have a broad impact on the design of hot-carrier optoelectronic devices and plasmon-driven photocatalytic systems.

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

Computational and Experimental Characterization of the Ligand Environment of a Ni-Oxo Catalyst Supported in the Metal–Organic Framework NU-1000

Heterogeneous catalysts exhibit significant changes in composition due to the influence of operating conditions, and these compositional changes can have dramatic effects on catalytic performance. For traditional bulk metal heterogeneous catalysts, relationships between composition and catalytic operating conditions are well documented. However, the influence of operating conditions on the compositions of single-site heterogeneous catalysts remains largely unresolved. To address this, we report a combined computational and experimental characterization of a Ni oxo catalyst under catalytic hydrogenation conditions. Specifically, pair distribution function (PDF) analysis is combined with ab initio thermodynamic modeling to investigate ligand environments present on a Ni oxo cluster supported in the metal–organic framework NU-1000. Comparisons of the experimentally observed and simulated Ni–O coordination numbers and Ni–O, Ni···Ni, and Ni···Zr distances provide insight into the Ni ligand environment under H 2 (g). These comparisons suggest significant OH and H 2 O content and, further, that different Ni ions within the cluster and/or NU-1000 structure may comprise subtly different numbers of these ligands. Further, the observation of significant H 2 O content under H 2 (g) suggests that the NU-1000 support supplies H 2 O to the cluster. Examples of ligand environments that could lead to the observed PDFs are provided. Furthermore, the combination of simulations and experiments provides new insights into the ligand environment for Ni-NU-1000 catalysts that will be useful for understanding the ligand environments of other single-site Ni catalysts as well.

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

Evaluating the evidence for exponential quantum advantage in ground-state quantum chemistry

Due to intense interest in the potential applications of quantum computing, it is critical to understand the basis for potential exponential quantum advantage in quantum chemistry. Here we gather the evidence for this case in the most common task in quantum chemistry, namely, ground-state energy estimation, for generic chemical problems where heuristic quantum state preparation might be assumed to be efficient. The availability of exponential quantum advantage then centers on whether features of the physical problem that enable efficient heuristic quantum state preparation also enable efficient solution by classical heuristics. Through numerical studies of quantum state preparation and empirical complexity analysis (including the error scaling) of classical heuristics, in both ab initio and model Hamiltonian settings, we conclude that evidence for such an exponential advantage across chemical space has yet to be found. While quantum computers may still prove useful for ground-state quantum chemistry through polynomial speedups, it may be prudent to assume exponential speedups are not generically available for this problem.

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

Machine learning for deep elastic strain engineering of semiconductor electronic band structure and effective mass

Abstract The controlled introduction of elastic strains is an appealing strategy for modulating the physical properties of semiconductor materials. With the recent discovery of large elastic deformation in nanoscale specimens as diverse as silicon and diamond, employing this strategy to improve device performance necessitates first-principles computations of the fundamental electronic band structure and target figures-of-merit, through the design of an optimal straining pathway. Such simulations, however, call for approaches that combine deep learning algorithms and physics of deformation with band structure calculations to custom-design electronic and optical properties. Motivated by this challenge, we present here details of a machine learning framework involving convolutional neural networks to represent the topology and curvature of band structures in k -space. These calculations enable us to identify ways in which the physical properties can be altered through “deep” elastic strain engineering up to a large fraction of the ideal strain. Algorithms capable of active learning and informed by the underlying physics were presented here for predicting the bandgap and the band structure. By training a surrogate model with ab initio computational data, our method can identify the most efficient strain energy pathway to realize physical property changes. The power of this method is further demonstrated with results from the prediction of strain states that influence the effective electron mass. We illustrate the applications of the method with specific results for diamonds, although the general deep learning technique presented here is potentially useful for optimizing the physical properties of a wide variety of semiconductor materials.

Tsymbalov, Evgenii↗
DOE OSTI2025

Extending Badger's rule. I. The relationship between energy and structure in hydrogen bonds

Here, we derive a new expression for the strength of a hydrogen bond (V HB ) in terms of the elongation of the covalent bond of the donor fragment participating in the hydrogen bond (Δr HB ) and the intermolecular coordinates R (separation between the heavy atoms) and θ (deviation of the hydrogen bond from linearity). The expression includes components describing the covalent D–H bond of the hydrogen bond donor via a Morse potential, the Pauli repulsion, and electrostatic interactions between the constituent fragments using a linear expansion of their dipole moment and a quadratic expansion of their polarizability tensor. We fitted the parameters of the model using ab initio electronic structure results for six hydrogen bonded dimers, namely, NH 3 –NH 3 , H 2 O–H 2 O, HF–HF, H 2 O–NH 3 , HF–H 2 O, and HF–NH 3 , and validated its performance for extended parts of their potential energy surfaces, resulting in a mean absolute error ranging from 0.07 to 0.31 kcal/mol. The derived expression describes the energy–structure relationship in terms of a single structural parameter, namely, the elongation of the donor’s covalent bond (Δr HB ), and suggests the novel relationship of 8.0 kcal/mol pm −1 (or 0.8 kcal/mol per 0.001 Å elongation). This structural parameter is easily obtained from theory and can serve as the single descriptor of the strength of individual hydrogen bonds.

Santis, Garrett D. [Univ. of Washington, Seattle, ↗
DOE OSTI2025

Effective field theories for collective excitations of atomic nuclei

Collective modes emerge as the relevant degrees of freedom that govern low-energy excitations of atomic nuclei. These modes—rotations, pairing rotations, and vibrations—are separated in energy from non-collective excitations, making it possible to describe them in the framework of effective field theory. Rotations and pairing rotations are the remnants of Nambu–Goldstone modes from the emergent breaking of rotational symmetry and phase symmetries in finite deformed and finite superfluid nuclei, respectively. The symmetry breaking severely constrains the structure of low-energy Lagrangians and thereby clarifies what is essential and simplifies the description. The approach via effective field theories exposes the essence of nuclear collective excitations and is defined with a breakdown scale in mind. This permits one to make systematic improvements and to estimate and quantify uncertainties. Effective field theories of collective excitations have been used to compute spectra, transition rates, and other matrix elements of interest. In particular, predictions of the nuclear matrix element for neutrinoless double beta decay then come with quantified uncertainties. This review summarizes these results and also compares the approach via effective field theories to well-known models and ab initio computations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Prediction of double-Weyl points in the iron-based superconductor CaKFe 4 As 4

Employing a combination of symmetry analysis, low-energy modeling, and ab initio simulations, we predict the presence of magnetic-field-induced Weyl points close to the Fermi level in CaKFe 4 As 4 . Depending on the relative strengths of the magnetic field and of the spin-orbit coupling, the Weyl fermions can carry a topological charge of ±1 or ±2, making CaKFe 4 As 4 a rare realization of a double-Weyl semimetal. Here, we further predict experimental manifestations of these Weyl points, both in bulk properties, such as the anomalous Hall effect, and in surface properties, such as the emergence of prominent Fermi arcs. Because CaKFe 4 As 4 displays unconventional fully gapped superconductivity below 30 K, our findings open a route to investigate the interplay between superconductivity and Weyl fermions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unraveling the magnetic ground state in the alkali-metal lanthanide oxide Na 2 ⁢Pr⁡O 3

Here, a comprehensive set of muon spin spectroscopy and neutron scattering measurements supported by ab initio and model Hamiltonian simulations have been used to investigate the magnetic ground state of Na 2 ⁢PrO 3 . μSR reveals a Néel antiferromagnetic order below T N ~ 4.9K, with a small static magnetic moment m static ≤ 0.22 μ B /Pr collinearly aligned along the c axis. Inelastic neutron measurements reveal the full spectrum of crystal field excitations and confirm that the Pr 4+ ground-state wave function deviates significantly from the Γ 7 limit that is relevant to the Kitaev model. Single- and two-magnon excitations are observed in the ordered state below T N =4.6K and are well described by nonlinear spin wave theory from the Néel state using a magnetic Hamiltonian with Heisenberg exchange J=1 meV and symmetric anisotropic exchange Γ/J=0.1, corresponding to an XY model. Intense two magnon excitations are accounted for by g-factor anisotropy g z /g ± = 1.29. A fluctuating moment δm 2 = 0.57 ⁢(22) ⁢μ$^2_B$/Pr extracted from the energy and momentum integrated inelastic neutron signal is reduced from expectations for a local J = 1/2 moment with average g factor g avg ≈ 1.1. Together, the results demonstrate that the small moment in Na 2 ⁢PrO 3 arises from crystal field and covalency effects and the material does not exhibit significant quantum fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coulomb excitation of the | T z | = 1 2 , A = 23 mirror pair

Background: Electric-quadrupole ( E 2 ) strengths relate to the underlying quadrupole deformation of a nucleus and present a challenge for many nuclear theories. Additionally, mirror nuclei in the vicinity of the line of N = Z represent a convenient laboratory for testing deficiencies in such models, making use of the isospin symmetry of the systems. Purpose: Uncertainties associated with literature E 2 strengths in Mg 23 are some of the largest in T z = | 1 2 | nuclei in the s d shell. The purpose of the present paper is to improve the precision with which these values are known, to enable better comparison with theoretical models. Methods: Coulomb-excitation measurements of Mg 23 and Na 23 were performed at the TRIUMF-ISAC facility using the TIGRESS spectrometer. They were used to determine the E 2 matrix elements of mixed E 2 / M 1 transitions. Results: Reduced E 2 transition strengths, B ( E 2 ) , were extracted for Mg 23 and Na 23 . Their precision was improved by factors of approximately 6 for both isotopes, while agreeing within uncertainties with previous measurements. Conclusions: A comparison was made with both shell-model and ab initio valence-space in-medium similarity renormalization group calculations. Valence-space in-medium similarity renormalization group calculations were found to underpredict the absolute E 2 strength, in agreement with previous studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

High-throughput determination of Hubbard $U$ and Hund $J$ values for transition metal oxides via the linear response formalism

DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 1000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an ATOMATE workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. Here, to demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-pU value results in a significantly improved agreement between the computed lattice parameters and experiment

36 MATERIALS SCIENCE↗