Synthesis and Structural Studies of Two Paramagnetic

May 24, 2019 - Copyright © 2019 American Chemical Society. RIGHTS & ... These studies illustrate that the electron spin can be used as a probe to ide...
3 downloads 0 Views 2MB Size
Download PDF
Article Cite This: Inorg. Chem. 2019, 58, 8162−8168

pubs.acs.org/IC

Synthesis and Structural Studies of Two Paramagnetic Metallofullerenes with Isomeric C72 Cage Chong Zhao,†,‡ Mingzhe Nie,†,‡ Haibing Meng,†,‡ Chunru Wang,*,† and Taishan Wang*,† †

Beijing National Laboratory for Molecular Sciences, Key Laboratory of Molecular Nanostructure and Nanotechnology, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China ‡ University of Chinese Academy of Sciences, Beijing 100049, China

Downloaded via UNIV OF SOUTHERN INDIANA on July 29, 2019 at 12:41:37 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: We synthesized and isolated two paramagnetic metallofullerenes of La@C72 and Y@C72 with different fullerene cages, which were characterized by electron paramagnetic resonance (EPR) spectroscopy and theoretical calculations. DFT calculations disclosed two possible isomers of La/Y@C72 with C72-C2 and C72-C2v cages, both of which have similar thermodynamic stability and one pair of fused pentagons. Their paramagnetic properties were then studied by EPR spectroscopy, and the obtained EPR signals were analyzed with very different hyperfine coupling constants, revealing distinct electron spin distributions for these two species. Furthermore, the experimental coupling constants were compared with those of calculated coupling constants, and comparison results revealed that the produced La@C72 has a C72-C2v cage and Y@C72 has a C72-C2 cage. These studies illustrate that the electron spin can be used as a probe to identify metallofullerene structure due to the susceptibility of spin-metal couplings. The successful isolation and characterizations of La@C72 and Y@C72 with such a small C72 cage reveal their stability that is important for application as paramagnetic molecule materials.

■

INTRODUCTION One of the intriguing properties of endohedral metallofullerene (EMF) is its paramagnetic character, which has potential in quantum information processing (QIP) and single molecular magnet (SMM) applications.1−7 The paramagnetic metallofullerene has an open-shell electronic structure, in which the unpaired spin is originated from the odd electron in the molecules.8−17 Spin-active metallofullerenes have been found to have high sensitivity to chemical modification, temperature, local magnetic field, host, and molecular rotation, which can greatly influence the spin distributions, spin relaxation, and spin couplings.18,19 Among them, the monometal endohedral fullerenes can provide valuable information about the nature of the single metal and have brought about distinct paramagnetic properties. Since the first discovery of monometal endohedral fullerene, La@C82,20 several La@C2n molecules and their derivatives have been synthesized. However, many pristine La@C2n monometal endohedral fullerenes have not been discovered due their high reactivity, and only their nonmagnetic derivatives can be isolated, such as the La@C72− C6H3Cl2.21 Likewise, many pristine Y@C2n have not been discovered except for [email protected]−26 Furthermore, due to the low yield of monometal endohedral fullerenes, it is difficult to identify their accurate structure. Therefore, it is meaningful to explore new method to characterize these paramagnetic monometal endohedral fullerenes. It is generally acknowledged that the couplings between the unpaired electron and magnetic nucleus always show energy splittings under an external magnetic field. The coupling parameters of A and g-factor values obtained by the electron © 2019 American Chemical Society

paramagnetic resonance (EPR) spectroscopy can be used to evaluate the chemical environment and kinetic information on electron spin as well as magnetic nucleus. Typically, the spinnucleus couplings in paramagnetic metallofullerenes provide related structural information. For example, for the paramagnetic Sc3C2@C80, the hyperfine couplings between spin and Sc3 indicate three equivalent Sc nuclei.1,4,27 Similar spinnucleus couplings are also found in Y2@C79N, and its EPR spectra illustrate the susceptible Y2 rotation inside the azafullerene C79N.3,9−11 Different Y2 and Y2C2 rotations within the same C82 cage were also observed by EPR spectroscopy.28 For the uncertain position of endohedral C and N for Sc3CN@ C80 in single crystal analysis, EPR spectroscopy can confirm them through analysis of hyperfine couplings.29 Therefore, the electron spin in metallofullerenes can be used as a probe to identify metallofullerene structure due to the susceptibility of spin-metal couplings. In this work, we successfully synthesized two monometal endohedral fullerenes La@C 72 and Y@C 72 by the Krätschmer−Huffman arc-discharging method and isolation by high-performance liquid chromatography (HPLC). Their structures and paramagnetic properties were studied by EPR spectroscopy combined with ab initio calculations. The experimental EPR coupling constants were compared with those of calculated coupling constants, and comparison results revealed that the produced La@C72 has a C72-C2v cage and Y@ C72 has a C72-C2 cage. These results show that different Received: April 2, 2019 Published: May 24, 2019 8162

DOI: 10.1021/acs.inorgchem.9b00954 Inorg. Chem. 2019, 58, 8162−8168

Article

Inorganic Chemistry

Figure 1. (a, b) HPLC separation profiles of La@C72 (20 × 250 mm Buckyprep column; flow rate of 6 mL/min; toluene as eluent). (c) The experimental MALDI-TOF mass spectra of La@C72. Insets show the experimental and simulated isotopic distributions. (d, e) HPLC separation profiles of Y@C72 (20 × 250 mm Buckyprep column; flow rate of 6 mL/min; toluene as eluent). (f) The experimental MALDI-TOF mass spectra of Y@C72. Insets show the experimental and simulated isotopic distributions. (g) UV−vis absorption spectra of Y@C72 and La@C72. (h−k) The optimized structures of La@C72-C2, La@C72-C2v, Y@C72-C2, and Y@C72-C2v. The yellow balls present the pair of fused pentagons. (Figure S3), and geometry optimizations are performed on the corresponding complex structures by B3LYP/6-31G* ∼ (lanl2dz for La and Y). The results are listed in Table S2 (the single point energy is by empirical dispersion correction).31 Frequency calculations were carried out to identify the nature of the optimized molecules. Rotational−vibrational partition functions are used to elaborate the relative concentration at the temperature span of 0−3500 K. The bonding features between La/Y and C72 were analyzed with a MO interaction diagram. The ESPs are mapped on the electron density isosurface (isovalue is 0.001 au). The total electron density is employed by a .fch file transferred by the formchk module in Gaussian 09. The ESP visualization is implemented through the Mutiwfn +VMD software. The above calculations were performed using the Gaussian 09 quantum chemical program package.32 The hfcc (hyperfine coupling constant) values of aiso are performed by the ORCA program at the level of BP86/Def-2 (La ∼ SARC-TZV for the DKH2 Hamiltonian).33 The EPR spectrum was simulated by easyspin codes (http://www.easyspin.org) used in MATLAB.

fullerene cages can greatly influence the spin distributions of paramagnetic metallofullerenes. Meanwhile, the EPR spectroscopy could be used as a tool to detect the structural characteristics.

■

EXPERIMENTAL SECTION

Synthesis of La@C72 and Y@C72. Metallofullerenes La@C72 and Y@C72 were prepared by the Krätschmer−Huffman arc discharge method.30 The as-prepared soot was Soxhlet-extracted with toluene for 24 h. The La@C72 and Y@C72 samples was isolated and purified by multistep high-performance liquid chromatography (HPLC), and the purity of the sample was confirmed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDITOF-MS). EPR Measurements. The EPR spectra were measured on a Bruker E500 instrument with a continuous-wave X band. The frequencies were 9.753 GHz for La@C72 and 9.859 GHz for Y@C72. Both spectra were obtained with the same power attenuation of 13.0 dB. All of the solution samples were dissolved in CS2 solution. Theoretical Calculations. The coordinates of C72 were generated by the revised CaGe software. There are 11190 classic isomers for C72. We selected the 42 isomers with zero to three B55 bonds (B55 means the shared C−C bond of fused pentagons) as the candidate cages. All the isomers of the selected C72 and their anions (C722−, C723−) are optimized at the PM3 and B3LYP/3-21G* level in turn. The results are listed in Table S1. Since the encaged cluster would transfer electrons to the parent cages, the stability of the anions is closely related to the corresponding metallofullerenes; the lowest energy isomers of anions are selected as candidate cages for encaging the La/ Y atom. La/Y is then put inside the cages at the typical seven sites

■

RESULTS AND DISCUSSION Structural Characterization. La@C72 and Y@C72 were both synthesized by the arc-discharging method under 200 Torr He. The produced fullerenes and metallofullerenes were isolated by multistep HPLC as shown in Figure 1a, b and Figure 1d, e. The retention time of target La@C72 or Y@C72 is similar to C60O and La@C74 on the Buckyprep column. Their MALDI-TOF mass spectra are shown in Figure 1c and 1f. The highest intensity of isotopic abundance ratio appears at 1003.91 (m/z) for La@C72 and at 953.91 (m/z) for Y@C72, 8163

DOI: 10.1021/acs.inorgchem.9b00954 Inorg. Chem. 2019, 58, 8162−8168

Article

Inorganic Chemistry

Table 1. Calculated Relative Energies,a Distance of La/Y−C,b and EEc of the Optimized Isomers of La@C72-C2/C2v and Y@ C72-C2/C2v isomer

△E-B3LYP/TPSSh

d(La−C)-B3LYP/TPSSh

EE-B3LYP/TPSSh

La@C72-C2 La@C72-C2v isomer

0 0.17/0.67 △E-B3LYP/TPSSh

2.64/2.62 2.59/2.58 d(Y−C)-B3LYP/TPSSh

3.43/5.75 3.00/5.17 EE-B3LYP/TPSSh

Y@C72-C2 Y@C72-C2v

0 0.69/1.28

2.39/2.38 2.37/2.36

3.53/5.55 2.95/5.24

a

Unit in kilocalories per mole (kcal/mol). bThe C represents one of the B55 bond atoms, with the unit in Angstroms (Å). cEE stands for encapsulation energy, with the unit in electronvolts (eV).

Figure 2. Relative concentrations of (a) La@C72 and (b) Y@C72 isomers computed with the TPSSh/6-31G* ∼ lanl2dz level using the RRHO approximation at a temperature in the range of 0−3500 K.

Figure 3. (a, b) Experimental and simulated EPR spectra of La@C72 and Y@C72 in CS2 solution at 298 K. (c−f) Spin density distributions of La@ C72-C2, La@C72-C2v, Y@C72-C2, and Y@C72-C2v. All isovalues of spin density distributions are the same value of 0.003.

C2v, La/Y@C72-2−2, La/Y@C72-2−3, and La/Y@C72-2−8 in turn at the B3LYP/6-31G* ∼ (lanl2dz for La and Y) level, whereas the energy differences change slightly at the B3LYP/631G* and TPSSh/6-31G* levels, as shown in Table 1. It is worth noting that the most stable isomers of C72-C2 and C72C2v are with small energy differences of 0.17/0.67 kcal/mol for La@C72 and 0.69/1.28 kcal/mol for Y@C72. In the previous report for La@C72,34 the C2 one is 0.26 kcal/mol more than the C2v one with a small difference at the B3LYP/6-31G* level. In fact, the two isomers show similar characteristics; for example, both of them have one pair of fused pentagons, similar La/Y sites as the optimized structures in Figure 1h−k

which are in good agreement with the simulated isotopic distribution patterns. In addition, the UV−vis absorption spectra were also characterized (Figure 1g), and the results show that there are obvious differences between the two spectra, indicating the different isomeric structures of prepared Y@C72 and La@C72. To elucidate the structures of La@C72 and Y@C72, DFT calculations were employed to select and optimize their isomers. The calculated relative energies of the first ten isomers of C72, C722−, and C723− at the B3LYP/3-21G* level are provided in Table S1. Table S2 shows that the five lowestenergy isomers of La/Y@C72 are La/Y@C72-C2, La/Y@C728164

DOI: 10.1021/acs.inorgchem.9b00954 Inorg. Chem. 2019, 58, 8162−8168

Article

Inorganic Chemistry

C72; therefore, the separated out La@C72 is determined as La@C72-C2v, which is different from reported La@C72− C6H3Cl2. For Y@C72, the results are the opposite. The calculated a(Y)iso for Y@C72-C2 and Y@C72-C2v are 0.75/1.31/ 0.95 G and 0.03/0.09/0.01 G (Table 2) at the three levels, respectively. The experimental hyperfine coupling constant of Y value (1.75 G) is closer to the calculated constant of Y@C72C2 but not that of Y@C72-C2v; therefore, the separated out Y@ C72 is determined as Y@C72-C2. Considering the natural paramagnetic property of La/Y@C72 and the sensitive hyperfine coupling splitting, the EPR results can be used as the “fingerprint” of the paramagnetic metallofullerenes. Figure 3c and d shows the spin density distributions of the La@C72 and Y@C72 four isomers. The spin of La@C72-C2 is mainly polarized on the hemisphere near the fused pentagons, whereas the spin of La@C72-C2v is polarized on the opposite hemisphere. Much the same thing was found for Y@C72, in that the spin of Y@C72-C2 is mainly polarized on the hemisphere near the fused pentagons, whereas the spin of Y@C72-C2v is polarized on the opposite hemisphere (Figure 3e, f). These results indicate that the hyperfine coupling values for a(La/Y) of La/Y@C72-C2 are bigger than that of La/Y@C72C2v, which coincides with the calculated hyperfine coupling constant values and experimental EPR results. MO Properties. For La@C72-C2v, Figure 4a displays the orbital interaction diagram of C72-C2v, La@C72-C2v, and La atom. The MOs of the La atom show the valence electrons of 5d16s2 orbits, where the upper two energy levels (green) show the 6s2 orbits (occupy the α and β orbits with the same energy of −4.11 eV) and the lower line shows the 5d1 orbit (occupy the α orbit), for La atom with the energy of −4.23 eV (red). Because the frontier MO energy levels of La are higher than that of the corresponding parent cage C72-C2v (−5.16 eV), the doublet La is prone to transfer electrons to the unoccupied molecular orbitals of the parent C72-C2v cage. Meanwhile, orbital hybridization between the cluster and the cage may result in the reduction of the energy level. Figure 4a clearly shows that metal La formally transfers 3-fold electrons to the cage which is marked by red and green lines from La. For Y@ C72-C2, the orbital interaction diagram of C72-C2, Y@C72-C2, and Y atom is also provided in Figure 4b. The MOs of Y atom show the valence electrons of 4d15s2 orbits, with the lower two (green) levels showing the 5s2 orbits with the similar energy of −4.64 eV (β orbit) and −4.74 eV (α orbit) and the upper energy level showing the 4d1 orbit (occupy the α orbit) with the energy of −4.36 eV (red). The frontier MO energy levels of Y are all higher than the corresponding parent cage C72-C2 (−5.10 eV); therefore, the doublet Y is prone to transfer three electrons to the unoccupied molecular orbitals of the parent C72-C2 cage like the La@C72. Furthermore, the pictures in Figure 4a and b show that the SOMOs of La@C72-C2v and Y@ C72-C2 are consistent with the spin density distributions as shown in Figure 3d and e, respectively. In addition, the surface ESPs (electrostatic potential) of La@C72-C2v and (C72-C2v)3− (Figure 4c) and of Y@C72-C2 and (C72-C2)3− (Figure 4d) show distinctly different distributions and values; although the ion model is widely accepted, the introduction of metal clusters leads to completely different surface electrostatic potentials as shown. In detail, the ESPs of the charged (C72-C2v/C2)3− cage are much more negative than that of EMFs, which is similar to the reported ESPs for C726− and La [email protected] Besides, the fused pentagons (metalcoordinated carbon atoms) are more negative than the isolated

(the bonds of La and the near B55 for La@C72 are 2.64 and 2.59 Å, respectively, as shown in Table 1, which is about 0.22− 0.25 Å longer than those of d(Y−C) for Y@C72 at the same method of B3LYP), and similar encapsulation energies (EEs). Therefore, it is difficult to distinguish the isomers by these indistinguishable indexes. Temperature-Dependent Stability. Relative concentrations (mole fractions) can be evaluated through their partition functions qi and the enthalpies at the absolute zero temperature and ground-state energies △H (i.e., the relative potential energies corrected for the vibrational zero-point energies) by the formula as reference (rotational−vibrational partition functions were using the rigid rotator and harmonic oscillator (RRHO) approximation).35 The relative concentrations (yields) of different isomers depend on temperature since they are formed at a very high temperature. The relative yields of the five lowest-energy isomers are examined to judge the temperature effect (Figure S4). Figure 2a shows the relative concentration of the two lowest-energy La@C72 isomers from 0 to 3500 K. The lowest-energy isomer La@ C72-C2 is prevailing at the temperatures exceeding 284 K, whereas La@C72-C2v has a slight advantage at the temperatures below 284 K. It is clear that La@C72-C2v and La@C72-C2 may coexist throughout the whole temperature range, which is similar to that as previously reported for Ca@C72;35 that is, the La@C72-C2 and La@C72-C2v are with similar thermodynamic stability, and these results also coincide with the reported calculation in ref 34. For Y@C72, the lowest-energy isomer of Y@C72-C2 is prevailing throughout the whole temperature region (Figure 2b), indicating that the C2 cage is more stable than the C2v cage. EPR Analysis. Experimentally, EPR spectroscopy characterizations of produced La@C72 and Y@C72 were performed. The EPR spectrum of La@C72 (Figure 3a) shows eight clear splitting peaks; the isotropic hyperfine coupling constant of a(La)iso is 0.16 G, and the g-factor is 2.0018. Among reported mono La-based metallofullerenes, this a(La)iso value is the smallest one,17 and the g value is similar to those of La@C82 derivatives36 and La@C86−II.37 The EPR spectrum of Y@C72 (Figure 3b) shows two clear splitting peaks; the isotropic hyperfine coupling constant of a(Y)iso is 1.75 G, and the gfactor is 2.0015. Moreover, we calculated the hyperfine coupling constants (the isotropic part, dipolar part, and second-order contribution part from SOC (spin orbit coupling) are considered) at the level of BP86/B3LYP/TPSS ∼ Def-2 with the ORCA program. For La@C72, the calculation results show that a(La)iso for La@C72-C2 and La@C72-C2v are 2.02/2.16/2.04 G and 0.13/0.17/0.09 G (Table 2) at the three levels, respectively. Obviously, the calculated hyperfine coupling constants of La@C72-C2v isomer are consistent with the experimental a(La)iso value (0.16 G) for produced La@ Table 2. DFT Calculated Hyperfine Coupling Constant Values of a(La) and a(Y) by the BP86, B3LYP, and TPSS Methods isomer

a(La)-BP86

a(La)-B3LYP

a(La)-TPSS

La@C72-C2 La@C72-C2v isomer

2.02 G 0.13 G a(Y)-BP86

2.16 G 0.17 G a(Y)-B3LYP

2.04 G 0.09 G a(Y)-TPSS

Y@C72-C2 Y@C72-C2v

0.75 G 0.03 G

1.31 G 0.09 G

0.95 G 0.01 G 8165

DOI: 10.1021/acs.inorgchem.9b00954 Inorg. Chem. 2019, 58, 8162−8168

Article

Inorganic Chemistry

Figure 4. (a) Molecular orbital interaction diagrams for C72-C2v, La@C72-C2v, and La. (b) Molecular orbital interaction diagram for C72-C2, Y@C72C2, and Y. The molecular orbitals of empty cages C72-C2v and C72-C2 in (a) and (b) are calculated with the fragments of the corresponding endohedral fullerenes without further optimization. (c) The calculated ESP for the optimized La@C72-C2v and (C72-C2v)3−. (d) The calculated ESP for the optimized Y@C72-C2 and (C72-C2)3−. Detailed values of ESPs are given in the Supporting Information.

■

pentagons fragments (IPR fragments) for both the EMFs and the −3 anion cages. In conclusion, the distribution of pentagons determines the relative intensity of the electrostatic potential, while the charge number of the carbon cage determines the absolute intensity of the surface electrostatic potential.

■

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.inorgchem.9b00954.

■

CONCLUSION

Two paramagnetic metallofullerenes La@C72 and Y@C72 have been discovered. Although monometal endohedral fullerenes have a high reactivity, the successful isolation and characterizations of La@C72 and Y@C72 reveal their stability even with such a small C72 cage. The two most advantageous isomers of La/Y@C72-C2 and La/Y@C72-C2v are with similar thermodynamics stability via DFT calculations. However, their different EPR coupling constants revealed varied electron spin distributions and structures for these two species. The electron spin in paramagnetic metallofullerenes provides much structural information as the spin-metal couplings are sensitive to the subtle structural differences. The experimental EPR coupling constants were compared with those of calculated coupling constants, and comparison results revealed that the produced La@C72 has a C72-C2v cage and Y@C72 has a C72-C2 cage. The MO analysis further proves the ionic model as (La)3+@(C72)3− and (Y)3+@(C72)3−. Therefore, this study further confirms that the electron spin in metallofullerenes can be used as a probe to identify metallofullerene structures. The successful isolation and characterizations of La@C72 and Y@ C72 reveal their stability, which is important for their application as paramagnetic molecule materials.

Experimental details, HPLC chromatogram, and calculation details (PDF)

AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected] (T.W.). *E-mail: [email protected] (C.W.). ORCID

Haibing Meng: 0000-0002-9538-3352 Chunru Wang: 0000-0001-7984-6639 Taishan Wang: 0000-0003-1834-3610 Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China (51672281, 61227902, 51832008). T.W. particularly thanks the Youth Innovation Promotion Association of CAS (2015025).

■

REFERENCES

(1) Wang, T.; Wu, J.; Xu, W.; Xiang, J.; Lu, X.; Li, B.; Jiang, L.; Shu, C.; Wang, C. Spin divergence induced by exohedral modification: ESR study of Sc3C2@C80 fulleropyrrolidine. Angew. Chem., Int. Ed. 2010, 49, 1786−1789. (2) Plant, S. R.; Jevric, M.; Morton, J. J. L.; Ardavan, A.; Khlobystov, A. N.; Briggs, G. A. D.; Porfyrakis, K. A two-step approach to the

8166

DOI: 10.1021/acs.inorgchem.9b00954 Inorg. Chem. 2019, 58, 8162−8168

Article

Inorganic Chemistry synthesis of N@C60 fullerene dimers for molecular qubits. Chem. Sci. 2013, 4, 2971−2975. (3) Feng, Y. Q.; Wang, T. S.; Li, Y. J.; Li, J.; Wu, J. Y.; Wu, B.; Jiang, L.; Wang, C. R. Steering Metallofullerene Electron Spin in Porous Metal Organic Framework. J. Am. Chem. Soc. 2015, 137, 15055− 15060. (4) Wu, B.; Wang, T. S.; Feng, Y. Q.; Zhang, Z. X.; Jiang, L.; Wang, C. R. Molecular magnetic switch for a metallofullerene. Nat. Commun. 2015, 6, 6468. (5) Liu, Z.; Dong, B.-W.; Meng, H.-B.; Xu, M.-X.; Wang, T.-S.; Wang, B.-W.; Wang, C.-R.; Jiang, S.-D.; Gao, S. Qubit crossover in the endohedral fullerene Sc3C2@C80. Chem. Sci. 2018, 9, 457−462. (6) Meng, H.; Zhao, C.; Li, Y.; Nie, M.; Wang, C.; Wang, T. An implanted paramagnetic metallofullerene probe within a metal-organic framework. Nanoscale 2018, 10, 3291−3298. (7) Ungur, L.; Chibotaru, L. F. Strategies toward High-Temperature Lanthanide-Based Single-Molecule Magnets. Inorg. Chem. 2016, 55, 10043−10056. (8) Hu, S.; Liu, T.; Shen, W.; Slanina, Z.; Akasaka, T.; Xie, Y.; Uhlik, F.; Huang, W.; Lu, X. Isolation and Structural Characterization of Er@C2v(9)-C82 and Er@Cs(6)-C82: Regioselective Dimerization of a Pristine Endohedral Metallofullerene Induced by Cage Symmetry. Inorg. Chem. 2019, 58, 2177−2182. (9) Zuo, T. M.; Xu, L. S.; Beavers, C. M.; Olmstead, M. M.; Fu, W. J.; Crawford, D.; Balch, A. L.; Dorn, H. C. M-2@C79N (M = Y, Tb): Isolation and characterization of stable endohedral metallofullerenes exhibiting M-M bonding interactions inside aza[80]fullerene cages. J. Am. Chem. Soc. 2008, 130, 12992−12997. (10) Ma, Y. H.; Wang, T. S.; Wu, J. Y.; Feng, Y. Q.; Jiang, L.; Shu, C. Y.; Wang, C. R. Susceptible electron spin adhering to an yttrium cluster inside an azafullerene C79N. Chem. Commun. 2012, 48, 11570−11572. (11) Zhao, C.; Meng, H.; Nie, M.; Jiang, L.; Wang, C.; Wang, T. Anisotropic Paramagnetic Properties of Metallofullerene Confined in a Metal−Organic Framework. J. Phys. Chem. C 2018, 122, 4635− 4640. (12) Fu, W.; Zhang, J.; Fuhrer, T.; Champion, H.; Furukawa, K.; Kato, T.; Mahaney, J. E.; Burke, B. G.; Williams, K. A.; Walker, K.; Dixon, C.; Ge, J.; Shu, C.; Harich, K.; Dorn, H. C. Gd2@C79N: isolation, characterization, and monoadduct formation of a very stable heterofullerene with a magnetic spin state of S = 15/2. J. Am. Chem. Soc. 2011, 133, 9741−9750. (13) Hu, Z.; Dong, B. W.; Liu, Z.; Liu, J. J.; Su, J.; Yu, C.; Xiong, J.; Shi, D. E.; Wang, Y.; Wang, B. W.; Ardavan, A.; Shi, Z.; Jiang, S. D.; Gao, S. Endohedral Metallofullerene as Molecular High Spin Qubit: Diverse Rabi Cycles in Gd2@C79N. J. Am. Chem. Soc. 2018, 140, 1123−1130. (14) Cimpoesu, F.; Frecus, B.; Oprea, C. I.; Ramanantoanina, H.; Urland, W.; Daul, C. On exchange coupling and bonding in the Gd2@C80 and Gd2@C79N endohedral dimetallo-fullerenes. Mol. Phys. 2015, 113, 1712−1727. (15) Zhang, Z.; Wang, T.; Xu, B.; Wang, C. Paramagnetic and theoretical study of Y(2)@C(8)(1)N: an endohedral azafullerene radical. Dalton Trans. 2014, 43, 12871−12875. (16) Hummelen, J. C.; Knight, B.; Pavlovich, J.; Gonzalez, R.; Wudl, F. Isolation of the Heterofullerene C59N as Its Dimer (C59N)2. Science 1995, 269, 1554−1556. (17) Popov, A. A.; Yang, S.; Dunsch, L. Endohedral fullerenes. Chem. Rev. 2013, 113, 5989−6113. (18) Wang, T.; Wang, C. Endohedral metallofullerenes based on spherical I(h)-C(80) cage: molecular structures and paramagnetic properties. Acc. Chem. Res. 2014, 47, 450−458. (19) Popov, A. A. Endohedral Fullerenes: Electron Transfer and Spin 2017, 1. (20) Akasaka, T.; Wakahara, T.; Nagase, S.; Kobayashi, K.; Waelchli, M.; Yamamoto, K.; Kondo, M.; Shirakura, S.; Okubo, S.; Maeda, Y.; Kato, T.; Kako, M.; Nakadaira, Y.; Nagahata, R.; Gao, X.; Van Caemelbecke, E.; Kadish, K. M. La@C82Anion. An Unusually Stable Metallofullerene. J. Am. Chem. Soc. 2000, 122, 9316−9317.

(21) Wakahara, T.; Nikawa, H.; Kikuchi, T.; Nakahodo, T.; Rahman, G. M.; Tsuchiya, T.; Maeda, Y.; Akasaka, T.; Yoza, K.; Horn, E.; Yamamoto, K.; Mizorogi, N.; Slanina, Z.; Nagase, S. La@C72 having a non-IPR carbon cage. J. Am. Chem. Soc. 2006, 128, 14228−14229. (22) Suzuki, S.; Kawata, S.; Shiromaru, H.; Yamauchi, K.; Kikuchi, K.; Kato, T.; Achiba, Y. Isomers and carbon-13 hyperfine structures of metal-encapsulated fullerenes M@C82 (M = Sc, Y, and La). J. Phys. Chem. 1992, 96, 7159−7161. (23) Lu, X.; Nikawa, H.; Feng, L.; Tsuchiya, T.; Maeda, Y.; Akasaka, T.; Mizorogi, N.; Slanina, Z.; Nagase, S. Location of the yttrium atom in Y@C82 and its influence on the reactivity of cage carbons. J. Am. Chem. Soc. 2009, 131, 12066−12067. (24) Shinohara, H.; Sato, H.; Saito, Y.; Ohkohchi, M.; Ando, Y. Mass spectroscopic and ESR characterization of soluble yttrium-containing metallofullerenes YC82 and Y2C82. J. Phys. Chem. 1992, 96, 3571− 3573. (25) Misochko, E. Y.; Akimov, A. V.; Belov, V. A.; Tyurin, D. A.; Bubnov, V. P.; Kareev, I. E.; Yagubskii, E. B. EPR spectrum of the Y@ C82 metallofullerene isolated in solid argon matrix: hyperfine structure from EPR spectroscopy and relativistic DFT calculations. Phys. Chem. Chem. Phys. 2010, 12, 8863−8869. (26) Yamada, M.; Feng, L.; Wakahara, T.; Tsuchiya, T.; Maeda, Y.; Lian, Y.; Kako, M.; Akasaka, T.; Kato, T.; Kobayashi, K.; Nagase, S. Synthesis and characterization of exohedrally silylated M@C82 (M = Y and La). J. Phys. Chem. B 2005, 109, 6049−6051. (27) Meng, H.; Zhao, C.; Nie, M.; Wang, C.; Wang, T. Optically Controlled Molecular Metallofullerene Magnetism via an Azobenzene-Functionalized Metal-Organic Framework. ACS Appl. Mater. Interfaces 2018, 10, 32607−32612. (28) Ma, Y.; Wang, T.; Wu, J.; Feng, Y.; Li, H.; Jiang, L.; Shu, C.; Wang, C. Electron Spin Manipulation via Encaged Cluster: Differing Anion Radicals of Y2@C82-Cs, Y2C2@C82-Cs, and Sc2C2@C82Cs. J. Phys. Chem. Lett. 2013, 4, 464−467. (29) Wang, T. S.; Feng, L.; Wu, J. Y.; Xu, W.; Xiang, J. F.; Tan, K.; Ma, Y. H.; Zheng, J. P.; Jiang, L.; Lu, X.; Shu, C. Y.; Wang, C. R. Planar quinary cluster inside a fullerene cage: synthesis and structural characterizations of Sc(3)NC@C(80)-I(h). J. Am. Chem. Soc. 2010, 132, 16362−16364. (30) Kroto, H. W.; Heath, J. R.; O’Brien, S. C.; Curl, R. F.; Smalley, R. E. C60: Buckminsterfullerene. Nature 1985, 318, 162−163. (31) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (32) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Keith, T.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian 09, revision A.02; Gaussian, Inc.: Wallingford, CT, 2016. (33) Pantazis, D. A.; Neese, F. All-Electron Scalar Relativistic Basis Sets for the Lanthanides. J. Chem. Theory Comput. 2009, 5, 2229− 2238. (34) Slanina, Z.; Uhlik, F.; Lee, S. L.; Adamowicz, L.; Akasaka, T.; Nagase, S. Stability computations for isomers of La@C(n) (n = 72, 74, 76). Molecules 2012, 17, 13146−13156. (35) Slanina, Z.; Lee, S.-L.; Uhlík, F.; Adamowicz, L.; Nagase, S. Computing relative stabilities of metallofullerenes by Gibbs energy treatments. Theor. Chem. Acc. 2007, 117, 315−322. 8167

DOI: 10.1021/acs.inorgchem.9b00954 Inorg. Chem. 2019, 58, 8162−8168

Article

Inorganic Chemistry (36) Maeda, Y.; Sato, S.; Inada, K.; Nikawa, H.; Yamada, M.; Mizorogi, N.; Hasegawa, T.; Tsuchiya, T.; Akasaka, T.; Kato, T.; Slanina, Z.; Nagase, S. Regioselective exohedral functionalization of La@C82 and its 1,2,3,4,5-pentamethylcyclopentadiene and adamantylidene adducts. Chem. - Eur. J. 2010, 16, 2193−2197. (37) Okubo, S.; Kato, T. ESR parameters of a series of La@C n isomers. Appl. Magn. Reson. 2003, 23, 481−493. (38) Popov, A. A.; Kästner, C.; Krause, M.; Dunsch, L. Carbon Cage Vibrations of M@C82and M2@C2n(M = La, Ce; 2n = 72, 78, 80): The Role of the Metal Atoms. Fullerenes, Nanotubes, Carbon Nanostruct. 2014, 22, 202−214.

8168

DOI: 10.1021/acs.inorgchem.9b00954 Inorg. Chem. 2019, 58, 8162−8168