Pooria Farahani*a,
Marcelo A. Oliveiraa,
Ignacio Fdez. Galvánbc and
Wilhelm J. Baader*a
aInstituto de Química, Departamento de Química Fundamental, Universidade de São Paulo, C. P. 05508-000, São Paulo, SP, Brazil. E-mail: fpooria@iq.usp.br; wjbaader@iq.usp.br
bDepartment of Chemistry – Ångström Laboratory, Theoretical Chemistry Programme, Uppsala University, P. O. Box 518, SE-751 20, Uppsala, Sweden
cUppsala Centre for Computational Chemistry, Uppsala University, P. O. Box 518, SE-751 20, Uppsala, Sweden
First published on 20th March 2017
1,2-Dioxetanones have been considered as model compounds for bioluminescence processes. The unimolecular decomposition of these prototypes leads mainly to the formation of triplet excited states whereas in the catalysed decomposition of these peroxides singlet states are formed preferentially. Notwithstanding, these cyclic peroxides are important models to understand the general principles of chemiexcitation as they can be synthesised, purified and characterised. We report here results of experimental and theoretical approaches to investigating the decomposition mechanism of spiro-adamantyl-1,2-dioxetanone. The activation parameters in the unimolecular decomposition of this derivative have been determined by isothermal kinetic measurements (30–70 °C) and the chemiluminescence activation energy calculated from the correlation of emission intensities. The activation energy for peroxide decomposition proved to be considerably lower than the chemiluminescence activation energy indicating the existence of different reaction pathways for ground and excited state formation. These experimental results are compared with the calculations at the complete active space second-order perturbation theory (CASPT2), which reveal a two-step biradical mechanism starting by weak peroxide bond breakage followed by carbon–carbon elongation. The theoretical findings also indicate different transition state energies on the excited and ground state surfaces during the C–C bond cleavage in agreement with the experimental activation parameters.
unimolecular, as well as catalysed decomposition, has been considered by several research groups.1–5 More recently, we reported experimental data on the catalysed decomposition of two new 1,2-dioxetanone derivatives, spiro-adamantyl-1,2-dioxetanone (4) and spiro-cyclopentyl-1,2-dioxetanone (5).6
The unimolecular decomposition of 1,2-dioxetanes and 1,2-dioxetanones leads, mainly, to the formation of triplet excited
carbonyl products, with quantum yields of up to 30%, whereas the singlet excited state formation is inefficient, reaching rarely up to 1%.1–4 Peroxide derivatives with higher degree of substituents are more stable and possess higher quantum yields.2,4,7
Two extreme pathways were suggested for this transformation; a concerted mechanism, first proposed by Turro,8 as well as a biradical path, initially formulated by Richardson.9,10 The concerted mechanism is a single-step reaction process in which the C–C and O–O bonds dissociate simultaneously. The two-step biradical mechanism, however, implies that once the O–O bond is cleaved, the system enters a so-called entropic trapping region11,12 of biradical nature where four singlet and triplet states are degenerated, then, the C–C rupture comes into action. Some of the obtained experimental data appear to be more compatible with the biradical mechanism. Although, the others can be better rationalised with a concerted transformation. Consequently, a “merged” reaction pathway, in which the O–O elongation is not fully completed when the C–C cleavage occurs, has been suggested maintaining the biradical nature of the process without formation of a biradical intermediate in a potential energy surface minimum.7 This mechanistic proposal appears to accommodate all (or at least most) of the available experimental data up to now. Although, generally no charge or electron transfer processes are considered in the unimolecular decomposition of cyclic peroxides, some recent work has pointed out a feasible importance of charge and electron spin transfer actions.13
During the last 10 years, significant efforts have been performed to understand the CL phenomenon with computational methodologies, by studying the decomposition mechanism of 1 and 2, in detail.12,14,15 These ab initio findings, as well, confirm a step-wise biradical mechanism for 1, taking into account the time it takes to pass the entropic trapping region before complete thermal decomposition. According to these efforts, in contrast to singlet excited state formation, triplet excitation requires almost the same energy as the thermal decomposition to ground state products. However, the dissociation mechanism of 2, is not supported yet neither by ab initio15 nor DFT studies, due to the lack of molecular dynamics simulations.16
An open mechanistic question in the unimolecular 4-membered ring peroxide decomposition is if there are distinct reaction pathways for the formation of ground and excited state products, which possess different activation energies? Only a few experimental studies have addressed this question; Turro and co-workers have investigated the possibility of different transition states for the formation of ground state, triplet and singlet excited state products, measuring the isothermal and CL activation energies for the decomposition of tetramethyl-1,2-dioxetane and concluded that ground as well as excited state reaction occur all by a common transition state.17 In a study on the unimolecular decomposition of dimethyl-1,2-dioxetanone (3), Schmidt and Schuster have reported different activation energies for the formation of ground and excited state products.18 However, these results have been questioned by Turro and Chow, which report similar values for the ground and excited state activation parameters in the decomposition of this 1,2-dioxetanone derivative.19
On the basis of the previous theoretical efforts with 1,2-dioxetane12,14,20 and 1,2-dioxetanone,15 a two-step biradical mechanism in which the C1–C2 rupture occurs after terminating the O1–O2 fragmentation, can be expected for the decomposition of spiro-adamantyl-1,2-dioxetanone to adamantanone and CO2 (see Fig. 1). The first step in the biradical mechanism is the O1–O2 cleavage (TSS0), after that the system enters a biradical region in which four singlet and four triplet manifolds are degenerated. The next step of the reaction mechanism is the C1–C2 bond breakage (TSS1,T1) which basically occurs along the lowest-lying excited states, since the TS corresponding the C1–C2 rupture is quite shallow on the ground state and hence difficult to compute.
The aim of the present article is to report a qualitative and quantitative study on the electronically excited state formation of spiro-adamantyl-1,2-dioxetanone as a prototype. The spiro-adamantyl-1,2-dioxetanone was chosen as the model peroxide derivative for the experimental studies, as this has been shown before to be relatively stable and hence, can be conveniently purified by low-temperature recrystallization. The chemiexcitation quantum yields and the activation parameters of the unimolecular decomposition of this 1,2-dioxetanone derivative were determined and the results obtained compared to those extracted from theoretical calculations applying multiconfigurational strategies. Especially, the possible existence of different reaction pathways for ground and excited state formation was addressed by experimental and theoretical approaches. It is noteworthy that a theoretical study on a similar but larger system, 3-(2-spiro-adamantyl)-4-methoxy-4-(3-phosphoryloxy)-phenyl-1,2-dioxetane has been recently reported proposing a reversible charge transfer-initiated luminescence mechanism.21
For the determination of the activation parameters for the unimolecular decomposition of 4, the decomposition kinetics in toluene were followed by the intensity of the direct CL emission, corresponding to 2-adamantanone fluorescence, in a temperature range of around 30 to 70 °C (Fig. 2). The kinetic curves give rise to the observed rate constants (kobs) and the initial emission intensities (I0), obtained by the first-order exponential fitting. The correlation of the decomposition rate constants with the temperature, according to Eyring plots, showed to be linear and allows the determination of the activation enthalpies (Δ≠H) and entropies (Δ≠S) (Fig. 3, Table 1). Using the I0 values, instead of kobs, in the Eyring plot, the activation enthalpy for the CL process (Δ≠HCL) can be obtained. This approach is based on the fact that the emission intensities correspond directly to a reaction rate, as these intensities depend on the decomposition rate constant, the 1,2-dioxetanone concentration and the emission quantum yield. Therefore, the emission intensities are proportional to the rate constant of the decomposition process leading to excited state formation and can be utilized to obtain the activation enthalpy of this process. However, as no rate constants are utilized for the correlation, it is not possible to obtain values for the activation entropy (Fig. 3, Table 1). The reported activation parameters are obtained from at least three independent rate constant measurements at each temperature. The error bars come from linear regression of rate constant values according to the Eyring equation, using error progression analysis for the calculation of the deviation in the Δ≠S and Δ≠G values.
Fig. 2 Kinetic emission intensity curves for the direct emission in the unimolecular decomposition of spiro-adamantyl-1,2-dioxetanone (4) (1.33 mmol L−1) at different temperatures. |
Δ≠H (kcal mol−1) | Δ≠S (cal mol−1 K−1) | Δ≠G (25 °C) (kcal mol−1) | Δ≠HCL (kcal mol−1) | |
---|---|---|---|---|
Unimol. decomp. | 21.7 ± 0.6 | −1.1 ± 0.8 | 22.0 ± 0.5 | 28 ± 2 |
The activation data obtained clearly indicate that the CL activation enthalpy, obtained from the correlation of the initial emission intensities, is significantly higher than the “normal” isothermal activation enthalpy, obtained from the correlation of the observed decay rate constants. Interestingly, the activation entropy is close to zero, in agreement with the occurrence of a unimolecular process, thereby excluding significant participation of dark peroxide decomposition (Table 1).
In the present study, MS-CASPT2(10in8)/ANO-VDZP geometry optimizations were carried out and the most relevant geometrical parameters for the CL reaction of 4 are compiled in Table 2. In contrast to 1,2-dioxetane14 and other alkyl-substituted 1,2-dioxetanes,25 where the O1–C1–C2–O2 dihedral angles are between 19 to 22 degrees, the starting geometry (Reac) has zero degree of the dihedral angle, whereas the C1–C2 and O1–O2 are more or less the same. The has the similar dihedral angle (6.7°) which reveals a stretching of O1–O2 bond instead of a torsion of the O1–C1–C2–O2 dihedral angle for 1,2-dioxetanes that was stated to be responsible for the formation of an entropic trapping process and hence producing excited state species.12 This difference is due to the conjugation effects caused by the C2O3, that restricts the twisting around the C1–C2 bond. It is noteworthy that the found TS dihedral angle at the CASSCF method is close to zero (0.77°). These findings are consistent with the study of 1,2-dioxetanone by Liu et al. that characterized the optimized geometries to be effectively planner.15
Reac | ||||
---|---|---|---|---|
Bond lengths (Å) | ||||
C1–C2 | 1.506 | 1.549 | 1.654 | 1.876 |
O1–O2 | 1.531 | 2.482 | 2.825 | 2.559 |
C2–O3 | 1.200 | 1.221 | 1.186 | 1.181 |
C1–O1 | 1.475 | 1.383 | 1.361 | 1.360 |
C2–O2 | 1.392 | 1.336 | 1.296 | 1.224 |
Angles (degree) | ||||
C1–C2–O3 | 139.4 | 128.8 | 110.6 | 116.2 |
O1–C1–C2 | 87.5 | 108.3 | 105.9 | 107.0 |
O2–C2–O3 | 127.1 | 119.5 | 127.9 | 140.4 |
Dihedral angle (degree) | ||||
O1–C1–C2–O2 | 0.00 | 6.76 | 70.1 | 1.82 |
Table 3 also compiles the relative energies of the optimized stationary points at the CASPT2 level of increasing accuracy as well as ZPVE and Gibbs free energy correction. The comparison of the relative energies showed that the CASSCF energies are affected by the dynamical correlation. For instance, in the structure the CASSCF underestimates the energy barrier by 12.8 kcal mol−1, this is caused by the lack of dynamical correlation. This huge difference between CASSCF and CASPT2 energies convinced us to perform CASPT2 geometry optimization on the transition structures, in order to provide quantitatively accurate activation energies. However, the final results revealed that CASSCF method is accurate enough to optimize geometries of spiro-adamantyl-1,2-dioxetanone and it, therefore, was used for the IRC calculations. Unless stated otherwise, the energy values shown in this study are the MS-CASPT2(16in13)/ANO-L-VTZP.
The ground state dissociation begins by surmounting a barrier height of 31 kcal mol−1 (the rate-limiting step). After that, 23.9 kcal mol−1 is required to cleave the C–C bond and lead to the products. Taking ZPVE corrections into account, however, decreases these values to 30.4 and 22.1 kcal mol−1, respectively, due to bond dissociation. The bond dissociation causes a change in entropy of the reaction as well. The computed CL activation energy () is 28.7 kcal mol−1, which is more or less the same as the one obtained experimentally by kinetic studies (∼28 kcal mol−1). These theoretical findings, both qualitatively and quantitatively support the manifestation of two different pathways; one with lower energy () leading to the products, and one with higher energy () producing chemiexcitation products.
Based on the MEP, the mechanism for the thermal decomposition can be described. First, stretching the O1–O2 bond leads the molecule from the reactant to formation of the . After that, the C1–C2 stretch begins and decomposes molecule to the products. At the of the MEP, the MS-CASPT2 energy gap between S0, S1 and T1 manifolds are less than 4 kcal mol−1, implying that there is a crossing where the adiabatic states strongly interact. The energetics at the reveals that the T1 is slightly higher in energy than the S0 (0.5 kcal mol−1), which indicates the presence of an intersystem crossing (ISC) on top of the internal conical intersection (CIx) between S0 and S1. This is proven by the large SOC values obtained at MS-CASPT2(16in13) level, see ESI Table S5.† From these crossing points adiabatic transformations take place on the MEPs from 1(σ,σ*) to 1,3(n,σ*). This transformation does not break the planer symmetry and after the crossing point, the 1,3(n,σ*) states more or less maintain the planer structure. However, once the molecule leaves the entropic trapping region, it twists around the O–C–C–O dihedral angle to form the which breaks the plane. This TS is the responsible for the C–C dissociation along the S0 manifold.
The MS-CASPT2//SA-CASSCF protocol did not significantly change the shape and relative positions of the ground and excited state MEPs optimized at the SA-CASSCF level. Although the SA-CASSCF(10in8) approach is not the most accurate method to describe this reaction mechanism in detail, it is quite safe to qualitatively interpret its reaction mechanism and to some extent quantitatively.
1,2-Dioxetanone | Δ≠H (kcal mol−1) | Δ≠S (cal mol−1 K−1) | Δ≠Ga (kcal mol−1) | Δ≠HCL (kcal mol−1) |
---|---|---|---|---|
a At 25 °C. | ||||
Unimolecular decomposition | ||||
4this work | 21.7 ± 0.6 | −1.1 ± 0.8 | 22.0 ± 0.5 | 28 ± 2 |
3 (ref. 3) | 19.2 ± 0.2 | −8.2 ± 0.5 | 21.6 ± 0.3 | — |
3 (ref. 4) | 21.7 ± 0.3 | 0 ± 1 | 21.7 | 25.0 ± 0.1 |
3 (ref. 5) | 21.7 ± 0.6 | — | — | 21.0 ± 0.8 |
6 (ref. 6) | 18.8 | −8.9 | 21.5 | — |
7 (ref. 5) | 17.5 ± 0.1 | −12.0 ± 0.2 | 21.1 ± 0.2 | 16.6 ± 1.3 |
10 (ref. 3) | 13.7 ± 0.3 | −24.1 ± 1.0 | 20.9 ± 0.5 | — |
9 (ref. 5) | 19.3 ± 1.4 | −4.85 ± 0.52 | 20.8 ± 1.4 | 16.6 ± 2.0 |
8 (ref. 7) | 11.2 ± 0.4 | −23.2 ± 0.8 | 18.2 ± 0.6 | — |
The activation entropies for the unimolecular decomposition of 1,2-dioxetanones vary from very negative (e.g., 10) to near-zero values (e.g., 4) (Table 4). For a unimolecular process, where two molecules of product are generated from one reactant molecule, Δ≠S values should be around zero or positive, due to a higher degree of disorder in the TS. Negative Δ≠S values for the unimolecular thermolysis of 1,2-dioxetanes have been interpreted as related to impurities-induced dark decomposition.1,2,27 This rationale may be applied to explain the considerably low values for Δ≠H of the derivatives 7 and 10, for whom highly negative activation entropy values were reported (Table 4); even so, the Δ≠G values of these compounds fall within the expected range for the series.
For the chemiexcitation reaction the Δ≠HCL values are in some cases similar to the reaction enthalpies Δ≠H (for the 1,2-dioxetanones 7 and 9) or significantly higher, as seen for peroxides 3 and 4 (Table 4). This Δ≠HCL > Δ≠H observation may be associated with the existence of an impurity-catalysed dark decomposition (with no light emission), which possess a rate constant kimp embedded to the experimentally determined kobs value, a reaction path possibly with lower activation energy. This may be the case for the unimolecular decomposition of 1,2-dioxetanone 3, which shows highly negative Δ≠S values. However, the observed Δ≠S ≈ 0 for the spiro-adamantyl derivative 4 (Table 4) indicates that this supposed impurity-catalysed process is absent or non-significant, showing that the unimolecular decomposition of this peroxide occurred in a “cleaner” system. Therefore, this 5.8 kcal mol−1 difference between Δ≠HCL and Δ≠H for 4 may point towards the existence of two distinct pathways for unimolecular peroxide decomposition: one with a lower energy, leading to ground state products, and the other with a higher energy barrier producing chemiexcited products. This is in agreement with earlier observations made by Schmidt and Schuster,28 who found a higher activation energy for the CL pathway in the unimolecular decomposition of 4,4-dimethyl-1,2-dioxetanone (3). However, in a subsequent work, Turro and Chow determined very similar activation parameters for the isothermal and CL pathway in the decomposition of the same cyclic peroxide derivative.29
These contradictive results indicate the high inherent difficulties to work with this kind of highly unstable peroxide and to obtain reproducible and solid experimental data. The results reported here are more trustable as they were obtained with a derivative which is more stable and could be purified by low-temperature recrystallization and characterized by NMR spectroscopy. Furthermore, the low activation entropy measured for the unimolecular decomposition indicates that the reaction system does not contain any impurities like heavy metal atoms, which could induce dark catalysis. Therefore, we believe to have obtained experimental evidence, for the first time, that the two concurrent reaction pathways may exist in the unimolecular decomposition of spiro-adamantyl-1,2-dioxetanone, one leading to excited state formation and the other giving directly to ground state products. This important experimental observation led us to study the process also by theoretical means.
The theoretical evidences indicate that the CL mechanism of spiro-adamantyl-1,2-dioxetanone occurs as follows. First, an activation energy of 30.4 kcal mol−1 corresponding to the O1–O2 cleavage is essential to enter the biradical region in which energy degeneracy takes place between lowest-lying four singlet and triplet states of the system. In this region, the singlet excited state is populated. Taking the significant amount of SOC into account, the population can also be transferred to the triplet excited state of the molecule. After that, a second activation energy related to the C1–C2 rupture (TSC–C) is required to dissociate the molecule. The essential energy values to surmount this barrier height for the ground and first singlet excited states are different. Our CASPT2 calculations clearly support a lower activation energy on the ground state (22.1 kcal mol−1) than the singlet manifold (27 kcal mol−1).
In summary, the theoretical calculations indicate that: (i) there is an extended biradical region in the decomposition of 4 on the ground state PES, where S0, S1 and T1 states are degenerated, allowing the formation of singlet and triplet excited states; (ii) the transition state for the rate limiting O–O bond cleavage is the same for S0, S1 and T1; (iii) on the ground state energy surface, an additional activation energy exists for C–C bond cleavage, therefore the decomposition process of 4 should be nominated as a biradical mechanism; (iv) an activation energy for C–C bond cleavage exists on the S1 and T1 excited state surfaces; indicating that singlet quantum yields (and also triplet quantum yields) should change with the temperature; (v) the S1 PES is higher in energy than the T1 PES, indication preferential formation of triplet state products.
These facts should now be related to the obtained experimental result: (a) formation of singlet as well as triplet excited states is observed in the thermal decomposition of 4, in agreement with (i); (b) the different activation energies for thermal decomposition and the chemiexcitation pathway would indicate different transition states for ground and excited state formation, contrarily to the theoretical results outlined in (ii), however, the observation of a higher energy PES for S1 formation during C–C bond cleavage, shown in (iv), can explain the experimental results, as higher quantum yields should be expected at higher temperatures, leading to the measurement of a higher CL activation energy; (c) the difference in the experimental and calculated activation energies for thermal decomposition might be due to the fact that the former were obtained in toluene and the later calculated in vacuo;30 (d) the experimentally observed (ΦT/ΦS) ratio of about 300 can be understood by the fact that the calculated triplet PES for C–C bond cleavage is lower in energy than the corresponding excited singlet PES.
Therefore, we can conclude that the theoretical calculations performed on this system are in excellent qualitative agreement with the experimental results obtained and also contribute to a more precise description of the decomposition mechanism of the model 1,2-dioxetanone derivative 4.
For the determination of chemiexcitation quantum yields, the decomposition of the peroxide 4 was performed in the presence of different concentrations of the energy acceptor ([PPO] = 5.0–50 mmol L−1, ([DBA] = 0.50–5.0 mml L−1) and the CL quantum yields determined at each energy acceptor concentration from the area under the emission intensity curves, obtained in absolute einstein (E) units by luminol calibration.23 The quantum yields at infinite energy acceptor concentration (Φ∞S and Φ∞T), which correspond to the yield of singlet and triplet excited carbonyl compounds initially formed in peroxide decomposition, were determined from double reciprocal plots between the respective quantum yields and the energy acceptor concentration (1/ΦS vs. [PPO] and 1/ΦT vs. [DBA]) as described before.7
In order to analyse state-interaction effects of the CASPT2 wave functions, the multi-state (MS) approach36 of the CASPT2 method was also used. The MS- and SS-CASPT2 approaches give rise to energy values with differences of lower than 1 kcal mol−1. This is caused by the small off-diagonal elements of the effective Hamiltonian matrix in the MS calculations. Throughout the CASPT2 calculations, freezing the core orbitals of non-hydrogen atoms as well as a modification of the zeroth-order Hamiltonian, excluding the standard ionization potential electron affinity (IPEA) (with 0.0 value) were applied.37 On the other hand, an imaginary level shift of 0.2 a.u. was used, in order to solve weakly interacting intruder-states problems.38
The selection of a reasonable active space in the both CASSCF and CASPT2 calculations was done based on the described two-step mechanism. The required amount of electrons and orbitals expected for the CL process of spiro-adamantyl-1,2-dioxetanone are 16 electrons distributed in 13 orbitals, similar active space to the study of parent 1,2-dioxetanone applying CASSCF method,15 proved by 2-RDM method.39 This selection corresponds to C1–C2, O1–O2, C1–O1, C2–O2 σ bonding and σ* antibonding, C2O3 π bonding and π* antibonding as well as the three oxygen lone pairs orbitals (see Fig. 6). This active space was employed for computing the final energies of the stationary points at MS-CASPT2 approach (MS-CASPT2(16in13)). However, the geometry optimizations applying this active space, would be time- and CPU-demanding and therefore unaffordable, with respect to size of the system. On the other hand, not all of the above orbitals are significantly important throughout the reaction mechanism. For instance, for the first step of the reaction, the six orbitals corresponding to C1–O2, C2–O2 σ bonding and σ* antibonding as well as the O3 lone pair are doubly occupied and still intact. Hence, for the geometry optimizations and the MEP calculations, it is safe to exclude these orbitals from the active space. The accuracy of this approach was consequently confirmed by performing two test calculations using CAS(10in8) and CAS(16in13) in both reactant and .
Four equally weighted roots were considered in the state average (SA)-CASSCF procedure for all the CASSCF computations. The final reported MEPs in the present study correspond to SA-CASSCF computations with the lowest-lying four singlet (S0, S1, S2 and S3) and four triplet (T1, T2, T3, T4) states.
Zero-point vibrational energy (ZPVE) and Gibbs free energy corrections were computed numerically at the SA(4)-CASSCF(10in8)/ANO-L-VDZP level of theory. Spin–orbit coupling (SOC) terms between singlet and triplet manifolds were calculated within the AMFI and CASSI frameworks40,41 at the MS-CASPT2(16in13)/ANO-L-VTZP wave function over four singlet and four triplet states. A development version of the MOLCAS-8 quantum chemistry package suite was employed for all the calculations.42
Consequently, the experimentally obtained data can be rationalized with the theoretical calculations, as these predict the temperature dependence of the quantum yields, with higher quantum yields expected for higher temperatures. This situation is equivalent to the observation of a higher activation enthalpy for the CL reaction, obtained from the temperature correlation of emission intensities.
Finally, the theoretical calculations indicate the occurrence of a biradical mechanism in the decomposition of 1,2-dioxetanone derivative 4 and explain the preferential formation of triplet excited carbonyl product.
Footnote |
† Electronic supplementary information (ESI) available: Cartesian coordinates for the involved structures and their corresponding relative electronic energies, as well as, dissociation along the lowest-lying triplet manifold. See DOI: 10.1039/c6ra26575h |
This journal is © The Royal Society of Chemistry 2017 |