Bin
Jiang
ab,
Daiqian
Xie
*a and
Hua
Guo
*b
aInstitute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China. E-mail: dqxie@nju.edu.cn
bDepartment of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA. E-mail: hguo@unm.edu
First published on 23rd October 2012
The bond selectivity in dissociative chemisorption of HOD on Cu(111) is investigated using a six-dimensional quantum model. It includes all vibrational modes of the impinging molecule on a density functional theory based interaction potential between the molecule and metal surface. It is shown that excitations in the HOD local stretching modes selectively enhance cleavage of the excited bond. This pronounced bond selectivity is attributed to a “late” or “product-like” barrier on the potential energy surface for the dissociative chemisorption and the slow intramolecular vibrational energy redistribution in the water molecule. The existence of mode and bond selectivities also underscores the inadequacy of statistical based transition-state theory in describing this industrially important surface reaction.
One such example is methane (CH4) dissociative chemisorption on metal surfaces.4 This reaction, producing adsorbed CH3 and H species, is the rate-limiting step in the steam reforming of methane, which is a key heterogeneous catalytic process in the production of H2.5 It has recently been shown in several landmark experiments that vibrational excitation in methane is very efficient in promoting the reaction, sometimes more so than the same amount of translational energy.6–12 In addition, clear bond selectivity was found in dissociative chemisorption of partially deuterated methane, in which vibrational excitation controls the products by selectively cleaving a specific bond.13 Similar mode and bond selectivities have also been demonstrated in gas phase bimolecular reactions involving methane.14–17 The selectivities can be attributed to two major factors. First, these reactions have “late” or “product-like” transition states, for which the appropriate vibrational coordinate(s) are often more efficient in promoting barrier crossing, according to Polanyi's rules.18 Second, methane has relatively inefficient IVR due to its local mode nature of vibration,19 so that the energy deposited into a single vibrational mode does not easily leak out to other degrees of freedom. However, there have been several recent reports on departures from Polanyi's rules,20,21 underscoring the complexity of multidimensional dynamics.22–24
The mode and bond specificities found in these systems cannot be readily explained by the widely used transition-state theory,25 in which all types of energy contribute equally. Instead, a dynamical model is needed. However, a complete theoretical understanding of methane dissociative chemisorption dynamics has been hampered by the necessarily large number of degrees of freedom (fifteen if surface is assumed rigid) involved. The commonly used reduced-dimensional models are often difficult to correctly describe the vibrational modes in CH4, thus inadequate in addressing the mode and bond selectivity in this process.26–32 Recently, full-dimensional models for this system have been reported,33–35 some with approximate dynamical treatments.34,35 In addition, a quantum dynamic treatment, which is necessary because of the large zero-point energy and tunneling involved in the reaction, is still formidable if all dimensions are included,36 and the difficulties are compounded by the complexity in developing an accurate high-dimensional potential energy surface (PES).33–35,37
The mode and bond selectivities are not restricted to methane. With similar properties, H2O has, for example, been shown to exhibit similar behaviors in gas phase bimolecular reactions.38–42 Its dissociative chemisorption on most transition metal surfaces, which produces adsorbed H and OH species, should also behave in a similar fashion, because it has slow IVR due to its local mode nature,19 and the reaction also has a “late” barrier.43 This reaction is ubiquitous and often obligatory in many industrially important heterogeneous catalytic processes5 and represents the rate-limiting step in low-temperature water–gas shift (WGS) reaction on copper catalysts.44 Hence, the ability to control such a reaction may have important implications in many heterogeneous processes involving water. While there has not yet been any report of experimental exploration, we have recently demonstrated, using a reliable quantum mechanical model on a realistic PES based on ∼25000 density functional theory (DFT) points, that the dissociative chemisorption of H2O on Cu(111) exhibits strong mode selectivity.45 Specifically, vibrational energy in all three modes was found to be more effective in promoting the reaction than translational energy. In this publication, we address a different but related question: can vibrational excitation in the singly deuterated water (HOD) control the bond cleavage in its dissociative chemisorption on the same copper surface?
Since a detailed description of the DFT calculations of the PES and the six-dimensional (6D) dynamic model used for H2O dissociative chemisorption on the rigid Cu(111) surface has been presented in our recent publication,45 only a brief discussion concerning modifications for HOD will be given here. More details can be found in the ESI.† The Jacobi coordinates used in the 6D surface-triatom model is shown in Fig. 1 and the Hamiltonian can be expressed as follows (ℏ = 1 hereafter):
![]() | (1) |
![]() | ||
Fig. 1 Coordinates used in the six-dimensional model for HOD dissociative chemisorption on Cu(111). |
The Hamiltonian was represented with a mixed discrete variable representation (DVR) and finite basis representation (FBR).47 The z and r2 degrees of freedom were discretized with the sine DVR,48 while r1 was treated as non-reactive with a potential optimized DVR (PODVR).49,50 A non-direct product FBR basis (Kj(cos θ1)
KJ(cos θ2)ΦK(ϕ)) was used to evaluate the rotational kinetic energy terms, where j, J, and K are the angular momentum quantum numbers associated with OH(OD), HOD, and the projection of Ĵ on the molecule-fixed z-axis which is along r2.
Kj(cos θ1) and
KJ(cos θ2) are normalized associated Legendre polynomials, and ΦK(ϕ) is the exponential Fourier function (
). We further took advantage of the inversion symmetry of the potential with respect to φ in our model, so as to expand the wave function of even/odd parities in terms of cos(Kφ)/sin(Kφ) basis functions. In practice, three one-dimensional pseudo-spectral transformations were performed in sequence to transform the wavefunction from FBR to DVR, in order to evaluate the action of the potential energy operator on a grid.51
The dissociation probabilities were calculated with a modified version of the transition-state wave packet approach52 with the Chebyshev propagator,53 as discussed in detail in our recent work,45,54 and in ESI.† Briefly, two dividing surfaces are required. An initial wave packet is constructed on the first dividing surface placed in the asymptote as a direct product of the positive one-dimensional eigenstate of the flux operator and the specific rovibrational eigenstate for the reactant. This wavepacket is then propagated by the Chebyshev propagator toward the second dividing surface placed behind the transition state. The energy-dependent total reaction flux is calculated at the second dividing surface via a cosine Fourier transform of the Chebyshev propagation states, yielding the dissociation probability on a pre-specified energy grid.
While the O–H and O–D bond cleavages can in principle be treated within a single wave packet propagation, the grid needed to accommodate both processes is rather large. We chose to carry out two separate calculations for each initial state, one for the HOD → OH + D reaction and the other for the HOD → OD + H reaction. In both cases, the non-reactive bond can be described with a relatively small grid. This is possible because the two transition states for the dissociative chemisorption of HOD are well separated and have negligible interference, as in the gas phase HOD + H reaction.42
Numerous convergence tests have been performed and the final parameters used in the calculations are listed in Table 1. It is interesting to note that the number of grid points in θ1 for the OH bond cleavage is significantly larger than that for the OD bond cleavage, presumably due to the larger moment of inertia of the OD moiety. Similar observations have been made for the HOD + H reaction.42
Parameters | HOD/Cu(111) | |
---|---|---|
OD + H | OH + D | |
Translational coordinate (z) | 260 sine-DVR in [2.0,14.5] | 280 sine-DVR in [2.0,14.5] |
Internal coordinate (r2) | 25 sine-DVR in [1.2,5.5] | 25 sine-DVR in [1.2,5.5] |
Nonreactive coordinate (r1) | 5 PODVR | 5 PODVR |
Quadrature points in (θ1, θ2, φ) = (jmax + 1, Jmax + 1, Kmax + 1) | (28, 30, 28) | (18, 32, 18) |
Dividing surface S1 | z flux = 12.4 | z flux = 12.4 |
Dividing surface S2 | r 2flux = 3.2 | r 2flux = 3.2 |
Absorbing potential in z |
![]() |
![]() |
Absorbing potential in r2 |
![]() |
![]() |
Potential cut-off (eV) | 5.0 | 5.0 |
Propagation step | 4000 | 4000 |
Since the O–H and O–D stretches of HOD are well-defined local modes,19 the three vibrational modes of HOD are most conveniently labeled as the OD stretching, bending, and OH stretching modes, namely HOD(vOD, vb, vOH). Their calculated frequencies are 2561, 1353, and 3486 cm−1 on our PES, in reasonable agreement with the experimental values 2727, 1402, and 3707 cm−1.55 Experimentally, the excitation of individual quantum states in HOD can be readily achieved with IR lasers. The initial states in our wave packet calculations were chosen to be associated with the ground and first excited state of rotationless HOD in each of the three vibrational modes.
In Fig. 2, we first compare the dissociation probabilities for the ground (000) vibrational state of H2O and HOD as a function of the translational energy. The inset of the figure plots the dissociation probabilities in logarithmic scale, which shows an exponential increase of the reaction probabilities at low energies due apparently to tunneling. While the tunneling in the OD + H channel below the classical barrier of 0.96 eV is substantial, the OH + D channel shows much less reactivity due apparently to the heavier mass of the deuterium. It is interesting to note that the reaction threshold for the OD + H channel in the HOD(000) dissociation is somewhat lower than that for the other channel. This is in accord with the lower vibrational frequency of OD than that of OH. However, the dissociation probability in this channel becomes smaller than the OH + D channel at higher energies. Overall, the total dissociation probability for H2O(000) is much larger than the combined dissociation probabilities in both the OH and OD channels for HOD(000), which was also found in the gaseous H + HOD reaction.42 A clear explanation for this behavior is still lacking.
![]() | ||
Fig. 2 Comparison of dissociative chemisorption probabilities for HOD and H2O on Cu(111) as a function of the translational energy. The probabilities are plotted in logarithmic scale in the inset. |
The influence of initial vibration excitations on the dissociation probabilities for the dissociative chemisorption of HOD is shown in Fig. 3 for the OH + D and OD + H channels. To facilitate a fair comparison, they are given as a function of the total energy relative to the reactant asymptote. As expected, vibrational excitation in the cleaving bond always enhances the reaction probability. For example, the first excited OD stretching state, HOD(100), has a much larger dissociation probability than HOD(000) for producing OH + D. On the other hand, the first OH stretching excited state, HOD(001), enhances the OD + H channel with a comparable magnitude. Interestingly, excitation of the vibration in the non-reactive bond typically inhibits the cleavage of the reactive bond. At 1.8 eV, for example, the dissociation probability for the OH bond cleavage into the OD + H channel in the dissociative chemisorption of HOD(100) (the OD stretching mode excited) is 0.00516, much lower than 0.0153 for HOD(000). In addition, the excitation of the bending mode also enhances moderately the reaction probabilities in both channels, relative to the vibrational ground state. This is presumably due to the geometry of the transition-state which has a different H–O–H angle than free H2O.45 The vibrational enhancement effects observed here bear many similarities with those found in the gas phase H + HOD reaction.42
![]() | ||
Fig. 3 Comparison of dissociative chemisorption probabilities for HOD into the OH + D channel (upper panel) and OD + H (lower panel) for several low-lying vibrational states of HOD as a function of the total energy (relative to the energy of H2O + Cu asymptote). The probabilities are plotted in logarithmic scale in the inset. |
In Table 2, the OH/OD branching ratios are listed for the translational energies of 1.0 eV and 1.3 eV. At both energies, the HOD(100) state shows significant selectivity towards the OH product, while the HOD(001) state has a strong preference for the OD product. In Fig. 4, the dissociation probabilities are shown as a function of translational energy in the logarithmic scale. The enhancement effects of the vibrationally excited HOD are quite apparent. To make the comparison realistic, we plot in Fig. 5 the same reaction probabilities assumed Boltzmann weighting factors for these initial HOD states at 700 K, based on the nozzle temperature reported in the experiments on methane dissociative chemisorption.4 At this temperature, it is seen that the dissociative chemisorption is dominated at low translational energies by the vibrationally excited HOD in which the cleaving bond is excited.
Vibrational state | 1.0 eV | 1.3 eV | ||||
---|---|---|---|---|---|---|
OH + D | OD + H | OH/OD | OH + D | OD + H | OH/OD | |
HOD(000) | 1.03 × 10−4 | 5.98 × 10−4 | 1![]() ![]() |
1.62 × 10−2 | 1.21 × 10−2 | 1.3![]() ![]() |
HOD(010) | 8.73 × 10−3 | 2.06 × 10−2 | 1![]() ![]() |
8.34 × 10−2 | 6.91 × 10−2 | 1.2![]() ![]() |
HOD(100) | 1.04 × 10−1 | 4.84 × 10−3 | 21.5![]() ![]() |
2.91 × 10−1 | 2.22 × 10−2 | 13.1![]() ![]() |
HOD(001) | 4.65 × 10−3 | 2.47 × 10−1 | 1![]() ![]() |
2.61 × 10−2 | 3.44 × 10−1 | 1![]() ![]() |
![]() | ||
Fig. 4 Comparison of dissociative chemisorption probabilities for HOD into the OH + D (upper panel) and OD + H channels (lower panel) for several low-lying vibrational states of HOD as a function of the translational energy. |
![]() | ||
Fig. 5 Dissociative chemisorption probabilities for HOD into the OH + D (upper panel) and OD + H channels (lower panel) for several low-lying vibrational states of HOD weighted by Boltzmann factors for the initial vibrational states of HOD at 700 K, which was the nozzle temperature in a recent molecular beam experiment for methane dissociative chemisorption.4 |
Most chemical reactions can be reasonably described by the transition-state theory, in which the probability to overcome a reaction barrier is solely determined by temperature or the total energy. In transition-state theory, the system is assumed to be ergodic and all modes have the same efficacy in promoting the reaction. Mode-selective reactions, in which the barrier crossing is preferentially promoted by energy in certain vibrational coordinate(s), suggest that all forms of energy are not equal. In other words, dynamics is essential in understanding the reactivity. In this work, we demonstrate, using a high-dimensional quantum dynamical model, the strong non-statistical nature of an industrially important heterogeneous catalytic process, namely the dissociative chemisorption of water.
Specifically, quantum dynamics studies on a DFT based PES have shown strong bond selectivity in the dissociative chemisorption of HOD on Cu(111), underscoring the central role of vibrational energy in the reaction. The mode and bond selectivities in this system can be attributed to the “late” transition state and the slow IVR rate in water. Despite the large body of work on the adsorption and dissociation of water on metal surfaces,56 there has not been any experimental report on water dissociative chemisorption dynamics. As our theoretical studies have demonstrated, significant control of the dissociative chemisorption of water can be achieved with vibrational excitation. Hence, it is hoped that our work will stimulate experimental explorations, particularly those with molecular beams and internal state resolution, of the mode and bond-selectivities in this important system in heterogeneous catalysis.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c2sc21393a |
This journal is © The Royal Society of Chemistry 2013 |