Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional 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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Role of Eu-Doping in the Electron Transport Behavior in the Zintl Thermoelectric Ca 5-x-y Yb x Eu y Al 2 Sb 6 System

A series of Eu-doped Zintl compounds belonging to theCa 5-x-y Yb x Eu y Al 2 Sb 6 (x = 0, 1.12; 0 ≤ y ≤ 0.63(2)) system have been successfully synthesized by both the arc-melting and the molten Pb-flux methods. All of the five title compounds initially crystallized in the Ca 5 Ga 2 As 6 -type phase (space group Pbam, Z = 2, Pearson code oP26) and maintained their original structure even after the post-heat treatment, unlike the recently reported n-type Zintl analogues in the Ca 5-x-y Yb x RE y Al 2 Sb 6 (RE = Pr, Nd, Sm) systems, which underwent a phase transition from the Ca 5 Ga 2 As 6 -type to the Ca 5 Al 2 Bi 6 -type phase after annealing. This research aimed to understand the origin of the structural preference of the title Ca 5-x-y Yb x Eu y Al 2 Sb 6 system, whether it was affected by the valence electron count or the cationic size. Electrical transport property measurements showed an increase in electrical conductivities and a decrease of Seebeck coefficients for Ca 4.89(1) Eu 0.11 Al 2 Sb 6 , Ca 4.82(1) Eu 0.18 Al 2 Sb 6 , and Ca 4.62(1) Eu 0.38 Al 2 Sb 6 , compared to the parental compound Ca 5 Al 2 Sb 6 . Hole effect measurements proved that these changes should be attributed to the reduced carrier concentration and enhanced carrier mobility. The comprehensive density functional theory calculations including electron density map analysis for the hypothetical model Ca 4.5 Eu 0.5 Al 2 Sb 6 revealed that the polarity between Al and Sb forming the anionic frameworks decreased as the Eu-dopants were introduced, which eventually affected the carrier mobility in the anionic frameworks. Thermal conductivity measurements proved that the Eu-doping successfully lowered the lattice thermal conductivity because of the enhanced atomic disordering. In conclusion, the magnetization measurements for Ca 4.37(2) Eu 0.63 Al 2 Sb 6 showed a typical Curie–Weiss behavior with weak antiferromagnetic nearest-neighbor interactions with θ p = -5.07 K.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Application of machine learning to discover new intermetallic catalysts for the hydrogen evolution and the oxygen reduction reactions

The adsorption energies for hydrogen, oxygen, and hydroxyl were calculated by means of density functional theory on the lowest energy surface of 24 pure metals and 332 binary intermetallic compounds with stoichiometries AB, A 2 B, and A 3 B taking into account the effect of biaxial elastic strains. This information was used to train two random forest regression models, one for the hydrogen adsorption and another for the oxygen and hydroxyl adsorption, based on 9 descriptors that characterized the geometrical and chemical features of the adsorption site as well as the applied strain. All the descriptors for each compound in the models could be obtained from physico-chemical databases. The random forest models were used to predict the adsorption energy for hydrogen, oxygen, and hydroxyl of ≈2700 binary intermetallic compounds with stoichiometries AB, A 2 B, and A 3 B made of metallic elements, excluding those that were environmentally hazardous, radioactive, or toxic. This information was used to search for potential good catalysts for the HER and ORR from the criteria that their adsorption energy for H and O/OH, respectively, should be close to that of Pt. Further, this investigation shows that the suitably trained machine learning models can predict adsorption energies with an accuracy not far away from density functional theory calculations with minimum computational cost from descriptors that are readily available in physico-chemical databases for any compound. Moreover, the strategy presented in this paper can be easily extended to other compounds and catalytic reactions, and is expected to foster the use of ML methods in catalysis.

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

Nuclear quantum effects of metal surface-mediated C–H activation

The nuclear quantum effects of surface-mediated C–H activation of surface CH 3 are considered for the pristine Pt(111) and Au(111) surfaces at 300 K. The kinetic barriers without nuclear quantum effects are calculated using both static density functional theory calculations and ab initio molecular dynamics. Static calculations are performed using the harmonic approximation while the free energy pathway is calculated using enhanced sampling molecular dynamics. Machine learning potentials are trained using generated datasets and validated against the ab initio molecular dynamics generated free energy pathways. The machine learning potentials are used to perform centroid molecular dynamics to consider the nuclear quantum effects of C–H activation. Nuclear quantum effects are found to have a very significant effect on the free energy pathway, with reduced importance at higher temperatures and in the CD 3 case.

Bunting, Rhys J. [Lawrence Livermore National Labo↗
DOE OSTI2021

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

Coexisting structural disorder and robust spin-polarization in half-metallic FeMnVAl

We report half-metallic ferromagnets (HMFs) are one of the most promising materials in the field of spintronics due to their unique band structure, consisting of one spin subband having metallic characteristics along with another subband with semiconductor-like behavior. In this work, we report the synthesis of a novel quaternary Heusler alloy—FeMnVAl—and we have studied the structural, magnetic, transport, and electronic properties complemented with first-principles calculations. Among different possible structurally ordered arrangements, the optimal structure is identified by theoretical energy minimization. The corresponding spin-polarized band-structure calculations indicate the presence of a half-metallic ferromagnetic ground state. A detailed and careful investigation of the x-ray diffraction data, Mössbauer spectra, and nuclear magnetic resonance spectra suggests the presence of site disorder between the Fe and Mn atoms in the stable ordered structure of the system. The magnetic susceptibility measurement clearly establishes a ferromagnetic-like transition below ~213 K. The Fe Mössbauer spectrometry measurements suggest that only the Mn spins could be responsible for the magnetic order, which is consistent with our theoretical calculation. Surprisingly, the density-functional-theory calculations reveal that the spin-polarization value is almost immunized (92.4% → 90.4%) from the Mn-Fe structural disorder, even when nonmagnetic Fe and moment-carrying Mn sites are entangled inseparably. The robustness of spin polarization and half-metallicity in the studied FeMnVAl compound comprising structural disorder is thus quite interesting and could provide a new direction to investigate and understand the exact role of disorders on spin polarization in this class of materials, beyond the available knowledge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Effect of the electronic charge gap on LO bond-stretching phonons in undoped La 2 CuO 4 calculated using LDA + U

Typical density-functional theory calculations that wrongly predict undoped cuprates to be metallic also predict Cu-O half- and full-breathing phonon energies that are significantly softer than observed, presumably because of weak on-site Coulomb repulsion on the Cu 3d orbitals. We used DFT+U calculations with antiferromagnetic supercells of La 2 CuO 4 to establish correlation between the on-site repulsion strength, tuned via adjusting the value of U, and phonon dispersions. Here, we find that breathing and half-breathing phonons reach experimental values when U is tuned to obtain the correct optical gap and magnetic moments. We demonstrate that using distorted supercells within DFT+U is a promising framework to model phonons in undoped cuprates and other perovskite oxides with complex, interrelated structural and electronic degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Al rich (111) and (110) surfaces of LaAlO 3

Lanthanum aluminate (LaAlO 3 ) has garnered significant interests since the discovery of a 2D electron gas at the LaAlO 3 /SrTiO 3 interface. Here in this work, the (110) surface of LaAlO 3 was studied using a combination of transmission electron microscopy plan and profile imaging, image simulation, electron diffraction, X-ray photoelectron spectroscopy, and density functional theory calculations. Profile imaging of the (110) surface was performed on a (111) oriented single crystalline sample with an Al rich surface. Our results indicate that high annealing temperatures (1000–1300°C) commonly used in preparing well ordered, flat surfaces of LaAlO 3 result in Al rich surfaces. The (110) surface was found to be single layered AlO x with a (2 × 1) reconstruction. Angle resolved X-ray photoelectron spectroscopy indicated that the (2 × 1) surface structure was not hydrated. Convex hull construction of surface enthalpies as a function of excess AlO 1.5 on the (110) surface was done using density functional theory calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Tight-binding approximations to time-dependent density functional theory — A fast approach for the calculation of electronically excited states

We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.

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

G-centers in irradiated silicon revisited: A screened hybrid density functional theory approach

Electronic structure calculations employing screened hybrid density functional theory are used to gain fundamental insight into the interaction of carbon interstitial (C{sub i}) and substitutional (C{sub s}) atoms forming the C{sub i}C{sub s} defect known as G-center in silicon (Si). The G-center is one of the most important radiation related defects in Czochralski grown Si. We systematically investigate the density of states and formation energy for different types of C{sub i}C{sub s} defects with respect to the Fermi energy for all possible charge states. Prevalence of the neutral state for the C-type defect is established.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Stoichiometry control and electronic and transport properties of pyrochlore B i 2 I r 2 O 7 thin films

Synthesizing stoichiometric and epitaxial thin films of pyrochlore iridates is an essential step toward the experimental realization of unusual topological and magnetic states that are theoretically predicted in this unique spin-orbit coupled material system. Here, we report on the stoichiometry control and electronic and transport properties of pyrochlore iridate Bi 2 Ir 2 O 7 thin films grown by pulsed laser deposition. The as-grown films form a bilayerlike structure, in which the top surface is highly Ir deficient while the bottom layer is mainly composed of iridium metal. By postannealing the as-deposited films in Ir O 2 + O 2 atmosphere, we improved the stoichiometry and homogeneity through the film thickness with the lattice constant close to the bulk value. Density functional theory calculation in the bulk limit shows a fourfold degenerate Dirac node slightly below the Fermi energy at the X point, along with trivial bands around the Γ point. The projected partial density of states suggests that the states in the vicinity of the Fermi energy (–3 to 0 eV) mainly consist of highly hybridized Ir 5d and O 2p with minor contributions from Bi 6s and 6p,while those far below the Fermi energy (–9 to –3 eV) are contributed primarily by the O bands. Transport measurements revealed a weakly metallic behavior at higher temperatures transitioning to a weakly insulating behavior below 150 K, and a low-temperature magnetoresistance qualitatively ascribed to multicarrier and band-structural effects. Here, the transport features are influenced by a density of states sharply peaked at the Fermi energy, and by the coexistence of trivial bands with the Dirac node, as revealed by the density functional theory calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2016

First-principles study of inert gas incorporation and migration in zirconium nitride

Uranium-Molybdenum alloy based fuels are the most promising fuel for the future high performance research reactors using low enriched uranium (LEU). However, when used in dispersion with aluminum, significant amount of interaction products form at the interface between the UMo fuel particle and the surrounding Al during in-pile tests which leads to non-acceptable swelling. To reduce the swelling in U-Mo/Al dispersion fuel at high burnups, different diffusion barrier coatings, e.g., ZrN, were deposited on the U-Mo particles to avoid the interaction between fuel particles and the Al matrix. These diffusion barrier coatings are very effective in reducing the formation of interaction layer and therefore the fuel swelling. So far, ZrN is selected as the candidate coating material for the European Dispersion Fuel program due to the superior thermophysical properties of ZrN. The thermodynamic and kinetic stability of ZrN coating with respect to Al and U-Mo has been studied by the current authors using density functional theory calculations. So far there are very few studies of the diffusion behaviors of inert gases in ZrN. A recent publication on irradiated (Pu,Zr)N reveals that the majority of fission-induced Xe is retained in the material, while He is release to the fuel pin to a large extent. This raises the questions on the underlying differences of inert gas behavior in nitride fuel. Therefore, the purpose of this study is to understand the incorporation and migration behavior of inert gases He, Kr and Xe in ZrN using density functional theory (DFT) based first-principles calculations. Using density-functional theory calculations we provided a systematic study of the incorporation and migration of fission gases in ZrN. We investigated the formation energies of several defects (vacancies, interstitials, divacancies, Frenkel pairs, and Schottky defects) and the incorporation energies of inert gases in these defects. Results show that the formation energies of point defects are strongly affected by the chemical potentials of N and Zr species. The energetically most favorable incorporation site for fission gases is found to be the first nearest neighbor Schottky defect, in which the fission gas atom sits at somewhere between the Zr and N vacancies. The predicted solution energies of rare gases, i.e, He, Kr and Xe, in ZrN, is almost linearly dependent on the atom size of inert gases. We compared the solution energy Xe in ZrN with several other nuclear materials with the same B1 structure (MX), e.g., TiN, UN and UC. The solution energies of Xe in MX are closely related to the size of vacancy sites. The potential migration mechanisms of Xe in ZrN include interstitial and vacancy-assisted diffusion. The calculations show that vacancies essentially act as Xe trapping sites. The migration behavior of Xe in N and Zr divacancy and Schottky defects were compared with that through tetrahedral interstitials. It is found that the migration barrier is much lower with the interstitial mechanism than those with the vacancy mechanism. Therefore, at low temperatures Xe diffusion in ZrN is mostly through the interstitial mechanism or existing grain boundaries while at high temperature volume diffusion could be preferred, i.e., diffusion with Zr and N vacancies. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Calculation of the Physicochemical Characteristics of a New Orthorhombic Form of Diamond

A geometrically optimized structure and some properties of a new diamond polymorph have been assessed using density-functional theory calculations in the local-density approximation. The new phase has a primitive orthorhombic lattice (Pmma). The calculation results indicate that the density, sublimation energy, bulk modulus, Vickers hardness, and band gap of this phase are lower than those of cubic diamond by 5–18%. The orthorhombic diamond polymorph can be obtained by compressing graphite in the [001] direction at a pressure of 59.4 GPa. For this phase transition, we calculated theoretical X-ray diffraction patterns of the new phase and intermediate states, which can be used for their identification in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Electronic state spectroscopy of diiodomethane (CH{sub 2}I{sub 2}): Experimental and computational studies in the 30 000–95 000 cm{sup −1} region

The electronic absorption spectrum of diiodomethane in the 30 000–95 000 cm{sup −1} region is investigated using synchrotron radiation; the spectrum in the 50 000–66 500 cm{sup −1} region is reported for the first time. The absorption bands in the 30 000–50 000 cm{sup −1} region are attributed to valence transitions, while the vacuum ultraviolet (VUV) spectrum (50 000–95 000 cm{sup −1}) is dominated by several Rydberg series converging to the first four ionization potentials of CH{sub 2}I{sub 2} at 9.46, 9.76, 10.21, and 10.56 eV corresponding to the removal of an electron from the outermost 3b{sub 2}, 2b{sub 1}, 1a{sub 2}, and 4a{sub 1} non-bonding orbitals, respectively. Rydberg series of ns, np, and nd type converging to each of the four ionization potentials are assigned based on a quantum defect analysis. Time dependent density functional theory calculations of excited states support the analysis and help in interpretation of the Rydberg and valence nature of observed transitions. Density functional theory calculations of the neutral and ionic ground state geometries and vibrational frequencies are used to assign the observed vibronic structure. Vibronic features accompanying the Rydberg series are mainly due to excitation of the C-I symmetric stretch (ν{sub 3}) and CH{sub 2} wag (ν{sub 8}) modes, with smaller contributions from the C-H symmetric stretch (ν{sub 1}). UV absorption bands are assigned to low lying valence states 1{sup 1}B{sub 2}, 1{sup 1}B{sub 1}, 2{sup 1}A{sub 1}, 3{sup 1}A{sub 1}, 2{sup 1}B{sub 1}, and 2{sup 1}B{sub 2} and the unusually high underlying intensity in parts of the VUV spectrum is attributed to valence states with high oscillator strength. This is the first report of a comprehensive Rydberg series and vibronic analysis of the VUV absorption spectrum of CH{sub 2}I{sub 2} in the 50 000–85 000 cm{sup −1} region. The VUV absorption spectrum of CD{sub 2}I{sub 2} which serves to verify and consolidate spectral assignments is also reported here for the first time.

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

First principles calculation of lattice thermal conductivity of metals considering phonon-phonon and phonon-electron scattering

The effect of phonon-electron (p-e) scattering on lattice thermal conductivity is investigated for Cu, Ag, Au, Al, Pt, and Ni. We evaluate both phonon-phonon (p-p) and p-e scattering rates from first principles and calculate the lattice thermal conductivity (κ{sub L}). It is found that p-e scattering plays an important role in determining the κ{sub L} of Pt and Ni at room temperature, while it has negligible effect on the κ{sub L} of Cu, Ag, Au, and Al. Specifically, the room temperature κ{sub L}s of Cu, Ag, Au, and Al predicted from density-functional theory calculations with the local density approximation are 16.9, 5.2, 2.6, and 5.8 W/m K, respectively, when only p-p scattering is considered, while it is almost unchanged when p-e scattering is also taken into account. However, the κ{sub L} of Pt and Ni is reduced from 7.1 and 33.2 W/m K to 5.8 and 23.2 W/m K by p-e scattering. Even though Al has quite high electron-phonon coupling constant, a quantity that characterizes the rate of heat transfer from hot electrons to cold phonons in the two-temperature model, p-e scattering is not effective in reducing κ{sub L} owing to the relatively low p-e scattering rates in Al. The difference in the strength of p-e scattering in different metals can be qualitatively understood by checking the amount of electron density of states that is overlapped with the Fermi window. Moreover, κ{sub L} is found to be comparable to the electronic thermal conductivity in Ni.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Atomically Precise PdSe 2 Pentagonal Nanoribbons

We report atomically precise pentagonal PdSe 2 nanoribbons (PNRs) fabricated on a pristine PdSe 2 substrate with a hybrid method of top-down and bottom-up processes. The PNRs form a uniform array of dimer structure with a width of 2.4 nm and length of more than 200 nm. In situ four-probe scanning tunneling microscopy (STM) reveals metallic behavior of PNRs with ballistic transport for at least 20 nm in length. Density functional theory calculations produce a semiconducting density of states of isolated PNRs and find that the band gap narrows and disappears quickly once considering coupling between PNR stacking layers or interaction with the PdSe 2 substrate. The coupling of PNRs is further corroborated by Raman spectroscopy and field-effect transistor measurements. The facile method of fabricating atomically precise PNRs offers an air-stable functional material for dimensional control.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗