JCHEMJournal of Chemistry2090-90712090-9063Hindawi10.1155/2017/54940385494038Research ArticleFive Regioisomers of Dimethyl Dodecahedrane Derivatives: A Hybrid DFT B3LYP StudyLeeSeol1ChoYongJae1http://orcid.org/0000-0003-3105-9753LeeKee Hag1Espinosa FeraoArturoDepartment of Chemistry and Nanoscale Sciences and Technology InstituteWonkwang UniversityIksanJeonbuk 570-749Republic of Koreawku.ac.kr
Dedicated to Professor Roald Hoffmann on the occasion of his 80th birthday
The hybrid density functional (B3LYP/6-31G(d, p)) method was used to understand why 1,2-dimethyl dodecahedrane has not been reported yet. From the direct dimethyl substitution of the dodecahedrane cage, we could have five C20H18(CH3)2 derivatives which were geometrically optimized without constraints. The results suggest that 1,16-dimethyl dodecahedrane derivative is the most stable, whereas the 1,2-dimethyl derivative is the most unstable; this implies that the distortion due to steric hindrance interactions would be a crucial effect on the relative energies of the dimethyl dodecahedranes. It would be disadvantaged thermodynamically for 1,2-dimethyl derivative that, in the experimental results, was not synthesized yet. The LUMO of each of the derivatives was equivalently delocalized over the void within the cage, implying that dimethyl derivatives are able to encapsulate atoms at the center of the cage. The HOMO was limitedly delocalized on the cage. The characteristics of the HOMO of derivatives show three patterns, implying that each derivative might undergo one of three entirely different sets of characteristic chemical reactions with electrophilic reagents.
Wonkwang University1. Introduction
The geometry of fullerenes obeys the isolated-pentagon rule. However, this rule cannot be satisfied for the smallest fullerene C20, which contains twelve 5-membered rings (MRs) and zero 6-MRs. Therefore, the C20 fullerene would be highly reactive owing to the fusion of strained 5-MRs. Furthermore, this strain could be removed through reactions of C20 that form sp3 bonds. By using density functional theory method, the structures and electronic states of alkyl radicals (R, CnH2n+1, n=1–4) added to C20 fullerenes were investigated [1]. Also it was proposed that thiyl radical-promoted polycyclization to form dodecahedrane would be exergonic [2].
Density functional calculations of the reactions CX3 radical + C20H20 (X = H, F, Cl, and Br) based on two pathways (H-displacement and H-abstraction from C20H20) showed that CH3 radical prefers H-abstraction from C20H20 while F-displacement is favorable for the reaction of CF3 radical with the fullerene. Exothermic characters of H-abstraction with the lower potential barrier indicate that the H-abstraction would dominate the reaction CCl3∙ + C20H20 while endothermic characters with high potential barrier heights for two reaction pathways of C20H20 with CBr3∙ indicate unfavorable thermodynamically and kinetically pathways for these reactions [3].
Using the PBE0/cc-pVDZ quantum chemical method, dipole moments, polarizabilities, the first hyperpolarizabilities for 1,2-diethynyldodecahedrane C20H18(CCH)2 and diethynylmethanododecahedrane C20H18C(CCH)2 including two neighboring ethynyl groups in vacuum and in tetrahydrofuran (THF) were reported [4].
A dodecahedrane derivative, namely, the D3d symmetric 1,16-dimethyl derivative, has been synthesized experimentally [5–8]. Since then, the synthesis of methyl-substituted dodecahedrane was also confirmed by X-ray analysis. Moreover, Pólya’s theorem was applied to the enumeration of substituted dodecahedrane isomers [9]. The total synthesis of dodecahedrane was confirmed by the single-peak nature of the 1H and 13C NMR spectra [10], as well as the X-ray spectra [11, 12]. Methyl dodecahedrane, 1,4-dimethyl (Cs), 1,6-dimethyl (C2v), 1,16-dimethyl (D3d) dodecahedranes, and trimethyl dodecahedrane have subsequently been reported [13–15]. However, at present no 1,2-dimethyl (C2v) and 1,7-dimethyl (C2) dodecahedranes have been reported, yet.
The structures of dodecahedrane and its derivatives were studied at the level of molecular mechanics calculations [16, 17], INDO calculations [18], and hybrid density functional B3LYP calculations [19]. Recently the relative energies, electronic properties, and atomic structures of disubstituted C20H18X2 regioisomers (X = F, Cl, Br, or OH) at the level of B3LYP/6-31G(d,p) have been determined [20]. Here the predominant electronic configurations of the regioisomers with all substituents X which are higher electron negativity than carbon are not so different, but the HOMO and LUMO maps for five regioisomers are not similar, implying that different exohedral complexes may undergo a distinct set of characteristic chemical reactions.
However, to the best of our knowledge, no the first principle calculations of the full optimization of the five dimethyl dodecahedranes, C20H18(CH3)2 derivatives using B3LYP/6-31G(d,p), have been reported. Such calculations would be required in order to understand, in the point of energetics, why only three derivatives have been synthesized. Accordingly, in this study we determine the relative energies, electronic properties, and atomic structures of five derivatives with the same spin state, following the methods of Hwang et al. 2012 [20]. Also it is interesting to see the effect of the methyl substituents on the frontier orbitals, because carbon atoms have the same electronegativity in the substituents and the cage.
2. Calculations
In this study, hybrid density functional theory (DFT) with Becke’s three-parameter hybrid method and the Lee-Yang-Parr exchange-correlation functional theory (B3LYP) [21–28] were used to optimize the geometries of the C20H20 and five dimethyl derivatives C20H18(CH3)2, as shown in Figure 1. Furthermore, the electron basis set 6-31G(d,p) was used [25]. We fully optimized the geometries of all the C20H18(CH3)2 derivatives using the Gaussian 2003 B.04 and Gaussian 2009 D01 package suites [26, 27]. We obtain highly accurate geometries, by using the convergence criterion with tight optimization and an ultrafine pruned (99,590) grid. Here, we analyzed the relative energies as well as the HOMO and LUMO orbitals of the derivatives [20, 28].
(a) The dodecahedrane and its numbering system (hydrogen atoms deleted for clarity). (b) The dimethyl substituted dodecahedrane derivatives C20H18(CH3)2.
To attest the veracity of the B3LYP/6-31G(d,p) relative energy values, as shown in Supplementary Material (see Supplementary Material available online at https://doi.org/10.1155/2017/5494038), we calculated single-point energies at the levels of M06-2X [29], PBE [30], B98 [31], B97D3 [32, 33], B3LYPD3 [21, 22, 33], PBED3 [30, 33], and M06-2XD3 [29, 33] with 6-31G(d,p) basis set on the optimized geometries and using the Zero Point Energy-Correction (ZPEC) already obtained at the B3LYP/6-31G(d) level, using Gaussian 2009 D01 package suites [27]. Here D3 represents the dispersion-corrected functional with the D3 version of Grimme’s dispersion with Becke-Johnson damping [33]. Also using the post-HF MP2 [34]/def2-TZVPP [35]//B3LYP/6-31G(d) level, we computed the energies as single-point calculations with the ZPEC already obtained at the B3LYP/6-31G(d) level.
3. Results and Discussion
By performing calculations for the full optimization of the atomic structures of the five C20H18(CH3)2 derivatives using B3LYP/6-31G(d,p) without any constraints, we determined the relative energies of the derivatives in Table 1 with the symmetries as shown in Figure 1.
The relative energies (rel. E), relative ZPEC energies (rel. ZPEC), relative enthalpy (rel. ΔH0), and relative Gibbs free energies (rel. ΔG0) in units of meV and distance (d, Å) between the two methyls of the five C20H18(CH3)2 derivatives using B3LYP/6-31G(d, p) calculationsa.
Derivative
Rel. E
Rel. ZPEC
ΔH0
ΔG0
d
1,2-cis-1
112.7
121.6
3
9
2.88
1,4-cis-2
2.9
2.5
−3
−5
4.32
1,7-eq
0.8
0.4
−2
−2
6.11
1,6-trans-2
1.9
0.2
−4
−7
6.99
1,16-trans-1
0.0
0.0
0
0
7.48
aThe relative value of each derivative is the energy minus the energy of trans-1 isomer.
The geometrical structures and the frontier orbitals (HOMOs and LUMOs) of the five C20H18(CH3)2 derivatives are shown in Figure 2. The calculated C-C bond length in the cage of neutral C20H20 with Ih symmetry is 1.556 Å, which is in the range of the experimental result (i.e., 1.555±0.001 Å) [36].
The (a) HOMOs and (b) LUMOs of the fully optimized geometries of the five C20H18(CH3)2 derivatives using B3LYP/6-31G(d,p) calculations. Here X-, Y-, and Z-axis represent the frontier orbitals, looking from x, y, and z direction, respectively. The HOMOs of C20H18(CH3)2 derivatives (isovalue = 0.05). trans-2 (HOMO: −7.08 eV) and the others (HOMO: −7.05); C20H20 (HOMO: −7.09 eV). The LUMOs of C20H18(CH3)2 derivatives (isovalue = 0.05). Five dimethyl derivatives (LUMO: 0.98 eV); C20H20 (LUMO: 0.90 eV).
In the five dimethyl derivatives whose geometries were fully optimized, the relative Zero Point Energy-Correction (ZPEC) energies of dimethyl dodecahedranes are arranged in the following order: 1,16-dimethyl derivative (trans-1, D3d) ≈ 1,7-dimethyl derivative (equatorial, C2) ≈ 1,6-dimethyl derivative (trans-2, C2v) ≈ 1,4-dimethyl derivative (cis-2, Cs) < 1,2-dimethyl derivative (cis-1, C2v). The distance between dimethyl groups in the 1,16-dimethyl derivative with the lowest energy is the longest among these. In 1,2-dimethyl derivative, the two substituents are located on the surface of the C20H20 fullerene within a distance that is shorter than the van der Waals distance (i.e., 3.4 Å).
The increasing order of relative ZPEC energies in them is the same to the decreasing order of the distance between the two methyl groups of derivatives. As shown in Figure 3 and Supplementary Materials, the veracity of the B3LYP/6-31G(d,p) relative energy values is attested at the levels of M06-2X, PBE, B98, B97D3, B3LYPD3, PBED3, and M06-2XD3, and MP2 in the point of view that the relative energy of cis-1 is about 10 kJ/mol higher than the others. The difference of relative energies of four isomers except cis-1 is less than 0.7 kJ/mol in the levels of DFT but within 2 kJ/mol in the MP2 method.
Trends in relative energies (unit: kJ/mol) for five isomers computed as single-point calculations at all different levels of DFT/6-31G(d, p)//B3LYP/6-31G(d, p) and MP2/def2-TZVPP//B3LYP/6-31G(d, p).
Our result suggests that both the strength of the cage distortion and the van der Waals force between the pair of methyl groups of derivatives are important factors that affect the position of exohedral dimethyl substitution on the dodecahedrane cage. And it could be suggested that in 1,2-dimethyl derivative, steric impedance brought on by the nearby two methyl groups is very crucial. It is interesting to compare the results above with the case of C20H18X2 (X = F, Cl, Br, OH). In this case, except for the dihydroxy derivative, 1,2-disubstituted derivatives have the highest energy whereas 1,16-derivatives have the lowest energy; however, in the dihydroxy derivatives, the increasing order is reversed [20, 37].
Here our calculations concerning the dimethyl dodecahedrane derivatives could explain the fact why three dimethyl dodecahedrane derivatives were reported in X-ray and NMR experiments [7, 15], but the 1,2-dimethyl dodecahedrane derivative has not been synthesized yet. And it could be suggested that the more stable derivative in the viewpoint of energetics could be synthesized easier than others.
In the following substitution reaction, the thermal reaction energies are endothermic:(1)C20H20+2CH3⟶C20H18CH32+2Hwhere 2.48, 2.86, 2.56, 3.78, and 5.79 eV of the thermal reaction energies at 298.15 K and 1 atm are absorbed for 1,16-dimethyl derivative (trans-1, D3d), 1,7-dimethyl derivative (equatorial, C2), 1,6-dimethyl derivative (trans-2, C2v), 1,4-dimethyl derivative (cis-2, Cs), and 1,2-dimethyl derivative (cis-1, C2v), respectively. The Gibbs free energies are 10.285, 10.278, 10.283, 10.280, and 10.294 eV for 1,16-dimethyl derivative (trans-1, D3d), C 1,6-dimethyl derivative (trans-2, C2v), 1,7-dimethyl derivative (equatorial, C2), 1,4-dimethyl derivative (cis-2, Cs), and 1,2-dimethyl derivative (cis-1, C2v), respectively.
Even though the very small variation of Gibbs free energy is not sufficient to decide which structure should be totally experimentally disfavored, it could suggest that, in the point of thermodynamics, a synthesis of cis-1 isomer might be a little difficult compared to the others. In the experimental results, dodecahedrane, methyl, and dimethyl derivatives are synthesized [12], but cis-1 and equatorial isomers are not synthesized yet.
Figure 2 shows the HOMOs and LUMOs for the directed dimethyl derivatives of the dodecahedrane cage. When we calculate the energy gaps from the HOMO/LUMO energy levels of the derivatives, C20H20 and C20H18(CH3)2 have slightly different energy gaps. Here the HOMOs of the C20H18(CH3)2 derivatives with the same number of symmetry elements are very similar, but the HOMOs of derivatives with the different number of symmetry elements are different, despite the fact that the LUMOs among different geometric structures are identical. This suggests that each derivative undergoes one of three distinct sets of characteristic chemical reactions only for the HOMO. It is interesting to see that they are different to the case of disubstituted C20H18X2 regioisomers (X = F, Cl, Br, or OH) which have five distinct set of reactions for both the HOMO and the LUMO [20].
Figure 2 shows the effect of the disubstituents on the frontier orbitals of C20H20. The HOMO and LUMO energies of five C20H18(CH3)2 derivatives are higher than the HOMO and LUMO levels of C20H20 (HOMO: −7.09 eV; LUMO: 0.90 eV). Also the band gaps of dimethyl substituted derivatives are wider than the gaps of C20H20 and C20 (HOMO: −5.03 eV; LUMO: −3.10 eV). Moreover, the LUMO energies are positive and higher than that of C20H20, which implies that these could function less efficiently as electron acceptors than C20H20.
4. Conclusion
We obtained the fully optimized geometries of the C20H18(CH3)2 derivatives (the smallest fullerene disubstituents) without any constraints using the hybrid B3LYP/6-31G(d,p) density functional method. We observed that the relative energy of the C20H20 disubstituted derivatives is mostly dependent on the distance between the pair of two methyl groups. It is interesting to note that the 1,2-dimethyl derivatives is less stable than the 1,16-dimethyl derivative. This implies that, in the C20H18(CH3)2 derivatives, the effect of distortion due to the steric hindrance interactions on the relative energy would be crucial. Therefore, our calculations in the point of thermodynamic view could support the experimental results, wherein three dimethyl dodecahedrane derivatives except the 1,2-dimethyl dodecahedrane derivative were preparative-synthetic and suggest that 1,7-dimethyl derivative (equatorial, C2) would be synthesized in the near future.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was supported by Wonkwang University in 2015.
AbeS.KawanoS.ToidaY.NakamuraM.InoueS.SanoH.YoshidaY.KawabataH.TachikawaH.Electronic states of alkyl-radical-functionalized C20 fullerene using density functional theory20165532-s2.0-8495989303710.7567/JJAP.55.03DD0303DD03WangS. C.TantilloD. J.Viability of dodecahedrane-forming radical polycyclizations20171591976197910.1039/C6OB02300BGhafouriR.EktefaF.Theoretical study on the mechanism of reactions of CX3 radicals (X = H, F, Cl and Br) with C20H20 and C20F20 fullerenes201711272963022-s2.0-8498268438110.1016/j.molstruc.2016.07.112SemenovS. G.MakarovaM. V.A quantum chemical study of diethynyl derivatives of dodecahedrane and buckminsterfullerene in vacuum and in tetrahydrofuran20151181464910.1134/S0030400X1501021XPaquetteL. A.BaloghD. W.UshaR.KountzD.ChristophG. G.Crystal and molecular structure of a pentagonal dodecahedrane198121144825755762-s2.0-000060672610.1126/science.211.4482.575PaquetteL. A.BaloghD. W.An expedient synthesis of 1,16-dimethyldodecahedrane198210437747832-s2.0-000000197110.1021/ja00367a021ChristophG. G.EngelP.UshaR.BaloghD. W.PaquetteL. A.Comparative X-ray crystal structure analysis of C2v-dimethylmonosecododecahedrane and D3d-dimethyldodecahedrane. The remarkable consequences of a missing framework bond on molecular topology198210437847912-s2.0-034263204310.1021/ja00367a022PaquetteL. A.Dodecahedrane—The chemical transliteration of Plato's universe (A Review)198279144495450010.1073/pnas.79.14.4495PaquetteL. A.TernanskyR. J.BaloghD. W.TaylorW. J.Total synthesis of a monosubstituted dodecahedrane. The methyl derivative198310516544154462-s2.0-3384555193810.1021/ja00354a042PaquetteL. A.TernanskyR. J.BaloghD. W.KentgenG.Total synthesis of dodecahedrane198310516544654502-s2.0-3384555097910.1021/ja00354a043PrinzbachH.FessenerW.-D.1987Oxford, UKBlackwellBertauM.WahlF.WeilerA.ScheumannK.WörthJ.KellerM.PrinzbachH.From pagodanes to dodecahedranes - Search for a serviceable access to the parent (C20H20) hydrocarbon1997532910029100402-s2.0-003078962810.1016/S0040-4020(97)00345-1PaquetteL. A.MiyaharaY.DoeckeC. W.Preparatively useful dehydrogenative method for dodecahedrane synthesis19861087171617182-s2.0-034478768310.1021/ja00267a067SantasI.BaloghD. W.DoeckeC. W.MarshallA. G.PaquetteL. A.Gas-phase basicities of cubane, dodecahedrane, and methyl-and 1,16-dimethyldodecahedrane as measured by Fourier transform ion cyclotron resonance mass spectrometry198610826818381852-s2.0-000100935610.1021/ja00286a011PaquetteL. A.MiyaharaY.Directed alkyl substitution of the dodecahedrane nucleus. The 1,4-dimethyl, 1,6-dimethyl, and 1,4,16-trimethyl derivatives. Indanododecahedrane by stepwise dehydrogenation of a benzylated seco derivative19875271265127210.1021/jo00383a017AllingerN. L.GeiseH. J.PyckhoutW.PaquetteL. A.GallucciJ. C.Structures of norbornane and dodecahedrane by molecular mechanics calculations (MM3), X-ray crystallography, and electron diffraction19891113110611142-s2.0-000142281510.1021/ja00185a048WahlF.WörthJ.PrinzbachH.The Pagodane Route to Dodecahedranes: An Improved Approach to the C20H20 Parent Framework; Partial and Total Functionalizations—Does C20-Fullerene Exist?199332121722172610.1002/anie.199317221SchulmanJ. M.DischR. L.Theoretical studies of dodecahedrane. 2. Dodecahedrane, inclusion compounds, and fluorine derivatives197810018567756812-s2.0-000041143710.1021/ja00486a016MoranD.StahlF.JemmisE. D.SchaeferH. F.IIISchleyerP. V. R.Structures, stabilities, and ionization potentials of dodecahedrane endohedral complexes200210620514451542-s2.0-003716217610.1021/jp014471zHwangY. G.LeeS.LeeK. H.DFT Study for Substitution Patterns of C20H18X2 Regioisomers (X = F, Cl, Br, or OH)201233264164610.5012/bkcs.2012.33.2.641BeckeA. D.Density-functional thermochemistry. III. The role of exact exchange1993987564810.1063/1.464913LeeC.YangW.ParrR. G.Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density19883727857892-s2.0-034549110510.1103/PhysRevB.37.785StephenP. J.DevlinF. J.ChablowskiC. F.FrischM. J.Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields19949845116231162710.1021/j100096a001HerteingR. H.KochW.On the parameterization of the local correlation functional. What is Becke-3-LYP?19972685-634535110.1016/S0009-2614(97)00207-8HehreW. J.DitchfieldK.PopleJ. A.Self-consistent molecular orbital methods. XII. Further extensions of gaussian-type basis sets for use in molecular orbital studies of organic molecules1972565225722612-s2.0-034717000510.1063/1.1677527FrischM. J.TrucksG. W.SchlegelH. B.2003Pittsburgh, Pa, USAGaussian, Inc.FrischM. J.TrucksG. W.SchlegelH. B.2009Pittsburgh, Pa, USAGaussian, Inc.ClarkT.KochR.1999Berlin, GermanySpringer Berlin Heidelberg10.1007/978-3-662-13150-3ZhaoY.TruhlarD. G.The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals2008120121524110.1007/s00214-007-0310-xPerdewJ. P.BurkeK.ErnzerhofM.Generalized gradient approximation made simple [Phys. Rev. Lett. 77, 3865 (1996)]1997787139610.1103/PhysRevLett.78.1396SchmiderH. L.BeckeA. D.Optimized density functionals from the extended G2 test set199810823962410.1063/1.476438BeckA. D.Density-functional thermochemistry. V. Systematic optimization of exchange-correlation functionals199710720855410.1063/1.475007GrimmeS.EhrlichS.GoerigkL.Effect of the damping function in dispersion corrected density functional theory2011327145614652-s2.0-7995294355910.1002/jcc.21759Head-GordonM.PopleJ. A.FrischM. J.MP2 energy evaluation by direct methods198815365035062-s2.0-2504445521710.1016/0009-2614(88)85250-3WeigendF.AhlrichsR.Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: design and assessment of accuracy20057183297330510.1039/b508541aHudsonB. S.AllisD. G.ParkerS. F.Ramirez-CuestaA. J.HermanH.PrinzbachH.Infrared, Raman, and inelastic neutron scattering spectra of dodecahedrane: An I h molecule in T h site symmetry200510915341834242-s2.0-1814439083510.1021/jp0503213WahlF.WeilerA.LandenbergerP.SackersE.VossT.HaasA.LiebM.HunklerD.WorthJ.KnotheL.PrinzbachH.Towards perfunctionalized dodecahedranes-en route to C20 fullerene20061224625562672-s2.0-3374833572910.1002/chem.200501618