The spin significance in the capture and activation of N2O by small Rh nanoparticles

https://doi.org/10.1016/j.molcata.2013.03.025Get rights and content

Highlights

  • Rhn + N2O reaction study must consider the geometry of the reactants and spin.

  • N2O reduction mechanism catalysed by Rhn clusters ends up breaking the N2single bondO bond.

  • The N2O activation never occurs through the Nsingle bondN bond rupture.

  • Most of Rh clusters capture N2O and there are some spin states that also activate it.

  • The spin states transfer enough charge from Rh to N2O to break the N2single bondO bond.

Abstract

This work presents many possible theoretical reaction pathways of N2O reduction to N2 and O on Rhn nanoparticles (n = 1–4) using density functional theory (DFT) method and the zero order regular approximation (ZORA), which explicitly considers the scalar relativistic corrections. The Rh spin multiplicity is an essential condition to dissociate N2O, because can promote or inhibit the electron back donation from the metal. Rh activates N2O by exothermic and spontaneous reactions. For each case presented, the optimized geometry adsorption site, reaction energy, spin multiplicity and Voronoi charges are calculated. On a single Rh atom in the ground and low-lying excited states, the N2O is captured only. On the Rh2 quintet ground state, N2O is also chemisorbed and dissociation occurs for the next two excited states (triplet and septet). In the Rh3 case, there are N2O adsorption and dissociation for Rh3 quartet ground state and for most excited states. Several dissociation cases take place when N2O is parallel to the plane and parallel to a Rh3 bond. On the Rh4, however, there are two optimal geometries: tetrahedral and square; for both cases there are N2O adsorption and dissociation.

Introduction

One of the real problems of the last decade is how to reduce atmospheric pollution caused by nitrogen oxides (NOx), specially the nitrous oxide (N2O). This pollutant is a major cause of ozone formation at the biosphere, and simultaneously, is the dominant substance that is depleting the ozone layer in the stratosphere [1]. At the present, the N2O amount in the atmosphere is measured and has been considered by the Montreal Protocol, but in earlier times it was considered harmless to atmosphere. For this reason, only recently has it become considered harmful and the study interest has grown. Industry has had to fabricate catalytic devices to reduce NOx emissions of stationary sources and automotive vehicles. Pt, Pd and Rh metals are commonly used in three-way catalyst (TWC) systems; Rh metal is recognized as the most appropriate metal to reduce NOx to N2, since Rh readily provides electrons for nitrogen oxide reduction. In fact, Rh has the lowest ionization energy of the three metals in the TWC (7.46 [2], 8.34 [2] and 8.96 eV [3] ionization energies for Rh, Pd and Pt atoms, respectively). N2O activation is complex, at certain temperatures N2O is an intermediate and its contents increases at temperatures below ~250 °C, i.e., up to 80% of NO is reduced to N2O when engine is warming up [4].

The need to reduce pollutant emissions has increased the interest in best experimental and theoretical knowledge of NOx reductions, especially the N2O intermediate. Many experimental and theoretical studies have been reported [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16] about the mechanisms of NOx activation reaction considering the best reduction conditions and catalytic systems.

Theoretical studies on these processes are suitable for understanding the NOx reduction produced by industries and motor vehicles, because these processes are quite complex and difficult to reproduce in laboratories using normal chemical devices. Due to current state development of the quantum chemistry methods it is also complicated to carry out some theoretical studies, but it is still possible to perform partial ab initio studies of the processes, although the transition metal atoms number is restricted. The density functional theory (DFT) represents a practical method to study systems with transition metal atoms. A further difficulty in theoretical studies is related to Rh nanoparticles, since they have many degenerate states with several structural isomers [17], [18], [19], [20], [21], [22]. Ghosh et al. [20] found that NO binds more strongly to Rh clusters than to Rh surfaces. Bae et al. [21] showed that open structures are more stable than icosahedral structures, which are normally employed. Moreover, many of these isomers are very close in energy and likely to be present in experiments, unless this are carried out at very low temperatures. In other metals, such as Au, the ground state is unable to bind and decompose N2O, whereas excited states capture and activate it [22].

This article is organized as follows: In Section 2, the computational method is described summarily; density functional theory has been used in all calculations taking into account relativistic effects included explicitly in the ZORA method, which are absolutely necessary to study transition metal clusters. In Section 3, results and discussion are shown. First, Rh clusters and their isomers with the lowest energy are shown. Subsequently, the N2O with Rhn interactions yielding to RhnN2O adduct are discussed. This is followed by the analysis of N2O adsorption and decomposition catalysed by Rhn particles with n = 1–4 on gas phase, considering the most stable Rh nanoclusters and some Rhn low-lying excited states, considering various N2O approaches to Rhn clusters with different spin multiplicities. In the cases where Rhn catalysed the N2O dissociation, the reduction mechanism and the role played by charge transfer processes were analysed through Voronoi charges. Finally, the conclusions are presented in Section 4.

Section snippets

Computational method

All calculations were performed using DFT method [23], [24] within the generalized gradient approximation (GGA) by means of the exchange-correlation functional free of empirical parameters proposed by Perdew–Burke–Ernzerhof (PBE) [25] using the ADF 2009.01 package [26]. A triple zeta basis set with two polarization functions (TZ2P) was used with the frozen core approximation for Rh atom, in which only 18 valence electrons are considered explicitly in the self-consistent-field cycles. For oxygen

Rh clusters

To explain the N2O activation catalysed by Rh small clusters, the first step is to define the catalyst geometrical factors and spin multiplicity. For small clusters, there are few experimental data due to the difficulty to determine the geometrical properties and energy states directly from small metal clusters [30]. Moreover, Rh has more degenerate states than any other transition metal; hence there is a significant dispersion in the available theoretical and experimental results.

Table 1

Conclusions

The quantum chemical theoretical PBE/ZORA method is adequate for studying Rhn + N2O reaction, because it describes quite well the reaction electronic states order considering different geometrical approaches and spin multiplicity.

The Rhn clusters calculations performed are in good agreement with calculations reported by other authors, most included in this article, particularly in spin multiplicity, optimal geometry, binding energies, charge distribution and Rhsingle bondRh distances.

According to our

Acknowledgements

The authors are indebted to the UAM-Iztapalapa Computer Centre by the computer time granted to accomplish this work in the Aitzaloa Supercomputer and with UAM-Azcapotzalco FAMA area by the computing time in the Molphys Supercomputer. R.A. is grateful to Instituto de Educación Media Superior of México City by the leave with pay grant. O.O.-N. thanks the PROMEP – México financial support project No. 22310793.

References (59)

  • P. Granger et al.
  • C.-L. Hu et al.

    Chem. Phys. Lett.

    (2007)
  • B.R. Wood et al.

    J. Catal.

    (2002)
  • F. Kapteijn et al.

    Appl. Catal. B: Environ.

    (1996)
  • R. Burch et al.

    Appl. Catal. B: Environ.

    (2002)
  • E. Poulain et al.

    Chem. Phys. Lett.

    (2010)
  • L. Wang et al.

    Chem. Phys. Lett.

    (2002)
  • R.J. Van Zee et al.

    Chem. Phys. Lett.

    (1992)
  • D. Dai et al.

    Chem. Phys. Lett.

    (1996)
  • D. Dai et al.

    Chem. Phys. Lett.

    (1992)
  • Y. Li et al.

    Surf. Sci.

    (1996)
  • F.G. Barker et al.

    Surf. Sci.

    (1973)
  • A.R. Ravishankara et al.

    Science

    (2009)
  • C.L. Callender et al.

    J. Opt. Soc. Am. B

    (1988)
  • A. Marijnissen et al.

    Phys. Rev. A

    (1995)
  • A. Yamada et al.

    Phys. Chem. Chem. Phys.

    (2012)
  • A.L. Yakovlev et al.

    Catal. Lett.

    (2001)
  • H. Orita et al.

    J. Phys. Chem. C

    (2010)
  • J.A. Ryder et al.

    J. Phys. Chem. B

    (2002)
  • T.M. Miller et al.

    J. Am. Chem. Soc.

    (1995)
  • E.J. Karlsen et al.

    J. Phys. Chem. A

    (2002)
  • P. Pannopard et al.

    ChemPhysChem

    (2012)
  • A.S. Wörz et al.

    J. Am. Chem. Soc.

    (2003)
  • V. Bertin et al.

    Int. J. Quantum Chem.

    (2010)
  • D. Harding et al.

    Phys. Chem. Chem. Phys.

    (2007)
  • P. Ghosh et al.

    J. Chem. Phys.

    (2008)
  • Y.-C. Bae et al.

    Phys. Rev. B

    (2004)
  • O. Olvera-Neria et al.

    J. Chem. Phys.

    (2010)
  • P. Hohenberg et al.

    Phys. Rev.

    (1964)
  • Cited by (10)

    • The π back-donation influence in CO oxidation on small and oxidized Au–Ag clusters

      2022, Chemical Physics
      Citation Excerpt :

      The exchange–correlation energy was calculated with the Perdew-Burke-Ernzerhof (PBE) functional [35] free of empirical parameters. The PBE functional adequately describes reactions that involve transition metal atoms [36–39]. We have used the Van Wüllen [40] implementation of the ZORA method in the ORCA3.0.1 program [41] to consider the scalar relativistic effects [30–32], which must be included for chemical reactions involve 4d or 5d transition metals.

    • Dissociation of N<inf>2</inf>O promoted by Rh<inf>6</inf> clusters. A ZORA/DFT/PBE study

      2015, Journal of Molecular Catalysis A: Chemical
      Citation Excerpt :

      As seen in Table 3B, without the inclusion of the relativistic effects, through ZORA, in the computational procedure, Rh6 is not able to reduce N2O. Considering all possible active sites of Rh6, approximately 47, they are close to the sum of the adsorption sites obtained for all the smaller Rhn particles previously studied [29]. That is, the number of active sites increases with the number of Rh atoms.

    View all citing articles on Scopus
    View full text