Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT calculation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Half-Metallic and Half-Semiconductor Gaps in Cr-Based Chalcogenides: DFT + U Calculations

Using full-potential linear-augmented plane waves plus local orbitals (FP-LAPW) method, within the generalized gradient approximation (GGA) and GGA + U (Hubbard Coulomb onsite correction) as well as the modified Becke-Johnson approach (mBJ-GGA and mBJ-GGA + U), we have carried out a systematic investigation of the structural, electronic, and magnetic properties of CrS, CrSe, and CrTe compounds in two competing phases: the zinc-blende (ZB) and wurtzite (WZ) structures. We obtained half-metallic (HM) gaps of 1.35 (2.32), 1.65 (2.49), and 1.52 (1.80) for the CrS, CrSe, and CrTe compounds, respectively, in the ZB phase within the GGA + U (mBJ-GGA + U) approach. The HM gaps are 1.78, 2.11, and 1.53 for CrS, CrSe, and CrTe compounds, respectively, in the WZ phase using the GGA + U. Notably, the CrSe and CrTe using the mBJ-GGA + U approach are a half-semiconductor (HSC) in the WZ phase. The HSC gaps are 0.94 and 0.25 eV for CrSe and CrTe compounds, respectively, while CrS has a HM gap of 2.30 eV. We observed a HSC-to-HM transition for the CrSe and CrTe crystal structures on application of strain. The CrS, CrSe, and CrTe compounds have a total magnetic moment of 4.0 µ{sub B} per formula unit with the magnetization predominantly from the Cr ion. The main features observed from the density of states and strain evolution should motivate further experimental exploration on the possible application of the CrS, CrSe, and CrTe compounds for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Elucidating the Impact of Cis – Trans Organic Structure Directing Agent Isomer Ratios on the Aluminum Distribution Within SSZ-39

Despite their widespread use, the mechanisms governing the synthesis of zeolite catalysts are still poorly understood. A notable example of this problem is the uncertainty surrounding the influence of synthesis conditions on the placement of Al atoms in the zeolite framework which determines the active sites available for catalytic species. In this work, the role of the cis to trans isomer ratio of the OSDA N,N-dimethyl-3-5-dimethylpiperidinium on the energetics of 26 distinct Al pair distributions in SSZ-39 is examined both in the presence and absence of Na using density functional theory calculations. The initial orientation of the OSDA was found to have a significant impact on the final energies present, necessitating the screening of a large number of initial orientations with force field calculations and single point DFT calculations. Ground state energies were found to vary significantly with the ratio of cis to trans OSDAs with a Boltzmann distribution revealing the most likely Al pair distributions shift from sharing the same 8 membered rings to sharing the same double six membered rings to having no shared subunits as one increases the amount of cis OSDA present within the framework. The presence of Na was found to favor Al pair distributions where both Als occupied the same 6-membered ring. When an implicit solvent model was used to evaluate ground state energies the ideal Na sites shifted from 6-membered rings to empty SSZ-39 cages while OSDA positions and orientations remained largely the same. To provide insight on how kinetic factors may influence Al distributions, formation energies we calculated for connected double six membered rings. Further, these formation energies revealed a preference for Al pairs to occupy the same 4-membered ring which indicates kinetic and thermodynamic control may lead to different Al distributions in SSZ-39.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2018

Noncollinear Relativistic DFT + U Calculations of Actinide Dioxide Surfaces

A noncollinear relativistic PBEsol + U study of low-index actinide dioxides (AnO 2 , An = U, Np, or Pu) surfaces has been conducted. The importance of magnetic vector reorientation relative to the plane of the surface is highlighted; this has often been ignored in collinear nonrelativistic models. The use of noncollinear relativistic methods is key to the design of reliable computational models. The ionic relaxation of each surface is shown to be confined to the first three monolayers, and we have explored the configurations of the terminal oxygen ions on the reconstructed (001) surface. The reconstructed (001) surfaces are ordered as (001)αβ < (001)α < (001)β in terms of energetics. Electrostatic potential isosurface and scanning tunneling microscopy images have also been calculated. In conclusion, by considering the energetics of the low-index AnO 2 surfaces, an octahedral Wulff crystal morphology has been calculated.

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

Combined hybrid functional and DFT+U calculations for metal chalcogenides

In the density-functional studies of materials with localized electronic states, the local/semilocal exchange-correlation functionals are often either combined with a Hubbard parameter U as in the LDA+U method or mixed with a fraction of exactly computed (Fock) exchange energy yielding a hybrid functional. Although some inaccuracies of the semilocal density approximations are thus fixed to a certain extent, the improvements are not sufficient to make the predictions agree with the experimental data. Here, we put forward the perspective that the hybrid functional scheme and the LDA+U method should be treated as complementary, and propose to combine the range-separated Heyd-Scuseria-Ernzerhof (HSE) hybrid functional with the Hubbard U. We thus present a variety of HSE+U calculations for a set of II-VI semiconductors, consisting of zinc and cadmium monochalcogenides, along with comparison to the experimental data. Our findings imply that an optimal value U{sup *} of the Hubbard parameter could be determined, which ensures that the HSE+U{sup *} calculation reproduces the experimental band gap. It is shown that an improved description not only of the electronic structure but also of the crystal structure and energetics is obtained by adding the U{sup *} term to the HSE functional, proving the utility of HSE+U{sup *} approach in modeling semiconductors with localized electronic states.

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

Electronic structure of ferromagnetic heavy fermion, YbPdSi, YbPdGe, and YbPtGe studied by photoelectron spectroscopy, x-ray emission spectroscopy, and DFT + DMFT calculations

Here, the electronic structures of ferromagnetic heavy fermion Yb compounds of YbPdSi, YbPdGe, and YbPtGe are studied by photoelectron spectroscopy around the Yb 4d–4f resonance, resonant x-ray emission spectroscopy at the Yb L 3 absorption edge, and density functional theory combined with dynamical mean field theory calculations. These compounds all have a temperature-independent intermediate Yb valence with large ${\rm Yb}^{3+}$ and small ${\rm Yb}^{2+}$ components. The magnitude of the Yb valence is evaluated to be YbPtGe $<$ YbPdGe $\lesssim $ YbPdSi, suggesting that YbPtGe is the closest to the quantum critical point among the three Yb compounds. Our results support the scenario of the coexistence of heavy fermion behavior and ferromagnetic ordering which is described by a magnetically-ordered Kondo lattice where the magnitude of the Kondo effect and the RKKY interaction are comparable.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2019

DFT + DMFT calculations of the complex band and tunneling behavior for the transition metal monoxides MnO, FeO, CoO, and NiO

Here, we report complex band structure (CBS) calculations for the four late transition metal monoxides MnO, FeO, CoO, and NiO in their paramagnetic phase. The CBS is obtained from density functional theory plus dynamical mean field theory (DMFT) calculations to take into account correlation effects. The so-called β parameters, governing the exponential decay of the transmission probability in the nonresonant tunneling regime of these oxides, are extracted from the CBS. Different model constructions are examined in the DMFT part of the calculation. The calculated β parameters provide theoretical estimation for the decay length in the evanescent channel, which would be useful for tunnel junction applications of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

de Haas–van Alphen study of role of 4 f electrons in antiferromagnetic CeZn 11 as compared to its nonmagnetic analog LaZn 11

Here we present a de Haas–van Alphen study of the Fermi surface of the low-temperature antiferromagnet CeZn 11 and its nonmagnetic analog LaZn 11 , measured by torque magnetometry up to fields of 33T and at temperatures down to 320 mK . Both systems possess similar de Haas–van Alphen frequencies, with three clear sets of features—ranging from 50 T to 4 kT —corresponding to three bands of a complex Fermi surface, with an expected fourth band also seen weakly in CeZn 11 . The effective masses of the charge carriers are very light (<1 m e ) in LaZn 11 but a factor of 2–4 larger in CeZn 11 , indicative of stronger electronic correlations. We perform detailed density functional theory (DFT) calculations for CeZn 11 and find that only DFT+ U calculations with U = 1.5 eV , which localize the 4 f states, provide a good match to the measured de Haas–van Alphen frequencies, once the presence of magnetic breakdown orbits is also considered. Finally, our study suggests that the Fermi surface of CeZn 11 is very close to that of LaZn 11 being dominated by Zn 3d , as the Ce 4 f states are localized and have little influence on its electronic structure, however, they are responsible for its magnetic order and contribute to enhance electronic correlations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Recent research on stishovite: Hugoniot and partial release Z experiments and DFT EOS calculations

We have conducted a series of ride-along experiments on the Z facility to ascertain the Hugoniot of silica centered in the stishovite phase over a range 0.4 - 1.0 TPa, together with partial release states produced at the interface between the sample and a fused silica window. The stishovite samples were synthesized in a large-volume multi-anvil press at 15 GPa and 1773 K, with an initial density of 4.29 gm/cc. The new Z experiments on stishovite fill in a gap between gas gun experiments and NIF experiments. The states are compared with the Hugoniots of quartz and fused silica for inferences as to EOS. They are generally consistent with Sesame 7360 predictions. Sound speed constraints from these data are discussed. The new Hugoniot data cross over the melting curve of stishovite; together with the partial-release data and predictions from density-functional theory modeling, they provide insights into the properties of solid and liquid under extreme conditions. These data are fundamentally important for understanding the interior of silicate-based super-Earths.

58 GEOSCIENCES↗
DOE OSTI2023

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
DOE OSTI2025

Kohn anomalies and phonon anharmonicity in iridium

Elemental iridium presents surprising challenges for both inelastic neutron scattering (INS) and theoretical thermal transport calculations due to its high neutron absorption cross-section and strong electron-phonon interactions, respectively. Here, in this study, we overcome these challenges to measure temperature-dependent phonon dispersion curves, compare these with calculations based on density functional theory (DFT), and ultimately examine the electron-phonon limited transport behaviors of this material. Our DFT calculations demonstrate Kohn anomalies, near the 𝐾 point of the iridium Brillouin zone, indicating coupling between electrons and phonons. Strong electron-phonon coupling can compete with anharmonic effects to determine electrical and thermal transport behaviors and make the Kohn anomalies challenging to observe. Nonetheless, our INS measurements map these anomalies and other dispersion features over the Brillouin zone from 100 to 700 K. These measurements also uncover unexpectedly large mode specific Grüneisen parameters obtained from the temperature-dependent phonon energies, highlighting strong anharmonicity in iridium. DFT-based Boltzmann transport calculations demonstrate how anharmonicity and electron-phonon couplings determine electronic and lattice transport behaviors. Furthermore, we correlate the Kohn anomalies with calculated electron-phonon nesting functions, Fermi surfaces, and DFT-derived coupling strengths. This study provides detailed insights into the temperature-dependent mode-resolved lattice dynamics and anharmonicity, transport behavior, and electron-phonon interactions.

DFT↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2020

Bypassing the Energy Functional in Density Functional Theory: Direct Calculation of Electronic Energies from Conditional Probability Densities

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP DFT calculations can overcome the convergence problems of Kohn-Sham DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗