Next Article in Journal
Semisynthesis and Biological Evaluation of Oleanolic Acid 3-O-β-d-Glucuronopyranoside Derivatives for Protecting H9c2 Cardiomyoblasts against H2O2-Induced Injury
Next Article in Special Issue
Pd-Catalyzed, Highly Selective C(sp2)-Br Bond Coupling Reactions of o-(or m-, or p-) Chloromethyl Bromobenzene with Arylboronic Acids
Previous Article in Journal
New Thiazoline-Tetralin Derivatives and Biological Activity Evaluation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monitoring Reaction Paths Using Vibrational Spectroscopies: The Case of the Dehydrogenation of Propane toward Propylene on Pd-Doped Cu(111) Surface

1
School of Chemistry and Materials Science, University of Science and Technology of China, iChEM (Collaborative Innovation Center of Chemistry for Energy Materials), Hefei 230026, China
2
Department of Physics and Institute of Theoretical Physics and Astrophysics, Xiamen University, Xiamen 361005, China
*
Author to whom correspondence should be addressed.
Molecules 2018, 23(1), 126; https://doi.org/10.3390/molecules23010126
Submission received: 23 November 2017 / Revised: 4 January 2018 / Accepted: 5 January 2018 / Published: 10 January 2018
(This article belongs to the Special Issue Palladium Catalysts 2018)

Abstract

:
Monitoring reaction paths is not only a fundamental scientific issue but also helps us to understand and optimize the catalytic process. Infrared (IR) and Raman spectroscopies are powerful tools for detecting particular molecules or intermediate products as a result of their ability to provide the molecular “finger-print”. However, theoretical modeling for the vibrational spectra of molecular adsorbates on metallic surfaces is a long-standing challenge, because accurate descriptions of the electronic structure for both the metallic substrates and adsorbates are required. In the present work, we applied a quasi-analytical IR and Raman simulation method to monitor the dehydrogenation of propane towards propylene on a Pd-doped Cu(111) surface in real-time. Different Pd ensembles were used to construct the single-atom catalyst (SAC). We found that the number of sublayer Pd atoms could only affect the intensity of the peak rather than the peak position on the vibrational spectra. However, with the dehydrogenation reaction proceeding, both IR and Raman spectra were changed greatly, which indicates that every reaction step can be distinguished from the point of view of vibrational spectroscopies. Additionally, we found that the catalytic process, which starts from different initial states, shows different spectral profiles. The present results suggest that the vibrational spectroscopies obtained by the high-precision simulations pave the way for identifying different catalytic reaction paths.

Graphical Abstract

1. Introduction

Studying reaction paths is not only a fundamental scientific issue, but it also helps in understanding and optimizing the catalytic process. Theoretically, transition state theory is usually employed to evaluate the activation energies; then the most likely reaction path can be determined accordingly. It is noted that the theoretical evaluations of activation energies are highly reliant on the simulation environments [1,2,3,4,5,6,7,8]. As a result, the relative values of activation energies rather than the absolute values are more useful in predicting reaction paths. Compared to the activation energies, spectroscopies usually provide more reliable information. Many technologies, such as nuclear magnetic resonance (NMR), X-ray photoelectron spectroscopy (XPS) and vibrational spectroscopies have been employed to detect specific intermediate products. On the basis of the detection of the intermediate product, the reaction path can be determined accordingly. Among these spectroscopic technologies, infrared (IR) and Raman spectra have advantages in detecting the intermediate products because of their in situ operability [1,2,3,4,9,10,11,12,13].
Up to now, there have been two approaches to model the adsorption of molecules on metallic surfaces and to simulate the corresponding IR and Raman spectra at the ab initio level: the cluster model and the periodic boundary condition (PBC) model. It is well known that the PBC model could provide a more reliable electronic structure of metallic substrates and nanoparticles compared to the cluster model. Therefore, many efforts have been made to develop IR and Raman simulation approaches using PBC models [14,15,16]. For example, Ding et al. developed a method that combines the geometrical optimization of the PBC model and Raman calculation of the cluster model [17]. They have tested the method by simulating the Raman spectra CO and ethylene adsorbed on Rh(111). Many efforts have also been made by our group in the previous works [18,19,20]. Our quasi-analytical method has also been applied in many specific systems, and the accuracy and stability have been tested [18,19,20].
Light olefins are important basic raw materials in the chemical industry. With their growing demand in the global chemical market, the catalytic production of olefins has attracted increased interest in the past decades [21]. It has been widely adopted that the dehydrogenation of light alkanes is a dominant commercial means to produce olefins. However, the highly endothermic process limits its industrial efficiency. Therefore, efficient catalysts are deeply desired. Recently, we designed a kind of single-atom catalyst (SAC) based on a Pd-doped Cu(111) surface. The theoretical simulations show that the embedded Pd atoms significantly improve the catalytic reactivity, whereas the presence of the Cu surface is beneficial for the diffusions of the detached H atoms. Therefore, we conclude that the high selectivity toward propylene dehydrogenation is realized on the Pd-doped Cu(111) surface. However, if these designed catalysts can be used in reality, how can we know the true reaction path from experiments at the in situ level?
In the present work, we attempt to apply our quasi-analytical method to systematically investigate the IR and Raman spectra of a whole reaction process. Specifically, the dehydrogenation of propane on Pd-doped Cu(111) is studied. By studying the changes in the IR and Raman spectra, the reaction paths are monitored in real-time. Most of all, different reaction paths are characterized, which provide us with a potential application of vibrational spectroscopies in tracing reaction paths.

2. Results

To investigate the effect of multiple doping on the catalytic performance, we designed four different Pd-doped Cu(111) surfaces. As shown in Figure 1a, the number of doped Pd atoms ranged from one to four. It is noted that the most stable configuration of propane adsorbed on the surface has been proven to be a parallel state (shown in Figure 1c, labelled as propane|| hereafter), which agrees well with the previous works [22,23,24,25,26]. On the other hand, Figure 1b shows the second stable configuration, in which a perpendicular adsorption is adopted. Our previous studies show that both parallel and perpendicular adsorbed propane could be catalyzed to propylene [27]. To explore the transition states of the propane dehydrogenation process, we employed the widely used climbing-image nudged elastic band (CI-NEB) method [28]. Frequency calculations were also performed to characterize the obtained transition-state structures, and these intermediates have only one imaginary frequency according to our previous study [27]. The corresponding configurations for the initial reactants, transition states and final products for the first dehydrogenation reaction are shown in Figure 1b,c. We can see that the main difference between the two reaction paths is the order of the dehydrogenation process. For the parallel adsorbed configuration, the first dehydrogenation occurs on the methylene group, and the second occurs on the methyl group. However, for the perpendicular configuration, the order of the two dehydrogenation processes is reversed. For the sake of discussion, we label these stationary points as propane|| (propane||), 1-propyl (2-propyl), 1-propyl-diff (2-propyl-diff), and 1-propylene (2-propylene). In particular, 1-propyl-diff (2-propyl-diff) indicates an extra stable state in which the first detached proton is diffused. Transition-state calculations show that the first dehydrogenation process needs the highest activation energies and that the parallel configuration shows better catalytic performance than the perpendicular configuration. To be specific, the activation energies for the parallel configuration with the 1–4 Pd-doped Cu(111) surfaces were found to be reduced by 5.1, 6.2, 8.0 and 9.7 kcal/mol compared to the pure Cu(111) surface, resulting in 20–400 times increases in the reaction rates [27].
In Table 1, we tabulate the predicted energy barriers for the first C–H bond cleavage from the methyl and methylene groups of propane on the pure and Pd-doped Cu(111) surfaces. From Table 1, one can see that the calculated activation energies varied only a little for these representative reaction paths according to our previous study. Taking the propane dehydrogenation over the single Pd-doped Cu(111) surface for instance, the first C–H activation from the CH3 group was 32.8 kcal/mol, while the dehydrogenation of CH3 required an energy barrier of 31.6 kcal/mol. The energy difference between these two reaction paths was only 1.2 kcal/mol; this tiny difference makes it very hard to determine the reaction process in the true reaction from transition-state calculations. Therefore, it is highly desirable to find an appropriate way to determine the reaction paths. Because the IR and Raman spectroscopies can provide the molecular “finger-print”, investigating the spectral changes during the whole reaction process became an ideal solution for identifying the true reaction path.
In this work, we employed our developed quasi-analytical method [18,19] to study the vibrational spectroscopies of the propane dehydrogenation reaction. Firstly, we evaluate the influence of the number of doped Pd atoms on the vibrational spectra. Figure 2 shows the IR adsorption and Raman scattering spectra for propane and propane|| adsorbed on the Pd-doped Cu(111) surfaces. We can see that the peak positions in the IR and Raman spectra were almost independent of the number of doped Pd atoms. For propane, the three most intense peaks in the IR spectra were located at 867, 1128 and 1324 cm−1, corresponding to C–C–C stretching, methyl CH3 rocking and methyl CH3 deforming, as shown in Table 2. Interestingly, with more Pd-doped atoms in the surface, the peaks located at around 867 and 1324 cm−1 became weaker. For the Raman spectra, there was only one dominant peak located at 867 cm−1 for propane. Compared to propane, propane|| showed somewhat different vibrational spectra. To be specific, four dominant peaks were found in the IR spectra, located at around 710, 863, 1136 and 1445 cm−1. These four peaks are attributed to methylene CH2 rocking, C–C–C stretching, methyl CH3 rocking and methyl CH3 deforming. Meanwhile, there were only two dominant peaks in the Raman spectra, located at 863 and 1136 cm−1. It is noted that for any other intermediate or final products, the IR and Raman spectra varied only a little, no matter how many Pd atoms were doped (all vibrational spectra are provided in Supplementary Information).
Special attention is paid to the change in the vibrational spectra during the reaction process. As shown in Figure 3, the IR spectra changed greatly as the reaction proceeded. Specifically, there were three intense peaks for propane, located at 867, 1128 and 1324 cm−1. However, after the first hydrogen dehydrogenated from the methyl group, only two dominant peaks (1013 and 1089 cm−1) appeared, no matter whether the detached proton was diffused or not. These two peaks are attributed to C–C–C stretching and dehydrogenated methyl twisting. It is noted that due to the dehydrogenation of CH3, the symmetry of the adsorbate changed. As a result, the peaks that were related to methyl vanished. When the two dehydrogenations finished and propylene formed, three peaks (877 and 937 and 1030 cm−1) dominated the IR spectrum, corresponding to dehydrogenated methylene C–H wagging, dehydrogenated methyl CH2 wagging and dehydrogenated methyl CH2 rocking. We can see from Figure 3a that the methyl- and methylene-related peaks diminished and the dehydrogenated groups-related peaks appeared.
A similar phenomenon can also be observed for the reaction path initiating from propane||. As shown in Figure 3b, four strong peaks dominated the IR spectra for propane||. However, after the first hydrogen dehydrogenated from the methylene group, there was only one predominant peak located at 1082 cm−1, which is attributed to methyl CH3 rocking. The peak located at 710 cm−1 (attributed to methylene CH2 rocking) diminished because the dehydrogenation occurred on methylene. The IR spectra varied only a little before and after the first detached proton diffused, although the CH3 rocking mode showed a 20 cm−1 red shift. Finally, the CH3 rocking mode become inactive in the IR spectrum after the second H dehydrogenated from the methyl group. It is very interesting to see that the final products obtained from the two reaction paths showed very similar IR spectra, indicating similarity in the geometries.
Several interesting phenomena were observed during the whole catalytic process. For example, the C–C–C stretching mode showed a blue shift, that is, 867, 888, 891 and 937 cm−1 for propane, 1-propyl, 1-propyl-diff, 1-propylene, and 863, 871, 889, and 937 cm−1 for propane||, 2-propyl, 2-propyl-diff, and 2-propylene, respectively. On the other hand, the CH3 rocking mode showed red shifts during the dehydrogenation. Particularly, 1128, 1114, 1088 and 1031 cm−1 were observed for propane, 1-propyl, 1-propyl-diff and 1-propylene during the first reaction path. For the second reaction path, these vibrational modes were located at 1136, 1082, 1062, and 1030 cm−1 for propane||, 2-propyl, 2-propyl-diff, and 2-propylene, respectively.
The ability of Raman spectroscopy to monitor reaction processes has also been investigated. As shown in Figure 4, only one peak located at 867 cm−1 was found in the Raman spectrum for propane, corresponding to C–C–C stretching. After the first proton detached from the methyl group, the most intense peak in the propyl species was located at 1010 cm−1, which is attributed to dehydrogenated methyl (–CH2) wagging. The diminishing and appearance of Raman peaks revealed the symmetry changes. Additionally, the proton diffusion hardly affected the Raman spectra. Finally, when the second proton dehydrogenated from 1-propyl, only one peak located at 1615 cm−1 dominated the Raman spectra, which is attributed to dehydrogenated C=C stretching. For propane||, two peaks appeared in the Raman spectra, 863 and 1136 cm−1, which are attributed to C–C–C stretching and methyl –CH3 rocking. After one proton detached from the methylene group, a small peak split from the main peak. Moreover, a red shift appeared for the –CH3 rocking peak, which shifted from 1136 to 1082 cm−1. Interestingly, the final product from the two different reaction paths showed similar Raman spectra. From Figure 4, we can also conclude that using Raman spectra, every dehydrogenation step can be easily monitored, and the two different reaction paths can also be distinguished.

3. Computational Methods and Details

The spin-polarized DFT calculations were performed using the Vienna ab initio simulation package (VASP) [29,30]. The interactions between the core ions and the valence electrons were described by using the projector augmented wave (PAW) pseudopotentials [31]. The exchange correlation effects were described by the GGA-PBE functional [32]. A plane-wave energy cutoff of 400 eV was used in the calculations. To mimic the flat (111) surface, a four-layer slab of a (5 × 5) unit cell and a 10 Å thick vacuum region were employed here. The Brillouin zone sampling was carried out by using the (5 × 5 × 1) Monkhorst-Pack grids. During the optimization, the lowermost two layers were fixed, while the atoms in the uppermost two layers were allowed to relax until the maximum force became less than 0.02 eV/Å.
After the optimization of the surface system, we fixed all metal atoms and allowed all adsorbate atoms to relax to obtain the normal modes using density functional perturbation theory (DFPT) [33,34]. Then the IR intensity was expressed as
I α k = i l i k Z α i m i
and the Raman tensors were calculated as
R α β k = i l i k m i Z α i F β
where α and β are Cartesian coordinates, l is the Cartesian displacement vector, m is the mass of the atom, and Z α i are the Born effective charges, which can be computed analytically. The derivative in Equation (2) is computed by the first-order finite-difference method, that is,
Z α i F β = Z α i ( Δ F β ) Z α i ( Δ F β ) 2 Δ F β
It is noted that the scaling factors are usually needed to make the simulated frequencies more accurate and more comparable to the experimental observations. However, we have noticed that Martin et al. have proven that the scaling factor for the PBE functional is 0.98–0.99 [35]. Furthermore, in the present work, the relative positions and intensities rather than the absolute values played an important role in monitoring the reaction path. As a result, no scaling factor was used in the present work. With the calculated Raman tensor, the differential Raman cross-section could then be calculated using the standard expression [36,37,38]:
d σ d Ω k = π 2 ϵ 0 2 ν ˜ in ν ˜ k 4 h 8 π 2 c ν ˜ k | ε i · R ^ · ε s | 2 1 1 exp h c ν ˜ k / k B T
where ϵ 0 is the electric permittivity of free space, ν ˜ in is the wavenumber of incident light (set to be 514 nm), ν ˜ k is the wavenumber of the corresponding frequency, h is the Planck constant, c is the speed of light in a vacuum, k B is the Boltzmann constant, and T is the temperature (set to be 298.15 K). Meanwhile, all the calculated Raman cross-sections were convoluted with a Lorentzian function with the full width at half maximum of 15 cm−1. Only R z z k needed to be taken into account, in view of the electric field being only along the z direction [16], from results based on PBC models.

4. Conclusions

In the present work, we propose a method to monitor reaction paths. We systematically studied the IR adsorption and Raman scattering spectra of all reactants, intermediates and final products for the dehydrogenation of propane towards propylene. Two different reaction paths initiating from perpendicular and parallel adsorbed propane were considered. The first-principle calculations of transition-state calculations could hardly distinguish between the two reaction paths. However, using vibrational spectra, every step during the reaction path could be distinguished, and the two different reaction paths could also be easily identified. As a result, the predicted IR and Raman spectra of these reaction intermediates and products provide a basis for their experimental characterization and for the determination of reaction paths.

Supplementary Materials

The supplementary materials are available online.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (21703223 and 21703178), the Fundamental Research Funds for the Central Universities (WK2060030026) and the China Postdoctoral Science Foundation (2017M612078). The Swedish National Infrastructure for Computing (SNIC) is acknowledged for computer time.

Author Contributions

Wei Hu suggested the project and carried out the theoretical simulations. Wei Hu and Xinrui Cao discussed the results. Wei Hu wrote the paper, and Xinrui Cao commented on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, R.; Zhang, Y.; Dong, Z.; Jiang, S.; Zhang, C.; Chen, L.; Zhang, L.; Liao, Y.; Aizpurua, J.; Luo, Y.; et al. Chemical Mapping of a Single Molecule by Plasmon-enhanced Raman Scattering. Nature 2013, 498, 82–86. [Google Scholar] [CrossRef] [PubMed]
  2. Aikens, C.M.; Madison, L.R.; Schatz, G.C. Raman spectroscopy: The effect of field gradient on SERS. Nat. Photonics 2013, 7, 508–510. [Google Scholar] [CrossRef]
  3. Cong, S.; Yuan, Y.; Chen, Z.; Hou, J.; Yang, M.; Su, Y.; Zhang, Y.; Li, L.; Li, Q.; Geng, F. Noble metal-comparable SERS enhancement from semiconducting metal oxides by making oxygen vacancies. Nat. Commun. 2015, 6, 7800. [Google Scholar] [CrossRef] [PubMed]
  4. Palonpon, A.F.; Ando, J.; Yamakoshi, H.; Dodo, K.; Sodeoka, M.; Kawata, S.; Fujita, K. Raman and SERS microscopy for molecular imaging of live cells. Nat. Protoc. 2013, 8, 677–692. [Google Scholar] [CrossRef] [PubMed]
  5. Pramanik, A.; Chavva, S.R.; Nellore, V.; Priya, B.; May, K.; Matthew, T.; Jones, S.; Vangara, A.; Ray, P.C. Development of a SERS Probe for Selective Detection of Healthy Prostate and Malignant Prostate Cancer Cells Using ZnII. Chem. Asian J. 2017, 12, 665–672. [Google Scholar] [CrossRef] [PubMed]
  6. Jia, P.; Chang, J.; Wang, J.; Zhang, P.; Cao, B.; Geng, Y.; Wang, X.; Pan, K. Fabrication and Formation Mechanism of Ag Nanoplate-Decorated Nanofiber Mats and Their Application in SERS. Chem. Asian J. 2016, 11, 86–92. [Google Scholar] [CrossRef] [PubMed]
  7. Hu, Y.; Zhao, T.; Zhu, P.; Zhu, Y.; Liang, X.; Sun, R.; Wong, C.P. Tailoring Size and Coverage Density of Silver Nanoparticles on Monodispersed Polymer Spheres as Highly Sensitive SERS Substrates. Chem. Asian J. 2016, 11, 2428–2435. [Google Scholar] [CrossRef] [PubMed]
  8. Philips, D.S.; Sreejith, S.; He, T.; Menon, N.V.; Anees, P.; Mathew, J.; Sajikumar, S.; Kang, Y.; Stuparu, M.C.; Sun, H.; et al. A Three-Photon Active Organic Fluorophore for Deep Tissue Ratiometric Imaging of Intracellular Divalent Zinc. Chem. Asian J. 2016, 11, 1523–1527. [Google Scholar] [CrossRef] [PubMed]
  9. Auer, B.; Skinner, J. IR and Raman spectra of liquid water: Theory and interpretation. J. Chem. Phys. 2008, 128, 224511. [Google Scholar] [CrossRef] [PubMed]
  10. Albrecht, A.C. On the theory of Raman intensities. J. Chem. Phys. 1961, 34, 1476–1484. [Google Scholar] [CrossRef]
  11. Kurouski, D.; Zaleski, S.; Casadio, F.; Van Duyne, R.P.; Shah, N.C. Tip-enhanced Raman Spectroscopy (TERS) for in situ Identification of Indigo and Iron Gall Ink on Paper. J. Am. Chem. Soc. 2014, 136, 8677–8684. [Google Scholar] [CrossRef] [PubMed]
  12. Molesky, B.P.; Guo, Z.; Cheshire, T.P.; Moran, A.M. Perspective: Two-dimensional resonance Raman spectroscopy. J. Chem. Phys. 2016, 145, 180901. [Google Scholar] [CrossRef] [PubMed]
  13. Girard, A.; Lermé, J.; Gehan, H.; Margueritat, J.; Mermet, A. Mechanisms of resonant low frequency Raman scattering from metallic nanoparticle Lamb modes. J. Chem. Phys. 2017, 146, 194201. [Google Scholar] [CrossRef] [PubMed]
  14. Karhánek, D.; Bučko, T.; Hafner, J. A Density Functional Study of the Adsorption of Methane-thiol on the (111) Surfaces of the Ni-group Metals: I. Molecular and Dissociative Adsorption. J. Phys. Condens. Matter 2010, 22, 265005. [Google Scholar] [CrossRef] [PubMed]
  15. Karhánek, D.; Bučko, T.; Hafner, J. A Density-functional Study of the Adsorption of Methane-thiol on the (111) Surfaces of the Ni-group Metals: II. Vibrational Spectroscopy. J. Phys. Condens. Matter 2010, 22, 265006. [Google Scholar] [CrossRef] [PubMed]
  16. Zayak, A.T.; Hu, Y.S.; Choo, H.; Bokor, J.; Cabrini, S.; Schuck, P.J.; Neaton, J.B. Chemical Raman Enhancement of Organic Adsorbates on Metal Surfaces. Phys. Rev. Lett. 2011, 106, 083003. [Google Scholar] [CrossRef] [PubMed]
  17. Ding, Z.B.; Tommasini, M.; Maestri, M. First-principles simulation of Raman Spectra of Adsorbates on Metal Surfaces. ChemPlusChem 2017, 82, 924–932. [Google Scholar] [CrossRef]
  18. Hu, W.; Duan, S.; Zhang, G.; Ma, Y.; Tian, G.; Luo, Y. Quasi-Analytical Approach for Modeling of Surface-Enhanced Raman Scattering. J. Phys. Chem. C 2015, 119, 28992–28998. [Google Scholar] [CrossRef]
  19. Hu, W.; Duan, S.; Luo, Y. Theoretical modeling of surface and tip-enhanced Raman spectroscopies. WIREs Comput. Mol. Sci. 2017, 7, e1293. [Google Scholar] [CrossRef]
  20. Hu, W.; Duan, S.; Zhang, Y.; Ren, H.; Jiang, J.; Luo, Y. Identifying the structure of 4-chlorophenyl isocyanide adsorbed on Au(111) and Pt(111) surfaces by first-principles simulations of Raman spectra. Phys. Chem. Chem. Phys. 2017, 19, 32389–32397. [Google Scholar]
  21. Sattler, J.J.H.B.; Ruiz-Martinez, J.; Santillan-Jimenez, E.; Weckhuysen, B.M. Catalytic Dehydrogenation of Light Alkanes on Metals and Metal Oxides. Chem. Rev. 2014, 114, 10613–10653. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, M.L.; Zhu, Y.A.; Zhou, X.G.; Sui, Z.J.; Chen, D. First-Principles Calculations of Propane Dehydrogenation over PtSn Catalysts. ACS Catal. 2012, 2, 1247–1258. [Google Scholar]
  23. Lo, J.M.H.; Premji, Z.A.; Ziegler, T.; Clark, P.D. First-Principles Investigation of Selective Oxidation of Propane on Clean and Sulfided V2O5 (010) Surfaces. J. Phys. Chem. C 2013, 117, 11258–11274. [Google Scholar] [CrossRef]
  24. Yang, M.L.; Zhu, Y.A.; Fan, C.; Sui, Z.J.; Chen, D.; Zhou, X.G. DFT Study of Propane Dehydrogenation on Pt Catalyst: Effects of Step Sites. Phys. Chem. Chem. Phys. 2011, 13, 3257–3267. [Google Scholar] [CrossRef] [PubMed]
  25. Vajda, S.; Pellin, M.J.; Greeley, J.P.; Marshall, C.L.; Curtiss, L.A.; Ballentine, G.A.; Elam, J.W.; Catillon-Mucherie, S.; Redfern, P.C.; Mehmood, F.; et al. Subnanometre Platinum Clusters as Highly Active and Selective Catalysts for the Oxidative Dehydrogenation of Propane. Nat. Mater. 2009, 8, 213–216. [Google Scholar] [CrossRef] [PubMed]
  26. Fu, H.; Liu, Z.P.; Li, Z.H.; Wang, W.N.; Fan, K.N. Periodic Density Functional Theory Study of Propane Oxidative Dehydrogenation over V2O5 (001) Surface. J. Am. Chem. Soc. 2006, 128, 11114–11123. [Google Scholar] [CrossRef] [PubMed]
  27. Cao, X. Insight into mechanism and selectivity of propane dehydrogenation over the Pd-doped Cu(111) surface. RSC Adv. 2016, 6, 65524–65532. [Google Scholar] [CrossRef]
  28. Henkelman, G.; Uberuaga, B.P.; Jónsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901–9904. [Google Scholar] [CrossRef]
  29. Kresse, G.; Furthmüller, J. Efficiency of Ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. J. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  30. Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B. 1996, 54, 11169. [Google Scholar] [CrossRef]
  31. Blöchl, P.E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef]
  32. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  33. Gajdoš, M.; Hummer, K.; Kresse, G.; Furthmüller, J.; Bechstedt, F. Linear Optical Properties in the Projector-augmented Wave Methodology. Phys. Rev. B 2006, 73, 045112. [Google Scholar] [CrossRef]
  34. Sternheimer, R. Electronic Polarizabilities of Ions from the Hartree-Fock Wave Functions. Phys. Rev. 1954, 96, 951. [Google Scholar] [CrossRef]
  35. Kesharwani, M.K.; Brauer, B.; Martin, J.M. Frequency and zero-point vibrational energy scale factors for double-hybrid density functionals (and other selected methods): Can anharmonic force fields be avoided? J. Phys. Chem. A 2014, 119, 1701–1714. [Google Scholar] [CrossRef] [PubMed]
  36. Neugebauer, J.; Reiher, M.; Kind, C.; Hess, B.A. Quantum Chemical Calculation of Vibrational Spectra of Large Molecules—Raman and IR Spectra for Buckminsterfullerene. J. Computat. Chem. 2002, 23, 895–910. [Google Scholar] [CrossRef] [PubMed]
  37. Long, D.A. The Raman effect : A Unified Treatment of the Theory of Raman Scattering by Molecules; Wiley Online Library: Hoboken, NJ, USA, 2002; Volume 1. [Google Scholar]
  38. Le Ru, E.; Etchegoin, P. Principles of Surface-Enhanced Raman Spectroscopy: And Related Plasmonic Effects; Elsevier: Amsterdam, The Netherlands, 2008. [Google Scholar]
Sample Availability: Not available.
Figure 1. (a) Different Cu(111) surfaces with one to four doped Pd atoms; (b,c) Geometries of the initial reactants and intermediate and final products during the dehydrogenation process of the perpendicular and parallel adsorbed propane towards propylene.
Figure 1. (a) Different Cu(111) surfaces with one to four doped Pd atoms; (b,c) Geometries of the initial reactants and intermediate and final products during the dehydrogenation process of the perpendicular and parallel adsorbed propane towards propylene.
Molecules 23 00126 g001
Figure 2. Infrared (IR) adsorption and Raman scattering spectra of propane adsorbed on the Pd-doped Cu(111) surface; (a,b) indicate two adsorption configurations: perpendicular and parallel adsorbed propane.
Figure 2. Infrared (IR) adsorption and Raman scattering spectra of propane adsorbed on the Pd-doped Cu(111) surface; (a,b) indicate two adsorption configurations: perpendicular and parallel adsorbed propane.
Molecules 23 00126 g002
Figure 3. Infrared (IR) adsorption spectra of all reactants, intermediate and final products, during the dehydrogenation of propane towards propylene on the Pd-doped Cu(111) surface; (a,b) indicate the two reaction paths initiating from parallel and perpendicular adsorbed propane.
Figure 3. Infrared (IR) adsorption spectra of all reactants, intermediate and final products, during the dehydrogenation of propane towards propylene on the Pd-doped Cu(111) surface; (a,b) indicate the two reaction paths initiating from parallel and perpendicular adsorbed propane.
Molecules 23 00126 g003
Figure 4. Raman scattering spectra of all intermediate products during the dehydrogenation of propane towards propylene on the Pd-doped Cu(111) surface; (a,b) represent the two reaction paths initiating from perpendicular and parallel adsorbed propane.
Figure 4. Raman scattering spectra of all intermediate products during the dehydrogenation of propane towards propylene on the Pd-doped Cu(111) surface; (a,b) represent the two reaction paths initiating from perpendicular and parallel adsorbed propane.
Molecules 23 00126 g004
Table 1. Energy barriers for the first C–H bond cleavage from methyl (path 1) and methylene group (path 2) of propane on pure and Pd-doped Cu(111) surfaces [27].
Table 1. Energy barriers for the first C–H bond cleavage from methyl (path 1) and methylene group (path 2) of propane on pure and Pd-doped Cu(111) surfaces [27].
Surfacesabcde
Path 1 (kcal/mol)37.932.831.729.928.2
Path 2 (kcal/mol)36.731.629.628.626.3
Table 2. Calculated infrared (IR) and Raman frequencies assignment of propane (propane||), 1-propyl (2-propyl), 1-propyl-diff (2-propyl-diff), and 1-propylene (2-propylene). Only the surface doped by one Pd atom is listed here. Modes 1–7 correspond to methylene –CH2– rocking, dehydrogenated methyl –CH2 wagging, C–C–C stretching, methylene –CH2 wagging, dehydrogenated methyl –CH2 twisting, methyl –CH3 rocking and methyl –CH3 deforming.
Table 2. Calculated infrared (IR) and Raman frequencies assignment of propane (propane||), 1-propyl (2-propyl), 1-propyl-diff (2-propyl-diff), and 1-propylene (2-propylene). Only the surface doped by one Pd atom is listed here. Modes 1–7 correspond to methylene –CH2– rocking, dehydrogenated methyl –CH2 wagging, C–C–C stretching, methylene –CH2 wagging, dehydrogenated methyl –CH2 twisting, methyl –CH3 rocking and methyl –CH3 deforming.
Mode1234567
Propane729881867101011281324
1-Propyl7458948881013108911141359
1-Propyl-diff7728898911006108010881389
1-Propylene876937108810311400
Propane||71086311361445
2-Propyl87110821443
2-Propyl-diff88910621449
2-Propylene877937108710301444

Share and Cite

MDPI and ACS Style

Hu, W.; Cao, X. Monitoring Reaction Paths Using Vibrational Spectroscopies: The Case of the Dehydrogenation of Propane toward Propylene on Pd-Doped Cu(111) Surface. Molecules 2018, 23, 126. https://doi.org/10.3390/molecules23010126

AMA Style

Hu W, Cao X. Monitoring Reaction Paths Using Vibrational Spectroscopies: The Case of the Dehydrogenation of Propane toward Propylene on Pd-Doped Cu(111) Surface. Molecules. 2018; 23(1):126. https://doi.org/10.3390/molecules23010126

Chicago/Turabian Style

Hu, Wei, and Xinrui Cao. 2018. "Monitoring Reaction Paths Using Vibrational Spectroscopies: The Case of the Dehydrogenation of Propane toward Propylene on Pd-Doped Cu(111) Surface" Molecules 23, no. 1: 126. https://doi.org/10.3390/molecules23010126

APA Style

Hu, W., & Cao, X. (2018). Monitoring Reaction Paths Using Vibrational Spectroscopies: The Case of the Dehydrogenation of Propane toward Propylene on Pd-Doped Cu(111) Surface. Molecules, 23(1), 126. https://doi.org/10.3390/molecules23010126

Article Metrics

Back to TopTop