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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Block Mott insulating state induced by next-nearest-neighbor hopping in the S = 3 2 zigzag chain BaCoTe 2 O 7

Quasi-one-dimensional correlated electronic multiorbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-field splitting, the electronic correlations (Hubbard repulsion U and Hund coupling J H ), and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe 2 ⁢ O 7 , with electronic density n=7, was prepared experimentally. Here, we systematically study the electronic and magnetic properties of this quasi-one-dimensional compound from the theoretical perspective. Based on first-principles density functional theory calculations, strongly anisotropic one-dimensional electronic Co 3⁢d bands were found near the Fermi level. By evaluating the relevant hopping amplitudes, we provide the magnitude and origin of the nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping matrices in BaCoTe 2 ⁢ O 7 . With this information, we constructed a three-orbital electronic Hubbard model for this zigzag chain system, and studied two cases: with only a NN hopping matrix, and with NN plus NNN hopping matrices. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group method, we constructed the ground-state phase diagram. A robust staggered ↑ - ↓ - ↑ - ↓ antiferromagnetic (AFM) region was found when only the NN hopping matrix in the chain direction was employed. However, for the realistic case where the NNN hopping matrix is also included, the dominant state becomes instead a block AFM ↑ - ↑ - ↓ - ↓ order, in agreement with experiments. The system displays Mott insulator characteristics with three half-filled orbitals, when the block AFM order is stable. In conclusion, our results for BaCoTe 2 ⁢O 7 provide guidance to experimentalists and theorists working on this zigzag one-dimensional chain and related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Cr is not an acceptor in 𝛽⁢−G⁢a 2 ⁢O 3

The intensity of red C⁢r 3+ photoluminescence (PL) in monoclinic gallium oxide (𝛽⁢−G⁢a 2 ⁢O 3 ) is suppressed by 𝑛-type conductivity, an effect that has been attributed to a Cr deep acceptor level in the bandgap. In 𝑛-type material, such an acceptor level would be occupied, resulting in the C⁢r 2+ oxidation state and the absence of C⁢r 3+ PL. To test this model, 𝑛-type 𝛽⁢−G⁢a 2 ⁢O 3 crystals co-doped with Cr and Zr (a donor) were grown from the melt. The samples show C⁢r 3+ optical absorption bands and a high free-electron concentration of 4 × 10 18 c⁢m −3 . If Cr were an acceptor, then it would be fully compensated and therefore would not exhibit the C⁢r 3+ optical signature. Hybrid functional calculations indicate that Cr occupies the substitutional octahedral Ga(II) site and that the C⁢r 2+ state is energetically unfavorable, i.e., Cr is not an acceptor. Weak C⁢r 3+ PL was observed in the Cr/Zr co-doped samples. In conclusion, the quenching of PL may be caused by a transfer of energy to free electrons, a nonradiative process that would reduce the emission intensity.

carrier generation & recombination↗
DOE OSTI2021

Automated Adsorption Workflow for Semiconductor Surfaces and the Application to Zinc Telluride

Surface adsorption is a crucial step in numerous processes, including heterogeneous catalysis, where the adsorption of key species is often used as a descriptor of efficiency. We present here an automated adsorption workflow for semiconductors which employs density functional theory calculations to generate adsorption data in a high-throughput manner. Starting from a bulk structure, the workflow performs an exhaustive surface search, followed by an adsorption structure construction step, which generates a minimal energy landscape to determine the optimal adsorbate-surface distance. An extensive set of energy-based, charge-based, geometric, and electronic descriptors tailored toward catalysis research are computed and saved to a personal user database. Finally, the application of the workflow to zinc telluride, a promising CO 2 reduction photocatalyst, is presented as a case study to illustrate the capabilities of this method and its potential as a material discovery tool.

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

Linking the HOMO-LUMO gap to torsional disorder in P3HT/PCBM blends

The electronic structure of [6,6]-phenyl C{sub 61} butyric acid methyl ester (PCBM), poly(3-hexylthiophene) (P3HT), and P3HT/PCBM blends is studied using soft X-ray emission and absorption spectroscopy and density functional theory calculations. We find that annealing reduces the HOMO-LUMO gap of P3HT and P3HT/PCBM blends, whereas annealing has little effect on the HOMO-LUMO gap of PCBM. We propose a model connecting torsional disorder in a P3HT polymer to the HOMO-LUMO gap, which suggests that annealing helps to decrease the torsional disorder in the P3HT polymers. Our model is used to predict the characteristic length scales of the flat P3TH polymer segments in P3HT and P3HT/PCBM blends before and after annealing. Our approach may prove useful in characterizing organic photovoltaic devices in situ or even in operando.

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

Anion-Activated Bases and Nucleophiles Characterized by Photoelectron Spectroscopy

In this study, negative ion photoelectron spectra at 20 K along with ab initio (CCSD(T)) and M06-2X density function theory calculations are reported for a series of six basic and nucleophilic pyridine derivatives with an anionic substituent (i.e., 3- and 4-PyrBX 3 –, where X = F, 4-t-BuC 6 H 4 , 4- MeOC 6 H 4 , and 3,5-(MeO) 2 C 6 H 3 ). Vertical detachment energies spanning from 4.50 – 5.85 eV are observed and well reproduced by M06-2X/aug-cc-pVTZ and CCSD(T)/maug-cc-pVTZ computations. Surprisingly, the VDEs are found to correlate with the S N 2 reactivity of the PPh 4 + salts of the anions with 1-iodooctane. Spin densities of the neutral structures formed from the anionic species are also reported and indicate where the excess electron density in the anions resides and likely cation-interaction sites in their corresponding salts.

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

The Rise of Neural Networks for Materials and Chemical Dynamics

Machine learning (ML) is quickly becoming a premier tool for modeling chemical processes and materials. ML-based force fields, trained on large data sets of high-quality electron structure calculations, are particularly attractive due their unique combination of computational efficiency and physical accuracy. This Perspective summarizes some recent advances in the development of neural network-based interatomic potentials. Designing high-quality training data sets is crucial to overall model accuracy. One strategy is active learning, in which new data are automatically collected for atomic configurations that produce large ML uncertainties. Another strategy is to use the highest levels of quantum theory possible. Transfer learning allows training to a data set of mixed fidelity. A model initially trained to a large data set of density functional theory calculations can be significantly improved by retraining to a relatively small data set of expensive coupled cluster theory calculations. These advances are exemplified by applications to molecules and materials.

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

Electronic transport, thermal transport, thermal expansion, and magnetization in the strongly correlated metal LaNi⁢O 3

Perovskite structured LaNiO 3 is a strongly correlated metal with intriguing thermal and magnetic properties. The volume dependence of calculated and measured physical properties can add additional critical information to develop a more in-depth understanding of this strongly correlated phenomenon. Taking advantage of recent single crystal LaNiO 3 growth using the floating-zone method, we have measured the thermal expansion, the magnetostriction, and the pressure dependence of the magnetic susceptibility, which then allows derivation of the Grüneisen parameters γ e = $\frac{dlnN(E_F)}{d lnV}$, γ χ = $\frac{dlnχ}{d lnV}$, as well as of electric and thermal transport properties. We simulate the volume dependence of structural and magnetic properties using Density Functional Theory calculations at the Generalized Gradient Approximation level. A large discrepancy between experimental values and calculated ones suggests that strong correlations are likely to be dynamic in nature. This study also provides a side-by-side comparison of measurements in single crystal and polycrystalline samples of LaNiO 3 to elucidate intrinsic materials properties. A broad hump at high temperatures in the temperature dependence of magnetization found in the single crystal sample of LaNiO 3 has been rationalized by a model that includes the influence of electron correlations on the Landau diamagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Vapor liquid solid-hydride vapor phase epitaxy (VLS-HVPE) growth of ultra-long defect-free GaAs nanowires: Ab initio simulations supporting center nucleation

High aspect ratio, rod-like and single crystal phase GaAs nanowires (NWs) were grown by gold catalyst-assisted hydride vapor phase epitaxy (HVPE). High resolution transmission electron microscopy and micro-Raman spectroscopy revealed polytypism-free zinc blende (ZB) NWs over lengths of several tens of micrometers for a mean diameter of 50 nm. Micro-photoluminescence studies of individual NWs showed linewidths smaller than those reported elsewhere which is consistent with the crystalline quality of the NWs. HVPE makes use of chloride growth precursors GaCl of which high decomposition frequency after adsorption onto the liquid droplet catalysts, favors a direct and rapid introduction of the Ga atoms from the vapor phase into the droplets. High influxes of Ga and As species then yield high axial growth rate of more than 100 μm/h. The diffusion of the Ga atoms in the liquid droplet towards the interface between the liquid and the solid nanowire was investigated by using density functional theory calculations. The diffusion coefficient of Ga atoms was estimated to be 3 × 10{sup −9} m{sup 2}/s. The fast diffusion of Ga in the droplet favors nucleation at the liquid-solid line interface at the center of the NW. This is further evidence, provided by an alternative epitaxial method with respect to metal-organic vapor phase epitaxy and molecular beam epitaxy, of the current assumption which states that this type of nucleation should always lead to the formation of the ZB cubic phase.

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

CO 2 Activation on Ni(111) and Ni(100) Surfaces in the Presence of H 2 O: An Ambient-Pressure X-ray Photoelectron Spectroscopy Study

Nickel-based catalysts play an important role in the chemical transformation of CO 2 . A fundamental understanding of the interaction between CO 2 and Ni surfaces at atomic level is necessary. In this study, the interfacial reactions of CO 2 and CO 2 + H 2 O mixture on Ni(111) as well Ni(100) surfaces were investigated using ambient pressure X-ray photoelectron spectroscopy as well as theoretical calculations. The results indicate that the distributions of dissociation products are very different on the Ni(111) and Ni(100) surfaces. Carbonate, CO, and graphitic carbon formed on both Ni(111) and Ni(100) surfaces in the presence of 0.2 Torr CO 2 . However, more than 90% adsorption species on the Ni(111) surface is the carbonate, whereas the Ni(100) surface is mainly covered by adsorbed CO* and graphitic carbon. The co-adsorption of H 2 O weakens the influence of structure effect on the interfacial interaction between CO 2 and Ni surfaces. Furthermore, the density functional theory calculations suggest that the activation of CO 2 on Ni(111) and Ni(100) tends to follow different reaction paths, which are consistent with the experimental results.

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

Evidence of Ba-substitution induced spin-canting in the magnetic Weyl semimetal EuCd 2 As 2

Recently EuCd 2 As 2 was predicted to be a magnetic Weyl semimetal with a lone pair of Weyl nodes generated by A-type antiferromagnetism and protected by a rotational symmetry. However, it was soon discovered that the actual magnetic structure broke the rotational symmetry and internal pressure was later suggested as a route to stabilize the desired magnetic state. In this work we test this prediction by synthesizing a series of Eu 1-x BaxCd 2 As 2 single crystals and studying their structural, magnetic, and transport properties via both experimental techniques and first-principles calculations. We find that small concentrations of Ba (~3%–10%) lead to a small out-of-plane canting of the Eu moment. However, for higher concentrations this effect is suppressed and a nearly in-plane model is recovered. Studying the transport properties we find that all compositions show evidence of an anomalous Hall effect dominated by the intrinsic mechanism as well as large negative magnetoresistances in the longitudinal channel. A nonmonotonic evolution of the transport properties is seen across the series which correlates to the proposed canting suggesting canting may enhance the topological effects. Careful density functional theory calculations using an all-electron approach revise prior predictions finding a purely ferromagnetic ground state with in-plane moments for both the EuCd 2 As 2 and Eu 0.5 Ba 0.5 Cd 2 As 2 compounds, corroborating our experimental findings. This work suggests that Ba substitution can tune the magnetic properties in unexpected ways which correlate to changes in measures of topological properties, encouraging future work to locate the ideal Ba concentration for Eu moment canting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Influence of surface relaxation on solute atoms positioning within atom probe tomography reconstructions

Highlights: • APT measurement of Al/Cu alloy is simulated with a combined electrostatics and molecular dynamics (MD) approach. • Results show segregation of dilute Cu at [111] poles in correspondence to experimental measurements. • No temperature is involved. The effect originates in both disparate evaporation field and bond-energies in the alloy. While thermal or field-induced surface migration has been invoked often in past literature to rationalize the experimental observations of solute and impurity atoms segregation to specific poles and zones lines in atom probe tomography (APT) reconstructions, it does not satisfactorily explain all experimental results. Using a combined finite-element and molecular dynamics simulation approach that models the field evaporation process from mechanically relaxed emitter structures, we show that an additional mechanism, athermal relaxation, can also drive significant solute migration. The mechanism is illustrated for Cu in solid solution in Al, which exhibits segregation of Cu to 111 poles and along [110] zone lines in reconstructed APT data. The higher bonding energy of Cu atoms on the surface, confirmed using density-functional theory calculations, causes longer retention on the surface and preferential evaporation of neighboring Al atoms with subsequent athermal relaxation of Cu into higher-neighbor-count sites left vacant by evaporating Al neighbors. This in turn results in barrier-free long-range migration of the Cu atoms by a series of relaxation events that follow the receding terrace ledges towards high-neighbor count poles and zone lines. This mechanism fundamentally limits spatial resolution in the limit of zero temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Electron interactions and Dirac fermions in graphene-Ge{sub 2}Sb{sub 2}Te{sub 5} superlattices

Graphene based superlattices have been attracted worldwide interest due to the combined properties of the graphene Dirac cone feature and all kinds of advanced functional materials. In this work, we proposed a novel series of graphene-Ge{sub 2}Sb{sub 2}Te{sub 5} superlattices based on the density functional theory calculations. We demonstrated the stability in terms of energy and lattice dynamics for such kind of artificial materials. The analysis of the electronic structures unravels the gap opening nature at Dirac cone of the insert graphene layer. The Dirac fermions in the graphene layers are strongly affected by the electron spin orbital coupling in the Ge{sub 2}Sb{sub 2}Te{sub 5} layers. The present results show the possible application in phase-change data storage of such kind of superlattice materials, where the Ge{sub 2}Sb{sub 2}Te{sub 5} layers exhibit as the phase-change data storage media and the graphene layer works as the electrode, probe, and heat conductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Multistability of isolated and hydrogenated Ga–O divacancies in β–Ga 2 O 3

This work systematically explores 19 unique configurations of the close-associate Ga–O divacancies (V Ga V O ) in β–Ga 2 O 3 , including their complexes with H impurities, using hybrid functional calculations. Interestingly, most configurations are found to retain the negative-U behavior of V O , as they exhibit a thermodynamic (–/3–) charge-state transition level energetically located in the upper part of the band gap, where the 3– charge state is associated with the formation of a Ga–Ga dimer. The energy positions of the thermodynamic (–/3–) charge-state transition levels divide the divacancy configurations into three different groups, which can be understood from the three possible Ga–Ga dimerizations resulting from the tetrahedral and octahedral Ga sites. The relative formation energies of the different divacancy configurations, and hence the electrical activity of the divacancies, is found to depend on the Fermi-level position, and the energy barriers for transformation between different divacancy configurations are explored from nudged elastic band calculations. Hydrogenation of the divacancies is found to either passivate their negative-U charge-state transition levels or shift them down in Fermi level position, depending on whether the H resides at V O or forms an O–H bond at V Ga , respectively. Finally, the divacancy is discussed as a potential origin of the so-called E$^{*}_{2}$ center previously observed by deep-level transient spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kohn anomaly and elastic softening in body-centered cubic molybdenum at high pressure

Transition metals in body-centered cubic (bcc) structures under compression can display several novel physical properties because of their complex electronic structures and electron-phonon interactions. Here, we used inelastic x-ray scattering experiments in a diamond-anvil cell up to ∼45 GPa and density-functional theory calculations up to 210 GPa to investigate the phonon dispersions, and electronic and elastic properties of single-crystal molybdenum (Mo). Our results show a pressure-induced Kohn anomaly at 𝑞∼0.5 along the [ξ00] direction in the longitudinal acoustic mode at ∼45 GPa; this anomaly is triggered by the pressure-enhanced Fermi-surface nesting effect. Theoretical calculations show that electron redistributions in the 𝑠-to-𝑑 orbitals of bcc-Mo contribute to the shear modulus anomaly at ∼50 GPa. In contrast, the Young's modulus anomaly in bcc-Mo at ∼210 GPa results from a Lifshitz-type electronic topological transition. In conclusion, our results shed light on the complex electronic behaviors that are associated with macroscopic elastic properties in typical bcc 𝑑-block transition metals under compression.

3-dimensional systems↗
DOE OSTI2019

Influence of Atomic Surface Structure on the Activity of Ag for the Electrochemical Reduction of CO 2 to CO

This work was undertaken to elucidate the facet-dependent activity of Ag for the electrochemical reduction of CO 2 to CO. To this end, CO2 reduction was investigated over Ag thin films with (111), (100), and (110) orientations prepared via epitaxial growth on single-crystal Si wafers with the same crystallographic orientations. Thus, this preparation technique yielded larger area electrodes than can be achieved using single-crystals, which enabled the electrocatalytic activity of the corresponding Ag surfaces to be quantified in the Tafel regime. The Ag(110) thin films exhibited higher CO evolution activity compared to the Ag(111) and Ag(100) thin films, consistent with previous single-crystal studies. Density functional theory calculations suggest that CO 2 reduction to CO is strongly facet-dependent, and that steps are more active than highly coordinated terraces. This is the result of both a higher binding energy of the key intermediate COOH as well as an enhanced double-layer electric field stabilization over undercoordinated surface atoms located at step edge defects. As a consequence, step edge defects likely dominate the CO 2 reduction activity observed over the Ag(111) and Ag(100) thin films. The higher activity observed over the Ag(110) thin film is then related to the larger density of undercoordinated sites compared to the Ag(111) and Ag(100) thin films. Our conclusion that undercoordinated sites dominate the CO 2 reduction activity observed over close-packed surfaces highlights the need to consider the contribution of such defects in studies of single-crystal electrodes.

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

AL4GAP: Active learning workflow for generating DFT-SCAN accurate machine-learning potentials for combinatorial molten salt mixtures

Machine learning interatomic potentials have emerged as a powerful tool for bypassing the spatiotemporal limitations of ab initio simulations, but major challenges remain in their efficient parameterization. We present AL4GAP, an ensemble active learning software workflow for generating multicomposition Gaussian approximation potentials (GAP) for arbitrary molten salt mixtures. The workflow capabilities include: (1) setting up user-defined combinatorial chemical spaces of charge neutral mixtures of arbitrary molten mixtures spanning 11 cations (Li, Na, K, Rb, Cs, Mg, Ca, Sr, Ba and two heavy species, Nd, and Th) and 4 anions (F, Cl, Br, and I), (2) configurational sampling using low-cost empirical parameterizations, (3) active learning for down-selecting configurational samples for single point density functional theory calculations at the level of Strongly Constrained and Appropriately Normed (SCAN) exchange-correlation functional, and (4) Bayesian optimization for hyperparameter tuning of two-body and many-body GAP models. Here, we apply the AL4GAP workflow to showcase high throughput generation of five independent GAP models for multicomposition binary-mixture melts, each of increasing complexity with respect to charge valency and electronic structure, namely: LiCl–KCl, NaCl–CaCl 2 , KCl–NdCl 3 , CaCl 2 –NdCl 3 , and KCl–ThCl 4 . Our results indicate that GAP models can accurately predict structure for diverse molten salt mixture with density functional theory (DFT)-SCAN accuracy, capturing the intermediate range ordering characteristic of the multivalent cationic melts.

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

Elucidation of the Active Site for the Oxygen Evolution Reaction on a Single Pt Atom Supported on Indium Tin Oxide

Single-atom catalysts (SACs) have attracted attention for their high catalytic activity and selectivity, but the nature of their active sites under realistic reaction conditions, involving various ligands, is not well understood. In this study, we use Density Functional Theory calculations and grand canonical basin hopping to theoretically investigate the active site for the oxygen evolution reaction (OER) on single Pt atom supported on Indium Tin Oxide, including the influence of the electrochemical potential. We show that the ligands on the Pt atom change from Pt-OH in the absence of electrochemical potential to PtO(OH) 4 in electrochemical conditions. This change of the chemical state of Pt is associated to a decrease of 0.3V for the OER overpotential. This highlights the importance of accurately identifying the nature of the active site under reaction conditions and the impact of adsorbates on the electrocatalytic activity. This theoretical investigation enhances our understanding of SACs for OER.

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

High-throughput dataset of impurity adsorption on common catalysts in biomass upgrading applications

Abstract An extensive dataset consisting of adsorption energies of pernicious impurities present in biomass upgrading processes on common catalysts and support materials has been generated. This work aims to inform catalyst and process development for the conversion of biomass-derived feedstocks to fuels and chemicals. A high-throughput workflow was developed to execute density functional theory calculations for a diverse set of atomic (Al, B, Ca, Cl, Fe, K, Mg, Mn, N, Na, P, S, Si, Zn) and molecular (COS, H 2 S, HCl, HCN, K 2 O, KCl, NH 3 ) species on 35 unique surfaces for transition-metal (Ag, Au, Co, Cu, Fe, Ir, Ni, Pd, Pt, Re, Rh, Ru) and metal-oxide (Al 2 O 3 , MgO, anatase-TiO 2 , rutile-TiO 2 , ZnO, ZrO 2 ) catalysts and supports. Approximately 3,000 unique adsorption geometries and corresponding adsorption energies were obtained.

09 BIOMASS FUELS↗