Understanding the Reactivity of Endohedral Metallofullerenes: C78 versus Sc3N@C78

The physical factors behind the reduced Diels–Alder reactivity of the Sc₃N@C₇₈ metallofullerene as compared with free C₇₈ have been investigated in detail by means of computational tools. To this end, the reactions between 1,3‐butadiene and free C₇₈ and endohedral Sc₃N@C₇₈ have been analysed in terms of regioselectivity and reactivity by using the activation strain model of reactivity in combination with the energy decomposition analysis method. Additional factors such as the molecular orbital overlap or the aromaticity of the corresponding transition states have been also explored. Our results indicate that the lower reactivity of the metallofullerene finds its origin mainly in the less stabilizing interaction between the deformed reactants along the reaction coordinate induced by the triscandium nitride moiety.     

Endohedral Metallofullerenes: New Structures and Unseen Phenomena

Endohedral metallofullerenes (EMFs), namely fullerenes with metallic species encapsulated inside, represent an ideal platform to investigate metal–metal or metal–carbon interactions at the sub-nanometer scale by means of single-crystal X-ray diffraction (XRD) crystallography. Herein, recent progress in the identification of new structures and unprecedented properties are discussed according to the categories of monometallofullerenes, dimetallofullerenes, carbide clusterfullerenes, and nitride clusterfullerenes. In particular, the dimerization and the cage-isomer dependent oxidation state of the inner metal atom are summarized in terms of pristine monometallofullerenes. Metal–metal bonds involving lanthanide–lanthanides or actinide–actinides are discussed based on both experimental and theoretical studies. The cluster–cage matching and/or mutual selections, as well as the rarely seen M=C double bonds, are discovered in M₂C₂@C_2n, U₂C@C₈₀, M₂TiC@C₈₀, and Ti₃C₃@C₈₀. Subsequently, the geometries of different M₃N clusters in various cages are discussed, revealing size-matching between the internal M₃N cluster and the outer cage induced by the planarity of the cluster. Finally, an outlook regarding the future developments of the molecular structures and applications of EMFs is presented.     

A Complete Guide on the Influence of Metal Clusters in the Diels–Alder Regioselectivity of Ih‐C80 Endohedral Metallofullerenes

The chemical functionalization of endohedral metallofullerenes (EMFs) has aroused considerable interest due to the possibility of synthesizing new species with potential applications in materials science and medicine. Experimental and theoretical studies on the reactivity of endohedral metallofullerenes are scarce. To improve our understanding of the endohedral metallofullerene reactivity, we have systematically studied with DFT methods the Diels–Alder cycloaddition between s‐cis‐1,3‐butadiene and practically all X@I_h‐C₈₀ EMFs synthesized to date: X=Sc₃N, Lu₃N, Y₃N, La₂, Y₃, Sc₃C₂, Sc₄C₂, Sc₃CH, Sc₃NC, Sc₄O₂ and Sc₄O₃. We have studied both the thermodynamic and kinetic regioselectivity, taking into account the free rotation of the metallic cluster inside the fullerene. This systematic study has been made possible through the use of the frozen cage model (FCM), a computationally cheap approach to accurately predicting the exohedral regioselectivity of cycloaddition reactions in EMFs. Our results show that the EMFs are less reactive than the hollow I_h‐C₈₀ cage. Except for the Y₃ cluster, the additions occur predominantly at the [5,6] bond. In many cases, however, a mixture of the two possible regioisomers is predicted. In general, [6,6] addition is favored in EMFs that have a larger charge transfer from the metal cluster to the cage or a voluminous metal cluster inside. The present guide represents the first complete and exhaustive investigation of the reactivity of I_h‐C₈₀‐based EMFs.     

Endohedral Hydrogen Bonding Templates the Formation of a Highly Strained Covalent Organic Cage Compound**

A highly strained covalent organic cage compound was synthesized from hexahydroxy tribenzotriquinacene (TBTQ) and a meta-terphenyl-based diboronic acid with an additional benzoic acid substituent in 2’-position. Usually, a 120° bite angle in the unsubstituted ditopic linker favors the formation of a [4+6] cage assembly. Here, the introduction of the benzoic acid group is shown to lead to a perfectly preorganized circular hydrogen-bonding array in the cavity of a trigonal-bipyramidal [2+3] cage, which energetically overcompensates the additional strain energy caused by the larger mismatch in bite angles for the smaller assembly. The strained cage compound was analyzed by mass spectrometry and ¹H, ¹³C and DOSY NMR spectroscopy. DFT calculations revealed the energetic contribution of the hydrogen-bonding template to the cage stability. Furthermore, molecular dynamics simulations on early intermediates indicate an additional kinetic effect, as hydrogen bonding also preorganizes and rigidifies small oligomers to facilitate the exclusive formation of smaller and more strained macrocycles and cages.     

The Regioselectivity of Bingel–Hirsch Cycloadditions on Isolated Pentagon Rule Endohedral Metallofullerenes

In this work, the Bingel–Hirsch addition of diethylbromomalonate to all non‐equivalent bonds of Sc₃N@D_3h‐C₇₈ was studied using density functional theory calculations. The regioselectivities observed computationally allowed the proposal of a set of rules, the predictive aromaticity criteria (PAC), to identify the most reactive bonds of a given endohedral metallofullerene based on a simple evaluation of the cage structure. The predictions based on the PAC are fully confirmed by both the computational and experimental exploration of the Bingel–Hirsch reaction of Sc₃N@D_5h‐C_80, thus indicating that these rules are rather general and applicable to other isolated pentagon rule endohedral metallofullerenes.     

Bingel–Hirsch Addition on Endohedral Metallofullerenes: Kinetic Versus Thermodynamic Control

An extensive theoretical study of the Bingel–Hirsch addition of bromomalonate on scandium nitride endohedral fullerenes has been carried out. The prototypical and highly symmetrical Sc₃N@I_h‐C₈₀, with a structure that satisfies the isolated pentagon rule (IPR), and the non‐IPR Sc₃N@D₃(6140)‐C₆₈ fullerene show analogous reaction paths despite the distinct topology of the carbon networks and different rotation freedom of the internal nitride cluster. For the two metallofullerenes, our results predict that the reaction takes place under kinetic control yielding open‐cage fulleroids on [6,6] bonds, which is in good agreement with experimental data. The theoretical studies also show that predicting the reactivity of endohedral metallofullerenes is not straightforward and often an accurate analysis of the potential energy surface is required.     

Contrast of Bonding and Characters for C6, B3N3, C6H6 and B3N3H6 Molecules

Through integrative consideration of NICS, MO, MOC and NBO, we precisely investigated delocalization and bonding characters of C₆, C₆H₆, B₃N₃ and B₃N₃H₆ molecules. Firstly, we originally discovered and testified that C₆ cluster was sp² hybridization. Negative NICS values in 0 and 1 Å indicated that C₆ had δ and Π aromaticity. Secondly, B₃N₃ with sp² hybridization had obvious δ aromaticity. Finally, WBI values approved that there were delocalization in C₆, C₆H₆ and B₃N₃ molecules, but B₃N₃H₆ structure did not have delocalization with the WBI 1.0. Moreover, total WBI values of carbon, boron and nitrogen atoms were four, three and three, respectively. Namely, the electrons of B₃N₃H₆ and B₃N₃ were localized in nitrogen atoms and they did not form delocalized bonding. In a word, bonding characters of carbon, boron and nitrogen atoms were dissimilar although the molecules composed of carbon, boron and nitrogen were regarded as isoelectronic structures.     

C68 Fullerene Isomers,Anions, and their Metallofullerenes: Charge‐Stabilizing Different Isomers

The complete set of 6332 classical isomers of the fullerene C₆₈ as well as several non‐classical isomers is investigated by PM3, and the data for some of the more stable isomers are refined by the DFT‐based methods HCTH and B3LYP. C₂:0112 possesses the lowest energy of all the neutral isomers and it prevails in a wide range of temperatures. Among the fullerene ions modeled, C₆₈^2?, C₆₈^4? and C₆₈^6?, the isomers C₆₈^2?(C_s:0064), C₆₈^4?(C_2v:0008), and C₆₈^6?(D₃:0009) respectively, are predicted to be the most stable. This reveals that the pentagon adjacency penalty rule (PAPR) does not necessarily apply to the charged fullerene cages. The vertical electron affinities of the neutral C_s:0064, C_2v:0008, and D₃:0009 isomers are 3.41, 3.29, and 3.10 eV, respectively, suggesting that they are good electron acceptors. The predicted complexation energy, that is, the adiabatic binding energy between the cage and encapsulated cluster, of Sc₂C₂@C₆₈(C_2v:0008) is ?6.95 eV, thus greatly releasing the strain of its parent fullerene (C_2v:0008). Essentially, C₆₈ fullerene isomers are charge‐stabilized. Thus, inducing charge facilitates the isolation of the different isomers. Further investigations show that the steric effect of the encaged cluster should also be an important factor to stabilize the C₆₈ fullerenes effectively.     

Packed to the Rafters: Filling up C60 with Rare Gas Atoms

How many rare gas atoms can be placed into a fullerene cage until the pressure becomes large enough to break the C₆₀ framework? The answer given by density functional and ab initio computations is surprising and underlines the high stability of this unique carbon structure.     

Boosting Activation of Oxygen Molecules on C60 Fullerene by Boron Doping

The activation of oxygen molecules on boron‐doped C₆₀ fullerene (C₅₉B) and the subsequent water formation reaction are systematically investigated by using hybrid density functional calculations. Results indicate that C₅₉B shows a favorable ability to activate oxygen molecules both kinetically and thermodynamically. The oxygen molecule is first adsorbed on the boron atom, which is identified to be the most reactive site in C₅₉B for O₂ adsorption because of its high positive charge and spin density. The adsorption structure C₅₉B?O₂ can further isomerize to form two products with small reaction barriers. Water formation reactions upon these two structures are energetically favorable and suggest a four‐electron mechanism for the oxygen reduction reaction catalyzed by C₅₉B. This work provides a reliable theoretical insight into the catalytic properties of boron‐doped fullerene, which is believed to be helpful to explore fullerene catalysts.     

Is This a Chemical Bond? A Theoretical Study of Ng2@C60 (Ng=He,Ne, Ar,Kr, Xe)

Quantum-chemical calculations using DFT (BP86) and ab initio methods (MP2, SCS-MP2) have been carried out for the endohedral fullerenes Ng2@C60 (Ng=He-Xe). The nature of the interactions has been analyzed with charge- and energy-partitioning methods and with the topological analysis of the electron density (Atoms-in-Molecules (AIM)). The calculations predict that the equilibrium geometries of Ng2@C60 have D3d symmetry when Ng=Ne, Ar, Kr, while the energy-minimum structure of Xe2@C60 has D5d symmetry. The precession movement of He2 in He2@C60 has practically no barrier. The Ng--Ng distances in Ng2@C60 are much shorter than in free Ng2. All compounds Ng2@C60 are thermodynamically unstable towards loss of the noble gas atoms. The heavier species Ar2@C60, Kr2@C60, and Xe2@C60 are high energy compounds which are at the BSSE corrected SCS-MP2/TZVPP level in the range 96.7-305.5 kcal mol(-1) less stable than free C60+2 Ng. The AIM method reveals that there is always an Ng--Ng bond path in Ng2@C60. There are six Ng--C bond paths in (D3d) Ar2@C60, Kr2@C60, and Xe2@C60, whereas the lighter D3d homologues He2@C60 and Ne2@C60 have only three Ng--C2 paths. The calculated charge distribution and the orbital analysis clearly show that the bonding situation in Xe2@C60 significantly differs from those of the lighter homologues. The atomic partial charge of the [Xe2] moiety is +1.06, whereas the charges of the lighter dimers [Ng2] are close to zero. The a2u HOMO of (D3d) Xe2@C60 in the 1A1g state shows a large mixing of the highest lying occupied sigma* orbital of [Xe2] and the orbitals of the C60 cage. There is only a small gap between the a2u HOMO of Xe2@C60 and the eu LUMO and the a2u LUMO+1. The calculations show that there are several triplet states which are close in energy to each other and to the 1A1g state. The bonding analysis suggests that the interacting species in Xe2@C60 are the charged species Xe2q+ and C60q-, where 1相似文献   

The Structural and Electronic Properties of Tin-Based Heteroatom Clusters Studied by the Density Functional Theory

The focus of this study is on compound clusters and, due to the existence of many phases with different structural properties, tin-based materials have been chosen as the reference case. The clusters considered below are of two types: in the first case the clusters have the skeleton of the pure tin clusters and are doped with oxygen and aluminum atoms with composition Sn _x Y _y with Y = Al, O, x = 1, 10 and y = 1, 2. In the second case the clusters have a rutile lattice with a columnar or a spherical shape and a size up to 80 atoms and are doped with a number of aluminum atoms up to 20. The calculations are based on the Density Functional Theory (DFT) and the results describe the cluster structure, its binding energy and the density of states (DOS). The general indication of the calculations is that the additive coordinates outside, rather than inside, the pristine skeleton with the formation of hybrid bonds with properties similar to the ones of the host atoms. Also conspicuous effects of hybridization are observed in the electronic structure and, due to these effects, the binding energy may decrease with respect to the one in the undoped clusters.     

Diels–Alder Reaction on Free C68 Fullerene and Endohedral Sc3N@C68 Fullerene Violating the Isolated Pentagon Rule: Importance of Pentagon Adjacency

The reaction mechanism and regioselectivity of the Diels–Alder reactions of C₆₈ and Sc₃N@C₆₈, which violate the isolated pentagon rule, were studied with density functional theory calculations. For C₆₈, the [5,5] bond is the most favored thermodynamically, whereas the cycloaddition on the [5,6] bond has the lowest activation energy. Upon encapsulation of the metallic cluster, the exohedral reactivity of Sc₃N@C₆₈ is reduced remarkably owing to charge transfer from the cluster to the fullerene cage. The [5,5] bond becomes the most reactive site thermodynamically and kinetically. The bonds around the pentagon adjacency show the highest chemical reactivity, which demonstrates the importance of pentagon adjacency. Furthermore, the viability of Diels–Alder cycloadditions of several dienes and Sc₃N@C₆₈ was examined theoretically. o‐Quinodimethane is predicted to react with Sc₃N@C₆₈ easily, which implies the possibility of using Diels–Alder cycloaddition to functionalize Sc₃N@C₆₈.     

Bonding in Barium Boryloxides,Siloxides, Phenoxides and Silazides: A Comparison with the Lighter Alkaline Earths

Barium complexes ligated by bulky boryloxides [OBR₂]⁻ (where R=CH(SiMe₃)₂, 2,4,6-ⁱPr₃-C₆H₂ or 2,4,6-(CF₃)₃-C₆H₂), siloxide [OSi(SiMe₃)₃]⁻, and/or phenoxide [O-2,6-Ph₂-C₆H₃]⁻, have been prepared. A diversity of coordination patterns is observed in the solid state for both homoleptic and heteroleptic complexes, with coordination numbers ranging between 2 and 4. The identity of the bridging ligand in heteroleptic dimers [Ba(μ₂-X₁)(X₂)]₂ depends largely on the given pair of ligands X₁ and X₂. Experimentally, the propensity to fill the bridging position increases according to [OB{CH(SiMe₃)₂}₂)]⁻<[N(SiMe₃)₂]⁻<[OSi(SiMe₃)₃]⁻<[O(2,6-Ph₂-C₆H₃)]⁻<[OB(2,4,6-ⁱPr₃-C₆H₂)₂]⁻. This trend is the overall expression of 3 properties: steric constraints, electronic density and σ- and π-donating capability of the negatively charged atom, and ability to generate Ba ⋅ ⋅ ⋅ F, Ba ⋅ ⋅ ⋅ C(π) or Ba ⋅ ⋅ ⋅ H−C secondary interactions. The comparison of the structural motifs in the complexes [Ae{μ₂-N(SiMe₃)₂}(OB{CH(SiMe₃)₂}₂)]₂ (Ae = Mg, Ca, Sr and Ba) suggest that these observations may be extended to all alkaline earths. DFT calculations highlight the largely prevailing ionic character of ligand-Ae bonding in all compounds. The ionic character of the Ae-ligand bond encourages bridging coordination, whereas the number of bridging ligands is controlled by steric factors. DFT computations also indicate that in [Ba(μ₂-X₁)(X₂)]₂ heteroleptic dimers, ligand predilection for bridging vs. terminal positions is dictated by the ability to establish secondary interactions between the metals and the ligands.     

Mechanism of the Asymmetric Autocatalytic Soai Reaction Studied by Density Functional Theory

The mechanism of the Soai reaction has been thoroughly investigated at the M05‐2X/6‐31G(d) level of theory, by considering ten energetically distinct paths. The study indicates the fully enantioselective catalytic cycle of the homochiral dimers to be the dominant mechanism. Two other catalytic cycles are shown to both be important for correct understanding of the Soai reaction. These are the catalytic cycle of the heterochiral dimer and the non‐enantioselective catalytic cycle of the homochiral dimers. The former has been proved to be not really competitive with the principal cycle, as required for the Soai reaction to manifest chiral amplification, whereas the latter, which is only slightly competitive with the principal one, nicely explains the experimental enantioselectivity observed in the reaction of 2‐methylpyrimidine‐5‐carbaldehyde. The study has also evidenced the inadequacy of the B3LYP functional for mechanistic investigations of the Soai reaction.     

Crystallographic Evidence for Direct Metal–Metal Bonding in a Stable Open‐Shell La2@Ih‐C80 Derivative

Endohedral metallofullerenes (EMFs) have novel structures and properties that are closely associated with the internal metallic species. Benzyl radical additions have been previously shown to form closed‐shell adducts by attaching an odd number of addends to open‐shell EMFs (such as Sc₃C₂@I_h‐C₈₀) whereas an even number of groups are added to closed‐shell EMFs (for example Sc₃N@I_h‐C₈₀). Herein we report that benzyl radical addition to the closed‐shell La₂@I_h‐C₈₀ forms a stable, open‐shell monoadduct instead of the anticipated closed‐shell bisadduct. Single‐crystal X‐ray diffraction results show the formation of a stable radical species. In this species, the La?La distance is comparable to the theoretical value of a La?La covalent bond and is shorter than reported values for other La₂@I_h‐C₈₀ derivatives, providing unambiguous evidence for the formation of direct La?La bond.