Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Investigating the Elusive Nature of Atomic O from CO 2 Dissociation on Pd(111): The Role of Surface Hydrogen

CO 2 dissociation is a key step in CO 2 conversion reactions to produce value-added chemicals typically through hydrogenation. In many cases, the atomic O produced from CO 2 dissociation can potentially block adsorption sites or change the oxidation state of the catalyst. Here, we used ambient pressure X-ray photoelectron spectroscopy (AP-XPS) and density functional theory (DFT) calculations to investigate the presence of surface species from the dissociation of CO 2 on Pd(111). AP-XPS results show that CO 2 was dissociated to produce adsorbed CO, but dissociated atomic O was not observed at room temperature. We were only able to observe atomic O when CO 2 was introduced at 500 K. Further investigations of O-covered Pd(111) revealed that chemisorbed O could be easily removed by low pressures of CO and H 2 . Notably, the effect of H 2 is quite prominent since it could react with chemisorbed O at a pressure as low as 2 × 10 -9 Torr, and the presence of H2 at ambient pressure prevented CO 2 dissociation. DFT calculations showed that in the presence of background H 2 , facile CO 2 dissociation took place via the reverse water–gas shift (rWGS) reaction, which resulted in the formation of adsorbed CO and removal of O by H 2 . DFT also identified the possible variation of surface species on simultaneous exposure of CO 2 and H 2 over Pd(111) depending on temperature and pressure, which opens alternative opportunities to tune the CO 2 hydrogenation catalysis by controlling the reaction conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Point Defects Control Guest Molecule Diffusion in the 1D Pores of Zn(tbip)

Molecular diffusion plays a critical role in metal-organic frameworks (MOFs) within the application of kinetic chemical separations. We carefully study in this work the unexpected role of point defects for short-chain alkanes diffusing in Zn(tbip), an MOF with rigid one-dimensional (1D) channels. Inside a defect-free Zn(tbip), guest molecules are expected to follow single-file diffusion along 1D channels. It has been found previously by Heinke et al. that these parallel 1D channels are connected for molecular diffusion. Our density functional theory (DFT) calculations suggest that linker vacancy defects could arise under experimentally relevant conditions by removing a pair of linkers. Further climbing-image nudged elastic band (cNEB) DFT calculations indicate that hopping of short-chain alkanes between adjacent 1D channels over defect windows can occur at moderate temperatures. In addition to providing insights into connected adjacent 1D channels in Zn(tbip), Heinke et al. also inferred that most 1D pores are blocked from a microkinetic model. Additionally, we explored the influence of hydrolyzed linker created by the formation of linker removal inside 1D pores. Our DFT calculations show that the linkers can effectively block the pores and the linker diffusion in 1D channels is slow. Our results, for the first time, offer a mechanistic explanation of the unexpected molecular diffusion behavior in this MOF with 1D channels.

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

Determining Catalytically Relevant Surfaces through Coverage-Dependent Lattice Gas Models: Carbon Adsorption on Fe(100)

Here, we have quantified the C–C lateral interactions on Fe(100) using a density functional theory (DFT)-parameterized lattice gas cluster expansion (LG CE) model trained using 265 unique configurations spanning a C coverage from 0 to 1 monolayer (ML). Our LG CE model shows high predictive accuracy with a leave-multiple-out cross-validation score of 10.2 and 16.6 meV/site for systems with and without the top two layers of Fe atoms fixed, respectively. Electronic ground-state structures identified from the lattice gas model (including the structures at 0 and 1 monolayers) were further used to generate ab initio phase diagrams under a range of temperatures and pressures. At low temperatures (<400 K), we found that the 1.0 monolayer structure is dominant, whereas at higher temperatures (>500 K), the 0.88 ML structure is most likely to form on the Fe surface. Interestingly, our model identified a c (2 × 2) ordered structure at 1/2 ML, which correlates well with previous DFT studies for carbon adsorption on iron surfaces and matches with the experimentally observed low-energy electron diffraction structure. Overall, the DFT-parameterized energies for the C/Fe system including effects of coverage and configurational space can further help in developing multiscale models for various heterogeneous reactions involving C–C and C–Fe interactions.

08 HYDROGEN↗
DOE OSTI2023

Ab Initio Molecular Dynamics Spectra for Characterization of Hydrated Al 2 O 3 Supported MoO x

The structure of supported MoO x /Al 2 O 3 catalysts is investigated using ab initio molecular dynamics (AIMD) and density-functional theory (DFT). Phase diagrams were computed to understand the hydroxylation coverage as a function of the temperature, H 2 O partial pressure, and MoO x loading. We relate the shifts in experimental Raman vibrational frequencies under hydrated conditions to hydrogen bonding interactions and the MoO x anchoring location. We showcase that the use of AIMD as a benchmark for DFT can provide insight into how, under certain hydrated conditions, DFT excels as an inexpensive method for computing vibrational frequencies, while, under dehydrated conditions, it is susceptible to the largest errors. Additionally, to facilitate the analysis of the hydroxyl region of experimental infrared spectra, we compute the power spectra of individual hydroxyls and see a strong relationship between OH stretching frequencies and the surrounding coordination environment, which are also impacted by anchored MoO x . Here, we compare computational and experimental IR and Raman spectra of catalysts synthesized herein under the same conditions.

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

Unveiling the Role of Surface Ir-Oxo Species in O 2 Evolution at IrO 2 Electrocatalysts via Embedded Cluster Multireference Calculations

Understanding the mechanisms driving the oxygen evolution reaction (OER) on iridium oxide (IrO 2 )-based catalysts is essential to improving their performance and enabling an actual scale-up of water-splitting photoelectrochemical cells. The mechanistic pathways at IrO 2 interfaces have been extensively investigated computationally using density functional theory (DFT), which predicts a high-energy barrier for the last step of the OER of molecular oxygen detachment and release from the catalyst surface. Nevertheless, surface O 2 over- and under-binding results by standard generalized gradient approximation and hybrid density functionals, respectively, call for further analysis of this crucial step via multireference methods. Aiming at unveiling the nature of such a barrier, we hereby address the formation of O 2 from the most-stable IrO 2 (110) surface with both periodic DFT and an electrostatic embedded cluster approach at the n-electron valence-state perturbation theory. With this multireference approach, we find a value for the aforementioned energy barrier that is much closer to experimental indications than DFT ones. Here, an in-depth analysis of the involved molecular orbitals suggests that the origin of this barrier is related to the breaking of a π interaction between O 2 and Ir surface atom and to a significant additional O 2 interaction with adjacent electrophilic Ir-oxo species, which is present under experimental operating conditions. Besides shedding light on the mechanism of the OER on IrO 2 , these findings point out the importance of multireference methods for dissecting complex reactions at electrocatalytic interfaces and pave the route for further investigations with effective embedding approaches.

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

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

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

Efficient Modeling of Structural, Electronic, and Optical Properties of Silver and Gold Metal Nanoclusters and Alloys Using Optimized SCC-DFTB Parameters

Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.

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

Quantum Monte Carlo Benchmarking of Molecular Adsorption on Graphene-Supported Single Pt Atom

The precise understanding of adsorption energetics and molecular geometry at catalytic sites is fundamental for advancing catalysis, particularly under the constraints of resource efficiency and environmental sustainability. Here, this study benchmarks the performance of density functional theory (DFT) calculations against diffusion Monte Carlo (DMC) calculations for adsorption properties of small gas molecules relevant to CO oxidation—namely O 2 , CO, CO 2 , and atomic oxygen—on a single Pt atom supported by pristine graphene. Our findings reveal that DMC calculations provide a significantly different landscape of adsorption energetics compared to DFT results. Notably, DFT predicts different lowest-energy configurations and spin states, particularly for O 2 , which suggests potential discrepancies in predicting the catalytic behavior. Furthermore, this study identifies the critical issue of CO poisoning, highlighted by the large disparity between the DMC adsorption energies of O 2 (−1.23(2) eV) and CO (−3.37(1) eV), which can inhibit the catalytic process. These results emphasize the necessity for more sophisticated computational approaches in catalysis research, aiming to refine the prediction accuracy of reaction mechanisms and to enhance the design of more effective catalysts.

Ahn, Jeonghwan [University of Illinois at Urbana-C↗
DOE OSTI2021

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Comparison of Exact and Model Exchange–Correlation Potentials for Molecules

Accurate exchange–correlation (XC) potentials for three-dimensional systems–via solution of the inverse density functional theory (DFT) problem–are now available to test the quality of DFT approximations. Herein, the exact XC potential for seven molecules–dihydrogen at four different bond-lengths, lithium hydride, water, and ortho-benzyne–are computed from full configuration interaction reference densities. These are compared to model XC potentials from nonlocal (B3LYP, HSE06, SCAN0, and M08-HX) and semilocal/local (SCAN, PBE, and PW92) XC functionals. Whereas for most systems, relative errors in the ground-state densities are O(10 –3 –10 –2 ) , the model XC potentials have much higher errors of O(10 –1 –10 0 ) . Among the model XC functionals, SCAN0 offers the best agreement with the exact XC potential, underlining the significance of satisfying exact conditions as well as including nonlocal effects in XC functionals. This work indicates that tests against the exact XC potential will provide a promising new direction for building more accurate XC functionals for DFT.

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

Spectroscopic Signature of Metal-hydroxo and Peroxo Species in K-edge X-ray Absorption Spectra

Metal-dioxygen species are important intermediates formed during dioxygen activations by metalloenzymes in various biological processes, by catalysts in fuel cells, and prior to O 2 evolution by photosystem II. In this work, we focus on manganese-porphyrin complexes using tetramesitylporphyrin ligand (TMP) to explore changes in Mn K-edge X-ray absorption spectroscopy (XAS) associated with the formation of Mn-hydroxide and Mn–O 2 peroxide species. With limited spectroscopic characterization of these compounds, Mn K β X-ray emission spectroscopy (XES), XAS, density functional theory (DFT), and time-dependent DFT (TD-DFT) analysis will enhance our understanding of their complex electronic structure. We show that the shape of the pre-edge in the K-edge Mn X-ray absorption near-edge structure (XANES) can serve as a spectroscopic signature of the Mn III -peroxo formation and thus can be used to track the presence of the side-on peroxide as an intermediate in time-resolved or in situ experiments. In conclusion, our results will help to further summarize the spectroscopic fingerprints for peroxo and hydroxo species, addressing the challenge of identifying the reactive metal species in catalytic reactions.

Extended X-ray absorption fine structure↗
DOE OSTI2025

Resolution of Selectivity Steps of CO Reduction Reaction on Copper by Quantum Monte Carlo

Electrochemical reduction of carbon monoxide to valuable fuels and chemicals on copper surfaces remains a challenging area in catalysis due to a limited understanding of adsorption mechanisms and reaction pathways. Although density functional theory (DFT)-based studies have investigated these processes, their accuracy varies across different functionals. Here, in this study, we present the application of fixed-node diffusion Monte Carlo (FNDMC) to benchmark the adsorption energies of CO*, H*, and key CO reduction reaction (CORR) intermediates, COH* and CHO* on the Cu(111) surface. Our results for CO* and H* adsorption energies closely align with experimentally measured chemisorption reactions, highlighting the limitations of DFT and providing site-specific energy comparisons that are often not available experimentally. Additionally, we explore the effect of explicit solvation, demonstrating how water stabilizes the COH* over CHO*, thus suggesting a critical role of COH* in CORR. Finally, we release our high-accuracy FNDMC benchmarks for testing and developing new DFT functionals for electrocatalysis. Overall, this study underscores the potential of FNDMC for detailed surface chemistry studies and offers new insights into catalytic processes.

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

Conjugation in Isomeric Cyclosilane Thioethers

Conjugation in thioether-functionalized cyclosilanes depends on both conformation and connectivity. Borane-catalyzed dehydrocoupling of thiophenols with partially hydrogenated cyclosilanes proceeded with high chemoselectivity for Si–H bonds in the presence of Si–Si bonds. Isomeric 1,3- and 1,4-functionalized cyclosilanes afforded isomeric cyclosilane thioethers that also differed in preference for chair or twist-boat conformations as demonstrated by both X-ray crystal structures and density functional theory (DFT) calculations. Here, the confounding effects of different conformations and connectivity could be distinguished with time-dependent DFT (TD-DFT) calculations and validated with experimentally obtained UV–vis spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Discharging Behavior of Hollandite α-MnO 2 in a Hydrated Zinc-Ion Battery

Hollandite, α-MnO 2 , is of interest as a prospective cathode material for hydrated zinc ion batteries (ZIBs); however, the mechanistic understanding of the discharge process remains limited. Herein, a systematic study on the initial discharge of α-MnO 2 cathode under hydrated environment was reported using density functional theory (DFT) in combination with complementary experiments, where the DFT predictions well described the experimental measurements on discharge voltages and manganese oxidation states. According to the DFT calculations, both protons (H + ) and zinc ions (Zn 2+ ) contribute to the discharging potentials of α-MnO 2 observed experimentally, where the presence of water plays an essential role during the process. This study provides valuable insights into the mechanistic understanding of the discharge of α-MnO 2 in hydrated ZIBs, emphasizing the crucial interplay among the H 2 O molecules, the intercalated Zn 2+ or H + ions, and the Mn 4+ ions on the tunnel wall to enhance the stability of discharged states and, thus, the electrochemical performances in hydrated ZIBs.

25 ENERGY STORAGE↗
DOE OSTI2025

Density Functional Tight Binding Insights into Plasmonic Silver–Platinum Nanoparticles and Alloys for Enhanced Photocatalysis

Developing accurate and efficient Slater-Koster (SK) tight-binding parameter sets is essential for quantum plasmonic studies of alloyed metal nanoparticles, as conventional time dependent density functional theory (TD-DFT) calculations are computationally prohibitive for larger clusters. In this work, we develop and validate density functional tight binding (DFTB) parameter sets for both ground state (GS-SK) and excited state (ES-SK) calculations to study the structural, electronic, and optical properties of silver (Ag), platinum (Pt), and Ag–Pt nanoalloys. Our investigation of the ground state properties demonstrates that the GS-SK parameters enable DFTB to closely reproduce the electronic structures of platinum clusters with diverse sizes and geometries – showing qualitative agreement with DFT for density of states (DOS) profiles and energy levels. The ES-SK parameters accurately describe excited-state properties compared to TD-DFT reference calculations, including the broad, featureless absorption profiles of Pt that are dominated by interband transitions. Using the ES-SK parameters within a real-time TD-DFTB framework, we compute size-dependent optical absorption spectra of Ag, Pt and Ag-Pt nanocubes containing up to 1099 atoms (size ∼4.18 nm). A detailed study of Ag–Pt and Pt-Ag core–shell nanoparticles shows quenching of the Ag plasmon resonance even at monolayer coverage for Ag-Pt, but not for Pt-Ag. We also show how to define submonolayer Ag-core Pt-shell cubic structures that have similar optical properties to those generated experimentally for much larger particles, which offers potential for describing plasmon-enhanced photocatalysis. Collectively, the GS-SK and ES-SK parameter sets provide an accurate, computationally efficient approach for modeling the complex optical and electronic behavior of noble–transition metal nanostructures and their alloys.

SPR↗
DOE OSTI2023

Mechanisms of Ethylene Epoxidation over Silver from Machine Learning-Accelerated First-Principles Modeling and Microkinetic Simulations

In this work we employed machine learning-augmented density functional theory (DFT) thermodynamic calculations to assess the stability of different AgO x structures under catalytic ethylene epoxidation reaction conditions. We found that there are multiple AgO x surface motifs that could co-exist under the relevant conditions. These included Ag surface oxides (e.g., AgO_p(4 × 4) and Ag 1.83 O) and atomic oxygen-covered Ag(111) surfaces. Furthermore, we employed DFT calculations to evaluate the energetics of different reaction mechanisms by which ethylene and oxygen can react on these surfaces. These studies revealed several energetically viable reaction pathways for ethylene epoxidation. Microkinetic modeling analyses, based on the DFT-calculated reaction pathways, showed that ethylene epoxidation can proceed on all surfaces and that multiple pathways, including those involving Langmuir–Hinshelwood and Eley–Rideal mechanisms, could be involved in selective and unselective reactions. The diversity of mechanisms that we discovered in the context of the relatively simple ethylene epoxidation reaction on Ag suggests that the richness and complexity of surface chemistry are most likely a rule rather than an exception in heterogeneous catalytic chemical transformations on metal surfaces and that the concept of a single or even a dominant mechanism and reaction intermediates might need to be revisited for many reactions.

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

Rapid Polyolefin Plastic Hydrogenolysis Mediated by Single-Site Heterogeneous Electrophilic/Cationic Organo-group IV Catalysts

A homologous series of cationic electrophilic group IV metal hydrides (M = Ti, Zr, Hf) created by chemisorption of the corresponding MNp 4 precursors on highly Brønsted acidic sulfated alumina (AlS) to yield loosely coordinated surface AlS/MNp 2 (Np = neopentyl) species are systematically characterized by ICP, EXAFS/XANES, DRIFTS, and solid-state high-resolution multidimensional NMR spectroscopy (SSNMR), as well as by energy span DFT computation. With effective stirring, these complexes readily undergo reaction with H 2 to yield AlS/M(alkyl)H species, which are highly active for the hydrogenolysis of diverse commercial polyethylenes, α-olefin-ethylene copolymers, isotactic polypropylene, and postconsumer polyolefins including high-density polyethylenes, yielding medium and small linear and branched hydrocarbons at turnover frequencies as high as 36,300 h –1 at 200 °C/17 atm H 2 for M = Zr. For given polyolefin and reaction conditions, turnover frequencies scale approximately as M = Zr > Hf > Ti, whereas catalyst thermal stability scales approximately as M = Hf ≈ Zr > Ti, and these trends are qualitatively understandable from the DFT analysis. These catalytic results reveal that the AlS/Hf(R)H-mediated hydrogenolysis favors wax-like and liquid products, whereas the AlS/Zr(R)H-mediated hydrogenolysis can be tuned between gases and liquids. Finally, DFT analysis identifies β-alkyl elimination as the turnover-limiting C–C scission process, which is particularly facile in these cationic d 0 complexes but not so in the neutrally charged analogues.

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

Selective Hydrogenation of Furfural Acetone over a Cu Catalyst: Combined Theoretical and Experimental Study

The selective hydrogenation of biomass-derived hydroxymethyl furfural (HMF)-acetone-HMF (HAH) presents an alternative route to producing highly functionalized polyesters and polyurethanes. While HAH undergoes furan ring hydrogenation over Pd, Ru, and Ni catalysts, furan ring hydrogenation is not observed over Cu. Herein, we combine reaction kinetics experiments and density functional theory calculations to elucidate the selective hydrogenation behavior of HAH over Cu catalysts. We identified furfural acetone (FAc) as a suitable surrogate for modeling HAH hydrogenation over Cu surfaces and performed reaction kinetics experiments between temperatures of 313–393 K and a H 2 partial pressure of 55 bar. Similar to the behavior of HAH, hydrogenation of FAc follows a consecutive two-step hydrogenation pathway over Cu and does not undergo furan ring hydrogenation. The apparent activation energy barriers for hydrogenation of the aliphatic double bond (0.58 eV) and carbonyl (0.43 eV) of FAc measured in a continuous flow reactor setup are consistent with those reported in prior studies for batch HAH hydrogenation. Reaction orders with respect to each reactant, including H 2 and FAc, were determined to be nearly one. Density functional theory (DFT; GGA-PBE-D3) calculations on Cu(111) showed that the hydrogenation of the aliphatic double bond of FAc is more facile than the hydrogenation of the furan ring, which displays weak interactions with the Cu surface. We determined an apparent activation energy barrier for FAc hydrogenation (0.58 eV) that agreed with our DFT predictions (highest barrier for FAc hydrogenation of 0.57 eV). Finally, our DFT calculations further show that weak interactions between the furan ring and Cu surface are responsible for the selective hydrogenation behavior.

Biomass↗