Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory calculations”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Emergent multiferroism with magnetodielectric coupling in EuTiO 3 created by a negative pressure control of strong spin-phonon coupling

Negative pressure has emerged as a powerful tool to tailor the physical properties of functional materials. However, a negative pressure control of spin-phonon coupling for engineering magnetism and multiferroicity has not been explored to date. Here, using uniform three-dimensional strain-induced negative pressure in nanocomposite films of (EuTiO 3 ) 0.5 :(MgO) 0.5 , we demonstrate an emergent multiferroicity with magnetodielectric coupling in EuTiO 3 , matching exactly with density functional theory calculations. Density functional theory calculations are further used to explore the underlying physics of antiferromagnetic-paraelectric to ferromagnetic-ferroelectric phase transitions, the spin-phonon coupling, and its correlation with negative pressures. The observation of magnetodielectric coupling in the EuTiO 3 reveals that an enhanced spin-phonon coupling originates from a negative pressure induced by uniform three-dimensional strain. Our work provides a route to creating multiferroicity and magnetoelectric coupling in single-phase oxides using a negative pressure approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding the Effect of Single Atom Cationic Defect Sites in an Al 2 O 3 (012) Surface on Altering Selenate and Sulfate Adsorption: An Ab Initio Study

Adsorption is a promising under-the-sink selenate remediation technique for distributed water systems. Recently it was shown that adsorption induced water network rearrangement control adsorption energetics on the α-Al 2 O 3 (012) surface. Here, we aim to elucidate the relative importance of the water network effects and surface cation identity on controlling selenate and sulfate adsorption energy using density functional theory calculations. Density functional theory (DFT) calculations predicted the adsorption energies of selenate and sulfate on nine transition metal cations (Sc–Cu) and two alkali metal cations (Ga and In) in the α-Al 2 O 3 (012) surface under simulated acidic and neutral pH conditions. We find that the water network effects had a larger impact on the adsorption energy than the cationic identity. However, cation identity secondarily controlled adsorption. Most cations decreased the adsorption energy, weakening the overall performance, the larger Sc and In cations enabled inner-sphere adsorption in acidic conditions because they relaxed outward from the surface, providing more space for adsorption. Additionally, only Ti induced Se selectivity over S by reducing the adsorbing selenate to selenite but not reducing the sulfate. Altogether, this study indicates that tuning water network structure will likely have a larger impact than tuning cation–selenate interactions for increasing adsorbate effectiveness.

Adsorption↗
DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Correct symmetry in density functional theory calculations of δ-Pu [Slides]

Density functional theory (DFT) calculations in the delta-phase of Pu have long relied on the antiferromagnetic (AFM) spin structure to approximate the material’s unique electronic structure. The resulting tetragonal distortion away from the correct cubic symmetry was believed to be a necessary consequence of allowing an ordered spin structure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Energetic Stability of Freestanding and Metal-Supported Borophenes: Quantum Monte Carlo and Density-Functional Theory Calculations

Quantum Monte Carlo (QMC) and density-functional theory (DFT) calculations were carried out to study cohesion energetics of two-dimensional (2D) sheets of boron atoms called borophenes. Our QMC calculations confirmed the polymorphism among free-standing borophenes that was reported in previous DFT studies. Although Perdew-Burke-Ernzerhof (PBE) calculations significantly overestimate the cohesive energies of 2D boron sheets, DFT-PBE relative energetics with respect to each other among various free-standing borophenes are found to be in quantitative agreement with the corresponding QMC results. This suggests that one can make reliable predictions for relative stability of different boron sheets through DFT-PBE calculations. Our analysis of PBE formation energies of borophenes on metal surfaces shows that the polymorphic range is extended for borophenes on the Ag(111) and the Au(111) surfaces beyond that of free-standing borophenes, reflecting recent experimental synthesis of beta(12) and chi(3) boron sheets on the Ag(111) surface. We have also found that a hexagonal borophene can be stabilized through charge transfer from a metal surface and is energetically favored on the AI(111) surface over other borophene structures. Finally, it is found that the bilayer formation could be energetically favored over its monolayer form for the borophene-Au system, especially for borophenes with hexagonal hole densities eta lower than 1/9. This leads to our prediction that in addition to its monolayer form, a bilayer eta = 1/12 borophene can be synthesized on the Au(111) surface, opening a new possibility for borophene-based electronic devices.

ab initio↗
DOE OSTI2025

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2022

Symmetry-correct bonding in density functional theory calculations for delta phase Pu

The long-held conclusion that magnetic order in delta phase Pu makes the structure mechanically unstable in density functional theory calculations is incorrect if the magnetic order is allowed to be non-collinear. The non-collinear 3Q spin structure is intrinsically cubic and makes all structurally equivalent bonds equivalent in their bonding character. Applied to density functional theory calculations on delta phase Pu, the 3Q spin structure results in elastic constants and phonons with the correct symmetry. Finally, the calculated phonon dispersion agrees better with experiment than previous calculations, and the calculated thermal expansion shows the unique behavior seen experimentally: it is negative.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A robust solver for wavefunction-based density functional theory calculations

A new iterative solver is proposed to efficiently calculate the ground state electronic structure in density functional theory calculations. This algorithm is particularly useful for simulating physical systems considered difficult to converge by standard solvers, in particular metallic systems. Here, the effectiveness of the proposed algorithm is demonstrated on various applications.

42 ENGINEERING↗
DOE OSTI2023

Structural Stability of Graphene-Supported Pt Layers: Diffusion Monte Carlo and Density Functional Theory Calculations

In this study, we used a combination of diffusion Monte Carlo and density functional theory calculations to investigate the stability and interlayer binding of various layered structures of Pt atoms adsorbed on graphene. Our findings show that vertically buckled Pt monolayer and bilayer with (111)-packing order are more energetically favorable than the corresponding buckled (100) or flat (100)-packing structures. This can be attributed to the significant lattice mismatch (>10%) between pristine graphene and a free-standing (100)-packing Pt layer. Additionally, our calculations reveal that among the (100)-packing Pt layers, an incommensurate structure with a Pt/C atomic ratio less than 1/2 may be more stable than the commensurate structures registered at the bridge sites, which aligns with recent experimental findings of incommensurate (100)-packing Pt layers on graphene. The interlayer binding between the Pt layer and graphene is found to be primarily driven by van der Waals interaction, except for the AA-stacked buckled-(100) Pt bilayer where the bottom Pt atoms show chemisorption to the graphene surface. In conclusion, this research offers a comprehensive examination of the stability and interlayer binding of metallic Pt layers, providing valuable insights for potential applications as a next-generation catalyst.

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

Challenges for density functional theory: calculation of CO adsorption on electrocatalytically relevant metals

Density Functional Theory (DFT) is currently the most tractable choice of theoretical model used to understand the mechanistic pathways for electrocatalytic processes such as CO 2 or CO reduction. In this work, we assess the performance of two DFT functionals designed specifically to describe surface interactions, RTPSS and RPBE, as well as two popular meta-GGA functionals, SCAN and B97M-rV, that have not been a priori optimized for better interfacial properties. We assess all four functionals against available experimental data for prediction of bulk and bare surface properties on four electrocatalytically relevant metals, Au, Ag, Cu, and Pt, and for binding CO to surfaces of these metals. To partially mitigate issues such as thermal and anharmonic corrections associated with comparing computations with experiments, molecular benchmarks against high level quantum chemistry are reported for CO complexes with Au, Ag, Cu and Pt atoms, as well as the CO–water complex and the water dimer. Overall, we find that the surface modified RPBE functional performs reliably for many of the benchmarks examined here, and the meta-GGA functionals also show promising results. Specifically B97M-rV predicts the correct site preference for CO binding on Ag and Au (the only functional tested here to do so), while RTPSS performs well for surface relaxations and binding of CO on Pt and Cu.

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

Nonconventional NMR Spin-Coupling Constants in Oligosaccharide Conformational Modeling: Structural Dependencies Determined from Density Functional Theory Calculations

Nonconventional NMR spin-coupling constants were investigated to determine their potential as conformational constraints in MA’AT modeling of the O-glycosidic linkages of oligosaccharides. Four ( 1 J C1',H1' , 1 J C1',C2' , 2 J C1',H2' , and 2 J C2',H1' ) and eight ( 1 J C4,H4 , 1 J C3,C4 , 1 J C4,C5 , 2 J C3,H4 , 2 J C4,H3 , 2 J C5,H4 , 2 J C4,H5 , and 2 J C3,C5 ) spin-couplings in methyl β-d-galactopyranosyl-(1→4)-β-d-glucopyranoside (methyl β-lactoside) were calculated using density functional theory (DFT) to determine their dependencies on O-glycosidic linkage C–O torsion angles, $\phi$ and ψ, respectively. Long-range 4 J H1',H4 was also examined as a potential conformational constraint of either $\phi$ or ψ. Secondary effects of exocyclic (hydroxyl) C–O bond rotation within or proximal to these coupling pathways were investigated. Based on the findings of methyl β-lactoside, analogous J-couplings were studied in five additional two-bond O-glycosidic linkages [βGlcNAc-(1→4)-βMan, 2-deoxy-βGlc-(1→4)-βGlc, αMan-(1→3)-βMan, αMan-(1→2)-αMan, and βGlcNAc(1→2)-αMan] to determine whether the coupling behaviors observed in methyl β-lactoside were more broadly observed. Of the 13 nonconventional J-couplings studied, 7 exhibit properties that may be useful in future MA’AT modeling of O-glycosidic linkages, none of which involve coupling pathways that include the linkage C–O bonds. The findings also provide new insights into the general effects of exocyclic C–O bond conformation on the magnitude of experimental spin-couplings in saccharides and other hydroxyl-containing molecules.

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

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Effects of substitutional Ga on the phonons of δ -phase Pu based on density functional theory calculations

The delta phase of Pu is stabilized by Ga doping, but the mechanism of this stabilization remains an open question. Density functional theory calculations focused on how Ga doping affects the phonons sheds some light on the phonons' contribution to the stabilization. Here, the calculated phonon modes of Ga-doped delta phase Pu fall into two distinct types: localized, high frequency Ga-dominated phonon modes, and Pu-dominated modes at lower frequencies. Increasing the Ga concentration has an effect on the Pu-dominated phonon modes opposite to that of compression: higher-frequency modes soften, and lower-frequency modes stiffen. The latter provides an indication that the stabilization mechanism is not due to a thermodynamic contribution from the phonons. Furthermore, the stiffened phonon modes include candidate modes that describe possible pathways into low-temperature phases, suggesting that doping with Ga could impede such pathways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Density functional theory calculation of helium segregation and decohesion effect in W110/Ni111 interphase boundary

Density functional theory (DFT) is employed to study the preferential distribution and decohesion effect of He in a W-NiFe composite consisting of W particles embedded in an Ni-based solid solution matrix. A slab containing {110}<100>W//{111}<110>Ni interface is employed as a surrogate model for the W-NiFe system. Firstly, the fracture energy of the W/Ni interphase boundary (IB) (4.37 J/m 2 ) is higher than the cleavage energy of Ni{111} (3.82 J/m^2) and lower than the cleavage energy of W{110} (6.60 J/m 2 ). The comparison shows that the cohesion of the IB is stronger than the Ni{111} planes that are away from the IB. However, the cohesion between the Ni{111} planes adjacent to the IB is found to be the weakest in this system, with a cleavage energy of 3.11 J/m 2 . Subsequently, formation energy of He is calculated in the Ni slab, W slab, and various interstitial sites in the IB. The DFT calculations show that He is significantly more stable in Ni than in W, by about 1.75 eV. Interestingly, He does not prefer to segregate at the IB as compared to bulk Ni. Nevertheless, it prefers to segregate to the region between the Ni{111} planes adjacent to the IB, and decreases the cohesion of the already weakest region. Based on an estimated amount of He gas production in 5 years under first wall neutron irradiation (neutron flux of 1.04x10 15 n/cm 2 /s), the He will decrease the cleavage energy of the weakest region by 21.2%, 15.4%, and 12.2% at 800, 1000, and 1200 C, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Analysis and visualization of energy densities. II. Insights from linear-response time-dependent density functional theory calculations

Inspired by the analysis of Kohn–Sham energy densities by Nakai and coworkers, we extended the energy density analysis to linear-response time-dependent density functional theory (LR-TDDFT) calculations. Using ethylene–tetrafluoroethylene and oxyluciferin–water complexes as examples, distinctive distribution patterns were demonstrated for the excitation energy densities of local excitations (within a molecular fragment) and charge-transfer excitations (between molecular fragments). It also provided a simple way to compute the effective energy of both hot carriers (particle and hole) from charge-transfer excitations via an integration of the excitation energy density over the donor and acceptor grid points.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
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 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↗