Abstract
Recent state-to-state experiments of methane scattering from Ni(111) and graphene-covered Ni(111) combined with quantum mechanical simulations suggest an intriguing correlation between the surface-induced vibrational energy redistribution (SIVR) during the molecule/surface scattering event and the catalytic activity for methane dissociation of the target surface (Werdecker, Phys. Rev. Res., 2020, 2, 043251). Herein, we report new quantum state and angle-resolved measurements for methane scattering from Ni(111) and Au(111) probing the extent of antisymmetric-to-symmetric conversion of methane stretching motion for two surfaces with different catalytic activities. Consistent with the expectations, the extent of SIVR occurring on the more catalytically active Ni(111) surface, as measured by the scattered population ratio, is found to be several times stronger than that on the more inert Au(111) surface. We also present additional insights on the rovibrational scattering dynamics contained in the angle- and state-resolved data. The results together highlight the power of state-resolved scattering measurements as a tool for investigating methane–surface interactions.
1 Introduction
State-to-state scattering experiments provide a powerful and well-established technique to elucidate the atomic-level details of what happens when molecules collide with other gas-phase molecules or with surfaces (Auerbach et al., 2021). Often, the connection between the measured final state distributions and the property under study is direct and obvious. For example, the degree of rotational excitation in a collision is directly related to the angular anisotropy of the potential energy surface (PES) of the system under study and can be used to test theoretical methods for calculating the PES. In other cases, the connection of the scattering data to the underlying question of what is going on is more subtle. For example, highly vibrationally excited NO or CO molecules scatter almost vibrationally elastically from insulator surfaces but undergo facile multi-quantum vibrational relaxation if they strike a metal surface. Qualitatively, these observations show that the presence of a continuum of low-lying electronically excited states in the metal is being excited by molecular vibrations-a strong violation of the Born–Oppenheimer approximation. More detailed theoretical work established that charge transfer in the collisions is the key driving factor and that the propensity for charge transfer therefore correlates with the degree of vibrational relaxation (Wodtke et al., 2008; Shenvi et al., 2009; Wagner et al., 2019).
An intriguing example of this more subtle connection of observations to dynamical propensities emerges from the first state-to-state measurements of methane scattering from surfaces. When CH4 (, that is, CH4 prepared with a single vibrational quantum in the antisymmetric C-H stretch normal mode and one quantum of vibrational angular momentum, strikes a clean Ni(111) surface, nearly 40% of the molecules scatter into a state with one quantum of the symmetric C-H stretch normal mode (Werdecker et al., 2018). In this study, we refer to this phenomenon of surface collision-mediated conversion of vibrational energy by the term surface-induced vibrational energy redistribution, or SIVR for short. We note that a similar SIVR process for scattering of H2O from Cu(111) was predicted by state-to-state quantum scattering calculations by Zhang and Jiang (2019) and Zhang and Jiang (2020).
Surprisingly, SIVR was not observed for CH4 () scattering from Gr/Ni(111), the graphene-covered Ni(111) surface (Werdecker et al., 2020) (for convenience we refer to reference [Werdecker et al., 2020] from now on as “Werdecker et al.”). Quantum dynamics calculations by Bret Jackson and coworkers propose an explanation for this remarkable difference. Their calculations show that as the incident CH4 () molecules approach the reactive Ni(111) surface, the C-H bond closest to the surface starts to elongate, introducing a coupling between the two C-H stretching normal modes and which differ only by the relative phases of the C-H stretching motions. Inelastic scattering calculations for collisions with different surface sites on Ni(111) reveal a correlation between the calculated C-H bond elongation and the catalytic activity of the surface site, resulting in an increasing probability for SIVR with increasing catalytic activity of the impact site. Based on these results, the authors hypothesized that the extent of SIVR in a molecule–surface collision is correlated with the catalytic activity of the surface.
In this paper, we present new measurements to test this hypothesis. We extend the range of surfaces studied to include the closely packed Au(111) surface which, based on density function calculations1 of the barrier to dissociative adsorption, is expected to have a catalytic activity intermediate to that of Ni(111) (Jackson and Nave, 2013) and Gr/Ni(111) (Li et al., 2018).
In addition to providing data on a new system with catalytic activity intermediate between Ni(111) and Gr/Ni(111), we have also improved the method of determining the “branching ratio” , which quantifies the extent of SIVR through the ratio of the total and scattering flux following the surface collision of an incident CH4 () molecule. We use a rotatable bolometer detector in this work, which allows us to measure quantum state-resolved angular distributions of scattered molecules. We base our reported values of on the scattering angle- and rotational state-integrated scattered flux of CH4 in the and states. Previous determinations of were based on measurements at a single scattering angle due to limitations of the instrument employed in those studies. Indeed, as shown in Section 3.2, the angular distributions differ for the different scattered CH4 rovibrational states probed. Therefore, for the measured branching ratio to reflect the total and scattering probabilities accurately, it is advantageous to measure over a wide range of scattering angles and integrate the and scattered intensities over the scattering angle.
Our measurements show that the value of measured for Au(111) is intermediate to those of Ni(111) and Gr/Ni(111). The literature values of the activation barriers to dissociative chemisorption increase when proceeding from Ni(111) (Jackson and Nave, 2013) to Au(111) (Jackson) to (free-standing) graphene (Li et al., 2018). Thus, our results are in accordance with the hypothesis that SIVR correlates with the catalytic activity of the surface.
2 Materials and methods
Figure 1 shows a schematic illustration of the experimental setup. Full details can be found in Reilly et al. (2023). Briefly, a continuous supersonic expansion from a temperature-controlled nozzle into a vacuum is used to create a molecular beam of methane with the mean kinetic energy of 100 meV and a FWHM spread of 30%. Approximately 10% of the methane molecules in the molecular beam are prepared in the meta-CH4 (i.e., nuclear spin ) rotational state of the antisymmetric stretch fundamental (i.e., with one quantum of vibration) by infrared laser excitation with a tunable, single-mode, continuous-wave (CW) infrared (IR) optical parametric oscillator (OPO) (Argos Lockheed Martin-Aculight, “Pump laser” in Figure 1) (Chadwick et al., 2014). The remaining of the molecules populate rotational states in the ground vibrational state not addressed by the pump laser.
FIGURE 1
Schematic of the experimental setup used for state-to-state surface scattering of methane. A molecular beam is generated in vacuum chamber P1 and traverses the differential pumping stages P2 and P3. Quantum state-specific preparation of methane in the molecular beam occurs in P2 which also contains a chopper disk for beam modulation. The state-prepared molecular beam collides in P4 with the single crystal target surface at an incident angle with respect to the surface normal (marked by the dotted line). Molecules scattered at an angle are excited with quantum state resolution by the tagging laser and detected by the rotatable bolometer.
The state-prepared molecular beam collides with a single crystal surface sample in the ultrahigh vacuum surface science chamber with a base pressure of 1 × 10−10 mbar. A 4-axis manipulator permits displacement the surface in three dimensions and rotation of the incident angle made between the surface normal and the molecular beam.
Scattered molecules are detected and resolved by the quantum state and scattering angle using the bolometer infrared laser tagging (BILT) technique (Reilly et al., 2023). Here, a second, tunable, single-mode, CW IR OPO (TOPO TOPTICA Photonics, “tagging laser” in Figure 1) laser beam is chopped at 237 Hz with 50% duty cycle using an optical chopper wheel (not shown) and crosses the scattering plane perpendicularly 16 mm from the target surface. By tuning the tagging laser into resonance with different rovibrational transitions, we selectively excite the methane molecules that have scattered from the initially prepared state into the lower state of the resonant transition, transferring to those molecules one quantum of vibrational energy. The approximately 370 meV/molecule of energy absorbed by the molecules from the tagging laser is then transferred to a cryogenic bolometer detector when the molecules adsorb on a cold 4-mm-diameter diamond absorber to which the bolometer is attached.
Using a lock-in amplifier, we record the response of the detector to the chopping of the tagging laser. To the raw lock-in output, we apply a series of corrections to obtain a “scattering intensity” which is intended to provide a relative measure of the fraction of molecules scattering into the range of solid angles subtended by the bolometer detector centered at an angle in the scattering plane and occupying the state being tagged. These corrections are described in detail in Section S1 of the Supplementary Material of Reilly et al. (2023).
Synchronized, computer-controlled rotation of the bolometer and the tagging laser beam enables variation of the measured final angle , enabling, in turn, measurement of the state-resolved angular distributions of scattered methane molecules, as shown in Figure 2 in the following section.
FIGURE 2
State-resolved in-plane angular distributions for CH4() scattering from Ni(111) (blue markers) and Au(111) (orange markers). The incoming molecules have a mean kinetic energy and strike the surface () at an incident angle . The points correspond to measured data, and the dashed lines correspond to the associated fits (see Eq. 1). The best-fit parameter values are listed in the boxes next to their associated plots. The main text shows the explanation of quantum state identifiers (, , , , and ). Reported scattering intensities are obtained from raw detector signals after the application of corrections for tagging efficiency and time-dependent bolometer sensitivity and represent the differential scattered flux (in molecules per steradian per second) within a multiplicative constant independent of the surface and rovibrational state. The absolute value of the angle between the surface projection of the molecule’s initial velocity and the crystallographic direction of the Ni(111) and Au(111) surface samples is and , respectively.
For preparation of atomically flat Ni(111) and Au(111) surfaces, we use Ar+ sputtering followed by annealing to 973 K and 773 K, respectively. Crystallographic order of prepared surfaces is verified using low-energy electron diffraction (LEED), and the absence of surface contaminants is verified using Auger electron spectroscopy. After conducting all experiments, the contamination levels were verified to be below the detection threshold of approximately 1% ML.
3 Results and discussion
For both the Ni(111) and Au(111) surfaces, we recorded 14 state-resolved angular distributions. To identify the rovibrational level being tagged, we use the quantum numbers (, , , , and ) which, for brevity, we sometimes condense down to a single compound index . specifies the symmetry ( for symmetric and for antisymmetric) of the stretch fundamental (or “stretch” for simplicity), denotes the total angular momentum, denotes the Coriolis label () characterizing the coupling of rotational and vibrational angular momentum (di Lauro, 2020), denotes the wavefunction symmetry (), and denotes the polyad-level identifier (Brown et al., 1992). The rovibrational level prepared by the pump laser, for example, has quantum numbers (). Note that for the subset of levels tagged in our study, the first four identifiers (i.e., , , , and ) suffice to uniquely identify a level.
Figure 2 illustrates a selected subset of the measured angular distributions. The complete set of measured distributions is available in Supplementary Section S1. The rotational states tagged include all meta-CH4 states in the and stretch fundamentals with total angular momentum quantum number along with one rotational level for both stretches. States with nuclear spin were not probed because interconversion among the different nuclear spin isomers in direct molecule–surface collisions is presumed to occur with negligible probability. Our own measurements using a double-resonance pump–probe technique showed no evidence of such interconversion in CH4/Au(111) collisions.
In the following sections, we present an analysis and discussion of the measured distributions. First, in Section 3.1, we treat the central topic of SIVR before moving on to a detailed look at aspects of scattering dynamics relating to the scattering angle (Section 3.2) and rotational state (Section 3.3).
3.1 Surface-induced vibrational energy redistribution
The principle objective of the current work is to test the hypothesis advanced in Werdecker et al. and discussed in the Introduction of the current work. We remind the reader that this hypothesis asserts a positive correlation between the efficacy of conversion induced in a methane–surface collision and the catalytic activity of the surface for the dissociative chemisorption of methane. To quantify the efficacy, we compute a branching ratio using the following procedure. For each tagged rovibrational state , we integrate the measured scattering intensity over the range of angles probed and then sum the resulting integrated fluxes over the different tagged rotational states of each stretch. We obtain a measure of the total scattered flux of each stretch fundamental . From and , we obtain the branching ratio .
The branching ratios obtained from the state-resolved angular distributions are shown in Table 1, along with the branching ratio of zero implied from the results of Werdecker et al. in their fixed scattering angle measurements for the graphene-covered nickel surface. Tabulated along with the branching ratios are the associated barrier heights for dissociative chemisorption of CH4 as determined by density functional theory (DFT) calculations. We note that the barrier height quoted for Gr/Ni(111) was in fact computed for free-standing graphene, although the addition of a Ni(111) substrate is not expected to alter this figure significantly1. Taking a low barrier height as an indication of high catalytic activity, it is clear that the order of the branching ratios and barrier heights for these three surfaces agrees with the stated hypothesis.
TABLE 1
| Surface | (%) | Barrier height (eV) |
|---|---|---|
| Ni(111) | 1.110 | |
| Au(111) | 1.99 | |
| Gr/Ni(111) | 0.0† | 4.1‡11 |
Scattering angle-integrated branching ratios computed from the state-resolved angular distributions measured in this study (Figure 2) along with theoretically calculated barrier heights for methane dissociation with associated references.
†From reference [Werdecker et al., 2020]. Measured under conditions , , , and fixed .
‡Calculated for free-standing graphene.
To assess the repeatability of the branching ratio measurements the full set of measurements for the Au(111) surface was performed twice. The data from the second run showed a somewhat better signal-to-noise ratio and so we present the branching ratio of 6.6% obtained from this run in Table 1. The first run yielded a branching ratio of 7.9%. The percentage deviation of between runs is of the same order of magnitude as those measured for other quantities derived from the same runs (e.g., the mean rotational energies and presented in Section 3.3). As there is a lack of necessary data for a more rigorous statistical analysis, we take as a reasonable estimate of our branching ratio uncertainty for the Ni(111) surface as well as the Au(111) surface, with the experimental procedures and data quality for the two surfaces being largely identical. From our own analysis of the signal-to-noise ratio of the measurements presented in Werdecker et al. on Gr/Ni(111), we arrive at a rough estimate of the upper limit of the branching ratio for this surface of 2%. We therefore consider the branching ratio order presented in Table 1 to be reliable.
3.2 Trends in state-resolved angular distributions
In computing the branching ratios, we integrate over the scattering angle and sum over rotational states and thus discard any dynamical information contained in the angular and rotational state dependence of the scattering intensities. To explore in finer detail the dynamics of the molecule–surface collision, we return in this section to the state-resolved angular distributions. In the following section (Section 3.3), we narrow our attention to the gross rotational dynamics via the scattering angle-integrated state fluxes .
We have found it useful to summarize individual state-resolved angular distributions by three parameters, namely, , , and , which characterize a distribution’s width, angle of peak intensity, and peak amplitude, respectively. The parameters are determined by performing a non-linear least squares fit of the measured data using a fit function of the following form:
That the chosen fit function gives a good description of our observations can be verified by inspection of Figure 2 which shows measured data overlaid with their associated best-fit curves (see also the Supplementary Materials for full set of angular distributions and fits).
Figure 3 shows the best-fit parameter values , , and for all 28 state-resolved measurements (although we include in Figure 3 the peak amplitude parameter for completeness, we omit discussion of this parameter, deferring discussion of the related scattering angle-integrated state populations for Section 3.3). From the figure, one can readily deduce that for both surfaces and for all lines tagged the scattering mechanism appears to be “direct” or “impulsive” in that the measured distributions strongly deviate from the behavior expected of a “trapping/desorption” scattering process characterized by an extended molecular residence time at the surface (Rettner and Auerbach, 1994). Such a lack of equilibration between the surface and the scattering methane molecule has been seen before both in experiment and calculations for Ni(111) (Milot et al., 2001; Al Taleb and Farías, 2017) and Pt (111) (Kondo et al., 2018).
FIGURE 3
Fitting parameters obtained from fitting Eq. 1 to the measured state-resolved in-plane angular distributions . The rovibrational level prepared by the pump laser is indicated by the overbar on the label. Tagged levels are grouped by the total angular momentum quantum number and displaced slightly left, right, or not at all depending on whether their associated Coriolis stack is , , or , respectively. Marker symbols denote the wavefunction symmetry. The red curve in the plot of for the levels is the predicted for a rigid methane molecule with scattering from a flat and rigid surface. The two black dotted horizontal lines in the plots mark the angle of specular scattering.
In addition to our classification of the scattering as direct/impulsive, we also identify three distinct trends from the analysis of the angular distributions. First, the exponents (Figure 3, top row) are systematically higher for scattering from Ni(111) than those from Au(111), reflecting a lower degree of angular dispersion caused by the former. Second, for both Ni(111) and Au(111), the widths of the scattering peaks increase with the rotational quantum number , suggesting a correlation between rotational excitation and momentum transfer either between the molecule and the surface or between the parallel and perpendicular components of the molecular momentum.
Third, the peak angles (Figure 3, middle row) for scattering from Au(111) are shifted further away from surface normal compared to those for scattering from Ni(111). In a recent publication, we reported (Reilly et al., 2023) the same finding for CH4/Ni(111) scattering with CH4 prepared in the vibrational ground state and noted in addition a positive trend between and total angular momentum quantum number . From Figure 3, it is clear that this correlation holds for the vibrationally excited molecules as well, the trend being somewhat more pronounced for Au(111). Note that a positive – correlation is what one anticipates for the scattering of a rigid body from a flat and rigid surface based on simple considerations of conservation of energy and parallel momentum (see the red curve of Figure 3). One might therefore conclude that the CH4/Au (111) system more closely approaches this ideal than the CH4/Ni(111) system. However, in the limit of perfect flatness and rigidity, one also expects a unique scattering angle for each final state (i.e., . On this account, the broader peaks (characterized by smaller exponents ) observed for scattering from Au(111) would then imply the opposite conclusion. Our results do not therefore admit a conclusive determination of which surface appears more “mirror-like” to a scattering methane molecule. Of course, a vibrating molecule cannot in the first place be said to be a rigid body, and certainly for the scattering channel, the situation is further complicated by the change in the vibrational mode structure and accompanying energy release.
Going beyond the trends just discussed, we can also make a few additional qualitative observations. We first point out the anomalously large exponent for the elastic channel on Ni(111) is 38% larger than the next-largest channel exponent and 61% larger than the largest exponent for the same surface. In a practical sense, this fact, combined with the fact that the elastic channel on Ni(111) is relatively intense, establishes that the angular distributions for the vibrationally elastic and inelastic channels can in fact differ and that it is, therefore, essential to measure overall scattered angles to obtain the most accurate estimate of the total scattering probability. Indeed, when we compute “fixed angle branching ratios” from our Ni(111) measurements, we find that their values vary by more than 50% over the range of scattering angles probed (Supplementary Section S2).
In another sense, that the vibrationally and rotationally elastic scattering peak is so much narrower than the others and should imply something about the molecule–surface interaction potential and associated scattering dynamics. We offer the following simple microscopic explanation for a strong and anomalously narrow elastic scattering distribution peaking near specular (). Of the full range of values of impact parameters characterizing the possible initial conditions of the scattering molecule (e.g., impact site, initial molecular orientation, and internal/external axis of rotation), there is a significant fraction which results in trajectories where the molecule is turned around at distances relatively far from the surface. At these distances, the interaction between the molecule and the surface manifests primarily as a mutually repulsive force directed along the surface normal. That is, for these trajectories, the molecules do not penetrate far enough for the molecule–surface interaction to either exhibit significant lateral surface corrugation or effect significant disruption of the molecule’s or surface’s internal structure. Thus, the various mechanisms available to induce a change in the rovibrational state or deflect the scattering molecule away from specular scattering (parallel-to-perpendicular momentum transfer, surface phonon creation/annihilation, and rovibrational inelasticity) would act only weakly, leading to a narrow elastic distribution peaking near specular scattering. On the other hand, trajectories permitting deep enough penetration for the interaction to effect change in the rovibrational state will also tend to disperse angularly due to the additional effects of lateral corrugation and phonon creation/annihilation, resulting in broader scattering distributions for the rovibrationally inelastic channels. Interestingly, the strength of these dispersive forces seems to be no greater for the trajectories undergoing vibrational state change than those simply changing rotational state.
To qualitatively summarize the explanation just presented, our observations can be interpreted as implying a certain order in which the different aspects of the molecule–surface interaction “turn on” as the molecule approaches the surface, with a two-body repulsive force preceding the onset of dispersive forces, which in turn precedes the onset of forces acting to distort or rotate the molecule.
We close the section with a discussion regarding the differences in peak widths and peak angles , particularly for the Ni(111) data, among the levels of equal J in the stretch. These levels have the same nominal rotational energy (within 20) but are distinguished by their associated Coriolis stack and wavefunction symmetry . The Coriolis stack of a level indicates, loosely speaking, the relative orientation of its vibrational angular momentum and its total angular momentum (parallel, perpendicular, or antiparallel for , , and , respectively) (di Lauro, 2020). Centrifugal coupling between rotation and vibration splits a Coriolis level ( worst case (Hecht, 1961a)) into “sublevels” of definite symmetry with respect to the point group of the CH4 molecules (Hecht, 1961b), of which those of and symmetries belong to the meta-CH4 nuclear spin isomer. The sublevels of a Coriolis stack are distinguished by the internal alignment of their total angular momentum, i.e., by the different molecule-fixed axes about which the molecule is more and less likely to be found rotating (again, loosely speaking).
Compare for example the best-fit and values for the () and () levels scattering from Ni(111). The peak of the level has a 25% larger best-fit exponent and a peak angle closer to normal than that of the level. Similar differences occur between the and sublevels in the same ( Coriolis stack. These rather subtle details of the internal molecular motion have a significant impact on the scattering distributions and highlight the “dynamic” nature of the CH4/Ni(111) interaction in that the form of the state-resolved distributions are not solely dictated by the state’s energy. In a recent publication published by the authors on surface scattering of CH4 in the vibrational ground state (Reilly et al., 2023), we found strong evidence for a more extensive molecule–surface collisional energy exchange (i.e., equilibration) after oxidation of a Ni(111) surface than before oxidation. We also found that the forms (and amplitudes) of the sublevel angular distributions (the ground vibrational state has only Coriolis stacks) produced upon scattering from the oxidized surface were essentially indistinguishable, while scattering from the unoxidized surface produced differences in the peak widths and peak angles for the different sublevels comparable to those observed in the vibrationally excited state. We argued that the lack of imbalances observed in the collisions with the oxidized surface was consistent with the other reported observations in that they reflected an equilibration among the sublevels caused by a prolonged molecule–surface energetic coupling acting to scramble any delicate internal rotational motion.
Concerning the Au(111) surface data, while the form of the distributions among levels of equal is similar, we will see in the following section (Section 3.3) that the level populations (related indirectly to the peak amplitudes ) differ dramatically depending on the internal alignment of the molecule’s total and vibrational angular momentum, so that the scattering for this system is, in this sense, no less “dynamic.”
3.3 Rotational dynamics
Figure 4 shows the scattering angle-integrated rovibrational state populations . In addition to illustrating what is already evident from Table 1—namely, that the overall conversion is much weaker for the Au(111) surface—the plotted distributions also reveal important additional features of the rotational dynamics of the molecular–surface collisions for the systems studied here.
FIGURE 4
Distribution of scattering angle-integrated rovibrational state fluxes , where . The rovibrational level prepared by the pump laser is indicated by the overbar on the label. The Coriolis stack (, , and ) of a tagged level is indicated by the hatching pattern of its associated bar and the wavefunction symmetry ( and ) by the bar color. All states are of the type. The weakly populated (, , ) level for the Au(111) data is of the type (same as for Ni(111)).
One feature that is directly observable is the differences in the forms of the distributions of the vibrationally elastic () channels for the two surfaces, with the Au(111) distribution having a relatively narrow distribution concentrated at low , suggesting a “colder” rotational population. To quantify the degree of rotational excitation in the different stretch fundamentals, we use two different measures. The first measure is the mean rotational energy obtained by taking the -weighted average of the rotational energies of the different tagged levels of the associated stretch . The rotational energy of a level is calculated by taking the energy difference between the level and the rovibrational ground state as determined by high resolution spectroscopy (Gordon et al., 2022) and subtracting the theoretically determined energy gap separating the vibrationally excited state from the ground vibrational state of a hypothetical non-rotating methane molecule (Boudon et al., 2006). The second measure is the effective rotational temperature obtained by fitting the of the associated stretch to a thermal (Boltzmann) distribution of the form , where is the Boltzmann constant and is the level degeneracy.
The extracted mean rotational energies and effective rotational temperatures (Figure 5 for fits) are presented in Table 2. The differences in the distributions mentioned previously are reflected in the mean energies, with the value for Ni(111) being 39% larger than that of Au(111). The difference in effective temperatures is less pronounced, due to the greater influence of the enhanced level population scattering from the Au(111) surface (level populations contribute linearly to , while the fitting procedure minimizes the sum of the squared logarithmic population deviations). That the effective rotational temperatures for all four combinations of surfaces and stretching modes lie well below the surface temperature of and implies a lack of equilibration between molecular rotation and the surface degrees of freedom, offering further evidence against a trapping-desorption scattering mechanism.
FIGURE 5
Boltzmann fits to the rotational-level population. Plotted are the scattering angle-integrated state-resolved fluxes weighted by the level degeneracy vs. level energy (within a constant -dependent shift ). The fluxes for the Au(111) surface data have been scaled down by a factor 20 to avoid overlap with the Ni(111) data. Solid (dashed) lines indicate fits of the () data to the Boltzmann form . Best-fit effective rotational temperatures are printed next to the associated best-fit line. The rovibrational level prepared by the pump laser corresponds to the left-most solid markers. The levels of a given stretch and total angular momentum are nominally isoenergetic (), with the energy ordering for the different Coriolis stacks of a given combination being . For both surfaces, the sublevel of the Coriolis stack is more populated than the sublevel.
TABLE 2
| Surface | [cm-1] | [cm-1] | [K] | [K] |
|---|---|---|---|---|
| Ni(111) | 69.8 | |||
| Au(111) |
Extracted mean rotational energies and effective rotational temperatures from the rotational distributions presented in Figure 4. Rotational distributions were measured twice for Au(111), and we report the values of and obtained from the run with the higher signal-to-noise ratio. The uncertainty of reported for the Au(111) values was obtained by taking the larger of the two deviations in measured for the two runs. The reported uncertainties for are the standard errors reported by the least-squares fitting procedure implemented in the optimize.curve_fitfunction of the SciPyPython library (Virtanen et al., 2020).
Perhaps the most interesting result emerging from this analysis is that the ratio of the degree of rotational excitation for molecules scattered into the and states differs dramatically for scattering from Au(111) versus Ni(111). The ratio of the mean rotational energies is 71% higher for Au(111) than for Ni(111). Similarly, the ratio is 72% higher for Au(111) than for Ni(111).
Why do the ratios and differ so dramatically for the two surfaces? We begin by pointing out that the fundamental lies in energy above the fundamental (Boudon et al., 2006). A more general question thus poses itself: where does this extra vibrational energy go during an SIVR process? Since we cannot measure the speed of the scattered molecules or directly detect energy transfer to the surface, a complete experimental accounting of the energy balance is not possible. Nonetheless, our observations indicate that a greater portion of this excess energy is transferred to rotation in a collision with Au(111) ) than is transferred in a collision with Ni(111) ().
In the reduced-dimensional CH4/Ni(111) SIVR calculations presented in Werdecker et al., the surface was held frozen and the molecule was forced to assume at all times the minimum energy orientation so that energy released in the conversion was constrained to be transferred completely into translational kinetic energy of the scattering molecule. In light of our results, it would be of keen interest to learn how this energy is partitioned in a dynamical simulation that includes both molecular rotation and surface vibrations as degrees of freedom. It would be equally interesting to learn to what extent the catalytic C-H bond elongation discussed in Introduction is also operative on Au(111). Perhaps the enhanced rotational excitation on Au(111) reflects a mechanism for SIVR on Au(111) that is fundamentally different from the catalytic effect thought to be operating on Ni(111). SIVR on Au(111) might, for example, be a purely “mechanical” effect arising from the violent recoil of the molecule experiences upon rebounding from the much more massive gold surface atoms. An influence of the surface atom mass on CH4/metal scattering inelasticity has been observed in recent dynamical calculations (Gerrits et al., 2019; Jackson, 2022) which finds Baule-like molecule–surface energy transfer in CH4/metal surface collisions. A measurement of CH4/Cu (111) SIVR efficiency might shed light on the role of a mechanical mechanism. The catalytic activity of Cu(111) with a calculated barrier height of 1.32 eV (Bhati et al., 2022) falls in between that of the Ni(111) and the Au(111) surface. In light of both this and the ability of the lighter copper surface atoms to better absorb the shock of impact, a Cu(111) branching ratio weaker than that of Au(111) would stand as strong evidence for a mechanical SIVR mechanism operating on Au(111).
The last observation we wish to make about the rotational distributions concerns the levels of equal , a subject we began discussing at the end of the previous section (Section 3.2) in the context of angular distributions. Figures 4, 5 clearly show that scattering from the Au(111) produces a highly imbalanced distribution of populations among sublevels of equal . These imbalances imply a scattered molecular flux with a high degree of both internal alignment of the total angular momentum (indicated by the imbalance of the and sublevel populations within the Coriolis stack) and relative alignment of the total and vibrational angular momentum (indicated by imbalances among the Coriolis stacks of equal for ). From the restricted subset of levels tagged in this experiment, there appears to be a strong tendency for preservation in the CH4/Au(111) scattering event of the parallel (i.e., ) alignment of total and vibrational angular momentum prepared by the pump laser. Although the population imbalances within a given are not as pronounced in the CH4/Ni(111) scattering data, it is nonetheless interesting to note that the order of the imbalances for the two surfaces are identical for the levels probed.
At the present time, the authors have no explanation for the form and extent of this internal alignment of molecular rotation observed in the scattering. We plan to present more extensive results on these interesting alignment phenomena in CH4 surface scattering in a future publication.
4 Summary and future directions
A number of promising avenues exist for further testing the connection between SIVR and catalytic activity. One extension of the work presented here would be a study of SIVR efficiency for a range of different incident kinetic energies . In one sense, one expects SIVR to correlate positively with since, with the increasing incident energy, the molecule can get closer to the surface and gain access to the catalytic forces that distort the molecule’s equilibrium geometry and mediate vibrational mode coupling. Indeed, a positive -SIVR correlation is observed in the CH4/Ni (111) dynamical calculations of Werdecker et al. discussed in Introduction. On the other hand, depending on the form of the repulsive molecule–surface interaction, a high-incident velocity may result in the brief molecule–surface interaction time, which might tend to weaken the SIVR produced on scattering.
Another direction, one which would offer a more direct validation of the theoretical results reported in Werdecker et al., would be to measure the variation in SIVR across the different microscopic impact sites of the same macroscopic surface. Although it is of course not feasible to experimentally control the impact site of a molecule, one can vary the relative concentration of different sites. In the calculations of Werdecker et al., the authors find an increase in SIVR at sites where a surface Ni atom is displaced from equilibrium. Such displacements are also expected to increase a site’s catalytic activity (Nave et al., 2010). By raising/lowering the surface temperature, one can increase/decrease the density of thermally generated lattice distortions and test for the expected increase/decrease in SIVR.
Another more controlled means of varying the density of different surface sites is to vary the exposed crystal plane, characterized by Miller indices (ijk) (Kolasinski, 2012). For fcc metals (e.g., Ni, Au, and Pt), the atoms of the (111) (or “close-packed”) surface are highly coordinated and therefore present sites of relatively low catalytic activity. The stepped fcc(211) surface, however, exposes under-coordinated atoms at step sites, offering lower barrier sites for methane dissociative chemisorption, as observed in experiments and calculations for CH4 sticking experiments on platinum (Chadwick et al., 2018). A measurement of the branching ratios for different surface planes of the same crystal could thus serve as a compelling experimental verification of the theoretically observed connection between a site’s SIVR efficiency and its catalytic activity.
We conclude with a brief summary. For the Ni(111) and Au(111) surfaces, we measured rovibrational state-resolved angular distributions for 10 (4) rotational levels of the () stretch fundamentals produced by scattering of incident CH4() molecules. From these distributions, we extracted for each surface a branching ratio characterizing the efficiency of collision-induced conversion. The branching ratio measured for the Au(111) surface is found to be more than five times smaller than that measured for the Ni(111) surface but significantly larger than that measured by Werdecker et al. at a fixed scattering angle for the Gr/Ni(111) surface. The ordering of the branching ratios for the three surfaces is found to match the ordering of their catalytic activities, which is in accordance with the hypothesis advanced in Werdecker et al.
Analysis of the shape of angular distributions reveals for both surfaces a broadening of the scattering distributions with increasing total angular momentum quantum number , with the broadening stronger for scattering from the Au(111) surface. In addition, we observe a tendency for super-specular scattering, and this tendency increases with and is more pronounced for the Au(111) surface. For the Ni(111) surface, the elastic channel is strong and anomalously narrow. The rovibrational levels of equal in the stretch show measurable differences in width and peak angle in scattering from the Ni(111) surface.
From the distribution of scattering angle-integrated state fluxes, the Ni(111) surface produces more highly rotationally excited molecules in the stretch than the Au(111) surface. Molecules in the stretch are, however, found to be dramatically more rotationally excited in scattering from Au(111) than from Ni(111), prompting a questioning into the nature of the flow of vibrational energy for the two surfaces. The rovibrational levels of equal in the stretch have strongly imbalanced populations upon scattering from Au(111), suggesting strong alignment of the internal angular momentum in the scattered flux and a preservation of the parallel Coriolis coupling of the total and vibrational angular momentum prepared in the incident molecules.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
RB conceptualized the project, secured funding, and supervised the project. DA designed the experimental setup with support from RB. PF and CR performed the experiments, and PF analyzed the data. CR wrote the manuscript with contributions from PF, DA, and RB. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors acknowledge financial support by the Swiss National Science Foundation through Grant No. 200497 and by the Max Planck-EPFL Center for Molecular Nanoscience and Technology. The authors thank the EPFL ISIC machine shop and Bo-Jung Chen for their invaluable contributions to the construction of the apparatus used to perform the experiments presented in this manuscript. The authors would also like to thank Bret Jackson for his calculations of the CH4/Au(111) dissociation barrier and valuable discussions. DA gratefully acknowledges support from the International Center for Advanced Energy Conversion in Göttingen, Germany.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2023.1238711/full#supplementary-material
Footnotes
1.^Jackson, B. (2023) Private communication.
References
1
Al Taleb A. Farías D. (2017). Coherent quantum scattering of CH 4 from Ni(111). Phys. Chem. Chem. Phys.19, 21267–21271. 10.1039/c7cp04559j
2
Auerbach D. J. Tully J. C. Wodtke A. M. (2021). Chemical dynamics from the gas-phase to surfaces. Nat. Sci.1, e10005. 10.1002/ntls.10005
3
Bhati M. Dhumal J. Joshi K. (2022). Lowering the C–H bond activation barrier of methane by means of SAC@Cu(111): Periodic DFT investigations. New J. Chem.46, 70–74. 10.1039/d1nj04525c
4
Boudon V. Rey M. Loëte M. (2006). The vibrational levels of methane obtained from analyses of high-resolution spectra. J. Quant. Spectrosc. Radiat. Transf.98, 394–404. 10.1016/j.jqsrt.2005.06.003
5
Brown L. R. Margolis J. Champion J. Hilico J. Jouvard J. Loete M. et al (1992). Methane and its isotopes: Current status and prospects for improvement. J. Quant. Spectrosc. Radiat. Transf.48, 617–628. 10.1016/0022-4073(92)90126-o
6
Chadwick H. Guo H. Gutiérrez-González A. Menzel J. P. Jackson B. Beck R. D. (2018). Methane dissociation on the steps and terraces of Pt(211) resolved by quantum state and impact site. J. Chem. Phys.148, 014701. 10.1063/1.5008567
7
Chadwick H. Hundt P. M. van Reijzen M. E. Yoder B. L. Beck R. D. (2014). Quantum state specific reactant preparation in a molecular beam by rapid adiabatic passage. J. Chem. Phys.140, 034321. 10.1063/1.4861054
8
di Lauro C. (2020). “Spectra of spherical top molecules,” in Rotational structure in molecular infrared spectra (China: Elsevier), 223–242.
9
Gerrits N. Chadwick H. Kroes G.-J. (2019). Dynamical study of the dissociative chemisorption of CHD3 on Pd(111). J. Phys. Chem. C123, 24013–24023. 10.1021/acs.jpcc.9b05757
10
Gordon I. E. Rothman L. Hargreaves R. Hashemi R. Karlovets E. Skinner F. et al (2022). The HITRAN2020 molecular spectroscopic database. J. Quant. Spectrosc. Radiat. Transf.277, 107949. 10.1016/j.jqsrt.2021.107949
11
Hecht K. T. (1961b). The vibration-rotation energies of tetrahedral XY4 molecules: Part I. Theory of spherical top molecules. J. Mol. Spectrosc.5, 355–389. 10.1016/0022-2852(61)90102-3
12
Hecht K. T. (1961a). Vibration-rotation energies of tetrahedral XY4 molecules: Part II. The fundamental ν3 of CH4. J. Mol. Spectrosc.5, 390–404. 10.1016/0022-2852(61)90103-5
13
Jackson B. Nave S. (2013). The dissociative chemisorption of methane on Ni(111): The effects of molecular vibration and lattice motion. J. Chem. Phys.138, 174705. 10.1063/1.4802008
14
Jackson B. (2022). Quantum studies of methane-metal inelastic diffraction and trapping: The variation with molecular orientation and phonon coupling. Chem. Phys.559, 111516. 10.1016/j.chemphys.2022.111516
15
Kolasinski K. (2012). “Surface and adsorbate structure,” in Surface science (China: John Wiley and Sons, Ltd), 9–49. 10.1002/9781119941798.ch1
16
Kondo T. Al Taleb A. Anemone G. Farías D. (2018). Low-energy methane scattering from Pt(111). J. Chem. Phys.149, 084703. 10.1063/1.5044744
17
Li K. Li H. Yan N. Wang T. Zhao Z. (2018). Adsorption and dissociation of CH4 on graphene: A density functional theory study. Appl. Surf. Sci.459, 693–699. 10.1016/j.apsusc.2018.08.084
18
Milot R. Kleyn A. W. Jansen A. P. J. (2001). Energy dissipation and scattering angle distribution analysis of the classical trajectory calculations of methane scattering from a Ni(111) surface. J. Chem. Phys.115, 3888–3894. 10.1063/1.1388224
19
Nave S. Tiwari A. K. Jackson B. (2010). Methane dissociation and adsorption on Ni(111), Pt(111), Ni(100), Pt(100), and Pt(110)-(1 x 2): Energetic study. J. Chem. Phys.132, 054705. 10.1063/1.3297885
20
Reilly C. S. Floß P. Chen B.-J. Auerbach D. J. Beck R. D. (2023). Quantum state-resolved methane scattering from Ni(111) and NiO(111) by bolometer infrared laser tagging: The effect of surface oxidation. J. Chem. Phys.158, 214202. 10.1063/5.0150009
21
Rettner C. T. Auerbach D. J. (1994). Distinguishing the direct and indirect products of A gas-surface reaction. Science263, 365–367. 10.1126/science.263.5145.365
22
Shenvi N. Roy S. Tully J. C. (2009). Nonadiabatic dynamics at metal surfaces: Independent-electron surface hopping. J. Chem. Phys.130, 174107. 10.1063/1.3125436
23
Virtanen P. Gommers R. Oliphant T. E. Haberland M. Reddy T. Cournapeau D. et al (2020). SciPy 1.0: Fundamental algorithms for scientific computing in Python. Nat. Methods17, 261–272. 10.1038/s41592-019-0686-2
24
Wagner V. Park G. B. Wallrabe M. Wodtke A. M. Schäfer T. (2019). Electron transfer mediates vibrational relaxation of CO in collisions with Ag(111). Phys. Chem. Chem. Phys.21, 1650–1655. 10.1039/c8cp06041j
25
Werdecker J. Chen B. J. Van Reijzen M. E. Farjamnia A. Jackson B. Beck R. D. (2020). State-to-state methane-surface scattering as a probe of catalytic activity. Phys. Rev. Res.2, 043251. 10.1103/physrevresearch.2.043251
26
Werdecker J. van Reijzen M. E. Chen B.-J. Beck R. D. (2018). Vibrational energy redistribution in a gas-surface encounter: State-to-State scattering of CH 4 from Ni(111). Phys. Rev. Lett.120, 053402. 10.1103/physrevlett.120.053402
27
Wodtke A. M. Matsiev D. Auerbach D. J. (2008). Energy transfer and chemical dynamics at solid surfaces: The special role of charge transfer. Prog. Surf. Sci.83, 167–214. 10.1016/j.progsurf.2008.02.001
28
Zhang L. Jiang B. (2019). Efficient vibrational energy redistribution between stretching modes: State-to-State quantum scattering of H 2 O from Cu(111). Phys. Rev. Lett.123, 106001. 10.1103/physrevlett.123.106001
29
Zhang L. Jiang B. (2020). State-to-state quantum dynamics of H 2 O/HOD scattering from Cu(111): Mode- and bond-selective vibrational energy transfer. J. Chem. Phys.153, 214702. 10.1063/5.0030490
Summary
Keywords
methane dissociation, state-to-state scattering, angular distributions, surface-induced vibrational energy redistribution, optothermal spectroscopy, bolometer infrared laser tagging, heterogeneous catalysis
Citation
Floß P, Reilly CS, Auerbach DJ and Beck RD (2023) Surface-induced vibrational energy redistribution in methane/surface scattering depends on catalytic activity. Front. Chem. 11:1238711. doi: 10.3389/fchem.2023.1238711
Received
12 June 2023
Accepted
05 July 2023
Published
25 July 2023
Volume
11 - 2023
Edited by
Marco Sacchi, University of Surrey, United Kingdom
Reviewed by
Bret Jackson, University of Massachusetts Amherst, United States
Hirokazu Ueta, Japan Atomic Energy Agency, Japan
Updates
Copyright
© 2023 Floß, Reilly, Auerbach and Beck.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rainer D. Beck, rainer.beck@epfl.ch
†These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.