Density Functional Theory Study on the Decay of Fullerenyl Radicals

Article pubs.acs.org/JPCA

Density Functional Theory Study on the Decay of Fullerenyl Radicals RC60•, ROC60•, and ROOC60• (R = tert-Butyl and Cumyl) and Polarizability of the Formed Fullerene Dimers Denis Sh. Sabirov,* Ralia R. Garipova, and Ramil G. Bulgakov Institute of Petrochemistry and Catalysis, Russian Academy of Sciences, 450075 Ufa, Russia ABSTRACT: Currently, there are no comparative studies on the possible routes of decay of fullerenyl radicals. Such information is required for development of new synthetic approaches to the fullerene derivatives, obtained via radical reactions, and understanding mechanisms of oxidative destruction of fullerenecontaining materials. In the present work, we have performed a theoretical study on the possible reactions of the selected fullerenyl radicals, generated in the wellknown experimental systems. We consider three main routes for the radical decay: formation of fullerene bisadducts via XC60• + X• reactions, dimerization of XC60•, and their interaction with molecular oxygen. The calculated heat effects explain the experimental regularities of fullerene radical reactions (e.g., the reversibility of the dimerization and low reactivity of fullerenyls toward molecular oxygen). We have found that the heat effects are maximal for the dimerization of ROOC60•, so the interaction of C60 with peroxy radicals may provide more effective synthetic routes to the single-bonded C60 dimers compared to the known ones. Moreover, the exothermicity of ROO• addition to C60 is increased during the subsequent addition that allows reconsidering the reasons underlying C60 antioxidant activity. Additionally, polarizability and its nonadditivity of fullerene dimers XC60−C60X have been found. Guided by the calculated polarizabilities, we assume a new bistable system “fullerenyl ↔ dimer” for molecular machinery. Depending on the conditions, such a system can exist in two states differing by their response to external electric fields.

■

radicals and oxygen. Generally, each of R•, RO•, or ROO• are able to attach to fullerenes (though with different rate constants),2,18−22 so it is experimentally difficult to study these radical processes separately. Experimental investigations23−25 have shown that formation of fullerenyl radicals XC60• is the first step of the C60 radical reactions:

INTRODUCTION Since the discovery of fullerenyl radicals by Krusic et al. in 1991,1 radical chemistry of fullerenes has been rapidly developed (see the recent exhaustive review2 and several key works in this field3−9). Being an electron-deficient molecule with a strong system of the conjugated double bonds, the C60 fullerene readily attaches diverse usual (nonfullerene) radicals, including C-, Si-, O-, S-, P- and metal-centered species. Numerous reaction sites in the C60 molecule cause its high effectiveness as a radical scavenger (e.g., it is able to capture up to 34 methyls1 or 16 perfluoroalkyl radicals9). Therefore, the C60 fullerene is metaphorically called a radical sponge.2 Moreover, radical reactions of fullerenes underlie syntheses of novel functional derivatives C60Xn.2,10 Type, number, and positional relationship of the radicals, tethered to the C60 core, open wide opportunities for production of fullerene derivatives with the desirable physicochemical properties. Another important issue of fullerene radical chemistry deals with ability of C60 and its derivatives to inhibit radical processes in chemical and biochemical systems.2,6 For example, effects of fullerene on the liquid-phase oxidation of hydrocarbons11−13 and lipids,14 radical polymerization (see ref 15 and references therein), and thermo-oxidative destruction of the fullerenecontaining polymers16,17 have been previously studied. Different radical particles can act as intermediates of the processes listed above: these are alkyl R•, alkoxy RO•, and peroxy radicals ROO•. The last one is a result of the reaction between alkyl © 2013 American Chemical Society

C60 + X• → XC60•

(1)

Due to delocalization of the unpaired electron on the carbon framework, addition of the next X• to any of 59 carbon atoms, remaining available in XC60•, is possible. However, the maximal spin density is concentrated on carbon atoms in positions 2 and 4, relative to the site of the first X• addition.23−26 It leads to 1,2and 1,4-bisadducts, the simplest molecular products of the fullerene radical reactions. Their formation has been experimentally detected.2 The bisadducts are formed according to the following equation: XC60• + X• → C60X 2

(2)

•

Addition of the next X to fullerenyl radicals is not the only way of XC60• decay. Their reaction with oxygen presents another alternative. As is known, fullerene is not oxidized by Received: October 3, 2013 Revised: November 21, 2013 Published: November 22, 2013 13176

dx.doi.org/10.1021/jp409845q | J. Phys. Chem. A 2013, 117, 13176−13183

The Journal of Physical Chemistry A

Article

oxygen (in the absence of light irradiation or a catalyst27); i.e., the fullerene core is rather stable toward O2. Usually, C60 reactivity is changed essentially when the unpaired electron appears on the fullerene cage28 but not in the case of fullerenyl radicals and especially with respect to oxygen. As shown in nonnumerous works,29,30 the reactivity of fullerenyls toward oxygen is surprisingly low. Indeed, ESR signals of (tBuO)nC60• radicals decay slowly in the presence of oxygen due to the reaction leading to O-centered fullerene-based radicals:29 XC60• + O2 → XC60OO•

In the present theoretical work, we have studied three main routes of decay for RC60•, ROC60•, and ROOC60• fullerenyls (R = tBu• and Ph(CH3)2C•). This set of radicals has been chosen due to the following reasons. Radicals tBuC60•, tBuOC60•, and t BuOOC60• are intermediates of the reactions underlying the modern synthetic strategies for production of oxygen-rich fullerene derivatives.21−23,40 The analogous cumyl-containing fullerenyl radicals (R = Ph(CH3)2C) are generated in the model systems of the liquid-phase oxidation of hydrocarbons.12,13,18,19 In addition, various fullerenyls ROC60• and ROOC60• are generated in different reactions in fullerene-containing systems, such as photolysis of peroxides ROOR20 and dialkoxy disulfides ROSSOR22 in the presence of C60 or catalytic fullerene alkoxylation.41

(3)

Our previous ESR study on the interaction between EtC60• and O2 has shown that the ESR signal intensity is not changed in air atmosphere (compared to its intensity in argon) and only a broadening of the ESR line is observed. It means that the complex EtC60•···O2 is formed whereas no chemical reaction between EtC60• and O2 takes place. As is known,31 hydrocarbon radicals facilely react with oxygen (e.g., in the case of alkyl radicals, the rate constants are ∼109 L/(mol s)). In this aspect, low reactivity of fullerenyl radicals toward O2 seems unusual. Dimerization of fullerenyl radicals is the third route of their decay. The reversible formation of the single-bonded dimers is typical for alkyl- and fluoroalkylfullerenyl radicals:2,4,9,32,33 30

XC60• ↔ XC60 − C60X

■

COMPUTATIONAL DETAILS The potential energy surface scanning and all optimizations have been performed by density functional theory method PBE/3ζ42,43 implemented in the Priroda program.44 The 3ζ basis set describes electronic configurations of molecular systems by the orbital basis sets of contracted Gaussian-type functions (5s,1p)/[3s,1p] for H, (11s,6p,2d)/[6s,3p,2d] for C and O, which have been used in combination with the densityfitting basis sets of uncontracted Gaussian-type functions (5s,2p) for H, (10s,3p,3d,1f) for C and O atoms. The PBE/ 3ζ method reproduces structures and physicochemical characteristics of fullerenes and their derivatives with high accuracy, as shown in numerous theoretical works.26,35,37,38,45−49 In addition, it correctly reproduces the experimental data on the distribution of spin density in hydro-, fluoro-, and methylfullerenyl radicals.26 After DFT optimizations and vibration modes solving (to prove that all the stationary points, respective to the radicals/ molecules under study, are minima of the potential energy surfaces) by standard techniques, the heat effects of the reactions have been calculated as the differences between the total energies E of the products and the reactants with inclusion of zero-point vibrational energy corrections εZPV and the temperature corrections Hcorr (T = 298 K):

(4)

Dimers XC60−C60X with X = R have been studied experimentally but the formation of the analogous structures with RO or ROO moieties attached has not been discussed. The interest in such compounds is caused by the opportunity of their use for generation of free radicals in fullerenecontaining systems.34 In addition, theoretical studies predict that C60 dimers should demonstrate outstandingly high values of the first35 and the second36 polarizabilities. It makes these chemicals prospective for optical applications and nanodevices. Thus, there are at least three alternative routes of fullerenyl radical decay. For simplification, we have excluded reactions of cross-addition, when different radicals are attached to one fullerene core, and reactions of fullerenyl radicals with the solvent (this is close to the experimental fullerene-containing systems with nonreactive medium, e.g., benzene11). We can propose that formation of 1,2- and 1,4-bisadducts in “C60 + X•” systems should be thermodynamically more favorable than the dimerization due to the observed reversibility of the last one. However, the comparative thermodynamic advantages of the aforementioned reactions have not been studied both theoretically and experimentally. At the same time, such data on the reactivity of fullerenyl radicals may be important to developing the effective synthetic methods, directed to the obtainment of only one desired product (only XC60−C60X or only C60X2). Thus, thermochemical parameters of these reactions can be applied to a study on the mechanism of oxidative destruction of the fullerene-containing materials. Reactions of radical addition usually have low activation barriers, so they are successfully described by heat effects, facilely calculated in terms of the modern DFT methods. Such an approach to estimation of fullerene reactivity has been widespread.37−39 For example, a DFT study on H•, F•, and CH3• addition to the diverse fullerenes (from C20 to C84) allowed us to obtain a linear correlation between the heat effects of the addition and the curvature indices of reaction sites in the fullerene molecules.26

ΔHr° =

∑

(Etot + εZPV + Hcorr)

products

−

∑

(Etot + εZPV + Hcorr)

reactants

(5)

Components of the polarizability tensors have been calculated in terms of the finite field approach as the secondorder derivatives of the total energy E with respect to the homogeneous external electric field. They have been calculated in the arbitrary coordinate system and then diagonalized. Eigenvalues of polarizability tensors αxx, αyy, αzz have been used for the calculation of the mean polarizabilities of the molecules: α=

1 (αxx + αyy + αzz) 3

(6)

The PBE/3ζ-based finite field approach allows reproducing the measured mean polarizabilitites of C60 and C70 fullerenes.35,49−51 The mathematical operations on polarizability tensors have been performed by POLARIZ program,52 worked out by us for the fast data processing upon calculations of polarizability. 13177



■

Article

Table 2. Spin Densities of Different Reaction Sites in XC60• Fullerenyls (Mulliken Population Analysis)

RESULTS AND DISCUSSION Formation of fullerenyl radicals XC60• at reaction 1 is exothermic, regardless of the type of X• (Table 1). The formed

X

Table 1. Heat Effects of R•, RO•, and ROO• Addition to C60 (Reaction 1, kJ/mol) R X• •

R RO• ROO•·

t

Bu

Ph(CH3)2C

−64.4 −59.4 −8.4

−11.9 −59.0 −13.0

position

R

RO

ROO

2 4 16

0.444 (0.444)a 0.211 (0.210) 0.090 (0.090)

0.430 (0.432) 0.220 (0.217) 0.087 (0.086)

0.419 (0.423) 0.214 (0.216) 0.088 (0.088)

a

Values for R = tBu are given; values for R = Ph(CH3)2C are given in parentheses.

equiprobability. For example, the small difference (less than 10 kJ/mol) in the highly exothermic heat effects of 1,2- and 1,4addition of tBuOO•, predicted here theoretically, well corresponds to the experimental identification of both 1,4C60(OOtBu)2 and 1,2-C60(OOtBu)2 among the products of interaction between C60 and tert-butylperoxy radicals.40 The heat effects of the second X• addition to the C60 core (reaction 2) are significantly higher than the respective values for the first step (reaction 1). Thus, in the case of 1,4-C60tBu2, 1,4-C60(OtBu)2, and 1,4-C60(OOtBu)2, the enthalpies of reactions 2 equal to −104.9, −144.1, and −102.5 kJ/mol whereas formation of the respective tBuC60•, tBuOC60•, and t BuOOC60• radicals are accompanied by the lower heat effects, being −64.4, −59.4, and −8.4 kJ/mol, respectively (Table 1). It is explained by the fact that decay of one radical occurs on the first step whereas on the next one, two radicals (X• and XC60•) cease to exist. Obviously, all the subsequent odd steps are characterized by the lower heat effects than the even ones. The latest statement has been demonstrated for the case of the subsequent addition of tert-butylperoxy radicals to C60, a wellknown experimental system, studied by Gan et al.40 (Figure 2). Here, the heat effects of the odd steps describe reactivity of C60 (n = 1) and its adducts C60(OOtBu)2 (n = 3) and C60(OOtBu)4 (n = 5) toward tBuOO• radicals. Figure 2 allows elucidating a trend to the enhancement of the reactivity during the addition of peroxy radicals to C60. Indeed, the exothermicity of the radical addition rises almost 6 times from the intact fullerene to C60(OOtBu)4. This information may be important

fullerenyl radicals XC60• are characterized by almost the same distribution of spin density: its maximal values correspond to carbon atoms in positions 2 and 4 (relative to the site of the first X added) (Figure 1, Table 2). That is why we have scrutinized the possibility of reactions 2−4 (bisadduct formation, O2 addition or dimerization) for the mentioned reaction sites. Additionally, we have considered the 1,16bisadducts and 1,16,1′,16′-dimers as the products of reactions 2 and 4. In spite of the lower spin density in this position compared to sites 2 and 4, the 1,16-bisadducts have been identified among the products of addition of bulky radicals to C60.10 Heat Effects of Fullerene Bisadducts Formation via Reactions 2. Interaction of fullerenyls XC60• with the same X• is a paramount route to the diverse fullerene bis- and polyadducts.1,2 In the case of alkyl radicals (X = R), formation of 1,16- and 1,4-bisadducts is thermodynamically more favorable according to the calculations performed (Table 3). Note that absolute values of the heat effects are decreased two times from tBu• to Ph(CH3)2C•, a more voluminous radical. Due to the sterical hindrances, formation of 1,2-C60R2 is endothermic in both tert-butyl and cumyl cases. In the case of RO• and ROO•, formation of bisadducts is characterized by exothermic heat effects for all the studied reaction sites. However, the reactions resulting in both 1,2- and 1,4-bisadducts have the largest ΔHr° values that proves their

Figure 1. Fullerenyl radicals RC60•, ROC60•, and ROOC60• with R = tBu. Reaction sites of the further additions are designated; key bond lengths are shown in Å for R = tBu and R = Ph(CH3)2C (in parentheses). 13178



Article

Table 3. Heat Effects of Isomeric Bisadducts C60X2 Formation via Reaction 2a (kJ/mol) product reaction

1,2-C60X2

1,4-C60X2

1,16-C60X2

Bu + BuC60 Ph(CH3)2C • + Ph(CH3)2CC60• t BuO• + tBuOC60• Ph(CH3)2CO• + Ph(CH3)2COC60• t BuOO• + tBuOOC60• Ph(CH3)2COO• + Ph(CH3)2COOC60•

+15.7 +68.5 −136.6 −120.66 −102.6 −108.4

−104.9 −52.9 −144.1 −145.98 −94.8 −102.5

−114.9 −64.6 −108.4 −110.3 −58.9 −60.9

t

a

•

t

•

The highest ΔHr° values in each row are shown in bold.

why RC60• and ROC60• have demonstrated the different reactivity toward oxygen in experiments. Indeed, the slow reaction has been detected in the case of alkoxyfullerenyls (tBuO)nC60• 29 whereas in the case of alkylfullerenyl EtC60•, no chemical reaction has been found, and only the reversible formation of intermolecular complex has been observed.30 We propose that reactivity of XC60• toward oxygen can be enhanced by the use of its excited state (singlet) instead of the ground one. As the calculations indicate, 1O2 addition to carbon atoms 2 and 4 in fullerenyls is highly exothermic. Herewith, 1,4RC60OO•, 1,2-ROC60OO•, and 1,2-ROOC60OO• are the thermodynamically more favorable products in reactions 3. In the case of RC60•, addition to position 2 is unfavorable due to the steric hindrances: bulky tert-butyl and cumyl addends in RC60• screen this reaction site. In the case of ROOC60• and ROC60OO•, the voluminous R moieties are more remote from the site that facilitates oxygen addition to the mentioned position. Wave functions of oxygen and peroxides are multiconfigurational, so investigation of the reasons underlying low reactivity of fullerenyls toward oxygen demands multiconfigurational quantum-chemical methods.53 Unfortunately, these are currently not applicable to such big systems as fullerenes and their derivatives. Nevertheless, we have theoretically justified possibility of a new class of radicals XC60OO•: all the mentioned structures are minima on potential energy surface. Such radicals may be precursors for synthesis of fullerenecontataining polymers, in which C60 cores are linked by peroxide bridges. Radicals XC60OO• have been almost unstudied both theoretically and experimentally. The generation of (tBuO)nC60OO• has been confirmed by the evolution of the ESR signal in aerobic conditions only in ref 29. Usually, interaction of C-centered radicals with molecular oxygen occurs intensively.31 To compare, we have calculated heat effects of the analogous reactions of other radicals: 1,3-butadien-2-yl, phenyl (these radicals have the delocalized unpaired electrons and can be considered as distant analogues of XC60•), and Gomberg radical (a classic example of a radical with steric hindrances) (Scheme 1). The PBE/3ζ-calculated heat effects are exothermic in all these cases. Thus, XC60• radicals essentially differ from their hydrocarbon analogues in the aspect of their reactivity toward molecular oxygen. Most likely, it is caused by spreading the unpaired electron over the fullerene cage. Indeed, in the mentioned hydrocarbon radicals, the free spins are concentrated mainly on the one reactive site, so the spin densities are much higher (0.629−0.979, Scheme 1) in comparison with those of fullerenyl radicals (less than 0.450, Table 2). Heat Effects of Fullerenyl Dimerization and Polarizability of the Dimers Formed in Reactions 4.

Figure 2. Heat effects of the subsequent addition of tBuOO• to C60. Sites and the order of addition are shown on the structural formula.

to understanding the reasons of high efficiency of C60 as an antioxidant and inhibitor of radical processes because most of the known explanations deal with the large number of the reaction sites in its molecule.2 On the basis of the calculations performed, we can assume that its efficiency, in addition to the reason above, may be explained by the exothermicity increasing with each radical captured. Heat Effects of XC60• Reactions with Molecular Oxygen. We pay attention to the reactions between XC60• and O2 because molecular oxygen is always present in both in vitro and in vivo studies of antioxidant activity of C60 and its derivatives. The O-centered radicals XC60OO• should be products of the mentioned interactions. The calculated heat effects of reactions 3 with O2 in the ground state (triplet) are endothermic for all the XC60• radicals, regardless of the nature of X (Table 4). It explains the experimental observation of very slow interaction of fullerenyls with oxygen.29,30 Endothermicity of O2 addition to the position 2 decays in the series RC60• > ROC60• > ROOC60•. It explains Table 4. Heat Effects of the Reactions between Fullerenyls and Triplet Molecular Oxygen (Data for the Analogous Reactions with Singlet Oxygen in Parentheses) (kJ/mol) product •

reaction

1,2-XC60OO

BuC60• + O2 t BuOC60• + O2 t BuOOC60• + O2 Ph(CH3)2CC60• + O2 Ph(CH3)2COC60• + O2 Ph(CH3)2COOC60• + O2

+29.6 (−128.1) +5.5 (−152.2) +0.4 (−157.2) +30.3 (−127.3) a +2.6 (−155)

t

1,4-XC60OO• +10 (−147.7) +13.7 (−144) +11.7 (−146) +8.4 (−149.3) +26.8 (−130.8) +10.4 (−147.2)

a

Localization of the stationary point corresponding to this radical has been failed. 13179



Article

calculation) (Table 6). To find out how polarizability changes upon the dimerization, we have calculated the deviations from the additive scheme for each [1 + 1]-dimer:

Scheme 1. Heat Effects of Reactions between Hydrocarbon Radicals and Molecular Oxygena

Δα = αDFT − αadd

a

Referred to the basic definitions, we have expected that the deviations Δα, calculated according to eq 8, would be negative because highly polarizable unpaired electrons of two XC60• radicals vanish upon dimerization. The analogy with the hydrocarbon chemistry is suitable here. For example, according to PBE/3ζ calculations, Δα is negative (−1.5 Å3) for the dimer of the CH2C•CHCH2 radical. In contrast, Δα values for the dimers of fullerenyl radicals are positive, lying in the range 5.8−11.4 Å3. It means that exaltation of polarizability takes place in the case of [1 + 1]-dimers, analogous to their [2 + 2]counterparts. The maximal exaltation Δα has been found typical for RC60−C60R dimers. It is larger than 10 Å3 that is enough to be detected by the modern experimental techniques for polarizability measurements, e.g., interferometry, previously applied to measuring polarizability of polyfluorofullerenes.55 Probably, the magnitude of Δα of XC60−C60X can be regulated by variation of X, attached to the C60 core. From a fundamental point of view, the positive deviation of polarizability from the additive scheme is unusual in terms of the minimum polarizability principle.56,57 According to this principle, the thermodynamically more stable state of a molecular system should demonstrate lower polarizability. In the systems “dimer (A)−fullerenyl (B)”, state A is more stable and, consequently, the deviation Δα should be negative. However, the contrary situation is observed in the computations. Thus, the single-bonded fullerene dimers are one of the compounds violating the minimum polarizability principle (along with the noted in the previous work58). The calculated polarizabilities of fullerenyl radicals and the respective dimers allow proposing such a system “dimer (A)− fullerenyl (B)” (Scheme 2) for application in molecular machinery, similar to the molecular switch concept.59 The dimerization is reversible and XC60−C60X molecules facilely dissociate (e.g., under the thermal treatment), resulting in the respective radicals, which can be again dimerized under the appropriate conditions. At the same time, states A and B of such a fullerene-containing system significantly differ in mean polarizabilities; i.e., the states differ in the response to external electric fields, and this difference is measurable. We believe that the described systems may find an application in nanotechnologies.

Spin densities of the reaction sites are shown in red.

Dimerization of fullerenyl radicals is another possible process in “fullerene + radical” systems. The calculated heat effects of the dimerization are small (Table 5), which reflects their reversibility, previously found in experimental studies.2,4,9,32,33 Comparison of ΔHr° shows that the dimerization should result in 1,4,1′,4′-dimers because their formation via reaction 4 is the only exothermic mode. This perfectly agrees with the detection of such dimers among the products of the radical reactions.2,32 Noteworthy, the dimerization is thermodynamically more favorable for peroxyfullerenyls ROOC60• than for the respective RC60• and ROC60•. This fact can be useful for development of new synthetic strategies to obtain fullerene dimers with higher yields compared to the known methods (current synthetic approaches are commonly based only on the dimerization of RC60• radicals or their derivatives2,32). Moreover, peroxy groups of ROOC60−C60OOR dimers can be facilely involved in the further chemical transformations, similar to one-cage C60 peroxides,54 that provide opportunities for the functionalization of such bicage fullerene derivatives. Previously, the exaltation of dipole polarizability has been found for [2 + 2]-dimer (C60)2 and its derivatives.35 Caused by the interaction of π-electronic systems of two C60 cores, this phenomenon is proposed for the application to nanomaterial and nanodevice design. Because the polarizability (α) of [1 + 1]-dimers, formed at radical reactions, has not been studied, we have calculated the polarizability tensors and the mean polarizabilities of XC60−C60X. Mean polarizabilities of the dimers have been calculated by quantum-chemical means (αPBE/3ζ) and in terms of the additive scheme based on the polarizability of the respective fullerenyl radicals: αadd(XC60 − C60X) = 2α(XC60•)

(8)

■

(7)

CONCLUSION Alternative routes of decay for RC60•, ROC60•, and ROOC60• fullerenyl radicals, generated in “R• + C60” systems under

We have found that the mean polarizabilities of XC60−C60X are much higher compared to the respective values for the original radicals or the C60 fullerene (82.7 Å3, PBE/3ζ

Table 5. Heat Effects of Dimerization of Fullerenyl Radicals (kJ/mol) product radical

1,2,1′,2′−XC60−C60X

1,4,1′,4′−XC60−C60X

1,16,1′,16′−XC60−C60X

BuC60• t BuOC60• t BuOOC60• Ph(CH3)2CC60• Ph(CH3)2COC60• Ph(CH3)2COOC60•

a +38.1 −10.4 a a +4.9

−3.4 −52.9 −65.5 −2.7 −40.6 −58.2

+6.4 +9 +8.6 +2.7 +20.5 +3.9

t

a

Such products do not exist due to steric hindrances. 13180



Article

Table 6. Components of the Diagonalized Polarizability Tensors (αxx, αyy, and αzz), Additive Polarizability (αadd), Mean Polarizability (αPBE/3ζ), and Its Deviation from the Additive Scheme (Δα) for 1,4,1′,4′-XC60−C60X dimers (Å3) dimer

αxx

αyy

αzz

αPBE/3ζ

BuC60−C60tBu t BuOC60−C60OtBu t BuOOC60−C60OOtBu Ph(CH3)2CC60−C60C(CH3)2Ph Ph(CH3)2COC60−C60OC(CH3)2Ph Ph(CH3)2COOC60−C60OOC(CH3)2Ph

157.47 158.77 160.93 186.96 169.48 174.18

177.23 184.00 189.79 194.43 202.05 211.05

258.00 255.90 256.37 269.37 272.70 271.14

197.57 199.56 202.36 216.92 214.74 218.79

t

αadd 2 2 2 2 2 2

× × × × × ×

93.63 95.37 96.96 102.77 104.48 105.33

Δα 10.31 8.82 8.44 11.38 5.78 8.13

Scheme 2. Proposed Fullerene-Containing Molecular System with Two States Differing by Their Response to External Electric Fields

in the radical reactions of C60. First, these are compounds with outstandingly high mean polarizabilities, i.e., high response to electric fields. Second, the polarizability of such dimers is more than two times higher than the polarizability of the respective fullerenyl radicals. Taking into account the well-known reversibility of the dimerization, we have proposed a fullerene-containing molecular system that can exist in two reversible states (dimer and radical) with different responses to external electric field.

aerobic conditions, have been studied by the density-functional theory method PBE/3ζ. The obtained thermochemical data explain the observed experimental regularities of fullerene radical reactions and allow predicting new chemical and physical properties that may be useful for synthetic organic chemistry and nanotechnology applications. Formation of bisadducts C60X2 is the most exothermic route of decay for the radicals under study. A comparison of the heat effects shows that the most thermodynamically favorable products in these reactions are 1,16- and 1,4-C60R2 in the case of alkyls; 1,4- and 1,2-C60(OR)2 in the case of alkoxyls; and 1,2- and 1,4-C60(OOR)2 in the case of peroxy radicals. Here, we have found that exothermicity of addition of peroxy radicals to the fullerene framework is increased during the subsequent addition. We will apply this interesting fact to a study on the mechanisms of C60 antioxidant activity. Much lower exothermic heat effects are typical for the dimerization of fullerenyls resulting in 1,4,1′,4′-dimers (the other variants of dimerization are endothermic according to the calculations performed). Herewith, the heat effect is maximal for the dimerization of ROOC60• among the other cases. Therefore, the interaction of C60 with peroxy radicals can be recommended as one of the way to increase the yields of the single-bonded C60 dimers instead of the common strategy based on the RC60• dimerization. The interaction with molecular oxygen is the less probable route of decay for fullerenyl radicals due to its endothermicity. Their reactivity toward O2 rises in the series RC60• < ROC60• < ROOC60• and can be enhanced by the use of the singlet molecular oxygen. This way, the O-centered radicals 1,4RC60OO•, 1,2-ROC60OO•, and 1,2-ROOC60OO• should be generated by highly exothermic addition of 1O2 to XC60•. In the context of new materials and nanotechnology, we pay particular attention to the polarizability of the dimers, formed

■

AUTHOR INFORMATION

Corresponding Author

*D. Sh. Sabirov: e-mail, [email protected]. Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS This work was supported by the Presidium of the Russian Academy of Sciences (program No. 24 “Foundations of Basic Research of Nanotechnologies and Nanomaterials”).

■

REFERENCES

(1) Krusic, P. J.; Wasserman, E.; Keizer, P. N.; Morton, J. R.; Preston, K. F. Radical Reactions of C60. Science 1991, 254, 1183−1185. (2) Tzirakis, M. D.; Orfanopoulos, M. Radical Reactions of Fullerenes: From Synthetic Organic Chemistry to Materials Science and Biology. Chem. Rev. 2013, 113, 5262−5321. (3) Morton, J. R.; Preston, K. F.; Krusic, P. J.; Hill, S. A.; Wasserman, E. ESR Studies of the Reaction of Alkyl Radicals with C60. J. Phys. Chem. 1992, 96, 3576−3578. (4) Tumanskii, B. L.; Bashilov, V. V.; Bubnov, N. N.; Solodovnikov, S. P.; Sokolov, V. I. ESR Study of the Reversible Dimerization of Phosphonylfullerenyl Radicals. Bull. Russ. Acad. Sci. Div. Chem. Sci. 1992, 41, 1519−1520. 13181



Article

Theoretical Study of HC60 and H3C60 Radicals. Can. J. Chem. 1994, 72, 776−782. (24) Morton, J. R.; Negri, F.; Preston, K. F. Addition of Free Radicals to C60. Acc. Chem. Res. 1998, 31, 63−69. (25) Morton, J. R.; Preston, K. F.; Krusic, P. J.; Wasserman, E. Electron Paramagnetic Resonance Spectra of R-C60 Radicals. Evidence for RC60C60R Dimers. J. Chem. Soc., Perkin Trans. 2 1992, 1425−1429. (26) Sabirov, D. Sh.; Bulgakov, R. G. Reactivity of Fullerenes Family towards Radicals in Terms of Local Curvature. Comput. Theor. Chem. 2011, 963, 185−190. (27) Dresselhaus, M. S.; Dresselhaus, G.; Eklund, P. C. Science of Fullerenes and Carbon Nanotubes; Elsevier Inc.: Amsterdam, 1996; pp 965. (28) Reed, C. A.; Bolskar, R. D. Discrete Fulleride Anions and Fullerenium Cations. Chem. Rev. 2000, 100, 1075−1120. (29) Zhang, Y.-K.; Janzen, E. G.; Kotake, Y. Stabilities and Reactivities of Buckminsterfullerene Radicals, (ButO)nC60•, Towards Dioxygen, Nitric Oxide and Spin Trapping Agents. J. Chem. Soc., Perkin Trans. 2 1996, 1191−1196. (30) Bulgakov, R. G.; Ponomareva, Yu. G.; Sadykov, R. A. Chemiluminescence upon the Dimerization of EtC60• Fullerenyl Radicals. Russ. Chem. Bull., Int. Ed. 2008, 57, 2028−2029. (31) Denisov, E. T.; Afanas’ev, I. B. Oxidation and Antioxidants in Organic Chemistry and Biology; CRC Press: Boca Raton, FL, 2005; pp1024. (32) Segura, J. L.; Martin, N. Fullerene Dimers. Chem. Soc. Rev. 2000, 29, 13−25. (33) Morton, J. R.; Preston, K. F.; Krusic, P. J.; Hill, S. A.; Wasserman, E. The Dimerization of Fullerene RC60 Radicals [R = alkyl]. J. Am. Chem. Soc. 1992, 114, 5454−5455. (34) Lu, S.; Jin, T.; Bao, M.; Yamamoto, Y. NaOH-Catalyzed Dimerization of Monofunctionalized Hydrofullerenes: TransitionMetal-Free, General, and Efficient Synthesis of Single-Bonded [60]Fullerene Dimers. Org. Lett. 2012, 14, 3466−3469. (35) Sabirov, D. Sh. Polarizability of C60 Fullerene Dimer and Oligomers: The Unexpected Enhancement and Its Use for Rational Design of Fullerene-Based Nanostructures with Adjustable Properties. RSC Adv. 2013, 3, 19430−19439. (36) Fang Ma, F.; Li, Zh.-R.; Zhou, Zh.-J.; Wu, D.; Li, Y.; Wang, Y.F.; Li, Z.-Sh. Modulated Nonlinear Optical Responses and Charge Transfer Transition in Endohedral Fullerene Dimers Na@C60C60@F with n-Fold Covalent Bond (n = 1, 2, 5, and 6) and Long Range Ion Bond. J. Phys. Chem. C 2010, 114, 11242−11247. (37) Shestakov, A. F. Reactivity of Fullerene C60. Russ. J. Gen. Chem. 2008, 78, 811−821. (38) Sabirov, D. Sh.; Bulgakov, R. G.; Khursan, S. L. Indices of the Fullerenes Reactivity. ARKIVOC 2011, No. 8, 200−224. (39) Osuna, S.; Swart, M.; Solà, M. The Reactivity of Endohedral Fullerenes. What Can Be Learnt From Computational Studies? Phys. Chem. Chem. Phys. 2011, 13, 3585−3603. (40) Gan, L.; Huang, Sh.; Zhang, X.; Zhang, A.; Cheng, B.; Cheng, H.; Li, X.; Shang, G. Fullerenes as a tert-Butylperoxy Radical Trap, Metal Catalyzed Reaction of tert-Butyl Hydroperoxide with Fullerenes, and Formation of the First Fullerene Mixed Peroxides C60(O)(OOtBu)4 and C70(OOtBu)10. J. Am. Chem. Soc. 2002, 124, 13384− 13385. (41) Bulgakov, R. G.; Kinzyabaeva, Z. S. A New Application of Cerium(IV) Ammonium Nitrate [(NH4)2Ce(NO3)3] in Fullerene Chemistry − Free-Radical Alkoxylation of C60 Fullerene on Reaction with Primary Alcohols. Tetrahedron Lett. 2012, 53, 5781−5783. (42) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (43) Laikov D. N. Ph.D. Dissertation, Moscow State University, 2000. (44) Laikov, D. N.; Ustynyuk, Yu. A. PRIRODA-04: a QuantumChemical Program Suite. New Possibilities in the Study of Molecular Systems with the Application of Parallel Computing. Russ. Chem. Bull., Int. Ed. 2005, 54, 820−826. (45) Misochko, E. Ya.; Akimov, A. V.; Belov, V. A.; Tyurin, D. A.; Laikov, D. N. High Resolution EPR Spectroscopy of C60F and C70F in

(5) Tumanskii, B. L.; Bashilov, V. V.; Solodovnikov, S. P.; Sokolov, V. I. EPR Investigation of the Adducts of Element-Centered Radicals with Polyhedral Carbon Clusters (Fullerenes). Bull. Russ. Acad. Sci. Div. Chem. Sci. 1992, 41, 1140−1141. (6) Sokolov, V. I. The Problem of Fullerenes. The Chemical Aspect. Russ. Chem. Bull. Int. Ed. 1993, 42, 1−11. (7) Krusic, P. J.; Roe, D. C.; Johnston, E.; Morton, J. R.; Preston, K. F. EPR Study of Hindered Internal Rotation in Alkyl-C60 Radicals. J. Phys. Chem. 1993, 97, 1736−1738. (8) Lyakhovetsky, Y. I.; Shilova, E. A.; Bashilov, V. V.; Sokolov, V. I.; Nekrasov, Y. S.; Tumanskii, B. L. Homolytic Reactive Mass Spectrometry of Fullerenes: Interaction of C60 and C70 with Organoand Organoelement Mercurials in the Electron Impact Ion Source of a Mass Spectrometer; EPR, CIDEP, and MS Studies of Several Analogous Reactions of C60 Performed in Solution. J. Phys. Chem. A 2009, 113, 13700−13710. (9) Fagan, P. J.; Krusic, P. J.; McEwen, C. N.; Lazar, J.; Parker, D. H.; Herron, N.; Wasserman, E. Production of Perfluoroalkylated Nanospheres from Buckminsterfullerene. Science 1993, 262, 404−407. (10) Troshin, P. A.; Lyubovskaya, R. N. Organic Chemistry of Fullerenes: The Major Reactions, Types of Fullerene Derivatives and Prospects for Their Practical Use. Russ. Chem. Rev. 2008, 77, 305−349. (11) Bulgakov, R. G.; Ponomareva, Yu. G. Effect of the Medium on the Reactivity of Fullerene C60 with Respect to the Peroxide Radicals RO2•.Chemiluminescence in the C60−AIBN−O2−C2H5Ph−PhH System. Russ. Chem. Bull., Int. Ed. 2009, 58, 786−790. (12) Galimov, D. I.; Bulgakov, R. G.; Gazeeva, D. R. Reactivity of Fullerene C60 Towards Peroxy Radicals Generated by Liquid-Phase Oxidation of Cumene and Ethylbenzene with Oxygen. Russ. Chem. Bull., Int. Ed. 2011, 60, 2107−2109. (13) Bulgakov, R. G.; Galimov, D. I.; Gazeeva, D. R. Addition of the Peroxyl Radicals RO2• to C70, C60 Molecules − the Dominant Reaction Inhibition by Fullerenes of the Hydrocarbons Oxidation. Fullerenes Nanotubes Carbon Nanostruct. 2013, DOI: 10.1080/ 1536383X.2012.684182. (14) Cataldo, F. Interaction of C60 Fullerene with Lipids. Chem. Phys. Lipids 2010, 163, 524−529. (15) Yumagulova, R. Kh.; Biglova, Yu. N.; Kuznetsov, S. I.; Ionova, I. A.; Kolesov, S. V.; Monakov, Yu. B. Fullerene C60 in Copolymerization of Allyl Chloride with Methyl Methacrylate. Polymer Sci., Ser. A 2008, 50, 260−264. (16) Troitskii, B. B.; Troitskaya, L. S.; Yachnov, A. S.; Dmitriev, A. A.; Denisova, V. N.; Novikova, M. A.; Anikina, L. I.; Domrachev, G. A. A Study of Fullerenes as High-Temperature Antioxidants of Poly(methyl methacrylate) and Polystyrene. Dokl. Phys. Chem. 1998, 363, 367−369. (17) Ginzburg, B. M.; Shibaev, L. A.; Ugolkov, V. L. Effect of Fullerene C60 on Thermal Oxidative Degradation of Polymethyl Methacrylate Prepared by Radical Polymerization. Russ. J. Appl. Chem. 2001, 74, 1329−1337. (18) Zeynalov, E. B.; Allen, N. S.; Salmanova, N. I. Radical Scavenging Efficiency of Different Fullerenes C60−C70 and Fullerene Soot. Polym. Degrad. Stab. 2009, 94, 1183−1189. (19) Gasanov, R. G.; Kalina, O. G.; Bashilov, V. V.; Tumanskii, B. L. Addition of Carbon-Centered Radicals to C60. Determination of the Rate Constants by the Spin Trap Method. Russ. Chem. Bull., Int. Ed. 1999, 48, 2344−2346. (20) Cremonini, M. A.; Lunazzi, L.; Placucci, G.; Krusic, P. J. Addition of Alkylthiyl and Alkoxy Radicals to C60 Studied by ESR. J. Org. Chem. 1993, 58, 4735−4738. (21) Borghi, R.; Guidi, B.; Lunazzi, L.; Placucci, G. ESR Detection of the Regioisomers Due to Addition of Methoxy and Methylthio Radicals to Fullerene C70. J. Org. Chem. 1996, 61, 5667−5669. (22) Borghi, R.; Lunazzi, L.; Placucci, G.; Cerioni, G.; Plumitallo, A. Photolysis of Dialkoxy Disulfides: A Convenient Source of Alkoxy Radicals for Addition to the Sphere of Fullerene C60. J. Org. Chem. 1996, 61, 3327−3331. (23) Morton, J. R.; Negri, F.; Preston, K. F. Addition of Alkyl Radicals to C60. Part 3. The EPR Spectra of R3C60 Radicals and a 13182



Article

Solid Argon: Reassignment of C70F Regioisomers. J. Chem. Phys. 2007, 127, 084301. (46) Sabirov, D. Sh.; Khursan, S. L.; Bulgakov, R. G. Quantum Chemical Modeling of Ozone Addition to C60 Fullerene. Fullerenes Nanotubes Carbon Nanostruct. 2008, 16, 534−537. (47) Sabirov, D. Sh.; Khursan, S. L.; Bulgakov, R. G. 1,3-Dipolar Addition Reactions to Fullerenes: The Role of the Local Curvature of Carbon Surface. Russ. Chem. Bull., Int. Ed. 2008, 57, 2520−2525. (48) Pankratyev, E. Yu.; Tulyabaev, A. R.; Khalilov, L. M. How Reliable are GIAO Calculations of 1H and 13C NMR Chemical Shifts? A Statistical Analysis and Empirical Corrections at DFT (PBE/3z) Level. J. Comput. Chem. 2011, 32, 1993−1997. (49) Sabirov, D. Sh.; Bulgakov, R. G. Polarizability of OxygenContaining Fullerene Derivatives C60On and C70O with Epoxide/ Oxidoannulene Moieties. Chem. Phys. Lett. 2011, 506, 52−56. (50) Bulgakov, R. G.; Galimov, D. I.; Sabirov, D. Sh. New Property of the Fullerenes: the Anomalously Effective Quenching of Electronically Excited States Owing to Energy Transfer to the C70 and C60 Molecules. JETP Lett. 2007, 85, 632−635. (51) Sabirov, D. Sh.; Bulgakov, R. G. Polarizability Exaltation of Endofullerenes X@Cn (n = 20, 24, 28, 36, 50, and 60; X is a Noble Gas Atom). JETP Lett. 2010, 92, 662−665. (52) Sabirov, D. Sh.; Saitgalina, A. D.; Garipova, R. R.; Bulgakov, R. G. Molecular Simulation on Fullerene Derivatives: Algorithms and Programs for Calculation of Their Curvature, Polarizability and Volume. Book of Abstracts “Computer Simulation of Advanced Materials”, 16−21 July 2012; Lomonosov Moscow State University: Moscow, Russia, 2012; p 35. (53) Antonovskii, V. L.; Khursan, S. L. Physical Chemistry of Organic Peroxides; CRC Press: Boca Raton, FL, 2005; p 549. (54) Huang, Sh.; Xiao, Z.; Wang, F.; Gan, L.; Zhang, X.; Hu, X.; Zhang, S.; Lu, M.; Pan, Q.; Xu, L. Selective Preparation of OxygenRich [60]Fullerene Derivatives by Stepwise Addition of tertButylperoxy Radical and Further Functionalization of the Fullerene Mixed Peroxides. J. Org. Chem. 2004, 69, 2442−2453. (55) Hornberger, K.; Gerlich, S.; Ulbricht, H.; Hackermüller, L.; Nimmrichter, S.; Goldt, I. V.; Boltalina, O.; Arndt, M. Theory and Experimental Verification of Kapitza−Dirac−Talbot−Lau Interferometry. New J. Phys. 2009, 11, 043032. (56) Ghanty, T. K.; Ghosh, S. K. A Density Functional Approach to Hardness, Polarizability, and Valency of Molecules in Chemical Reactions. J. Phys. Chem. 1996, 100, 12295−12298. (57) Holm, U. Is There a Minimum Polarizability Principle in Chemical Reactions? J. Phys. Chem. A 2000, 104, 8418−8423. (58) Blair, S. A.; Thakkar, A. J. How Often Is the Minimum Polarizability Principle Violated? Chem. Phys. Lett. 2013, 556, 346− 349. (59) Mateo-Alonso, A.; Guldi, D. M.; Paolucci, F.; Prato, M. Fullerenes: Multitask Components in Molecular Machinery. Angew. Chem., Int. Ed. 2007, 46, 8120−8126.

13183


Density Functional Theory Study on the Decay of Fullerenyl Radicals

Recommend Documents