Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT 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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

In-situ Synthesis of PuCl 3 and Corresponding Raman and Density Functional Theory Vibrational Modes

The experimental and calculated Raman spectrum of PuCl 3 has been reported for the first time. PuCl 3 is a primary species found in plutonium metal refinement, specifically in pyrochemical salt processes including multicycle direct oxide reduction, metal chlorination, and electrorefining. As such, Raman signatures of PuCl 3 could serve as potential forensic indicators of material process history. A novel technique for synthesizing PuCl 3 from the in situ chlorination of plutonium metal with HCl was developed to establish these signatures. Cerium metal surrogates were utilized to ensure optimization of the plutonium experiments and to minimize personnel exposure, and all experiments were carried out in a Raman reaction chamber designed for air-tight, high vacuum environments. In situ Raman spectroscopy was employed in conjunction with density functional theory (DFT) to investigate the vibrational modes of PuCl 3 . Associated mixed oxy and hydroxyl phases are also reported. The combined Raman and DFT results have eliminated inconsistencies in Raman mode assignments for the MX 3 family of metal chlorides having P6 3 /m symmetry, and IR modes derived from the DFT calculations are additionally presented. The data observed in this study are of potential interest to nuclear forensic analyses, nuclear sample aging, nuclear energy, plutonium processing, and stockpile stewardship.

36 MATERIALS SCIENCE↗
DOE OSTI2023

In-situ Synthesis of PuCl 3 and Corresponding Raman and Density Functional Theory Vibrational Modes

The experimental and calculated Raman spectrum of PuCl 3 has been reported for the first time. PuCl 3 is a primary species found in plutonium metal refinement, specifically in pyrochemical salt processes including multicycle direct oxide reduction, metal chlorination, and electrorefining. As such, Raman signatures of PuCl 3 could serve as potential forensic indicators of material process history. A novel technique for synthesizing PuCl 3 from the in-situ chlorination of plutonium metal with HCl was developed to establish these signatures. Cerium metal surrogates were utilized to ensure optimization of the plutonium experiments and to minimize personnel exposure, and all experiments were carried out in a Raman reaction chamber designed for air-tight, high vacuum environments. In-situ Raman spectroscopy was employed in conjunction with Density Functional Theory (DFT) to investigate the vibrational modes of PuCl 3 . Associated mixed oxy and hydroxyl phases are also reported. The combined Raman and DFT results have eliminated inconsistencies in Raman mode assignments for the MX 3 family of metal chlorides having P6 3 /m symmetry, and IR modes derived from the DFT calculations are additionally presented. The data observed in this study is of potential interest to nuclear forensic analyses, nuclear sample aging, nuclear energy, plutonium processing, and stockpile stewardship.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DFT-guided flux synthesis of a family of layered titanates crystallizing in the lepidocrocite structure type

K 0.72 Fe 0.72 Ti 1.28 O 4 , a layered mixed-metal oxide belonging to the lepidocrocite structure type, was grown out of a molten KCl-KF mixture at 800 °C. The compound crystallizes in orthorhombic space group Cmcm and exhibits two-dimensional layers, composed of mixed Fe/Ti octahedra, separated by K ions that occupy the interlayer layer space and maintain charge balance. The targeted flux synthesis of the Ni, Cu, and Zn analogs was carried out after density functional theory (DFT) calculations predicted their formation to be energetically favorable, yielding K 0.77 Ni 0.38 Ti 1.62 O 4 , K 0.77 Cu 0.38 Ti 1.62 O 4 , and K 0.73 Zn 0.36 Ti 1.64 O 4 . Finally, all four compositions can be prepared via conventional solid-state technique by mixing and heating stoichiometric amounts of KNO 3 , appropriate transition metal oxide precursors, and TiO 2 .

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

Catalytic N–H Bond Formation Promoted by a Ruthenium Hydride Complex Bearing a Redox-Active Pyrimidine-Imine Ligand

The synthesis of a piano-stool ruthenium hydride, [(η 5 -C 5 Me 5 )Ru(PmIm)H] (PmIm = (N-(1,3,5-trimethylphenyl)-1-(pyrimidin-2-yl)ethan-1-imine) for the dual purpose of catalytic dihydrogen activation and subsequent hydrogen atom transfer for the formation of weak chemical bonds is described. Introduction of a neutral, potentially redox-active PmIm supporting ligand was designed to eliminate the possibility of deleterious C(sp 2 )–H reductive coupling and elimination that has been identified as a deactivation pathway with related rhodium and iridium catalysts. Treatment of [(η 5 -C 5 Me 5 )RuCl 2 ] n with one equivalent PmIm ligand in the presence of zinc and sodium methoxide resulted in isolation of the diruthenium complex, [(η 5 -C 5 Me 5 )Ru(PmIm)] 2 , arising from C–C bond formation between two PmIm chelates. Addition of H 2 to the ruthenium dimer under both thermal and blue light irradiation conditions furnished the targeted hydride, [(η 5 -C 5 Me 5 )Ru(PmIm)H] which has a relatively weak DFT-calculated Ru–H bond dissociation free energy (BDFE) of 47.9 kcal/mol. Addition of TEMPO to [(η 5 -C 5 Me 5 )Ru(PmIm)H] generated the 17-electron metalloradical, [(η 5 -C 5 Me 5 )Ru(PmIm)], which was characterized by EPR spectroscopy. The C–C bond forming process was reversible as irradiation of [(η 5 -C 5 Me 5 )Ru(PmIm)] 2 generated [(η 5 -C 5 Me 5 )Ru(PmIm)H] and a piano-stool ruthenium complex containing an enamide ligand derived from H-atom ab-straction from the PmIm chelate. Equilibration studies were used to establish an experimental estimate of the effective Ru–H BDFE and a value of 50.8 kcal/mol was obtained, in agreement with the observed loss of H 2 and the DFT-computed value. Here, the ruthenium hydride was an effective catalyst for the thermal catalytic hydrogenation of TEMPO, acridine and a cobalt-imido complex and for the selective reduction of azobenzene to diphenylhydrazine using H 2 as the terminal reductant, highlighting the role of this complex is catalytic weak bond formation using H 2 as the stoichiometric reductant.

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

Degradation and fate of 2,4-dinitroanisole (DNAN) and its intermediates treated with Mg/Cu bimetal: Surface examination with XAS, DFT, and LDI-MS

In this article, a novel Mg-based bimetal reagent (Mg/Cu) was used as an enhanced reductive system to degrade insensitive munition 2,4-dinitroanisole (DNAN), a contaminant found in energetic-laden waste. Degradation of DNAN was significantly impacted by dissolved oxygen and studied in anoxic and oxic bimetal systems (i.e., purging with N 2 , air, or O 2 gas). Degradation occurred through sequential nitroreduction: first one nitro group was reduced (ortho or para) to form short-lived intermediates 2-amino-4-nitroanisole or 4-amino-2-nitroanisole (2-ANAN or 4-ANAN), and then subsequent reduction of the other nitro group formed 2,4-diaminoanisole (DAAN). The nitro-amino intermediates demonstrated regioselective reduction in the ortho position to 2-ANAN; Regioselectivity was also impacted by the anoxic/oxic environment. Under O 2 -purging DNAN degradation rate was slightly enhanced, but most notably O 2 significantly accelerated DAAN generation. DAAN also further degraded only in the oxygenated Mg/Cu system. Adsorption of DNAN byproducts to the reagent occurred regardless of anoxic/oxic condition, resulting in a partition of carbon mass between the adsorbed phase (27%-35%) and dissolved phase (59%-72%). Additional surface techniques were applied to investigate contaminant interaction with Cu. Density functional theory (DFT) calculations identified preferential adsorption structures for DNAN on Cu with binding through two O atoms of one or both nitro groups. X-ray absorption spectroscopy (XAS) measurements determined the oxidation state of catalytic metal Cu and formation of a Cu-O-N bond during treatment. Laser desorption ionization mass spectrometry (LDI-MS) measurements also identified intermediate 2-ANAN adsorbed to the bimetal surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Accelerating defect predictions in semiconductors using graph neural networks

First-principles computations reliably predict the energetics of point defects in semiconductors but are constrained by the expense of using large supercells and advanced levels of theory. Machine learning models trained on computational data, especially ones that sufficiently encode defect coordination environments, can be used to accelerate defect predictions. Here, we develop a framework for the prediction and screening of native defects and functional impurities in a chemical space of group IV, III–V, and II–VI zinc blende semiconductors, powered by crystal Graph-based Neural Networks (GNNs) trained on high-throughput density functional theory (DFT) data. Using an innovative approach of sampling partially optimized defect configurations from DFT calculations, we generate one of the largest computational defect datasets to date, containing many types of vacancies, self-interstitials, anti-site substitutions, impurity interstitials and substitutions, as well as some defect complexes. We applied three types of established GNN techniques, namely crystal graph convolutional neural network, materials graph network, and Atomistic Line Graph Neural Network (ALIGNN), to rigorously train models for predicting defect formation energy (DFE) in multiple charge states and chemical potential conditions. We find that ALIGNN yields the best DFE predictions with root mean square errors around 0.3 eV, which represents a prediction accuracy of 98% given the range of values within the dataset, improving significantly on the state-of-the-art. We further show that GNN-based defective structure optimization can take us close to DFT-optimized geometries at a fraction of the cost of full DFT. The current models are based on the semi-local generalized gradient approximation-Perdew–Burke–Ernzerhof (PBE) functional but are highly promising because of the correlation of computed energetics and defect levels with higher levels of theory and experimental data, the accuracy and necessity of discovering novel metastable and low energy defect structures at the PBE level of theory before advanced methods could be applied, and the ability to train multi-fidelity models in the future with new data from non-local functionals. The DFT-GNN models enable prediction and screening across thousands of hypothetical defects based on both unoptimized and partially optimized defective structures, helping identify electronically active defects in technologically important semiconductors.

Rahman, Md Habibur (ORCID:000000027705984X)↗
DOE OSTI2025

Formation of U(VI) peroxide nanoclusters from cascade reactions with a persulfate radical initiator

Radiolysis of water in high radiation fields generates a variety of reactive oxygen species that influence the chemical behavior and complexation of hexavalent uranium. This study investigates the behavior of interaction of a uranyl cation (UO 2 2+ (VI)) with a series of free radicals that are formed in situ via activation of the free radical initiator persulphate (S 2 O 8 2− ), which releases both SO 4 ˙ − and ˙OH species in the solution. Electron Paramagnetic Resonance (EPR) and Raman spectroscopy were used to evaluate the presence of the hydroperoxyl radical (HO 2 ˙) and superoxide radicals (O 2 ˙ − ) that are formed within the solution through radical cascade reactions. In addition, a uranyl peroxide cluster solid (NaU 24 ) was crystallized and characterized using single crystal X-ray diffraction (SCXRD), vibrational spectroscopy, and EPR spectroscopy. The presence of the hydroperoxyl radical (HO 2 ˙) and superoxide radicals (O 2 ˙ − ) was also observed in the solid-state compound, but spectroscopic evidence suggests that it was associated with the Na + network and not the cluster itself. Density functional theory (DFT) calculations were also utilized to further confirm the radical species produced and determine the potential stabilization of radicals detected within the cluster and lattice.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Elucidating facet dependent electronic and electrochemical properties of Cu 2 O nanocrystals using AFM/SCEM and DFT

Cuprous oxide (Cu 2 O) has extensively been studied owing to its excellent optical, magnetic, and catalytic properties. Many of these properties are facet-dependent and have not been well elucidated. This work synthesized cubic, cuboctahedral, octahedral, and rhombic dodecahedral shaped Cu 2 O nanocrystals of ~300 nm in size to evaluate the facet-dependent electrochemical activities. Here, cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) were firstly used to reveal the average electrochemical activities at the ensemble level. Atomic force microscopy-scanning electrochemical microscopy (AFM-SECM) was further used to assess the electrochemical activities of different Cu 2 O nanocrystals at the facet level. Hexaammineruthenium (III) chloride ({Ru(NH 3 ) 6 }Cl 3 ) was employed as the probe molecules that reacted with four different Cu 2 O nanocrystals under —400mV and yielded ~300 pA current between the probing tip and the nanocrystal surface. The tip-current mapping results indicate that rhombic dodecahedral Cu 2 O exhibits higher electrocatalytic activity than other shaped Cu 2 O, due to the presence of dominant exposed facet of {110} as indicated by the relatively high tip current. Density-functional theory (DFT) calculations confirmed the facet dependence of local surface energy and electronic structure of Cu 2 O nanocrystals. Besides electrochemical activity, the surface work function and adsorptive properties were both observed to vary with the shape and dominant exposed facets of Cu 2 O. This study presented a unique experimental and computational chemistry approach to analyze surface electrochemical properties of Cu 2 O crystals at a crystalline facet level.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Electronic Structure of Zirconium and Hafnium Monochalcogenides

High-level ab initio CCSD(T) and spin–orbit icMRCI+Q calculations were used to predict potential energy curves (PECs) for the lowest-lying states of ZrO, ZrS, HfO, and HfS. The prediction of the ground state is basis set dependent at the icMRCI+Q level for ZrO and ZrS due to the small singlet–triplet splitting between the lowest 1 Σ + and 3 Δ states. CCSD(T) with a spin orbit correction predicted the 1 Σ + ground state in agreement with experiment. New all-electron basis sets were developed for Hf to improve the results over those predicted by use of effective core potentials (ECPs) that subsume the 4f electrons into the definition of the core. The use of the new DK-4f basis sets rather than ECPs became more important for HfO and HfS where there is a lack of a good core–valence separation. icMRCI+Q, CCSD(T), and DFT calculations for the spectroscopic parameters of ZrO, ZrS, HfO, and HfS were benchmarked with available experimental data. Bond dissociation energies (BDEs) of these four systems were calculated at the Feller–Peterson–Dixon (FPD) level to be 762.1 (ZrO), 543.5 (ZrS), 803.8 (HfO), and 575.1 kJ/mol (HfS), in excellent agreement with experiment. The HfS BDE was remeasured using the R3PI method, providing an updated experimental measurement of D 0 (HfS) = 5.978 ± 0.002 eV = 576.8 ± 0.2 kJ/mol. This experimental value, combined with experimental measurements of the ionization energies of Hf and HfS, gives the cationic BDE of D 0 (Hf + -S) = 5.124 ± 0.002 eV = 494.4 ± 0.2 kJ/mol.

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

Gas-phase negative ion photoelectron spectroscopy and reactivity of phenylacetylide

Arylacetylides such as phenylacetylide (PhCC–, 1) are important nucleophiles used in synthetic chemistry yet rarely have they been studied as bare carbanions. In this work, the phenylacetylide anion was formed via electrospray ionization (ESI) or multistage mass spectrometry (MSn) experiments and subsequently examined in the gas phase by negative ion photoelectron spectroscopy (NIPES), ion-molecule reactions (IMR), alongside theoretical calculations to probe its fundamental structure and reactivity. Photoelectron spectra of PhCC– (1) revealed vertical (VDE) and adiabatic detachment energies (ADE), both of 3.220 eV. The latter value is also the electron affinity (EA) of the phenylethynyl radical (PhCC•) from which a bond dissociation energy (BDE) of phenylacetylene (PhCCH) was derived to be ca. 131.3 ± 1.7 kcal mol-1 using a gas-phase thermochemical cycle. Detachment of an electron from PhCC– at 266 nm likely results in resonant autodetachment, supported by Franck-Condon Factor (FCF) simulations. The phenylacetylide anion reacts with methyl iodide (CH3I) and allyl iodide (C3H5I) via SN2 nucleophilic displacement with measured rate coefficients of 4.29 and 3.66 × 10-10 cm3 molecule-1 s-1, respectively. The mechanisms associated with these displacement reactions are explained and understood in terms of the reaction kinetics, natural bond orbital (NBO) theory, and Density Functional Theory (DFT) calculations.

Ma, Howard Z.↗
DOE OSTI2023

Determining the twist angle of stacked MoS 2 layers using machine learning-assisted low-frequency interlayer Raman fingerprints

The investigation of twisted stacked few-layer MoS 2 has revealed novel electronic, optical, and vibrational properties over an extended period. For the successful integration of twisted stacked few-layer MoS 2 into a wide range of applications, it is crucial to employ a noninvasive, versatile technique for characterizing the layered architecture of these complex structures. Here, in this work, we introduce a machine learning-assisted low-frequency Raman spectroscopy method to characterize the twist angle of few-layer stacked MoS 2 samples. A feedforward neural network (FNN) is utilized to analyze the low-frequency breathing mode as a function of the twist angle. Moreover, using finite difference method (FDM) and density functional theory (DFT) calculations, we show that the low-frequency Raman spectra of MoS 2 are mainly influenced by the effect of the nearest and second nearest layers. A new improved linear chain model (TA-LCM) with taking the twist angle into the consideration is developed to understand the interlayer breathing modes of stacked few-layer MoS 2 . This approach can be extended to other 2D materials systems and provides an intelligent way to investigate naturally stacked and twisted interlayer interactions.

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

Tuning the Molecular Structure and Reaction Mechanism of Olefin Metathesis by Model Bilayered Supported MoO x /AlO x /SiO 2 Catalysts

The molecular structure and activity of supported MoO x olefin metathesis catalysts are heavily impacted by the choice of catalyst support. In this study, surface modification of the SiO 2 support with AlO x and selective anchoring of the MoO x on the surface AlO x sites were used to tune the structure, activation, and reactivity of the resulting surface MoO x sites. Extensive in situ molecular characterization, chemical probe studies, and density functional theory (DFT) calculations reveal that the enhanced activity of the supported MoO x /AlO x /SiO 2 catalyst over the MoO x / SiO 2 catalyst is associated with more favorable activation and kinetics of surface MoO x anchored at AlO x sites.

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

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗
DOE OSTI2021

Radiation-Induced Oxidation Reactions of 2-Selenouracil in Aqueous Solutions: Comparison with Sulfur Analog of Uracil

One-electron oxidation of 2-selenouracil (2-SeU) by hydroxyl (●OH) and azide (●N3) radicals leads to various primary reactive intermediates. Their optical absorption spectra and kinetic characteristics were studied by pulse radiolysis with UV-vis spectrophotometric and conductivity detection and by the density functional theory (DFT) method. The transient absorption spectra recorded in the reactions of ●OH with 2-SeU are dominated by an absorption band with an λmax = 440 nm, the intensity of which depends on the concentration of 2-SeU and pH. Based on the combination of conductometric and DFT studies, the transient absorption band observed both at low and high concentrations of 2-SeU was assigned to the dimeric 2c-3e Se-Se-bonded radical in neutral form (2●). The dimeric radical (2●) is formed in the reaction of a selenyl-type radical (6●) with 2-SeU, and both radicals are in equilibrium with Keq = 1.3 × 104 M−1 at pH 4 (below the pKa of 2-SeU). Similar equilibrium with Keq = 4.4 × 103 M−1 was determined for pH 10 (above the pKa of 2-SeU), which admittedly involves the same radical (6●) but with a dimeric 2c-3e Se-Se bonded radical in anionic form (2●−). In turn, at the lowest concentration of 2-SeU (0.05 mM) and pH 10, the transient absorption spectrum is dominated by an absorption band with an λmax = 390 nm, which was assigned to the ●OH adduct to the double bond at C5 carbon atom (3●) based on DFT calculations. Similar spectral and kinetic features were also observed during the ●N3-induced oxidation of 2-SeU. In principle, our results mostly revealed similarities in one-electron oxidation pathways of 2-SeU and 2-thiouracil (2-TU). The major difference concerns the stability of dimeric radicals with a 2c-3e chalcogen-chalcogen bond in favor of 2-SeU.

2-selenouracil↗
DOE OSTI2025

Mechanistic and kinetic relevance of hydrogen and water in CO 2 hydrogenation on Cu-based catalysts

Here, we ally steady-state kinetics, kinetic isotope effects, and density functional theory (DFT) calculations to illustrate that Cu-based catalysts remain saturated by H-adatoms (H*) and molecular formic acid (HCOOH**) during CO 2 hydrogenation. High H* coverage under methanol synthesis conditions is evidenced by reverse water-gas shift (RWGS) rates that exhibit positive H 2 reaction orders only at P H2 ≲ 0.5 bar, above which methanol synthesis and RWGS rates exhibit first and zeroth order dependence on P H2 , respectively. HCOOH** also accumulates on the surface with increasing P CO2 as informed by the Langmuir-type dependence on P CO2 (0.25-23 bar) for both methanol synthesis and RWGS. As both HCOOH** and H* have one H-atom per site occupied, the two species share the same P H2 dependence and give rise to CO 2 reaction orders that are independent of P H2 . Surface coverages determined based on kinetic analyses are further corroborated with DFT-derived adsorption energies that show favorable HCOOH** adsorbate-adsorbate interactions as well as repulsive interactions for bidentate formate (HCOO**) on H*-saturated surfaces. Methanol selectivity remains invariant with P CO2 and P CO despite CO inhibiting reaction rates, thereby demonstrating methanol synthesis and RWGS occur on the same active site. In contrast, water preferentially inhibits methanol synthesis rates, increases methanol synthesis H 2 reaction order from 1.0 to 1.5, and alters the methanol synthesis H 2 /D 2 kinetic isotope effect; the inhibitory effect of H 2 O thus cannot be attributed to competitive adsorption alone and instead reflects a change in the rate-determining step for methanol synthesis. The disparate kinetics of methanol synthesis and RWGS evince a branching pathway where methanol is formed from formates and CO is formed from carboxylates. The presented work thus identifies the relevant surface species, underscores the distinct catalytic role of water in branching methanol synthesis and RWGS pathways, and, in doing so, details a mechanistic picture that yields predictable rates and reaction orders for both methanol synthesis and RWGS on Cu-based CO 2 hydrogenation catalysts.

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

Electronic tuning of confined sub-nanometer cobalt oxide clusters boosting oxygen catalysis and rechargeable Zn–air batteries.

Reasonable design of robust bifunctional oxygen catalysts from an electronic structure perspective is intriguing and challenging for the development of high active rechargeable zinc-air batteries (ZABs). In this study, the favorable regulation of the electronic structure of the cobalt oxide nanoclusters was firstly predicted by density functional theory (DFT) simulation, and then experimentally verified by confining sub-nanometer CoOx clusters (0.86 nm) into the small pore of ZIF-8 derived N-doped nanomaterials (PNC) using a microporous MOFs confinement strategy. The confined effect of the MOF micropores not only enhanced the stability of the subnanometer cobalt oxide clusters, but also make it coupled with Co-Nx to further regulate the electronic structure of the former, synergistic resulting in enhanced ORR/OER actives. As a result, the optimized 0.05CoOx@PNC catalyst demonstrates outstanding bifunctional oxygen performance with a smaller potential gap of 0.67 V. Moreover, the rechargeable Zn-air batteries integrated 0.05CoOx@PNC air cathode displays encouraging performance with a peak power density of 157.1 mW cm(-2), a specific capacity of 887 mAh g(Zn)(-1)at 10 mA cm(-2) and long-term cyclability for over 200 h, significantly outperforming the benchmark electrode couple consisted of Pt/C/RuO2. DFT calculation further revealed that reducing particle size and coupling with Co-N could effectively regulate the charge distribution of CoOx nanoclusters and downshift the D-band center of Co adsorption sites in CoOx nanoclusters, which reduced the reaction barrier of intermediate O-2* and OH* and ORR/ OER over potential, thus accelerating the overall ORR/OER kinetic process. This work offers a novel reference for the construction of a robust sub-nanometer cluster catalysts in the field of ZABs.

Bifunctional oxygen electrocatalysts↗
DOE OSTI2021

Automated discovery of a robust interatomic potential for aluminum

Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.

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

Effects of water on the kinetics of acetone hydrogenation over Pt and Ru catalysts

We employed an approach combining reaction kinetics measurements at steady state conditions, electronic structure calculations employing density functional theory, and microkinetic modeling for acetone hydrogenation to provide insights into the effects of water on metal catalyst surfaces for the hydrogenation of oxygenates over a wide range of reaction conditions. Elucidation of the repulsive interactions due to adsorbed water molecules at various reaction conditions provides a basis to formulate rate expressions for heterogeneous catalytic processes of biomass oxygenates. Reaction kinetics experiments were carried out at partial pressures of H 2 , acetone, water and helium in the range of 0.51-0.79, 0.02-0.13, 0.08-0.23, 0-0.28 atm, respectively. We show that the addition of water enhances the hydrogenation rate at 353 K and 1 atm on oxophilic metal catalysts such as Ru/C, whereas the same promotional effect of water is not observed for Pt-based catalysts. Microkinetic model predictions for the hydrogenation of acetone on Ru in the absence and presence of water, using enthalpies and entropies obtained from DFT calculations, were in agreement with the experimentally observed reaction orders and activation barriers. The model shows that a water-assisted hydroxypropyl path is expected to be the favored path on Ru with a rate-determining step of H-OH-mediated hydrogenation of C 3 H 6 OH (i.e., the hydroxypropyl intermediate formed by H 2 O-mediated initial hydrogenation of acetone) to produce isopropyl alcohol (IPA). Furthermore, hydrogen, acetone, hydroxypropyl intermediate and hydroxyl species were predicted to be abundant on the Ru surface with a high coverage of nearly 85%. Lastly, the combined studies of computational and experimental catalysis on hydrogenation reactions help to elucidate the mechanistic role of water on metal catalyzed reactions for producing chemical building blocks from biomass-derived oxygenates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗