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

Tailoring Bandgap of Perovskite BaTiO 3 by Transition Metals Co-Doping for Visible-Light Photoelectrical Applications: A First-Principles Study

The physical and chemical properties of V-M" and Nb-M" (M" is 3d or 4d transition metal) co-doped BaTiO 3 were studied by first-principles calculation based on density functional theory. Our calculation results show that V-M" co-doping is more favorable than Nb-M" co-doping in terms of narrowing the bandgap and increasing the visible-light absorption. In pure BaTiO 3 , the bandgap depends on the energy levels of the Ti 3d and O 2p states. The appropriate co-doping can effectively manipulate the bandgap by introducing new energy levels interacting with those of the pure BaTiO 3 . The optimal co-doping effect comes from the V-Cr co-doping system, which not only has smaller impurity formation energy, but also significantly reduces the bandgap. Detailed analysis of the density of states, band structure, and charge-density distribution in the doping systems demonstrates the synergistic effect induced by the V and Cr co-doping. The results can provide not only useful insights into the understanding of the bandgap engineering by element doping, but also beneficial guidance to the experimental study of BaTiO 3 for visible-light photoelectrical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Marcus–Hush–Chidsey kinetics at electrode–electrolyte interfaces

Electrochemical kinetics at electrode–electrolyte interfaces limit the performance of devices including fuel cells and batteries. While the importance of moving beyond Butler–Volmer kinetics and incorporating the effect of electronic density of states of the electrode has been recognized, a unified framework that incorporates these aspects directly into electrochemical performance models is still lacking. In this work, we explicitly account for the density functional theory-calculated density of states numerically in calculating electrochemical reaction rates for a variety of electrode–electrolyte interfaces. We first show the utility of this for two cases related to Li metal electrodeposition and stripping on a Li surface and a Cu surface (anode-free configuration). The deviation in reaction rates is minor for cases with flat densities of states such as Li, but is significant for Cu due to nondispersive d-bands creating large variation. Finally, we consider a semiconducting case of a solid-electrolyte interphase consisting of LiF and Li2CO3 and note the importance of the Fermi level at the interface pinned by the redox reaction occurring there. We identify the asymmetry in reaction rates as a function of discharge/charge naturally within this approach.

Kurchin, Rachel (ORCID:0000000221474809)↗
DOE OSTI2014

Ab initio calculations on magnetism induced by composite defects in magnesium oxide

The local magnetic state induced by the composite defects, composed of an oxygen vacancy and a nitrogen substituting oxygen, in magnesium oxide has been studied by using ab initio calculation based on density functional theory. The calculated results show that local magnetic moment can be induced by the composite defects around the oxygen vacancy, when the exchange split of the oxygen vacancy is enhanced either by the hybridization between the N-p and nearest neighbor O-p orbitals or by applying on-site Coulomb repulsion (U) and exchange interaction (J). We show that the magnetic state induced by the composite defect is energetically more stable than the non-magnetic state. In addition, we show that the U and J applied on the p-orbitals of N and O atoms may significantly impact the calculated magnetic state of the composite defect, resulting in magnetic state for a configuration that is non-magnetic by generalized gradient approximation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Point defects and doping in wurtzite LaN

Wurtzite LaN (wz-LaN) is a semiconducting nitride that has piezoelectric and ferroelectric properties, making it promising for applications in electronics, either as a binary compound or in alloys such as LaAlN. The prospects for wz-LaN in devices are influenced by the properties of point defects and impurities; here, we use first-principles density functional theory with a hybrid functional to calculate their formation energies, as well as their atomic and electronic structures. Among native point defects, we find that nitrogen-related defects, both vacancies ($V^+_N$) and interstitials ($N^-_i$), are energetically most favorable under most relevant chemical potentials and positions of the Fermi level; $V^0_N$ may additionally be observed under N-poor conditions, and $N^0_i$ may be prominent under N-rich conditions. We also investigate the incorporation of oxygen and hydrogen, which will likely be present as unintentional impurities. We find that the $O^+_N$ substitutional species readily forms, but oxygen will not lead to n-type conductivity due to formation of DX centers and compensation by interstitial defects. Similarly, substitutional HN and interstitial H i can compensate both p- and n-type dopants. Our results provide detailed, microscopic guidance for the development of electronic devices based on wz-LaN.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Trivial to nontrivial topology transition in rare-earth pnictides with epitaxial strain

The combination of magnetotransport and topological properties has brought great attention to rare-earth monopnictides semimetals. For some of them, like LaSb, it is unclear whether they show nontrivial topology or not based on density functional theory calculations and angular resolved photoemission spectroscopy measurements. Here we use hybrid density functional theory to demonstrate that LaSb is in fact a trivial topological semimetal, in agreement with experiments, but on the verge of a transition to a topological phase. In this work, we show that under compressive epitaxial strain, the La $\textit{d}$ band crosses the Sb $\textit{p}$ band near the $\textit{Z}$ point in the Brillouin zone, stabilizing a topologically nontrivial phase, opening unique opportunities to probe the inter-relation between magnetotransport properties and the effects of band topology by examining epitaxially strained and unstrained thin films of the same material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Accuracy of finite-difference harmonic frequencies in density functional theory

Analytic Hessians are often viewed as essential for the calculation of accurate harmonic frequencies, but the implementation of analytic second derivatives is nontrivial and solution of the requisite coupled-perturbed equations engenders a sizable memory footprint for large systems, given that these equations are not required for energy and gradient calculations in density functional theory. Here, we benchmark the alternative approach to harmonic frequencies based on finite differences of analytic first derivatives, a procedure that is amenable to large-scale parallelization. Not only for absolute frequencies but also for isotopic and conformer-dependent frequency shifts in flexible molecules, we find that the finite-difference approach exhibits mean errors < 0.1 cm -1 as compared to results based on an analytic Hessian. For very small frequencies corresponding to nonbonded vibrations in noncovalent complexes (for which the harmonic approximation is questionable anyway), the finite-difference error can be larger, but even in these cases the errors can be reduced below 0.1 cm -1 by judicious choice of the displacement step size and a higher-order finite-difference approach. Finally, the surprising accuracy and robustness of the finite-difference results suggests that availability of the analytic Hessian is not so important in today's era of commodity processors that are readily available in large numbers.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Simulating the Nanomechanical Response of Cyclooctatetraene Molecules on a Graphene Device

We investigate the atomic and electronic structures of cyclooctatetraene (COT) molecules on graphene and analyze their dependence on external gate voltage using first-principles calculations. The external gate voltage is simulated by adding or removing electrons using density functional theory calculations. This allows us to investigate how changes in carrier density modify the molecular shape, orientation, adsorption site, diffusion barrier, and diffusion path. For increased hole doping, COT molecules gradually change their shape to a more flattened conformation and the distance between the molecules and graphene increases while the diffusion barrier drastically decreases. For increased electron doping, an abrupt transition to a planar conformation at a carrier density of -8 × 10 13 e/cm 2 is observed. Furthermore, these calculations imply that the shape and mobility of adsorbed COT molecules can be controlled by externally gating graphene devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Theoretical Investigations into Defected Graphene for Electrochemical Reduction of CO 2

Here, despite numerous experimental efforts that have been dedicated to studying carbon-based materials for electrochemical reduction of CO 2 , a rationalization of the associated trends in the intrinsic activity of different active motifs has so far been elusive. In the present work, we employ density functional theory calculations to examine a variety of different active sites in N-doped graphene to give a comprehensive outline of the trends in activity. We find that adsorption energies of COOH* and CO* do not follow the linear scaling relationships observed for the pure transition metals, and this unique scaling is rationalized through differences in electronic structure between transition metals and defected graphene. This finding rationalizes most of the experimental observations on the carbon-based materials which present promising catalysts for the two-electron reduction of CO 2 to CO. With this simple thermodynamic analysis, we identify several active sites that are expected to exhibit a comparable or even better activity to the state-of-the-art gold catalyst, and several configurations are suggested to be selective for CO 2 RR over HER.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nature of native atomic defects in ZrTe 5 and their impact on the low-energy electronic structure

Over the past decades, investigations of the anomalous low-energy electronic properties of ZrTe 5 have reached a wide array of conclusions. An open question is the growth method's impact on the stoichiometry of ZrTe 5 samples, especially given the very small density of states near its chemical potential. In this work, we report on high-resolution scanning tunneling microscopy and spectroscopy measurements performed on samples grown via different methods. Using density functional theory calculations, we identify the most prevalent types of atomic defects on the surface of ZrTe 5 , namely, Te vacancies and intercalated Zr atoms. Finally, we precisely quantify their density and outline their role as ionized defects in the anomalous resistivity of this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Kinetics of Silicon field evaporation

Highlights: • A theory of field evaporation developed for metals is applied to Silicon. • Potential energy curves and evaporation field strengths are calculated with density functional theory. • Good agreement with experimental results is obtained for Silicon. The theory of field evaporation, previously developed for metals and adsorbates, is applied to semiconductors and insulators, specifically Si. Potential energy curves and evaporation field strengths of Si{sub 6} and Si{sub 7} are calculated with density functional theory. The evaporation field strength is found to be proportional to the coordination. Good agreement with experimental results is obtained for Silicon as well as an estimate of the experimental field strength.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ti-6Al-4V to over 1.2 TPa: Shock Hugoniot experiments, ab initio calculations, and a broad-range multiphase equation of state

Titanium alloys are used in a large array of applications. In this work we focus our attention on the most used alloy, Ti-6Al-4V (Ti64), which has excellent mechanical and biocompatibility properties with applications in aerospace, defense, biomedical, and other fields. Here we present high-fidelity experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for Ti64 to more than threefold compression, up to over 1.2 TPa. We use the data to validate our ab initio molecular dynamics simulations and to develop a highly reliable, multiphase equation of state (EOS) for Ti64, spanning a broad range of temperature and pressures. The first-principles simulations show very good agreement with Z data and with previous three-stage gas gun data from Sandia’s STAR facility. The resulting principal Hugoniot and the broad-range EOS and phase diagram up to 10 TPa and 10 5 K are suitable for use in shock experiments and in hydrodynamic simulations. The high-precision experimental results and high-fidelity simulations demonstrate that the Hugoniot of the Ti64 alloy is stiffer than that of pure Ti and reveal that Ti64 melts on the Hugoniot at a significantly lower pressure and temperature than previously modeled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Direct Measure of the Dense Methane Phase in Gas Shale Organic Porosity by Neutron Scattering

Here, we report the first direct measurements of methane density in shale gas using small-angle neutron scattering. At a constant pressure, the density of methane in the inorganic pores is similar to the gas bulk density of the system conditions. Conversely, the methane density is 2.1 ± 0.2 times greater in the organic mesopores. Furthermore, classical density functional theory calculations show that this excess density in the organic pores persists to elevated temperatures, typical of shale gas reservoir conditions, providing new insight into the hydrocarbon storage mechanisms within these reservoirs.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

In Situ Studies of Methanol Decomposition Over Cu(111) and Cu2O/Cu(111): Effects of Reactant Pressure, Surface Morphology, and Hot Spots of Active Sites

The dissociative adsorption of methanol was investigated on Cu(111) and ultrathin Cu 2 O films. We employed synchrotron-based Ambient Pressure X-ray Photoelectron Spectroscopy (AP-XPS) and Scanning Tunneling Microscopy (STM) to study the dynamics of gas–solid interactions, and calculations based on Density Functional Theory (DFT) were used to examine the reaction path. C 1s XPS spectra revealed that methanol underwent dissociative adsorption on plain Cu(111) to form methoxy (CH 3 O), formaldehyde (H 2 CO), and formate (HCOO) at a pressure range of 0.5–10 mTorr, with these species remaining on the surface after evacuation. This was accompanied by the appearance of a low coverage (~0.05 ML) of O ads in the O 1s which can be considered a highly active site for methanol activation. The high activity is apparent by a coverage of 0.8 ML of methoxy at room temperature. STM was unable to image these species at room temperature as they were highly mobile on metallic copper. In contrast, for CH 3 OH on Cu 2 O/Cu(111), STM showed clear hot spots for reaction and a complex array of adsorption structures. On the oxide substrate, there was decomposition of methanol to H 2 CO, CH 3 O, HCOO, and hydrocarbon species (CH x ) due to the subsequent interactions of methanol with lattice oxygen. Cu(111) remained entirely saturated with decomposition products under 10 mTorr of methanol (θ ≈ 0.97 ML), whereas the Cu 2 O overlayer was saturated at a much lower coverage (θ ≈ 0.30 ML). STM revealed rows and step edges of Cu 2 O decorated with decomposition products and metallic Cu islands ~5 nm in size. The difference in activity between Cu(111) and Cu 2 O/Cu(111) is attributed to the significant amount of O present on the oxide surface. Additionally, Density Functional theory (DFT) calculations described the XPS measurements well, showing a likely methanol dissociation to *CH 3 O and therefore a surface reduction. More importantly, the DFT results revealed that it was the chemisorbed oxygen on Cu 2 O/Cu(111) which oxidized the dissociated *CH 3 O to *HCOO and eventually CO 2 , while the reaction only involving upper oxygen on the Cu 2 O hexagonal ring led to the formation of H 2 CO.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pressure-induced charge orders and their postulated coupling to magnetism in hexagonal multiferroic LuFe 2 O 4

Hexagonal LuFe 2 O 4 is a promising charge order (CO) driven multiferroic material with high charge and spin-ordering temperatures. The coexisting charge and spin orders on Fe 3+ /Fe 2+ sites result in magnetoelectric behaviors, but the coupling mechanism between the charge and spin orders remains elusive. Here, by tuning external pressure, we reveal three charge-ordered phases with suggested correlation to magnetic orders in LuFe 2 O 4 : (i) a centrosymmetric incommensurate three-dimensional CO with ferrimagnetism, (ii) a non-centrosymmetric incommensurate quasi-two-dimensional CO with ferrimagnetism, and (iii) a centrosymmetric commensurate CO with antiferromagnetism. Experimental in situ single-crystal X-ray diffraction and X-ray magnetic circular dichroism measurements combined with density functional theory calculations suggest that the charge density redistribution caused by pressure-induced compression in the frustrated double-layer [Fe 2 O 4 ] cluster is responsible for the correlated spin-charge phase transitions. The pressure-enhanced effective Coulomb interactions among Fe-Fe bonds drive the frustrated (1/3, 1/3) CO to a less frustrated (1/4, 1/4) CO, which induces the ferrimagnetic to antiferromagnetic transition. Our results not only elucidate the coupling mechanism among charge, spin, and lattice degrees of freedom in LuFe 2 O 4 , but also provide a new way to tune the spin-charge orders in a highly controlled manner.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Collisional stopping power of ions in warm dense matter

Here, a model for the collisional stopping of ions on free electrons in warm dense matter is developed and explored. It is based on plasma kinetic theory, but with modifications to address the warm dense matter regime. Specifically, it uses the Boltzmann-Uehling-Uhlenbeck kinetic equation to incorporate effects of Fermi degeneracy of electrons. The cross section is computed from quantum scattering of electrons and ions occurring via the potential of mean force derived from an average atom model, which incorporates effects of strong Coulomb correlations. Predictions from this model show comparable accuracy to results from time-dependent density functional theory calculations for deuterium near solid density and a temperature of several electronvolts, at a fraction of the computational cost. Further, the model captures the transition of a plasma from the classical limit to the degenerate limit, including qualitative behaviors of solid state theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Perfect separation of intraband and interband excitations in PdCoO 2

The temperature dependence of the optical properties of the delafossite PdCoO 2 has been measured in the a-b planes over a wide frequency range. The optical conductivity due to the free-carrier (intraband) response falls well below the interband transitions, allowing the plasma frequency to be determined from the f-sum rule. Drude-Lorentz fits to the complex optical conductivity yield estimates for the free-carrier plasma frequency and scattering rate. Here, the in-plane plasma frequency has also been calculated using density functional theory. The experimentally determined and calculated values for the plasma frequencies are all in good agreement; however, at low temperature the optically determined scattering rate is much larger than the estimate for the transport scattering rate, indicating a strong frequency-dependent renormalization of the optical scattering rate. In addition to the expected in-plane infrared-active modes, two very strong features are observed that are attributed to the coupling of the in-plane carriers to the out-of-plane longitudinal optic modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hydrogen in disordered titania: connecting local chemistry, structure, and stoichiometry through accelerated exploration

Hydrogen incorporation in native surface oxides of metal alloys often controls the onset of metal hydriding, with implications for materials corrosion and hydrogen storage. A key representative example is titania, which forms as a passivating layer on a variety of titanium alloys for structural and functional applications. These oxides tend to be structurally diverse, featuring polymorphic phases, grain boundaries, and amorphous regions that generate a disparate set of unique local environments for hydrogen. Here, we introduce a workflow that can efficiently and accurately navigate this complexity. First, a machine learning force field, trained on ab initio molecular dynamics simulations, was used to generate amorphous configurations. Density functional theory calculations were then performed on these structures to identify local oxygen environments, which were compared against experimental observations. Second, to classify subtle differences across the disordered configuration space, we employ a graph-based sampling procedure. Finally, local hydrogen binding energies and hopping kinetics are computed using exhaustive density functional theory calculations on representative configurations. Here, we leverage this methodology to show that hydrogen binding energetics are described by local oxygen coordination, which in turn is affected by stoichiometry, and form the basis of hopping kinetics and diffusion. Together these results imply that hydrogen incorporation and transport in TiO x can be tailored through compositional engineering, with implications for improving the performance and durability of titanium-derived alloys in hydrogen environments.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Geometry and electronic structure of iridium adsorbed on graphene

We report investigation of the geometry and electronic structure of iridium atoms adsorbed onto graphene through a combined experimental and theoretical study. Ir atoms were deposited onto a flake of graphene on a Pt(111) surface and found to form clusters even at low temperatures. The areal density of the observed clusters on the graphene flake suggests the clusters are most likely pairs of Ir atoms. Theoretical ab initio density functional (DFT) calculations indicate that these Ir dimers are oriented horizontally, near neighboring "bridge" sites of the graphene lattice, as this configuration has the strongest adsorption energy of all high-symmetry configurations for the Ir dimer. A large peak in the local density of states (LDOS) at the Dirac point energy was measured via scanning tunneling spectroscopy, and this result is reproduced by a DFT calculation of the LDOS. The peak at the Dirac point energy is found to be from the Ir s and p states. Finally, the LDOS in the monomer case was also calculated, and is found to significantly differ from the experimentally determined data, further supporting the hypothesis of low-temperature clustering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗