Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory”

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

Support of the Third Summer School and Workshop Focused on Theory and Applications of Time-Dependent Density Functional Theory, Electronic Excited States, and Excited State Dynamics

To create a community of scholars with a deep understanding of TDDFT and excited state methods, both their capabilities and limitations, and how they can be used to model the excited state properties and charge transfer in materials and molecules, we are planning the third iteration of a summer school and workshop focusing on theory and applications of TDDFT, Excited States, and Excited State Dynamics to be offered in the summer of 2023. The two iterations of the US-based school that we previously organized have been very successful with overwhelmingly positive feedback. There is a high level of engagement, enthusiasm, and active participation and networking amongst the students of the school and the scientists of the workshop.

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

Density functional theory calculations of the mixing enthalpy of ternary uranium carbide compounds

The high melting point of uranium-zirconium carbides (U,Zr)C makes them an ideal fuel for nuclear thermal propulsion (NTP) reactors. Gaps remain in the current understanding of the U-Zr-C system due to the difficulty of conducting thermodynamic experiments at NTP operation conditions. Density functional theory calculations using the Hubbard U model (DFT+U) were performed using orbital matrix occupation (OMC) to obtain the mixing enthalpy for UC and ZrC for (U,Zr)C ternary compounds. Similarly, DFT+U calculations were also carried out for the (U,Nb)C and (U,Ta)C systems. In conclusion, the DFT results are envisioned to be used in thermodynamic assessments of the uranium carbide systems based on the CALPHAD approach to supplement the lack of experimental data for the mixing thermodynamics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Density Functional Theory Study of Epitaxially Strained Monolayer Transition Metal Chalcogenides for Piezoelectricity Generation

Two-dimensional transition metal chalcogenides (2D TMCs) are known for their wide range of bandgaps, flexibility, and high strength. Recent synthesis and data mining efforts indicate that 56 2D TMCs have low exfoliation energies and are relatively stable in monolayer form. Under epitaxial strain, we predict using density functional theory (DFT) calculations that the majority of these 2D TMCs can accommodate ±10% strain without breaking their crystal symmetry. The elastic and piezoelectric tensors indicate that 22 of 56 candidates are piezoelectric, and we derive their in-plane piezoelectric coefficient d 11 . The epitaxial strain is further predicted to enhance the d 11 by over 100% at 10% tensile epitaxial strain for most of these piezoelectric 2D TMCs. ReSe 2 at pristine state and Au 2 Se 2 at +5% epitaxial strain are predicted to obtain the extreme d 11 coefficients at -120 and 326 pm/V, respectively. Lastly, these findings have implications for the use of high-performance 2D piezoelectric materials in devices.

2D materials↗
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 OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2020

Molecular Structure and Electron Affinity of Metal-Solvent Complexes: Insights from Density Functional Theory Simulations

A molecular level understanding of the structure and energetics of the monovalent and divalent metal ion complexes is of great importance for development of next-generation batteries. Here, Density Functional Theory (DFT) simulations at the ωb97xD/6-31 + G(d,p) level of theory are performed to investigate the interaction of metal ions (Li + , Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ ) with 26 organic solvent molecules. The reduction energetics (electron affinity and reduction potential) and structural responses of the solvent molecules and the molecular complexes are discussed. The DFT calculations are carried out to investigate the structure, energetics and electron affinities of chelated complexes of water (H 2 O), tetrahydrofuran (THF) and di-methoxy ethane (DME) solvent molecules. Additionally, ab initio dynamic simulations (AIMD) at 298 K using atom centered density matrix propagation (ADMP) formalism are performed to understand the spontaneous structure formation upon electron attachment of the metal ion-solvent complexes. The ADMP simulations indicate the decomposition of Mg + -(DME) 3 complex via cleavage of C–O bond of one of the three DME molecules indicating irreversible decomposition of DME in the presence of the Mg + radical. We believe that the data collected as part of this investigation serves as a library of fundamental knowledge towards a deeper understanding of the electrode-electrolyte interfacial reactions.

25 ENERGY STORAGE↗
DOE OSTI2020

Radionuclide incorporation in negative thermal expansion α-Zr(WO 4 ) 2 : A density functional theory study

The incorporation of uranium, plutonium and technetium in the negative thermal expansion (NTE) α-Zr (WO 4 ) 2 has been investigated within the framework of density functional theory (DFT). It is found that the vacancy formation energies of the charged vacancies are overall larger than that of its counterpart neutral Frenkel defects and Schottky defects. DFT calculations suggest that U and Pu substitutions for the Zr site are preferred in α-Zr (WO 4 ) 2 . In case of Tc substitution, both Tc(IV) for the Zr site and Tc(VII) for the W site are considered under oxygen-poor and oxygen-rich conditions, while Tc(VII) substitution can be improved significantly by including Y 2 O 3 (charge compensation).

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

Toward a better understanding of Ni coarsening in solid oxide cells: NiH on Ni (111) examined at the level of density-functional theory (and KMC)

Ni coarsening in the fuel electrode of solid oxide cells (SOCs) is known to be significantly faster under a humid atmosphere. In this talk, density-functional theory and kinetic Monte Carlo methods are used to explore the hypothesis that the surface diffusion of NiH on Ni may promote Ni coarsening in the SOC operated in electrolysis cell mode. Defining the surface diffusivity as the product of the fractional surface coverage and single-molecule surface diffusivity, the surface diffusivity of NiH on Ni (111) is found to be large enough under a significant overpotential to support the above hypothesis. However, as NiH could dissociate on Ni (111), more work is needed to show that NiH may promote Ni coarsening.

Mantz, Yves↗
DOE OSTI2025

Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory

Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.

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 OSTI2025

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

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

Irreversible oxygen poisoning: Modeling the initial water dissociation kinetics on δ-Pu(111) and δ-Pu(100) through density functional theory

In this work, the initial kinetics of water dissociation on two facets of δ-plutonium, δ-Pu(111) and δ-Pu(100), are explored through density functional theory in order to understand how water dissociation occurs on these facets. We explored the dissociation of water via the formation of hydroxyls, atomic hydrogen, and atomic oxygen species on each facet. We calculate low energetic barriers for water to split to adsorbed hydrogen and hydroxyl species at 0.19 eV for δ-Pu(111) and 0.07 eV for δ-Pu(100). The hydroxyl has a barrier of 0.64 and 0.37 eV to cleave the hydrogen–oxygen bond on δ-Pu(111) and δ-Pu(100), respectively. Due to the highly exergonic adsorption free energy of atomic oxygen of −2.10 eV, the metallic surfaces are found to be fully covered in oxygen, even with the inclusion of oxygen lateral interactions. When combined with the reaction thermodynamics, this free energy of adsorption forms a molecular oxygen desorption barrier greater than 9.51 eV (918 kJ/mol). These results, combined with simulated temperature programed desorption spectra, indicate that oxygen formed via water dissociation on δ-Pu (111) and (100) facets induces an irreversible poison that prevents further reaction of water directly on metallic plutonium surfaces, most likely due to the strong hybridization of the Pu and O valence states. Therefore, these results imply that another mechanism is responsible for the continuously experimentally measured water dissociation, which is most likely due to the Pu oxide.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crystal growth, density functional theory, and scintillation properties of Tl 3 LnCl 6 :Ce 3+ and TlLn 2 Cl 7 :Ce 3+ (Ln = Y, Gd)

In this paper we report on the crystal growth, density functional theory (DFT) calculations and scintillation properties of Tl 3 LnCl 6 :Ce and TlLn 2 Cl 7 :Ce(Ln = Y, Gd). Crystals were grown by the Vertical Bridgman technique up to 16 mm in diameter and 25 mm long. Crystals of Tl 3 YCl 6 :Ce and Tl 3 GdCl 6 :Ce belong to the family of the Nesohalides which have the monoclinic crystal structure. Crystals of TlY 2 Cl 7 :Ce and TlGd 2 Cl 7 :Ce belong to the family of the Phyllohalides which have either the monoclinic or the orthorhombic crystal structure. We report the light yields of these Nesohalides and Phyllohalides are typically on the order of 35,000–40,000 ph/MeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2025

Activation, Dehydrogenation, and Carbon–Carbon Coupling of Methane by Iridium Cations Studied by Infrared Multiple Photon Dissociation Spectroscopy and Density Functional Theory

Products resulting from the sequential activation of one, two, three, and four methane molecules by atomic iridium cations were characterized by gas-phase infrared multiple photon dissociation spectroscopy and density functional theory (DFT) calculations. Iridium cations were generated using a laser ablation source and reacted with methane in a linear radiofrequency ion trap before mass analysis and spectroscopic interrogation in a Fourier transform ion cyclotron resonance mass spectrometer coupled to the free-electron laser for intracavity experiments (FELICE) beamline. Product ions were irradiated using infrared light over the 250–1500 cm –1 range. Comparisons between the experimental and DFT-calculated spectra enabled structural determination of the products formed. The observed products include HIrCH + , Ir(CH 2 ) 2 + , H s Ir(C 3 H 5 ) + , and Ir(CH 3 ) s (C 3 H 5 ) + , where the subscript s denotes a syn orientation of the two ligands. Furthermore, formation of the latter two products provides evidence for efficient C–H bond activation and subsequent C–C coupling on the atomic iridium cation.

Chemical reactions↗