Printer Friendly

Prediction of tetraoxygen reaction mechanism with sulfur atom on the singlet potential energy surface.

1. Introduction

The oxygen atom exists naturally in four basic forms including free atomic particle, diatomic oxygen [O.sub.2], ozone [O.sub.3], and tetraoxygen [O.sub.4]. Tetraoxygen is typically an unstable, rare, nonmagnetic pale blue gas. Interest in tetraoxygen molecules dates to a 1924by Lewis [1]; since then, studies of [O.sub.4] known as "oxozone" have enjoyed along history of investigation [2-10].

Since its discovery, [O.sub.4] has been a subject of interest from both experimental and theoretical point of views. The first experimental study of molecular oxygen by using spectroscopic method came back to Arnold and coworker's research [11]. They suggested the existence of [O.sub.4] molecule as van der Waals complex. In addition, Helm and Walter [12], are observed formation of [O.sub.4] molecules when electron transferred to form of [O.sub.4.sup.+]. The kinetic energy release and the nature of dissociation have been investigated using electronic excited state. The metastable structure of [O.sub.4] molecule has been recently observed in mass spectroscopy experiments by Cacace et al. [13]. Then, Bevsek et al. [14] described 1:1 (or 1+1) resonant photoionization spectra of an energetic metastable [O.sub.4] species produced in a dc discharge. In another study, Adamantides et al. defined a long-predicted covalent forms, either cyclic, [D.sub.2d], or pinwheel, [D.sub.3h], structures [15].

In theoretical approach, Adamantides et al. (1980) [15] predicted the existence of metastable covalent [O.sub.4] molecule completely different from the van der Waals structure of [O.sub.2]-[O.sub.2] molecule which was experimentally detected. Then, Roeggn and Nilssen (1989) [16] have shown a branch form of [O.sub.4] analogues with the stable isoelectronic molecules N[O.sub.3.sup.-] and O[O.sub.3]. A central atom and three equivalent ligand atoms characterized it, which can exist in metastable [D.sub.3h] form. In 1990, Dunn et al. [6] predicted a cyclostructure for [O.sub.4] molecule with the bond length of 1.437 A. In addition, dissociation of cyclostructure [O.sub.4] molecule is investigated by Seidl and Schaefer [9] based on new high-energy materials. An ab initio calculation has been done on the ground and low-lying electronic states of [O.sub.4] along the minimum energy to explain the reactivity and electronic state of tetraoxygen. Obtained results provided a solid basis to establish the stability of the [O.sub.4] chemically bound molecule. In addition, surface crossings between singlet and triplet states were found and further characterized by evaluating their spin-orbit coupling matrix elements [7]. Peterka [5] present rotationally resolved photoionization spectra, photoelectron spectra, and ab initio calculations for providing strong evidence for the identity of this species as a novel complex between a ground state molecule and one excited state. In summary, theoretical calculations have predicted the existence of metastable [O.sub.4] molecules with two different shapes: cyclobutane and trigonal planar similar to boron trifluoride. The existence of [O.sub.4] was further confirmed in some other theoretical and experimental research [14,17-20].

On the other hand, processes such as volcanic eruptions, biogenic activity, and fossil fuel combustion all contribute to the emission of sulfur gases into the atmosphere. Sulfur has long been recognized as an important constituent of atmospheric aerosols and sulfur compounds are major pollutants of the environment [21]. The sixteen molecular species of type [S.sub.n][O.sub.m] (S[O.sub.2], S[O.sub.2], S[O.sub.3], S[O.sub.4], etc.) have been detected experimentally [22]. Among these species, the sulfur dioxide (S[O.sub.2]) molecule plays an important role in the atmospheric formation of sulfuric acid.

Sulfur tetraoxide (S[O.sub.4]) as an interesting atmospheric compound is potentially relevant to atmospheric chemistry [23]. McKee [24] considered several structures for the S[O.sub.4] species and predicted that the lowest-energy structure is to have a three-membered SOO ring and [C.sub.2v] symmetry. Schriver et al. [25] also reported IR bands at 1,442 and 1,272 [cm.sup.-1], which they attribute to S[O.sub.4]. Kugel and Taube [26] give a discussion of the synthesis and characterization of the monomeric S[O.sub.4]. Previously, sulfur tetraoxide has been proposed as intermediate in a variety of sulfur oxide reaction systems, as the active oxidizing species in per-sulfate oxidations [27]. Therefore, the experimental and theoretical studies reveal that the most stable structure for S[O.sub.4] is a three-membered ring (SOO) with two terminal oxygen atoms attached to ring sulfur with [C.sub.2v] symmetry (S[O.sub.4] ([C.sub.2v])). S[O.sub.4] ([C.sub.2v]) molecule has been confirmed in both the experimental and theoretical studies, while S[O.sub.4] ([C.sub.3v]) molecule has been only confirmed in the theoretical studies.

To our knowledge, there is not any theoretical or experimental work on the reaction mechanism of S+[O.sub.4] reaction. Hence, in this work the reaction mechanism of sulfur atom reaction with [O.sub.4] molecule on the singlet potential energy surface are investigated. Then, thermodynamic parameters of the reaction are obtained and results the reaction is spontaneous with the releasing of the energy. As well as, the structure of [O.sub.4] molecule and possible products are reported on the singlet potential energy surface.

2. Computational Details

All the calculations are performed with the Gaussian 03 system of codes [28]. The geometries of all species involved in the [sup.1]S+[sup.1][O.sub.4] reaction are fully optimized at the B3LYP [29] level in conjunction with the 6-311+G (3df) basis set. The 6-311+G (3df) basis set has been used because of the lack of hydrogen atoms in structures. These energies are corrected by zero point energies of B3LYP level. In present study the coupled cluster singles and doubles (CCSD) method has been employed [30, 31]. This method has the advantage of being able to evaluate the correlation surface in a size consistent manner. Hence, the CCSD calculations are carried out to compute accurate relative energies. The counterpoise procedure has been used to correct the interaction energy for basis set superposition error, BSSE, for the structures. So, the obtained energies are corrected by BSSE energies.

At the motioned level of treatment, we computed the vibrational frequencies to determine the nature of the stationary points according to the number of the eigenvalues of the Hessian matrix. The results show that one of reactants ([O.sub.4]), products, and intermediates possesses all real frequencies, and any transition state has only one imaginary frequency. The imaginary frequency that characterizes each transition state is included as well. Another reactant, S, has no vibrational frequencies. To verify that the TS structure is the right saddle point connecting the corresponding reactants and products of interest, the intrinsic reaction coordinate, IRC [32], calculations were performed at the B3LYP level. These calculations show that all TS structures of title reactions are valid. At the level of this method, it is possible to generate a wave function in a suitable form to execute a topological analysis of atoms in molecules [33] using the AIM2000 series programs [34]. Moreover, thermodynamic data have been calculated using the statistical mechanics.

3. Results and Discussion

The reaction mechanisms of S with [O.sub.4] have been investigated on singlet PES. One stable complex has been found between the reactants. Hence, all of the elementary reactions considered in the title reaction begin with the stable initial complex after passing corresponding transition states to produce final adducts. The intermediates prefix is IN followed by a number to differentiate one from each other. The transition states are designated by the prefix TS and a number. The structures of the reactants, products, intermediates, and transition states involved in the [sup.1]S+[sup.1][O.sub.4] reaction are shown in Figure 1. Their Cartesian coordinates are available in the Supplementary Material, available online at http://dx.doi.org/10.1155/2014/912391. By means of the TSs and their connected intermediates (or products), a schematic PES for the reaction is plotted in Figure 2. Table 1 gives the total energies for the B3LYP and CCSD levels and corresponding relative energies in comparison with the reactants in the singlet state. To simplify our discussion, the energy of reactant R, [sup.1]S+[sup.1][O.sub.4], is set to be zero for reference. The AIM theory topological parameters for [O.sub.4] as a reactant on the B3LYP level are exhibited in Figure 3. In addition, the topological characteristics in the bond critical points in the mentioned molecules are tabulated in Table 2. The vibrational frequencies of all species are summarized in Table 3. In addition, the values of the AIM theory topological parameters for selected bonds of products on the B3LYP level are shown in Figure 3. Then, the topological characteristics in the bond critical points, BCPs, in the molecules are tabulated in Table 4. The highest occupied molecular orbital, HOMO, and the lowest unoccupied molecular orbital, LUMO, energies and band gap energy for HOMO to LUMO transition of the selected species for products are shown in Table 5. In addition, frontier molecular orbitals of the same molecules are shown in Figure 4. Finally, the standard thermodynamic functions at room temperature and atmospheric pressure are tabulated in Table 6.

3.1. [O.sub.4] Molecular Property. Up to date, the experimental and theoretical studies show that [O.sub.4] has different structures with various group points. In addition, its bond strengths and characteristic (van der Waals or covalentbonds) are discussed [7, 12, 16]. In our study, we applied the new structure of neutral tetraoxygen with group point of [D.sub.2h]. Of course, in study of Helm and Walter, [D.sub.2h] group point of [O.sub.4] with positive charge has been investigated. [O.sub.4] ([D.sub.2d]) has stability energy of -299.0706 au while its structure with the [D.sub.3h] group point shows stability energy of -300.0949 au. In addition, in the mentioned study of Helm and Walter, the bond length of O-O in the [O.sub.4] ([D.sub.2d]) is reported 1.433 [Angstrom] and angle of O1O2O3 atoms is 87.0[degrees] while in the [O.sub.4] with [D.sub.3h] point group, the O-O bond length is obtained 1.330 [Angstrom] and angle of three mentioned atoms is 120.0[degrees]. In [O.sub.4]+ ([D.sub.2h]), the bond lengths are 1.160 and 2.396 [Angstrom], which shows rectangular shape of molecule [12].

In the [O.sub.4] ([D.sub.2h]) on the singlet PES, our results confirmed that the bond lengths are 1.934 and 1.201 [Angstrom] with rectangular structure of molecule. The stability energy of that is -300.6897 au, which is stable than the other structures with diverse group points on the singlet PES.

The bond formation between the oxygen atoms in the tetraoxygen molecule on the singlet and triplet states is confirmed by atoms in molecules topological analysis of wave function. This analysis reveals the presence of bond critical points located between [O.sub.1]-[O.sub.2] atoms in [sup.1][O.sub.4] ([D.sub.2h]) with density of the wave function: [rho]([r.sub.bcp]) = 0.0817 au and the Laplacian of electronic density at critical point, [[nabla].sup.2][rho]([r.sub.bcp]) = 0.3479 au. The results confirmed that van der Waals interaction between the [O.sub.1]-[O.sub.2] atoms. Although, the AIM analysis reveals covalent bond formation between the [O.sub.2]-[O.sub.4] with the Laplacian of electronic charge density of [[nabla].sup.2][rho]([r.sub.bcp]) = -0.7198 au. As shown in Figure 3, furthermore, we are analyzed the structure of tetraoxygen on the triplet PES. AIM analyses of [sup.3][O.sub.4] ([D.sub.2h]) confirm the presence of one BCP located between [O.sub.1]-[O.sub.2] with [rho]([r.sub.bcp]) = 0.0527 au, and [[nabla].sup.2][rho]([r.sub.bcp]) = 0.2142 au. The results display the van der Waals interaction between [O.sub.1]-[O.sub.2] atoms. In addition, another BCP is located between [O.sub.2]-[O.sub.4] atoms in [sup.3][O.sub.4] ([D.sub.2h]) molecule with the electronic density [rho]([r.sub.bcp]) = 0.5507 au, and the Laplacian of electronic density [[nabla].sup.2][rho]([r.sub.bcp]) = -0.7601 au. The AIM analysis shows the covalent bond formation between the [O.sub.2]-[O.sub.4] atoms. So, this work suggestion is two van der Waals interaction and two covalent bonds in the structure of [O.sub.4] ([D.sub.2h]) molecule on the singlet and triplet states. The relative energies and Laplacian of electronic charge density analysis confirmed that the tetraoxygen molecule on the triplet state with [D.sub.2h] point group is more realable than the other suggested structures.

3.2. Initial Interactions of Reactants. Since the ground state of sulfur atom is triplet, the mechanism of S + [O.sub.4] reaction is considered on the triplet PES. In spite of numerous attempts, we could not find any transition state for the reaction on the triplet PES. Only two complexes between the reactants are found on the triplet PES. Therefore, our attempt focused on the singlet PES for finding some intermediate, transition state, products, and suggested mechanism for the reaction pathway at the B3LYP level. Therefore, in the first step, the sulfur atom in triplet state, [sup.3]S, excited to the singlet state, and the reaction can continue in the singlet state This is explained as follows.

One stable reactant complex, IN1 ([O.sub.2]-S-[O.sub.2]), has been considered between the reactants on the singlet PES. The intermediate IN1 is formed when the sulfur atom bridged to the oxygen atoms of [O.sub.4] using the formation of four van der Waals bonds and rupture two [O.sub.2]-[O.sub.2] van der Waals interactions in [O.sub.4] molecule. In the intermediate IN1, the newly formed S-O bonds are two types. One, bond lengths are 1.525 [Angstrom] and two they are 1.702 [Angstrom] at the B3LYP level with the formation relative energy of -109.524 kcal/mol. No transition states have been found for the formation of IN1 intermediate. Therefore, the results show that the formation of IN1 is a barrier-less and exothermic process. In the complex, four oxygen atoms of [O.sub.4] are associated with sulfur atom and two three-membered rings are formed (see IN1 structure in Figure 1). In the first ring formation the lone pair of the [O.sub.2] and [O.sub.4] atoms of tetraoxygen interacts with the sulfur atom and in the second ring [O.sub.1] and [O.sub.3] atoms associate with the atomic sulfur. Molecular geometry indicates that the IN1 complex has C1 symmetry.

The stable reactant complex formation is confirmed by AIM topological analysis of wave function. Our analysis reveals the presence of bond critical points; one located between two oxygen atoms of [O.sub.4] molecule and sulfur atom with density of wave function, ([rho]([r.sub.bcp]) = 0.1698 au and [[nabla].sup.2][rho]([r.sub.bcp]) = 0.0245 au) and the other located between two other oxygen and sulfur atoms ([rho]([r.sub.bcp]) = 0.2267 au and [[nabla].sup.2][rho]([r.sub.bcp]) = 0.0196 au). Therefore, all the O ... S interactions are of the van der Waals types. In addition, the topological analysis of the wave function shows the existence of two bond critical points; one is located between [O.sub.1] and [O.sub.3] ([rho]([r.sub.bcp]) = 0.2300au and [[nabla].sup.2][rho]([r.sub.bcp]) = 0.2334 au), and the other is located between [O.sub.2] and [O.sub.4] with the same electronic density and Laplacian. Calculated results allow us to classify the O ... O interactions as being of the van der Waals type. The rings formation between [O.sub.2][S.sub.5][O.sub.4] and [O.sub.1][S.sub.5][O.sub.3] atoms also is lightened using AIM analysis. The topological analysis of the wave functions and their Laplacians are [rho]([r.sub.bcp]) = 0.1449 au and [[nabla].sup.2][rho]([r.sub.bcp]) = 0.4256 au, respectively.

3.3. Formation Pathways of P1 (OOSOO). There is one possible pathway for production of P1 as follows:

Path 1: R [right arrow] IN1 [right arrow] TS1 [right arrow] IN2 [right arrow] TS2 [right arrow] P1 (OOSOO).

As is shown in path 1, the singlet reactants of S and [O.sub.4] directly transformed to IN1 without entrancing any transition state. Then IN1 transformed to IN2 after passing TS1 with the energy barrier of 35.429 kcal/mol at the B3LYP level. In this step, one O-S van der Waals bond of three-membered ring structure in IN1 molecule ruptures and IN2 appears. IN2 is 86.279 kcal/mol lower than the[sup. 1]S+[sup.1][O.sub.4] original reactant. TS1 has one imaginary frequency at 378i [cm.sup.-1] that is vibrated in the reaction coordinate. TS1 has a three-membered ring structure (see Figure 1). The electronic charge density and its Laplacian in ring critical points are 0.1712 and 0.4378 au, respectively. The intermediate IN2 converted into final product after passing TS2 with the energy barrier of 35.155 kcal/mol and imaginary frequency of 163i [cm.sup.-1] in the reaction coordinate. The S-O bond of three-membered ring complex, IN2, is ruptured and converts to the final product P1. The conversion is an endothermic process but the overall reaction from singlet reactants is exothermic with the reaction enthalpy -68.286 kcal/mol and is spontaneous with the Gibbs free energy -62.177 kcal/mol. All steps of reaction as formation of ring structure and its collapse are confirmed by electronic charge density analysis using AIM calculation results. Formation channel of product P1 with two transition states and lower level of energy barriers is favorable process from kinetic point of view.

3.4. Formation Pathways of P2 ([O.sub.2] ([sup.1][DELTA] g) + S[O.sub.2]). There is one possible pathway for production of P2 as follows:

Path 2: R [right arrow] IN1 [right arrow] TS1 [right arrow] IN2 [right arrow] TS3 [right arrow] IN3 [right arrow] TS4 [right arrow] IN4 [right arrow] P2 ([sup.1][O.sub.2] ([sup.1][[DELTA].sub.g]) + S[O.sub.2]).

As shown in Figure 2, formation of the intermediate IN2 is similar to path P1; then, it transforms into IN3 using TS3 with energy barrier of 11.048 kcal/mol and imaginary frequency of 591i [cm.sup.-1]. The [O.sub.2]-[O.sub.4] bond in ring structure of TS3 is 51% longer in comparison with the corresponding bond in IN2. Four-membered ring complex, IN3, (see Figure 1) undergoes elongation of [O.sub.3]-S bond and collapse of four-membered ring structure turn to IN4 after entrancing TS4 with the energy barrier of 11.892 kcal/mol and imaginary frequency of 152i [cm.sup.-1]. The elongation of the O3-S bond of TS4 almost 0.379 corresponds to a relative increase of 22% in TS4 compared to the IN3 molecule. Obtained specie is known as a complex between singlet molecular oxygen and sulfur dioxide. Then the product complex IN4 can directly transform into the product [sup.1][O.sub.2]+S[O.sub.2] through O-O bond rupture without any transition state. Singlet molecular oxygen due to spin multiplicity changes and transforms to triplet oxygen through a relaxation process. The reaction thermodynamic parameters indicate that P2 complex, [sup.1][O.sub.2] ([sup.11][[DELTA].sub.g]) + S[O.sub.2], formation is exothermic with enthalpy of reaction, -170.708 kcal/mol, in comparison with the original reactants. In all steps of path P2, the formation of products from reactants is confirmed using related IRC curves.

3.5. Formation Pathways ofP3 (S[O.sub.4] ([C.sub.2v])). For production of P3, there is one possible pathway as follows:

Path 3: R [right arrow] IN1 [right arrow] TS1 [right arrow] IN2 [right arrow] TS3 [right arrow] IN3 [right arrow] TS4 [right arrow] IN4 [right arrow] TS5 [right arrow] P3.

In this path, the formation of the intermediate IN4 is similar to path P2. Then, the IN4 exchange to final product through TS5 by formation of S-[O.sub.4] and [O.sub.3]-S bonds. For this step, barrier height is 33.603 kcal/mol. The optimized structure of transition state TS5 shows that [O.sub.3]-S bond length is 45% longer than corresponding bond in S[O.sub.4]. On the other hand, newly formed bond, [O.sub.4]-S, is 51% greater than the corresponding bond in IN4 molecule. Imaginary frequency of TS5 is 455i [cm.sup.-1].

The reaction thermodynamic parameters indicate that P3 formation process is exothermic by 188.042 kcal/mol releasing heat and spontaneous with -179.147 kcal/mol in Gibbs free energy at atmospheric pressure and 298.15 K temperature. The formation of products along the path is confirmed using AIM results and IRC curve.

3.6. Formation Pathways of P4 (S[O.sub.4] ([C.sub.3v])). For production of P4, there is one possible pathway as follows:

Path 4: R [right arrow] IN1 [right arrow] TS1 [right arrow] IN2 [right arrow] TS3 [right arrow] IN3 [right arrow] TS4 [right arrow] IN4 [right arrow] TS5 [right arrow] P3 [right arrow] TS6 [right arrow] P4.

In this path, the formation of the P3 is similar to path P3. Then P3 can be converted into the P4 by angle-variation through TS6 transition state. The angle of [O.sub.3]S[O.sub.4] in P3 is 59.6[degrees] while it is 106.0[degrees] in P4. The breaking bond of [O.sub.3]-[O.sub.4] is 1.576 in P3 that extend to 2.176 in TS6. The energy height of this process is 18.231 kcal/mol and imaginary frequency of TS6 is 281i [cm.sup.-1].

The reaction thermodynamic parameters show that P4 formation process is exothermic by 164.857 kcal/mol releasing heat and spontaneous with -156.922 kcal/mol in Gibbs free energy. The S[O.sub.4] ([C.sub.2v]) molecule has been confirmed in both the experimental and theoretical studies, while S[O.sub.4] ([C.sub.3v]) molecule has been only confirmed in the theoretical studies. In this work, we are obtained both of them that molecule with [C.sub.2v] symmetry is more stable than another. In summary, the [sup.1]S+[sup.1][O.sub.4] reaction is energetically feasible to formation of all four products all over corresponding paths which all step of reaction confirm using IRC and AIM calculations.

3.7. Topological Analysis of Electronic Density of the Products. AIM topological analysis has been shown to provide important information about many different chemical systems by an analysis of molecular electron density distribution. The electron density, [rho](r), density Laplacian, [[nabla].sup.2][rho](r), and bond ellipticity, [epsilon] = ([[lambda].sub.1]/ [[lambda].sub.2]) - 1, at bond critical points (BCP) have been extracted using AIM topological analysis. The A1, A2, and [[lambda].sub.3] are the eigenvalues of the Hessian matrix. The first two eigenvalues correspond to the perpendicular curvatures and the latter provides curvatures along the internuclear axis. According to the AIM topological analysis, electron density, [rho](r), and Laplacian of the electron density, [[nabla].sup.2][rho](r), are used to describe the strength and the characteristic of the bond, respectively. According to the theory of AIM, the positive value of the Laplacian ([[nabla].sup.2][rho](r)) indicates a weak interaction or an ionic bond; however, the negative value of the Laplacian shows a strong covalent bond between the atoms. The ellipticity shows that the ratio of the rate of density decrease in two directions perpendicular to the bond path at the bond critical point.

The topological analysis of the electronic density were performed for all bonds in the OOSOO, S[O.sub.2] and S[O.sub.4] ([C.sub.2v]), and S[O.sub.4] ([C.sub.3v]) species as stable products of our studied reaction. The topological characteristics in the BCP's of S-O and O-O bonds are tabulated in Table 4. In addition, the [[nabla].sup.2][rho](r) values in Table 4 indicate that the every bond in OOSOO has covalent character ([[nabla].sup.2][rho](r) < 0). Their electron density values show that the bond strength manner O-O > S-O is corresponding to the electronic charge densities 0.3606 and 0.2283 au, respectively.

In S[O.sub.4] ([C.sub.2v]) molecule, the S-O bond has covalent character with Laplacian of -0.0738 au, whereas the S=O bond has electrostatic interaction nature with Laplacian of the electron density value 0.9675 and O-O bond with 0.3346 au in [[nabla].sup.2][rho](r) has van der Waals origin. The S=O bond in S[O.sub.4] ([C.sub.2v]), in comparison with corresponding bonds in S[O.sub.2] and S[O.sub.4] ([C.sub.3v]) molecules, has high electronic density. The results have confirmed that mentioned bond in S[O.sub.4] ([C.sub.2v]) is stronger than the corresponding bonds in other species. The Laplacian of S=O bond in S[O.sub.2] and S[O.sub.4] ([C.sub.2v]) is very similar and its absolute value is greater than Laplacian of corresponding bond in S[O.sub.4] ([C.sub.3v]). The results suggested that S=O bonds in three mentioned molecules have electrostatic nature and can be ionic origin.

The values for the ellipticity of various bonds for three products are listed in Table 4. The high ellipticity ([epsilon]) for S-O bond in OOSOO molecule shows nonsymmetrical distribution of density about the bond path, whereas low value belongs to the O-O bond which is corresponding with symmetrical distribution of density between two atoms. In S[O.sub.4] ([C.sub.3v]) molecule, all of the bonds have the same ellipticity (0.0462 au). This manner shows that bonds have the same symmetry in distribution of [pi] bonding character. Comparison of the S-O bond ellipticities in all of species shows that the S-O electron density in OOSOO builds up in a particular orientation and it reveals nonsymmetrical distribution of density between two atoms.

3.8. Molecular Orbital Analyses of the Products. The features of the HOMO and LUMO for products are calculated by B3LYP/6-311+G (3df) level in gaseous phase. The HOmO energy characterizes the ability of electron to donate, and the LUMO energy characterizes the ability of electron to accept [36]. The energy gap between HOMO and LUMO characterizes whether the molecules are chemically stable or not. The energy band gap ([DELTA][epsilon]) (transition from HOMO to LUMO), ionization potential, and electron affinity of the OOSOO, S[O.sub.2], S[O.sub.4] ([C.sub.2v]), and S[O.sub.4] ([C.sub.3v]) species can be seen in Table 5. The HOMO-LUMO energy gap results predict that the S[O.sub.2] (0.2073 au) and S[O.sub.4] ([C.sub.2v]) (0.2027 au) with large energy gaps ([DELTA][epsilon]) are more stable than the OOSOO (0.0.610 au) and S[O.sub.4] ([C.sub.3v]) (0.0539 au) molecules. The frontier molecular orbitals of four products are depicted in Figure 4. In OOSOO molecule, the LUMO is localized on the bond between oxygen-sulfur atoms. The HOMO is mainly composed on the p-type atomic orbitals of oxygen and sulfur atoms. In S[O.sub.2] molecule, the HOMO orbital is localized on the bond between oxygen-sulfur atoms while the LUMO is composed in the [pi] type orbitals of tree atoms. In addition, results show that HOMO orbital is localized on the [pi] type orbitals of oxygen atoms in S[O.sub.4] ([C.sub.2v]) while the LUMO is composed mainly just in the n type orbitals of two oxygen atoms. Finally, in S[O.sub.4] ([C.sub.3v]), the HoMO and LUMO orbitals that are completely similar are localized in the p-type atomic orbitals of oxygen atoms.

4. Conclusion

In the percent work, the reaction mechanism of sulfur atom with [O.sub.4] molecule is investigated on the singlet PES. The experimental and theoretical studies show that [O.sub.4] has different structures with various points group. In the neutral structure of tetraoxygen with [D.sub.2h] point group and triplet state are more stable than singlet state and other structures. In spite of numerous attempts, no transition state has been found for the S+[O.sub.4] reaction on the triplet PES. Therefore, the details of theoretical investigation on the singlet PES of S+[O.sub.4] reaction have been carried out at the B3LYP/6-311+G (3df) and CCSD levels. Therefore, in the first step, the sulfur atom on the triplet state, [sup.3]S+[sup.1][O.sub.4], is excited to the singlet state, [sup.1]S+[sup.1][O.sub.4], and the reaction continues in the singlet state. One stable reactant complex, INI ([O.sub.2]-S-[O.sub.2]), has been considered between the reactants on the singlet PES. Through variety of INI transformations, four kinds of products OOSOO, S[O.sub.2] + [O.sub.2] ([sup.1][DELTA] g), S[O.sub.4] ([C.sub.2v]), and S[O.sub.4] ([C.sub.3v]) are obtained which have enough thermodynamic stability. In conclusion, S[O.sub.4] ([C.sub.2v]) is the most stable product and its production is spontaneous and exothermic reaction with -188.042 and -179.147 kcal/mol in Gibbs free energy and enthalpy of reaction at the B3LYP level, respectively. In kinetic viewpoint, it is expected that the formation channel of P1 with two transition states and low level energy barriers is favorable process.

http://dx.doi.org/10.1155/2014/912391

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] G. N. Lewis, "The magnetism of oxygen and the molecule [O.sub.4]," Journal of the American Chemical Society, vol. 46, no. 9, pp. 2027-2032, 1924.

[2] A. J. C. Varandas and J. L. Llanio-Trujillo, "On triplet tetraoxygen: ab initio study along minimum energy path and global modelling," Chemical Physics Letters, vol. 356, no. 5-6, pp. 585-594, 2002.

[3] O. Prasad, L. Sinha, N. Misra et al., "Study of electrostatic potential surface and molecular orbitals of [O.sub.4] nano cluster by first principles," Scholars Research Library, vol. 1, no. 2, pp. 79-85, 2009.

[4] C. A. Long and G. E. Ewing, "Spectroscopic investigation of van der waals molecules. I. The infrared and visible spectra of [([O.sub.2]).sub.2]," The Journal of Chemical Physics, vol. 58, pp. 4824-4834, 1973.

[5] D. S. Peterka, "Unraveling the mysteries of metastable [O.sub.4]," Journal of Chemical Physics, vol. 110, no. 13, pp. 6095-6098, 1999.

[6] K. M. Dunn, G. E. Scuseria, and H. F. Schaefer III, "The infrared spectrum of cyclotetraoxygen, [O.sub.4]: a theoretical investigation employing the single and double excitation coupled cluster method," The Journal of Chemical Physics, vol. 92, no. 10, pp. 6077-6080, 1990.

[7] R. Hernandez-Lamoneda and A. Ramirez-Solis, "Reactivity and electronic states of [O.sub.4] along minimum energy paths," Journal of Chemical Physics, vol. 113, no. 10, pp. 4139-4145, 2000.

[8] D. Pillay, Y. Wang, and G. S. Hwang, "Prediction of tetraoxygen formation on rutile Ti[O.sub.2] (110)," Journal of the American Chemical Society, vol. 128, no. 43, pp. 14000-14001, 2006.

[9] E. T Seidl and H. F. Schaefer III, "Is there a transition state for the unimolecular dissociation of cyclotetraoxygen ([O.sub.4])?" The Journal of Chemical Physics, vol. 96, no. 2, pp. 1176-1182, 1992.

[10] V. Aquilanti, D. Ascenzi, M. Bartolomei et al., "Molecular beam scattering of aligned oxygen molecules. The nature of the bond in the [O.sub.2]--[O.sub.2] dimer," Journal of the American Chemical Society, vol. 121, no. 46, pp. 10794-10802, 1999.

[11] S. J. Arnold, E. A. Ogryzlo, and H. Witzke, "Some new emission bands of molecular oxygen," The Journal of Chemical Physics, vol. 40, no. 6, pp. 1769-1770, 1964.

[12] H. Helm and C. W. Walter, "Observation of electronically excited states of tetraoxygen," The Journal of Chemical Physics, vol. 98, no. 7, pp. 5444-5449, 1993.

[13] F. Cacace, G. de Petris, and A. Troiani, "Experimental detection of tetraoxygen," Angewandte Chemie International Edition, vol. 40, pp. 4062-4065, 2001.

[14] H. M. Bevsek, M. Ahmed, D. S. Peterka, F. C. Sailes, and A. G. Suits, "Direct detection and spectroscopy of [O.sub.4]," Faraday Discussions, vol. 108, pp. 131-138, 1997.

[15] V. Adamantides, D. Neisius, and G. Verhaegen, "Ab initio study of the [O.sub.4] molecule," Chemical Physics, vol. 48, no. 2, pp. 215-220, 1980.

[16] I. Roeggen and E. W. Nilssen, "Prediction of a metastable [D.sub.3h] form of tetra oxygen," Chemical Physics Letters, vol. 157, no. 5, pp. 409-414, 1989.

[17] F. A. Gorelli, L. Ulivi, M. Santoro, and R. Bini, "The [epsilon] phase of solid oxygen: evidence of an [O.sub.4] molecule lattice," Physical Review Letters, vol. 83, no. 20, pp. 4093-4096, 1999.

[18] L. F. Lundegaard, G. Weck, M. I. McMahon, S. Desgreniers, and P. Loubeyre, "Observation of an O8 molecular lattice in the e phase of solid oxygen," Nature, vol. 443, no. 7108, pp. 201-204, 2006.

[19] M. Hotokka and P Pyykko, "An ab initio study of bonding trends in the series B[O.sup.3-.sub.3], C[O.sup.2-.sub.33], N[O.sup.-.sub.3] and [O.sub.4]([D.sub.3h])," Chemical Physics Letters, vol. 157, no. 5, pp. 415-418, 1989.

[20] A. H. Jubert and E. L. Varetti, "On the possible existence of the [O.sub.4] molecule with D3h symmetry," Anales de Quimica, vol. 82, pp. 227-230, 1986 (Spanish).

[21] H. R. Pruppacher and J. D. Klett, Microphysics of Clouds and Precipitation, Reidel, Boston, Mass, USA, 1980.

[22] M. W. Wong, Y. Steudel, and R. Steudel, "Structures and vibrational spectra of the sulfur-rich oxides [S.sub.n]O (n = 4-9): the importance of [[pi].sup.*]- [[pi].sup.*] interactions," Chemistry A, vol. 13, no. 2, pp. 502-514, 2007.

[23] A. Jacob and C. A. Winkler, "Kinetics of the reactions of oxygen atoms and nitrogen atoms with sulphur trioxide," Journal of the Chemical Society, Faraday Transactions, vol. 68, pp. 2077-2082, 1972.

[24] M. L. McKee, "Computational studies on S[O.sub.4] and [S.sub.2][O.sub.3]," Journal of the American Chemical Society, vol. 115, no. 20, pp. 9136-9142, 1993.

[25] L. Schriver, D. Carrere, A. Schriver, and K. Jaeger, "Matrix-isolation photolysis of S[O.sub.2], [O.sub.3] and [H.sub.2]O: evidence for the [H.sub.2]O:S[O.sub.3] complex," Chemical Physics Letters, vol. 181, no. 6, pp. 505-511, 1991.

[26] R. Kugel and H. Taube, "Infrared spectrum and structure of matrix-isolated sulfur tetroxide," Journal of Physical Chemistry, vol. 79, no. 20, pp. 2130-2135, 1975.

[27] L. S. Levitt, "A new mechanism for persulphate oxidations," Canadian Journal of Chemistry, vol. 31, no. 10, pp. 915-922, 1953.

[28] M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Pople JA Gaussian 03, Revision B.03, Gaussian, Pittsburgh, Pa, USA, 2003.

[29] A. D. Becke, "A new mixing of hartree-fock and local density-functional theories," The Journal of Chemical Physics, vol. 98, no. 2, pp. 1372-1377, 1993.

[30] G. E. Scuseria and H. F. Schaefer III, "Is Coupled Cluster Singles and Doubles (CCSD) more computationally intensive than Quadratic Configuration Interaction (QCISD)?" The Journal of Chemical Physics, vol. 90, no. 7, pp. 3700-3703, 1989.

[31] G. E. Scuseria, C. L. Janssen, and H. F. Schaefer III, "An efficient reformulation of the closed-shell Coupled Cluster Single and Double Excitation (CCSD) equations," The Journal of Chemical Physics, vol. 89, no. 12, pp. 7382-7387, 1988.

[32] C. Gonzalez and H. B. Schlegel, "Reaction path following in mass-weighted internal coordinates," Journal of Physical Chemistry, vol. 94, no. 14, pp. 5523-5527, 1990.

[33] R. G. Parr and W. Yang, Density-Functional Theory of Atoms and Molecules, Oxford University Press, Oxford, UK, 1989.

[34] F. Biegler-Konig and J. Schoenbohm, AIM2000, Buro fur Innovative Software, Bielefeld, Germany, 2nd edition, 2002.

[35] R. D. Johnson III and NIST, Eds., "Computational chemistry comparison and benchmark database," NIST Standard Reference Database, Number 69, April, 2013, http://cccbdb.nist.gov/.

[36] P. W. Atkins, Physical Chemistry, Oxford University Press, Oxford, UK, 2001.

Ashraf Khademzadeh, Morteza Vahedpour, and Fereshte Karami

Chemistry Department, University of Zanjan, Zanjan 45371-38791, Iran

Correspondence should be addressed to Morteza Vahedpour; vahed@znu.ac.ir

Received 17 August 2013; Accepted 11 October 2013; Published 23 January 2014

Academic Editors: D. Dondi and B. Fernandez

TABLE 1: The total energies (in Hartree) and relative energies (in
parentheses and kcal/mol) of the reactants, products, and
intermediates in the [sup.1]S+[sup.1][O.sub.4] reaction on the
singlet PES.

Species                           B3LYP                  CCSD

[sup.1][O.sub.4]                -300.6780             -300.1169
[sup.1]S+[sup.1][O.sub.4]   -698.7512 (0.000)     -697.7038 (0.000)
[sup.3]S+[sup.1][O.sub.4]  -698.8125 (-38.466)   -697.7576 (-33.760)
IN1                        -698.9258 (-109.524)  -697.8998 (-122.962)
IN2                        -698.8887 (-86.279)   -697.8493 (-91.305)
IN3                        -698.9700 (-137.288)  -697.9428 (-149.967)
IN4                        -699.0211 (-169.363)  -698.0045 (-188.706)
TS1                        -698.8693 (-74.095)   -697.8344 (-81.980)
TS2                        -698.8327 (-51.124)   -697.7908 (-54.616)
TS3                        -698.8711 (-75.231)   -697.8172 (-71.159)
TS4                        -698.9510 (-125.396)  -697.9046 (-125.983)
TS5                        -698.9676 (-135.760)  -697.9385 (-147.258)
TS6                        -699.0077 (-160.916)  -697.9687 (-166.210)
OOSOO                      -698.8598 (-68.119)   -697.7884 (-53.090)
S[O.sub.2]+[O.sub.2]
  ([sup.1][[DELTA]
  .sub.g])                 -699.0234 (-170.792)  -698.0023 (-187.280)
S[O.sub.4] ([C.sub.2v])    -699.0489 (-186.769)  -698.0316 (-205.672)
S[O.sub.4] ([C.sub.3v])    -699.0127 (-164.089)  -697.9665 (-164.835)

TABLE 2: Topological analyses of [O.sub.4] molecules on the B3LYP
level; values are in atomic units.

                     Bond
                   critical
Species             point     [rho](r)   [[lambda].sub.1]

[sup.1][O.sub.4]   Point 1     0.0817        -0.1465
                   Point 2     0.5523        -1.4969

[sup.3][O.sub.4]   Point 1     0.0527        -0.0813
                   Point 2     0.5507        -1.4961

Species            [[lambda].sub.2]   [[lambda].sub.3]

[sup.1][O.sub.4]       -0.1199             0.6143
                       -1.4626             2.2397

[sup.3][O.sub.4]       -0.0736             0.3691
                       -1.4713             2.2321

Species            [[nabla].sup.2][rho](r)    Nature of bond

[sup.1][O.sub.4]            0.3479            van der Waals
                           -0.7198               Covalent

[sup.3][O.sub.4]            0.2142            van der Waals
                           -0.7601               Covalent

TABLE 3: The vibrational frequencies ([cm.sup.-1]) of the reactants,
products, and transition states calculated at the B3LYP level.
The value in parentheses refers to experimental data [35].

Species                              Frequencies (B3LYP)

[sup.1][O.sub.4]                     198, 519, 547, 810,1432,1624
[sup.3][O.sub.4]                     248, 279, 398, 626, 1471, 1637
IN1                                  246, 332, 414, 521, 567, 759, 776,
                                     984, 1041
IN2                                  156, 184, 370, 546, 554, 666,791,
                                     842, 1010
IN3                                  73, 254, 411, 602, 652, 775, 921,
                                     930, 1301
IN4                                  17, 47, 52, 57, 121, 521, 1177,
                                     1371, 1633
TS1                                  -378, 202, 242, 410, 546, 608,
                                     696, 929,1042
TS2                                  -163, 93, 174, 251, 369, 696,
                                     759, 921, 1012
TS3                                  -591, 93, 369, 399, 520, 608,
                                     767, 970, 1052
TS4                                  -152, 265, 355, 420, 576, 795,
                                     874, 992, 1233
TS5                                  -455, 188, 290, 424, 488, 690,
                                     906, 1250, 1442
TS6                                  -281, 304, 334, 384, 479, 979,
                                     986, 1110, 1130
OOSOO                                74, 80, 235, 317, 343, 694, 706,
                                     925, 1005
S[O.sub.2]                           519 (518), 1178 (1151),
                                     1377 (1362)
[O.sub.2] ([sup.1][[DELTA].sub.g])   1634 (1484)
S[O.sub.4] ([C.sub.2v])              310, 478 (490), 485 (498), 497,
                                     678 (611), 737 (777), 967 (927),
                                     1279 (1267), 1447 (1437)
S[O.sub.4] ([C.sub.3v])              243, 244, 298, 308, 477, 994,
                                     1035, 1088, 1089

TABLE 4: Topological analyses of products in atomic units.

                                    OOSOO              S[O.sub.2]

                               O-O       S-O        S-O         S=O
[[lambda].sub.1]             -0.8486   -0.4070    -0.6145     -0.6264
[[lambda].sub.2]             -0.8424   -0.2703    -0.5457     -0.5500
[[lambda].sub.3]             1.6457    0.5719      2.1282     2.0939
[rho](r)                     0.3606    0.2283      0.3080     0.3283
[[nabla].sup.2][rho](r)      -0.0453   -0.1054     0.9680     0.9175
[epsilon]                    0.0073    0.5057      0.1261     0.1389

                              S[O.sub.4]                S[O.sub.4]
                             ([C.sub.2v])      O-O     ([C.sub.3v])

                                  S-O          O-O         S=O
[[lambda].sub.1]                -0.4161      -0.3970     -0.6069
[[lambda].sub.2]                -0.3708      -0.2838     -0.5801
[[lambda].sub.3]                0.7131       1.0154       1.6849
[rho](r)                        0.2373       0.2038       0.3015
[[nabla].sup.2][rho](r)         -0.0738      0.3346       0.4979
[epsilon]                       0.1222       0.3989       0.0462

TABLE 5: HOMO and LUMO energies (in atomic unit) of OOSOO,
S[O.sub.2], S[O.sub.4] ([C.sub.2v]), and S[O.sub.4]
([C.sub.3v]) molecules.

                                                  S[O.sub.4]
                           OOSOO    S[O.sub.2]   ([C.sub.2v])

[E.sub.HOMO]              -0.2763    -0.3472       -0.3640
[E.sub.LUMO]              -0.2153    -0.1399       -0.1613
[[DELTA].sub.[epsilon]]   -0.0610    -0.2073       -0.2027

                           S[O.sub.4]
                          ([C.sub.3v])

[E.sub.HOMO]                -0.3772
[E.sub.LUMO]                -0.3233
[[DELTA].sub.[epsilon]]     -0.0539

TABLE 6: The thermodynamic data of the S+[O.sub.4] reaction on the
singlet PES at the B3LYP method (kcal/mol).

Reaction                      [DELTA][E.sup.0]   [DELTA][H.sup.0]

[sup.1]S+[sup.1][O.sub.4]         -67.694            -68.286
  [right arrow] OOSOO
[sup.1]S+[sup.1][O.sub.4]         -170.708           -170.708
  [right arrow] S[O.sub.2]+
  [sup.1][O.sub.2] ([sup.1]
  [[DELTA].sub.g])
[sup.1]S+[sup.1][O.sub.4]         -187.450           -188.042
  [right arrow] S[O.sub.4]
  ([C.sub.2v])
[sup.1]S+[sup.1][O.sub.4]         -164.265           -164.857
  [right arrow] S[O.sub.4]
  ([C.sub.3v])

Reaction                      [DELTA][G.sup.0]   T[DELTA][S.sup.0]

[sup.1]S+[sup.1][O.sub.4]         -62.177             -6.108
  [right arrow] OOSOO
[sup.1]S+[sup.1][O.sub.4]         -173.338             2.629
  [right arrow] S[O.sub.2]+
  [sup.1][O.sub.2] ([sup.1]
  [[DELTA].sub.g])
[sup.1]S+[sup.1][O.sub.4]         -179.147            -8.896
  [right arrow] S[O.sub.4]
  ([C.sub.2v])
[sup.1]S+[sup.1][O.sub.4]         -156.922            -7.935
  [right arrow] S[O.sub.4]
  ([C.sub.3v])
COPYRIGHT 2014 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Khademzadeh, Ashraf; Vahedpour, Morteza; Karami, Fereshte
Publication:The Scientific World Journal
Article Type:Report
Date:Jan 1, 2014
Words:7254
Previous Article:Study of water absorption in Raffia vinifera Fibres from Bandjoun, Cameroon.
Next Article:Pulpal thermal changes following Er-YAG laser debonding of ceramic brackets.
Topics:

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters