Measuring transient reaction rates from non-stationary catalysts

Up to now, the methods available for measuring the rate constants of reactions taking place on heterogeneous catalysts require that the catalyst be stable over long measurement times. But catalyst are often non-stationary, they may become activated under reaction conditions or become poisoned through use. It is therefore desirable to develop methods with high data acquisition rates for kinetics, so that transient rates can be measured on non-stationary catalysts. In this work, we present velocity resolved kinetics using high repetition rate pulsed laser ionization and high-speed ion imaging detection. The reaction is initiated by molecular beam pulses incident at the surface and the product formation rate is observed by a sequence of laser pulses at a high repetition rate. Ion imaging provides the desorbing product flux (reaction rate) as a function of reaction time for each laser pulse. We demonstrate the method using a 10 Hz pulsed CO molecular beam pulse train to initiate CO desorption from Pd(332) - desorbing CO is detected every millisecond by non-resonant multiphoton ionization using a 1-kHz Ti:Sapphire laser. This approach overcomes the time-consuming scanning of the delay between CO and laser pulses needed in past experiments and delivers a data acquisition rate that is 10-1000 times higher. We also apply this method to CO oxidation on Pd(332) - we record kinetic traces of CO$_2$ formation while a CO beam titrates oxygen atoms from an O-saturated surface. This provides the reaction rate as a function of O-coverage in a single experiment. We exploit this to produce controlled yet inhomogeneously mixed reactant samples for measurements of reaction rates under diffusion-controlled conditions.


Introduction
Methods to measure the kinetics of surface reactions are crucial to improving our understanding of heterogeneous catalysis. Traditionally, temperature programmed reaction, molecular beam relaxation spectrometry and phase-lag detection have been available to experimentalists 1-4 . Recently, the kinetic trace was obtained using velocity resolved methods 5 based on ion imaging [6][7][8] . This is essentially a pump-probe technique where a molecular beam pump-pulse initiates the reaction and pulsed laser ionization probes the desorbing products. Varying the delay between the two pulses provides the time base of the reaction kinetics. The ionized products are recorded with ion imaging providing product velocity information with every detection pulse. This allows measured product densities to be converted to product flux, which is by definition the reaction rate for a surface reaction. Furthermore, flight-times irrelevant to the reaction time can be subtracted from the experimental time-axis 9 . Like all pump-probe measurements, during the time that the delay between pump and probe is being scanned, the reacting system under study must not change.
However, catalysts are often dynamic. Catalyst composition can change dramatically under reactive conditions 10 -living catalyst [11][12][13] -and catalytic use can lead to poisoning 14 . For such systems, we need methods that can rapidly obtain kinetic information while the catalyst is changing.
In this work, we demonstrate velocity resolved kinetics with high repetition rate detection. The reaction starts when a pulse of molecules arrives at the surface and ion images are recorded for each pulse of a high repetition-rate laser that ionizes desorbing products. The ion images preserve the velocity information from which the rate of reaction is derived. The inverse repetition rate of the laser sets the temporal resolution.
We demonstrate a duty cycle that is one to three orders of magnitude higher than previous methods 5 , allowing measurements on a changing catalyst. The present experiments use a 1-kHz Ti:Sapphire laserfuture experiments with Yb-fiber lasers operating at 10 kHz provide a perspective for improvement.

Experiment
We previously described the apparatus in detail elsewhere [6][7][8] . Briefly, we produce two molecular beams in two vacuum chambers, each equipped with piezo-electrically actuated pulsed valves. The valves' repetition rates are variable up to 500 Hz. The pulse durations can be as low as 30 μs. Each beam passes through two differential pumping chambers, before entering an ultrahigh vacuum (UHV) chamber with a base pressure of 2×10 −10 mbar, where they intersect one another and collide with a Pd(332) surface. One beam collides at normal incidence dosing the sample with oxygen. The second beam, incident at 30° to the normal, initiates the reaction with a pulse of CO. The CO beam can either be used alone to study CO trapping/desorption or with an oxidized surface to initiate CO2 formation. A single crystal of Pd cut and polished to expose the (332) surface is mounted on a 5-axis manipulator and can be heated to 1150 K using electron bombardment.
The instrument is equipped with an Ar + sputtering source for cleaning the surface as well as an Auger electron spectrometer to check its cleanliness.
A homogeneous electric field oriented parallel to the surface is formed by two parallel flat meshes (repeller and extractor), between which both molecular beams pass. After ionization of the reaction products by a non-resonant multiphoton process, using an ultra-short Ti:Sapphire laser (Coherent Astrella, 800 nm, 35 fs, 0.5 mJ, 1 kHz) focused with a 150 mm plano-convex lens, a 3 kV pulse applied to the repeller of the ion imaging system directs the ions to the imaging detector. This maps the products' density and in-plane velocity vectors, which is used to create a flux image. A region of the flux image is then integrated to provide the rate of reaction at a specific time. We record ion images with a 56 mm Chevron MCP detector coupled to a P43-phosphor screen, whose phosphorescence detected by a high-frame rate CMOS camera (Vision Research Phantom VEO 710). We took advantage of commercial data acquisition software (DaVis LaVision GmbH) and a software-controlled timing unit (PTUX, LaVision GmbH). The timing unit is triggered both at 10 Hz-synchronized with the pulsed nozzle-and at 1−kHz-synchronized with the laser. Several thousand images are recorded over several seconds and stored on the camera's internal memory, only to be transferred later to a computer's hard disk. Figure. 1 shows a comparison to methods requiring the delay between pulsed molecular beam and laser, , to be scanned. In that case (Fig. 1a), an ion image is measured for a fixed value and an ion image is accumulated over many (typically 50) molecular beam pulses. is then incremented and the process is repeated. Here, one ion image is recorded for every molecular beam pulse, whose repetition rate is typically 10 to 100 Hz. Using high repetition-rate detection, the ion image is recorded every millisecond. Each pulse of the laser (points in Fig. 1b) corresponds to a point in the temporal evolution of the reaction. The P43 phosphor screen decays over %→ % = 1.3 ms, while the time between laser pulses is only 1 ms. Hence, after downloading the image sequence to the computer we subtracted from each image the "afterglow background" remaining from the previous image.

Proof of principle: application to CO desorption from Pd(332)
As a proof of principle, we performed measurements on CO trapping/desorption from Pd(332) between 583-623 K. Here, a clean Pd(332) crystal is exposed to a pulsed molecular beam of pure CO operating at 10 Hz. The surface temperature is controlled so that a 1−kHz detection rate is sufficiently rapid to follow the desorption kinetics, while also ensuring that all CO molecules desorb between molecular beam pulses.
Following Ref.'s 6-7 , we extract the kinetic trace by integrating flux images between 300 and 900 m/s and ±4° form the surface normal. This captures most of the desorbing molecules, while suppressing signal from directly scattered (higher velocity) and background (lower velocity) CO. Figure 2 shows data from a typical 5-second experiment, requiring 1% the measurement time needed for delay-scanning. 50 kinetic traces result, one from each of 50 CO molecular beam pulses. The inset shows three kinetic traces in detail. We filter to the raw data (blue) with a periodic Savitzky-Golay filter 15 applied by first sorting the data according to and then employing a moving linear fit to a single data point and 10 of its neighboring data points-all with the same . The value of the fitted line then replaces the data point and the process is repeated on the next data point. This leads to the filtered output (black). The CO desorption rate constant, , is determined by fitting each pulsed decay with a function that convolves the incident CO beam's temporal profile with an exponential decay-red line in Fig. 2 inset decay 5 . In this way, we derive fifty independent values of , from which we obtain an average value and a standard deviation.

Duty-cycle analysis
We consider now the quantitative duty-cycle improvements possible with high rep-rate detection, within the specific context of desorption rates near zero-coverage. We first define a characteristic desorption time, , which is the inverse of the desorption rate constant, = . This imposes an upper limit of the molecular beam's repetition rate (' ( ) * ) and therefore a minimum repeat-time, )+, = 1 ' ( ) * ⁄ , needed to maintain the low coverage condition. While there is some ambiguity involved, we set )+, = 5 , the time at which a 1 st -order decay has reached 0.7% of its initial value. Data obtained within )+, are most important to the fitting-we label this data "relevant".
The number of relevant data obtained from each molecular beam pulse used in the high rep-rate approach, / 011 is given by: where ' is the detection laser repetition rate and data acquisition rate, /2 011 , is given by: where ' ( is the repetition rate of the pulsed molecular beam.
The number of relevant data per molecular beam pulse in a conventional delay-scanning experiments, / 34 , is 1 and the data acquisition rate is then: Taking the ratio of these two data acquisition rates, we find that the theoretical improvement in duty cycle is given by: The data acquisition rate using the delay-scanning approach is limited by )+, and, of course, technical limitations to the rep-rate of pulse beams (in our experience ~500 Hz), whereas ' is the only limiting factor to the data acquisition rate for the high rep-rate method. We emphasize that ' can be improved dramatically. This work used a Ti:Sapphire laser, ' = 1 kHz; newly available Yb-fiber lasers achieve repetition rates of 10 − 10 kHz, while still providing pulse energies and peak intensities sufficient for non-resonant multiphoton ionization.
The analysis so far neglects the number of ions produced in each experiment, which is equally important as the rate of data acquisition. All velocity resolved kinetics signals are proportional to the rate of product formation 6 . Hence in the desorption experiments presented above, the number of ions detected per laser pulse is proportional to 1/ . The dependence on reflects the temporal dilution seen for slow reactions.
Each molecular beam pulse deposits the same number of CO molecules on the surface; so, the observed density is diluted greatly over time for slow reactions and less so for fast reactions. Taking this into account, we may define the "count acquisition rate" (CAR).
This quantity determines the signal-to-noise (S/N) of the data obtained in any experiment. These equations point out that experiments using delay-scans exhibit decreasing S/N as increases, but this is not the case for high rep-rate measurements. This of course, mirrors the implications of Eq. 4. This also means that comparing different data acquisition methods should be done as a function of . Figure 4 shows calculated values of CAR vs. for a few different experimental configurations. Here, we only consider values larger than the shortest molecular beam pulse, which defines the kinetic resolution (black vertical line). To ensure that the CAR results only from relevant data, the molecular beam repetition rate should be matched to )+, = 5 = 1 ' ( ⁄ . This is true for either delay-scanning or high rep-rate detection. This gives rise to CAR plots for optimized delay-scanning (blue dashed line) and optimized 1−kHz detection (blue solid line) in Fig. 4. The red solid line shows CAR when using optimized 100−kHz detection. We also show in Fig. 4 Fig. 4 through the slope of CAR vs. for delay-scan measurements, which is steeper than that of high rep-rate experiments. Furthermore, increasing ' further increases CAR. This shows that the high rep-rate method becomes extremely attractive for measuring slow rates. From our experience, the feasibility limit in an optimized delay scanning experiment is reached for ~10-40 ms. Delay scanning measurements under these conditions take on the order of 1h. The same limit is reached in a 1−kHz measurement when ~5 s, which can be extended to 500 s with 100−kHz detection. This shows that the high rep-rate detection approach can be applied to measure values over ~7 orders of magnitude, whereas, delay scanning is limited to at most 3 orders of magnitude. High rep-rate detection thus enables measurements over a wider temperature range, providing more accurate Arrhenius parameters and greater sensitivity to non-Arrhenius behavior.
We also compare experimentally observed CARs obtained from our actual CO desorption experiments. The vertical gray dash-dotted line of Fig. 4 (marked with @ 593 K) represents the temperature at which the CO desorption experiments presented in Fig. 2 were carried out. Here, delay-scanning required 20 min to obtain ~250 relevant data, while 1−kHz detection provided ~70 relevant data in 10s. The derived rate constants were of similar accuracy for both methods. Normalizing to the number of relevant data points obtained, we find that CAR increased by a factor of ~30 for 1−kHz detection compared to delay scanning.
Seen at the -value at 593 K, the theoretical CAR plots (magenta line) and (green dashed line) show a theoretical enhancement factor that is close the observed enhancement.

A real time titration experiment for CO oxidation at Pd(332)
The velocity resolved kinetics experiment carried out with delay scanning provides time-resolved information by recording a signal arising from two pulses with a variable delay. Such experiments require coverage on the surface. We compare the integral from oxygen saturated surfaces with those obtained from a steady-state oxygen covered surface to determine the fraction of the total oxygen coverage that remains under steady state conditions. Clearly, this procedure is not optimal; ideally, one would like to know the kinetic trace at each point in the titration. While this is tremendously tedious and time consuming to perform with delay scanning, it is easily achieved with high rep-rate detection. Figure 5 shows such a measurement carried out on Pd(332) at H 4 = 503 K. Here, we first saturated the surface with oxygen by dosing with a 500 Hz O2 molecular beam pulse for 5 minutes (total exposure 300 ± 80 ML). The high rep-rate raw data (blue lines) results from a CO pulsed beam operating at 50 Hz. With each CO titrant pulse, a certain amount of oxygen is removed from the surface, so that each kinetic trace probes a different O-atom surface coverage. Using a 51 point periodic Savitzky-Golay filter as described above we filtered the raw data (blue lines of Fig. 5) to yield the filtered data (black lines of Fig. 5). Two insets in Fig. 5 show representative kinetic traces at early and late times in the titration. We find that the signal amplitude decreases and the rate slows with increasing titration time, reflecting the consumption of oxygen with each subsequent CO pulse. The filtered data can be represented by a first order decay for the entire 30 s titration time and for temperatures between 473 and 533 K.

Reaction rate analysis at high oxygen coverages
This approach provides new information about the nature of the kinetics and improves the performance of the velocity resolved kinetics methods. One advantage is the ability to obtain rates of reactions at saturated oxygen coverage, where the absolute oxygen coverage is unambiguously defined. To demonstrate this, we apply a simple model previously suggested by Engel and Ertl to describe CO oxidation kinetics on Pd 18 .
The model incorporates four processes: and the CO2 formation rate is given by: We take advantage of the fact that the velocity resolved kinetics signal is directly proportional to the rate of CO ,M formation and that during the first few CO pulses arriving at the surface, we probe a well-defined oxygen coverage [O * ] F G = 0.292 ML 6 . We present examples of such kinetic traces in Fig. 6 Fig. 7 as an Arrhenius plot. Note that in Fig. 6 for H 4 = 533 K, the temporal resolution is insufficient to provide a reliable fit. This problem could be solved by repetitive measurements of this type at a variety of delays of the CO molecular beam pulse -says interleaved in 0.1 ms steps. One could also use a higher repetition rate laser. Alternatively, ^ can be obtained simply from the amplitude of the kinetic trace, shown as arrows in Fig. 6, which is proportional to the initial rate of product formation. This allows the rate constants at all three temperatures to be placed on the same scale using the Arrhenius law. In Fig. 7, the values derived from initial rates (×) are placed on an absolute scale by comparison to the rate constants obtained by exponential fitting at the lower two temperatures (o). The best fit Arrhenius parameter for the reaction rate constant ^ are = 0.76 ± 0.02 eV and $ = 10 . ± .b s ML . These results are consistent with independently obtained results using delay scanning.

Diffusion limited surface reaction rates
In the course of studying the behavior high rep-rate detection, we made what were, at first, surprising observations. We found that near the end of titrations when the rate of CO2 formation had nearly vanished, oxygen coverage remained on the surface in regions outside the crossing region of the two molecular beams.
By translating the crystal in a direction perpendicular to the surface normal, we could observe a sudden increase of the CO2 production rate. These observations indicated that a successful modelling of these experiments would require characterizing both the spatial as well as the temporal evolution of the reactant coverages.
In this section, we describe such modelling show that titration experiments often produce conditions where CO diffusion effects the rates of reactions.
We first imagine dividing the reacting surface into e spatial elements.
The total rate is the sum of reactive (Eq. 8 and 9) and diffusive (Eq. 10 and 11) contributions. In Eq. 10 and 11, l n and ∆[/ * ] f are the species specific and concentration independent diffusion coefficients and concentration gradients, respectively, used in application of Fick's second law of diffusion. I J,f is the time dependent incoming flux to the spatial element e, produced by the molecular beams.
The dosing function I J,f is described with a periodic function (in time) that reassembles the spatial, oKp f N, and temporal shape, 'q r, of our molecular beams. Specifically, we modelled it with: 'q r = cos n Kt uu q − rN, where uu is the repetition rate of the nozzle, the reference timing and / is an integer chosen to best represent the temporal shape of the beam. Using ion imaging, we experimentally determined the spatial intensity profile of each molecular beam, from which we deduced their radial profiles. Both molecular beams have a nominal projected diameter of 2 mm. For I J,f we use a flattop Gaussian that resembles the experimentally determined radial profile, oqp f r, of the beam, which is given by: where ~ and z are parameters representing the shape of the experimental beam profile. The combined and normalized dosing function is then given by: where • is the normalization to define the observed molecular beam flux.
We made sure that both molecular beams overlap on the surface and checked this by ensuring that oxygen coverage remained symmetrically distributed around the molecular beams crossing point at the end of the titration measurement. Hence, we conclude that our experiments approximately preserve radial symmetry, which allows us to solve the diffusion equations, Eq.'s 10 and 11, in polar coordinates. The diffusion formalism is derived in the Appendix to this paper. The rate equations including diffusion and reaction are solved numerically using LSODA from the Fortran ODEPACK library 19 . The concentrations of CO* and O* in each spatial element e are propagated in time.
To simulate measurements like those of Fig. 5, we initiate the model calculations with adsorbed oxygen produced by many pulses of the O2 beam. This requires an initial O * spatial profile (black line of Fig. 8) that is much broader than the nominal O2 beam profile (thick gray line of Fig. 8), as O* coverage quickly saturates near the center of the beam and after that only the wings of the O2 beam add additional O*. We simulated the spatial evolution of concentrations within a radial extent of 3 mm and with each spatial element, j, being 5 µm in size. The corresponding total CO2 formation rate is given by summing the rate of each spatial element e and weighting it by the respective area $ f , in the following manner: where the area of e G… spatial element is given by: The simulation accounts for the influence of reactions R -R as well as CO desorption and diffusion. The reaction rate constants were determined previously (see Sec. 3

.1 and 3.3). Oxygen desorption is unimportant
at these surface temperatures [20][21] . Oxygen diffusion is found to be unimportant under our conditions 22 . We estimated the diffusion coefficient for CO using an activation energy of 0.12eV from Ref. 23 and the fitted pre-factor for CO diffusion needed to obtain agreement with our measurements. The optimized pre-factor was 10 . †± . cm s . The CO diffusion rates we obtain in this way are consistent with previous measurements on Pd(111) -see the Appendix.
We derived the absolute incident beam fluxes from measurements of the steady state pressures of CO and O2 in the UHV chamber combined with a knowledge of the chamber pumping speed. The model results were insensitive to the O2 flux, but highly sensitive to the assumed CO flux. We found best agreement with experiment when using a CO flux ~30% smaller than that derived from our experimental estimate.
We used a 2 nd order Langmuir expression for the coverage dependent sticking coefficient of O2, with Z k Q , = 0.4 16 . Best agreement with the experiment is achieved when an oxygen coverage independent sticking coefficient of 0.6 ± 0.1 is used for CO. We assume that the sticking probability of CO decreases linearly with CO coverage. Fig. 9 (panel A) shows a comparison of this model to the titration experiment of Fig. 5. Note that the amplitude quickly decays over the first 300-400 CO beam pulses, thereafter decaying more slowly, behavior that is captured in the kinetic model. The transition between the fast and slow decay regions is accompanied by an increase of the baseline (shown in magnification in panel B of Fig. 9). This indicates a continuous production of CO2. The experimentally observed increase of the baseline is also present in the model.
Looking in more detail (panels C of Fig. 9), we find that the single pulse transient rate is decreasing with increasing titration time; furthermore, the transient rates are well reproduced by the kinetic model as is the continuous production of CO2 seen in the later stages of the titration.
The qualitative behavior can be understood by recalling that the amplitude of the titration curve reflects the initial rate of CO2 production, which is directly proportional to the oxygen coverage. With increasing titration time, the initial rate decreases indicating that the oxygen coverage is dropping. As a consequence of the reduced reaction rate at lower oxygen coverage, the lifetime of CO molecules on the surface increases, while the rate of CO adsorption remains constant. Since CO's desorption rate is slow at these temperatures, CO begins to build up from one molecular beam pulse to the next; this leads to quasicontinuous CO2 formation and to baseline increase at later times in the titration.
We also performed a sensitivity analysis of the fit to the titration kinetics. The degree of rate control 24 exhibited by the kinetic parameter, ‰ , is given by a sensitivity coefficient, Š ‰,J : where u J is the CO2 formation rate. A high absolute value of Š ‰,J indicates importance of the process to the reaction rate. A positive (negative) value of Š ‰,J means that an increase of the rate parameter produces an increase (decrease) of the CO2 formation rate. In Fig. 9 (panels D) we plot Š ‰,J for reaction (purple, dotted), CO desorption (blue, dashed) and CO diffusion (green, dash-dotted).
The reaction between CO* and O* dominates the rate of product formation up to a titration time of about 7 sec, thereafter CO desorption and diffusion become increasingly important. Between 12 and 24 seconds, where the three processes are of similar importance, their influence appears at different points in the kinetic trace. Consider the kinetic traces found at ~23 seconds. Here, the beginning of the kinetic trace is dominated by the influence of the reaction, whereas diffusion and desorption influence later times in the trace. Note that desorption decreases while diffusion increases the rate of CO2 production. This can be understood by realizing that at later stages of the titration, O* has been depleted near the center of the CO beam. Each new CO pulse produces higher CO* concentration in the "doughnut hole" of O* concentration (see Fig. 8).
These are the conditions where the quasi-continuous CO2 formation rate (i.e. the CO2 being produced prior to the next pulse) can appear as it is due to a diffusion-controlled reaction between CO* and O*. With our validated kinetic model we can also estimate the associated reaction front speed (see Appendix) which is a characteristic property that can be measured for spatio-temporal pattern formation. In Fig. 10 the reaction front speed at 503 K is shown as a function of titration time. Prior to 6 seconds after the beginning of the titration, no reaction front is formed. However, from 6 to 8 seconds titration time a reaction front forms and its speed accelerated to 175 µm/s. With increasing titration time, the front speed decreases reaching speeds of around 10 µm/s at titration times longer than 15 s. We emphasize that the derived values of front speed are similar to those obtained in previous work for CO oxidation on Pt(110), which ranged from 1 to 100 µm/sec [25][26] . The fact that we derive front speeds nearly a factor of two higher than that work, probably results from faster thermal diffusion for CO on Pd compared to Pt 23,27 .
It is important to highlight that we have modeled the real time titration experiment without coverage dependent rate constants. Since we achieve good correspondence with the experiment, we claim that the rate constants have weak dependence on oxygen coverage in CO oxidation on Pd(332). However, this is in contradiction to the findings that were previously made on Pd(111) by Engel and Ertl 18 . We also find that our reaction rate constant is about 4-8 times higher than those reported from Pd(111). We think that steps lead to a higher reaction rate, consistent with previous observations on Pd 28 and Pt 6 . The reason why we have not taken reaction at steps and terraces explicitly into account is that we have not needed it for good match with the experiment. This is probably due to a rather fast exchange of CO and O atoms between terraces and steps which leads to an effective reaction rate composed of both reactions at steps and terraces.
We plan to investigate the details of the kinetic mechanism of CO oxidation at steps and terraces of Pd further in future.

Conclusions
This work shows how high-repetition rate lasers and ion imaging detection can be used to obtain the kinetic traces of catalytic processes from a single molecular beam pulse, overcoming the need for delay scans that are typical for pump-probe methods. The new approach provides an increased duty cycle resulting in rates of acquisition for kinetic data that are 10-1000 times faster than conventional delay scanning methods. The new method can measure rates over 5-7 orders of magnitude, dramatically better than when using delay scanning. The method is particularly attractive for measuring slow processes where temporal dilution would make delay scanning impossible. Our results are also consistent with an oxygen coverage independent sticking coefficient of CO of 0.6 ± 0.1. The desorption and diffusion rate constants of CO agree well with the parameters determined earlier from Pd(111), indicating that CO has no energetic preference for steps and that they are not influencing its mobility on the surface. The reaction rate constant is found to be approximately a factor 4-8 higher than previous reports for Pd (111), indicating that steps are more reactive for CO oxidation on Pd than terraces.
While in this work we were limited to a detection rep-rate of 1−kHz due to the fact that we used a Ti:Sapphire laser, we plan to extend our capabilities to a detection rate of 100 kHz and study reaction rates at changing catalysts conditions in more detail using a Yb-Fiber laser. We think that this method offers the possibility to accurately study catalytic reaction rates and kinetic mechanisms at the intersection between the well-defined conditions that are desirable for surface science and the more dynamic conditions relevant to industrial catalysis. Fig. 1: Comparison of delay scanning versus high rep-rate detection employed in velocity resolved kinetics measurements. (a) Delay scanning involves the acquisition of many (e.g. 50) images at each time delay between the initiating molecular beam pulse and the laser ionization pulse. Points in the kinetic trace recorded by scanning the delay between a molecular beam pulse that initiates the reaction and a laser ionization pulse that detects the products. The catalytic system must be stable throughout the course of the delay scanning procedure. (b) High rep-rate detection with high-speed imaging records many points in the kinetic trace for each molecular beam pulse. Here, the molecular beam initates the reaction every 0.1 s and points in the kinetic trace are recorded by each pulse of a 1−kHz detection laser. The duty cycle of this method can be much higher than delay scanning. Furthermore, the kinetics can be recorded while the catalyst is changing.    High rep-rate detection of velocity resolved kinetics for a non-stationery catalyst. The kinetics of CO oxidation on Pd(332) are recorded starting with saturated oxygen coverage. Adsorbed oxygen is removed during the experiment and the kinetics change accordingly. The surface temperature was 503 K and the CO beam operated at 50 Hz. The CO beam clean up a pre-oxidized surface that had been exposed to 300 ± 80 ML of O2. The raw data are shown as blue lines, the Savitzky-Golay filtered data is shown as black lines. Kinetic fits (first order decay convoluted with incident beam shape) are shown as red dashed lines in the insets. The insets are indicated by colored bars and borders. The gray dash-dotted line in the insets indicates the reaction time at which the reaction is initiated by the pulsed CO beam.  Temperature dependence of CO oxidation rate constants at saturated oxygen coverage. The first order rate constants for CO2 formation determined from the data of Fig. 6. The circles are the rate constants determined from the shape of the single pulse kinetic trace and crosses are initial rates determined from their amplitude. The initial rates are scaled to match the first order rate constants at low temperature. The dotted line is the limit above which the rate constants cannot be derived from the shape of the kinetic trace. The dash-dotted line is the kinetic resolution (in this experiment around 110 µs) above which no kinetic information can be derived from transient kinetics. The full and dashed curves are Arrhenius fits to the initial rate and the first order rate constant, respectively. Fig. 8: Spatial distributions of adsorbed oxygen atoms during a CO oxidation titration. The distributions are assumed cylindrically symmetrical about the CO beam axis. The radial distance from the CO beam center-line is shown on the '-axis. The solid black line indicates the initial oxygen coverage distribution produced by long exposure with a molecular beam of O2. The radial distribution of the CO beam (gray thick line) peaks near 0 and preferentially removes O-atoms there. As time progresses, a "doughnut hole" reaction develops, where the CO is concentrated along the CO beams center-line and adsorbed oxygen atoms form a ring around the CO beam. In later stages of titration, the reaction forms a front where the CO and O concentrations overlap. Diffusion of CO from the center of the doughnut hole to the oxygen ring also influences the reaction rate.