Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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

Improving the Precision of First-Principles Calculation of Parton Physics from Lattice Quantum Chromodynamics

Large momentum effective theory (LaMET) provides a general framework for computing the multi-dimensional partonic structure of the proton from first principles using lattice quantum chromodynamics (QCD). In this effective field theory approach, LaMET predicts parton distributions through a power expansion and perturbative matching of a class of Euclidean observables—quasi-distributions—evaluated at large proton momenta. Recent advances in lattice renormalization, such as the hybrid scheme with leading renormalon resummation, together with improved matching kernel that incorporates higher-loop corrections and resummations, have enhanced both the perturbative and power accuracy of LaMET, enabling a reliable quantification of theoretical uncertainties. Moreover, the Coulomb-gauge correlator approach further simplifies lattice analyses and improves the precision of transverse-momentum-dependent structures, particularly in the non-perturbative region. State-of-the-art LaMET calculations have already yielded certain parton observables with important phenomenological impact. In addition, the recently proposed kinematically enhanced lattice interpolation operators promise access to unprecedented proton momenta with greatly improved signal-to-noise ratios, which will extend the range of LaMET prediction and further suppress the power corrections. The remaining challenges, such as controlling excited-state contamination in lattice matrix elements and extracting gluonic distributions, are expected to benefit from emerging lattice techniques for ground-state isolation and noise reduction. Thus, lattice QCD studies of parton physics have entered an exciting stage of precision control and systematic improvement, which will have a broader impact for nuclear and particle experiments.

Zhao, Yong [Argonne National Laboratory (ANL), Arg↗
DOE OSTI2023

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗
DOE OSTI2024

First-principles calculation of Hubbard U for Terbium metal under high pressure

Abstract Using density functional theory (DFT) and linear response approaches, we compute the on-site Hubbard interaction U of elemental Terbium (Tb) metal in the pressure range ∼ 0–65 GPa. The resulting first-principles U values with experimental crystal structures enable us to examine the magnetic properties of Tb using a DFT+U method. The lowest-energy magnetic states in our calculations for different high-pressure Tb phases—including hcp, α -Sm, and dhcp—are found to be compatible with the corresponding magnetic ordering vectors reported in experiments. The result shows that the inclusion of Hubbard U substantially improves the accuracy and efficiency in modeling correlated rare-earth materials. Our study also provides the necessary U information for other quantum many-body techniques to study Tb under extreme pressure conditions.

36 MATERIALS SCIENCE↗
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 OSTI2025

Atomic-Scale Insights into Carbon Dissolution in α-, γ-, and θ-Al 2 O 3 : Phase-dependent Transport Dynamics from First-Principles Calculations

α-Al 2 O 3 exhibits superior carburizing corrosion resistance compared to metastable γ-Al 2 O 3 and θ-Al 2 O 3 phases in high-temperature CO 2 environments, yet its atomic-scale origins remain unclear. Using first-principles density functional theory, we systematically investigate carbon dissolution and diffusion in α-Al 2 O 3 , γ-Al 2 O 3 , and θ-Al 2 O 3 , including the effects of oxygen (O) and aluminum (Al) vacancies. Our results show that α-Al 2 O 3 consistently exhibits higher carbon solution enthalpies than γ-Al 2 O 3 and θ-Al 2 O 3 in both pristine and defective structures, indicating lower intrinsic carbon solubility in α-Al 2 O 3 . Vacancies significantly enhance carbon incorporation: O vacancies reduce solution enthalpy, while Al vacancies further amplify this effect, with a strong preference for carbon at Al vacancy sites. Carbon diffusion barriers are also highest in α-Al 2 O 3 , reflecting slower carbon mobility. Al vacancies increase diffusion barriers across all phases, while O vacancies raise barriers in α- and γ-Al 2 O 3 but slightly lower them in θ-Al 2 O 3 . These results reveal a dual mechanism behind the carburizing resistance of α-Al 2 O 3 : reduced carbon solubility and elevated diffusion barriers. Furthermore, this work provides atomic-scale insights to guide the design of alumina-based materials with improved carburizing resistance through phase selection and defect engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exciton Lifetime and Optical Line Width Profile via Exciton–Phonon Interactions: Theory and First-Principles Calculations for Monolayer MoS 2

Exciton dynamics dictates the evolution of photoexcited carriers in photovoltaic and optoelectronic devices. However, interpreting their experimental signatures is a challenging theoretical problem due to the presence of both electron–phonon and many-electron interactions. Here, we develop and apply here a first-principles approach to exciton dynamics resulting from exciton–phonon coupling in monolayer MoS 2 and reveal the highly selective nature of exciton–phonon coupling due to the internal spin structure of excitons, which leads to a surprisingly long lifetime of the lowest-energy bright A exciton. Moreover, we show that optical absorption processes rigorously require a second-order perturbation theory approach, with photon and phonon treated on an equal footing, as proposed by Toyozawa and Hopfield. Such a treatment, thus far neglected in first-principles studies, gives rise to off-diagonal exciton–phonon self-energy, which is critical for the description of dephasing mechanisms and yields exciton line widths in excellent agreement with experiment.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Vibrational and electronic properties of Np 2 O 5 from experimental spectroscopy and first principles calculations

High-valence actinide oxides are critical to understanding the behavior of 5f-electrons, yet their structural and electronic properties remain poorly understood due to challenges in synthesis and handling. We report the first Raman spectroscopic study of single-crystalline Np 2 O 5 and the first scanning tunneling spectroscopy (STS) measurement on any neptunium-containing material. Hydrothermally synthesized crystals were structurally verified by X-ray diffraction. Raman spectra revealed sharply resolved vibrational features, including previously unreported low-frequency modes. STS measurements revealed a band gap of 1.5 eV. Density functional theory (DFT) enables vibrational mode assignments, reveals neptunium-dominated low-frequency phonons, oxygen-dominated high-frequency modes, and predicts an indirect band gap of 1.68 eV. This predicted value is in excellent agreement with the experimentally measured STS gap. This combined Raman, DFT, and STS approach provides a robust framework for correlating lattice dynamics and electronic structure in actinide materials, providing benchmark data for Np 2 O 5 , and opening new avenues for probing structure–property relationships in complex f-electron materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

High-Temperature Gas Sensor Materials with Properties Predicted via First-Principles Calculations with Machine Learning Modeling and Experimental Corroboration

Understanding the temperature dependence of functional properties of sensing materials is vital for their applications in combustion environments. The electron-phonon coupling that derives the electronic structure change with temperatures is a key property of interest as it affects other sensing responses. Herein, we first assess the temperature dependence of band gap renormalization in sensing materials by employing Allen-Heine-Cardona (AHC) theory with density functional theory (DFT) simulations corroborated with experimental observation. As the AHC calculations are impractical for high-throughput screening of materials, we employ data-driven Gaussian process regression to predict the parameters employed in the O’Donnell empirical model from a set of physical features. To mitigate the reliability issues arising from the small size of the dataset, we apply a Bayesian technique to improve the generalizability of the data-driven models as well as to quantify the uncertainty associated with theoretical predictions. These models capture well the overall trend of the O’Donnell parameters with respect to a reduced feature set obtained by transforming the available physical features. Quantifying the associated uncertainty helps us understand the reliability of the predictions and, therefore, the variation of bandgap as a function of temperature for other novel materials. The predicted candidates from machine learning models are further validated by experiments and DFT calculations.

bandgap renormalization↗
DOE OSTI2025

Temperature-Dependent Mechanical Properties of Ni-Based Concentrated Alloys: Insights from First-Principles Calculations

The present work focuses on predicting temperature-dependent mechanical properties of Ni-based concentrated alloys Ni 18 Cr 10 Co 10 Fe 6 M 4 (abbreviated by X 44 M 4 , with M = Al, V, Mn, Fe, Nb, Mo, and W) using density functional theory (DFT). These predictions are based on shear (plastic) and elastic deformations, utilizing the special quasirandom structure (SQS), the phonon-based quasiharmonic approach (QHA), and the quasistatic approach. The resulting properties include coefficient of thermal expansion via QHA, ideal shear strength (τ IS ), and stable and unstable stacking fault energies (γ SF and γ US ) through pure alias shear deformation, and elastic constants (c ij ), bulk modulus (B 0 ), and shear modules (G 0 ) via elastic deformation. Notably, predicting accurate γ SF is challenging due to uncertainties that can exceed the γ SF values. τ IS and γ US exhibit a strong linear relationship, enabling the accurate prediction of γ US based on the precisely determined τ IS . All mechanical properties of X 44 M 4 decrease with increasing temperature, except for some γ SF cases such as X 44 M 4 with M = V, Mn, Fe, Mo, and W. Among the X 44 M 4 alloys, X 44 Nb 4 exhibits the lowest τ IS , γ US , and G 0 values, and the highest B 0 /G 0 ratio, while X 44 Mn 4 has the lowest B 0 and B 0 /G 0 ratio. We found that volume is a crucial descriptor for understanding and modeling mechanical properties (except B0 and maybe also γ SF ) affected by alloying elements and temperature. Ni-based dilute alloys (e.g., Ni 11 M 1 and Ni 31 M 1 ) and concentrated alloys (e.g., X 44 M 4 ) show similar trends in mechanical properties influenced by alloying elements and temperature, simplifying the analysis and design of Ni-based alloys.

Elastic properties↗
DOE OSTI2022

CALPHAD modeling of uranium nitride (UN) fabrication routes enabled by first-principles calculations

The thermochemical details of fabricating uranium nitride (UN) by ammonolysis of uranium tetraflouride (UF 4 ) were determined using density functional theory (DFT) and CALculation of PHAse Diagrams (CALPHAD) computational methods. The thermochemical data of all binary, ternary, and quaternary U-H-N-F phases were computed using DFT, and the data for the phases that have not been measured experimentally, including UN 2 and NH 4 F(g), were combined with existing experimentally-determined data for CALPHAD modeling. The DFT data were benchmarked using experimental Gibbs energy of reaction and experimental thermochemical data for individual species. Phase diagrams relevant to the ammonolysis reaction are depicted, showing regions of stability for solid U-N, U-F and U-N-F phases. An unidentified phase produced in a previous experiment was identified as UN 0.95 F 1.2 (UNF) by comparing its X-ray diffraction spectrum to the experimental spectrum, and its formation during the fabrication of UN from UF 4 is supported by the simulated phase diagram. Here, it is calculated that UN 2 can be produced by the ammonolysis of UF 4 , but requires elevated temperatures, high NH 3 (g) partial pressure, and large amounts of flowing NH 3 (g) to avoid solid flu oride impurities in the uranium nitride. Likewise, U 2 N 3 can be produced instead at temperatures greater than 980 K. The use of silane (SiH 4 ) gas was investigated as a potential additive in the ammonolysis fabrication route to speed removal of fluorine. The addition of SiH 4 (g) offers little advantage to the removal of fluorine, and adds the complication of Si 3 N 4 formation. The use of DFT to fill in missing data to perform CALPHAD calculations demonstrated here allows for the determination of more comprehensive and trustworthy phase diagrams than the use of existing experimental data alone.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Influence of External Conditions on the Black-to-Yellow Phase Transition of CsPbI 3 Based on First-Principles Calculations: Pressure and Moisture

All-inorganic CsPbI 3 perovskite has emerged as a promising candidate for next-generation solar cells. However, owing to the poor structural stability of black phase CsPbI 3 perovskite at room temperature, it spontaneously transforms to the yellow, photoinactive, nonperovskite phase at room temperature, limiting further development of CsPbI 3 perovskite solar cells. Here, to better understand the mechanism driving such undesirable phase transformation, we examine the thermodynamics and kinetics associated with γ-to-δ phase transformation using density functional theory calculation. A solid-state nudged elastic band method was employed to find minimum energy pathways assuming a concerted solid-state phase transformation between these two phases. The gas molecules representing atmospheric (H 2 O, O 2 , and N 2 ) and inert (Ar) ambient environments, as well as applied external pressure, were considered to examine the influence of external conditions on the black-to-yellow phase transition. Our calculation reveals that, with increasing pressure, the thermodynamic driving force for converting the γ-phase to the δ-phase increases, whereas a larger activation energy must be overcome for the phase transition to occur. We also investigate the moisture-induced phase transformation with an H 2 O molecule and its dissociated species (H + /OH – ) and demonstrate that the reaction energy barriers can be significantly lowered in the presence of H 2 O or OH – . On the other hand, other nonpolar species (O 2 , N 2 , and Ar) have a negligible effect on the phase transformation kinetics, while they might reduce the thermodynamic driving force for the phase transformation and suppress the undesirable phase transformation. Our theoretical prediction could support recent observations that the γ-to-δ phase transformation occurs rapidly and catalytically in the presence of moisture, whereas in dry argon and dry oxygen atmospheres, the γ-CsPbI 3 remains stable.

14 SOLAR ENERGY↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗