Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
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 OSTI2026

Point defect energetics in gallium arsenide, a comprehensive density functional theory study

In materials, point defects often control or modify functional properties. To predict the performance of materials intended for application in optoelectronic devices, it is imperative to understand the properties of those point defects. For the first time, all six intrinsic defects of GaAs, a key optoelectronics material, and their charge transition levels are calculated using density functional theory with the HSE06 functional. For comparison, both PBE and r 2 SCAN calculations are also carried out. The HSE06 results are found to be in better agreement with experimental data than previous calculations. In conclusion, the importance of using the exact electron exchange present in hybrid functionals and larger supercells to accurately determine defect levels and ground state defect configurations is demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Density functional theory study on the coupling and reactions of diferuloylputrescine as a lignin monomer

Diferuloylputrescine has been found in a variety of plant species, and recent work has provided evidence of its covalent bonding into lignin. Results from nuclear magnetic resonance spectroscopy revealed the presence of bonding patterns consistent with homo-coupling of diferuloylputrescine and the possibility of cross-coupling with lignin. In the present work, density functional theory calculations have been applied to assess the energetics associated with radical coupling, rearomatization, and dehydrogenation for possible homo-coupled dimers of diferuloylputrescine and cross-coupled dimers of diferuloylputrescine and coniferyl alcohol. Furthermore, the values obtained for these reaction energetics are consistent with those reported for monolignols and other novel lignin monomers. As such, this study shows that there would be no thermodynamic impediment to the incorporation of diferuloylputrescine into the lignin polymer and its addition to the growing list of non-canonical lignin monomers.

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

Binding of uranyl cations to a Zr-based metal-organic framework by density functional theory

Zr-based metal–organic frameworks (Zr-MOFs) have been widely used as ion adsorbents for the removal or extraction of toxic and/or radionuclide species from aqueous solutions. However, the mechanisms by which uranyl cations (UO 2 2+ ) interact with Zr-MOFs have not been established. In this work, the nature of the bonding of uranyl cations with a Zr-MOF was determined using density functional theory for nineteen structurally distinct candidate complexes. Further, the results showed that in all cases the binding energy was of the order of 1.5 eV, but depended on the specific bonding site. The most stable structure involved coordination of the uranyl cation and two structurally distinct oxygens in the Zr-MOF metal node. It was also found that higher degree of deprotonation in Zr-MOF correlated with higher binding energy between the Zr-MOF and uranyl cations. These insights can aid in the design of Zr-MOFs with optimized features for efficient capture of uranyl cations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

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

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Small bimetallic clusters Agn-1M (M = Au, Co, Cu, Ni, Pd, Pt; n = 3, 9, 15): Density functional theory and genetic algorithm

We investigated the effect of size and composition on the properties of bimetallic nanoclusters. The geometric structures, stabilities, and electronic properties of size-selected Ag n-1 M (M = Au, Co, Cu, Ni, Pd, Pt; n = 3, 9, 15) bimetallic nanoclusters are systematically analyzed using spin-polarized density functional theory (DFT) within the generalized gradient approximation (GGA). We determine the most stable geometries for these clusters using a genetic algorithm (GA) in combination with DFT. Our results show that doping pure silver clusters with an M atom (transition metal), referred to as a “guest atom”, increases the stability as compared to pure Ag n (n = 3, 9, 15) clusters. The results for various properties including formation energy per atom, electronic structure, magnetic moments, and vibrational density of states (VDOS) are evaluated as a function of both size and composition of the system. The adsorption of selected bimetallic clusters on hydroxylated alumina substrate shows weak binding and minor changes in geometric properties except for Ag 8 Pt.

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

Density Functional Theory Study of Iron–Oxygen Divacancies in Magnetite (Fe 3 O 4 ) and Hematite (Fe 2 O 3 )

Density functional theory (DFT) calculations are employed to investigate the formation energies, charge redistribution, and binding energies of iron–oxygen divacancies in magnetite (Fe 3 O 4 ) and hematite (Fe 2 O 3 ). For magnetite, we focus on the low-temperature phase to explore variations with local environments. Building on previous DFT calculations of the variations in formation energies for oxygen vacancies with local charge and spin order in magnetite, we extend this analysis to include octahedral iron vacancies before analyzing the iron–oxygen divacancies. We also assessed the relative stability of iron–oxygen divacancies by comparing their formation energies with those of individual vacancies. Our findings reveal a significant energetic driving force for the formation of divacancy clusters, particularly in magnetite, where divacancies in the +1 charge state exhibit formation energies comparable to those of neutral iron vacancies under oxidizing conditions. In hematite, the results indicate a strong tendency for oxygen vacancies to bind to iron vacancies. These results highlight the significance of iron–oxygen vacancy complexes in the transport properties of iron oxides, with particular relevance to diffusion mechanisms under irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polymorphism and phase transitions in Na 2 U 2 O 7 from density functional perturbation theory

Polymorphism and phase transitions in sodium diuranate, Na 2 U 2 O 7 , are investigated with density functional perturbation theory (DFPT). Thermal properties of crystalline α-, β- and γ-Na 2 U 2 O 7 polymorphs are predicted from DFPT phonon calculations, i.e., the first time for the high-temperature γ-Na 2 U 2 O 7 phase (R$\bar{3}$ with combining macron]m symmetry). The standard molar isochoric heat capacities predicted within the quasi-harmonic approximation are C p (298.15K)= 219.4 and 220.9 J K -1 mol -1 mfor P21/a α-Na 2 U 2 O 7 and C2/m β-Na 2 U 2 O 7 , respectively. Gibbs free energy calculations reveal that α-Na 2 U 2 O 7 (P2 1 /a) and β-Na 2 U 2 O 7 (C2/m) are almost energetically degenerate at low temperature, with β-Na 2 U 2 O 7 becoming slightly more stable than α-Na 2 U 2 O 7 as temperature increases. Here these findings are consistent with XRD data showing a mixture of α and β phases after cooling of γ-Na 2 U 2 O 7 to room temperature and the observation of a sluggish α → β phase transition above ca. 600 K. A recently observed α-Na 2 U 2 O 7 structure with P21 symmetry is also shown to be metastable at low temperature. Based on Gibbs free energy, no direct β → γ solid-solid phase transition is predicted at high temperature, although some experiments reported the existence of such phase transition around 1348 K. This, along with recent experiments, suggests the occurrence of a multi-step process consisting of initial β-phase decomposition, followed by recrystallization into γ-phase as temperature increases.

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

Electron dynamics in extended systems within real-time time-dependent density-functional theory

Due to a beneficial balance of computational cost and accuracy, real-time time-dependent density-functional theory has emerged as a promising first-principles framework to describe electron real-time dynamics. Here we discuss recent implementations around this approach, in particular in the context of complex, extended systems. Results include an analysis of the computational cost associated with numerical propagation and when using absorbing boundary conditions. Here we extensively explore the shortcomings for describing electron–electron scattering in real time and compare to many-body perturbation theory. Modern improvements of the description of exchange and correlation are reviewed. In this work, we specifically focus on the Qb@ll code, which we have mainly used for these types of simulations over the last years, and we conclude by pointing to further progress needed going forward.

2D materials↗
DOE OSTI2020

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

3DRISM↗
DOE OSTI2023

Density Functional Theory Guided Investigation of Ligand‐Induced Neptunyl‐Neptunyl Interactions

Abstract Actinyl‐actinyl interactions are particularly prevalent for the pentavalent neptunyl cation (Np(V)O 2 ) + where these interactions appear either as a T ‐ or D ‐shape (diamond‐shape). T ‐shaped interactions have been previously identified in high concentration Np(V) solutions containing simple anions (NO 3 − , ClO 4 − , Cl − ) whereas D ‐shaped have only been isolated in the solid‐state in the presence of carboxylate ligands. In this study, Density Functional Theory (DFT) calculations were paired with Raman spectroscopy to evaluate the formation of D ‐shaped interactions in the presence of aliphatic (R=H (formate), CH 3 (acetate), CH 2 CH 3 (propionate)) and aromatic (R=C 6 H 5 (benzoate), C 6 H 4 OH (4‐hydroxybenzoate), C 5 H 4 N (isonicotinate)) carboxylate ligands. DFT studies indicate that the Δ G to form hydrated T ‐ and D ‐shaped forms are not spontaneous but become so with the addition of the carboxylate ligands. Raman spectra of the Np(V) carboxylate solutions contained vibrational modes associated with the D ‐shaped interactions, but spectral changes observed over time indicate a dynamic system. Crystallization experiments from the Np(V) carboxylate systems confirmed the presence of D‐ shaped dimers for the aromatic carboxylates, suggesting that the choice of the anion in solution favors actinyl‐actinyl interactions even at low concentrations (≤20 mM) of Np(V).

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

Density functional theory-based surrogate kinetic models for heterogeneous reactions of hydrocarbon intermediates on silicon carbide

The increasing demand for high-performance materials in advanced technologies highlights the importance of achieving a fundamental understanding and potential control of silicon carbide (SiC) deposition processes. However, existing models often lack sufficient theoretical detail, relying heavily on empirical data and offering limited predictive capability. In particular, the complex surface chemistry governing SiC growth remains poorly understood. This study addresses these challenges by employing density functional theory (DFT) to investigate key heterogeneous reactions involving hydrocarbon intermediates on SiC surfaces, including dehydrogenation, hydrogenation, and carbon deposition. Transition state searches were conducted to identify reaction pathways and energy barriers. While first-principles calculations offer high accuracy, they are computationally intensive. To extend the utility of these first-principles results, vibrational analyses were performed using phonon-based statistical thermochemistry to compute temperature-dependent reaction rates which were used to develop Arrhenius-type surrogate kinetic models. Furthermore, the resulting framework provides a more rigorous, physically grounded basis for integrating atomistic insights into continuum-scale modeling, ultimately enabling improved prediction and optimization of SiC film growth in high-performance material systems.

Density Functional Theory↗
DOE OSTI2020

A combined machine learning and density functional theory study of binary Ti-Nb and Ti-Zr alloys: Stability and Young’s modulus

The multicomponent Ti alloys, specifically the -phase, have experienced a strong growth over the last decades, due to their outstanding properties of ultra-high strength and low Young’s modulus. These properties play a significant role in many aerospace and biomedical applications. Selection and optimization of multicomponent alloys is challenging due to the vast chemical and compositional space. Here we investigate the use of machine learning techniques informed by density functional calculations to guide the selection of Nb- and Zr-based Ti binary alloys. From the cubic structures obtained from high throughput calculations and literature, we identify several structures with Young’s moduli below 40 GPa. The multivariant decision tree methods provide efficient surrogate models to identify structure variables have high influences on the energetic stability and Young’s modulus. Furthermore, we implement a workflow of incorporating DFT provided results and machine learning method to explore the chemical and composition space of other binary and multicomponent alloys, to eventually accelerate the material design via taking advantages of identified key variables.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Competitive effects of glucan’s main hydrolysates on biochar formation: A combined experiment and density functional theory analysis

We report the complexity of polysaccharide hydrothermal products increases the difficulty of exploring the formation of biochar, limiting the development of biochar. This work clarifies the completive effects of glucan’s main hydrolysates on biochar formation from three aspects: experimental, thermodynamic, and kinetic. The products distribution illustrates that 5-HMF, FA, and LA are mainly involved in the formation of biochar. Biochar mainly includes furan ring, ether group, and ester group by the analysis of magic-angle-spinning nuclear magnetic resonance, X-ray photoelectron spectroscopy, and elemental analysis. Combined experiments and density functional theory analysis, the etherification reaction of 5-HMF itself is most likely to occur and is key to form biochar, followed by the esterification of FA with 5-HMF, and then the etherification of 5-HMF and LA. The further verified experiments also manifest these results. This work will develop a foundation for exploring the complex formation mechanism of cellulose-based biochar.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Accurate frozen core approximation for all-electron density-functional theory

We implement and benchmark the frozen core approximation, a technique commonly adopted in electronic structure theory to reduce the computational cost by means of mathematically fixing the chemically inactive core electron states. The accuracy and efficiency of this approach are well controlled by a single parameter, the number of frozen orbitals. Explicit corrections for the frozen core orbitals and the unfrozen valence orbitals are introduced, safeguarding against seemingly minor numerical deviations from the assumed orthonormality conditions of the basis functions. A speedup of over twofold can be achieved for the diagonalization step in all-electron density-functional theory simulations containing heavy elements, without any accuracy degradation in terms of the electron density, total energy, and atomic forces. This is demonstrated in a benchmark study covering 103 materials across the Periodic Table and a large-scale simulation of CsPbBr3 with 2560 atoms. Our study provides a rigorous benchmark of the precision of the frozen core approximation (sub-meV per atom for frozen core orbitals below −200 eV) for a wide range of test cases and for chemical elements ranging from Li to Po. The algorithms discussed here are implemented in the open-source Electronic Structure Infrastructure software package.

Yu, Victor Wen-zhe↗
DOE OSTI2025

Design guidance for ferrites: Insights from density functional theory on magnetic properties

A grand challenge in materials research is identifying the relationship between composition and performance. Herein, we explore this relationship for magnetic properties, specifically magnetic saturation (M s ) and magnetocrystalline anisotropy energy (K) of ferrites. Ferrites are materials derived from magnetite (chemical formula = Fe 3 O 4 ) that comprise metallic elements such as Fe, Mn, Ni, Co, Cu and Zn. Experimentally, synthesizing and characterizing ferrites is time consuming. Further, selection of compositions to achieve particular magnetic properties currently relies on intuition. To address this, in this work, density functional theory (DFT) is used to predict M s and K for 571 ferrite structures. These structures are primarily double-substituted non-stoichiometric ferrites with formulae M1 x M2 y Fe 3–x–y O 4 , where M1 and M2 can be Mn, Ni, Co, Cu and/or Zn and 0 ≤ x ≤ 1 and y = 1–x. Calculated magnetic properties for the structures in this study vary from 0.04 × 10 5 to 9.6 × 10 5 A m −1 for M s and from 0.02 × 10 5 to 14.08 × 10 5 J m −3 for K. All structures are made publicly available in a FAIR database.

36 MATERIALS SCIENCE↗