Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density Functionals”

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 OSTI2026

Mott vs Kondo: Influence of various density functional based methods on the Ce isostructural phase transition mechanism

The cerium iso-structural phase transition (γ–α) is dominated by f-electron localization changes that results in a magnetic ordering change and a volume collapse. Generally, these physics are difficult to capture with ab initio and first principles methods. However, previous works have shown various methods to be successful in predicting at least some of the physics of the γ to α phase transition. Therefore, here, we perform a broad survey of density functional based methods across three levels of theory and types of functions [GGA (generalized gradient approximation), MetaGGA, and hybrid functionals] and compare the results, focusing on hydrostatic compression across the phase boundary at zero kelvin. For the methods that best reproduce experimental results, we directly probe the predicted mechanisms and frame the results in the Mott/Kondo debate, assessing how the underlying methods and assumptions of different functionals can assess the physical drivers in the phase transition, providing insight into the governing processes of this unique phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗
DOE OSTI2025

Doping Effects on Multivalence States, Electronic Structure, and Optical Band Gap in LaCrO 3 under Varied Atmospheres: An Integrated Experimental and Density Functional Theory Study

Doping effects on the valence state, electronic structure, and optical band and the effects on electrical conductivity were studied on the doped lanthanum chromite (LaCrO 3 ) system. The specific compositions studied were La 1–x Ca x CrO 3 (LCCx), La 1–x Sr x CrO 3 (LSCx), and La 0.8 Sr 0.2 Cr 1–x Mn x O 3 (LSCMx) (0.1 ≤ x ≤ 0.4). The powders were synthesized using a modified Pechini sol–gel method, and the ceramic samples were densified using a reactive sintering method resulting in densities >96% theoretical. X-ray photoelectron spectroscopy (XPS) was completed to characterize the defect states and cationic valence compensation as a result of divalent (Ca 2+ or Sr 2+ ) and trivalent (Mn 3+ ) substitutions. XPS was completed for samples tested in oxidizing and reducing atmospheres (up to 1500 °C), which provided insights into the oxidation state transitions induced by the Ca 2+ and Sr 2+ dopants. The work notably demonstrated, for the first time, the oxidation/reduction transitions of Cr 4+ to Cr 3+ in Sr 2+ /Mn 3+ co-doped samples under reducing atmospheres. Reflectance UV–vis spectrophotometry optical band gap measurements were also completed for the same materials; a decrease in the optical band gap (2.81–3.12 eV) was shown with increased substitution, suggesting electronic structure modifications in the LaCrO 3 perovskite. Density functional theory calculations validated experimental trends, predicting a diminishing band gap with a rising dopant concentration. The transition in Cr oxidation states was attributed to the presence of divalent/trivalent cations. These findings contribute some insights into methods to tune the LaCrO 3 electrical properties for various low- and high-temperature applications.

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

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electronic structure of the kagome compound CaTi 3⁢ Bi 4 using high-field torque magnetometry and density functional theory

Here, we report systematic torque magnetometry measurements to investigate the electronic properties of the newly discovered kagome compound CaTi 3 ⁢Bi 4 . Electrical transport, magnetic susceptibility, and thermal measurements reveal no evidence of a magnetic ground state in this material. Torque data obtained in magnetic fields up to 41.5 T exhibit clear de Haas–van Alphen (dHvA) oscillations, with nine distinct frequencies ranging from 13 to 6164 T. Angular-dependent dHvA measurements show that, with the exception of the lowest frequency (13 T), all observed frequencies nearly follow a 1/cos ⁡𝜃 dependence, where 𝜃 is the angle between the crystallographic 𝑐 axis and the magnetic field direction. This behavior is characteristic of quasi-two-dimensional Fermi-surface sheets with nearly circular cross sections. To further elucidate the electronic structure, we performed density functional theory (DFT) calculations of the band structure and Fermi surface. The calculated bands reveal the presence of multiple Dirac points (DP), flat bands (FB), and van Hove singularities (VHS) near the Fermi level. The resulting Fermi surface consists of several quasi-two-dimensional cylindrical sheets, consistent with the experimentally observed 1/cos⁡ 𝜃 dependence. Notably, the theoretical dHvA frequencies, derived from extremal Fermi-surface cross-sectional areas, agree well with the experimental values and reproduce their angular dependence. Remarkably, all experimentally observed frequencies are captured by the DFT predictions. Pressure-dependent calculations up to 10 GPa show that the electronic features—DP, FB, and VHS—evolve systematically with pressure. In particular, the VHS shifts closer to the Fermi level, demonstrating that pressure acts as an effective tuning parameter in this material. These combined experimental and theoretical results provide a comprehensive understanding of the electronic structure of CaTi 3 ⁢Bi 4 and demonstrate how pressure can be used to tune its key electronic features, offering valuable guidance for exploring related kagome materials.

Shtefiienko, Kyryl [West Texas A & M University, C↗
DOE OSTI2024

Antimony stable isotope fractionation during adsorption onto birnessite: A molecular perspective from X-ray absorption spectroscopy and density functional theory

Sorption of antimony (Sb) onto birnessite significantly influences the fate of Sb in oceanic and terrestrial environments and fractionates Sb isotopes. Nevertheless, little is known about Sb isotopic fractionation during its adsorption on birnessite. Here, in this study, we show the value of Δ 123 Sb adsorbed-aqueous increases from −0.398 to −0.332 ‰ in 1 h and then decreases and stabilizes at −0.384 ‰ in 72 h. The enrichment of the light Sb isotope is predominantly due to the distortion of the octahedral symmetry. X-ray absorption spectroscopy results indicate Sb first forms a double-corner-sharing complex on birnessite and then transforms to a double-edge-sharing complex during adsorption. The optimized bond distances for double-corner-sharing (3.37 Å) and double-edge-sharing (2.90 Å) complexes calculated using density functional theory (DFT) fits well with the structure (3.41 and 3.00 Å) revealed by X-ray absorption spectroscopy, respectively. The fractionation derived from reduced partition function ratios calculated using DFT aligns well with the experimental results. Therefore, the variation in Sb isotopic fractionation during adsorption is attributed to the evolving structure of Sb complexes on birnessite. Our results demonstrate the isotopic fractionation of Sb during adsorption on birnessite and provide a molecular-scale understanding of Sb behavior, contributing to the correct reconstruction of the Sb isotope composition of ancient seawater using ferromanganese crusts and nodules, and efforts to trace Sb migration in epigenetic mining environments.

Adsorption↗
DOE OSTI2025

Evaluation of probability density function descriptions for three-component Rayleigh–Taylor mixing

Results from simulations of a three-component Rayleigh–Taylor (RT) mixing problem are presented. These simulations are conducted in heavy–light–heavy and heavy–intermediate–light configurations, and each of these configurations are further considered in high- and low-Reynolds-number regimes. This results in RT-unstable flow with one or both interfaces initially unstable, permitting the influence of problem configuration on the statistical description of three-component RT-driven mixing to be considered. Mass fraction covariances are observed to undergo a sign change through the mixing layer in all four configurations considered. This appears to be unique to the multi-component case and represents another way in which multi-component RT mixing differs from the two-component case. Qualitative and quantitative comparisons of joint and marginal probability density function (PDF) descriptions of species concentration are made. Three-, five-, and six-parameter model PDFs are compared against simulation data to assess how accurately they describe the mixing, and it is found that three-component mixing requires at least a five-parameter model PDF to accurately describe the mixing. Notably, the marginal distributions of three-component mixing do not appear to conform to a beta distribution, representing a departure from the classical two-component RT case. In conclusion, statistical neutrality also appears to influence the optimal choice of model PDF, which is found to be a function of problem configuration.

Large-eddy simulation↗
DOE OSTI2024

Simple self-consistent method for excited states in density functional theory to characterize defect-derived behavior in wide-band-gap-based microelectronic materials

This final report summarizes the results of the Laboratory Direct Research and Development (LDRD) Project Number 229740. Wide band gap semiconductors such as gallium nitride (GaN) have features highly desirable for multiple mission electronic applications. Realization of their potential requires atomic-scale understanding of electronic behavior. The principal experimental tools for electronically probing defects in GaN are chemically undifferentiating and lack a practical theoretical counterpart needed to identify and characterize specific defects. This project investigated whether a simple idea for modeling defect excited states and their associated photoluminescence (PL) energies is viable, as a path to accelerate the understanding of defect behavior and gain valuable insights into engineering new electronic materials and devices. The research implemented a non-self-consistent total-energy evaluation of a Koopmans-type estimation of an excited electronic state energy in density functional theory (DFT) calculations, and proceeded to design, implement, and assess a self-consistent method for computing excited states based upon an OCcupation-Constrained-DFT (occ-DFT). The occ-DFT was verified in test calculations of defect excited states and validated against well-characterized PL data for 3d transition metal defects in GaN. The method proved stable and robust in computing excited states and gave accurate predictions compared to experimental PL data. The combined ground state/excited-state capability proved capable of chemically differentiating defect species in GaN. In application to 3d dopants in GaN, we reinterpreted extensive experimental literature, proposed new defects as prospective candidates for use in quantum information applications, and outlined design strategies to create and exploit these potentially useful functional defects in GaN.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First principles density functional theory study of tritium species adsorption on Ni(111) surface and diffusion in nickel-sublayer for tritium storage

The nickel-plated zircaloy-4 is used as a tritium ( 3 H) getter in the tritium-producing burnable absorber rods (TPBARs) to capture 3 H produced in the 6 Li-riched annular γ-LiAlO 2 pellet under neutron irradiation. The experimental data and our previous theoretical results showed that the 3 H species produced from the γ-LiAlO 2 pellet were mainly 3 H 2 and 3 H 2 O. These 3 H species diffuse from the surface of the LiAlO 2 pellet across vacuum to the nickel-plated zircaloy-4 getter and then further diffuse into the getter to chemically form metal hydrides. While a number of studies show that oxygen binds strongly as compared to 3 H on the nickel (Ni) layer, the detailed mechanism of 3 H species absorption and diffusion across the Ni plate and Ni/Zr interface are still unclear. By employing density functional theory calculations, here we explored the 3 H 2 and 3 H 2 O species adsorption and dissociation on the Ni(111) surface and diffusion into the Ni sublayer. Our results indicated that the 3 H 2 and 3 H 2 O dissociate on the Ni(111) surface. The NiO x and Ni(O 3 H) x could be formed in the Ni layer due to the higher oxygen (O) diffusion energy barrier and formation of Ni vacancy defects. The oxygen was found to be retained in the Ni layer from diffusing across the Ni–Zr interface. This was revealed by comparing the diffusion barriers for 3 H with O. 3 H was found to have nearly three times smaller diffusion barrier than for O, making 3 H comparatively easier to diffuse through the Ni layer. In conclusion, the obtained results provide guidelines for experimental measurements on 3 H retention behavior in TPBARs and may open further avenues to explore the impurity effects on 3 H diffusion and storage at the Ni/zircaloy interfaces.

Tafen, De Nyago [National Energy Technology Lab. (↗
DOE OSTI2025

Evidence for Ga clusters in β-Ga 2 O 3 from Raman spectroscopy and density functional theory

Monoclinic gallium oxide (β-Ga 2 O 3 ) single crystals have a Raman mode at ∼250 cm −1 that is strongly correlated with free-electron density. Prior work attributed this peak to an electronic excitation of a shallow donor impurity band. However, heavily n-type thin films grown by metalorganic chemical vapor deposition or molecular beam epitaxy do not have the peak. In the present work, an alternate model is proposed: the 250 cm −1 Raman peak arises from Ga clusters, defined as two or more Ga atoms that form Ga–Ga bonds. Raman mapping reveals variations in the frequency that are consistent with a distribution of cluster sizes. The intensity of the peak decreases as the temperature is raised, attributed to melting of the Ga clusters. First-principles calculations indicate that the 250 cm −1 mode is due to Ga–Ga bond-stretching vibrations. As the Fermi energy is raised, the formation of Ga–Ga dimers becomes energetically favorable, explaining the correlation between n-type conductivity and the appearance of the Raman peak.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interactions involved in the adsorption of ethylene glycol and 2-hydroxyethoxide on the Au(111) surface: a Density Functional Theory study

Abstract Context:The monolayers of ethylene glycol and 2-hydroxyethoxide on gold surfaces have been used in hybrid materials as biosensors. In this article, the adsorption of ethylene glycol and 2-hydroxyethoxide on the Au(111) surface was analyzed. For the first system, ethylene glycol on Au(111), there are Au$$\cdot \cdot \cdot $$ · · · O and Au$$\cdot \cdot \cdot $$ · · · H interactions. To the best of our knowledge, the Au$$\cdot \cdot \cdot $$ · · · H interaction has been overlooked until now. However, in this work, there is strong evidence that this interaction is important to stabilize the system. For the second system, the atomic interactions mentioned previously are also predicted, although there is an additional interaction between 2-hydroxyethoxide molecules. Such an interaction induces the link -O-H-O-, with high values of the electron density at the critical points of the corresponding bond path of the O-H interaction. These links suggest the forming of ethylene glycol chains.Methods:The calculations were performed using two exchange-correlation functionals: BEEF-vdW and C09$$_{x}$$ x -vdW; both functionals incorporate dispersion effects within the Kohn-Sham approach in Density Functional Theory as implemented in GPAW code and ASE computational packages. The contacts between the molecules considered in this article and the Au(111) surface were analyzed through the Quantum Theory of Atoms in Molecules implemented in GPUAM code.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Range-separated thermal hybrid exchange-correlation density functional for accurate band-gap calculations of warm dense matter

Here, in this study, we present a range-separated thermal hybrid exchange-correlation (XC) functional, termed RS-KDT0, which provides highly accurate band-gap predictions at elevated temperature while reducing the computational cost relative to thermal global hybrid functionals. At zero temperature, RS-KDT0 recovers the widely used ground-state range-separated hybrid functional HSE06. To assess its performance across temperature regimes, we conducted static calculations of electronic pressure, band gaps, and band structures. RS-KDT0 exhibits the same level of accuracy as HSE06 in the low-temperature regime and demonstrates improved agreement with finite-temperature GW calculations at elevated temperatures when compared to other thermal XC functionals such as the generalized gradient approximation (GGA)-level KDT16 and the global hybrid KDT0.

density functional theory↗
DOE OSTI2024

Putting error bars on density functional theory dataset

This dataset contains submission files and raw output files from high-throughput DFT simulations to analyze the systemic errors in lattice constant, bulk moduli and formation energy predictions for a range of binary and ternary oxides using four exchange correlation functionals (LDA, PBE, PBEsol and vdW-DF-C09). This data was then used as the basis for employing materials informatics methods to predict the expected errors in the lattice constants of the studied compounds. Predicted errors were also used to better the DFT-predicted lattice parameters. Our results emphasize the link between the computed errors and the electron density and hybridization errors of a functional. In essence, these results provide “error bars” for choosing a functional for the creation of high-accuracy, high-throughput datasets as well as avenues for the development of XC functionals with enhanced performance, thereby enabling the accelerated discovery and design of new materials.

36 MATERIALS SCIENCE↗