Synthesis,Isolation, and Trifluoromethylation of Two Isomers of C84‐Based Monometallic Cyanide Clusterfullerenes: Interplay between the Endohedral Cluster and the Exohedral Addends

As an emerging member of endohedral fullerenes, metal cyanide clusterfullerenes (CYCF) are unique in terms of the encapsulation of a monometallic cluster. To date the reported carbon cages of CYCFs are limited to C₈₂ and C₇₆, and little is known about the chemical reactivity of CYCFs. Herein, two isomers of the first C₈₄‐based CYCFs, YCN@C₈₄, were isolated as trifluoromethyl derivatives, including YCN@C₈₄(23)(CF₃)₁₈ and three isomers of YCN@C₈₄(13)(CF₃)₁₆, which are based on a unique chiral C ₂‐C₈₄(13) cage. As a common feature of the CF₃ addition patterns, the YCN@C₈₄(CF₃)_16/18 compounds are stabilized by the formation of isolated C=C bonds and benzenoid rings on the carbon cages. The interplay between the endohedral YCN cluster and the exhohedral CF₃ addends was unveiled according to single‐crystal X‐ray diffraction studies, thus offering new insight into the chemical reactivity of CYCFs.     

Synthesis,Isolation, and Trifluoromethylation of Two Isomers of C84‐Based Monometallic Cyanide Clusterfullerenes: Interplay between the Endohedral Cluster and the Exohedral Addends

As an emerging member of endohedral fullerenes, metal cyanide clusterfullerenes (CYCF) are unique in terms of the encapsulation of a monometallic cluster. To date the reported carbon cages of CYCFs are limited to C₈₂ and C₇₆, and little is known about the chemical reactivity of CYCFs. Herein, two isomers of the first C₈₄‐based CYCFs, YCN@C₈₄, were isolated as trifluoromethyl derivatives, including YCN@C₈₄(23)(CF₃)₁₈ and three isomers of YCN@C₈₄(13)(CF₃)₁₆, which are based on a unique chiral C ₂‐C₈₄(13) cage. As a common feature of the CF₃ addition patterns, the YCN@C₈₄(CF₃)_16/18 compounds are stabilized by the formation of isolated C=C bonds and benzenoid rings on the carbon cages. The interplay between the endohedral YCN cluster and the exhohedral CF₃ addends was unveiled according to single‐crystal X‐ray diffraction studies, thus offering new insight into the chemical reactivity of CYCFs.     

Monometallic cyanide cluster fullerene YCN@C78: A theoretical prediction

Stimulated by the recent experimental success in production and characterization of YCN@C_s(6)‐C₈₂, the possibility of encapsulating YCN cluster in the C₇₈ fullerene has been performed using density functional theory. Six isomers of YCN@C₇₈ are considered based on six lowest energy C₇₈^2? isomers. The results reveal that YCN@D_3h(24109)‐C₇₈ and YCN@C_2v(24107)‐C₇₈, both of which satisfy the isolated‐pentagon rule, present excellent thermodynamic stability with very small energy differences. Moreover, the large HOMO‐LUMO gaps (1.55 and 1.47 eV for YCN@D_3h(24109)‐C₇₈ and YCN@C_2v(24107)‐C₇₈, respectively) indicate their high kinetic stabilities. Significantly, in both the structures, the encapsulated YCN cluster is triangular, similar to the cases of YCN@C_s(6)‐C₈₂ and TbCN@C₂(5)‐C₈₂. In addition, electronic absorption spectra, infrared spectra, and ¹³C nuclear magnetic resonance spectra of two stable structures have also been explored to further disclose the molecular structures and properties. © 2015 Wiley Periodicals, Inc.     

Search for the global minimum structures of AlB3H2n (n = 0 − 6) clusters

The global minimum structures of AlB₃H_2n (n = 0–6) clusters are determined using the stochastic search method at the B3LYP/6–31G level of theory. These initially specified geometries are recalculated using B3LYP and CCSD(T) methods using the 6–311++G^** basis set. The structural and electronic properties of the two lowest‐lying isomers are presented. The structural parameters obtained for aluminum borohydride are compared with the experimental and theoretical results. The H₂ fragmentation energies of the most stable isomers are investigated. Chemical bonding analyses for the global minimum of AlB₃H_2n (n = 0–6) clusters are performed using the adaptive natural density partitioning method. © 2014 Wiley Periodicals, Inc.     

The Evolution of Bonding and Thermodynamic Properties of Boron‐Doped Small Carbon Clusters: An Ab Initio Study

Theoretical studies on BC_n (n=1–6) clusters are carried out using density functional theory, Møller–Plesset second‐order perturbation theory (MP2), coupled‐cluster calculations including up to triple excitations (CCSD(T)), and higher‐level approaches. All possible isomers depending on the positions of the boron atom are generated and the lowest‐energy isomers are determined for doublet and quartet electronic states. The three potential evolution paths of the clusters are determined as a function of their size. The energetic and electronic consequences for the increased size of structures differ significantly, which leads to representatives of the ground electronic state from different structural groups. The ab initio calculated thermal functions allow enhancements to the available atomization energies and improve the agreement between the calculated and experimental heat content.     

Structures and electronic properties of the small rubidium‐doped silicon RbSin (n = 1–12) clusters

The geometries, relative stabilities, and electronic properties of small rubidium‐doped silicon clusters RbSi_n (n = 1–12) have been systematically investigated using the density functional theory at the B3LYP/GENECP level. The optimized structures show that lowest‐energy isomers of RbSi_n are similar with the ground state isomers of pure Si_n clusters and prefer the three‐dimensional for n = 3–12. The relative stabilities of RbSi_n clusters have been analyzed on the averaged binding energy, fragmentation energy, second‐order energy difference, and highest occupied molecular orbital‐lowest unoccupied molecular orbital energy gap. The calculated results indicate that the doping of Rb atom enhances the chemical activity of Si_n frame and the magic number is RbSi₂. The Mulliken population analysis reveals that the charges in the corresponding RbSi_n clusters transfer from the Rb atom to Si atoms. The partial density of states and chemical hardness are also discussed. © 2014 Wiley Periodicals, Inc.     

A high‐accuracy theoretical study of the CHnP Systems n = 1–3

We have performed high‐level electronic structure computations on the most important species of the CH_nP systems n = 1–3 to characterize them and provide reliable information about the equilibrium and vibrationally averaged molecular structures, rotational constants, vibrational frequencies (harmonic and anharmonic), formation enthalpies, and vertical excitation energies. Those chemical systems are intermediates for several important reactions and also prototypical phosphorus‐carbon compounds; however, they are often elusive to experimental detection. The present results significantly complement their knowledge and can be used as an assessment of the experimental information when available. The explicitly correlated coupled‐cluster RCCSD(T)‐F12 method has been used for geometry optimizations and vibrational frequency calculations. Vibrational configuration interaction theory has been used to account for anharmonicity effects. Basis‐set limit extrapolations have been carried out to determine accurate thermochemical quantities. Electronic excited states have been calculated with coupled‐cluster approaches and also by means of the multireference configuration interaction method. © 2013 Wiley Periodicals, Inc.     

Structural interconnections and the role of heptagonal rings in endohedral trimetallic nitride template fullerenes

Recent experiments indicate that fullerene isomers outside the classical definition can also encapsulate metallic atoms or clusters to form endohedral metallofullerenes. Our systematic study using DFT calculations, suggests that many heptagon‐including nonclassical trimetallic nitride template fullerenes are similar in stability to their classical counterparts, and that conversion between low‐energy nonclassical and classical parent cages via Endo–Kroto insertion/extrusion of C₂ units and Stone–Wales isomerization may facilitate the formation of endohedral trimetallic nitride fullerenes. Close structural connections are found between favored isomers of trimetallic nitride template fullerenes from C₇₈ to C₈₂. It appears that the lower symmetry and local deformations associated with introduction of a heptagonal ring favor encapsulation of intrinsically less symmetrical mixed metal nitride clusters. © 2016 Wiley Periodicals, Inc.     

A theoretical study on the stability difference of the borane and carborane C2Bn−2Hn (5 ≤ n ≤ 7) clusters

In order to study the electronic structure and structural stability of borane and carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters, especially the stability difference between the borane and carborane C₂B₃H₅. The frontier orbital energy levels of the borane and carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters are calculated at CCSD(T)/aug‐cc‐pVXZ//B3LYP/def2‐TZVPP level. The results are further analyzed by qualitative frontier orbital method based on the cap–ring interaction. The results reveal that: (1) the larger E_gap(HOMO‐LUMO energy gap) of carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters than borane (5 ≤ n ≤ 7) clusters originates from the more effective cap–ring orbital overlap of carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters than that of borane (5 ≤ n ≤ 7) clusters; (2) the smallest E_gap of the borane results from the highest energy level of the ring symmetry‐adapted linear combination orbital of cluster; and (3) the largest E_gap of the carborane C₂B₃H₅ is induced by the most effective cap–ring orbital interaction of C₂B₃H₅ cluster. © 2014 Wiley Periodicals, Inc.     

Theoretical study of silicon–sulfur clusters (SiS2)n− (n=1–6)

The possible geometrical structures and relative stability of silicon–sulfur clusters (SiS₂) (n=1–6) are explored by means of density functional theory (DFT) quantum chemical calculations. We also compare DFT with second‐order Møller–Plesset (MP2) and Hartree–Fock (HF) methods. The effects of polarization functions, diffuse functions, and electron correlation are included in MP2 and B3LYP quantum chemical calculations, and B3LYP is effective in larger cluster structure optimization, so we can conclude that the DFT approach is useful in establishing trends. The electronic structures and vibrational spectra of the most stable geometrical structures of (SiS₂)_n⁻ are analyzed by B3LYP. As a result, the regularity of the (SiS₂)_n⁻ cluster growing is obtained, and the calculation may predict the formation mechanism of the (SiS₂)_n⁻ cluster. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 280–290, 2001     

Electronic and magnetic properties of small RhnCa (n = 1–9) clusters: A DFT study

The equilibrium geometries, relative stabilities, electronic and magnetic properties of small Rh_nCa (n = 1–9) clusters have been investigated by DFT calculations. The obtained results show that the three‐dimensional geometries are adopted for the lowest‐energy Rh_nCa clusters, and the doped Ca atom prefers locating on the surface of the cluster. Based on the analysis of the second‐order difference of energies, fragmentation energies and the HOMO‐LUMO energy gaps, we identify that the Rh₄Ca, Rh₆Ca, and Rh₈Ca clusters are relatively more stable than their neighboring clusters, and the doping of Ca enhances the chemical reactivity of the pure Rh_n clusters, suggesting that the Rh_nCa clusters can be used as nanocatalysts in many catalytic reactions. The magnetic moment for these clusters is mostly localized on the Rh atoms, and the doping Ca atom has no effect on the total magnetic moment of Rh_nCa clusters. The partial density of states, VIP, VEA, and η of these clusters in their ground‐state structures were also calculated and discussed. © 2015 Wiley Periodicals, Inc.     

Far‐Infrared Spectra of Yttrium‐Doped Gold Clusters AunY (n=1–9)

The geometric, spectroscopic, and electronic properties of neutral yttrium‐doped gold clusters Au_nY (n=1–9) are studied by far‐infrared multiple photon dissociation (FIR‐MPD) spectroscopy and quantum chemical calculations. Comparison of the observed and calculated vibrational spectra allows the structures of the isomers present in the molecular beam to be determined. Most of the isomers for which the IR spectra agree best with experiment are calculated to be the energetically most stable ones. Attachment of xenon to the Au_nY cluster can cause changes in the IR spectra, which involve band shifts and band splittings. In some cases symmetry changes, as a result of the attachment of xenon atoms, were also observed. All the Au_nY clusters considered prefer a low spin state. In contrast to pure gold clusters, which exhibit exclusively planar lowest‐energy structures for small sizes, several of the studied species are three‐dimensional. This is particularly the case for Au₄Y and Au₉Y, while for some other sizes (n=5, 8) the 3D structures have an energy similar to that of their 2D counterparts. Several of the lowest‐energy structures are quasi‐2D, that is, slightly distorted from planar shapes. For all the studied species the Y atom prefers high coordination, which is different from other metal dopants in gold clusters.     

Theoretical investigation of LaC3n+ (n = 0, 1, 2) clusters by density functional theory

LaC₃ⁿ⁺ (n=0, 1, 2) clusters have been studied using B3LYP (Becke 3-parameter–Lee-Yang-Parr) density functional method. The basis set is Dunning/Huzinaga valence double zeta for carbon and [2s2p2d] for lanthanum, denoted LANL1DZ. Four isomers are presented for each cluster; two of them are edge binding isomers with C_2v symmetry, the other two are linear chains with C_∞v symmetry. Meanwhile, two spin states for each isomer, that is, singlet and triplet for LaC₃⁺, doublet and quartet for LaC₃ and LaC₃²⁺, respectively, are also considered. Geometries, vibrational frequencies, infrared intensities, and other quantities are reported and discussed. The results indicate that at some spin states; the C_2v symmetry isomers are the dominant structures, while for the other spin states, linear isomers are energetically favored. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 301–307, 1998     

Nonclassical fullerenes with a heptagon violating the pentagon adjacency penalty rule

Nonclassical fullerenes with heptagon(s) and their derivatives have attracted increasing attention, and the studies on them are performing to enrich the chemistry of carbon. Density functional theory calculations are performed on nonclassical fullerenes C_n (n = 46, 48, 50, and 52) to give insight into their structures and stability. The calculated results demonstrate that the classical isomers generally satisfy the pentagon adjacency penalty rule. However, the nonclassical isomers with a heptagon are more energetically favorable than the classical ones with the same number of pentagon–pentagon bonds (B₅₅ bonds), and many of them are even more stable than some classical isomers with fewer B₅₅ bonds. The nonclassical isomers with the lowest energy are higher in energy than the classical ones with the lowest energy, because they have more B₅₅ bonds. Generally, the HOMO–LUMO gaps of the former are larger than those of the latter. The sphericity and asphericity are unable to rationalize the unique stability of the nonclassical fullerenes with a heptagon. The pyramidization angles of the vertices shared by two pentagons and one heptagon are smaller than those of the vertices shared by two pentagons and one hexagon. It is concluded that the strain in the fused pentagons can be released by the adjacent heptagons partly, and consequently, it is a common phenomenon for nonclassical fullerenes to violate the pentagon adjacent penalty rule. These findings are heuristic and conducive to search energetically favorable isomers of C_n, especially as n is 62, 64, 66, and 68, respectively. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Computational studies of stable hexanuclear CulAgmAun (l + m + n = 6; l,m, n > 0) clusters

A DFT study was carried out on the ground state structures of ternary Cu_lAg_mAu_n (l + m + n = 6) clusters, with the aim of investigating changes of thermal and kinetic stabilities as an effect of composition, as well as the composition dependence of the electrostatic potential, of stable planar structures. DFT optimizations were performed using the PBE functional and the SDD basis set. All the optimized structures adopt planar geometries with bent triangular structures. Calculated binding energy values are in the range 1.5–1.9 eV/atom, which shows their thermal stability. The predicted HOMO‐LUMO energy gap values are in the semiconductor region, providing a qualitative indication of a moderate kinetic stability. NBO analyses indicate the existence of two mechanisms promoting planar structural stability, one due to bonding‐antibonding orbital interaction, and the other one due to the well‐known spd hybridization. Wiberg indices were obtained showing interatomic bonding. Electrostatic potential calculations show the existence of nucleophilic attack regions preferentially around silver and copper atoms located at the vertices while electrophilic attack regions are found in the vicinity of gold atoms over the cluster plane. Apparently, charge transfer occurs toward gold from silver and copper atoms when the concentration is favorable in the proximity of gold atoms. In particular, if the small ternary clusters discussed here contain only one gold atom, then a high electron density is observed at the site of this gold atom. © 2016 Wiley Periodicals, Inc.     

Electronic,bonding, and optical properties of 1d [CuCN]n (n = 1–10) chains, 2d [CuCN]n (n = 2–10) nanorings,and 3d [Cun(CN)n]m (n = 4, m = 2, 3; n = 10, m = 2) tubes studied by DFT/TD‐DFT methods

The electronic, bonding, and photophysical properties of one‐dimensional [CuCN]_n (n = 1–10) chains, 2‐D [CuCN]_n (n = 2–10) nanorings, and 3‐D [Cu_n(CN)_n]_m (n = 4, m = 2, 3; n = 10, m = 2) tubes are investigated by means of a multitude of computational methodologies using density functional theory (DFT) and time‐dependent‐density‐functional theory (TD‐DFT) methods. The calculations revealed that the 2‐D [CuCN]_n (n = 2–10) nanorings are more stable than the respective 1‐D [CuCN]_n (n = 2–10) linear chains. The 2‐D [CuCN]_n (n = 2–10) nanorings are predicted to form 3‐D [Cu_n(CN)_n]_m (n = 4, m = 2, 3; n = 10, m = 2) tubes supported by weak stacking interactions, which are clearly visualized as broad regions in real space by the 3D plots of the reduced density gradient. The bonding mechanism in the 1‐D [CuCN]_n (n = 1–10) chains, 2‐D [CuCN]_n (n = 2–10) nanorings, and 3‐D [Cu_n(CN)_n]_m (n = 4, m = 2, 3; n = 10, m = 2) tubes are easily recognized by a multitude of electronic structure calculation approaches. Particular emphasis was given on the photophysical properties (absorption and emission spectra) of the [CuCN]_n chains, nanorings, and tubes which were simulated by TD‐DFT calculations. The absorption and emission bands in the simulated TD‐DFT absorption and emission spectra have thoroughly been analyzed and assignments of the contributing principal electronic transitions associated to individual excitations have been made. © 2015 Wiley Periodicals, Inc.     

1‐Aryl‐4‐Silylmethyl[60]fullerenes: Synthesis,Properties, and Photovoltaic Performance

The efficient nucleophilic addition of aryl Grignard reagents (aryl=4‐MeOC₆H₄, 4‐Me₂NC₆H₄, Ph, 4‐CF₃C₆H₄, and thienyl) to C₆₀ in the presence of DMSO produced 1,2‐arylhydro[60]fullerenes after acid treatment. The reactions of the anions of these arylhydro[60]fullerenes with either dimethylphenylsilylmethyl iodide or dimethyl(2‐isopropoxyphenyl)silylmethyl iodide yielded the target compounds, 1‐aryl‐4‐silylmethyl[60]fullerenes. The properties and structures of these 1‐aryl‐4‐silylmethyl[60]fullerenes (aryl=4‐MeOC₆H₄, thienyl) were examined by electrochemical studies, X‐ray crystallography, flash‐photolysis time‐resolved microwave‐conductivity (FP‐TRMC) measurements, and electron‐mobility measurements by using a space‐charge‐limited current (SCLC) model. Organic photovoltaic devices with a polymer‐based bulk heterojunction structure and small‐molecule‐based p–n and p–i–n heterojunction configurations were fabricated by using 1‐aryl‐4‐silylmethyl[60]fullerenes as an electron acceptor. The most efficient device exhibited a power‐conversion efficiency of 3.4 % (short‐circuit current density: 8.1 mA/ cm², open‐circuit voltage: 0.69 V, fill factor: 0.59).     

Theoretical study on monometallic cyanide cluster fullerenes YCN@C72

The geometries, stabilities, and electronic properties of new endohedral fullerene YCN@C₇₂ have been investigated by the B3LYP and PBE1PBE density functional (DFT) methods. The C_2v(11188)‐C₇₂ cage, which violates the isolated pentagon rule (IPR) with a pair of fused pentagons, is predicted to be the lowest energy isomer for both empty and YCN@C₇₂. The relatively large HOMO‐LUMO gap (B3LYP: 1.48 eV, PBE1PBE: 1.68 eV) for YCN@C_2v(11188)‐C₇₂ reveals this structure kinetic stability. Significantly, the encased YCN cluster adopts a triangular structure inside the C_2v(11188)‐C₇₂ cage, similar to the results reported on YCN@C_s(6)‐C₈₂ and TbCN@C₂(5)‐C₈₂. Furthermore, the vertical ionization potential and electron affinity, UV‐vis‐NIR and IR spectra of YCN@C_2v(11188)‐C₇₂ have been predicted to facilitate future experimental characterization. © 2015 Wiley Periodicals, Inc.     

Synthesis and Isolation of the Titanium–Scandium Endohedral Fullerenes—Sc2TiC@Ih‐C80, Sc2TiC@D5h‐C80 and Sc2TiC2@Ih‐C80: Metal Size Tuning of the TiIV/TiIII Redox Potentials

The formation of endohedral metallofullerenes (EMFs) in an electric arc is reported for the mixed‐metal Sc–Ti system utilizing methane as a reactive gas. Comparison of these results with those from the Sc/CH₄ and Ti/CH₄ systems as well as syntheses without methane revealed a strong mutual influence of all key components on the product distribution. Whereas a methane atmosphere alone suppresses the formation of empty cage fullerenes, the Ti/CH₄ system forms mainly empty cage fullerenes. In contrast, the main fullerene products in the Sc/CH₄ system are Sc₄C₂@C₈₀ (the most abundant EMF from this synthesis), Sc₃C₂@C₈₀, isomers of Sc₂C₂@C₈₂, and the family Sc₂C_{2 n} (2 n=74, 76, 82, 86, 90, etc.), as well as Sc₃CH@C₈₀. The Sc–Ti/CH₄ system produces the mixed‐metal Sc₂TiC@C_{2 n} (2 n=68, 78, 80) and Sc₂TiC₂@C_{2 n} (2 n=80) clusterfullerene families. The molecular structures of the new, transition‐metal‐containing endohedral fullerenes, Sc₂TiC@I_h‐C₈₀, Sc₂TiC@D_5h‐C₈₀, and Sc₂TiC₂@I_h‐C₈₀, were characterized by NMR spectroscopy. The structure of Sc₂TiC@I_h‐C₈₀ was also determined by single‐crystal X‐ray diffraction, which demonstrated the presence of a short Ti=C double bond. Both Sc₂TiC‐ and Sc₂TiC₂‐containing clusterfullerenes have Ti‐localized LUMOs. Encapsulation of the redox‐active Ti ion inside the fullerene cage enables analysis of the cluster–cage strain in the endohedral fullerenes through electrochemical measurements.