Engineering topics

Heyden, Andreas

Publications and source records attributed to Heyden, Andreas.

7 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↗