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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Assessment of exchange-correlation functionals on oxygen vacancy formation energies of metal oxides

Oxygen vacancies are a common defect in oxides and play vital roles in many technological applications such as oxygen separation, catalytic reactors, solid oxide fuel cells, and solar thermochemical water splitting. The oxygen vacancy formation energy is directly related to the oxide reduction enthalpy and the vacancy concentration and, therefore, is a key quantity underlying these applications. The complexity of measuring oxygen vacancy formation energies experimentally suggests the utility of calculations based on density functional theory (DFT). However, the calculated results are strongly dependent on the exchange-correlation functionals and other parameters such as the Hubbard U . In this work, we compare the performance of the meta-GGA strongly constrained and appropriately normed (SCAN) functional and the commonly used semilocal generalized gradient approximation (GGA) functionals with experimental values of structural parameters, band gaps, magnetic structures, and oxygen vacancy formation energies for six representative oxides, i.e., CaMnO 3 , SrMnO 3 , LaMnO 3 , YMnO 3 , LaFeO 3 , and CeO 2 . Further, our results show that SCAN usually has better agreement with the experimental lattice constants and band gaps and larger magnetic moments than the commonly used GGA functionals. Although SCAN overestimates the oxygen vacancy formation energies of transition metal oxides and therefore requires unusually large Hubbard U values to reproduce the experimental reduction enthalpies, it does predict the correct oxygen vacancy formation energy of CeO 2 , which challenges the commonly used GGA and hybrid functionals. Our results underscore the challenges that exist in describing these complex oxides by DFT and may shed light on developing more accurate exchange-correlation functionals for oxygen vacancy formation energy calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Selective activation of methane C–H bond in the presence of methanol

Direct methane to methanol (MTM) conversion over heterogeneous catalysts is a promising route for valorization of methane. The methane C-H bond activation is considered as the key step of the MTM and is the focus of considerable research activity. However, the formed methanol typically suffers from overoxidation largely due to the cleavage of a methanol C-H bond, whose bond dissociation energy is ca. 0.5 eV lower than that of the methane C-H bond, which usually translates to a transition state energy of the methanol C-H bond cleavage that is ca. 0.55 eV lower than that of methane whenever the reactions proceed through a radical mechanism. Furthermore, we propose a general approach for decreasing the transition state energy difference between the CH 4 and CH 3 OH C-H bond dissociation. When a metal-oxide supported cationic transition metal atom and a neighboring oxygen on the oxide surface serve as the active site, the transition state energy difference through a surface-stabilized pathway can be noticeably narrowed as compared with that of a radical pathway. For Ir, Pt and Rh-doped anatase TiO 2 (101), the CH 4 C-H bond activation can be preferred over that of CH 3 OH at significant methanol mole fraction. Also, for PdAu alloys containing adsorbed oxygen and positively charged Pd, calculations suggest in agreement with recent experiments (Science 367 (6474), 193-197) that the CH 4 C-H bond can be selectively activated in the presence of CH 3 OH.

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

Phase stability and coercivity in La 2 Fe 14 B magnet

Critical rare-earth free La 2 Fe 14 B (2:14:1) has the potential to be a gap permanent magnet. However, La 2 Fe 14 B decomposes into La, α-Fe, and La 4 Fe 4 phases below 1067 K. The phase stability and coercivity have been studied in La 2 Fe 14 B magnet using first principles DFT (density functional theory) calculation and micromagnetic simulation. For a perfect La 2 Fe 14 B cube (edge length of 256 nm) without any structural defects and soft magnetic secondary phases, the coercivity (8.5 kOe) is reduced to ~40% of its magnetocrystalline anisotropy field ( H A = 20 kOe). Further, the coercivity sharply reduces to 3.2 kOe upon forming a thin layer (2 nm) of α-Fe on the surface of the La 2 Fe 14 B cube particle. The DFT calculations indicate that a partial replacement of La by other rare-earth (RE) elements can enhance the structural stability of 2:14:1. The gains in cohesive energy are 0.75, 0.10, and 0.33 eV per formula unit in (La 0.5 RE 0.5 ) 2 Fe 14 B with RE = Ce, Pr, and Nd, respectively. Stabilizing the 2:14:1 structure and mitigating the formation of soft magnetic structural defects or impurity phases such as α-Fe is necessary to develop La 2 Fe 14 B based magnet, which can be moderately achieved via partial substitution of La by other rare earth elements such as Ce, Pr, and Nd.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Salt-flux synthesis, crystal structure and theoretical characterization of Rb 0.74 Ga 6.62 Ti 0·38 O 11

Here single crystals of Rb 0.74 Ga 6.62 Ti 0·38 O 11 (RGTO) were grown from a mixed RbCl–RbF flux at 850 °C. The compound crystallizes in the RbGa 7 O 11 structure type, which is reminiscent of the hollandite and β-Ga 2 O 3 structure types. RGTO crystallizes in the monoclinic space group P2/m with lattice parameters a = 8.3355 (8) Å, b = 3.0286 (3) Å, c = 9.5028 (9) Å, and β = 114.620 (3)°. The crystal structure of RGTO is comprised of GaO 6 and mixed (Ga/Ti)O 6 octahedra and GaO 4 tetrahedra connected in a complex three-dimensional, anionic framework exhibiting eight-sided channels that are occupied by disordered Rb cations required for charge balance. First-principles calculations in the form of density functional theory were performed, which indicated the complex to be a charge transfer semiconductor.

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

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗
DOE OSTI2023

Elucidating Biomass-Derived Pyrolytic Lignin Structures from Demethylation Reactions through Density Functional Theory Calculations

Pyrolytic lignin is a fraction of pyrolysis oil that contains a wide range of phenolic compounds that can be used as intermediates to produce fuels and chemicals. However, the characteristics of the raw lignin structure make it difficult to establish a pyrolysis mechanism and determine pyrolytic lignin structures. Herein this study proposes dimer, trimer, and tetramer structures based on their relative thermodynamic stability for a hardwood lignin model in pyrolysis. Different configurations of oligomers were evaluated by varying the positions of the guaiacyl (G) and syringyl (S) units and the bonds βO4 and β5 in the hardwood model lignin through electronic structure calculations. The homolytic cleavage of βO4 bonds is assumed to occur and generate two free radical fragments. These can stabilize by taking hydrogen radicals that may be in solution during the intermediate liquid (pathway 1) formation before the thermal ejection. An alternative pathway (pathway 2) could occur when the radicals use intramolecular hydrogen, turning themselves into stable products. Subsequently, a demethylation reaction can take place, thus generating a methane molecule and new oligomeric lignin-derived molecules. The most probable resulting structures were studied. We used FTIR and NMR spectra of selected model compounds to evaluate our calculation approach. Thermophysical properties were calculated using group contribution methods. The results give insights into the lignin oligomer structures and how these molecules are formed. They also provide helpful information for the design of pyrolysis oil separation and upgrading equipment.

09 BIOMASS FUELS↗
DOE OSTI2020

Thermodynamics of Metal Carbonates and Bicarbonates and Their Hydrates for Mg, Ca, Fe, and Cd Relevant to Mineral Energetics

The heats of formation of carbonate, bicarbonate and bicarbonate/hydroxide metal complexes and hydrated versions these complexes of Mg 2+ , Ca 2+ , Fe 2+ , and Cd 2+ are predicted from atomization energies using correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit following the Feller-Peterson-Dixon (FPD) approach. Using the calculated gas phase values and the available experimental solid state values, the cohesive energies of selected minerals were calculated. The gas phase decomposition energies into MO, CO 2 and H 2 O follow the order Mg ~ Ca > Cd ~ Fe, and correlate with the hardness of the metal +2 ions. Gas phase hydration energies show that the order is Mg > Fe > Ca ~ Cd. There are a number of bulk hydrated Mg and Ca complexes that occur as minerals but there are few if any for Fe and Cd, suggesting that a number of factors are important in determining the stability of the bulk mineral hydrates. The FPD heats of formation were used to benchmark a range of density functional theory exchange-correlation functionals, including those commonly used in solid state mineral calculations. None of the functionals provided chemical accuracy agreement (± 1 kcal/mol) with the FPD results The best functionals are ωB97X and ωB97X-D with mean average unsigned errors of 10 kcal/mol.

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

Thermally stable manganese(Ⅲ) peroxido complexes with hindered N3 tripodal ligands: Structures and their physicochemical properties

Mononuclear manganese(III) peroxido complexes are candidates for the reaction intermediates in manganese containing proteins, such as manganese superoxide dismutase (Mn-SOD) etc. In this study, manganese(III) peroxido complexes [Mn(O 2 )(L3)] and [Mn(O 2 )(L10)] ligated by anionic N3 type ligands with sterically hindered substituents, hydrotris(3-tertiary butyl-5-isopropyl-1-pyrazolyl)borate (L3 - ) and hydrotris(3-adamantyl-5-isopropyl-1-pyrazolyl)borate (L10 - ), respectively, were structurally characterized. These complexes are the first examples of structurally characterized five-coordinate manganese(III) peroxido complexes. Their characteristic ν(O—O) and ν(Mn—O) stretchings were determined by using H 2 18 O 2 for the first time. Theoretical calculations were performed to obtain further insight into their structural parameters. The decomposed products were obtained as [{Mn III (μ-O)(L3)} 2 Mn IV ] and [Mn III (OH){L10(O)}] from [Mn(O 2 )(L3)] and [Mn(O 2 )(L10)], respectively.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

17 O NMR Spectroscopy Reveals CO 2 Speciation and Dynamics in Hydroxide-Based Carbon Capture Materials

Carbon dioxide capture technologies are set to play a vital role in mitigating the current climate crisis. Solid-state 17 O NMR spectroscopy can provide key mechanistic insights that are crucial to effective sorbent development. In this work, we present the fundamental aspects and complexities for the study of hydroxide-based CO 2 capture systems by 17 O NMR. We perform static density functional theory (DFT) NMR calculations to assign peaks for general hydroxide CO 2 capture products, finding that 17 O NMR can readily distinguish bicarbonate, carbonate and water species. However, in application to CO 2 binding in two test case hydroxide-functionalised metal-organic frameworks (MOFs) – MFU-4l and KHCO 3 -cyclodextrin-MOF, we find that a dynamic treatment is necessary to obtain agreement between computational and experimental spectra. We therefore introduce a workflow that leverages machine-learning force fields to capture dynamics across multiple chemical exchange regimes, providing a significant improvement on static DFT predictions. In MFU-4l, we parameterise a two-component dynamic motion of the bicarbonate motif involving a rapid carbonyl seesaw motion and intermediate hydroxyl proton hopping. For KHCO 3 -CD-MOF, we combined experimental and modelling approaches to propose a new mixed carbonate-bicarbonate binding mechanism and thus, we open new avenues for the study and modelling of hydroxide-based CO 2 capture materials by 17 O NMR.

NMR spectroscopy↗
DOE OSTI2022

Metal-Coordinated Phthalocyanines as Platform Molecules for Understanding Isolated Metal Sites in the Electrochemical Reduction of CO 2

Single-atom catalysts (SACs) of non-precious transition metals (TMs) often show unique electrochemical performance, including the electrochemical carbon dioxide reduction reaction (CO 2 RR). However, the inhomogeneity in their structures makes it difficult to directly compare SACs of different TM for their CO 2 RR activity, selectivity, and reaction mechanisms. Here, the comparison of isolated TMs (Fe, Co, Ni, Cu, and Zn) is systematically investigated using a series of crystalline molecular catalysts, namely TM-coordinated phthalocyanines (TMPcs), to directly compare the intrinsic role of the TMs with identical local coordination environments on the CO 2 RR performance. The combined experimental measurements, in situ characterization, and density functional theory calculations of TM-Pc catalysts reveal a TM dependent CO 2 RR activity and selectivity, with the free energy difference of ΔG(*HOCO) - ΔG(*CO) being identified as a descriptor for predicting the CO 2 RR performance.

25 ENERGY STORAGE↗
DOE OSTI2025

Electronic and magnetic properties of hole-doped topological kagome Fe 1−𝑥 ⁢Mn 𝑥 ⁢Sn thin films

We have investigated the electronic and magnetic structures of topological kagome Fe 1−𝑥 ⁢Mn 𝑥 ⁢Sn (0 ≤ 𝑥 ≤ 0.3) thin films via neutron diffraction, electronic transport measurements, and ab initio density functional theory (DFT) to understand the interplay between hole doping, magnetism, and the electronic structures. Temperature-dependent neutron diffraction measurements on parent FeSn reveal the Néel temperature to be 𝑇 N ∼ 355 K and the underlying A-type antiferromagnetic ordering is associated with a wave vector 𝒒 = (001/2). Upon Mn doping to 𝑥 = 0.15, 𝑇 N decreases slightly while the magnetic ordering vector remains the same. Resistivity measurements show metallic characteristics and in-plane anisotropy down to 10 K for all the investigated samples. The effects of hole doping are mapped in terms of electronic ground state calculations via DFT which show that the Dirac point is moved closer to the Fermi level (𝐸 F ) and the flat bands get pushed away from 𝐸 F upon hole doping. However, a comparison between hole-doped Fe 1−𝑥⁢ Mn 𝑥⁢ Sn and electron-doped Fe 1−𝑥 ⁢Co 𝑥 ⁢Sn indicates that the Néel temperature does not scale with the position of 𝐸 F relative to the flat band. Furthermore, our results establish the antiferromagnetic state of FeSn and Fe 1−𝑥 ⁢Mn 𝑥⁢ Sn films at room temperature, laying the groundwork for future studies of magnetism in kagome heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Discovery of Stable Surfaces with Extreme Work Functions by High‐Throughput Density Functional Theory and Machine Learning

Abstract The work function is the key surface property that determines the energy required to extract an electron from the surface of a material. This property is crucial for thermionic energy conversion, band alignment in heterostructures, and electron emission devices. This work presents a high‐throughput workflow using density functional theory (DFT) to calculate the work function and cleavage energy of 33,631 slabs (58,332 work functions) that are created from 3,716 bulk materials. The number of calculated surface properties surpasses the previously largest database by a factor of ≈27. Several surfaces with an ultra‐low (<2 eV) and ultra‐high (>7 eV) work function are identified. Specifically, the (100)‐Ba‐O surface of BaMoO 3 and the (001)‐F surface of Ag 2 F have record‐low (1.25 eV) and record‐high (9.06 eV) steady‐state work functions. Based on this database a physics‐based approach to featurize surfaces is utilized to predict the work function. The random forest model achieves a test mean absolute error (MAE) of 0.09 eV, comparable to the accuracy of DFT. This surrogate model enables rapid predictions of the work function (≈ 10 5 faster than DFT) across a vast chemical space and facilitates the discovery of material surfaces with extreme work functions for energy conversion and electronic device applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Enhancement of van der Waals Interlayer Coupling through Polar Janus MoSSe

Interlayer coupling plays essential roles in the quantum transport, polaritonic, and electrochemical properties of stacked van der Waals (vdW) materials. In this work, we report the unconventional interlayer coupling in vdW heterostructures (HSs) by utilizing an emerging 2D material, Janus transition metal dichalcogenides (TMDs). In contrast to conventional TMDs, monolayer Janus TMDs have two different chalcogen layers sandwiching the transition metal and thus exhibit broken mirror symmetry and an intrinsic vertical dipole moment. Such a broken symmetry is found to strongly enhance the vdW interlayer coupling by as much as 13.2% when forming MoSSe/MoS 2 HS as compared to the pristine MoS 2 counterparts. Our noncontact ultralow-frequency Raman probe, linear chain model, and density functional theory calculations confirm the enhancement and reveal the origins as charge redistribution in Janus MoSSe and reduced interlayer distance. Furthermore, our results uncover the potential of tuning interlayer coupling strength through Janus heterostacking.

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

Tetragonal Kondo Insulator EuCd 2 Sb 2 Discovered via High Pressure High Temperature Synthesis

Abstract Magnetic and electronic properties of quantum materials heavily rely on the crystal structure even in the same chemical compositions. In this study, it is demonstrated that a layered tetragonal EuCd 2 Sb 2 structure can be obtained by treating bulk trigonal EuCd 2 Sb 2 under high pressure (6 GPa) and high temperature (600 °C). Magnetization measurements of the newly formed layered tetragonal EuCd 2 Sb 2 confirm an antiferromagnetic ordering with Neel temperature ( T N ) around 16 K, which is significantly higher than that ( T N ≈ 7 K) of trigonal EuCd 2 Sb 2 , consistent with heat capacity measurements. Moreover, bad metal behavior is observed in the temperature dependence of the electrical resistivity and the resistivity shows a dramatic increase around the Neel temperature. Electronic structure calculations with local density approximation dynamic mean–field theory (LDA+DMFT) show that this material is strongly correlated with well‐formed large magnetic moments, due to Hund's coupling, which is known to dramatically suppress the Kondo scale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Computational and Experimental 1H-NMR Study of Hydrated Mg-Based Minerals

Magnesium oxide (MgO) can convert to different magnesium-containing compounds depending on exposure and environmental conditions. Many MgO-based phases contain hydrated species allowing 1H-nuclear magnetic resonance (NMR) spectroscopy to be used in the characterization and quantification of proton-containing phases; however, surprisingly limited examples have been reported. Here, 1H-magic angle spinning (MAS) NMR spectra of select Mg-based minerals are presented and assigned. These experimental results are combined with computational NMR density functional theory (DFT) periodic calculations to calibrate the predicted chemical shielding results. This correlation is then used to predict the NMR shielding for a series of different MgO hydroxide, magnesium chloride hydrate, magnesium perchlorate, and magnesium cement compounds to aid in the future assignment of 1H-NMR spectra for complex Mg phases.

1H-NMR↗
DOE OSTI2020

Probing Conformational Evolution and Associated Dynamics of Mg(N(SO 2 CF 3 ) 2 ) 2 ·Dimethoxyethane Adduct Using Solid-State 19 F and 1 H NMR

Bis(trifluoromethanesulfonimide) or TFSI is widely used as a counter anion in electrolyte design due to its structural flexibility and chemical stability. Here we studied the conformational variations and associated dynamics of TFSI in adduct of Mg(TFSI) 2 with dimethoxyethane (DME), a solvate crystalline material using solid-state 1 H and 19 F NMR. TFSI molecular motion in this solvate structure falls within the timescale of the 19 F NMR experiment, yielding spectroscopic signatures for unique TFSI conformers under the coordination environment of Mg 2+ cation. Within the temperature range of -5 to 82 °C, we observe nine distinct TFSI sites in both crystalline and disordered regions using 19 F NMR, reflecting complexity of structural and dynamics of TFS1 anions within solvate structure. The four distinguishable sites in the disordered region for the two CF 3 groups of the same TFSI molecule are identified using chemical shift analysis. The exchange rate constants from site to site are calculated through variable temperature 19 F NMR and two-dimensional (2D) exchange spectroscopy (EXSY) experiments, along with respective activation enthalpies using Eyring's formulation. The flip rate of CF 3 around the S-C bond is estimated as ~ 15 s -1 at 8 °C with ΔH ≠ ~ 22 kJ/mol, but the rotation of the entire TFSI is 4.8 s -1 at 8 °C with a significantly greater ΔH ≠ = 98 ± 10 kJ/mol. Furthermore, the slow conversion of trans to cis conformers at a lower temperature (T ≤ 1 °C) in the crystalline region is monitored, with a conversion rate of ~ 2 x 10 -5 s -1 at -5 °C. Density functional theory (DFT)-based calculations were performed to support further the assignment of experimental chemical shifts, and the activation energy E a = 21.1 kJ/mol obtained for the cis to trans transition is consistent with experimental values. The combined set of 19 F and 1 H under both one-dimensional (1D) and 2D NMR methods demonstrated here can be further used for examining electrode-electrolyte interfaces to probe the motions of various constituents that can enable detailed studies of interfacial processes and dynamics. Ultimately, such studies will aid in the design and discovery of interfacial constructs in which directed defect chemistry, chemical moiety distribution, and nanostructure are employed to drive efficient charge transport.

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

Kinetics and Reaction Mechanisms of Acetic Acid Hydrodeoxygenation over Pt and Pt–Mo Catalysts

In this work, kinetic measurements for silica-supported Pt and Pt-Mo catalysts were collected in vapor-phase acetic acid hydrodeoxygenation by varying the hydrogen partial pressure between 18 and 72 kPa and acetic acid partial pressure between 7 and 18 kPa at 423-473 K. At all testing conditions, the Pt-Mo catalyst was more active and selective. In addition, the apparent activation energy for Pt-Mo of 76±3 kJ/mol was lower than that of 84±4 kJ/mol for Pt. The apparent reaction orders were also different. The order in hydrogen of 0.8±0.1 for Pt-Mo changed to zero at higher hydrogen partial pressures while that for Pt remained constant at 0.6±0.1. A near-zero order in acetic acid for Pt-Mo changed to -2.1 at higher acetic acid pressures while that for Pt remained constant at -2.9±0.3. These differences in reaction kinetics as well as in selectivity trends with changes in temperature and feed composition indicated a change in the reaction mechanism for Pt-Mo. The catalysts were characterized with hydrogen temperature programmed desorption, oxygen temperature programmed oxidation and transmission electron microscopy with energy-dispersed X-ray spectroscopy elemental mapping. Mo was present in the form of subnanometer-size clusters on the surface of Pt nanoparticles. Both Pt and Pt-Mo catalysts were stable under the reaction conditions for 10 h, and the size and structure of Pt and Pt-Mo particles remained mostly unchanged, without coke accumulation. Density functional theory calculations show that surface acetate is not a major reaction intermediate on both Pt and Pt-Mo and, instead, C-OH bond splitting with the formation of acetyl is the first reaction step in hydrodeoxygenation. The activation energy for this step is dramatically lower on Pt-Mo. Furthermore, the activity of acetyl on Pt-Mo is different. Unlike on pure Pt, the reaction of C-O bond splitting becomes exothermic with a lower activation energy on Pt-Mo, and the reaction of C-C bond splitting, in contrast, becomes endothermic with a higher activation energy, explaining the experimentally observed higher activities and selectivities of Pt-Mo. In addition, the calculations demonstrate that a pair of Pt-Mo surface atoms acts as a single active site where Mo serves as a preferential binding anchor for O atoms. The presence of Mo atoms changes the structure and energy of adsorbed and reacting surface species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗