Engineering topics

Mavrikakis, Manos

Publications and source records attributed to Mavrikakis, Manos.

At least 19 records

DOE OSTI2025

Operando probing dynamic migration of copper carbonyl during electrocatalytic CO2 reduction

Single crystals and shape-controlled nanocrystals are well known to exhibit facet-dependent catalytic properties. However, few studies have investigated how those nanocrystals evolve and (de)activate during reactions, calling for the development of nanoscale time-resolved operando methods. In this context, we have designed Cu nanocubes as a model system to elucidate the underlying driving force of dynamic nanocatalyst reconstruction during the CO2 reduction reaction (CO2RR). Operando electrochemical liquid-cell scanning transmission electron microscopy (EC-STEM) and synchrotron-based X-ray spectroscopy reveal the size- and potential-dependent complete transformation from (100)-oriented Cu@Cu2O nanocubes to polycrystalline metallic Cu nanograins under CO2RR conditions. In addition, machine learning-assisted operando four-dimensional STEM reveals that large Cu nanograins derived from nanocubes form mainly crystalline domains, while their smaller counterparts are more amorphous due to faster evolution kinetics. In situ Raman spectroscopy and density functional theory calculations suggest that CO drives the ejection of single Cu atoms, resulting in few-nanometre Cu clusters and the surface migration of highly mobile copper carbonyl (Cu–CO) species. Combined, these multimodal operando methods and theoretical approaches pave the way for understanding the complex structural evolution of energy-related nanocatalysts under electrochemical conditions.

Yang, Yao↗

DOE OSTI2024

CO–induced roughening of Cu(111): formation and detection of reactive nanoclusters on metal surfaces

The formation of nanoclusters on metal surfaces in the presence of reactive environments is a phenomenon with important implications for catalysis. These nanoclusters are composed of atoms ejected from undercoordinated sites such as step edges, and their presence alters the catalytic properties of solid materials. We perform density functional theory (DFT) and kinetic Monte Carlo (KMC) simulations to investigate the formation and reactivity of copper clusters on Cu(111). Our results indicate a considerably higher reactivity of small copper nanoclusters, with up to seven atoms in size on roughened copper surfaces than on pristine Cu(111) and Cu(211). Regarding the restructuring events that give rise to nanoclusters under CO atmospheres, we determine that the ejection of Cu atoms from step edges and their migration therefrom to adjacent Cu(111) terraces are, by and large, driven by CO coverage effects. By means of KMC simulations, which account for CO–CO lateral interactions and CO–induced surface restructuring, we show that temperature programmed desorption (TPD) holds promise for the detection of highly reactive nanoclusters. Furthermore, our approach showcases how surface restructuring and surface–adsorbate bond breaking can be combined when modeling surface reactions and contributes to the development of an advanced understanding of the nature of active site under reaction conditions.

catalyst dynamic restructuring↗

DOE OSTI2024

Moiré Patterns in Pt Overlayers on Gold: A Graph Neural Network Interatomic Potential Study

Overlayer structures in bimetallic catalysts are relevant to a variety of catalytic reactions, particularly in electrocatalysis for fuel cell applications. Previous computational studies largely consider these overlayer structures to be those of a pseudomorphic overlayer, where there is a 1:1 atomic ratio between the overlayer and the support metal. Our previous work based on density functional theory (DFT) has shown that there exist nonstoichiometric overlayer structures that are more stable than the stoichiometric ones. Here, in this work, we developed a graph neural network interatomic potential (GNN-IP) to analyze structures formed in Pt overlayers on Au(111). The GNN-IP was used to explore properties of nonstoichiometric overlayers at length scales that are prohibitively expensive to pursue with planewave DFT calculations. In particular, we examined large Pt islands on top of a 48 × 48 Au(111) unit cell to explore the influence of the rotational angle (α) between the Pt overlayer and support Au(111) on both the stability of and the preferred atomic density in the overlayer. Island structures with smaller rotational angles between the Pt overlayer and Au(111) tend to be more stable. Further, smaller rotational angles tend to result in lower atomic density in the Pt overlayer.

chemical structure↗

DOE OSTI2024

The Effects of Polyolefin Structure and Source on Pyrolysis-Derived Plastic Oil Composition

Seven types of plastics were pyrolyzed in a fluidized bed reactor: post-consumer recycled (PCR) high-density polyethylene (HDPE), PCR polypropylene (PP), virgin resins of varying molecular weights of HDPE, virgin resins of low-density polyethylene (LDPE), linear low-density polyethylene (LLDPE), and (PP). Pyrolysis produced non-condensable gases (C1-C3), liquid phase products (C4-C40), and solids (C40+ and chars), with alkane, alkene, alkadiene, aromatic, and multi-cycloaromatics as the predominant compounds. Polymer structure had the greatest impact on product distribution, with minimal influence from molecular weight. Branches in polyethylene (PE) acted as thermal defects initiating degradation. Higher branch density in PE led to increased concentrations of aromatics, branched alkanes, and internal alkenes. PP and PE exhibited distinct degradation mechanisms, with PP requiring less energy for decomposition and yielding more oil. Here, pyrolysis oil from PCR HDPE and PCR PP contained a higher proportion of branched compounds. Additives in PCR plastics may promote isomerization during pyrolysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Adsorbate-induced adatom formation on Au-Cu bimetallic alloys and its possible consequences for CO 2 electroreduction

The adsorbate-induced formation of sub-nanometer clusters on transition-metal single crystals observed in previous high-pressure microscopic studies hinted at the in-situ formation of unique active sites even on large nanoparticle catalysts. We propose that the adatom formation energy can be used as an energetic descriptor for the initial step toward the adsorbate-induced metal-cluster formation process. This descriptor can be efficiently computed using density functional theory (DFT) calculations and applied for screening and identification of metal catalysts where this phenomenon may play an important role in generating active sites in-situ. As a proof of concept, here, we construct an adatom formation energy database for three Au x Cu y alloys (x:y = 3:1, 1:1, or 1:3) and eighteen adsorbates (H, C, N, O, F, S, Cl, Br, I, CH x , NH x (x = 1 – 3), CO, NO, and OH) commonly involved in catalytic reactions. The energetics of adatom formation were examined in all cases where the (111) terrace, (211) step-edge, and (874) kink were the sources of the adatom. We demonstrate that the presence of an adsorbate could alter not only the energetics for adatom formation but also the elemental nature of the preferred adatom being formed. Using our database, we identified promising systems which favor adsorbate-induced adatom formation under near-ambient conditions. Specifically, CO-induced adatom formation on all three Au-Cu alloy surfaces could occur under CO 2 electroreduction (CO 2 RR) conditions. This phenomenon offers a qualitative explanation for the experimentally observed CO 2 RR activity on Au-Cu alloy catalysts. As a result, our methodology offers an easily expandable and efficient approach for large-scale catalyst screening with regards to adatom/cluster formation under reaction conditions and provides insight into the possible nature of active sites on alloy catalysts from a novel perspective.

Active site↗

DOE OSTI2024

Origins of enhanced oxygen reduction activity of transition metal nitrides

Transition metal nitride (TMN-) based materials have recently emerged as promising non-precious-metal-containing electrocatalysts for the oxygen reduction reaction (ORR) in alkaline media. However, the lack of fundamental understanding of the oxide surface has limited insights into structure–(re)activity relationships and rational catalyst design. Here, in this work, we demonstrate how a well-defined TMN can dictate/control the as-formed oxide surface and the resulting ORR electrocatalytic activity. Structural characterization of MnN nanocuboids revealed that an electrocatalytically active Mn 3 O 4 shell grew epitaxially on the MnN core, with an expansive strain along the [010] direction to the surface Mn 3 O 4 . The strained Mn 3 O 4 shell on the MnN core exhibited an intrinsic activity that was over 300% higher than that of pure Mn 3 O 4 . A combined electrochemical and computational investigation indicated/suggested that the enhancement probably originates from a more hydroxylated oxide surface resulting from the expansive strain. This work establishes a clear and definitive atomistic picture of the nitride/oxide interface and provides a comprehensive mechanistic understanding of the structure–reactivity relationship in TMNs, critical for other catalytic interfaces for different electrochemical processes.

electrocatalysis↗

DOE OSTI2024

Catalytic production of δ-valerolactone (DVL) from biobased 2-hydroxytetrahydropyran (HTHP) – Combined experimental and modeling study

δ-Valerolactone (DVL) is a five-carbon (C5) cyclic ester that can undergo ring-opening polymerization to yield high-performance, biocompatible polyesters. But current market prices of C5 chemicals like DVL are very high due to poor availability of C5 feedstock in petroleum. Herein, we demonstrate a novel route to DVL synthesis via dehydrogenation of biomass-derived 2-hydroxytetrahydropyran (HTHP) over Cu/SiO 2 without the use of toxic reagents. Since HTHP exists in thermal equilibrium with 3,4-dihydropyran (DHP) via dehydration, and with 2,2’-oxybis(tetrahydropyran) and 5-(tetrahydropyran-2-yloxy)pentanal via acetalization, we have also determined the thermochemistry (ΔH rxn and ΔG rxn ) of each competing reaction using density functional theory (DFT) calculations at the M06-2X/cc-pVTZ level. Finally, by developing a kinetic model of all 8 reactions involved, we have achieved 84% selectivity to DVL at 150°C in a packed bed reactor for over 72 hours of time on stream.

42 ENGINEERING↗

DOE OSTI2024

Adsorption Properties of Au−Ni Surface Alloys with a Nonstoichiometric Moiré Structure: A Density Functional Theory Study

Due to the large lattice mismatch between gold and nickel, gold–nickel surface alloys can form unique nonstoichiometric overlayer structures characterized by a moiré pattern and subsurface defects. For this work, we performed density functional theory (DFT) calculations to study the adsorption of molecular oxygen, atomic hydrogen, and atomic carbon on a gold–nickel(111) surface alloy with 0.46 monolayer gold randomly distributed in the surface layer. We observed six distinct adsorption structures for molecular oxygen characterized by intramolecular stretching frequencies of <700, 729, 795, 857, 929, and 1004 cm –1 , which describe well the experimentally observed high-resolution electron energy-loss spectra. Surface atomic hydrogen adsorption is associated with adsorbate–surface modes in the ∼1000 cm –1 range, while subsurface hydrogen can have features as low as ∼400 cm –1 . We observed a unique adsorption structure for atomic carbon inside the surface dislocation loop defect, which explains the experimentally observed low carbon-surface mode at ∼340 cm –1 . Our study sheds light on the unique adsorption properties of the gold–nickel surface alloys and helps with rationalizing vibrational frequency experimental studies for this system.

adsorption↗

DOE OSTI2024

Exploring the Impact of Active Site Structure on the Conversion of Methane to Methanol in Cu‐Exchanged Zeolites

Abstract In the past, Cu‐oxo or ‐hydroxy clusters hosted in zeolites have been suggested to enable the selective conversion of methane to methanol, but the impact of the active site's stoichiometry and structure on methanol production is still poorly understood. Herein, we apply theoretical modeling in conjunction with experiments to study the impact of these two factors on partial methane oxidation in the Cu‐exchanged zeolite SSZ‐13. Phase diagrams developed from first‐principles suggest that Cu‐hydroxy or Cu‐oxo dimers are stabilized when O 2 or N 2 O are used to activate the catalyst, respectively. We confirm these predictions experimentally and determine that in a stepwise conversion process, Cu‐oxo dimers can convert twice as much methane to methanol compared to Cu‐hydroxyl dimers. Our theoretical models rationalize how Cu‐di‐oxo dimers can convert up to two methane molecules to methanol, while Cu‐di‐hydroxyl dimers can convert only one methane molecule to methanol per catalytic cycle. These findings imply that in Cu clusters, at least one oxo group or two hydroxyl groups are needed to convert one methane molecule to methanol per cycle. This simple structure–activity relationship allows to intuitively understand the potential of small oxygenated or hydroxylated transition metal clusters to convert methane to methanol.

Chemistry↗

DOE OSTI2024

Kinetic Study of Polyvinyl Chloride Pyrolysis with Characterization of Dehydrochlorinated PVC

In this paper, we study the kinetics of polyvinyl chloride (PVC) decomposition using a combination of experimental and computational approaches. Here, we develop a simplified kinetic model that contains only two steps: dehydrochlorination of PVC and further decomposition of the PVC residue. The model is consistent with density functional theory (DFT) calculations and experimental data. Dehydrochlorination is an autocatalytic reaction that starts with a tertiary chloride (Cl) and generates hydrogen chloride (HCl) and benzene as the main products. Benzene and HCl formation rates showed similar trends, indicating that HCl likely catalyzes a homolytic carbon-carbon (C-C) bond cleavage, which gives rise to benzene and an aliphatic fragment. We characterized the structure of dehydrochlorinated PVC (PVC residue) by using thermal gravimetric analysis (TGA), Fourier-transform infrared spectroscopy (FTIR), and nuclear magnetic resonance spectroscopy (NMR). FTIR and NMR results indicate that the PVC residue contains 20% quaternary carbon content, indicating a high concentration of cross-linked molecules. We predict that the most probable structure in the cross-linked centers of the PVC residue is cyclohexadiene, which is supported by DFT calculations, FTIR, and NMR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Actuating Liquid Crystals Rapidly and Reversibly by Using Chemical Catalysis

Abstract Microtubules and catalytic motor proteins underlie the microscale actuation of living materials, and they have been used in reconstituted systems to harness chemical energy to drive new states of organization of soft matter (e.g., liquid crystals (LCs)). Such materials, however, are fragile and challenging to translate to technological contexts. Rapid (sub‐second) and reversible changes in the orientations of LCs at room temperature using reactions between gaseous hydrogen and oxygen that are catalyzed by Pd/Au surfaces are reported. Surface chemical analysis and computational chemistry studies confirm that dissociative adsorption of H 2 on the Pd/Au films reduces preadsorbed O and generates 1 ML of adsorbed H, driving nitrile‐containing LCs from a perpendicular to a planar orientation. Subsequent exposure to O 2 leads to oxidation of the adsorbed H, reformation of adsorbed O on the Pd/Au surface, and a return of the LC to its initial orientation. The roles of surface composition and reaction kinetics in determining the LC dynamics are described along with a proof‐of‐concept demonstration of microactuation of beads. These results provide fresh ideas for utilizing chemical energy and catalysis to reversibly actuate functional LCs on the microscale.

Chemistry↗

DOE OSTI2024

Structure sensitivity in adsorbate-induced adatom formation on FCC transition-metal surfaces

Recent surface science discoveries reveal that adsorbates may induce in situ sub-nanometer cluster formation on transition-metal surfaces. To elucidate the structure sensitivity behind this phenomenon, we performed density functional theory calculations to construct an adatom formation energy database for eight fcc metals and 26 adsorbates commonly involved in catalytic reactions. We show that the adatom formation on (100) surfaces is generally easier than that on (111) surfaces. Many adsorbate/metal pairs exist, mostly on Pd, Ni, Rh, Pt, and Ir, for which adsorbates might induce adatom formation under near-ambient conditions on the (100) facet, but not on (111), highlighting the role of more open facets in adatom formation on metals intrinsically harder than Ag, Au, and Cu. Furthermore, our study offers a new perspective towards understanding structure sensitivity in heterogeneous thermal- and electro-catalytic systems such as methane steam reforming, Fischer-Tropsch synthesis, and ammonia decomposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Atomically Dispersed Zn/Co–N–C as ORR Electrocatalysts for Alkaline Fuel Cells

Hydrogen fuel cells have drawn increasing attention as one of the most promising next-generation power sources for future automotive transportation. Developing efficient, durable, and low-cost electrocatalysts, to accelerate the sluggish oxygen reduction reaction (ORR) kinetics, is urgently needed to advance fuel cell technologies. Herein, we report on metal–organic frameworks-derived nonprecious dual metal single-atom catalysts (SACs) (Zn/Co–N–C), consisting of Co–N 4 and Zn–N 4 local structures. These catalysts exhibited superior ORR activity with a half-wave potential (E 1/2 ) of 0.938 V versus RHE (reversible hydrogen electrode) and robust stability (ΔE 1/2 = –8.5 mV) after 50k electrochemical cycles. Moreover, this remarkable performance was validated under realistic fuel cell working conditions, achieving a record-high peak power density of ~1 W cm –2 among the reported SACs for alkaline fuel cells. Operando X-ray absorption spectroscopy was conducted to identify the active sites and reveal catalytic mechanistic insights. The results indicated that the Co atom in the Co–N 4 structure was the main catalytically active center, where one axial oxygenated species binds to form an O ads –Co–N 4 moiety during the ORR. In addition, theoretical studies, based on a potential-dependent microkinetic model and core-level shift calculations, showed good agreement with the experimental results and provided insights into the bonding of oxygen species on Co–N 4 centers during the ORR. Furthermore, this work provides a comprehensive mechanistic understanding of the active sites in the Zn/Co–N–C catalysts and will pave the way for the future design and advancement of high-performance single-site electrocatalysts for fuel cells and other energy applications.

25 ENERGY STORAGE↗

DOE OSTI2024

Thermal Stability of Au Rhombic Dodecahedral Nanocrystals Can Be Greatly Enhanced by Coating Their Surface with an Ultrathin Shell of Pt

Rhombic dodecahedral nanocrystals have been considered particularly difficult to synthesize because they are enclosed by {110}, a low-index facet with the greatest surface energy. Recently, we demonstrated the use of seed-mediated growth for the facile and robust synthesis of Au rhombic dodecahedral nanocrystals (AuRD). While the unique shape and surface structure of AuRD are desirable for potential applications in plasmonics and catalysis, respectively, their high surface energy makes them highly susceptible to thermal degradation. Here we demonstrate that it is feasible to greatly improve the thermal stability with some sacrifice to the plasmonic properties of the original AuRD by coating their surface with an ultrathin shell made of Pt. Our in situ electron microscopy analysis indicates that the ultrathin Pt coating can increase the thermal stability from 60 up to 450 °C, a trend that is also supported by the results from a computational study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Molecular-scale Insights into Cooperativity Switching of x TAB Adsorption on Gold Nanoparticles

Quantifying adsorption behaviors is crucial for various applications such as catalysis, separation, and sensing, yet it is generally challenging to access in solution. Here, we report a combined experimental and computational study of the adsorption behaviors of alkyl-trimethylammonium bromides (xTAB), a class of ligands important for colloidal nanoparticle stabilization and shape control, with various alkyl chain lengths x on Au nanoparticles. We use density functional theory (DFT) to calculate xTAB binding energies on Au{111} and Au{110} surfaces with standing-up and lying-down configurations, which provides insights into the adsorption affinity and cooperativity differences of xTAB on these two facets. We demonstrate the key role of van der Waals interactions in determining the xTAB adsorption behavior. These computational results predict and explain the experimental discovery of xTAB’s adsorption behavior switch from stronger affinity, negative cooperativity to weaker affinity, positive cooperativity when the concentration of xTAB increases in solution. We also show that in the standing-up configuration, bilayer adsorption may occur on both facets, which can lead to different differential binding energies and consequently adsorption crossover between the two facets when the ligand concentration increases. Our combined experimental and computational approaches demonstrate a paradigm for gaining molecular-scale insights into adsorbate–surface interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Fast and Non‐equilibrium Uptake of Hydrogen by Pd Icosahedral Nanocrystals

Abstract We report for the first time that Pd nanocrystals can absorb H via a “single‐phase pathway” when particles with a proper combination of shape and size are used. Specifically, when Pd icosahedral nanocrystals of 7‐ and 12‐nm in size are exposed to H atoms, the H‐saturated twin boundaries can divide each particle into 20 smaller single‐crystal units in which the formation of phase boundaries is no longer favored. As such, absorption of H atoms is dominated by the single‐phase pathway and one can readily obtain PdH x with anyx in the range of 0–0.7. When switched to Pd octahedral nanocrystals, the single‐phase pathway is only observed for particles of 7 nm in size. We also establish that the H‐absorption kinetics will be accelerated if there is a tensile strain in the nanocrystals due to the increase in lattice spacing. Besides the unique H‐absorption behaviors, the PdH x ( x =0–0.7) icosahedral nanocrystals show remarkable thermal and catalytic stability toward the formic acid oxidation due tothe decrease in chemical potential for H atoms in a Pd lattice under tensile strain.

Zhou, Siyu↗

DOE OSTI2023

Hydroformylation of pyrolysis oils to aldehydes and alcohols from polyolefin waste

Waste plastics are an abundant feedstock for the production of renewable chemicals. Pyrolysis of waste plastics produces pyrolysis oils with high concentrations of olefins (>50 weight %). The traditional petrochemical industry uses several energy-intensive steps to produce olefins from fossil feedstocks such as naphtha, natural gas, and crude oil. In this work, we demonstrate that pyrolysis oil can be used to produce aldehydes through hydroformylation, taking advantage of the olefin functionality. These aldehydes can then be reduced to mono- and dialcohols, oxidized to mono- and dicarboxylic acids, or aminated to mono- and diamines by using homogeneous and heterogeneous catalysis. This route produces high-value oxygenated chemicals from low-value postconsumer recycled polyethylene. Here, we project that the chemicals produced by this route could lower greenhouse gas emissions ~60% compared with their production through petroleum feedstocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Adsorbate-Induced Adatom Formation on Lithium, Iron, Cobalt, Ruthenium, and Rhenium Surfaces

Recent experimental and theoretical studies have demonstrated the reaction-driven metal–metal bond breaking in metal catalytic surfaces even under relatively mild conditions. Here, we construct a density functional theory (DFT) database for the adsorbate-induced adatom formation energy on the close-packed facets of three hexagonal close-packed metals (Co, Ru, and Re) and two body-centered cubic metals (Li and Fe), where the source of the ejected metal atom is either a step edge or a close-packed surface. For Co and Ru, we also considered their metastable face-centered cubic structures. We studied 18 different adsorbates relevant to catalytic processes and predicted noticeably easier adatom formation on Li and Fe compared to the other three metals. The NH 3 - and CO-induced adatom formation on Fe(110) is possible at room temperature, a result relevant to NH 3 synthesis and Fischer-Tropsch synthesis, respectively. There also exist other systems with favorable adsorbate effects for adatom formation relevant to catalytic processes at elevated temperatures (500–700 K). Our results offer insight into the reaction-driven formation of metal clusters, which could play the role of active sites in reactions catalyzed by Li, Fe, Co, Ru, and Re catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗