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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Quantitative Attachment of Bimetal Combinations of Transition-Metal Ions to the Surface of TiO 2 Nanorods

Here we report the sequential, quantitative loading of transition-metal ions (Cr 3+ , Mn 2+ , Fe 2+ , Co 2+ , Ni 2+ , and Cu 2+ ) onto the surface of rod-shaped anatase TiO 2 nanocrystals in bimetallic combinations ( 6 C 2 = 15) to form M,M'-TiO 2 nanocrystals. The materials were characterized with transmission electron microscopy (TEM), powder X-ray diffraction (XRD), elemental analysis, X-ray photoelectron spectroscopy (XPS), and UV–visible spectroscopy. TEM and XRD data indicate that the sequential adsorption of metal ions occurs with the retention of the phase and morphology of the nanocrystal. Atomistic models of the M,M'-TiO 2 nanocrystals were optimized with density functional theory calculations. Calculated UV–visible absorption spectra and partial charge density maps of the donor and acceptor states for the electronic transitions indicate the importance of metal-to-metal charge transfer (MMCT) processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structures and localized vibrational states of defects in graphite by tight-binding calculations

The structural and vibrational properties of pristine graphite and point defects in graphite are studied by tight-binding (TB) calculations using a three-center TB potential model. We showed that the three-center TB potential without ‘‘ad hoc’’ van der Waals interaction corrections can accurately describe the inter-layer distance of graphite and the lowest-energy structures and stabilities of typical point defects in graphite. The results from our TB calculations are in good agreement with those from density-functional theory calculations with van der Waals interaction corrections. We also investigated the vibrational properties to gain better understanding on the localization of vibrational states induced by the point defects. Our calculation results show that although localized or quasi-localized vibrational modes can be found in all defected graphite, the localization induced by Frenkel pair, dual-vacancy, and dual-interstitial defects is much stronger. Atomic displacements associated with the localized vibrational modes induced by these three point defects are also analyzed

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

Magnetic properties of the quasi-one-dimensional S = 1 spin chain antiferromagnet BaNi Te 2 O 7

We report a quasi-one-dimensional S=1 spin chain compound BaNiTe 2 O 7 . This magnetic system has been investigated by magnetic susceptibility, specific heat, and neutron powder diffraction. These results indicate that BaNiTe 2 O 7 develops a short-range magnetic correlation around T ~ 22K. With further cooling, an antiferromagnetic phase transition is observed at T N ~ 5.4K. Neutron powder diffraction revealed antiferromagnetic noncollinear order with a commensurate propagation vector k = (1/2,1,0). The refined magnetic moment size of Ni 2+ at 1.5 K is 1.84 μB , and its noncollinear spin texture is confirmed by first-principles calculations. In conclusion, inelastic neutron-scattering results and density functional theory calculations confirmed the quasi-one-dimensional nature of the spin systems.

1-dimensional spin chains↗
DOE OSTI2013

Atomistic force field for alumina fit to density functional theory

We present a force field for bulk alumina (Al{sub 2}O{sub 3}), which has been parametrized by fitting the energies, forces, and stresses of a large database of reference configurations to those calculated with density functional theory (DFT). We use a functional form that is simpler and computationally more efficient than some existing models of alumina parametrized by a similar technique. Nevertheless, we demonstrate an accuracy of our potential that is comparable to those existing models and to DFT. We present calculations of crystal structures and energies, elastic constants, phonon spectra, thermal expansion, and point defect formation energies.

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

Towards an Optimal Gradient-dependent Energy Functional of the PZ-SIC Form

Results of Perdew–Zunger self-interaction corrected (PZ-SIC) density functional theory calculations of the atomization energy of 35 molecules are compared to those of high-level quantum chemistry calculations. While the PBE functional, which is commonly used in calculations of condensed matter, is known to predict on average too high atomization energy (overbinding of the molecules), the application of PZ-SIC gives a large overcorrection and leads to significant underestimation of the atomization energy. The exchange enhancement factor that is optimal for the generalized gradient approximation within the Kohn-Sham (KS) approach may not be optimal for the self-interaction corrected functional. The PBEsol functional, where the exchange enhancement factor was optimized for solids, gives poor results for molecules in KS but turns out to work better than PBE in PZ-SIC calculations. The exchange enhancement is weaker in PBEsol and the functional is closer to the local density approximation. Furthermore, the drop in the exchange enhancement factor for increasing reduced gradient in the PW91 functional gives more accurate results than the plateaued enhancement in the PBE functional. A step towards an optimal exchange enhancement factor for a gradient dependent functional of the PZ-SIC form is taken by constructing an exchange enhancement factor that mimics PBEsol for small values of the reduced gradient, and PW91 for large values. The average atomization energy is then in closer agreement with the high-level quantum chemistry calculations, but the variance is still large, the F 2 molecule being a notable outlier.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

First principles, microkinetic, and experimental analysis of Lewis acid site speciation during ethanol dehydration on Sn-Beta zeolites

Density functional theory calculations are combined with kinetic measurements of ethanol dehydration to diethyl ether to identify the relative catalytic contributions of structurally distinct speciations of Sn sites in zeolite Beta frameworks. The structural complexities of the Beta framework require nonstandard techniques for entropy and energy calculations, including consideration of anharmonic effects in vibrational modes, employment of quasi-harmonic densities of states methods to evaluate entropies, and use of hybrid density functionals to evaluate binding energies. Calculated energies and entropies are used to construct a microkinetic model that is iteratively refined to identify all kinetically and thermodynamically sensitive reaction steps and intermediates which are subsequently treated with the higher-level methods. The rate and equilibrium constants obtained from this tiered approach agree well with measured reaction orders in ethanol and water. Site balances provide evidence for the interconversion of Sn sites between “closed” configurations that are tetra-coordinated to the framework, open configurations formed by hydrolysis that are tri-coordinated to the framework and contain a hydroxyl ligand and proximal silanol group (“hydroxy-open”), and open tri-coordinate configurations formed by reaction with ethanol, yielding an ethoxy ligand and proximal silanol group (“ethoxy-open”). Closed, hydroxy-open, and ethoxy-open Sn sites adsorb ethanol via distinct modes and react via distinct pathways, with the consequence that prevalent dehydration pathways depend on the speciation of Sn sites under reaction conditions. The kinetic modeling indicates that, under the conditions studied (404 K, 0.5–35 kPa ethanol, 0.1–50 kPa water), bimolecular dehydration on the closed Sn site is the sole kinetically-relevant step, and ethanol, ethanol-ethanol dimers, and ethanol-water dimers are the most abundant surface intermediates. In conclusion, these results highlight the importance of considering the distribution of heteroatom coordination modes to zeolite frameworks under reaction conditions, as well as pathways for their interconversion in the presence of reacting molecules, to obtain more robust mechanistic and kinetic interpretations of catalytic pathways in Lewis acid zeolites.

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

Computation of the Thermal Expansion Coefficient of Graphene with Gaussian Approximation Potentials

Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this work, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

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

Energetics of intrinsic defects in NiO and the consequences for its resistive random access memory performance

Energetics for a variety of intrinsic defects in NiO are calculated using state-of-the-art ab initio hybrid density functional theory calculations. At the O-rich limit, Ni vacancies are the lowest cost defect for all Fermi energies within the gap, in agreement with the well-known p-type behaviour of NiO. However, the ability of the metal electrode in a resistive random access memory metal-oxide-metal setup to shift the oxygen chemical potential towards the O-poor limit results in unusual NiO behaviour and O vacancies dominating at lower Fermi energy levels. Calculated band diagrams show that O vacancies in NiO are positively charged at the operating Fermi energy giving it the advantage of not requiring a scavenger metal layer to maximise drift. Ni and O interstitials are generally found to be higher in energy than the respective vacancies suggesting that significant recombination of O vacancies and interstitials does not take place as proposed in some models of switching behaviour.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Electron Density Determination, Bonding and Properties of Tetragonal Ferromagnetic Intermetallics

The project developed quantitative convergent-beam electron diffraction (QCBED) methods by energy-filtered transmission electron microscopy (EFTEM) and used them in combination with density functional theory (DFT) calculations to study the electron density distribution in metallic and intermetallic phases with different cubic and non-cubic crystal structures that comprise elements with d-electron shells. The experimental methods developed here focus on the bonding charge distribution as one of the quantum mechanical characteristics central for understanding of intrinsic properties and validation of DFT calculations. Multiple structure and temperature factors have been measured simultaneously from nano-scale volumes of high-quality crystal with sufficient accuracy and precision for comparison with electron density distribution calculations by DFT. The often anisotropic temperature factors for the different atoms and atom sites in chemically ordered phases can differ significantly from those known for relevant pure element crystals due to bonding effects. Thus they have been measured from the same crystal volumes from which the structure factors have been determined. The ferromagnetic ordered intermetallic phases FePd and FePt are selected as model systems for 3d-4d and 3d-5d electron interactions, while the intermetallic phases NiAl and TiAl are used to probe 3d-3p electron interactions. Additionally, pure transition metal elements with d-electrons have been studied. FCC metals exhibit well defined delocalized bonding charge in tetrahedral sites, while less directional, more distributed bonding charge attains in BCC metals. Agreement between DFT calculated and QCBED results degrades as d-electron levels fill in the elements, and for intermetallics as d-d interactions become prominent over p-d interactions. Utilizing the LDA+U approach enabled inclusion of onsite Coulomb-repulsion effects in DFT calculations, which can afford improved agreements with QCBED results. This implies that onsite Coulomb repulsion effects become non-negligible as the d-orbitals fill. The use of now easily measured low-order structure factors as an additional experimental metric in validation of DFT calculation of electronic structures of crystals offers potential to capture better both total energy related properties and details of the interatomic bonding in system with d-electron orbital contributions. This effort advanced the state of the art in quantitative TEM experimentation, provides original experimental data uniquely suited for new validation approaches of DFT calculations of d-electron affected transition metals and intermetallics.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A method of orbital analysis for large-scale first-principles simulations

An efficient method of calculating the natural bond orbitals (NBOs) based on a truncation of the entire density matrix of a whole system is presented for large-scale density functional theory calculations. The method recovers an orbital picture for O(N) electronic structure methods which directly evaluate the density matrix without using Kohn-Sham orbitals, thus enabling quantitative analysis of chemical reactions in large-scale systems in the language of localized Lewis-type chemical bonds. With the density matrix calculated by either an exact diagonalization or O(N) method, the computational cost is O(1) for the calculation of NBOs associated with a local region where a chemical reaction takes place. As an illustration of the method, we demonstrate how an electronic structure in a local region of interest can be analyzed by NBOs in a large-scale first-principles molecular dynamics simulation for a liquid electrolyte bulk model (propylene carbonate + LiBF{sub 4})

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

A First-Principles Calculation of Electronic Structure and dHvA Simulation of the Novel Superconductor Re{sub 6}Hf

In this paper, we performed a density functional theory calculation study on the newly discovered superconductor Re{sub 6}Hf containing anti-symmetric spin-orbit coupling (ASOC) properties. The calculated densities of electronic states demonstrate that the d–d interaction of Re{sub 6}Hf is much stronger than that in Re{sub 6}Zr. A fully relativistic DFT study shows a special performance of densities of states for the spin-up and spin-down electrons. Band structure calculations indicate that the lifting at several high symmetric k points and plenty of splits in band structures yield significant topological transitions in the Fermi surfaces of Re{sub 6}Hf. The energy splitting of bands caused by ASOC are estimated to be 34 and 44 meV, respectively. The de Hass-van Alphen effect simulation is demonstrated and analyzed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Mixed cerium-platinum oxides: Electronic structure of [CeO]Pt{sub n} (n = 1, 2) and [CeO{sub 2}]Pt complex anions and neutrals

The electronic structures of several small Ce–Pt oxide complexes were explored using a combination of anion photoelectron (PE) spectroscopy and density functional theory calculations. Pt and Pt{sub 2} both accept electron density from CeO diatomic molecules, in which the cerium atom is in a lower-than-bulk oxidation state (+2 versus bulk +4). Neutral [CeO]Pt and [CeO]Pt{sub 2} complexes are therefore ionic, with electronic structures described qualitatively as [CeO{sup +2}]Pt{sup −2} and [CeO{sup +}]Pt{sub 2}{sup −}, respectively. The associated anions are described qualitatively as [CeO{sup +}]Pt{sup −2} and [CeO{sup +}]Pt{sub 2}{sup −2}, respectively. In both neutrals and anions, the most stable molecular structures determined by calculations feature a distinct CeO moiety, with the positively charged Ce center pointing toward the electron rich Pt or Pt{sub 2} moiety. Spectral simulations based on calculated spectroscopic parameters are in fair agreement with the spectra, validating the computationally determined structures. In contrast, when Pt is coupled with CeO{sub 2}, which has no Ce-localized electrons that can readily be donated to Pt, the anion is described as [CeO{sub 2}]Pt{sup −}. The molecular structure predicted computationally suggests that it is governed by charge-dipole interactions. The neutral [CeO{sub 2}]Pt complex lacks charge-dipole stabilizing interactions, and is predicted to be structurally very different from the anion, featuring a single Pt–O–Ce bridge bond. The PE spectra of several of the complexes exhibit evidence of photodissociation with Pt{sup −} daughter ion formation. The electronic structures of these complexes are related to local interactions in Pt-ceria catalyst-support systems.

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

High-throughput determination of Hubbard $U$ and Hund $J$ values for transition metal oxides via the linear response formalism

DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 1000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an ATOMATE workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. Here, to demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-pU value results in a significantly improved agreement between the computed lattice parameters and experiment

36 MATERIALS SCIENCE↗
DOE OSTI2015

Atomic-scale properties of Ni-based FCC ternary, and quaternary alloys

The aim of our study is to characterize some atomic-scale properties of Ni-based FCC multicomponent alloys. For this purpose, we use Monte Carlo method combined with density functional theory calculations to study short-range order (SRO), atomic displacements, electronic density of states, and magnetic moments in equimolar ternary NiCrCo, and quaternary NiCrCoFe alloys. The salient features for the ternary alloy are a negative SRO parameter between Ni Cr and a positive between Cr Cr pairs as well as a weakly magnetic state. For the quaternary alloy we predict negative SRO parameter for Ni Cr and Ni Fe pairs and positive for Cr Cr and Fe Fe pairs. Atomic displacements for both ternary and quaternary alloys are negligible. In contrast to the ternary, the quaternary alloy shows a complex magnetic structure. The electronic structure of the ternary and quaternary alloys shows differences near the Fermi energy between a random solid solution and the predicted structure with SRO. Despite that, the calculated EXAFS spectra does not show enough contrast to discriminate between random and ordered structures. Finally, the predicted SRO has an impact on point-defect energetics, electron phonon coupling and thermodynamic functions and thus, SRO should not be neglected when studying properties of these two alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Predicted boron-carbide compounds: A first-principles study

By using developed particle swarm optimization algorithm on crystal structural prediction, we have explored the possible crystal structures of B-C system. Their structures, stability, elastic properties, electronic structure, and chemical bonding have been investigated by first-principles calculations with density functional theory. The results show that all the predicted structures are mechanically and dynamically stable. An analysis of calculated enthalpy with pressure indicates that increasing of boron content will increase the stability of boron carbides under low pressure. Moreover, the boron carbides with rich carbon content become more stable under high pressure. The negative formation energy of predicted B{sub 5}C indicates its high stability. The density of states of B{sub 5}C show that it is p-type semiconducting. The calculated theoretical Vickers hardnesses of B-C exceed 40 GPa except B{sub 4}C, BC, and BC{sub 4}, indicating they are potential superhard materials. An analysis of Debye temperature and electronic localization function provides further understanding chemical and physical properties of boron carbide.

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

Quantum Monte Carlo calculation of the binding energy of the beryllium dimer

The accurate calculation of the binding energy of the beryllium dimer is a challenging theoretical problem. In this study, the binding energy of Be{sub 2} is calculated using the diffusion Monte Carlo (DMC) method, using single Slater determinant and multiconfigurational trial functions. DMC calculations using single-determinant trial wave functions of orbitals obtained from density functional theory calculations overestimate the binding energy, while DMC calculations using Hartree-Fock or CAS(4,8), complete active space trial functions significantly underestimate the binding energy. In order to obtain an accurate value of the binding energy of Be{sub 2} from DMC calculations, it is necessary to employ trial functions that include excitations outside the valence space. Our best estimate DMC result for the binding energy of Be{sub 2}, obtained by using configuration interaction trial functions and extrapolating in the threshold for the configurations retained in the trial function, is 908 cm{sup −1}, only slightly below the 935 cm{sup −1} value derived from experiment.

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

A Noncentrosymmetric Polymorph of LuRuGe

In this work, we report a new polymorph of LuRuGe, obtained in indium flux. This phase exhibits the noncentrosymmetric ZrNiAl-type structure with the space group $Ρ\bar 62m$ as determined by single-crystal X-ray diffraction. This polymorph can convert into another centrosymmetric polymorph (TiNiSi-type structure, space group Pnma ) at high temperatures. In this work, we performed electrical transport, magnetization, and specific heat measurements on this new phase. It shows metallic behavior with a Hall sign change from negative at 2 K to positive at 125 K. LuRuGe exhibits Pauli paramagnetism as the ground state with no local magnetic moments from either the Ru or Lu site. The Debye temperature Θ = 348 K and electronic coefficient γ e = 3.6 mJ K –2 mol –1 are extracted from the low-temperature specific heat data in LuRuGe. We also carried out first-principles density functional theory calculations to map out the electronic band structure and density of states. There are several electronic bands crossing the Fermi level, supporting a multiband scenario consistent with the Hall sign change. The density of states around the Fermi level is mainly from Ru 4d and Ge 4p electrons, indicating a strong hybridization between those atomic orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗