Operando Characterization and Molecular Simulations Reveal the Growth Kinetics of Graphene on Liquid Copper During Chemical Vapor Deposition.

In recent years, liquid metal catalysts have emerged as a compelling choice for the controllable, large-scale, and high-quality synthesis of two-dimensional materials. At present, there is little mechanistic understanding of the intricate catalytic process, though, of its governing factors or what renders it superior to growth at the corresponding solid catalysts. Here, we report on a combined experimental and computational study of the kinetics of graphene growth during chemical vapor deposition on a liquid copper catalyst. By monitoring the growing graphene flakes in real time using in situ radiation-mode optical microscopy, we explore the growth morphology and kinetics over a wide range of CH4-to-H2 pressure ratios and deposition temperatures. Constant growth rates of the flakes' radius indicate a growth mode limited by precursor attachment, whereas methane-flux-dependent flake shapes point to limited precursor availability. Large-scale free energy simulations enabled by an efficient machine-learning moment tensor potential trained to density functional theory data provide quantitative barriers for key atomic-scale growth processes. The wealth of experimental and theoretical data can be consistently combined into a microkinetic model that reveals mixed growth kinetics that, in contrast to the situation at solid Cu, is partly controlled by precursor attachment alongside precursor availability. Key mechanistic aspects that directly point toward the improved graphene quality are a largely suppressed carbon dimer attachment due to the facile incorporation of this precursor species into the liquid surface and a low-barrier ring-opening process that self-heals 5-membered rings resulting from remaining dimer attachments.

Tripling of the scattering vector range of X-ray reflectivity on liquid surfaces using a double-crystal deflector.

The maximum range of perpendicular momentum transfer (q z) has been tripled for X-ray scattering from liquid surfaces when using a double-crystal deflector setup to tilt the incident X-ray beam. This is achieved by employing a higher-energy X-ray beam to access Miller indices of reflecting crystal atomic planes that are three times higher than usual. The deviation from the exact Bragg angle condition induced by misalignment between the X-ray beam axis and the main rotation axis of the double-crystal deflector is calculated, and a fast and straightforward procedure to align them is deduced. An experimental method of measuring scattering intensity along the q z direction on liquid surfaces up to q z = 7 Å-1 is presented, with liquid copper serving as a reference system for benchmarking purposes.

Characterization of CO Adsorbed to Clean and Partially Oxidized Cu(211) and Cu(111).

Copper-based catalysts gain activity through the presence of poorly coordinated Cu atoms and incomplete oxidation at the surface. The catalytic mechanisms can in principle be observed by controlled dosing of reactants to single-crystal substrates. However, the interconnected influences of surface defects, partial oxidation, and adsorbate coverage present a large matrix of conditions that have not been fully explored in the literature. We recently characterized oxygen and carbon monoxide coadsorption on Cu(111), a nominally defect-free surface, and now extend our study to the stepped surface Cu(211). Temperature-programmed desorption of CO adsorbed to bare metal surfaces confirms that two sites dominate desorption from a saturated layer: atop terrace atoms of local (111) character and atop step edge atoms with CO bound more strongly to the latter. At low coverage, discrete CO resonances in reflection adsorption infrared spectra can be assigned to these sites: 2077 cm-1 for extended (111) terraces, 2093 cm-1 for step sites, and additional kink-adsorbed molecules at 2110 cm-1. With increasing coverage, in contrast to Cu(111), the infrared spectral features on Cu(211) evolve and shift as a consequence of dipole-dipole coupling between differentially occupied types of sites. Auger electron spectroscopy shows that exposure to background O2 oxidizes the (211) surface at a rate nearly 1 order of magnitude greater than (111); we argue that the resulting surface is stoichiometric Cu2O, as previously found for Cu(111). This oxide binds CO less strongly than the bare metal and the underlying crystal cut continues to influence the adsorption sites available to CO. On oxidized (111) terraces, broad absorption peaks at 2115-2120 cm-1; on oxidized Cu(211), CO adsorbed to step sites appears as a resolved secondary peak at 2144 cm-1. This suite of spectroscopic signatures, obtained under carefully controlled conditions, will help to determine the origin and fate of adsorbed species in future studies of reaction mechanisms on copper.

Spectroscopic investigation of a Co(0001) model catalyst during exposure to H2 and CO at near-ambient pressures.

Cobalt catalysts, although already used industrially for Fischer-Tropsch synthesis, are prone to a number of deactivation mechanisms such as oxidation of the active phase, and the deposition of carbon and reaction products. We have performed near-ambient-pressure X-ray photoelectron spectroscopy on Co(0001) model catalysts during exposure to gases relevant to Fischer-Tropsch synthesis, i.e., CO and H2, at 0.25 mbar total pressure. At this pressure, CO seems to be more efficient at keeping the Co(0001) surface metallic than H2, which is the opposite behavior as reported in the literature for other pressure ranges. We offer an interpretation of these differences based on the preferred adsorption and dissociation sites of CO and H2 compared to the oxidizing agent water (present as impurity in the gas feed and one of the products of the reaction). Additionally, detailed carbon spectra measured at the HIPPIE beamline of MAX IV allow for the distinction of different adsorbed species: CO and COx species are present in correlation to the presence of oxygen on the surface. Carbidic carbon and graphitic carbon can both be removed by hydrogen, whereas adsorbed hydrocarbons possibly poison the surface.

Graphene at Liquid Copper Catalysts: Atomic-Scale Agreement of Experimental and First-Principles Adsorption Height.

Liquid metal catalysts have recently attracted attention for synthesizing high-quality 2D materials facilitated via the catalysts' perfectly smooth surface. However, the microscopic catalytic processes occurring at the surface are still largely unclear because liquid metals escape the accessibility of traditional experimental and computational surface science approaches. Hence, numerous controversies are found regarding different applications, with graphene (Gr) growth on liquid copper (Cu) as a prominent prototype. In this work, novel in situ and in silico techniques are employed to achieve an atomic-level characterization of the graphene adsorption height above liquid Cu, reaching quantitative agreement within 0.1 Å between experiment and theory. The results are obtained via in situ synchrotron X-ray reflectivity (XRR) measurements over wide-range q-vectors and large-scale molecular dynamics simulations based on efficient machine-learning (ML) potentials trained to first-principles density functional theory (DFT) data. The computational insight is demonstrated to be robust against inherent DFT errors and reveals the nature of graphene binding to be highly comparable at liquid Cu and solid Cu(111). Transporting the predictive first-principles quality via ML potentials to the scales required for liquid metal catalysis thus provides a powerful approach to reach microscopic understanding, analogous to the established computational approaches for catalysis at solid surfaces.

X-ray reflectivity from curved surfaces as illustrated by a graphene layer on molten copper.

The X-ray reflectivity technique can provide out-of-plane electron-density profiles of surfaces, interfaces, and thin films, with atomic resolution accuracy. While current methodologies require high surface flatness, this becomes challenging for naturally curved surfaces, particularly for liquid metals, due to the very high surface tension. Here, the development of X-ray reflectivity measurements with beam sizes of a few tens of micrometres on highly curved liquid surfaces using a synchrotron diffractometer equipped with a double crystal beam deflector is presented. The proposed and developed method, which uses a standard reflectivity θ-2θ scan, is successfully applied to study in situ the bare surface of molten copper and molten copper covered by a graphene layer grown in situ by chemical vapor deposition. It was found that the roughness of the bare liquid surface of copper at 1400 K is 1.25 ± 0.10 Å, while the graphene layer is separated from the liquid surface by a distance of 1.55 ± 0.08 Šand has a roughness of 1.26 ± 0.09 Å.

Investigation of Active Catalysts at Work.

Even after being in business for at least the last 100 years, research into the field of (heterogeneous) catalysis is still vibrant, both in academia and in industry. One of the reasons for this is that around 90% of all chemicals and materials used in everyday life are produced employing catalysis. In 2020, the global catalyst market size reached $35 billion, and it is still steadily increasing every year. Additionally, catalysts will be the driving force behind the transition toward sustainable energy. However, even after having been investigated for 100 years, we still have not reached the holy grail of developing catalysts from rational design instead of from trial-and-error. There are two main reasons for this, indicated by the two so-called "gaps" between (academic) research and actual catalysis. The first one is the "pressure gap", indicating the 13 orders of magnitude difference in pressure between the ultrahigh vacuum lab conditions and the atmospheric pressures (and higher) of industrial catalysis. The second one is the "materials gap", indicating the difference in complexity between single-crystal model catalysts of academic research and the real catalysts, consisting of metallic nanoparticles on supports, promoters, fillers, and binders. Although over the past decades significant efforts have been made in closing these gaps, many steps still have to be taken. In this Account, I will discuss the steps we have taken at Leiden University to further our fundamental understanding of heterogeneous catalysis at the (near-)atomic scale. I will focus on bridging the pressure gap, though we are also working on closing the materials gap. Over the past years, we developed state-of-the-art equipment that is able to investigate the (near-)atomic-scale structure of the catalyst surface during the chemical reaction using several surface-science-based techniques such as scanning tunneling microscopy, atomic force microscopy, optical microscopy, and X-ray-based techniques (surface X-ray diffraction, grazing-incidence small-angle X-ray scattering, and X-ray reflectivity, in collaboration with ESRF). Simultaneously with imaging the surface, we can investigate the catalyst's performance via mass spectrometry, enabling us to link changes in the catalyst structure to its activity, selectivity, or stability. Although we are currently investigating many industrially relevant catalytic systems, I will here focus the discussion on the oxidation of platinum during, for example, CO and NO oxidation, the NO reduction reaction on platinum, and the growth of graphene on liquid (molten) copper. I will show that to be able to obtain the full picture of heterogeneous catalysis, the ability to investigate the catalyst at the (near-)atomic scale during the chemical reaction is a must.

Real-Time Multiscale Monitoring and Tailoring of Graphene Growth on Liquid Copper.

The synthesis of large, defect-free two-dimensional materials (2DMs) such as graphene is a major challenge toward industrial applications. Chemical vapor deposition (CVD) on liquid metal catalysts (LMCats) is a recently developed process for the fast synthesis of high-quality single crystals of 2DMs. However, up to now, the lack of in situ techniques enabling direct feedback on the growth has limited our understanding of the process dynamics and primarily led to empirical growth recipes. Thus, an in situ multiscale monitoring of the 2DMs structure, coupled with a real-time control of the growth parameters, is necessary for efficient synthesis. Here we report real-time monitoring of graphene growth on liquid copper (at 1370 K under atmospheric pressure CVD conditions) via four complementary in situ methods: synchrotron X-ray diffraction and reflectivity, Raman spectroscopy, and radiation-mode optical microscopy. This has allowed us to control graphene growth parameters such as shape, dispersion, and the hexagonal supra-organization with very high accuracy. Furthermore, the switch from continuous polycrystalline film to the growth of millimeter-sized defect-free single crystals could also be accomplished. The presented results have far-reaching consequences for studying and tailoring 2D material formation processes on LMCats under CVD growth conditions. Finally, the experimental observations are supported by multiscale modeling that has thrown light into the underlying mechanisms of graphene growth.

Structural Characterization of a Novel Two-Dimensional Material: Cobalt Sulfide Sheets on Au(111).

Transition metal dichalcogenides (TMDCs) are a type of two-dimensional (2D) material that has been widely investigated by both experimentalists and theoreticians because of their unique properties. In the case of cobalt sulfide, density functional theory (DFT) calculations on free-standing S-Co-S sheets suggest there are no stable 2D cobalt sulfide polymorphs, whereas experimental observations clearly show TMDC-like structures on Au(111). In this study, we resolve this disagreement by using a combination of experimental techniques and DFT calculations, considering the substrate explicitly. We find a 2D CoS(0001)-like sheet on Au(111) that delivers excellent agreement between theory and experiment. Uniquely this sheet exhibits a metallic character, contrary to most TMDCs, and exists due to the stabilizing interactions with the Au(111) substrate.

In situ observations of an active MoS2 model hydrodesulfurization catalyst.

The hydrodesulfurization process is one of the cornerstones of the chemical industry, removing harmful sulfur from oil to produce clean hydrocarbons. The reaction is catalyzed by the edges of MoS2 nanoislands and is operated in hydrogen-oil mixtures at 5-160 bar and 260-380 °C. Until now, it has remained unclear how these harsh conditions affect the structure of the catalyst. Using a special-purpose high-pressure scanning tunneling microscope, we provide direct observations of an active MoS2 model catalyst under reaction conditions. We show that the active edge sites adapt their sulfur, hydrogen, and hydrocarbon coverages depending on the gas environment. By comparing these observations to density functional theory calculations, we propose that the dominant edge structure during the desulfurization of CH3SH contains a mixture of adsorbed sulfur and CH3SH.

Transferability of the Specific Reaction Parameter Density Functional for H2 + Pt(111) to H2 + Pt(211).

The accurate description of heterogeneously catalyzed reactions may require chemically accurate evaluation of barriers for reactions of molecules at the edges of metal nanoparticles. It was recently shown that a semiempirical density functional describing the interaction of a molecule dissociating on a flat metal surface (CHD3 + Pt(111)) is transferable to the same molecule reacting on a stepped surface of the same metal (Pt(211)). However, validation of the method for additional systems is desirable. To address the question whether the specific reaction parameter (SRP) functional that describes H2 + Pt(111) with chemical accuracy is also capable of accurately describing H2 + Pt(211), we have performed molecular beam simulations with the quasi-classical trajectory (QCT) method, using the SRP functional developed for H2 + Pt(111). Our calculations used the Born-Oppenheimer static surface model. The accuracy of the QCT method was assessed by comparison with quantum dynamics results for reaction of the ro-vibrational ground state of H2. The theoretical results for sticking of H2 and D2 on Pt(211) are in quite good agreement with the experiment, but uncertainties remain because of a lack of accuracy of the QCT simulations at low incidence energies and possible inaccuracies in the reported experimental incidence energies at high energies. We also investigated the nonadiabatic effect of electron-hole pair excitation on the reactivity using the molecular dynamics with the electron friction (MDEF) method, employing the local density friction approximation (LDFA). Only small effects of electron-hole pair excitation on sticking are found.

Structural Dynamics of Al2O3/NiAl(110) During Film Growth in NO2.

While continuum descriptions of oxide film growth are well established, the local structural dynamics during oxide growth are largely unexplored. Here, we investigate this using scanning tunneling microscopy (STM) and X-ray photoelectron spectroscopy (XPS) for the example of alumina film growth on NiAl(110) following NO2 exposure. To maintain a well-defined system, we have adopted a cyclic growth approach of NO2 adsorption and annealing. NO2 adsorption at 693 K results in the formation of a vacancy island pattern in the NiAl(110) substrate, which is filled with AlOx by diffusion of O through the alumina film. The patches of AlOx coalesce to form smooth terraces upon annealing to 1200 K. By repeated cycling, we have grown films of up to 0.9 nm thick. While peak shifts in the XPS spectra indicate that the film maintains its insulating character upon thickening, our STM data show that there is a finite density of states within the band gap. The thickening of the alumina film is accompanied by the formation of trenches in the surface, which we interpret to be the result of film stress relief.

Observing the oxidation of platinum.

Despite its importance in oxidation catalysis, the active phase of Pt remains uncertain, even for the Pt(111) single-crystal surface. Here, using a ReactorSTM, the catalytically relevant structures are identified as two surface oxides, different from bulk α-PtO2, previously observed. They are constructed from expanded oxide rows with a lattice constant close to that of α-PtO2, either assembling into spoked wheels, 1-5 bar O2, or closely packed in parallel lines, above 2.2 bar. Both are only ordered at elevated temperatures (400-500 K). The triangular oxide can also form on the square lattice of Pt(100). Under NO and CO oxidation conditions, similar features are observed. Furthermore, both oxides are unstable outside the O2 atmosphere, indicating the presence of active O atoms, crucial for oxidation catalysts.Improving platinum as an oxidation catalyst requires understanding its structure under catalytic conditions. Here, the authors discover that catalytically important surface oxides form only when Pt is exposed to high pressure and temperature, highlighting the need to study catalysts in realistic environments.

Simultaneous scanning tunneling microscopy and synchrotron X-ray measurements in a gas environment.

A combined X-ray and scanning tunneling microscopy (STM) instrument is presented that enables the local detection of X-ray absorption on surfaces in a gas environment. To suppress the collection of ion currents generated in the gas phase, coaxially shielded STM tips were used. The conductive outer shield of the coaxial tips can be biased to deflect ions away from the tip core. When tunneling, the X-ray-induced current is separated from the regular, 'topographic' tunneling current using a novel high-speed separation scheme. We demonstrate the capabilities of the instrument by measuring the local X-ray-induced current on Au(1 1 1) in 800 mbar Ar.

In Situ Optical Reflectance Difference Observations of CO Oxidation over Pd(100).

Using a home-built reflectometer, we have investigated the changes in the optical reflectivity of a Pd(100) model catalyst during CO oxidation under high-pressure, high-temperature conditions. We observe changes in optical contrast when exposing the surface to CO oxidation conditions at 200 mbar from room temperature up to 400 °C. These changes in reflectivity are a result both of the formation of a surface oxide layer and of a change in surface roughness because of gas exposure. However, the reflectivity is more sensitive to the presence of a thin, flat oxide layer than to surface roughness. CO oxidation plays an important role in the decrease of the reflectivity. Since adding a reducing agent to the gas mixture renders it unlikely that the oxide thickness increases, we conclude that the observed decrease in reflectivity is dominated by increased surface roughness because of the catalytic reaction. We contribute this observed surface roughening to a Mars-van Krevelen-type reaction mechanism.

Surface science under reaction conditions: CO oxidation on Pt and Pd model catalysts.

Platinum and palladium are frequently used as catalytic materials, for example for the oxidation of CO. This is one of the most widely studied reactions in the field of surface science. Although seemingly uncomplicated, it remains an active and interesting topic, which is partially explained by the push to conduct experiments on model systems under relevant reaction conditions. Recent developments in the surface-science methodology have allowed obtaining chemical and structural information on the active phase of model catalysts. Tools of the trade include near-ambient-pressure X-ray photoelectron spectroscopy, high-pressure scanning tunneling microscopy, high-pressure surface X-ray diffraction, and high-pressure vibrational spectroscopy. Interpretation is often aided by density functional theory in combination with thermodynamic and kinetic modeling. In this review, results for the catalytic oxidation of CO obtained by these techniques are compared. On several of the Pt and Pd surfaces, new structures develop in excess O2. For Pt, this requires a much larger excess of O2 than for Pd. Most of these structures also develop in pure O2 and are identified as (surface) oxides. A large body of evidence supports the conjecture that these oxides are more reactive than the corresponding O-covered metallic surfaces under similar conditions, although still debated in the literature. An outlook on this developing field, including directions that move away from CO oxidation towards more complex chemistry, concludes this review.

From dull to shiny: A novel setup for reflectance difference analysis under catalytic conditions.

We have developed an experimental setup for optically monitoring a catalytically active surface under reaction conditions. A flow reactor with optical access allows us to image the behavior of an active catalyst surface down to the millimeter length scale. We use reflectance difference measurements with 625 nm light to investigate CO oxidation on Pd(100) at 300 mbar and 320 °C. We conclude that the changes in visible contrast result from the formation of an oxide layer after surface oxidation.

Combined scanning probe microscopy and x-ray scattering instrument for in situ catalysis investigations.

We have developed a new instrument combining a scanning probe microscope (SPM) and an X-ray scattering platform for ambient-pressure catalysis studies. The two instruments are integrated with a flow reactor and an ultra-high vacuum system that can be mounted easily on the diffractometer at a synchrotron end station. This makes it possible to perform SPM and X-ray scattering experiments in the same instrument under identical conditions that are relevant for catalysis.

Evidence of stable high-temperature D(x)-CO intermediates on the Ru(0001) surface.

We demonstrate the formation of complexes involving attractive interactions between D and CO on Ru(0001) that are stable at significantly higher temperatures than have previously been reported for such intermediate species on this surface. These complexes are evident by the appearance of new desorption features upon heating of the sample. They decompose in stages as the sample temperature is increased, with the most stable component desorbing at >500 K. The D:CO ratio remaining on the surface during the final stages of desorption tends towards 1:1. The new features are populated during normally incident molecular beam dosing of D(2) on to CO pre-covered Ru(0001) surfaces (180 K) when the CO coverage exceeds 50% of the saturation value. The amount of complex formed decreases somewhat with increasing CO pre-coverage. It is almost absent in the case of dosing on to the fully saturated surface. The results are interpreted in terms of both local and long-range rearrangements of the overlayer that give rise to the observed CO coverage dependence and limit the amount of complex that can be formed.