Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Unveiling the Electrocatalytic Hydrogen Evolution Reaction Pathway on RuP 2 through Ab Initio Grand Canonical Monte Carlo

In this study, the high catalytic reactivity of ruthenium phosphide (RuP 2 ) has been identified by first-principles density functional theory (DFT) calculations for the electrocatalytic hydrogen evolution reaction (HER). Complex surface reconstructions are considered by applying the ab initio grand canonical Monte Carlo (ai-GCMC) algorithm, efficiently providing a sufficient phase-space exploration of possible surfaces. Combined with surface-phase Pourbaix diagrams, we are able to identify the actual surfaces that obtained under specific experimental environments, thus leading to a more accurate understanding of the nature of the active sites and the binding strength of adsorbates. Specifically, through hundreds of surface reconstructions and hydrogenation states generated with ai-GCMC, we identify the most favorable surface phases of RuP 2 under aqueous acidic conditions. We discover that the HER activity is determined by multiple surfaces with different stoichiometries within a narrow electrode potential window. Low HER overpotential (η) has been found for each of the identified surfaces, as low as 0.04 V. High H-coverage reconstructed surfaces have been discovered under acidic conditions, and the surface Ru sites introduced by additional Ru adatoms or exposed by P-vacancies serve as the active sites for HER based on their nearly reversible H binding. Furthermore, this work provides atomistic insights into the origin of high HER activity on RuP 2 by exploring the dynamic surface phases of electrocatalysts and features a generalizable method to explore the reconstructed/hydrogenated surface space as a function of experimental conditions.

25 ENERGY STORAGE↗
DOE OSTI2025

Effects of co-adsorbed water on different bond cleavages involved in acetic acid decomposition on Pt (111)

Acetic acid decomposition on Pt (111) in the presence of co-adsorbed water is a good model system for oxygenate decomposition on Pt (111) in aqueous phase, with application in hydrogen production and biomass conversion. Here, in this study, we present a density functional theory (DFT) theory calculation of how co-adsorbed water affects different bond cleavages of acetic acid decomposition on Pt (111). The presence of co-adsorbed water generally enhances O$-$H bond cleavage while inhibiting OC$-$O and OC$-$OH bond cleavage. The influence of co-adsorbed water on C$-$H bond cleavage varies the most and depends on the nature of the transition state and how co-adsorbed water stabilizes the initial and final state. Although these trends are useful as general guidance, they are not sufficient to predict the effect on a complex reaction network such as acetic acid decomposition on Pt (111) which has several parallel reaction paths with similar energies. In the absence of co-adsorbed water, the two lowest energy pathways are decarboxylation (DCX) and decarbonylation (DCN) pathways through a common CH 2 COO intermediate, in which the DCX pathway (C$-$H bond cleavage of CH 2 COO) is more favorable than DCN pathway (OC$-$O bond cleavage of CH 2 COO). In the presence of co-adsorbed water, the energy difference between Csingle bondH bond cleavage and OC$-$O bond cleavage of CH 2 COO increases, suggesting an increase in the favorability of acetic acid decarboxylation (formation of carbon dioxide) over acetic acid decarbonylation (formation of carbon monoxide) on Pt (111).

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

First-principles study of the structures and redox mechanisms of Ni-rich lithium nickel manganese cobalt oxides

To reduce the cobalt (Co) content in lithium-ion batteries, Ni-rich (high-Ni) lithium nickel manganese cobalt oxides (NMC) are pursued as one of the next-generation cathode materials. However, there is still debate on the crystal and electronic structures of the baseline, LiNiO 2 . Density Functional Theory (DFT) calculations were performed to provide a theoretical understanding of Ni-rich NMC. First, it was found that the commonly used $R\bar{3}m$ structure for LiNiO 2 is metallic, contrary to the experimentally reported mix-conducting behavior. Among the four different space groups, $R\bar{3}m$, C2/m, P2 1 /c, and P2/c, P2/c with charge disproportionation of Ni 2+ and Ni 4+ is the most energetically stable and semiconducting structure of LiNiO 2 . Therefore, the atomic structures of representative Ni-rich NMC were built by partially replacing Ni with Co or Mn in the P2/c LiNiO 2 to form Li x Ni y Mn z Co 1-y-z O 2 . In the fully lithiated (x=1.0) high Ni content NMC (y>0.5), the oxidation state of all Mn ions becomes 4+, while Co ions still maintain 3+, and part of the Ni ions become 3+ to compensate for the charge. Upon delithiation, the local environment shows more variation of the charge states on the transition metal (TM) ions. The average oxidation on each TM follows a sequence of losing electrons that starts from Ni 2+ to Ni 3+ , then oxidizing Ni 3+ and Co 3+ , while Mn 4+ remains electrochemically inactive till x=0. Finally, a general relationship for the oxidation state change in each TM as a function of x is derived and shows agreement with both modeling and experimental data.

25 ENERGY STORAGE↗
DOE OSTI2021

Predicting and Synthesizing Interface Stabilized 2D Layers

The compound (Pb 2 MnSe 3 ) 0.6 VSe 2 was predicted to be kinetically stable based on density functional theory (DFT) calculations on an island of Pb 2 MnSe 3 between layers of VSe 2 . This approach provides a high degree of freedom by not forcing interlayer lattice match, making it ideal to investigate the likelihood of formation of new incommensurate layer misfit structures. The free space around the island is critical, as it allows atoms to diffuse and hence exploring the local energy landscape around the initial configuration. (Pb 2 MnSe 3 ) 0.6 VSe 2 was synthesized via a near diffusionless reaction from precursors where a repeating sequence of elemental layers matches the local composition and layer sequence of the predicted compound. The VSe 2 layer consists of a Se–V–Se trilayer with octahedral coordination of the V atoms. The Pb 2 MnSe 3 layer consists of three rock-salt-like planes, with a MnSe layer between the planes of PbSe. The center MnSe plane stabilizes the puckering of the outer PbSe layers. Electrical properties indicate that (Pb 2 Mn 1 Se 3 ) 0.6 VSe 2 undergoes a charge density wave transition at ~100 K and orders ferromagnetically at 35 K. Overall, the combination of theory and experiment enables a faster convergence to new heterostructures than either approach in isolation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Model for Atomic Precision p-Type Doping with Diborane on Si(100)-2×1

Diborane (B 2 H 6 ) is a promising molecular precursor for atomic precision p-type doping of silicon that has recently been experimentally demonstrated. In this study, we use density functional theory (DFT) calculations to determine the reaction pathway for diborane dissociating into a species that will incorporate as electrically active substitutional boron after adsorbing onto the Si(100)-2×1 surface. Our calculations indicate that diborane must overcome an energy barrier to adsorb, explaining the experimentally observed low sticking coefficient (<1 × 10 –4 at room temperature) and suggesting that heating can be used to increase the adsorption rate. Upon sticking, diborane has an ≈50% chance of splitting into two BH 3 fragments versus merely losing hydrogen to form a dimer such as B 2 H 4 . As boron dimers are likely electrically inactive, whether this latter reaction occurs is shown to be predictive of the incorporation rate. The dissociation process proceeds with significant energy barriers, necessitating the use of high temperatures for incorporation. Using the barriers calculated from DFT, we parameterize a Kinetic Monte Carlo model that predicts the incorporation statistics of boron as a function of the initial depassivation geometry, dose, and anneal temperature. Our results suggest that the dimer nature of diborane inherently limits its doping density as an acceptor precursor and furthermore that heating the boron dimers to split before exposure to silicon can lead to poor selectivity on hydrogen and halogen resists. This suggests that, while diborane works as an atomic precision acceptor precursor, other non-dimerized acceptor precursors may lead to higher incorporation rates at lower temperatures.

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

Influence of Al location on formation of silver clusters in mordenite

Formation of zeolite supported Ag 0 clusters depends on a combination of thermodynamically stable atomic configurations, charge balance considerations, and mobility of species on the surface and within pores. Periodic density functional theory (DFT) calculations were performed to evaluate how the location of Al in the mordenite (MOR) framework and humidity control Ag 0 nanocluster formation. Four Al framework sites were studied (T1-T4) and the Al positions in the framework were identified by the shifts in the differential Al…Al pair distribution function (PDF). Furthermore, structural information about the Ag 0 nanoclusters, such as dangling bonds, can be identified by Ag…Ag PDF data. For Ag 0 formation in vacuum MOR structures with a Si:Al ratio of 5:1 with Al in the T1 position resulted in the most framework flexibility and the lowest Ag 0 nanocluster charge, indicating the best result for formation of charge neutral nanoclusters. When water is present, Al in the T3 and T4 positions results in the formation of the smallest average Ag 0 nanoclusters plus greater expansion of the O-T-O bond angle than in vacuum, indicating easier diffusion of the Ag 0 nanoclusters to the surface. Here, the presence of Al in 4-membered rings and in pairs indicates favorable MOR structures for formation of single Ag atoms, despite the existence of synthesis challenges. Therefore, Al in the T2 position is the least favorable for Ag 0 nanocluster formation in both vacuum and in the presence of water. Al in the T1, T3, and T4 positions provides beneficial effects through framework flexibility and changes in nanocluster size or charge that can be leveraged for design of zeolites for formation of metallic nanoclusters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Orthorhombic Structures as Inorganic Halide Perovskite Models for High-Throughput DFT Investigations

Although high-symmetry Pm-3m space group cubic models are computationally efficient for high-throughput density functional theory (DFT) calculations of inorganic ternary (ABX 3 ) halide perovskites (HPs), they frequently predict band gaps (E g ) that disagree with experiment. Conversely, while low-symmetry cubic polymorphous networks (PN) comprised of 160 to 320 atoms incur significantly greater computational cost, they predict E g ’s that are more closely aligned with experiment. Here, in this study, we compare the DFT total energies and E g ’s predicted by four high-symmetry structure models (Pnma orthorhombic, R3m trigonal, P4/mbm tetragonal, and Pm-3m cubic) to cubic PNs for 5 experimentally characterized ternary HPs and find that the orthorhombic model computes Eg’s with the smallest MAD of 0.23 eV relative to the PNs. Pair distribution functions and DFT-computed total energies show that octahedral tilting, which is present in the 20-atom orthorhombic and 160-atom cubic PN models but not in the 5-atom cubic models, stabilizes all 5 compositions in our benchmarking set. We also find that imposing PN constraints when generating and optimizing these orthorhombic structures by fixing the unit cell lattice vectors and displacing the atoms prior to ionic relaxation with DFT, which we call the orthorhombic surrogate model (OSM), lowers the MAD of E g predictions to 0.09 eV. Our OSM predicts the PN band gaps of an additional 95 theoretical inorganic ternary HPs with MAD of 0.08 eV, supporting its usage in high-throughput DFT investigations to closely estimate PN band gaps with much less computational expense.

14 SOLAR ENERGY↗
DOE OSTI2020

Kinetics of low-temperature methane activation on IrO 2 (1 1 0): Role of local surface hydroxide species

The ability of the IrO 2 (1 1 0) surface to promote CH 4 activation at low temperatures (~150 K) suggests possibilities for developing IrO 2 -based catalysts to selectively convert light alkanes to high-value chemicals. In this study, we present experimental results and microkinetic simulations, based on density functional theory (DFT) calculations, of temperature programmed reaction spectra (TPRS) obtained for CH 4 σ-complexes adsorbed on IrO 2 (1 1 0), and focus on clarifying how surface OH groups modify the branching between CH 4 desorption, activation and subsequent reactions during TPRS. Our DFT results predict that surface OH groups strongly destabilize CH 4 σ-complexes on IrO 2 (1 1 0) in addition to deactivating surface O-atoms that are needed to achieve CH 4 activation at low temperature. We demonstrate that a microkinetic model that incorporates the influence of surface OH groups on the CH 4 binding and reactivity reproduces experimental TPRS results which show that adsorbed CH 4 σ-complexes preferentially dissociate at low CH 4 coverage, but that an increasing fraction of the adsorbed CH 4 desorbs at low temperature (~120 K) with increasing initial CH 4 or H coverage. The simulations reveal that low-temperature CH4 desorption during TPRS arises primarily from CH 4 σ-complexes that are kinetically-trapped between adjacent surface OH groups, and are thus destabilized and unable to access reactive O-atoms. In addition, when the effect of adjacent OH groups is incorporated we find that the PBE-D3 functional incorporating dispersion provides CH 4 binding energies that agree more closely with that obtained from the TPRS experiments. Overall, our results demonstrate that the local effect of adjacent OH groups must be incorporated into any microkinetic models to properly capture the selectivity between extensive and partial oxidation in alkane conversion on IrO 2 (1 1 0) under reaction conditions.

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

Methanol tolerance of atomically dispersed single metal site catalysts: mechanistic understanding and high-performance direct methanol fuel cells

Proton-exchange membrane fuel cells (PEMFCs) and direct methanol fuel cells (DMFCs) are promising power sources from portable electronic devices to vehicles. The high-cost issue of these low-temperature fuel cells can be primarily addressed by using platinum-group metal (PGM)-free oxygen reduction reaction (ORR) catalysts, in particular atomically dispersed metal–nitrogen–carbon (M–N–C, M = Fe, Co, Mn). Furthermore, a significant advantage of M–N–C catalysts is their superior methanol tolerance over Pt, which can mitigate the methanol cross-over effect and offer great potential of using a higher concentration of methanol in DMFCs. In this work, we investigated the ORR catalytic properties of M–N–C catalysts in methanol-containing acidic electrolytes via experiments and density functional theory (DFT) calculations. FeN 4 sites demonstrated the highest methanol tolerance ability when compared to metal-free pyridinic N, CoN 4 , and MnN 4 active sites. The methanol adsorption on MN 4 sites is even strengthened when electrode potentials are applied during the ORR. The negative influence of methanol adsorption becomes significant for methanol concentrations higher than 2.0 M. However, the methanol adsorption does not affect the 4e - ORR pathway or chemically destroy the FeN 4 sites. The understanding of the methanol-induced ORR activity loss guides the design of promising M–N–C cathode catalyst in DMFCs. Accordingly, we developed a dual-metal site Fe/Co–N–C catalyst through a combined chemical-doping and adsorption strategy. Instead of generating a possible synergistic effect, the introduced Co atoms in the first doping step act as “scissors” for Zn removal in metal–organic frameworks (MOFs), which is crucial for modifying the porosity of the catalyst and providing more defects for stabilizing the active FeN 4 sites generated in the second adsorption step. The Fe/Co–N–C catalyst significantly improved the ORR catalytic activity and delivered remarkably enhanced peak power densities (i.e., 502 and 135 mW cm -2 ) under H 2 –air and methanol–air conditions, respectively, representing the best performance for both types of fuel cells. Notably, the fundamental understanding of methanol tolerance, along with the encouraging DMFC performance, will open an avenue for the potential application of atomically dispersed M–N–C catalysts in other direct alcohol or ammonia fuel cells.

25 ENERGY STORAGE↗
DOE OSTI2021

Mechanistic Investigation of Dimethylmercury Formation Mediated by a Sulfide Mineral Surface

Mercury (Hg) pollution is a global environmental problem, but many key transformations among 29 various Hg species are not well understood. Although anaerobic microorganisms are the primary 30 producers of monomethylmercury (MMeHg), the abiotic formation of dimethylmercury (DMeHg) 31 from MMeHg have been suggested to account for a large portion of DMeHg in the ocean. It has 32 been shown that abiotic formation of DMeHg from MMeHg can be facilitated by reduced sulfur 33 groups on mineral surfaces. However, a detailed mechanistic understanding of this reaction is 34 lacking. Here we perform density functional theory (DFT) calculations to explore mechanisms of 35 DMeHg formation on the 110 surface of a CdS(s) (hawleyite) nanoparticle. We show that 36 coordination of MMeHg substituents to adjacent reduced sulfur groups protruding from the surface 37 facilitate DMeHg formation and that the reaction proceeds through direct transmethylation from 38 one MMeHg substituent to another. Solvation effects were found to play an important role in the 39 surface reconstruction of the nanoparticle and in decreasing the energetic barrier for DMeHg 40 formation. We also investigate DMeHg formation from MMeHg mediated by organic dithiols and 41 compare the results to those mediated by the mineral surface. In contrast to the surface reaction, 42 the calculations suggest that dithiol-mediated formation of DMeHg instead proceeds through a 43 mechanism in which Hg–Hg interactions activate a methyl group from one MMeHg substituent 44 for intramolecular transfer to the adjacent MMeHg substituent. These findings refine previously 45 proposed mechanisms of abiotic DMeHg formation and fill gaps in our understanding of global 46 Hg cycling in the environment.

Lian, Peng↗
DOE OSTI2022

Proton-Coupled Electron Transfer in a Ruthenium(II) Bipyrimidine Complex in Its Ground and Excited Electronic States

In this work, proton-coupled electron transfer (PCET) was studied for the ground and excited electronic states of a [Ru(terpy)(bpm)(OH 2 )(PF 6 ) 2 ] complex, Ru-bpm. Cyclic voltammetry measurements show that the Ru(II)-aqua moiety undergoes PCET to form a Ru(IV)-oxo moiety in the anodic region, while the bpm ligand undergoes PCET to form bpmH 2 in the cathodic region. The photophysical behavior of Ru-bpm was studied using steady-state and femtosecond transient UV–vis absorption spectroscopy, coupled with density functional theory (DFT) calculations. The lowest-lying excited state of Ru-bpm is described as a (Ru → bpm) metal-to-ligand charge-transfer (MLCT) state, while the metal-centered (MC) excited state was found computationally to be close in energy to the lowest-energy bright MLCT state (MC state was 0.16 eV above the MLCT state). The excited-state kinetics of Ru-bpm were found via transient absorption spectroscopy to be short-lived and were fit well to a biexponential function with lifetimes $τ_1$ = 4 ps and $τ_2$ = 65 ps in aqueous solution. Kinetic isotope effects of 1.75 ($τ_1$) and 1.61 ($τ_2$) were observed for both decay components, indicating that the solvent plays an important role in the excited-state dynamics of Ru-bpm. Based on the pH-dependent studies and the results from prior studies of similar Ru-complexes, we hypothesize that the 3 MLCT state forms an excited-state hydrogen-bond adduct with the solvent molecules and that this process occurs with a 4 ps lifetime. The formation of such a hydrogen-bond complex is consistent with the electronic density accumulation at the peripheral N atoms of the bpm moiety in the 3 MLCT state. The hydrogen-bonded state 3 MLCT decays to the ground state with a 65 ps lifetime. Such a short lifetime is likely associated with the efficient vibrational energy transfer from the 3 MLCT state to the solvent.

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

Engineering bimetallic interfaces and revealing the mechanism for carbon dioxide electroreduction to C 3+ liquid chemicals

Reduction reaction of CO 2 (CO 2 RR) to liquid C 3+ chemicals is a net-zero-carbon process and can increase local resiliency to power outages and fuel consumption. The mechanism and the catalyst design rules of CO 2 RR-to-C 3+ are unknown. Engineering bimetallic interface (e.g., Pd/Au) to tune the intermediate adsorption is promising for promoting C 3+ formation. Our density functional theory (DFT) calculations find that *CH 2 could be the key intermediate and C 1 -CH 2 coupling could be the rate-limiting step to generate C 3+ . High CO surface coverages can promote the bimetallic interfacial sites, lower the energetics of the C 1 -CH 2 coupling step, and enhance the C 3+ formation. We further construct a volcano plot of C 1 -CH 2 kinetics as a function of the binding strength of key intermediate *CH 2 via engineering the d-band center of the interfacial site. Our findings could guide the rational design of bimetallic interfaces and their near-surface microenvironment for enhancing CO 2 RR-to-C 3+ .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Bimetallic-Derived Catalytic Structures for CO 2 -Assisted Ethane Activation

In recent years, the simultaneous upgrading of CO 2 and ethane has emerged as a promising approach for generating valuable gaseous (CO, H 2 , and ethylene) and liquid (aromatics and C 3 oxygenates) chemicals from the greenhouse gas CO 2 and large-reserved shale gas. The key challenges for controlling product selectivity lie in the selective C-H and C-C bond cleavage of ethane with the assistance of CO 2 . Bimetallic-derived catalysts likely undergo alloying or oxygen-induced segregation under reaction conditions, thus providing diverse types of interfacial sites, e.g., metal/support (M/M'O x ) interface and metal oxide/metal (M'O x /M) inverse interface, that are beneficial for selective CO 2 -assisted ethane upgrading. The alloying extent can be initially predicted by cohesive energy and atomic radius (or Wigner-Seitz radius), while the preference for segregation to form the on-top suboxide can be approximated using the work function, electronegativity, and binding strength of adsorbed oxygen. Furthermore, bimetallic-derived catalysts are typically supported on high surface area oxides. Modifying the reducibility and acidity/basicity of the oxide supports and introducing surface defects facilitate CO 2 activation and oxygen supplies for ethane activation. Using in situ synchrotron characterization and density functional theory (DFT) calculations, we found that the electronic properties of oxygen species influence the selective cleavage of C-H/C-C bonds in ethane, with electron-deficient oxygen over the metal (or alloy) surface promoting nonselective bond scission to produce syngas and electron-enriched oxygen over the metal oxide/metal interface enhancing selective C-H scission to yield ethylene. Here, we further demonstrate that the preferred structures of the catalyst surfaces, either alloy surfaces or metal oxide/metal inverse interfaces, can be controlled through the appropriate choice of metal combinations and their atomic ratios. Through a comprehensive comparison of experimental results and DFT calculations, the selectivity of C-C/C-H bond scission is correlated with the thermodynamically favorable bimetallic-derived structures (i.e., alloy surfaces or metal oxide/metal inverse interfaces) under reaction conditions over a wide range of bimetallic catalysts. These findings not only offer structural and mechanistic insights into bimetallic-derived catalysts but also provide design principles for selective catalysts for CO 2 -assisted activation of ethane and other light alkanes. This Account concludes by discussing challenges and opportunities in designing advanced bimetallic-derived catalysts, incorporating new reaction chemistries for other products, employing precise synthesis strategies for well-defined structures with optimized site densities, and leveraging time/spatial/energy-resolved in situ spectroscopy/scattering/microscopy techniques for comprehensive structural analysis. The research methodologies established here are helpful for the investigation of dynamic alloy and interfacial structures and should inspire more efforts toward the simultaneous upgrading of CO 2 and shale gas.

03 NATURAL GAS↗
DOE OSTI2026

Investigating trends in actinide covalency and magnetism with 35/37 Cl SSNMR spectroscopy and first-principles calculations

Solid-state NMR (SSNMR) spectroscopy is a powerful technique for studying actinide chemistry but has been significantly limited due to the complex paramagnetism, and radiological hazards presented by these materials. Lanthanide and actinide salts often feature magnetic ordering and can be paramagnetic, ferromagnetic, or antiferromagnetic depending on temperature and electronic structure. Paramagnetic interactions can manifest in SSNMR both as secular spectral shifts and/or couplings as well as contributions from non-secular relaxation. Both effects can be directly measured with NMR and used to extrapolate rich chemical information such as coordination environments, bonding characteristics, local molecular dynamics, and correlation times. Typically, these studies are carried out on high-γ and highly abundant NMR-active isotopes (e.g., 1 H, 6/7 Li, 19 F, 23 Na, etc.) or on enriched rare isotopes (e.g., 2 H and 17 O), which can be expensive. Herein, we present a facile methodology to measure the 35/37 Cl electric-field gradient (EFG) and paramagnetic shift anisotropy (SA) tensor components using static wideline SSNMR measurements of LaCl 3 , NdCl 3 , UCl 3 , and UCl 4 . The static powder spectra were measured with both 35 Cl and 37 Cl SSNMR to increase the fidelity of the extracted tensor parameters. Variable temperature NMR of a select case confirms the Curie–Weiss paramagnetism. Relaxation measurements of both nuclei further corroborate observations owing to the paramagnetic relaxation enhancement and reveal simultaneous quadrupolar relaxation mechanisms. Density functional theory (DFT) calculations using Hubbard U corrections to the uranium valence orbitals show excellent agreement with experimental EFG tensor parameters and help describe the bonding characteristics in these lanthanide and actinide systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Gaseous cyclodextrin- closo -dodecaborate complexes χCD·B 12 X 12 2– (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin- closo -dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo -dodecaborate dianions B 12 X 12 2– (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B 12 X 12 2– ) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B 12 X 12 2– prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e. , X = Cl, Br, and I, are found outside of α-CD, while B 12 F 12 2– is completely wrapped by γ-CD. Partial embedment of B 12 X 12 2– into CDs is observed for the other complexes via multipronged B–X … H–O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host–guest interactions for the X = F series partially due to the formation of a O/C–H … X–B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

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

Assessing the Role of Interfacial and Metal Sites in Pt/TiO2-Catalyzed Acetic Acid Hydrodeoxygenation

Thermochemical conversion of biomass to produce drop-in quality biofuels typically involves hydrodeoxygenation (HDO) steps following catalytic fast pyrolysis (CFP) to remove excess oxygen and create a more-stable bio-oil product. HDO reactions are performed by co-feeding the CFP vapor-phase product and H2 gas over a bi-functional catalyst. Noble-metal catalysts supported on reducible metal oxides (e.g., Pt/TiO2) are active and selective toward these HDO reactions. Griffin and co-workers showed that Pt/TiO2 catalysts promote the desired deoxygenation steps for m-cresol HDO while mitigating undesired C-C bond-breaking steps that reduce the overall value/energy density of the biofuel. Such model-compound studies to inform the design of improved catalysts for HDO chemistry are necessary to improve the overall process economics/efficiencies for biofuels production. One important class of bio-derived compounds that has not been studied extensively with respect to HDO chemistry, particularly at the atomic level, is carboxylic acids. Past experimental work indicates that Pt/C and Pt/TiO2 catalysts are selective toward C-C and C-O bond-dissociation products for acetic acid HDO, respectively. To understand the role of Pt/TiO2 active sites in this observed change in selectivity and guide catalyst development, the work in this presentation focuses on modeling the role of Pt-metal and Pt-TiO2-interfacial sites in promoting key C-C bond-breaking, C-O bond-breaking, and (de)hydrogenation steps in acetic acid HDO. Density functional theory (DFT) calculations were performed using the Vienna Ab initio Simulation Package (VASP). The exchange-correlation functional was approximated by the Perdew-Burke-Ernzerhof functional. Dispersion interactions were captured using the D3 method. Projector augmented-wave potentials described electron-ion interactions, and electron wavefunctions were expanded via a planewave basis with an energy cutoff of 500 eV. Activation barriers were calculated using the climbing image nudged elastic band method. Results and Discussion: To discern the role of Pt-metal and Pt-TiO2 interface sites in promoting acetic acid HDO chemistry, Pt(111) slab and anatase TiO2(101)-supported Pt-nanowire (PtNW/OH-TiO2) surface models were constructed, respectively. Because H2 is co-fed in HDO reactions, the anatase support was terminated with OH groups. Interfacial vacancies have been shown to facilitate Ru/TiO2-catalyzed phenol HDO; thus, an interfacial model with an OH vacancy was also considered (PtNW/OHv-TiO2). Pt-TiO2-interface sites stabilize adsorption of all studied acetic acid HDO surface intermediates relative to terrace Pt-metal sites, particularly when an interfacial-OH vacancy is present. Oxygenated species prefer to bind at the OH vacancy through the O atom, suggesting a preference for C-O over C-C bond cleavage at these sites. This hypothesis is supported by net-negative and net-positive average shifts in the reaction energy and activation energy barriers for C-O and C-C bond-breaking steps, respectively, at Pt-TiO2-interface sites relative to Pt-metal sites. Using the calculated energetics, the predicted minimum-energy pathway was determined for each surface. Pt(111) is predicted to follow decarboxylation, producing undesired methane and carbon dioxide. Conversely, PtNW/OH-TiO2 and PtNW/OHv-TiO2 are both predicted to produce desired acetaldehyde and ethane. The interfacial vacancy may also play a key role in facilitating the first C-O bond-breaking step in acetic acid HDO, lowering the barrier by 0.6 eV relative to the defect-free interface model. These results demonstrate the critical role of the Pt-TiO2 interface in the shift in acetic acid HDO selectivity experimentally observed on Pt/C and Pt/TiO2 catalysts. The results herein demonstrate the important role of the Pt-TiO2 interface and interfacial oxygen vacancies in improving the carbon efficiency for HDO reactions in CFP upgrading.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2025

Electronic density of states as the descriptor of elastic bond strength, ductility, and local lattice distortion in BCC refractory alloys

Although electronic density of states (DOS) is fundamental to materials properties, its general relationship to mechanical properties of alloys is not well established. In this paper, using density functional theory (DFT) calculations, we show that the electronic occupancy at the Fermi level, N(E f ), obtained from DOS is a key descriptor of alloy strength and ductility. Our comprehensive analysis of numerous body centered cubic (BCC) refractory high entropy alloys (RHEAs) shows an overwhelming correlation that low N(E f ) indicates strong bonds that have high stiffness resulting in high elastic constants. High bond stiffness indicates presence of covalent nature of bonds that are directional in nature resulting in resistance to deformation leading to high bulk (B) and shear (G) moduli. Consequently, N(E f ) provides a direct correlation to the tendency of alloy ductility evidenced in the Pugh ratio (G/B). As stiffer bonds result in lower local lattice distortion (LLD), N(E f ) are LLD are also found to be corelated which opens up a correlation to solid solution strengthening and yield strength. Thus, this work unveils fundamental correlations between N(E f ) and (1) elastic bond strength, (2) ductility, and (3) LLD. These correlations open opportunities for the design of high strength high ductile RHEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High spin polarization in the disordered quaternary Heusler alloy FeMnVGa

Here in this paper, we report the successful synthesis of a Fe-based highly spin-polarized quaternary Heusler alloy FeMnVGa and its structural, magnetic, and transport properties probed through different experimental methods and theoretical techniques. Density functional theory (DFT) calculations performed on different types of structures reveal that the structure with Ga at 4a, V at 4b, Mn at 4c, and Fe at 4d (space group F$\bar{4}$3m) possess minimum energy among all the ordered variants. Ab initio simulations in the most stable ordered structure show that the compound is a ferromagnet having a large spin-polarization (89.9%). Neutron diffraction reveals that the compound crystallizes in disordered type-2 structure (space group Fm$\bar{3}$m) in which Ga occupies at 4a, V 4b and Fe/Mn occupy 4c/4d sites with 50:50 proportions. The structural disorder is further confirmed by x-ray diffraction, extended x-ray absorption fine structure, 57 Fe Mössbauer spectrometry results, and DFT calculations. Magnetization studies suggest that the compound orders ferromagnetically below T C ~ 293 K and the saturation magnetization follows the Slater-Pauling rule. Mössbauer spectrometry, along with neutron diffraction, suggest that Mn is the major contributor to the total magnetism in the compound, consistent with the theoretical calculations, which also indicates that spin polarization remains high (81.3%), even in the presence of such large atomic disorder. The robustness of the half-metallic ferromagnetic (HMF) property in the presence of disorder is a quite unique characteristic over other reported HMF in literature and makes this compound quite promising for spintronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗