Search NASA⌕ Search

SEARCH · Search NASA

Results for “Calculation”

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

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2024

Exciton Enhanced Nonlinear Optical Responses in Monolayer h-BN and MoS 2 : Insight from First-Principles Exciton-State Coupling Formalism and Calculations

Excitons are vital in the photophysics of materials, especially in low-dimensional systems. The conceptual and quantitative understanding of excitonic effects in nonlinear optical (NLO) processes is more challenging compared to linear ones. Here, we present an ab initio approach to second-order NLO responses, incorporating excitonic effects, that employs an exciton-state coupling formalism and allows for a detailed analysis of the role of individual excitonic states. Taking monolayer h-BN and MoS 2 as two prototype 2D materials, we calculate their second harmonic generation (SHG) susceptibility and shift current conductivity tensor. We find strong excitonic enhancement in the NLO responses requires that the resonant excitons are not only optically bright themselves but also able to couple strongly to other bright excitons. Our results explain the occurrence of two strong peaks in the SHG of monolayer h-BN and why the A and B excitons of MoS 2 unexpectedly exhibit minimal excitonic enhancement in both SHG and shift current generation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Benchmarking Density Functional Theory Methods for Efficient Calculations of a Strongly Correlated Li 1– x Ni 1– y O 2−δ System

Transition metal oxides (TMOs), such as LiNiO 2 , are promising candidates for energy storage and electronic devices due to their unique electronic properties, exceptional physical and chemical characteristics, and ability to adopt multiple oxidation states. However, accurately predicting their properties using mean-field density functional theory (DFT) is challenging due to the presence of strongly correlated d-electrons and the complex interplay between their structural, electronic, and magnetic responses. These challenges are further exacerbated by the need to model defects, surfaces, and interfaces, which require computationally efficient, large-scale simulations. To address these issues, we carry out a benchmark study on the Li 1–x NiO 2 system, evaluating the performance of several popular functionals. Our findings demonstrate that combining SCAN functional relaxation with single-step HSE calculations provides a practical and scalable computational strategy. This approach balances accuracy and efficiency, enabling high-throughput simulations of strongly correlated TMOs and improved predictive modeling capability of TMOs for practical applications.

25 ENERGY STORAGE↗
DOE OSTI2024

First-Principles Calculations of the Electrical Conductivity of Carbon Nanotubes Functionalized with Copper and Nitrogen: Implications for Electronics, Energy Storage, and Nanodevices

In this work, we investigate the electrical conductivity of carbon nanotubes (CNTs), with a particular focus on the effects of doping. Using first-principles-based approaches, we study the electronic structure, phonon dispersion, and electron–phonon scattering to understand the finite-temperature electrical transport properties in CNTs. Our study covers both prototypical metallic and semiconducting CNTs, with special emphasis on the influence of typical defects such as vacancies and the incorporation of copper or nitrogen, such as pyridinic N, pyrrolic N, graphitic N, and oxidized N. Our theoretical study shows significant improvements in the electrical conduction properties of copper-CNT composites, especially when semiconducting CNTs are functionalized with nitrogen. Doping is found to cause significant changes in the electronic density of states near the Fermi level, which affects the electrical conductivity. Calculations show that certain types of functional groups, such as N-pyrrolic, result in more than 30-fold increase in the conductivity of semiconducting CNTs compared to Cu-incorporated CNTs alone. For metallic CNTs, the conductivity is in agreement with existing experimental data, and our prediction of significant increases in conductivity with N-pyrrolic functional group is consistent with recent experimental results, demonstrating the effectiveness of doping in modifying conductivity. In conclusion, our study provides valuable insight into the electronic properties of doped CNTs and contributes to the development of ultrahigh conductivity CNT composites.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From Oxo to Oxyl to Biradical: Systematic Multireference Calculations of Methane Activation at MOF Nodes

Methane C–H activation at transition-metal sites often involves electronic structures that challenge conventional single-reference electronic structure descriptions. Although Kohn–Sham density functional theory (DFT) is widely used to study catalytic trends, its reliability for reactions involving strongly correlated species remains uncertain. Here we present a systematic multireference investigation of methane activation at metal–organic framework (MOF) node catalysts across the 3d transition-metal series. We introduce an automated workflow for active space selection to enable consistent application of multireference methods, including multiconfiguration pair-density functional theory and n-electron valence state perturbation theory, to these catalytic systems. These calculations show substantial static correlation in the C–H activation reaction step and predict activation barriers that differ from DFT by 30–70 kJ mol–1, with DFT often qualitatively disagreeing in barrier height trends across transition metals. Analysis of multireference wave functions shows that reactivity is governed by the electronic structure of the M–O moiety along a continuum from metal–oxo to oxyl radical and O biradical character. Increased oxygen-centered spin density and weakened M–O bonding are identified as descriptors of catalytic activity which correlate with lower activation barriers.

Wardzala, Jacob↗
DOE OSTI2025

DFT calculations and theory do not support enantiospecificity in NMR J-coupling constants

Distinguishing the enantiomers of small organic molecules is an industrially relevant problem with important implications for the health of the population. In a recent publication, Bouchard and co-workers have suggested that large differences in indirect spin-spin (J) coupling constants between enantiomers are possible. A close inspection of their work revealed significant flaws in their density functional theory (DFT) calculations and that the reported effects disappear with appropriate care. We thus conclude that enantiospecificity in spin-spin coupling constants has not be demonstrated either experimentally or theoretically.

Perras, Frédéric A. [Ames Laboratory (AMES), Ames,↗
DOE OSTI2025

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2025

Importance of enforcing Hund’s rules in density functional theory calculations of rare earth magnetocrystalline anisotropy

Density functional theory (DFT) and its extensions, such as DFT+U and DFT+dynamical mean-field theory, are invaluable for studying magnetic properties in solids. However, rare-earth (R) materials remain challenging due to self-interaction errors and the lack of proper orbital polarization. We show how the orbital dependence of self-interaction error contradicts Hund’s rules and plagues magnetocrystalline anisotropy (MA) calculations, and how analyzing DFT states that respect Hund’s rules can mitigate this issue. We benchmark MA in RCo 5 , R 2 Fe 14 B, and RFe 12 , extending prior work on RMn 6 Sn 6 , achieving excellent agreement with experiments. Additionally, we illustrate a semi-analytical perturbation approach that treats crystal fields as a perturbation in the large spin-orbit coupling limit. Using Gd-4f crystal-field splitting, this method provides a microscopic understanding of MA and enables rapid screening of high-MA materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Database and deep-learning scalability of anharmonic phonon properties by automated brute-force first-principles calculations

Understanding the anharmonic phonon properties of crystal compounds—such as phonon lifetimes and thermal conductivities—is essential for investigating and optimizing their thermal transport behaviors. These properties also impact optical, electronic, and magnetic characteristics through interactions between phonons and other quasiparticles and fields. In this study, we develop an automated first-principles workflow to calculate anharmonic phonon properties and build a comprehensive database encompassing more than 6500 inorganic compounds. Utilizing this dataset, we train a graph neural network model to predict thermal conductivity values and spectra from structural parameters, demonstrating a scaling law in which prediction accuracy improves with increasing training data size. High-throughput screening with the model enables the identification of materials exhibiting extreme thermal conductivities—both high and low. The resulting database offers valuable insights into the anharmonic behavior of phonons, thereby accelerating the design and development of advanced functional materials.

Ohnishi, Masato [University of Tokyo (Japan); Inst↗
DOE OSTI2026

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗
DOE OSTI2025

Contextual subspace variational quantum eigensolver calculation of the dissociation curve of molecular nitrogen on a superconducting quantum computer

Abstract We present an experimental demonstration of the Contextual Subspace Variational Quantum Eigensolver on superconducting hardware. Calculating the potential energy curve of molecular nitrogen proves challenging for many conventional quantum chemistry techniques, since static correlation dominates in the dissociation limit. Our quantum simulations retain good agreement with the Full Configuration Interaction energy, outperforming all benchmarked single-reference wavefunction techniques in capturing the bond-breaking appropriately. Moreover, our methodology is competitive with multiconfigurational approaches but at a saving of quantum resource, meaning larger active spaces can be treated for a fixed qubit allowance. To achieve this result, we deploy an error mitigation/suppression strategy comprised of Dynamical Decoupling, Measurement-Error Mitigation and Zero-Noise Extrapolation. Circuit parallelization also provides passive noise-averaging and improves the effective shot yield to reduce the measurement overhead. Furthermore, we introduce a modified adaptive ansatz construction algorithm that incorporates hardware awareness into our variational circuits, minimizing the transpilation cost for the target qubit topology.

Physics↗
DOE OSTI2026

Low-cost calculation and analysis of 2D IR spectra of model diiron trinitrosyl complexes in the NO stretch region with vibrational perturbation theory

Two-dimensional infrared spectroscopy offers unique capabilities for probing vibrational coupling in complex metal–ligand systems. In this paper, we combine two-dimensional infrared spectroscopy with vibrational perturbation theory to investigate vibrational coupling in a diiron trinitrosyl complex across three stable redox states. Although these systems are challenging for electronic structure methods, we demonstrate that key features of experimental 2D IR spectra can be accurately reproduced using reduced-dimensional anharmonic calculations with a small harmonic frequency scaling. Analysis reveals that N–O stretching modes maintain high locality across all redox states, with coupling patterns that directly reflect variations in Fe–N bond strength. Using curvilinear coordinate analysis, we demonstrate these differences result from systematic changes in cubic anharmonic force constants rather than mode delocalization. Our results establish N–O stretches as sensitive probes of metal–ligand bonding strength, expanding the toolkit for studying biologically relevant nitrosyl complexes.

Moran, Hayden A. [Texas A&M University, College St↗
DOE OSTI2026

Investigating trends in actinide covalency and magnetism with 35/37 Cl SSNMR spectroscopy and first-principles calculations

Solid-state NMR (SSNMR) spectroscopy is a powerful technique for studying actinide chemistry but has been significantly limited due to the complex paramagnetism, and radiological hazards presented by these materials. Lanthanide and actinide salts often feature magnetic ordering and can be paramagnetic, ferromagnetic, or antiferromagnetic depending on temperature and electronic structure. Paramagnetic interactions can manifest in SSNMR both as secular spectral shifts and/or couplings as well as contributions from non-secular relaxation. Both effects can be directly measured with NMR and used to extrapolate rich chemical information such as coordination environments, bonding characteristics, local molecular dynamics, and correlation times. Typically, these studies are carried out on high-γ and highly abundant NMR-active isotopes (e.g., 1 H, 6/7 Li, 19 F, 23 Na, etc.) or on enriched rare isotopes (e.g., 2 H and 17 O), which can be expensive. Herein, we present a facile methodology to measure the 35/37 Cl electric-field gradient (EFG) and paramagnetic shift anisotropy (SA) tensor components using static wideline SSNMR measurements of LaCl 3 , NdCl 3 , UCl 3 , and UCl 4 . The static powder spectra were measured with both 35 Cl and 37 Cl SSNMR to increase the fidelity of the extracted tensor parameters. Variable temperature NMR of a select case confirms the Curie–Weiss paramagnetism. Relaxation measurements of both nuclei further corroborate observations owing to the paramagnetic relaxation enhancement and reveal simultaneous quadrupolar relaxation mechanisms. Density functional theory (DFT) calculations using Hubbard U corrections to the uranium valence orbitals show excellent agreement with experimental EFG tensor parameters and help describe the bonding characteristics in these lanthanide and actinide systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Investigating PVC polymer–plasticizer interactions with atomistic MD simulations and potential of mean force calculations

For this work, atomistic molecular dynamics (MD) simulations coupled with potential of mean force (PMF) calculations were employed to investigate the interactions between PVC polymer chains containing 6–20 repeating units and various plasticizers, with the goal of identifying potential replacements for the toxic plasticizer di(2-ethylhexyl)phthalate (DEHP) used in blood bags. The selected plasticizers belong to various chemical families, such as orthophthalates, citrates, adipates, and the terephthalate DEHT. Both the polymer and the plasticizers lack ionizable groups, and their interactions are primarily governed by van der Waals and electrostatic forces. A model correlating PMF profiles with interaction forces was developed and validated across polyvinyl chloride (PVC) polymers of different lengths and all investigated plasticizers. This model provides insight into how structural variations in plasticizers influence their respective PMF values. The study ranks the investigated plasticizers based on binding affinity, identifying TOTM (tris(2-ethylhexyl) trimellitate, TEHTM), BTHC (butyryl trihexyl citrate, Citroflex B-6), ATHC (acetyl trihexyl citrate, Citroflex A-6, CA-6), and DEHT (bis(2-ethylhexyl)terephthalate, DOTP/DEHTP) as promising alternatives for further investigation. This comprehensive study encompasses PVC polymers of four different lengths and 14 plasticizers, with all data averaged over 30 independent simulations. The approach provides a deeper understanding of molecular interactions, enabling the tailoring of polymer–plasticizer systems for diverse applications, including extractables and leachables, with relevance spanning materials science to biomedical engineering, and also serves as a basis for developing coarse-grained simulation protocols. Overall, this study provides valuable insights for designing safer and more efficient plasticizer substitutes and is well supported by other studies.

Shet, Sai Athmeeya G. [Sri Sathya Sai Institute of↗
DOE OSTI2024

Reconciling calculations and measurements of inverse bremsstrahlung absorption

It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

PyLRO: A Python calculator for analyzing long-range structural order

We present PyLRO, an open-source Python calculator designed to detect, quantify, and display long-range order in periodic structures. The program’s design methodology, workflow, and approach to order quantification are described and demonstrated using a simple toy model. Additionally, we apply PyLRO to a series of metastable AlPO 4 structural intermediates from a prior high-pressure study, demonstrating how to compute and visualize structural order in all directions on a Miller sphere. We further highlight the program’s capabilities through a high-throughput analysis of structural patterns in the pressure-induced amorphization of AlPO 4 , revealing atomistic insights into specific energy regions of massive amorphous structures. These results suggest that PyLRO can be a valuable tool for investigating crystal–amorphous transition in materials research.

36 MATERIALS SCIENCE↗