Next Article in Journal
Study of Ammonia Adsorption on Magnetite Surfaces with Molecular Dynamics Simulations
Next Article in Special Issue
Investigation of the Anti-Inflammatory Properties of Bioactive Compounds from Olea europaea: In Silico Evaluation of Cyclooxygenase Enzyme Inhibition and Pharmacokinetic Profiling
Previous Article in Journal
Efficient Transformation of Water Vapor into Hydrogen by Dielectric Barrier Discharge Loaded with Bamboo Carbon Bed Structured by Fibrous Material
Previous Article in Special Issue
Identification and Evaluation of Olive Phenolics in the Context of Amine Oxidase Enzyme Inhibition and Depression: In Silico Modelling and In Vitro Validation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinetic Study and Reaction Mechanism of the Gas-Phase Thermolysis Reaction of Methyl Derivatives of 1,2,4,5-Tetroxane

by
Alexander G. Bordón
1,
Mariela I. Profeta
1,
Jorge M. Romero
1,
María J. Jorge
1,
Lilian C. Jorge
1,
Nelly L. Jorge
1,
C. Ignacio Sainz-Díaz
2,*,
Juliana Cuéllar-Zuquin
3,
Daniel Roca-Sanjuán
3,
César Viseras Iborra
2,4,
André Grand
5 and
Alfonso Hernández-Laguna
2,*
1
Laboratorio de Investigaciones en Tecnología Ambiental, Área de Química Física Facultad de Ciencias Exactas y Naturales y Agrimensura, Universidad Nacional del Nordeste, Av. Libertad 5460, Corrientes 3400, Argentina
2
Instituto Andaluz de Ciencias de la Tierra, Consejo Superior de Investigaciones Científicas (CSIC), Av. Las Palmeras 4, 18100 Armilla, Granada, Spain
3
Instituto de Ciencia Molecular, Universitat de València, Apartado 22085, 46071 Valencia, Spain
4
Department of Pharmacy and Pharmaceutical Technology, Faculty of Pharmacy, University of Granada, Campus de Cartuja s/n, 18071 Granada, Spain
5
Unité de Formation et de Recherche, Université Grenoble Alpes, Commissariat á l’Énergy Atomique (CEA), Centre National de la Recherche Scientifique (CNRS), Institute for Nanoscience and Cryogenics-Systèmes Moléculaires et Nanomatériaux por l’Énergieset la Santé (INAC-SyMMES), 38000 Grenoble, France
*
Authors to whom correspondence should be addressed.
Molecules 2024, 29(14), 3274; https://doi.org/10.3390/molecules29143274
Submission received: 27 May 2024 / Revised: 24 June 2024 / Accepted: 3 July 2024 / Published: 11 July 2024

Abstract

:
Tetroxane derivatives are interesting drugs for antileishmaniasis and antimalaric treatments. The gas-phase thermal decomposition of 3,6,-dimethyl-1,2,4,5-tetroxane (DMT) and 3,3,6,6,-tetramethyl-1,2,4,5-tetroxane (acetone diperoxide (ACDP)) was studied at 493–543 K by direct gas chromatography by means of a flow reactor. The reaction is produced in the injector chamber at different temperatures. The resulting kinetics Arrhenius equations were calculated for both tetroxanes. Including the parent compound of the series 1,2,4,5-tetroxane (formaldehyde diperoxide (FDP)), the activation energy and frequency factors decrease linearly with the number of methyl groups. The reaction mechanisms of ACDP and 3,6,6-trimethyl-1,2,4,5-tetroxane (TMT) decomposition have been studied by means of the DFT method with the BHANDHLYP functional. Our calculations confirm that the concerted mechanism should be discarded and that only the stepwise mechanism occurs. The critical points of the singlet and triplet state potential energy surfaces (S- and T-PES) of the thermolysis reaction of both compounds have been determined. The calculated activation energies of the different steps vary linearly with the number of methyl groups of the methyl-tetroxanes series. The mechanism for the S-PES leads to a diradical O···O open structure, which leads to a C···O dissociation in the second step and the production of the first acetaldehyde/acetone molecule. This last one yields a second C···O dissociation, producing O2 and another acetone/acetaldehyde molecule. The O2 molecule is in the singlet state. A quasi-parallel mechanism for the T-PES from the open diradical to products is also found. Most of the critical points of both PES are linear with the number of methyl groups. Reaction in the triplet state is much more exothermic than the singlet state mechanism. Transitions from the singlet ground state, S0 and low-lying singlet states S1–3, to the low-lying triplet excited states, T1–4, (chemical excitation) in the family of methyl tetroxanes are also studied at the CASSCF/CASPT2 level. Two possible mechanisms are possible here: (i) from S0 to T3 by strong spin orbit coupling (SOC) and subsequent fast internal conversion to the excited T1 state and (ii) from S0 to S2 from internal conversion and subsequent S2 to T1 by SOC. From these experimental and theoretical results, the additivity effect of the methyl groups in the thermolysis reaction of the methyl tetroxane derivatives is clearly highlighted. This information will have a great impact for controlling these processes in the laboratory and chemical industries.

Graphical Abstract

1. Introduction

Organic peroxides are a series of compounds with many different applications, due to their therapeutic properties and explosive character. Some of them are used to fight malaria in subtropical regions, and as antibiotic; some others therapeutical properties are the intervention in biological metabolites, cancer illness, cellular aging, and herbicide use [1,2,3,4,5].
Furthermore, commercial vinyl polymers—such as polystyrene, polyethylene, and polyvinyl chloride—are produced industrially through free radical polymerization, where organic peroxides are commonly used as free radical initiators. The thermal decomposition of organic diperoxides through the homolytic cleavage of the peroxide bond generates a biradical intermediate, which would initiate the polymerization of vinyl monomers. This biradical species would be incorporated into the polymer structure, giving polymeric species that contain peroxide groups in their structures, which are then decomposed in the next processes. From the structure of the organic peroxide, a biradical is generated, which influences on the properties of the polymer. Therefore, it is important to clarify the decomposition reaction mechanism of the di- and triperoxide organic compounds, which can be considered multifunctional initiators due to the presence of more than one peroxide group [1]. However, the knowledge of the behavior and reaction mechanism of these compounds is limited. Hence, a further understanding of these processes would allow the control of many industrial processes with high impact in economic costs, and environmental effects.
The tetroxane-derivative antimalaric drugs are compounds with large molecular weight, low toxicity and no genetic toxicity. It is well known that the diperoxide compounds are more suitable as drugs than the equivalent ozonides and trioxane compounds [2]. Therefore, the diperoxide compounds are considered as new anti-malaric drugs and can be a good choice in front of artemisinin derivative drugs [2,3,4]. A deeper knowledge of the reactivity of these compounds close to receptor environment would help in the design of optimal antimalaric drugs.
Diperoxide compounds are also found with an important herbicide activity [6]. Many tetroxanes have got phyto-toxicity, inhibiting the growing of roots, especially cucumber and sorghum [6]. This activity is statistically comparable to the glyphosate and imazapyr [6]. One of the compounds of this family (3,3,6,6-tetramethyl-1,2,4,5-tetroxane—in short, either acetone diperoxide (ACDP) or tetramethyl-tetroxane), which is of interest in the present work, has shown especially high phytotoxicity.
In previous papers, the experimental kinetic and theoretical mechanism of the thermolysis reaction at gas phase of the parent compound of the series (1,2,4,5-tetroxane (formaldehyde diperoxide, FDP)) was accomplished [7]. Additionally, the theoretical mechanism of the thermolysis reaction of 3-methyl-1,2,4,5-tetroxane (MFDP) [8] and 3,6-dimethyl-1,2,4,5-tetroxanes (DMT) [9] was also determined. Previously, a concerted mechanism of the thermolysis reaction was proposed [10], but from quantum mechanics calculations, a stepwise mechanism was also determined. In these previous works [7,8,9], a common stepwise mechanism for the thermolysis reaction was found with lower activation energy than the previously proposed concerted mechanism. In the stepwise mechanism, the first step brings the reactant to a diradical open (o) intermediate in its singlet state (S0), and, afterwards two step-wise mechanisms follows: (i) a decomposition step by step from the o structure to two aldehydes molecules plus an oxygen molecule in its excited singlet state (S1); and (ii), a transition from the S-o diradical structure to an o structure in its triplet state (T1), and a subsequent triplet stepwise mechanism, producing two molecules of aldehydes plus one oxygen molecule in its ground triplet state (T1). This singlet-to-triplet-state non-adiabatic transition is deeply studied in this work. In general, the rate-limiting step was the second step. The thermal decomposition mechanism in dioxetane was also found to proceed in a stepwise manner and with a non-negligible nonadiabatic transition from the singlet to triplet state [11].
The methyl group substitution effect was analyzed previously [9], showing that the activation energies (Ea), reaction energies (Er), and reaction coordinate lengths of the different steps of the stepwise mechanism of the reaction behave linearly with the number of methyl groups.
Because of the peroxide bonds and the cyclic structure, these compounds become highly explosive. In general, they are compounds of an unstable and explosive nature due to the peroxide groups, which are intrinsically unstable and highly reactive [5], so the synthesis and manipulation of these materials require strict care to avoid accidents. The multiple incidents recorded in recent decades involving this type of species have raised the need for a constant and exhaustive physicochemical study. Methyl derivatives are especially dangerous because of their high explosivity, and their synthesis and manipulation should be performed in very small amounts. Additionally, the products formed during the thermolysis reaction of the tetroxane derivatives augment three times the number of moles with respect to the reactants, contributing also to their high-explosive character. The successive methylation of the tetroxanes yields different reaction heats. These features plus the surplus incentive of knowing the kinetic and mechanism of thermolysis reaction of the highest and lowest methyl-substituted members of the series and its comparison with the parent compound increase the interest of this work. The knowledge of these complex non-adiabatic processes will facilitate the disposition of control techniques for further synthesis and manufacturing scale processes. This will have a great economic and environmental impact in the chemical industry.
Herein, the following goals are pursued: (i) to carry out an experimental kinetic characterization of DMT and ACDP of the thermolysis reaction in the gas phase; (ii) to determine the theoretical thermolysis reaction mechanism of trimethyltetroxane (TMT) and ACDP; and (iii) to theoretically revisit tetroxane and its methyl derivatives in order to study the non-adiabatic transition occurring in the diradical open structure (chemiexcitation) by multiconfigurational quantum chemistry. The obtained data will be studied in order to know the methyl substitution effect in the gas phase thermolysis reaction of the methyl-tetroxanes.

2. Results and Discussion

2.1. Analytical Results

After the introduction of ACDP [10] into the injector chamber [12], only three compounds were detected: solvent (retention time = 1.3 min), acetone (2.84 min), and the remaining non-reacted ACDP (8.7 min). As a consequence of the very short residence time (40 s) of ACDP in the injector chamber, acetone was the only product detected, and no other products coming from possible secondary reactions were found.
The introduction of DMT into the injector chamber showed the detection of three compounds: acetaldehyde (1.4 min), the solvent (3.2 min), and the remaining non-reacted DMT (7.2 min). Acetaldehyde was the only product detected, a consequence of the very short residence time (40 s) of DMT in the injector chamber.

2.2. Kinetic Results

According to Arrhenius’s equation, the Ea was calculated from the slope of least-square fitting of ln[ln(C0/C)] as a function of 1/T and the intersect of this fitting yields the ln(A.t). Considering the retention time, which is constant in all measures, A can be obtained.
The rate constant of reaction, kexp, for DMT and ACDP at different temperatures (Table 1), gives the following Arrhenius equations:
ln kexp,DMT = (25.66 ± 0.80) (26.50 ± 1.00)/(RT); R2 = 0.999
ln kexpACDP = (22.48 ± 0.85) − (22.7 ± 1.00)/(RT); R2 = 0.999
Therefore, the Ea values are 26.5 ± 1.0 and 22.7 ± 1.0 kcal/mol for DMT and ACDP, respectively, and 1.4 × 1011 and 5.8 × 109 s−1 for the frequency factors of DMT and ACDP, respectively. In a previous work, Ea for ACDP turned out to be 39 ± 2.5 kcal/mol [13], coming from a very different technique, where many reactors (Pyrex vessels) were introduced at a constant temperature (at different temperatures), and the reaction was stopped at different times (up to 50 h) without controlling that all species were in the gas phase. Other cyclic peroxides, such as dioxetane (Ea = 22.7 kcal/mol) [11] dimeric diphenyl tetroxane [12] and FDP [7] showed Ea = 21.6 and 29.3 kcal/mol, respectively. From these results, it is possible to observe a substitution effect on the Ea of this series of compounds. The methyl substitution from FDP to dimethyl tetroxane (DMT) produced a decrease of Ea, where the gas phase of DMT is made of a mix of the different position isomers [9]. It is possible to infer that the methyl substitution would go on affecting the Ea of these derivatives. Ea values for these three derivatives (FDP, DMT, and ACDP) are fitted as a function of the number of methyl groups (n), resulting in the following linear function:
Ea = (29.28 ± 0.04) − (1.37 ± 0.01)n; R2 = 0.99
This indicates an average decrease of 1.4 kcal/molCH3 for Ea. The average value of the frequency factor of FDP was 5.2 × 1013 s−1. Our frequency factor shows a much smaller value, having a parallelism with the methyl substitution of the dioxetanes [14,15], which were in the range of 2 × 1012–1.3 × 1013 s−1. Nonetheless, a new linear relationship is also found with the number of methyl group:
lnA = (31.10 ± 1.1) − (2.3 ± 0.4)n; R2 = 0.97
Therefore, we have to emphasize the additivity of the methyl groups in the thermolysis reaction of the methyl tetroxane derivatives, which could bring further to think that the additivity of the methyl groups in the properties of these compounds would be strongly feasible. The hyperconjugative effect of methyl groups increases the delocalization of the electronic cloud stabilizing the transition states and decreasing the activation energy.
Previous experimental and theoretical works [7,8,9] indicated a stepwise mechanism. In this case, we are going to undertake a computational study in order to confirm and widen the previous proposed mechanisms.

2.3. Thermolysis Reaction Mechanism

The possible reaction mechanisms are given in Scheme 1 and analyzed with our calculations, in which we can describe a concerted mechanism (S-I) and a stepwise mechanism (S-II) via a diradical open ·O···O· structure.

2.3.1. Concerted Mechanism

The tri/tetra-methyl-substituted-tetroxane reactants are in a chair conformation (S-c). The H atoms of the methyl groups in geminal position show a symmetric eclipsed conformation two to two with respect to the average plane of the molecule, in which all carbons are included (Figure 1a). One H atom of each methyl group is in this plane, and two H atoms out of the plane of two geminal methyl groups are facing each other. The concerted mechanism has been calculated as one single step, where the transition state (TS) is called S-TSX (Table 2, Figure 1a), with a transition vector (TV) of opening and closing the peroxide O-O bonds and opening C···O bonds to form the O2 molecule and two molecules of acetaldehyde/acetone via the opening/closing of the C-O bonds. This TS and TV are depicted in Figure 1a. The Ea values are 66.3, 68.8, and 65.9 kcal/mol for TMT-axial (TMTax), TMT-equatorial (TMTeq), and ACDP, respectively, being lower than the less-substituted compounds of the series [7,8,9]. The TMTax has lower Ea than TMTeq probably due to the anomeric effect that stabilizes the axial isomer TS.
Taking into account all methyl derivatives, a correlation of the average values of Ea as a function of the number of methyl groups has been performed (Figure 1b). We observe that the Ea decreases linearly with the increasing number of methyl groups, with a slope of −2.04 kcal/molCH3 (Figure 1b; in this figure, average values of axial and equatorial isomers of TMT are also considered). This slope is larger than that of the experimental Ea as a function of the number of methyl group (Equation (3)). This fitting shows that the FDP Ea for a concerted mechanism would be 73.8 kcal/mol, being considerably higher than the experimental value (29.3 kcal/mol) [7]. Additionally, the value of Ea = 65.9 kcal/mol for ACDP is the lowest of the series, which is much higher than our experimental value but is still high enough to avoid the concerted mechanism in front of the stepwise mechanism (vide supra).
The reaction product (S-P) (Branch S-I, Scheme 1) shows an O2 molecule in the singlet state surrounded by hydrogen bonds of two acetaldehyde/acetone molecules depending on if the reactant is TMT or ACDP. Reaction energy (Er) is exothermic, being −31.1 and −28.0 kcal/mol for TMT axial and equatorial isomers, respectively, and −32.8 kcal/mol for ACDP. The series follows a linear equation as a function of the number of methyl groups: Er = −17.1 − 4.0n (R2 = 0.98) (FDP is not included in this equation).

2.3.2. Stepwise Mechanism in the Singlet State

The stepwise mechanism is depicted in Scheme 1 (Branch S-II). The mechanism has three steps: (i) the production of the open diradical structure (S-o), (ii) the formation of acetone/acetaldehyde (S-b′) and oxide-peroxide intermediates (S-b), and (iii) the generation of the end products (S-p), where the oxygen molecule is in the singlet state. Step (i) is common to all methyl derivatives of tetroxane. In Step (ii), only one possible reaction path is possible for ACDP, but in the case of TMT, two different secondary reaction pathways are possible (vide supra).
Let us start with the reactant, S-c, which goes throughout one TS (S-TSco) to the open diradical structure S-o (Figure 2). This S-TSco has a much lower Ea than the above S-TSX. The Ea for S-TSco of TMT turns out to be 17.7 and 18.3 kcal/mol for axial and equatorial methyl position isomers, respectively (Table 3 and Table 4); and for ACDP Ea = 18.5 kcal/mol (Figure 2, Table 5). Considering the Ea as an average value (Eaaver) of the axial and equatorial isomers, a function with the number (n) of methyl groups of the different derivatives of the series can be expressed by a linear equation Eaaver = 16.15 + 0.61n, R2 = 0.993. The most important distance in the TS, and so the simplified reaction coordinate, in the structure of S-TSco, is the O···O distance, which, in the average of the axial and equatorial isomers (dO···Oaver), also shows a linear function with the number of methyl groups: dO···Oaver = 1.90 + 0.014n, R2 = 0.96. Indeed, the Ea and dO···O as a function of the number methyl groups show qualitative parallel trends, increasing the distance as increasing Ea, but with a much lower slope of 0.014 Å/CH3, indicating that the main coordinate of the TS has a small variation on the effect of substituents in the same reaction.
The first intermediate in the reaction is the diradical open structure, S-o, with a distance between both O atoms of 3.22 Å and a spin density of −0.94 and 0.94 e (Figure 2), indicating that both electrons have different spins, maintaining the identity of the singlet state. This intermediate gives endothermic reaction energy, increasing with the number (n) of methyl groups, following also a linear function: Erco = 9.0 + 1.22n; R2 = 0.992. The slope of Erco is twice the slope of Ea-TSco, indicating that the methyl substitution effect is twice as large in the diradical open intermediates as in TSco. Indeed, the linear equation relating both Erco and Ea (TSco) is Erco = −23.5 + 2.0Ea (R2 = 0.996). The O···O distance of S-o follows the linear equation dO···O = 3.519 − 0.074n, R2 = 0.997. Additionally, the reaction coordinate O···O of S-o intermediate decreases with a slope of 4.6 Å with respect to this O···O distance in S-TSco, following the linear equation dO···Oo = 12.2 − 4.6d(O···O)TSco (R2 = 0.99). All these effects agree with the Leffler–Hammond postulate [16,17].
The next step is the TS of a broken C-O bond, S-TSob′, which depends on either the two-methyl- or single-methyl-substituted C atom at TMT (Figure 3). In the first case, acetone will be obtained as the product (S-b′, Branch SII-ii of Scheme 1), while in the second case, the mono-substituted carbon (R = H, CH3 in Scheme 1) will give acetaldehyde as the first product of reaction (S-b′, Branch SII-i of Scheme 1). Let us start with the mono-substituted carbon (Branch SII-i, Scheme 1). In Figure 3a the S-o intermediate goes throughout S-TSob′ to S-b′ (second intermediate product). The Ea of S-TSob′ is 17.1 and 17.5 kcal/mol for axial and equatorial isomers with respect to the S-o intermediate. In S-TSob′, the distinction between equatorial and axial isomers of CH3 positions is not clear. Nonetheless, some slightly structural differences are found in the potential energy surface (PES) critical points. For this reason, average values are also taken. The transition vectors are mainly formed by the vibrations of C···O distance—which are 1.704 and 1.707 Å for axial and equatorial isomers, respectively—and the spin density is shared by all oxygen atoms. The product is an acetaldehyde molecule plus an oxide–peroxide diradical, S-b′ (Figure 3), whose Er is exothermic with respect to the S-o intermediate. However, it is slightly endothermic with respect to S-c.
When acetone is left as the first product in TMT (Branch ii in Scheme 1), the corresponding Ea of TSob′ and Erob′ are 15.4 and −13.9 kcal/mol, respectively, for the axial isomer and 15.4 and −12.7 kcal/mol, respectively, for the equatorial isomer (Figure 3b). This reaction is more favorable than the above, where the aldehyde is left out, though the differences are not meaningful. ACDP yields acetone and a peroxyacetone diradical as reaction products in the second stage, giving 15.31 and −14.83 kcal/mol for Ea and Er, respectively.
Considering the Ea as an average value (Eaaver) of the axial and equatorial isomers and the two possible intermediates formed in this step of the different derivatives of the series, a linear function can be obtained. These Ea values low down to 1.3 kcal/molCH3, following the linear equation Ea(TSob′) = 20.6 − 1.3n (R2 = 0.93). The C···O distances have a linear function with respect to the number (n) of methyl groups: d(C···O) = 1.714 − 0.0013n (R2 = 0.99). The slope is quite small, indicating that this TS is quite constant with respect to the number of methyl groups. The average values of axial and equatorial isomers give an exothermicity of 2.5 kcal/molCH3 with respect to S-o, resulting from the linear equation Erob′ = −5.0 − 2.5n (R2 = 0.98). The correlation is low because the product’s disposition is not very regular. In general, the Ea of this reaction step decreases with the number of methyl groups, and Er is growing in exothermicity. These facts indicate the additivity of methyl groups in our methyl-tetroxane derivatives; that is, each methyl group contributes to decreasing the activation energy.
From now on, one of the products of the step (either acetaldehyde or acetone), S-b′, is removed, and the oxide–peroxide diradical intermediate (S-b) continues reacting. A second molecule of acetaldehyde or acetone and one oxygen molecule will be the products of this reaction step (Table 3, Table 4 and Table 5). When S-b is ·O-CH(CH3)-O-O·, TMT has the same Ea of 11.6 kcal/mol (as average value). This value is similar to MFDP, and DMT [8,9]. The Er of the acetaldehyde product (Scheme 1, Branch SII-ii) is −29.01 and −28.64 kcal/mol with respect to b for TMTax and TMTeq, respectively (Figure 4a, Table 4). However, when the acetone is the second product, S-b is ·O-C(CH3)2-O-O· (Scheme 1, Branch SII-i), where the Ea and Er are 10.73 and −31.10 kcal/mol, respectively, for TMT and ACDP (Figure 4b). Therefore, the reaction of TMT and ACDP in the singlet state is exothermic in the last step of reaction.
In the derivatives with three and four methyl groups, the largest TS is that of the first step—that is, S-TSco, which could be the rate-limiting step. On the contrary, in FDP, MFDP, and DMT, the rate-limiting step looks like the second step, where S-TSob′ has the highest energy barrier [7,8,9]. However, these energy differences are very small and both steps can be considered limiting step at experimental conditions. Our experimental Ea for ACDP obtained in this work is 22.7 kcal/mol, which is consistent with our theoretical value of 18.5 kcal/mol of S-TSco (Table 5).

2.3.3. Thermolysis Reaction as Triplet State (Scheme 2)

For ACDP, a quasi-chair excited structure (T-qc) in the PES of the triplet state has been found at 14.9 kcal/mol above S-c. This structure is not as symmetric as S-c and has a dO···O = 2.82 Å, much larger than the other peroxide bonds (1.41 Å). The methyl groups are in the same conformation as S-c. This excited structure has a spin density of 0.942 e on the two open oxygen atoms. The equivalent structure for TMT has not been found. The second structure found in the T-PES is the open diradical structure (T-o) close to S-o. This T-o presents an energy of 13.7 kcal/mol with respect to S-c and only −1.2 kcal/mol with respect to T-qc. Then, the energy of T-o is very close to T-qc, where the quasi-open structure of T-qc has a large dO···O distance similar to the open structure T-o. If a transition state exists from T-qc to T-o, we could not find it, but it should be with low Ea from T-o to T-qc.
The energy and structure of T-o are very close to the S-o, having a spin density of 0.939 e on each open oxygen of the peroxide bond, and an O···O distance of 3.222 Å for ACDP, and 3.286 and 3.330 Å for TMT axial and equatorial isomers, respectively (Figure 5). The O···O distances decrease as the number (n) of methyl groups following the linear equation dO···O = 3.527 − 0.075n (R2 = 0.998). This equation and distances are very close to the S-o-dO···O distances, indicating the similar o structures for the different derivatives and their very close variation with the number of methyl groups.
From T-o a transition state T-TSob′ is also found. This new structure is the TS for breaking and forming the first acetone/acetaldehyde molecule. In TMT, different values are found depending on if either acetone or acetaldehyde is the first product in b′. If acetaldehyde is the first product in b′ (Scheme 2, Branch i), Ea is 16.98 and 15.13 kcal/mol for the axial and equatorial isomers, respectively (Figure 5a). If acetone is the first product (T-b′, Scheme 2, Branch ii), Ea is 15.20 kcal/mol for both axial and equatorial isomers (Figure 5b). In ACDP, Ea is 15.13 kcal/mol with respect to T-o (Figure 5b). The Ea as a function of the number of methyl groups from FDP to DMT was previously calculated with a quadratic function [9]. With the addition of the three and four substituted derivatives to the system, the quadratic function is Ea = 25.3 − 6.2n + 0.9n2 (R2 = 0.97), similar to that found in Reference [9] (Table 6, Table 7 and Table 8).
C···O distances for T-TSob′ are 1.708 Å for ACDP and 1.704 and 1.715 Å for TMT as average values of the position axial/equatorial isomers when acetaldehyde and acetone are the first products, respectively. These distances are not quite different from those of S-TSob′. TV (Figure 5) is viewed as a vibration of the C···O breaking/forming the C···O bond. These distances vary with the number (n) of methyl groups with a linear function: dC···O = 1.714 − 0.0014n (R2 = 0.991). The slope of the function is very small, where the structure of the TS can be considered very stable as the number of methyl groups increases. This equation is very close to the dC···O equation of the S-TSob′, which indicates an identity of structures in the TSob′ of the diradical part of reaction independently of the multiplicity of the state. Additionally, the o diradicals also have very close structures for both multiplicity states, as already mentioned above. From this, we can extract that if both multiplicity states have equal reactants and TS’, the reaction path should be very close. The spin density of T-TSob′ is shared over all oxygen atoms in the structure (Figure 5).
T-TSob′ goes to T-b′, a super-molecular system with an acetone/acetaldehyde molecule and a peroxy-acetaldehyde/acetone diradical, with a Cacetone/acetaldehyde···Operoxy distance of 3.956 and 3.473 Å for TMTaver(ax-eq) and acetone as the first product of the reaction; 3.647 Å is the average of both position isomers when the acetaldehyde is the first product of the reaction. The reaction energy is −11.06 kcal/mol with respect to T-o for TMT as an average value of both position isomers when acetaldehyde is the first product of reaction (Figure 5a, Table 6); and −13.36 kcal/mol for the equivalent value when acetone is the first product of the reaction (Figure 5b, Table 7). The Er as a function of the number of methyl groups behaves decreasing with the linear function Er = −5.5 − 2.9n, R2 = 0.93. The spin density of the triplet state is shared on the three oxygen atoms of the peroxy-acetone radical. From now on, the acetone molecule is removed from the super-molecular system and the critical point is optimized from the peroxy-acetone radical. Nonetheless, the energy of acetone molecule is constantly added up to the energy of the critical points in order to preserve the whole energy of the system. Following this diradical intermediate, T-b, a T-TSbp yields the final product (T-p, Figure 6). The Ea is really the lowest of the T-PES and S-PES with values of 6.88–7.90 kcal/mol (Table 6, Table 7 and Table 8). The Ea as a function of the number of methyl groups is also linear: Ea = 9.48 − 0.67n, R2 = 0.98. In the T-PES all Eas low down with the increasing number of methyl groups, indicating the reaction rate is faster with the methyl substitution. Besides, this TSbp is the lowest TS of the S- and T-PES. The last TS gives rise to the second acetaldehyde/acetone molecule plus the 3O2 molecule. The TVs clearly show this formation (Figure 6a,b). The spin density also is shared over the three oxygen atoms of the system. The highest Ea on T-PES is that of T-TSob′, being the rate-limiting step of the reaction, which agrees with the other derivatives of the series.
Finally, the products are 3O2 plus the second acetone/acetaldehyde molecule. The reaction energy with respect to T-o is exothermic, with −58.2 kcal/mol for ACDP (Erob′+ Erbp, Table 8) and −54.2 and −53.9 kcal/mol for TMTaver(ax-eq) for acetone (Table 6) and acetaldehyde (Table 7) as the second product, respectively. A linear function is also found for this step of reaction, Erop = −41.4 − 4.3n, R2 = 0.998, describing the reaction energy of op, giving an exothermicity of 4.3 kcal/molCH3. Considering the first step of the reaction in singlet state for obtaining the biradical S-o, the reaction energies of the reactionin triplet state with respect to reactant S-c are −44.7 kcal/mol for ACDP and −41.4 and −40.7 kcal/mol for TMTaver(ax-eq) for acetone and acetaldehyde as final products, respectively. If the reaction would proceed in singlet state PES, the Er of the equivalent critical point for S-p, with respect to S-c is −24.0 kcal/mol and −29.9 kcal/mol for TMTaver(ax-eq) for acetone and acetaldehyde as final products, respectively and −32.4 kcal/mol with respect to S-c for ACDP (Table 3, Table 4 and Table 5). This indicates that the reaction energies on T-PES are much more exothermic than in the S-PES, as expected based on the fact that the singlet molecular oxygen is more unstable than triplet molecular oxygen. The results in the T-PES are more in agreement with the experimental reaction energies of the diperoxide compounds.
The influence of the methyl groups in the thermolysis reaction of tetroxane derivatives increases the exothermicity of the reaction, and they change the rate-limiting step from the second step to the first step, although the differences in Ea are very small, and both TSs could be viewed as limiting steps at experimental conditions; additionally, they also change the geometry of intermediates and transition states in an approximately linear way. The singlet and triplet states present very close structures at the different critical point intermediates of the reaction. These similarities of structures indicate a possible nonadiabatic transition from the singlet to triplet states in such a way that reaction would be more exothermic if this transition of different multiplicity states is produced, which would agree with the high exothermicity of the tetroxanes. Therefore, it is worth providing general details on the intersystem crossing process, as described in the next section.

2.4. Excited States and Non-Adiabatic Chemistry of Tetroxane and the Methylated Systems

The lowest-lying singlet and triplet excited states of tetroxane (FDP) at the reactant structure correspond to electronic excitations from the lone pairs of the oxygen atoms (ni) or the σ bonding orbitals of the -O-O- bonds (σi) to the σ anti-bonding orbitals of those bonds (σ*j)—that is, niσ*j or σiσ*j excitations (see orbitals in Figure 7). These electronic configurations are associated with the weakening of the peroxide bonds. The vertical energy position of the three first singlet (Si) excited states (all of them niσ*j in nature) relative to the energy of the ground state (S0) ranges from 5.09 to 5.78 eV (Table 9). The manifold of triplet states (also niσ*j in nature) is lower-lying than that of the singlet states, being in a range of 3.87–4.81 eV for the first four of them. Test computations show that σiσ*j configurations appear on the geometry of the reactant at higher energies (Table 9).
Upon -O-O- bond stretching, the mentioned electronic configurations become stabilized as is confirmed by the calculated energy of the excited Si and Ti states at the diradical geometry. Here, the four singlets and four triplet states appear in an energy range lower than 1 eV (Figure 8). T1 has the same energy as the ground state (S0), and other states (S1, S2, T2, and T3) are located only 0.3 eV above S0. Considering that the kinetic energy released after the transition state of the -O-O- bond breaking is of the same order; those excited states may be populated in the non-adiabatic dynamical process, which follows the peroxide bond decomposition.
To analyze the efficiency of the non-adiabatic transitions at the diradical region, SOCs were calculated between the different singlet and triple states (Table 10). The coupling between S0 and T1 is negligible, and the maximum value is obtained with T3, which is located at 0.29 eV. Additionally, a spin orbital coupling (SOC) between S2 and T1 at 0.29 eV also has a very high value. Both SOCs indicate two possibilities for transferring the population from S0 to T1 along the thermal decomposition: (i) S0 transfers the population initially to T3 via strong SOC and subsequently by an ultrafast internal conversion process T3 decays to T1, or (ii) S0 transfers the population initially to S2 via an internal conversion process and next the molecule changes to the T1 manifold through strong SOC between S2 and T1.
Analysis of the open-shell spin densities at the diradical region shall help to provide a rationale on the distinct efficiencies related to the population transfer between the singlet and triplet states (Figure 9). Such densities reflect the distinct occupation of the unpaired electrons of the diradical in the oxygen 2px and 2py atomic-like orbitals, which are perpendicular to the C-O bonds. As can be seen, high SOC values are obtained when singlets and triplets have one unpaired electron in a perpendicular 2p orbital, in accordance with El-Sayed’s rules [18].
Once we comprehend the electronic-structure features of the parent tetraoxane molecule (FDP), it is worth extending the analysis to the methylated molecules of the family, Figure 8 shows that ACDP has the electronic states closer in energy compared with FDP. Energies of S3 and T4 show a clear trend of decreasing energy upon methylation. Regarding the most accessible states (S1, S2, S3, T2, and T3), we see that an odd number of methyl group substitutions (MFDP and TMT) increases the energy gap, while those states are energetically slightly more accessible for FDP and the molecules with an even number of methyl groups. Comparison of SOCs in the distinct molecules of the family (Table 10) also do not show an efficient transition from S0 to T1 at the diradical region and a higher efficiency for S0 -> T3 and S2 -> T1. The highest values appear also for the molecules with an even methylation. ACDP, which shows the lowest energy position of the S2 and T3 states at the diradical region and high SOCs arises as the most efficient system in the family for non-adiabatic decomposition. The dynamical effects, not considered in this work, are expected to favor this triplet population even more, as inferred from the study by Vacher et al. comparing the dynamics of 1,2-dioxetane and the tetramethylated compound [19]. By increasing the mass of the substituents, the time spend at the diradical region increases, allowing a more effective re-distribution of the population among the near-degenerate singlet and triplet states and accordingly higher yields for the evolution on T1 [19]. These transitions are more favorable for FDP and even-methylated derivatives (the highest SOCs and the least ΔE), being less efficient at the odd-methylated compounds. For FDP and the even-methylated derivatives, SOC is linear with the number of methyl groups (SOCS2->T1 = 96.15–0.96·n; R2 = 0.999; SOC of DMT being the average of the three position isomers); SOC of S0 to T3 equally decreases with the number of methyl groups, but non-linearly. ΔE (Figure 8) of the tetramethyl derivative is the least of all derivatives. SOCs of both transitions of the two odd-methylated derivatives increase with the number of methyl groups, SOCs at S2 to T1 being the highest.

3. Materials and Methodology

3.1. Materials

DMT was prepared by reacting 0.6 mL of acetaldehyde, which is added to a stirred and cooled solution at 263.15 K containing 2.5 mL of ethanol (absolute; Merck, Darmstadt, Germany), 0.3 mL of hydrogen peroxide (56% v/v) and 1.5 mL of sulfuric acid (98%, Merck). The mixture was stirred, and sulfuric acid was added up to 3 mL, precipitating a white solid. Then, water was slowly added to complete a volume of 5 mL in order to increase the insolubility of the peroxide precipitate in the alcoholic medium. A white precipitate with a gelatinous characteristic was obtained, which was washed and centrifugated repeatedly with distilled water until neutral pH. Thus, the wet solid was dissolved directly in benzene (avoiding the incorporation into the organic phase of possible polar impurities). This dissolution was analyzed using the gas chromatographic (GC) method by using an Agilent 7890-A chromatograph (Santa Clara, CA, USA).
ACDP was synthesized by carefully adding acetone to a stirring solution of hydrogen peroxide in strong acid medium at 253 K, as it was previously proposed [7]. ACDP solid was characterized by IR spectroscopy by using a Nicolet spectrometer (Thermo Fisher Scientific, Waltham, MA, USA), and a solution was analyzed by UV spectroscopy with a CampSpec M330 spectrophotometer (Leeds, UK) and GC. THF was used as a solvent in the thermolysis reaction of ACDP.

3.2. Analytic and Kinetic Methods

The GC method was used to study the kinetic of the thermolysis reaction of DMT and ACDP. This technique is very suitable for small quantities of samples. The reactor in this method is the injection chamber of a gas chromatograph, where the flow and reaction time (40 s) are maintained constant, and the reaction temperature is a parameter. When the compound flows out of the reactor and goes to the chromatographic column, temperature lowers enough to stop the thermolysis reaction, but temperature is maintained high enough to produce the chromatographic separation [12] (ΔT = −150 K in the most adverse case). The method was also used for FDP, obtaining good results [7].
A solution (1 μL) of DMT (0.001 M) in benzene was introduced in the injection chamber. The reaction was controlled by the injector temperature within the range of 493.15–543.15 K, and the rate of ΔT = 10 K. As the flow and linear rate were kept constant, the reaction time is consequently kept also constant. The reaction products and reactant were analyzed by GC with an electronic device of constant flow (Q = 1.40 mL/min; VL = 20 cm/s) using nitrogen as carrier gas. A flame-ionization detector and a melted silica capillary column (25 m × 0.53 mm) coated with 5% of phenyl-methylpolysiloxane were used. Temperature of the chromatographic column was started at 313 K for 2 min, then programmed at a rate of 30 K/min up to 463 K, and held constant for 10 min.
A solution (1 μL) of ACDP (0.02 M) in THF was introduced in the injection chamber under the same experimental conditions as DMT. The reaction was controlled by the injector temperature.
Residence time, which coincides with the reaction time, in the injector chamber was calculated by Reference [12]:
t = l/v
where l (cm) is the liner length and v is the linear rate of the carrier gas. All previous works have found the first-order reaction model for this kind of thermolysis reaction [7,12]. Then, Ea was calculated by the Arrhenius equation in a first-order integral kinetic equation such as:
ln[ln(C0/C)] = −Ea/RT + ln(A·t)
where C0 and C are the initial and remaining concentrations of the reactant, respectively; T is the temperature; A is the frequency factor of the Arrhenius equation; and t is the residence time, Equation (4). ln[ln(C0/C)] is fitted as a function of 1/T to a linear equation.

3.3. Computational Methods

The electron density properties of the studied molecules (ACDP and TMT) were calculated according to references [7,8,9]. The geometry and energy of the reactant, intermediates, and products were computed by using the Gaussian 09 program [20], at the DFT/BHANDHLYP/6-311 + G** level [21,22,23]. Critical points (minimums and TS) of the S- and T-PES were calculated independently and optimized with the Berny method [24]. The topology of the critical point was characterized by means of the harmonic frequencies. Transition vectors (TVs) for the TSs were inspected and drawn with the help of the GaussView package [25]. The intrinsic reaction coordinate method was used in order to connect the reaction paths between the TS and reactant, products or intermediates [26,27]. The relative energies of the critical points of the PES were calculated by adding the total energy to the zero-point energy (ZPE).
The tetroxane series (FDP, MFDP, DMT, TMT, and ACDP) with their relevant axial and equatorial conformers, as determined in the current and previous related works, was revisited to obtain a high-level characterization of distinct nonadiabatic transitions between T and S states occurring along the thermal decomposition reaction. Four S states (S0–S3) and four T states (T1–T4) were energetically characterized accounting for the spin orbit coupling (SOC) and transition probability between all of them. The CASSCF/CASPT2 method [28,29] as implemented in OpenMolcas [30] was used to compute energies and SOCs at the DFT geometries. The valence double-ζ plus polarization atomic natural orbital L-type basis set (ANO-L-VDZP) was employed. State-average CASSCF wave functions, averaging 4 S and 4 T separately, were used as zeroth-order references for the CASPT2 calculations. The active space comprised 12 electrons distributed in 8 orbitals, CASSCF(12,8). In the CASPT2 computations, the ionization potential electron affinity parameter (IPEA) was set to 0.00 au, and an imaginary level shift of 0.2 au [31] was employed to minimize the presence of weakly intruder states. SOC matrix elements T x , m s = 1,0 , 1 H ^ S O S y were computed using the CASSCF wave functions. The SOC complex vector norms | H ^ S O ( T x , m s = 1,0 , 1 / S y ) | were obtained as the square root of the sum of the absolute square of each m s component of the triplet state T x with the corresponding singlet state S y . The effective one-electron spin–orbit Hamiltonian was derived from atomic mean field integrals, employing the RASSI method [32,33].

4. Conclusions

In this work, experimental activation energy and the pre-exponential factor of the thermolysis reaction of DMT and ACDP have been studied by means of gas chromatography. Thermolysis reactions were produced in the injector chamber at different temperatures. A first order velocity constant was found experimentally. The experimental activation energies are 26.5 and 22.7 kcal/mol for DMT and ACDP, respectively. Additionally, the reaction mechanism of this reaction was accomplished by means of theoretical methods in TMT and ACDP. All features of the mechanism were searched as a function of the number of methyl groups. Position isomers were also considered. Finally, the first S and T excited states of the series were also calculated, and the SOC between them was determined.
The calculated Ea of the stepwise mechanism was lower than the concerted mechanism, and it agrees with the experimental results of the tetroxane derivatives and other similar compounds. The first intermediate is a ·O···O· diradical structure, which is very close in geometry and energy and has similar structure to the triplet diradical. The next step is the C···O dissociation and the production of the first molecule of acetaldehyde/acetone in the singlet state. The transition state for this dissociation is also very similar to that of T state, and the geometry and energy of the open diradical and the transition states are linear with the number methyl groups. If it would be possible to extrapolate the reaction paths, they would be also very close for both multiplicity states.
The reaction energies of both states are exothermic, but those of the T-PES is much more exothermic. A transition from the S to T state in the open diradical could be produced by two possible mechanisms: (i) from the singlet ground state to S2, for internal conversion and from this one to T1 by one strong SOC; and (ii) from S0 to T3 by a strong SOC and from this one to T1, decaying for internal conversion. Such singlet to triplet non-adiabatic transition seems to be more hindered for the mono- and tri-methylated tetroxanes, being more favorable in FDP and the di- and tetra-methylated molecules, in particular the last one.
Therefore, our experimental and theoretical studies emphasize the additivity of the methyl groups in the thermolysis reaction of the methyl-tetroxane derivatives and its properties. This work is a clear example where theoretical calculations can explain the experimental results of the behavior of these tetroxanes, the reaction products, and the mechanisms. The knowledge acquired from this work can be applied to obtain better control of these processes at laboratory and industrial scales, with economic and environmental impact in the industry along with better therapeutic applications.

Author Contributions

Conceptualization: N.L.J. and A.G.; methodology: N.L.J., M.J.J., L.C.J., A.G.B., D.R.-S., A.H.-L., A.G. and C.I.S.-D.; software: D.R.-S., A.H.-L. and C.I.S.-D.; formal analysis: M.I.P., M.J.J. and L.C.J.; investigation: N.L.J., D.R.-S., J.C.-Z. and A.H.-L.; data curation: M.J.J. and A.H.-L.; writing—original draft: J.M.R., A.G.B., D.R.-S. and A.H.-L.; writing—review and editing: C.V.I., N.L.J., L.C.J., A.H.-L. and C.I.S.-D.; supervision: A.H.-L. and C.I.S.-D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Andalusian Government (project P18-RT-3786), the Spanish Science and Innovation Ministry and European Social Fund Plus (project PID2022-137603OB-I00, PID2021-127199NB-I00, PRE2022-101260), and Generalitat Valenciana (CIAICO/2022/121).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, and further inquiries can be directed to the corresponding authors.

Acknowledgments

We thank the Supercomputing Centers of UGR, CSIC, and CESGA for high-performance computing services.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cerna, J.; Morales, G.; Eyler, G.; Cañizo, A. Bulk polymerization of styrene catalyzed by biand tri-functional cyclic initiators. J. Appl. Polym. Sci. 2002, 83, 1–11. [Google Scholar] [CrossRef]
  2. Dong, Y.; McCulloug, K.J.; Wittlin, S.; Chollet, J.; Vennerstrom, J.L. The structure and antimalarial activity of dispiro-1, 2, 4, 5-tetraoxanes derived from (+)-dihydrocarvone. Bioorg. Med. Chem. Lett. 2010, 20, 6359–6361. [Google Scholar] [CrossRef]
  3. Atheaya, H.; Khan, S.I.; Mamgaina, R.; Rawata, D.S. Synthesis, thermal stability, antimalarial activity of symmetrically and asymmetrically substituted tetraoxanes. Bioorg. Med. Chem. Lett. 2008, 18, 1446–1449. [Google Scholar] [CrossRef]
  4. Creek, D.J.; Ryan, E.; Charman, W.N.; Chiu, F.C.; Prankerd, R.J.; Vennerstrom, J.L.; Charman, S.A. Stability of peroxide antimalarials in the presence of human hemoglobin. Antimicrob. Agents Chemother. 2009, 53, 3496–3500. [Google Scholar] [CrossRef]
  5. Oxley, J.C.; Smith, J.L.; Chen, H. Decomposition of a Multiperoxidic Compound: Triacetone Triperoxide (TATP). Prop. Explos. Pyrotech. 2002, 27, 209–216. [Google Scholar] [CrossRef]
  6. Antolínez, I.V.; Barbosa, L.C.A.; Maltha, C.R.A.; Pereira, G.A.M.M.; Silva, A.A. Synthesis and phytotoxic profile of a new tetraoxane designed from a commercial auxin. Quim. Nova 2020, 43, 1–8. [Google Scholar] [CrossRef]
  7. Jorge, N.L.; Romero, J.M.; Grand, A.; Hernández-Laguna, A. Gas-phase thermolysis reaction of formaldehyde diperoxide. Kinetic study and theoretical mechanisms. Chem. Phys. 2012, 39, 37–45. [Google Scholar] [CrossRef]
  8. Profeta, M.I.; Romero, J.M.; Jorge, N.L.; Grand, A.; Hernández-Laguna, A. Theoretical study of the gas-phase thermolysis of 3-methyl-1, 2, 4, 5-tetroxane. J. Mol. Model. 2014, 2, 2224–2226. [Google Scholar] [CrossRef] [PubMed]
  9. Bordón, A.G.; Pila, A.N.; Profeta, M.I.; Romero, J.M.; Jorge, L.C.; Jorge, N.L.; Sainz-Díaz, C.I.; Grand, A.; Hernández-Laguna, A. Theoretical study of the gas-phase thermolysis reaction of 3, 6-dimethyl-1, 2, 4, 5-tetroxane. Methyl and axial-equatorial substitution effects. J. Mol. Model. 2019, 25, 217. [Google Scholar] [CrossRef]
  10. Leiva, L.C.; Jorge, N.L.; Romero, J.M.; Cafferata, L.F.R.; Gómez Vara, M.E.; Castro, E.A. Thermal decomposition reaction of 3, 3, 6, 6- tetramethyl- 1, 2, 4, 5- tetroxane in 2-methoxy-ethanol solution. Chem. Heterocycl. Compd. 2010, 45, 1455–1459. [Google Scholar] [CrossRef]
  11. Farahani, P.; Roca-Sanjuán, D.; Zapata, F.; Lindh, R. Revisiting the nonadiabatic process in 1, 2-dioxetane. J. Chem. Theory Comput. 2013, 9, 5404–5411. [Google Scholar] [CrossRef] [PubMed]
  12. Schmigel, W.W.; Litt, F.A.; Cowan, D.O. A convenient gas chromatographic method for determining activation energies of first-order reactions. J. Org. Chem. 1968, 33, 3334–3345. [Google Scholar] [CrossRef]
  13. Cafferata, L.F.R.; Lombardo, J.D. Kinetics and Mechanism of the Thermal Decomposition Reaction of Acetone Cyclic Diperoxide in the Gas Phase. Int. J. Chem. Kin. 1994, 26, 503–509. [Google Scholar] [CrossRef]
  14. Richardson, W.H.; Montgomery, F.C.; Yelvington, M.B.; HO’Neal, E. Kinetics of the thermal decomposition of 3, 3-diphenyl- and 3, 3-dibenzyl-1, 2-dioxetane. Consideration of stepwise and concerted mechanisms. J. Am. Chem. Soc. 1974, 96, 7525–7532. [Google Scholar] [CrossRef]
  15. Adam, W.; Baader, W.J. Effects of methylation on the thermal stability and chemiluminescence properties of 1, 2-dioxetanes. J. Am. Chem. Soc. 1985, 107, 410–416. [Google Scholar] [CrossRef]
  16. Lefller, J.E. Parameters for the Description of Transition States. Science 1953, 117, 340–341. [Google Scholar] [CrossRef] [PubMed]
  17. Hammond, G.S. A Correlation of Reaction Rates. J. Am. Chem. Soc. 1955, 77, 334–338. [Google Scholar] [CrossRef]
  18. El-Sayed, M.A. Triplet state. Its radiative and nonradiative properties. Acc. Chem. Res. 1968, 1, 8–16. [Google Scholar] [CrossRef]
  19. Vacher, M.; Farahani, P.; Valentini, A.; Frutos, L.M.; Karlsson, H.O. How Do Methyl Groups Enhance the Triplet Chemiexcitation Yield of Dioxetane? J. Phys. Chem. Lett. 2017, 8, 3790–3794. [Google Scholar] [CrossRef]
  20. Frisch, M.J.; Truck, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision B.01; Gaussian Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  21. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  22. Slater, J.C.; Phillips, J.C. Quantum Theory of Molecules and Solids Vol. 4: The Self-Consistent Field for Molecules and Solids. Phys. Today 1974, 27, 49–50. [Google Scholar] [CrossRef]
  23. Lee, C.; Yang, L.W.; Par, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef] [PubMed]
  24. Schlegel, H.B. Optimization of Equilibrium Geometries and Transition Structures. J. Comput. Chem. 1982, 3, 214–218. [Google Scholar] [CrossRef]
  25. Frisch, A.; Nielsen, A.B.; Holder, A.J. Gauss View Molecular Visualization Program; User Manual; Gaussian Inc.: Pittsburgh, PA, USA, 2001. [Google Scholar]
  26. Gonzalez, C.; Schlegel, H.B. An improved algorithm for reaction path following. J. Chem. Phys. 1989, 90, 2154–2161. [Google Scholar] [CrossRef]
  27. Gonzalez, C.; Schlegel, H.B. Reaction Path Following in Mass-Weighted Internal Coordinates. J. Phys. Chem. 1990, 94, 5523–5529. [Google Scholar] [CrossRef]
  28. Andersson, K.; Malmqvist, P.A.; Roos, B.O.; Sadlej, A.J.; Wolinski, K. Second-order perturbation theory with a CASSCF reference function. J. Phys. Chem. 1990, 94, 5483–5488. [Google Scholar] [CrossRef]
  29. Roca-Sanjuán, D.; Aquilante, F.; Lindh, R. Multiconfiguration second-order perturbation theory approach to strong electron correlation in chemistry and photochemistry. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 585–603. [Google Scholar] [CrossRef]
  30. Aquilante, F.; Autschbach, J.; Carlson, R.K.; Chibotaru, L.F.; Delcey, M.G.; De Vico, L.; Fernández-Galván, I.; Ferre, N.; Frutos, L.M.; Gagliardi, L.; et al. MOLCAS 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table. J. Comput. Chem. 2016, 37, 506–541. [Google Scholar] [CrossRef]
  31. Forsberg, N.; Malmqvist, P.A. Multiconfiguration perturbation theory with imaginary level shift. Chem. Phys. Lett. 1997, 274, 196–204. [Google Scholar] [CrossRef]
  32. Malmqvist, P.A.; Roos, B.O. The CASSCF state interaction method. Chem. Phys. Lett. 1989, 155, 189–194. [Google Scholar] [CrossRef]
  33. Malmqvist, P.A.; Roos, B.O.; Schimmelpfennig, B. The restricted active space (RAS) state interaction approach with spin–orbit coupling. Chem. Phys. Lett. 2002, 357, 230–240. [Google Scholar] [CrossRef]
Scheme 1. Reaction paths of the singlet ground state ACDP and TMT by single-step/concerted (Branch S-I: from S-c to S-P) and stepwise (Branch S-II: from S-c to S-p) mechanisms. Arrow up and down means spin up and down, respectively.
Scheme 1. Reaction paths of the singlet ground state ACDP and TMT by single-step/concerted (Branch S-I: from S-c to S-P) and stepwise (Branch S-II: from S-c to S-p) mechanisms. Arrow up and down means spin up and down, respectively.
Molecules 29 03274 sch001
Figure 1. (a) Concerted mechanism of S-PES for TMT and ACDP, black color for ACDP, blue for axial TMT, and green for equatorial TMT, and blue arrows indicate the transition vector of TS. (b) Activation energy of TS of the concerted mechanism (TSX) as a function of the number of methyl groups in FDP. Values for FDP, MFDP, and DMT come from references [7,8,9].
Figure 1. (a) Concerted mechanism of S-PES for TMT and ACDP, black color for ACDP, blue for axial TMT, and green for equatorial TMT, and blue arrows indicate the transition vector of TS. (b) Activation energy of TS of the concerted mechanism (TSX) as a function of the number of methyl groups in FDP. Values for FDP, MFDP, and DMT come from references [7,8,9].
Molecules 29 03274 g001aMolecules 29 03274 g001b
Figure 2. First step of stepwise mechanism S-PES of TMT and ACDP. Black color for ACDP, blue for axial TMT, and green for equatorial TMT. The main non-bonding distances are shown in Å. The spin density is indicated in brackets. Blue arrows indicate the transition vector of TS.
Figure 2. First step of stepwise mechanism S-PES of TMT and ACDP. Black color for ACDP, blue for axial TMT, and green for equatorial TMT. The main non-bonding distances are shown in Å. The spin density is indicated in brackets. Blue arrows indicate the transition vector of TS.
Molecules 29 03274 g002
Figure 3. Second step of the stepwise mechanism: (a) when acetaldehyde is the first product for both axial and equatorial isomers of TMT; (b) for ACDP and for TMT when acetone is the first product for axial and equatorial isomers. The interatomic distance of the breaking bonds are shown, and the spin densities are in brackets. The values of ACDP are in black, and those of TMT are in blue for axial TMT and green for equatorial TMT. Blue arrows mean the transition vector of TS.
Figure 3. Second step of the stepwise mechanism: (a) when acetaldehyde is the first product for both axial and equatorial isomers of TMT; (b) for ACDP and for TMT when acetone is the first product for axial and equatorial isomers. The interatomic distance of the breaking bonds are shown, and the spin densities are in brackets. The values of ACDP are in black, and those of TMT are in blue for axial TMT and green for equatorial TMT. Blue arrows mean the transition vector of TS.
Molecules 29 03274 g003
Figure 4. Third step between S-b to S-p of the S-PES: (a) for TMT for axial and equatorial isomers when acetaldehyde and oxygen are the last products; (b) for ACDP and for TMT for axial and equatorial isomers when acetone and oxygen are the last products. The interatomic distance of the breaking bond is indicated in Å, and the electron densities are shown in brackets. Black color for ACDP, blue for axial TMT, and green for equatorial TMT. Blue arrows mean the transition vector of TS.
Figure 4. Third step between S-b to S-p of the S-PES: (a) for TMT for axial and equatorial isomers when acetaldehyde and oxygen are the last products; (b) for ACDP and for TMT for axial and equatorial isomers when acetone and oxygen are the last products. The interatomic distance of the breaking bond is indicated in Å, and the electron densities are shown in brackets. Black color for ACDP, blue for axial TMT, and green for equatorial TMT. Blue arrows mean the transition vector of TS.
Molecules 29 03274 g004aMolecules 29 03274 g004b
Scheme 2. Reaction paths of the triplet ground states of ACDP and TMT. The two arrows down means two parallel spins.
Scheme 2. Reaction paths of the triplet ground states of ACDP and TMT. The two arrows down means two parallel spins.
Molecules 29 03274 sch002
Figure 5. The first steps of the stepwise mechanism between T-o and T-b′ of the T-PES: (a) TMT when acetaldehyde is the first product for both axial (ax) and equatorial (eq) isomers; (b) TMT and ACDP when the acetone is the first product. Black color is for ACDP, blue for axial TMT, and green for equatorial TMT. Blue arrows indicate the transition vector of TS, the breaking bond distances are indicated in Å, and the spin density values are in brackets.
Figure 5. The first steps of the stepwise mechanism between T-o and T-b′ of the T-PES: (a) TMT when acetaldehyde is the first product for both axial (ax) and equatorial (eq) isomers; (b) TMT and ACDP when the acetone is the first product. Black color is for ACDP, blue for axial TMT, and green for equatorial TMT. Blue arrows indicate the transition vector of TS, the breaking bond distances are indicated in Å, and the spin density values are in brackets.
Molecules 29 03274 g005
Figure 6. The second step between T-b to T-p of the T-PES of (a) TMT and ACDP when the acetone and oxygen are the products; (b) TMT when the acetaldehyde and oxygen are the products. The distance of the breaking bond is indicated in Å, and the spin density values are in brackets. The TVs are indicated by blue arrows.
Figure 6. The second step between T-b to T-p of the T-PES of (a) TMT and ACDP when the acetone and oxygen are the products; (b) TMT when the acetaldehyde and oxygen are the products. The distance of the breaking bond is indicated in Å, and the spin density values are in brackets. The TVs are indicated by blue arrows.
Molecules 29 03274 g006
Figure 7. Bonding (green color) and anti-bonding (* and red color) natural orbitals used in the CASSCF/CASPT2(12-in-8)/ANO-L-VDZP computations of FDP.
Figure 7. Bonding (green color) and anti-bonding (* and red color) natural orbitals used in the CASSCF/CASPT2(12-in-8)/ANO-L-VDZP computations of FDP.
Molecules 29 03274 g007
Figure 8. Vertical energies of the four singlet and four triplet electronic states at the ground-state geometry of the diradical obtained with the CASSCF/CASPT2(12-in-8)/ANO-L-VDZP calculations. Green levels are excited triplet states (from T1 to T4), and blue levels are ground singlet state (S0) and their excited states (S1 to S3) of each derivative in increasing order of energy.
Figure 8. Vertical energies of the four singlet and four triplet electronic states at the ground-state geometry of the diradical obtained with the CASSCF/CASPT2(12-in-8)/ANO-L-VDZP calculations. Green levels are excited triplet states (from T1 to T4), and blue levels are ground singlet state (S0) and their excited states (S1 to S3) of each derivative in increasing order of energy.
Molecules 29 03274 g008
Figure 9. CASSCF (12-in-8)/ANO-L-VDZP electron spin density for the triplet states of the diradical of FDP. The corresponding densities for the singlets S0, S1, S2, and S3 are similar to those for T4, T2, T3, and T1, respectively.
Figure 9. CASSCF (12-in-8)/ANO-L-VDZP electron spin density for the triplet states of the diradical of FDP. The corresponding densities for the singlets S0, S1, S2, and S3 are similar to those for T4, T2, T3, and T1, respectively.
Molecules 29 03274 g009
Table 1. Temperature (K) and rate constant (k) (s−1) of DMT 0.001 M in solution of benzene and ACDP 0.02 M in solution of tetrahydrofuran (THF).
Table 1. Temperature (K) and rate constant (k) (s−1) of DMT 0.001 M in solution of benzene and ACDP 0.02 M in solution of tetrahydrofuran (THF).
Tkexp,DMTkexp,ACPD
4930.2450.461
5030.4190.732
5130.7041.142
5231.0721.749
5331.8672.642
5432.9593.926
Table 2. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies ν b, spin contamination S2, and relative energies (ΔE) c of critical points of the singlet-ground-state PES of the thermolysis reaction of TMT and ACDP. Concerted mechanism.
Table 2. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies ν b, spin contamination S2, and relative energies (ΔE) c of critical points of the singlet-ground-state PES of the thermolysis reaction of TMT and ACDP. Concerted mechanism.
StructureE aZPE aν bS2ΔE c
TMT d
C-6axC-6eqC-6axC-6eqC-6axC-6eqC-6axC-6eqC-6axC-6eq
c−497.124296−497.1292050.1604750.160432 0.0000.0000.000.00
TSx−497.010763−497.0119190.1526430.1479808477820.0000.00066.3368.79
P−497.162535−497.1627480.1491170.149286 0.0000.000−31.12−28.04
ACDP
E aZPE aν bS2ΔE c
c−536.4304630.188318 0.0000.00
TSx−536.3181980.1810898270.00065.91
P−536.4724920.178024 0.000−32.83
a Hartree. b Imaginary frequencies in cm−1. c Relative energies in kcal/mol, Δ E = [ E t o t a l + Z P E ] ( c ) E t o t a l + Z P E ( i ) , where i, represents any critical point of the species in the table. d In position C-6 axial (ax) and equatorial (eq) isomers.
Table 3. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the singlet-ground-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers. Acetaldehyde as the first intermediate product (Figure 3a) and acetone as the second product (Figure 4b).
Table 3. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the singlet-ground-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers. Acetaldehyde as the first intermediate product (Figure 3a) and acetone as the second product (Figure 4b).
StructureE aZPE aν bS2ΔE c
C-6axC-6eqC-6axC-6eqC-6axC-6eqC-6axC-6eqC-6axC-6eq
Steps 1–2
c−497.124296−497.1292050.1604750.160432 0.0000.0000.000.00
TSco−497.092215−497.0961800.1565380.1565932902350.0290.03017.6618.31
o−497.101741−497.1017410.1549810.154855 0.0790.07911.7813.73
TSob′−497.070352−497.0715630.1525740.1525009218800.4340.43828.89 (17.11) d31.19 (17.46) d
b′−497.114959−497.1149560.1515740.151577 0.0960.0960.27
(−11.50) d
3.38
(−10.35) d
Step 3
b−343.323732−343.3237320.0934650.093465 0.0930.0930.000.00
TSbp −343.303619−343.3036190.0904560.090456174817480.3310.33110.7310.73
p−343.370545−343.3705450.0907150.090715 0.000 0.000−31.10−31.10
a Hartrees. b Imaginary frequencies in cm−1. c Relative energies in kcal/mol, Δ E = [ E t o t a l + Z P E ] ( c / b ) [ E t o t a l + Z P E ] ( i ) where i, represents any critical point of the species in the table. d Relative energies in kcal/mol with respect critical point o.
Table 4. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the singlet-ground-state PES’ of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers. Acetone as intermediate product of step 2 (Figure 3b) and acetaldehyde as product of the last step (Figure 4a).
Table 4. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the singlet-ground-state PES’ of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers. Acetone as intermediate product of step 2 (Figure 3b) and acetaldehyde as product of the last step (Figure 4a).
StructureE aZPE aνbS2ΔE c
C-6axC-6eqC-6axC-6eqC-6axC-6eqC-6axC-6eqC-6axC-6eq
Steps 1–2
c−497.124296−497.1292050.1604750.160432 0.0000.0000.000.00
TSco−497.092215−497.0961800.1565380.1565932902350.0290.03017.6618.31
o−497.100033−497.1017410.1549810.154855 0.0790.07911.7813.73
TSob′−497.073225−497.0748720.1526870.1524998408500.4340.42927.16 (15.38) d29.11 (15.38) d
b′−497.119271−497.1188080.1520030.151632 0.0970.095−2.16
(−13.94) d
1.00
(−12.72) d
Step 3
b−304.018286−304.0188820.0650150.064871 0.0940.0990.000.00
TSbp−303.996727−303.9973980.0621880.061560389243960.3770.38911.7511.40
p−304.061264−304.0611380.0617600.061722 0.0000.000−29.01−28.64
a Hartrees. b Imaginary frequencies in cm−1. c Relative energies in kcal/mol, Δ E = [ E t o t a l + Z P E ] ( c / b ) [ E t o t a l + Z P E ] ( i ) , where i, represents any critical point of the species in the table. d Relative energies in kcal/mol with respect to critical point o.
Table 5. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the singlet-ground-state PES’ of the thermolysis reaction of 3,3,6,6-tetramethyl-1,2,4,5-tetroxane (ACDP).
Table 5. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the singlet-ground-state PES’ of the thermolysis reaction of 3,3,6,6-tetramethyl-1,2,4,5-tetroxane (ACDP).
StructureE aZPE aν bS2ΔE c
Step 1
c−536.4304630.188318 0.0000.00
TSco−536.3973110.1846532330.03118.50
o−536.4039120.183363 0.07713.55
TSob′−536.3771670.1810228460.42228.86 (15.31) d
b′−536.4246360.180440 0.096−1.29 (−14.83) d
Step 2
b−343.3237320.093465 0.0930.00
TSbp −343.3036190.09045617480.33110.73
p−343.3705450.090715 0.000−31.10
a Hartrees. b Imaginary frequencies in cm−1. c Relative energies in kcal/mol, Δ E = [ E t o t a l + Z P E ] ( c / b ) [ E t o t a l + Z P E ] ( i ) , where i, represents any critical point of the species in the table. d Relative energies in kcal/mol with respect to critical point o.
Table 6. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the triplet-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers. Acetaldehyde as the first intermediate product (Figure 5a) and acetone as the second product (Figure 6a).
Table 6. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the triplet-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers. Acetaldehyde as the first intermediate product (Figure 5a) and acetone as the second product (Figure 6a).
StructureE aZPE aν bS2ΔE c
C-6axC-6eqC-6axC-6eqC-6axC-6eqC-6ax C-6eqC-6axC-6eq
Step 1
o−497.100010497.1017060.1549870.154850 2.0092.0000.000.00
TSob′−497.070759−497.0718790.1528030.1525528888752.0002.00016.98 d15.13 d
b′−497.115218−497.1153400.1516760.151731 2.0212.000−11.62 e−10.51 e
Step 2
b−343.320980−343.3209800.0935080.093508 2.0002.0000.000.00
TSbp −343.310328−343.3103280.0906960.090696461246122.0002.0006.886.88
p−343.389853−343.3898530.0905780.090578 2.0002.000−43.10−43.10
a Hartrees. b Imaginary frequencies in cm−1. c Relative energies in kcal/mol, Δ E = [ E t o t a l + Z P E ] ( o / b ) [ E t o t a l + Z P E ] ( i ) , where i represents any critical point of the species in the table. d The relative energies with respect to Sc are 28.78 and 31.03 kcal/mol for the axial and equatorial isomers, respectively. e The relative energies with respect to Sc are −3.67 and −4.33 kcal/mol for the axial and equatorial isomers, respectively.
Table 7. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the triplet-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers. Acetone as the first intermediate product (Figure 5b) and acetaldehyde as the second product (Figure 6b).
Table 7. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of critical points of the triplet-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers. Acetone as the first intermediate product (Figure 5b) and acetaldehyde as the second product (Figure 6b).
StructureE aZPE aν bS2ΔE c
C-6axC-6eqC-6axC-6eqC-6axC-6eqC-6axC-6eqC-6axC-6eq
Step 1
o−497.100010−497.1017060.1549870.154850 2.0092.0000.000.00
TSob′−497.073526−497.0751790.1527270.1525518378482.0002.00015.20 d15.20 d
b′−497.121021−497.1166530.1519320.151264 2.0212.000−15.10 e−11.63 e
Step 2
b−304.018708−304.0198760.0650590.064864 2.0002.0000.000.00
TSbp−304.003558−303.0043280.0624160.061914378435032.0002.0007.907.84
p−304.080579−304.0806680.0617370.061720 2.0002.000−40.91−40.12
a Hartrees. b Imaginary frequencies in cm−1. c Relative energies in kcal/mol, Δ E = [ E t o t a l + Z P E ] ( o / b ) [ E t o t a l + Z P E ] ( i ) where i, represents any critical point of the species in the table. d The relative energies with respect to Sc are 27.16 and 28.95 kcal/mol for the axial and equatorial isomers, respectively. e The relative energies with respect to Sc are −2.16 and −2.12 kcal/mol for the axial and equatorial isomers, respectively.
Table 8. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of the critical points of the triplet PES for the thermolysis reaction of 3,3,6,6-tetramethyl-1,2,4,5-tetroxane (ACDP).
Table 8. Total energies (E) a, zero-point energies (ZPE) a, imaginary frequencies (ν) b, spin contamination (S2), and relative energies (ΔE) c of the critical points of the triplet PES for the thermolysis reaction of 3,3,6,6-tetramethyl-1,2,4,5-tetroxane (ACDP).
StructureE aZPE aΝ bS2ΔE c
qc−536.4017870.183464 2.0001.39 d
Step 1
o−536.4038950.183362 2.0000.00
TSob′ −536.3774830.1806908442.00015.13 e
b′−536.4251140.180509 2.000−15.12 f
Step 2
b−343.3209800.093508 2.0000.00
TSbp−343.2196310.09069646122.0006.88
p−343.3898430.090578 2.000−43.10
a Hartrees. b Imaginary frequencies in cm−1. c Relative energies in kcal/mol, Δ E = [ E t o t a l + Z P E ] ( o / b ) [ E t o t a l + Z P E ] ( i ) , where i represents any critical point of the species in the table. d With respect to T-o. e The relative energy with respect to Sc is 28.70 kcal/mol. f The relative energy with respect to Sc is −1.56 kcal/mol.
Table 9. Spectroscopic characterization of the reactant. EVA corresponds to the CASPT2 (12-in-8)/ANO-L-VDZP vertical absorption energy (eV) of FDP. Main nature of the transitions and the corresponding weight of the electronic configurations are also shown (see orbitals in Figure 7).
Table 9. Spectroscopic characterization of the reactant. EVA corresponds to the CASPT2 (12-in-8)/ANO-L-VDZP vertical absorption energy (eV) of FDP. Main nature of the transitions and the corresponding weight of the electronic configurations are also shown (see orbitals in Figure 7).
StateNatureWeight (%)EVA (eV)
S1n3σ1*785.09
S2n2σ2*745.14
S3n4σ1*805.78
T1n3σ1*393.87
n2σ2*26
T2n2σ1*343.93
n3σ2*30
T3n4σ1*544.42
n1σ2*27
T4n4σ2*454.81
n1σ1*35
Table 10. SOC complex vector norms (see Computational Details) in cm−1 for the tetraoxane series. Data obtained with the CASSCF(12,8)/ANO-L-VDZP method. Note that Ti encompasses the three ms components of the triplet state.
Table 10. SOC complex vector norms (see Computational Details) in cm−1 for the tetraoxane series. Data obtained with the CASSCF(12,8)/ANO-L-VDZP method. Note that Ti encompasses the three ms components of the triplet state.
FDPMFDP (ax)MFDP (eq)DMT (aa)DMT (ae)DMT (ee)TMT (aae)TMT (aee)ACDP
T1-S00.030.060.050.100.090.070.150.130.21
T2-S021.1265.2766.4623.7241.8127.7761.3057.9333.84
T3-S096.0476.5072.4295.2988.9194.2476.5879.2791.94
T4-S00.230.130.200.120.140.160.410.230.65
T1-S121.1362.6564.1623.6732.0827.7651.9546.5633.71
T2-S10.100.500.430.190.410.161.010.750.33
T3-S10.440.040.511.850.840.171.951.103.13
T4-S121.2862.2563.8123.9732.1526.7451.6946.4934.10
T1-S296.1275.9074.5795.5193.0094.3483.4686.6192.26
T2-S21.220.210.660.370.180.831.250.341.16
T3-S20.030.480.430.030.260.030.480.440.00
T4-S295.3975.0873.8894.3591.9793.3982.5385.6191.23
T1-S30.140.100.120.350.220.030.600.350.93
T2-S321.2664.9866.1823.9041.6627.9961.6057.7133.99
T3-S395.4772.8771.8794.5388.1993.4976.0878.6491.48
T4-S30.170.210.200.290.260.240.360.320.48
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bordón, A.G.; Profeta, M.I.; Romero, J.M.; Jorge, M.J.; Jorge, L.C.; Jorge, N.L.; Sainz-Díaz, C.I.; Cuéllar-Zuquin, J.; Roca-Sanjuán, D.; Viseras Iborra, C.; et al. Kinetic Study and Reaction Mechanism of the Gas-Phase Thermolysis Reaction of Methyl Derivatives of 1,2,4,5-Tetroxane. Molecules 2024, 29, 3274. https://doi.org/10.3390/molecules29143274

AMA Style

Bordón AG, Profeta MI, Romero JM, Jorge MJ, Jorge LC, Jorge NL, Sainz-Díaz CI, Cuéllar-Zuquin J, Roca-Sanjuán D, Viseras Iborra C, et al. Kinetic Study and Reaction Mechanism of the Gas-Phase Thermolysis Reaction of Methyl Derivatives of 1,2,4,5-Tetroxane. Molecules. 2024; 29(14):3274. https://doi.org/10.3390/molecules29143274

Chicago/Turabian Style

Bordón, Alexander G., Mariela I. Profeta, Jorge M. Romero, María J. Jorge, Lilian C. Jorge, Nelly L. Jorge, C. Ignacio Sainz-Díaz, Juliana Cuéllar-Zuquin, Daniel Roca-Sanjuán, César Viseras Iborra, and et al. 2024. "Kinetic Study and Reaction Mechanism of the Gas-Phase Thermolysis Reaction of Methyl Derivatives of 1,2,4,5-Tetroxane" Molecules 29, no. 14: 3274. https://doi.org/10.3390/molecules29143274

APA Style

Bordón, A. G., Profeta, M. I., Romero, J. M., Jorge, M. J., Jorge, L. C., Jorge, N. L., Sainz-Díaz, C. I., Cuéllar-Zuquin, J., Roca-Sanjuán, D., Viseras Iborra, C., Grand, A., & Hernández-Laguna, A. (2024). Kinetic Study and Reaction Mechanism of the Gas-Phase Thermolysis Reaction of Methyl Derivatives of 1,2,4,5-Tetroxane. Molecules, 29(14), 3274. https://doi.org/10.3390/molecules29143274

Article Metrics

Back to TopTop