Seven-dimensional microcanonical treatment of hydrogen dissociation dynamics on Cu(111): clarifying the essential role of surface phonons

A simple picture of the hydrogen dissociation/associative desorption dynamics on Cu(111) emerges from a two-parameter, full dimensionality microcanonical unimolecular rate theory (MURT) model of the gas-surface reactivity. Vibrational frequencies for the reactive transition state were taken from density functional theory calculations of a six-dimensional potential energy surface [Hammer et al., Phys. Rev. Lett. 73, 1400 (1994)]. The two remaining parameters required by the MURT were fixed by simulation of experiments. These parameters are the dissociation threshold energy, E(0)=79 kJmol, and the number of surface oscillators involved in the localized H(2)Cu(111) collision complex, s=1. The two-parameter MURT quantitatively predicts much of the varied behavior observed for the H(2) and D(2)Cu(111) reactive systems, including the temperature-dependent associative desorption angular distributions, mean translational energies of the associatively desorbing hydrogen as a function of rovibrational eigenstate, etc. The divergence of the statistical theory's predictions from experimental results at low rotational quantum numbers, J < or approximately 5, suggests that either (i) rotational steering is important to the dissociation dynamics at low J, an effect that washes out at high J, or (ii) molecular rotation is approximately a spectator degree of freedom to the dissociation dynamics for these low J states, the states that dominate the thermal reactivity. Surface vibrations are predicted to provide approximately 30% of the energy required to surmount the activation barrier to H(2) dissociation under thermal equilibrium conditions. The MURT with s=1 is used to analytically confirm the experimental finding that partial differential "E(a)(T(s))" partial differential E(t)= -1 for eigenstate-resolved dissociative sticking at translational energies E(t)相似文献   

Dissociative chemisorption and energy transfer for methane on Ir(111)

A 3-parameter local hot spot model of gas-surface reactivity is employed to analyze and predict dissociative sticking coefficients for CH(4) incident on Ir(111) under varied nonequilibrium and equilibrium conditions. One Ir surface oscillator and the molecular vibrations, rotations, and translational energy directed along the surface normal are treated as active degrees of freedom in the 14 dimensional microcanonical kinetics. The threshold energy for CH(4) dissociative chemisorption on Ir(111) derived from modeling molecular beam experiments is E(0) = 39 kJ/mol. Over more than 4 orders of magnitude of variation in sticking, the average relative discrepancy between the beam and theoretically derived sticking coefficients is 88%. The experimentally observed enhancement in dissociative sticking as beam translational energies decrease below approximately 10 kJ/mol is consistent with a parallel dynamical trapping/energy transfer channel that likely fails to completely thermalize the molecules to the surface temperature. This trapping-mediated sticking, indicative of specific energy transfer pathways from the surface under nonequilibrium conditions, should be a minor contributor to the overall dissociative sticking at thermal equilibrium. Surprisingly, the CH(4) dissociative sticking coefficient predicted for Ir(111) surfaces at thermal equilibrium, based on the molecular beam experiments, is roughly 4 orders of magnitude higher than recent measurements on supported nanoscale Ir catalysts at 1 bar pressure, which suggests that substantial improvements in catalyst turnover rates may be possible.     

Effusive molecular beam study of C2H6 dissociation on Pt(111)

The dissociative sticking coefficient for C2H6 on Pt(111) has been measured as a function of both gas temperature (Tg) and surface temperature (Ts) using effusive molecular beam and angle-integrated ambient gas dosing methods. A microcanonical unimolecular rate theory (MURT) model of the reactive system is used to extract transition state properties from the data as well as to compare our data directly with supersonic molecular beam and thermal equilibrium sticking measurements. We report for the first time the threshold energy for dissociation, E0 = 26.5 +/- 3 kJ mol(-1). This value is only weakly dependent on the other two parameters of the model. A strong surface temperature dependence in the initial sticking coefficient is observed; however, the relatively weak dependence on gas temperature indicates some combination of the following (i) not all molecular excitations are contributing equally to the enhancement of sticking, (ii) that strong entropic effects in the dissociative transition state are leading to unusually high vibrational frequencies in the transition state, and (iii) energy transfer from gas-phase rovibrational modes to the surface is surprisingly efficient. In other words, it appears that vibrational mode-specific behavior and/or molecular rotations may play stronger roles in the dissociative adsorption of C2H6 than they do for CH4. The MURT with an optimized parameter set provides for a predictive understanding of the kinetics of this C-H bond activation reaction, that is, it allows us to predict the dissociative sticking coefficient of C2H6 on Pt(111) for any combination of Ts and Tg even if the two are not equal to one another.     

Using effusive molecular beams and microcanonical unimolecular rate theory to characterize CH4 dissociation on Pt(111)

The dissociative sticking coefficient for CH4 on Pt(111) has been measured as a function of both gas temperature (Tg) and surface temperature (Ts) using effusive molecular beam and angle-integrated ambient gas dosing methods. The experimental results are used to optimize the three parameters of a microcanonical unimolecular rate theory (MURT) model of the reactive system. The MURT calculations allow us to extract transition state properties from the data as well as to compare our data directly to other molecular beam and thermal equilibrium sticking measurements. We find a threshold energy for dissociation of E0 = 52.5 +/- 3.5 kJ mol(-1). Furthermore, the MURT with an optimized parameter set provides for a predictive understanding of the kinetics of this C-H bond activation reaction, that is, it allows us to predict the dissociative sticking coefficient of CH4 on Pt(111) for any combination of Ts and Tg even if the two are not equal to one another, indeed, the distribution of molecular energy need not even be thermal. Comparison of our results to those from recent thermal equilibrium catalysis studies on CH4 reforming over Pt nanoclusters ( approximately 2 nm diam) dispersed on oxide substrates indicates that the reactivity of Pt(111) exceeds that of the Pt nanocatalysts by several orders of magnitude.     

Dynamically biased RRKM model of activated gas-surface reactivity: vibrational efficacy and rotation as a spectator in the dissociative chemisorption of CH4 on Pt(111)

The reactivity of CH(4) impinging on a Pt(111) surface was examined using a precursor-mediated microcanonical trapping model of dissociative chemisorption wherein the effects of rotational and vibrational energy could be explored. Dissociative sticking coefficients for a diverse range of non-equilibrium effusive beam, supersonic beam, and eigenstate-resolved experiments were simulated and an average relative discrepancy between theory and experiment of better than 50% was achieved by treating molecular rotations and translation parallel to the surface as spectator degrees of freedom, and introducing a dynamically-biased vibrational efficacy. The model parameters are {E(0) = 57.9 kJ mol(-1), s = 2, η(v) = 0.40} where E(0) is the apparent threshold energy for reaction, s is the number of surface oscillators participating in energy exchange within each gas-surface collision complex formed, and η(v) is the mean vibrational efficacy for reaction relative to normal translational energy which figures in the assembly of the active exchangeable energy which is available to surmount the activation barrier to dissociative chemisorption. GGA-DFT electronic structure calculations provided vibrational frequencies for the transition state for dissociative chemisorption. The asymmetry of the rotational state populations in supersonic and effusive molecular beam experiments allowed kinetic analysis to establish that taking rotation as a spectator degree of freedom is a good approximation. Surface phonons, rather than the incident molecules, are calculated to play the dominant role in supplying the energy required to overcome the activation barrier for dissociative chemisorption under the thermal equilibrium conditions relevant to high pressure catalysis. Over the temperature range 300 K ≤T≤ 1000 K, the thermal dissociative sticking coefficient is predicted to be well described by S(T) = S(0) exp(-E(a)/RT) where S(0) = 0.62 and E(a) = 62.6 kJ mol(-1).     

Microcanonical unimolecular rate theory at surfaces. II. Vibrational state resolved dissociative chemisorption of methane on Ni(100)

A three-parameter microcanonical theory of gas-surface reactivity is used to investigate the dissociative chemisorption of methane impinging on a Ni(100) surface. Assuming an apparent threshold energy for dissociative chemisorption of E(0)=65 kJ/mol, contributions to the dissociative sticking coefficient from individual methane vibrational states are calculated: (i) as a function of molecular translational energy to model nonequilibrium molecular beam experiments and (ii) as a function of temperature to model thermal equilibrium mbar pressure bulb experiments. Under fairly typical molecular beam conditions (e.g., E(t)>/=25 kJ mol(-1), T(s)>/=475 K, T(n)/=100 K the dissociative sticking is dominated by methane in vibrationally excited states, particularly those involving excitation of the nu(4) bending mode. Fractional energy uptakes f(j) defined as the fraction of the mean energy of the reacting gas-surface collision complexes that derives from specific degrees of freedom of the reactants (i.e., molecular translation, rotation, vibration, and surface) are calculated for thermal dissociative chemisorption. At 500 K, the fractional energy uptakes are calculated to be f(t)=14%, f(r)=21%, f(v)=40%, and f(s)=25%. Over the temperature range from 500 K to 1500 K relevant to thermal catalysis, the incident gas-phase molecules supply the preponderance of energy used to surmount the barrier to dissociative chemisorption, f(g)=f(t)+f(r)+f(v) approximately 75%, with the highest energy uptake always coming from the molecular vibrational degrees of freedom. The predictions of the statistical, mode-nonspecific microcanonical theory are compared to those of other dynamical theories and to recent experimental data.     

Nonequilibrium activated dissociative chemisorption: SiH4 on Si(100)

A three-parameter local hot spot model of gas-surface reactivity is employed to analyze and predict dissociative sticking coefficients for SiH4 incident on Si(100) under varied nonequilibrium conditions. Two Si surface oscillators and the molecular vibrations, rotations, and translational energy directed along the local surface normal are active degrees of freedom in the 15 dimensional microcanonical kinetics. The threshold energy for SiH4 dissociative chemisorption is found to be 19 kJ/mol, in quantitative agreement with recent GGA-DFT calculations. A simple scheme for increasing the rate of chemical vapor deposition of silicon from SiH4 at low surface temperatures is modeled.     

Microcanonical unimolecular rate theory at surfaces. III. Thermal dissociative chemisorption of methane on Pt(111) and detailed balance

A local hot spot model of gas-surface reactivity is used to investigate the state-resolved dynamics of methane dissociative chemisorption on Pt(111) under thermal equilibrium conditions. Three Pt surface oscillators, and the molecular vibrations, rotations, and the translational energy directed along the surface normal are treated as active degrees of freedom in the 16-dimensional microcanonical kinetics. Several energy transfer models for coupling a local hot spot to the surrounding substrate are developed and evaluated within the context of a master equation kinetics approach. Bounds on the thermal dissociative sticking coefficient based on limiting energy transfer models are derived. The three-parameter physisorbed complex microcanonical unimolecular rate theory (PC-MURT) is shown to closely approximate the thermal sticking under any realistic energy transfer model. Assuming an apparent threshold energy for CH(4) dissociative chemisorption of E(0)=0.61 eV on clean Pt(111), the PC-MURT is used to predict angle-resolved yield, translational, vibrational, and rotational distributions for the reactive methane flux at thermal equilibrium at 500 K. By detailed balance, these same distributions should be observed for the methane product from methyl radical hydrogenation at 500 K in the zero coverage limit if the methyl radicals are not subject to side reactions. Given that methyl radical hydrogenation can only be experimentally observed when the CH(3) radicals are kinetically stabilized against decomposition by coadsorbed H, the PC-MURT was used to evaluate E(0) in the high coverage limit. A high coverage value of E(0)=2.3 eV adequately reproduced the experimentally observed methane angular and translational energy distributions from thermal hydrogenation of methyl radicals. Although rigorous application of detailed balance arguments to this reactive system cannot be made because thermal decomposition of the methyl radicals competes with hydrogenation, approximate applicability of detailed balance would argue for a strong coverage dependence of E(0) with H coverage--a dependence not seen for methyl radical hydrogenation on Ru(0001), but not yet experimentally explored on Pt(111).     

Adsorption dynamics of water on Pt{110}-(1 x 2)

The dynamics of H(2)O adsorption on Pt{110}-(1 x 2) is studied using supersonic molecular beam and temperature programed desorption techniques. The sticking probabilities are measured using the King and Wells method at a surface temperature of 165 K. The absolute initial sticking probability s(0) of H(2)O is 0.54+/-0.03 for an incident kinetic energy of 27 kJmol. However, an unusual molecular beam flux dependence on s(0) is also found. At low water coverage (theta<1), the sticking probability is independent of coverage due either to diffusion in an extrinsic precursor state formed above bilayer islands or to incorporation into the islands. We define theta=1 as the water coverage when the dissociative sticking probability of D(2) on a surface predosed with water has dropped to zero. The slow falling H(2)O sticking probability at theta>1 results from compression of the bilayer and the formation of multilayers. Temperature programed desorption of water shows fractional order kinetics consistent with hydrogen-bonded islands on the surface. A remarkable dependence of the initial sticking probability on the translational (1-27 kJ/mol) and internal energies of water is observed: s(0) is found to be essentially a step function of translational energy, increasing fivefold at a threshold energy of 5 kJ/mol. The threshold migrates to higher energies with increasing nozzle temperature (300-700 K). We conclude that both rotational state and rotational alignment of the water molecules in the seeded supersonic expansion are implicated in dictating the adsorption process.     

Communication: angle-resolved thermal dissociative sticking of CH4 on Pt(111): further indication that rotation is a spectator to the gas-surface reaction dynamics

Effusive molecular beam measurements of angle-resolved thermal dissociative sticking coefficients for CH(4) impinging on a Pt(111) surface, at a temperature of 700 K, are reported and compared to theoretical predictions. The reactivity falls off steeply as the molecular angle of incidence increases away from the surface normal. Successful modeling of the thermal dissociative sticking behavior, consistent with existent CH(4) supersonic molecular beam experiments involving rotationally cold molecules, required that rotation be treated as a spectator degree of freedom.     

On the dynamics of H2 adsorption on the Pt(111) surface

The dynamics and kinetics of the dissociation of hydrogen over the hexagonal close packed platinum (Pt(111)) surface are investigated using Car–Parrinello molecular dynamics and static density functional theory calculations of the potential energy surfaces. The calculations model the reference energy‐resolved molecular beam experiments, considering the degrees of freedom of the catalytic surface. Two‐dimensional potential energy surfaces above the main sites on Pt(111) are determined. Combined with Car–Parrinello trajectories, they confirm the dissociative adsorption of H₂ as the only adsorption pathway on this surface at H₂ incindence energies above 5 kJ/mol. A direct determination of energy‐resolved sticking coefficients from molecular dynamics is also performed, showing an excellent agreement with the experimental data at incidence energies in the 5–30 kJ/mol range. Application of dispersion corrections does not lead to an improvement in the prediction of the H₂ sticking coefficient. The adsorption reaction rate obtained from the calculated sticking coefficients is consistent with experimentally derived literature values.     

Dynamics of analyte binding onto a metallophthalocyanine: NO/FePc

The gas-surface reaction dynamics of NO impinging on an iron(II) phthalocyanine (FePc) monolayer were investigated using King and Wells sticking measurements. The initial sticking probability was measured as a function of both incident molecular beam energy (0.09-0.4 eV) and surface temperature (100-300 K). NO adsorption onto FePc saturates at 3% of a monolayer for all incident beam energies and surface temperatures, suggesting that the final chemisorption site is confined to the Fe metal centers. At low surface temperature and low incident beam energy, the initial sticking probability is 40% and decreases linearly with increasing beam energy and surface temperature. The results are consistent with the NO molecule sticking onto the FePc molecules via physisorption to the aromatics followed by diffusion to the Fe metal center, or precursor-mediated chemisorption. The adsorption mechanism of NO onto FePc was confirmed by control studies of NO sticking onto metal-free H2Pc, inert Au111, and reactive Al111.     

Methylthiolate on Au(111): adsorption and desorption kinetics

Low energy electron diffraction, Auger electron spectroscopy, X-ray photoelectron spectroscopy and line of sight mass spectrometry have been used to study the adsorption and desorption of dimethyldisulfide (DMDS) on Au(111). At 300 K adsorption is dissociative, forming a chemisorbed adlayer of methylthiolate with a 1/3 ML, (sq rt 3 x sq rt 3)R30 degrees, structure. At 100 K adsorption is molecular, with dissociation to form the 1/3 ML (sq rt 3 x sq rt 3)R30 degrees methylthiolate structure occurring at 138-160 K. A physisorbed DMDS layer, with a coverage of 1/6 ML of DMDS, forms on top of the (sq rt 3 x sq rt 3)R30 degrees chemisorbed MT surface for T < or = 180 K, with multilayers forming for T < or = 150 K. In temperature programmed desorption, multilayers of DMDS desorbed with zero order kinetics and an activation energy of 41 kJ mol(-1); the physisorbed layer desorbed with first order kinetics, exhibiting repulsive lateral interactions with an activation energy which varied from 63 kJ mol(-1) (theta = 0) to 51 kJ mol(-1) (theta = 1); the chemisorbed methylthiolate layer desorbed associatively as DMDS via the physisorbed layer, the activation energy for the reaction, 2 methylthiolate --> physisorbed DMDS, exhibiting repulsive lateral interactions with an activation energy which varied from 65 kJ mol(-1) (theta = 0) to 61 kJ mol(-1) (theta = 1). The physisorbed disulfide layer explains the pre-cursor state adsorption kinetics observed in sticking probability measurement, while its relatively facile formation provides a mechanism by which thiolate self-assembled monolayers can become mobile at room temperature.     

TPEPICO spectroscopy of vinyl chloride and vinyl iodide: neutral and ionic heats of formation and bond energies

The 0 K dissociative ionization onsets of C2H3X --> C2H3(+) + X (X = Cl, I) are measured by threshold photoelectron-photoion coincidence spectroscopy. The heats of formation of C2H3Cl (Delta H(f,0K)(0) = 30.2 +/- 3.2 kJ mol(-1) and Delta(H f,298K)(0) = 22.6 +/- 3.2 kJ mol(-1)) and C2H3I (Delta(H f,0K)(0) = 140.2 +/- 3.2 kJ mol(-1) and Delta(H f,298K)(0) = 131.2 +/- 3.2 kJ mol(-1)) and C- X bond dissociation enthalpies as well as those of their ions are determined. The data help resolve a longstanding discrepancy among experimental values of the vinyl chloride heat of formation, which now agrees with the latest theoretical determination. The reported vinyl iodide heat of formation is the first reliable experimental determination. Additionally, the adiabatic ionization energy of C2H3I (9.32 +/- 0.01 eV) is measured by threshold photoelectron spectroscopy.     

平动能对激光诱导气-固表面蚀刻反应的增强效应

本文采用超声分子束和时间分辨质谱技术研究了入射分子平动能对激光诱导气-固表面反应的增强效应. 对于由可见激光(56 nm)诱导的Cl_2与Ge(111)、Si(111)和GaAS(100)表面蚀刻反应, 研究发现提高Cl_2分子的入射平功能将明显地增加反应产率, 而且都存在一个入射平动能的反应阈值, 其数值为5～7 kJ·mol~(-1). 此外, 从反应产物的飞行时间谱测得入射分子平动能对产物平动温度的影响. 这些结果可以通过平动能促进Cl_2分子在表面上解离化学吸附过程来解释。     

The role of molecular rotation in activated dissociative adsorption on metal surfaces

The role of molecular rotation in dissociative adsorption of H2 on the activated NiAl(110) metal surface is systematically investigated by means of classical dynamics calculations performed on ab initio six-dimensional potential energy surfaces. The calculations show that molecules rotate abruptly when they are close to the surface and that this rotation allows the molecules to adopt the orientation that is more convenient for dissociation (i.e., nearly parallel to the surface). Also, in reactive sectors of the NiAl(110) unit cell, there is an "angular threshold" below which molecules cannot dissociate. This angular threshold goes down as the incidence energy increases, which explains the rise of the dissociation probability and the fact that it reaches a value close to 1 at incidence energies of the order of 2 eV. The fact that switching on molecular rotation favors dissociation establishes a competition between dissociation and rotational excitation of reflected molecules above the dissociation threshold. Measurements on rotational excitation might thus bring indirect evidence on the dissociation dynamics. Sample calculations for nonactivated Pd(111) and activated Cu(110) metal surfaces suggest that some of these conclusions may be of general validity.     

Rotationally inelastic scattering of HD from Cu(100) and Pd(111)

Rotational excitation of HD scattered from Cu(100), Pd(111), and Pd(111):H(D) was measured using molecular beam and quantum-state-specific laser spectroscopy techniques. Greater than 91% of the incident HD population was in the v = 0, J = 0 state. The final rotational distributions from Cu(100), Pd(111), and Pd(111):H(D) were compared for a HD beam at an incident energy of 74 meV. For all the three surfaces studied, rotationally inelastic scattering probabilities were large. We find that the final HD rotational distributions are remarkably similar for the three surfaces even though Pd(111) is very reactive to dissociative adsorption of HD whereas Cu(100) and Pd(111):H(D) are chemically inert.     

Heat of adsorption of naphthalene on Pt(111) measured by adsorption calorimetry

The heat of adsorption of naphthalene on Pt(111) at 300 K was measured with single-crystal adsorption calorimetry. The heat of adsorption on the ideal, defect-free surface is estimated to be (300 - 34 - 199(2)) kJ/mol. From this, a C-Pt bond energy for aromatic hydrocarbons on Pt(111) of approximately 30 kJ/mol is estimated, consistent with earlier results for benzene on Pt(111). There is higher heat of adsorption at very low coverage, attributed to step sites where the adsorption heat is >/=330 kJ/mol. Saturation coverage, = 1 ML, corresponds to 1.55 x 10(14) molecules/cm(2). Sticking probability measurements of naphthalene on Pt(111) give a high initial value of 1.0 and a Kisliuk-type coverage dependence that implies precursor-mediated sticking. The ratio of the hopping rate to the desorption rate of this precursor is approximately 51. Naphthalene adsorbs transiently on top of chemisorbed naphthalene molecules with a heat of adsorption of 83-87 kJ/mol.     

Density functional theory study of H and H2 interacting with NiAl(110)

We present results of extensive density functional theory (DFT) calculations for H and H2 interacting with NiAl(110). Continuous representations of the full dimensional potential energy surface (PES) for the H/NiAl(110) and H2/NiAl(110) systems are obtained by interpolation of the DFT results using the corrugation reducing procedure. We find a minimum activation energy barrier of approximately 300 meV for dissociative adsorption of H2, which is consistent with the energy threshold obtained in molecular beam experiments for H2 (nu=0). We explain vibrational enhancement observed in experiments as the consequence of vibrational softening in the entrance channel over the most reactive surface site. The H2/NiAl(110) PES shows a high surface site selectivity: for energies up to 0.1 eV above threshold, H2 adsorption can only take place around top-Ni sites (within a circle of radius approximately 0.3 A). A strong energetic corrugation is observed: energy barriers for dissociation vary by more than 1 eV between the most and the least reactive sites. In contrast, geometric corrugation is much less pronounced and comparable to that of low index single metal surfaces like Cu or Pt.