Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso–meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4- phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi’s golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

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

Magnetoelastic interactions in SrCu 2 (BO 3 ) 2 studied by Raman scattering experiments and first principles calculations

Dynamic and static crystal lattice properties of SrCu 2 (BO 3 ) 2 are studied by means of Raman scattering, magnetostriction, and thermal expansion measurements in magnetic fields to 45 T. Raman experiments versus temperature reveal that some phonon modes show an unusual behavior: their frequencies soften (modes at 200 and 450 cm –1 ) while others harden (modes at 385 and 478 cm –1 ) when decreasing the temperature below 15 K. Magneto-Raman experiments show that their field dependence correlates with their respective temperature dependencies; e.g., modes that are hardened with increasing temperature also harden with applied magnetic fields and modes that become softer with temperature also soften with applied fields. We use density functional theory to successfully model and compute the energies of these modes, classifying them into two types: pantograph (modes that soften when decreasing the temperature) and nonpantograph. We conclude that the former involves the modification of the intradimer exchange interaction J and the latter the interdimer J'. Lastly, dilatometry is used to correlate field-dependent Raman modes to the closing of the spin gap as well as fractional-magnetization stripe states M = 1/4 M s and M = 1/3 M s , where M s is the saturation magnetization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Descriptors for Cu facets for CO2 reduction reaction activity

Computation screening is crucial for designing efficient electrochemical catalysts for carbon dioxide (CO2R) reduction to valuable hydrocarbons and oxygenates. Herein, leveraging density functional theory calculations of the CO adsorption energy ΔE_CO on seventeen Cu terminations, we discover a strong linear correlation between ΔE_CO and the recently experimentally measured CO2R electrochemical currents (ACS Catal. 2022, 12, 11, 6578–6588). Examining the ab initio thermodynamics of the early critical intermediates CO*, COH*, and CHO*, we find that CO* → CHO* is the thermodynamically preferred step, and notably shows a volcano trend with the experimental currents where the maximum CO2R current corresponds to the moderate CHO* formation energy. Importantly, we show that increasing the step and kink density of the Cu termination not only enhances CO adsorption strength but also modulates the CO* → CHO* pathway, as respectively exemplified in the (941) and (741) facets. We also explain why (741) is exceptional with high CO2R activity as measured experimentally due to its relatively low activity toward the hydrogen evolution reaction compared with the other Cu surfaces. Beyond the general CO adsorption energy that only shows a linear trend with CO2R activity, we show that the reaction CO* → CHO* free energy is a descriptor that displays a volcano relationship with the overall CO2R activity on Cu facets.

machine learning↗
DOE OSTI2023

Insights into the Oxygen Evolution Reaction on Graphene-Based Single-Atom Catalysts from First-Principles-Informed Microkinetic Modeling

Single-atom transition metals embedded in nitrogen-doped graphene have emerged as promising electrocatalysts due to their high activity and low material cost. These materials have been shown to catalyze a variety of electrochemical reactions, but their active sites under reaction conditions remain poorly understood. Using first-principles density functional theory calculations, we develop a pH-dependent microkinetic model to evaluate the relative performance of transition metal catalysts embedded in fourfold N-substituted double carbon vacancies in graphene for the oxygen evolution reaction. We find that reaction pathways involving intermediates co-adsorbed on the metal site are preferred on all transition metals. These pathways lead to enhancements in catalytic activity and broaden the activity peak when compared with purely thermodynamics-based predictions. Furthermore, these findings demonstrate the importance of investigating reaction pathways on graphene-based catalysts and other two-dimensional (2D) materials that involve metal active centers decorated by spectator intermediate species.

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

High‐Field NMR, Reactivity, and DFT Modeling Reveal the γ‐Al 2 O 3 Surface Hydroxyl Network**

Abstract Aluminas are strategic materials used in many major industrial processes, either as catalyst supports or as catalysts in their own right. The transition alumina γ‐Al 2 O 3 is a privileged support, whose reactivity can be tuned by thermal activation. This study provides a qualitative and quantitative assessment of the hydroxyl groups present on the surface of γ‐Al 2 O 3 at three different dehydroxylation temperatures. The principal [AlOH] configurations are identified and described in unprecedented detail at the molecular level. The structures were established by combining information from high‐field 1 H and 27 Al solid‐state NMR, IR spectroscopy and DFT calculations, as well as selective reactivity studies. Finally, the relationship between the hydroxyl structures and the molecular‐level structures of the active sites in catalytic alkane metathesis is discussed.

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

Silica-copper catalyst interfaces enable carbon-carbon coupling towards ethylene electrosynthesis

Membrane electrode assembly (MEA) electrolyzers offer a means to scale up CO 2 -to-ethylene electroconversion using renewable electricity and close the anthropogenic carbon cycle. To date, excessive CO 2 coverage at the catalyst surface with limited active sites in MEA systems interferes with the carbon-carbon coupling reaction, diminishing ethylene production. With the aid of density functional theory calculations and spectroscopic analysis, here we report an oxide modulation strategy in which we introduce silica on Cu to create active Cu-SiO x interface sites, decreasing the formation energies of OCOH* and OCCOH*—key intermediates along the pathway to ethylene formation. We then synthesize the Cu-SiO x catalysts using one-pot coprecipitation and integrate the catalyst in a MEA electrolyzer. By tuning the CO 2 concentration, the Cu-SiO x catalyst based MEA electrolyzer shows high ethylene Faradaic efficiencies of up to 65% at high ethylene current densities of up to 215 mA cm -2 ; and features sustained operation over 50 h.

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