Engineering topics

Heyden, Andreas

Publications and source records attributed to Heyden, Andreas.

11 records

DOE OSTI2026

Hydrogen Adsorption over Transition Metals in Water

The adsorption free energy of atomic hydrogen on Pt(111), Pd(111), Ni(111), Ru(0001), Cu(111), and Rh(111) in liquid water was computed using a quantum mechanical/molecular mechanical free-energy perturbation scheme. The Pt(111) computations indicate that the solvent effect on H adsorption on atop sites (+0.20 eV) is almost twice that on fcc (+0.12 eV), showing it is less likely to find adsorbed hydrogen in atop position in the presence of water than in the gas phase. The solvent effect for the fcc site, which is the most favorable site for adsorbed H on Pt(111), agrees qualitatively with experimental work by Lercher et al., who reported an effect of +0.2 eV. Overall, an endergonic solvent effect for hydrogen adsorption is observed for all metals, indicating a lower hydrogen coverage relative to free site coverage at metal-water interfaces compared to metal-gas interfaces, even when hydrogen transport effects through the fluid phase are negligible; a result with important implications for (de)hydrogenation catalysis.

Zare, Mehdi↗

DOE OSTI2025

Characterizing the surface compositions of supported bimetallic PtSn clusters: Effects of cluster-support interactions and surface adsorbates

PtSn bimetallic clusters on TiO2(110) and highly oriented pyrolytic graphite (HOPG) surfaces have been characterized by scanning tunneling microscopy, low energy ion scattering (LEIS), Xray photoelectron spectroscopy, and temperature programmed desorption (TPD); density functional theory (DFT) calculations have also been performed to better understand adsorption of CO and D2 on the PtSn surfaces. On TiO2 at coverages of 2 ML of Pt and 2 ML of Sn, exclusively bimetallic clusters are formed for both orders of deposition because clusters of the first metal completely cover the surface such that all atoms of the second metal are incorporated into the existing clusters. In contrast, on HOPG, the high mobility and weak cluster-support interactions on HOPG result in much larger 2 ML monometallic clusters (~30 Å high) that do not completely cover the surface, and deposition of the second metal produces larger clusters as well as smaller ones. Despite the difference in cluster morphologies for the different orders of deposition and supports, the LEIS experiments demonstrate that in all cases, the PtSn clusters are rich in Sn at the surface, as expected based on the lower surface free energy for Sn compared to Pt. Furthermore, the +0.2 eV shift in the Sn(3d5/2) binding energy observed on all surfaces in the presence of Pt is consistent with PtSn alloy formation. Deposition of 2 ML of Sn on TiO2 produces two-dimensional clusters with oxidation of Sn and reduction of titania at the clustersupport interface, but addition of Pt to the Sn clusters causes Sn to diffuse away from this interface, leaving Sn in the metallic state. TPD experiments on 2 ML Pt/TiO2 with increasing coverages of Sn show that the number of adsorption sites for D2 sharply decreases to nearly zero at 0.5 ML, while CO adsorption decreases to zero only at much higher Sn coverages of 2 ML. DFT studies for Sn modified Pt surfaces and bulk structures demonstrate that for CO adsorption at low Sn coverages (<0.25 ML), the strong Pt-CO interactions induce diffusion of Pt to the cluster surface and the formation of a bulk Pt3Sn alloy, whereas D2 adsorption does not lead to interactions with the Pt surface that are strong enough to induce alloy formation. A single Sn adatom prevents D2 adsorption on four neighboring Pt atoms via site-blocking and the donation of electron density to Pt.

Li, Fangliang↗

DOE OSTI2024

Investigation of Ethane Dehydrogenation and Hydrogenolysis on Pt(111), Pt(211), and Pt(100): Bayesian Quantification and Correction of DFT-Based Enthalpic and Entropic Uncertainties

Computational investigations of heterogeneously catalyzed reactions using density functional theory (DFT) are often inaccurate, largely due to uncertainties in the choice of DFT functional (enthalpic uncertainty) and approximations for modeling adsorbate movement along the catalyst surface (entropic uncertainty). This work illustrates that both uncertainties are significant in the investigation of ethane dehydrogenation (EDH) and hydrogenolysis on Pt catalysts by considering the complete deconstruction of ethane on Pt(111), Pt(211), and Pt(100) using microkinetic modeling (MKM). Hence, this work uses both noncalibrated and Bayesian-calibrated MKMs to quantify and correct inaccuracies in macroscopic properties due to both uncertainties. A Bayesian approach to the correction of entropic errors was introduced using a “Modified Fermi Function (MFF)” to calibrate between the two bounds of entropy represented by the harmonic oscillator (HO) and free translator (FT) approximations. Regardless of enthalpic and entropic uncertainties, all three surfaces are capable of ethane activation; however, Pt(211) was found to be the most active and is largely responsible for methane production. Next, Pt(111) is largely responsible for acetylene production, and Pt(100) has the highest ethylene selectivity but is most susceptible to coking. By comparison of different calibrated models, the FT entropy approximation was found to better describe EDH under typical experimental conditions. Statistical evidence was found to support Pt(111) as the active site for EDH, assuming that one single site is responsible for the chemistry. On the three surfaces, competing second dehydrogenations to CH 2 CH 2 and CH 3 CH were observed as well as isomerization of CH 3 CH back to CH 2 CH 2 and deeper dehydrogenation of CH 3 CH. In conclusion, C–C cleavage was found to largely proceed via the CH 3 C intermediate on Pt(100) and Pt(111), while on Pt(211), it was via both CHC and CH 3 C.

Bayesian model selection↗

DOE OSTI2024

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal–water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Here, due to the large system sizes needed for studying processes at liquid water–metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate–metal interactions. Next, for the water–water and water–adsorbate interactions, we used established classical force fields. Finally, for the water–Ru surface interaction, for which no reliable force fields have been published, we used Behler–Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Adsorption↗

DOE OSTI2024

Polyolefin melt-phase effects on alkane hydrogenolysis over Pt catalysts

Supported transition metal catalyzed, chemical upcycling of polyolefins by hydrogenolysis typically occurs in a polymer melt phase at elevated temperatures (T > 200 ºC). Currently, the impact of the melt phase on the catalytic activity and selectivity of the transition metal is largely unknown. Herein, we use a hybrid quantum mechanical/molecular mechanical (QM/MM) approach to investigate the melt-phase effects on the adsorption free energy (∆∆G$^{gas→liq}_{Adsorbate}$) of atomic hydrogen, 12 hydrocarbon molecules, and 4 transition states in the hydrogenolysis mechanism of butane on a Pt(111) catalyst surface at 573 K in the presence of a polyethylene surrogate melt consisting of C 36 H 74 chains. The smallest and largest endergonic melt phase effects, (∆∆G$^{gas→liq}_{Adsorbate}$), belong to hydrogen (0.045 eV) and butane (1.357 eV). Altogether, we find melt-phase effects are significant and change the activity of transition metal catalysts. Beyond an overall reduced adsorption strength, elementary surface reactions are also affected by the melt phase.

Pt catalysis↗

DOE OSTI2024

Machine Learning Accelerated First-Principles Study of the Hydrodeoxygenation of Propanoic Acid

The complex reaction network of catalytic biomass conversions often involves hundreds of surface intermediates and thousands of reaction steps, greatly hindering the rational design of metal catalysts for these conversions. Here, we present a framework of machine learning (ML)-accelerated first-principles studies for the hydrodeoxygenation (HDO) of propanoic acid over transition metal surfaces. The microkinetic model (MKM) is initially parametrized by ML-predicted energies and iteratively improved by identifying the rate-determining species and steps (RDS), computing their energies by density functional theory (DFT), and reparameterizing the MKM until all the RDS are computed by DFT. The Gaussian process (GP) model performs significantly better than the linear ridge regression model for predicting both the adsorption free energies and transition state free energies. Parameterized with energies from the GP model, only 5–20% of the full reaction network has to be computed by DFT for the MKM to possess DFT-level accuracy for the TOF and dominant reaction pathway. While the linear ridge regression model performs worse than the GP model, its performance is greatly improved when only transition states are predicted by the regression model and adsorption energies are computed by DFT. Overall, we find that a high accuracy in adsorption free energies is more important for a reliable MKM than a high accuracy in TS free energies. Lastly, based on the GP model with GOH and GCHCHCO as catalyst descriptors, we build two-dimensional volcano plots in activity and selectivity that can help design promising alloy catalysts for HDO reactions of organic acids.

adsorption↗

DOE OSTI2024

A redox-reversible A/B-site co-doped BaFeO 3 electrode for direct hydrocarbon solid oxide fuel cells

Solid oxide fuel cells (SOFCs) can directly convert the chemical energy in fuel to electrical energy with fuel flexibility; however, the conventional nickel-based anodes face great challenges due to coking upon direct oxidation of hydrocarbon fuels and redox instability. Thus, developing new anode materials which can provide high coking resistance as well as redox stability is crucial. In this work, Ba 0.6 La 0.4 Fe 0.8 Mo 0.1 Ni 0.1 O 3-δ (BLFMN) has been synthesized in air using a sol–gel combustion method, resulting in a dual phase consisting of a cubic BLFMN main phase and scheelite BaMoO 4 (BMO 4 ) secondary phase. By heat-treating the BLFMN dual phase in H 2 at 800 °C for 5 h, a metallic nanoparticle-decorated BLFMN triple phase compound comprising cubic BLFMN, cubic BaMoO 3 (BMO 3 ) and in situ exsolved FeNi 3 alloy was obtained. BLFMN was subsequently investigated as an electrode material for La 0.8 Sr 0.2 Ga 0.83 Mg 0.17 O 3-δ (LSGM) electrolyte (300 μm) supported SOFCs. Symmetrical cells using BLFMN as electrodes with the cell configuration of BLFMN//LSGM//BLFMN showed excellent redox reversibility and a peak power density (PPD) of 1.32 W cm -2 at 850 °C when using H 2 as fuel. Single cell with the cell configuration of BLFMN//LSGM//LSCF (La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ ) reached PPD of 1.61 and 0.41 W cm -2 at 850 °C when operating with H 2 and CH 4 fuel, respectively. Moreover, the single cell exhibit excellent stability (over 300 h) upon direct oxidation of hydrocarbon fuels of CH 4 and C 3 H 8 . This study indicates that BLFMN is a promising redox reversible and coking resistant anode for SOFCs.

08 HYDROGEN↗

DOE OSTI2024

Invariant Molecular Representations for Heterogeneous Catalysis

Catalyst screening is a critical step in the discovery and development of heterogeneous catalysts, which are vital for a wide range of chemical processes. In recent years, computational catalyst screening, primarily through density functional theory (DFT), has gained significant attention as a method for identifying promising catalysts. However, the computation of adsorption energies for all likely chemical intermediates present in complex surface chemistries is computationally intensive and costly due to the expensive nature of these calculations and the intrinsic idiosyncrasies of the methods or data sets used. This study introduces a novel machine learning (ML) method to learn adsorption energies from multiple DFT functionals by using invariant molecular representations (IMRs). To do this, we first extract molecular fingerprints for the reaction intermediates and later use a Siamese-neural-network-based training strategy to learn invariant molecular representations or the IMR across all available functionals. Our Siamese network-based representations demonstrate superior performance in predicting adsorption energies compared with other molecular representations. Notably, when considering mean absolute values of adsorption energies as 0.43 eV (PBE-D3), 0.46 eV (BEEF-vdW), 0.81 eV (RPBE), and 0.37 eV (scan+rVV10), our IMR method has achieved the lowest mean absolute errors (MAEs) of 0.18 0.10, 0.16, and 0.18 eV, respectively. These results emphasize the superior predictive capacity of our Siamese network-based representations. The empirical findings in this study illuminate the efficacy, robustness, and dependability of our proposed ML paradigm in predicting adsorption energies, specifically for propane dehydrogenation on a platinum catalyst surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Methane Partial Oxidation over Multifunctional 2-D Materials

The objective of this research is to design, synthesize, and evaluate highly selective, active, and stable multifunctional catalysts for the low temperature (< 500 Kelvin (K)) partial oxidation of methane to methanol (MTM) with molecular oxygen: CH 4 (g) + $\frac{1}{2}$O 2 (g) → CH 3 OH(g). Methane, the primary component of natural gas, is a source of energy and economic growth as well as an environmental concern. Recent developments in horizontal drilling and enhanced extraction methods have resulted in production of an estimated 62.4 trillion m 3 of ‘stranded’, or uneconomic, natural gas. Uneconomical natural gas is often flared or vented at remote oil production sites. Leaked, flared, and/or vented gas represents a "lost opportunity”, and this research project aims to maximize the value of the resource. Conventional catalysts for MTM suffer from low methanol selectivity since they exhibit ~0.55 eV higher barrier for C-H bond activation of methane compared to methanol. Without breaking these scaling relations, methanol oxidation is orders of magnitude faster than methane oxidation and it is very challenging to envision a process with economically viable single-pass yield. Here, we chose to investigate single-atom catalysts embedded and stabilized in two-dimensional materials such as graphene (GR) and "supported" on Group VIII and IB transition metals such as nickel. The electronic atomic monolayer-metal support interaction (EAMSI) present in these systems could promote methanol selectivity by breaking the scaling relations of the C-H bond activation of methane and methanol. A density functional theory (DFT) based computational study focused on predicting families of GR-based catalysts that could be active and selective for MTM. The catalyst systems predicted by the computational study were synthesized and evaluated for the gas phase MTM under relevant conditions. Unfortunately, the experimental activity and selectivity was lower than computationally predicted. The origin for the discrepancy is likely related to difficulties in synthesizing single atom catalysts in a threecomponent catalyst system at high density and with high selectivity. Future work in our groups is thus focused on reducing the system complexity to a two-component catalyst system. Finally, a techno-economic analysis (TEA) was also conducted to identify critical bottlenecks that inhibit future commercialization.

03 NATURAL GAS↗

DOE OSTI2023

Unlocking the Potential of A-Site Ca-Doped LaCo 0.2 Fe 0.8 O 3-δ : A Redox-Stable Cathode Material Enabling High Current Density in Direct CO 2 Electrolysis

Massive carbon dioxide (CO 2 ) emission from recent human industrialization has affected the global ecosystem and raised great concern for environmental sustainability. The solid oxide electrolysis cell (SOEC) is a promising energy conversion device capable of efficiently converting CO 2 into valuable chemicals using renewable energy sources. However, Sr-containing cathode materials face the challenge of Sr carbonation during CO 2 electrolysis, which greatly affects the energy conversion efficiency and long-term stability. Thus, A-site Ca-doped La1– x CaxCo 0.2 Fe 0.8 O 3-δ (0.2 ≤ x ≤ 0.6) oxides are developed for direct CO 2 conversion to carbon monoxide (CO) in an intermediate-temperature SOEC (IT-SOEC). With a polarization resistance as low as 0.18 O cm 2 in pure CO 2 atmosphere, a remarkable current density of 2.24 A cm –2 was achieved at 1.5 V with La 0.6 Ca 0.4 Co 0.2 Fe 0.8 O 3-δ (LCCF64) as the cathode in La 0.8 Sr 0.2 Ga 0.83 Mg 0.17 O 3-δ (LSGM) electrolyte (300 µm) supported electrolysis cells using La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ (LSCF) as the air electrode at 800 °C. Furthermore, symmetrical cells with LCCF64 as the electrodes also show promising electrolysis performance of 1.78 A cm –2 at 1.5 V at 800 °C. In addition, stable cell performance has been achieved on direct CO 2 electrolysis at an applied constant current of 0.5 A cm –2 at 800 °C. The easily removable carbonate intermediate produced during direct CO 2 electrolysis makes LCCF64 a promising regenerable cathode. The outstanding electrocatalytic performance of the LCCF64 cathode is ascribed to the highly active and stable metal/perovskite interfaces that resulted from the in situ exsolved Co/CoFe nanoparticles and the additional oxygen vacancies originated from the Ca 2 Fe 2 O 5 phase synergistically providing active sites for CO 2 adsorption and electrolysis. Here this study offers a novel approach to design catalysts with high performance for direct CO 2 electrolysis.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2023

Modeling the Effect of Surface Platinum–Tin Alloys on Propane Dehydrogenation on Platinum–Tin Catalysts

Uncertainty analysis, reported experimental literature data, and density functional theory were synthesized to model the effect of surface tin coverage on platinum-based catalysts for nonoxidative propane dehydrogenation to propylene. Here, this study tests four different platinum–tin skin surface models as potential catalytic sites, Pt 3 Sn/Pt(100), PtSn/Pt(100), Pt 3 Sn/Pt(111), and Pt 2 Sn/Pt(211), and compares them to the corresponding pure Pt surface sites using an uncertainty analysis methodology that uses BEEF-vdW with its ensembles (BMwE) to generate the uncertainty for the energies of the intermediates and transition states. One experimental data set with two experimental observations, selectivity to propylene and turnover frequency of propylene, was used as a calibration data set to evaluate the impact of the experimental data on informing the models. This study finds that the prior model for Pt 3 Sn/Pt(100) is the most active and Pt 2 Sn/Pt(211) is the most selective toward propylene. Active sites on the (100) facet have the highest probability of being responsible for C 1 and C 2 product formations (C–C bond cleavage). Increasing the Sn coverage on the (100) surface facet to a PtSn/Pt(100) active site leads to a significantly reduced rate and might explain the experimentally observed higher selectivity of Sn-doped catalysts relative to pure Pt catalysts. Next, this study finds that for all surfaces, except PtSn/Pt(100), the rate-controlling steps are the initial dehydrogenation steps alongside some partially rate-controlling second dehydrogenation steps. For PtSn/Pt(100), only the initial terminal dehydrogenation step to CH 3 CH 2 CH 2 * and second dehydrogenation steps are rate-controlling. Next, the calibrated models for all surfaces were found to be selective toward propylene production and model the reported turnover frequency successfully. Nevertheless, Pt 2 Sn/Pt(211) emerges as the active site with some (minor) evidence as the main active site based on Jeffreys’ scale interpretation of Bayes factors. This observation agrees with prior studies that also found step sites to be most likely the most relevant active sites for pure Pt catalysts. Overall, the results indicate that tin, in addition to affecting the binding strength of the adsorbed species, prevents deeper dehydrogenation (reducing coking) and cracking reactions through increasing activation barriers for unwanted side reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗