Abstract

Energy profiles linking the reactants M + SO2 (M = V(4F), Nb(6D;4F), and Ta(4F)) with the products observed for these reactions under matrix-isolation conditions, mainly the oxide complex OV(η2-SO) and the sulfide oxides SVO2, SNbO2, and STaO2, have been obtained from DFT and CASSCF-MRMP2 calculations. For each of these interactions, the radical fragments MO + SO can be reached from the lowest-lying quadruplet electronic states of the reactants. As the quadruplet and doublet radical asymptotes that vary only in the spin of the unpaired parallel electrons of the nonmetallic fragment are degenerated, a second reaction leading to the rebounding of the radical fragments can take place through both multiplicity channels. Reaction along the doublet pathway leads in each case to the most stable structure for the oxide SMO2. For the vanadium interaction, recombination of the radical species through the quadruplet channel explains for the oxide product OV(η2-SO).

1. Introduction

Besides being an uncontrollable air contaminant, sulfur dioxide is also a poison in many catalytic processes [13]. The study of the activation of the S-O bond by transition metal containing systems has been one of main lines followed to explore the factors that could favor the decomposition of this molecule. Most of these studies have been focused in analyzing the adsorption and decomposition processes of SO2 in pure metals, alloys, and metal oxides [49]. Kinetic investigations have also been carried out on the reactions of this molecule with different transition metal atoms. Results emerging from those studies show that SO2 is more reactive than other gases (O2, SO2, CO2, N2, N2O, and NO) in this kind of reactions [1014].

The interaction of SO2 with laser-ablated transition metal atoms in groups IV–VI of the periodic table has been investigated by IR spectroscopy under matrix-isolation conditions [1517]. The main products observed for these interactions are the oxide complex OM(η2-SO) and/or the sulfide oxide SMO2. Importantly, for most of these reactions, the spin multiplicity assigned to the products does not match the spin corresponding to the ground state of the reactants. For example, the frequencies detected in the IR-matrix spectra for the products emerging from the reactions M + SO2 (M = V(4F), Nb(6D), and Ta(4F)), mainly the complex OV(η2-SO) and the sulfide oxides SMO2 (M = V, Nb, Ta), fix better those determined through DFT calculations for the corresponding doublet structures [17]. Liu et al. rationalized the different spin multiplicity between the reactants and the products by invoking interactions between electronic states of different spin through intersystem-crossings. For instance, they proposed that the vanadium reaction undergoes through an intersystem crossing between the potential energy curves emerging from the quadruplet ground state and the first excited doublet state of the reactants. In our opinion, the use of spin-flip models could be inadequate for the theoretical description of reactions in which only rather light elements participate. For these reactions, it is not expected that the relativistic effects that lead to interactions of electronic states of different spin be relevant. In addition, doublet electronic states do not appear among the low-lying states for the vanadium atom, so it is unlikely that they take place in reactions proceeding under cryogenic conditions [18].

Recently, we studied theoretically the interactions of acetonitrile with group V transition metal atoms [19]. In accord with the results in reference 19, the products determined under matrix-isolation conditions for these reactions can be explained in terms of a nonspin flip scheme consisting of two sequential reactions involving the radical species M-NC + CH3. For the niobium and tantalum reactions, these radical fragments are energetically accessible; thus, they could be reached from the ground state of the reactants. As the radical moieties remain caged at the matrix, they can recombine themselves in a second reaction to yield the inserted structure CH3-M-NC. Like the quadruplet and doublet asymptotes M-NC + CH3 that vary only in the spin of the methyl group are degenerate, the rebounding reaction can take place along both multiplicity channels. Once the inserted structure CH3-M-NC has been formed, the low-multiplicity channel leads to the products detected for these interactions by migration of one or two hydrogen atoms toward the metal atom. According to this scheme, the lack of inserted structures in the product distribution determined for the vanadium reaction could be rationalized in terms of the high energy found for the radical species V-NC + CH3.

In this contribution, we use a similar two-step reaction scheme to describe the interactions of the SO2 molecule with the same group of transition metal atoms. It is shown that the results obtained for each of the investigated reactions allow to rationalize the product distributions observed for them in the IR-matrix spectra, without invoking interactions between electronic states of different spin multiplicity.

2. Computational Details

For all the investigated reactions, DFT-B3LYP-Def2-TZVP and CASSCF-MRMP2 calculations were performed to obtain the energy profiles that join the ground state of the reactants with the observed products in the IR-matrix spectra, mainly the oxide complex OM(η2-SO) and the sulfide oxide SMO2. Since the lowest-lying excited quadruplet state of niobium lies only 3.5 kcal/mol above the sextuplet ground state, the energy profile emerging from the fragments Nb(4F) + SO2 was also calculated [18]. All the located stationary points were characterized as energy minima or transition states through frequency analysis calculations at DFT level of theory. These calculations were done employing the functional B3LYP as used before by Liu et al. in a previous contribution on the investigated interactions for comparing the theoretical frequencies with the experimental ones [17]. Inclusion of a portion of exact exchange from Hartree–Fock theory using a hybrid functional can be important for the comparative description of systems involving electronic states with unpaired electrons, such as the quadruplet and doublets electronic states studied in this contribution. The energy of each of the stationary points was revaluated through single-point CASSCF-MRMP2 calculations.

For each of the investigated interactions, the energies of the lowest-lying quadruplet and doublet radical fragments MO + SO that vary only in the spin of the unpaired parallel electrons of the SO fragment were calculated at CASSCF-MRMP2 level of theory. Likewise, potential energy curves for the rebounding of the radical fragments to yield the metal oxide structure OM(η2-SO) were obtained from partial geometry calculations by approaching perpendicularly the metal fragment to the SO radical at different distances, without any additional symmetry restrictions.

For all the atoms, the Def2-TZVP basis sets optimized by Weigend and Ahlrichs were used [2022]. The niobium and tantalum atoms were described using the pseudopotentials proposed by Dolg et al., which include 28 and 60 electrons in the inner core, respectively [23].

The CASSCF calculations were carried out using an active space consisting of 7 electrons in 7 active orbitals, which expands around 700 CSF’s. The orbital space included in all the cases the d- and s-type functions of the transition metal atom as well as the outermost p-type functions of sulfur and oxygen. MRMP2 calculations were carried out using the same active space.

All the calculations were performed using the GAMESS and GAUSSIAN09 computational packages [24, 25].

3. Results and Discussion

3.1. V + SO2 Reaction

In Figure 1 are shown the two sequential radical reactions that lead from the ground state of the reactants V(4F) + SO2 to the products observed for this interaction under matrix-isolation conditions, mainly the oxide OV(η2-SO) and the sulfide oxide SVO2 compounds [17]. Geometrical parameters for these structures are provided in Table 1.

As it is seen in Figure 1, the radical fragments VO + SO lie 39 kcal/mol below the ground state reference. Hence, it is plausible that these radical species are yielded from the reactants V(4F) + SO2 in a first reaction. It has been proposed that a rhombic intermediate arising from the electrostatic interaction between the fragments could be formed in the first stage of the reaction [17].

However, the potential well calculated for this structure is only 1.7 kcal/mol below the energy corresponding to the radical species VO + SO. Hence, it could be expected that the reaction taking place from the ground state of the reactants passes by this adduct and reaches readily the radical species. Once the radical fragments are attained, they can recombine themselves in a second reaction to form the metal oxide structure OV(η2-SO). In Table 2 is shown the valence electronic configuration of the lowest-lying quadruplet and doublet states of the radical fragments that differ only in the spin of the unpaired parallel electrons located on the nonmetallic fragment. As these asymptotes are degenerated, the rebounding of the radical moieties can occur along both multiplicity channels (the low-multiplicity pathway is activated when the spin of the unpaired electrons in the SO moiety is opposite to that corresponding to the quadruplet electronic state). The potential energy curves for the recombination of the radical fragments are shown in Figure 2.

Stable structures for the oxide complex OV(η2-SO) are obtained from the radical rebounding along both spin channels. The doublet electronic state is more stable than the quadruplet by 17.5 kcal/mol. As it is seen in Figure 1, the reaction along the low-multiplicity channel can evolve to the sulfide oxide SVO2 after surmounting a relatively small energy barrier of 23.3 kcal/mol.

The potential well calculated for this product is 121.1 kcal/mol below the ground state reference and represents the global energy minimum for this interaction. For the quadruplet channel, the height of the barrier separating the OV(η2-SO) and the SVO2 species is 34.3 kcal/mol; it is unlikely that the sulfide oxide could be reached along this channel under cryogenic conditions. Thus, the reaction should stop once the quadruplet complex OV(η2-SO) is formed.

According to the proposed scheme, the products observed under confinement conditions for this reaction are obtained from different spin-multiplicity channels. The sulfide oxide SVO2 emerges from the rebounding reaction along the doublet channel, whereas the complex OV(η2-SO) is obtained when the reaction takes place through the quadruplet pathway.

3.2. Nb + SO2 and Ta + SO2 Reactions

In Figure 3, it is shown the two-step reaction scheme connecting the reactants with the sulfide oxide SNbO2. This is the only product observed for this reaction. The sextuplet sulfide oxide is predicted as unstable; it lies 69.9 kcal/mol above the ground state reference. For this reason, channels emerging from the sextuplet ground state of the reactants were not investigated. However, the lowest-lying quadruplet excited state of the reactants is located 5.7 kcal/mol above the ground state reference. It is plausible that the reaction taking place from this excited state plays an important role in determining the product distribution observed for this reaction by leading to the radical fragments NbO + SO, which lie 65.2 kcal/mol below the ground state reference.

As it is shown in Figure 3 and Table 2, the quadruplet and doublet radical asymptotes that vary only in the spin of the unpaired electrons located on the nonmetallic fragment (both parallel or antiparallel) are degenerated (−65.2 kcal/mol). Hence, recombination of the radical fragments can take place along both channels (Figure 2). The most stable structure for the complex ONb(η2-SO) is reached through the low-multiplicity channel. The energy minimum for this structure is located 120 kcal/mol below the ground state reference. The quadruplet state is roughly 14 kcal/mol above the doublet one. Unlike the vanadium interaction discussed before, reaction along both channels can reach the sulfide oxide SNbO2 after surmounting in each case a relatively small energy barrier (24.3 and 10.7 kcal/mol, resp.). The doublet pathway evolves by far to the most stable sulfide oxide structure. It lies 161.2 kcal/mol below the ground state of the original reactants Nb + SO2. Thus, the picture attained through the two-step reaction scheme is consistent with the experimental determinations made on this interaction under cryogenic conditions, as the sulfide oxide is the only detected product.

The two-reaction scheme shown in Figure 3 for the Ta + SO2 interaction is similar to that discussed before for the niobium reaction.

The degenerated quadruplet and doublet radical asymptotes TaO + SO are located around 30 kcal/mol below the ground state of the reactants (the small difference between the energies calculated for these asymptotes arises as consequence of the limited spaces used for the CASSCF expansions). The recombination of the radical fragments through both multiplicity channels leads to stable structures for the complex OTa(η2-SO). The energy potential well for the doublet state lies 144.7 kcal/mol below the ground state reference. Reaction proceeding along both the quadruplet and doublet channels could reach the sulfide oxide STaO2 (the height of the barriers separating the OTa(η2-SO) and the STaO2 species are 24.5 and 3.2 kcal/mol). As for the niobium reaction, this is the only detected product for this interaction under matrix-isolation conditions.

3.3. Spin-Flip Models versus the Two-Step Reaction Scheme

The picture attained for the investigated interactions in terms of two sequential reactions involving the radical fragments OM + SO is consistent with the product distribution determined for them in the IR-matrix spectra [17]. This scheme provides a plausible explanation for the low-spin products determined for these reactions without invoking unlikely interactions between electronic states of different spin multiplicity, as it has been proposed in the previous study by Liu et al. [17]. This could be of significance for the interactions of vanadium and niobium, as it is not expected that spin-orbit effects would be relevant in the description of their energy profiles. However, it is important stressing out that the corresponding reaction with tantalum is also suitably described by this model.

According to the proposed scheme, the recombination of the radical fragments through the low-multiplicity channel leads in all the cases to the products observed under matrix-isolation conditions. As the radical species formed in the first reaction would remain trapped, presumably the confinement imposed by the matrix could play an important role in determining the observed product distributions by enabling their recombination.

In a previous contribution, the results obtained from matrix-isolation determinations for the reactions of V, Nb, and Ta with CH3CN were described in terms of a similar reaction scheme involving the radical species CH3 + MNC [19]. As for the reactions with SO2 discussed here, the products determined for those reactions evolve from the doublet channel activated by the rebounding of radical fragments with opposite spin. It is striking that the same scheme explains for the main features determined for these groups of reactions.

Interestingly, the product distributions obtained under confinement conditions for other reactions have been analyzed by using a similar reaction scheme. The low-spin carbenes and carbynes detected in the IR-matrix spectra for the reactions of zirconium, iron, and ruthenium with fluoromethanes were rationalized by means of two sequential radical reactions [2628]. Likewise, in a recent contribution, we have extended the use of this scheme to the description of interactions involving transition-metal complexes. Particularly, we used this model to rationalize the irregularities observed in the product distribution determined for the reactions of triethyl-methane and silane molecules with the CpCOCo and CpCORh compounds [29]. Unlike previous studies on these interactions, in the above-mentioned investigations, we have not invoked interactions between the potential energy surfaces belonging to electronic states of different spin-multiplicity to explain for the low-spin species observed for these reactions [3035]. Hence, the use of this sequential radical reaction scheme overcomes the main disadvantage of the spin-flip models [3638]. The relativistic effects that could lead to interactions between electronic states of different spin multiplicity can be relevant in reactions involving heavy atoms, but not in those in which only lighter atoms participate.

In fact, important mechanistic aspects can be missed or hided when spin-flip models are applied to reactions involving only light atoms. As discussed before, the product distribution detected for the reaction of vanadium with SO2 is different from those observed for the corresponding reactions with niobium and tantalum. Whereas for the vanadium reaction are detected both the oxide OV(η2-SO) and the sulfide oxide SVO2 complexes, for the remaining reactions only the sulfide oxide is yielded [17]. According to the proposed scheme, the products observed for the vanadium interaction evolve from different spin channels. Reaction along the doublet channel reaches readily the sulfide oxide complex SVO2. However, due to the high energy barrier separating the OV(η2-SO) and the SVO2 compounds (34.2 kcal/mol), it is expected that the reaction through the quadruplet pathway ends when the oxide OV(η2-SO) is yielded. For the corresponding reactions with niobium and tantalum, the height of the energy barriers for the high-multiplicity channels is lower (around 24 kcal/mol in each case). For these interactions, both the quadruplet and doublet channels could attain the sulfide oxide. This interpretation is essentially different from that given by Liu et al. [17]. These authors explain the differences in the product distributions observed for these interactions in terms of the relative stability of the doublet OM(η2-SO) and the SMO2 complexes (the role played by the quadruplet state is missed once the intersystem crossing occurs).

4. Conclusions

The reactions of the SO2 molecule with the V, Nb, and Ta atoms have been studied theoretically through DFT and CASSCF-MRMP2 calculations. Energy profiles joining the ground state of the reactants with the products observed under cryogenic matrix-isolation conditions have been built for the three interactions by using a two-step reaction scheme involving the radical species MO + SO.

The description attained for each of the investigated interactions allows explaining for the low-spin products observed for them in the IR-matrix spectra without invoking interactions between electronic states of different spin multiplicity, as it has been proposed in a previous study. The use of this model also allows highlighting subtle aspects that can play an important role in determining the product distributions for these interactions.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the Universidad Nacional Autónoma de México (DGAPA-IN113515). The authors acknowledge the computational resources provided by DGCTIC-UNAM for the use of the Miztli supercomputer. C. Velásquez gratefully acknowledges CONACYT for his graduate scholarship (576688). A. E. Torres gratefully acknowledges CONACYT for the postdoctoral scholarship (364788/245467).