Computationally Designed Families of Flat,Tubular, and Cage Molecules Assembled with “Starbenzene” Building Blocks through Hydrogen‐Bridge Bonds

Using density functional calculations, we demonstrate that the planarity of the nonclassical planar tetracoordinate carbon (ptC) arrangement can be utilized to construct new families of flat, tubular, and cage molecules which are geometrically akin to graphenes, carbon nanotubes, and fullerenes but have fundamentally different chemical bonds. These molecules are assembled with a single type of hexagonal blocks called starbenzene (D_6h C₆Be₆H₆) through hydrogen‐bridge bonds that have an average bonding energy of 25.4–33.1 kcal mol^?1. Starbenzene is an aromatic molecule with six π electrons, but its carbon atoms prefer ptC arrangements rather than the planar trigonal sp² arrangements like those in benzene. Various stability assessments indicate their excellent stabilities for experimental realization. For example, one starbenzene unit in an infinite two‐dimensional molecular sheet lies on average 154.1 kcal mol^?1 below three isolated linear C₂Be₂H₂ (global minimum) monomers. This value is close to the energy lowering of 157.4 kcal mol^?1 of benzene relative to three acetylene molecules. The ptC bonding in starbenzene can be extended to give new series of starlike monocyclic aromatic molecules (D_4h C₄Be₄H₄^2?, D_5h C₅Be₅H₅^?, D_6h C₆Be₆H₆, D_7h C₇Be₇H₇⁺, D_8h C₈Be₈H₈^2?, and D_9h C₉Be₉H₉^?), known as starenes. The starene isomers with classical trigonal carbon sp² bonding are all less stable than the corresponding starlike starenes. Similarly, lithiated C₅Be₅H₅ can be assembled into a C₆₀‐like molecule. The chemical bonding involved in the title molecules includes aromaticity, ptC arrangements, hydrogen‐bridge bonds, ionic bonds, and covalent bonds, which, along with their unique geometric features, may result in new applications.     

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.     

Starlike aluminum-carbon aromatic species

Is it possible to achieve molecules with starlike structures by replacing the H atoms in (CH)_n^q aromatic hydrocarbons with aluminum atoms in bridging positions? Although D_4h C₄Al₄^2? and D₂ C₆Al₆ are not good prospects for experimental realization, a very extensive computational survey of fifty C₅Al₅^? isomers identified the starlike D_5h global minimum with five planar tetracoordinate carbon atoms to be a promising candidate for detection by photoelectron detachment spectroscopy. BOMD (Born–Oppenheimer molecular dynamics) simulations and high‐level theoretical computations verified this conclusion. The combination of favorable electronic and geometric structural features (including aromaticity and optimum C–Al–C bridge bonding) stabilizes the C₅Al₅^? star preferentially.     

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.     

Heavy-Atom Tunneling Through Crossing Potential Energy Surfaces: Cyclization of a Triplet 2-Formylarylnitrene to a Singlet 2,1-Benzisoxazole

Not long ago, the occurrence of quantum mechanical tunneling (QMT) chemistry involving atoms heavier than hydrogen was considered unreasonable. Contributing to the shift of this paradigm, we present here the discovery of a new and distinct heavy-atom QMT reaction. Triplet syn-2-formyl-3-fluorophenylnitrene, generated in argon matrices by UV-irradiation of an azide precursor, was found to spontaneously cyclize to singlet 4-fluoro-2,1-benzisoxazole. Monitoring the transformation by IR spectroscopy, temperature-independent rate constants (k≈1.4×10⁻³ s⁻¹; half-life of ≈8 min) were measured from 10 to 20 K. Computational estimated rate constants are in fair agreement with experimental values, providing evidence for a mechanism involving heavy-atom QMT through crossing triplet to singlet potential energy surfaces. Moreover, the heavy-atom QMT takes place with considerable displacement of the oxygen atom, which establishes a new limit for the heavier atom involved in a QMT reaction in cryogenic matrices.     

The Exohedral Diels–Alder Reactivity of the Titanium Carbide Endohedral Metallofullerene Ti2C2@D3h‐C78: Comparison with D3h‐C78 and M3N@D3h‐C78 (M=Sc and Y) Reactivity

The chemical functionalization of endohedral (metallo)fullerenes has become a main focus of research in the last few years. It has been found that the reactivity of endohedral (metallo)fullerenes may be quite different from that of the empty fullerenes. Encapsulated species have an enormous influence on the thermodynamics, kinetics, and regiochemistry of the exohedral addition reactions undergone by these species. A detailed understanding of the changes in chemical reactivity due to incarceration of atoms or clusters of atoms is essential to assist the synthesis of new functionalized endohedral fullerenes with specific properties. Herein, we report the study of the Diels–Alder cycloaddition between 1,3‐butadiene and all nonequivalent bonds of the Ti₂C₂@D_3h‐C₇₈ metallic carbide endohedral metallofullerene (EMF) at the BP86/TZP//BP86/DZP level of theory. The results obtained are compared with those found by some of us at the same level of theory for the D_3h‐C₇₈ free cage and the M₃N@D_3h‐C₇₈ (M=Sc and Y) metallic nitride EMFs. It is found that the free cage is more reactive than the Ti₂C₂@D_3h‐C₇₈ EMF and this, in turn, has a higher reactivity than M₃N@D_3h‐C₇₈. The results indicate that, for Ti₂C₂@D_3h‐C₇₈, the corannulene‐type [5, 6] bonds c and f , and the type B [6, 6] bond 3 are those thermodynamically and kinetically preferred. In contrast, the D_3h‐C₇₈ free cage has a preference for addition to the [6, 6] 1 and 6 bonds and the [5, 6] b bond, whereas M₃N@D_3h‐C₇₈ favors additions to the [6, 6] 6 (M=Sc) and [5, 6] d (M=Y) bonds. The reasons for the regioselectivity found in Ti₂C₂@D_3h‐C₇₈ are discussed.     

Heavy‐Atom Tunneling Through Crossing Potential Energy Surfaces: Cyclization of a Triplet 2‐Formylarylnitrene to a Singlet 2,1‐Benzisoxazole

Not long ago, the occurrence of quantum mechanical tunneling (QMT) chemistry involving atoms heavier than hydrogen was considered unreasonable. Contributing to the shift of this paradigm, we present here the discovery of a new and distinct heavy‐atom QMT reaction. Triplet syn‐2‐formyl‐3‐fluorophenylnitrene, generated in argon matrices by UV‐irradiation of an azide precursor, was found to spontaneously cyclize to singlet 4‐fluoro‐2,1‐benzisoxazole. Monitoring the transformation by IR spectroscopy, temperature‐independent rate constants (k≈1.4×10^?3 s^?1; half‐life of ≈8 min) were measured from 10 to 20 K. Computational estimated rate constants are in fair agreement with experimental values, providing evidence for a mechanism involving heavy‐atom QMT through crossing triplet to singlet potential energy surfaces. Moreover, the heavy‐atom QMT takes place with considerable displacement of the oxygen atom, which establishes a new limit for the heavier atom involved in a QMT reaction in cryogenic matrices.     

Tunneling Assists the 1,2‐Hydrogen Shift in N‐Heterocyclic Carbenes

At room temperature, 1,2‐hydrogen‐transfer reactions of N‐heterocyclic carbenes, like the imidazol‐2‐ylidene to give imidazole is shown to occurr almost entirely (>90 %) by quantum mechanical tunneling (QMT). At 60 K in an Ar matrix, for the 2, 3‐dihydrothiazol‐2‐ylidene→thiazole transformation, QMT is shown to increase the rate about 10⁵ times. Calculations including small‐curvature tunneling show that the barrier for intermolecular 1,2‐hydrogen‐transfer reaction is small, and QMT leads to a reduced rate of the forward reaction because of nonclassical reflections even at room temperature. A small barrier also leads to smaller kinetic isotope effects because of efficient QMT by both H and D. QMT does not always lead to faster reactions or larger KIE values, particularly when the barrier is small.     

Full Exploration of the Diels–Alder Cycloaddition on Metallofullerenes M3N@C80 (M=Sc,Lu, Gd): The D5h versus Ih Isomer and the Influence of the Metal Cluster

In this work a detailed investigation of the exohedral reactivity of the most important and abundant endohedral metallofullerene (EMF) is provided, that is, Sc₃N@I_h‐C₈₀ and its D_5h counterpart Sc₃N@D_5h‐C₈₀, and the (bio)chemically relevant lutetium‐ and gadolinium‐based M₃N@I_h/D_5h‐C₈₀ EMFs (M=Sc, Lu, Gd). In particular, we analyze the thermodynamics and kinetics of the Diels–Alder cycloaddition of s‐cis‐1,3‐butadiene on all the different bonds of the I_h‐C₈₀ and D_5h‐C₈₀ cages and their endohedral derivatives. First, we discuss the thermodynamic and kinetic aspects of the cycloaddition reaction on the hollow fullerenes and the two isomers of Sc₃N@C₈₀. Afterwards, the effect of the nature of the metal nitride is analyzed in detail. In general, our BP86/TZP//BP86/DZP calculations indicate that [5,6] bonds are more reactive than [6,6] bonds for the two isomers. The [5,6] bond D _5h ‐b , which is the most similar to the unique [5,6] bond type in the icosahedral cage, I _h ‐a , is the most reactive bond in M₃N@D_5h‐C₈₀ regardless of M. Sc₃N@C₈₀ and Lu₃N@C₈₀ give similar results; the regioselectivity is, however, significantly reduced for the larger and more electropositive M=Gd, as previously found in similar metallofullerenes. Calculations also show that the D_5h isomer is more reactive from the kinetic point of view than the I_h one in all cases which is in good agreement with experiments.     

Activation of Molecular Hydrogen by a Singlet Carbene through Quantum Mechanical Tunneling

Carbenes are among the few metal‐free molecules that are able to activate molecular hydrogen. Whereas triplet carbenes have been shown to insert into H₂ through a two‐step mechanism that at low temperature is assisted by quantum mechanical tunneling (QMT), singlet carbenes insert in concerted reactions with considerable activation barriers, and are thus unreactive towards H₂ at cryogenic temperatures. Here we show that 1‐azulenylcarbene with a singlet ground state readily inserts into H₂, and that QMT governs the insertion into both H₂ and D₂. This is the first example that shows that QMT can also be important for singlet carbenes inserting into dihydrogen.     

A classification of metaborate crystals based on Raman spectroscopy

The Raman spectra of crystals built solely of metaborate triangles provide fingerprint identification of three distinct network types. Classified according to increasing cation field strength these are: rings with degenerate intra-annular bonds and D_3h symmetry, distorted rings with alternating intra-annular bonds and C_3h symmetry, as well as chains. The occurrence of each network type has been associated with a characteristic range of cationic field strength. This approach led to the discovery of a hitherto unknown C_3h ring strontium metaborate crystal, with Sr²⁺ cations in 9- or 10-fold coordination to oxygen atoms. The Raman spectra of the mixed cation metaborates Ba₂Ca(B₂O₄)₃ and Ba₂Mg(B₂O₄)₃ confirm their C_3h-ring structure and clearly point to the fact that the synergetic effect of dissimilar cations to the metaborate network cannot be predicted by the additivity of their field strengths.     

A MNDO study of carbon clusters with specifically fitted parameters

Summary A MNDO method with new parameters for carbon clusters is presented. The parameters in the new sets are specifically tuned to fit the properties of small carbon clusters, C₂, C₃, C₅ and C₇–C₁₀, and buckminsterfullerene, C₆₀. The validity of these MNDO parameters is verified by experimental data. The calculated (with new parameters) IR spectra of C₆₀ and the heat of formation, geometry and IR spectra of C₇₀ agree satisfactorily with observed data. Heats of formation of other fullerenes, from C₂₀ to C₈₄, and C₆₀O are evaluated. The resulting heats of formation of the isomers of C₇₆ and C₈₄ are reliable and their relative stability is in excellent agreement with other reports. The predicted IR spectra of several fullerenes, C₂₄(C_6v ), C₂₈(T _d ), C₃₂(D₃), C₃₆(D_6h ), C₅₀(C_5h ) and C₈₀(D_5d ) are provided to aid assignments of experimental spectra.     

Mapping the Metal Positions inside Spherical C80 Cages: Crystallographic and Theoretical Studies of Ce2@D5h‐C80 and Ce2@Ih‐C80

The dynamic positions of the dimetallic cluster inside the mid‐sized spherical cages of C₈₀–C₈₂ have been seldom studied, despite the high abundance of M₂@C_2n (2n=80, 82) species among various endohedral metallofullerenes. Herein, using crystallographic methods, we first unambiguously map the metal positions for both Ce₂@D_5h‐C₈₀ and Ce₂@I_h‐C₈₀, showing how the symmetry or geometrical change in cage structure can influence the motional behavior of the cluster. Inside the D_5h cage, the primary cerium sites have been identified along a cage belt of the contiguous hexagons, which suggests the significant influence of such a cage motif on endohedral cluster motion. Further analysis revealed a distorted D_5h cage owing to the “punch‐out” effect of cerium atoms. The consequence is the presence of two localized electrostatic potential minima inside the cage of (D_5h‐C₈₀)^6?, thus reflecting the primary ionic cerium–cage interaction. In contrast, a different motional behavior of Ce₂ cluster was observed inside the I_h cage. With the major cerium sites, the molecule of Ce₂@I_h‐C₈₀ presented an approximate D_2h configuration. With the combined theoretical study, we propose that the additional unidentified influence of Ni^II(OEP) (OEP=octaethylporphyrin) might be also relevant for the location of cerium sites inside the I_h cage.     

Is the stereomutation of methane possible?

Large basis set ab initio calculations at correlated levels, including MP2, single reference, as well as multireference configuration interaction, carried out on the methane potential energy surface, have located and characterized a transition structure for stereomutation (one imaginary frequency). This structure is best described as a pyramidal complex between singlet methylene and a side-on hydrogen molecule with C_s symmetry. At the single reference CI level, it lies 105 kcal/mol above the methane T_d-ground state but is stable relative to dissociation into CH₂(¹A₁) and H₂ by 13 kcal/mol at 0 K (with harmonic zero point energy (ZPE) corrections for all structures). Dissociation of the transition state into triplet methylene and hydrogen also is endothermic (by 4 kcal/mol), but single bond rupture to give CH and H^. is 3 kcal/mol exothermic. Thus, it does not appear likely that methane can undergo stereomutation classically beneath the dissociation limit. Confirming earlier conclusions, side-on insertion of ¹A₁ CH₂ into H₂ in a perpendicular geometry occurs without activation energy. Planar (D_4h) methane (130.5 kcal/mol) has four imaginary frequencies. Two of these are degenerate and lead to equivalent planar C_2v structures with one three-center, two-electron bond and two two-electron bonds and two imaginary frequencies. The remaining imaginary frequencies of the D_4h form lead to tetrahedral (T_d) and pyramidal (C_4v) methane. The latter has three negative eigenvalues in the force-constant matrix; one of these leads to the T_d global minimum and the other to the C_s (parallel) stereomutation transition structure. Multireference CI calculations with a large atomic natural orbitals basis set produce similar results, with the electronic energy of the C_s stereomutation transition state 0.7 ± 0.5 kcal/mol higher than that of CH + H^. dissociation products, and a ZPE-corrected energy which is 5 ± 1 kcal/mol higher. Also considered are photochemical pathways for stereomutation and the possible effects of nuclear spin, inversion tunneling, and the parity-violating weak nuclear interaction on the possibility of an experimental detection of stereomutation in methane. © 1995 by John Wiley & Sons, Inc.     

Linear Group 13 E≡E Triple Bonds in E2Li62+

The triply bonded heavier main-group compounds have a textbook trans-bent geometry, in contrast to a familiar linear form found for the lightest analogues. Strikingly, the unexpected linear group 13 E≡E triple bonds were herein found in the D_4h-symmetry E₂Li₆²⁺ clusters, and they possess a large barrier (>18.0 kcal/mol) towards the dissociation of Li⁺. The perfectly surrounded Li₄ motifs and two linear coordinated Li atoms strongly suppress the increasing nonbonded electron density of heavier E atoms, making two degenerate π bonds and one multi-center σ bond in linear heavier main-group triple bonds. The surrounding Li₆ motifs not only creates an effective electronic structure to form a linear E≡E triple bond, but the resulting electrostatic interactions account for the highly stable global E₂Li₆²⁺ clusters.     

Characterization of the Retro-[3+2] Reaction of Protonated 2,3′-Bisindolylmethanes by Electrospray Ionization–Tandem Mass Spectrometry

Twelve 2,3′-bisindolylmethanes with various substituents were investigated using electrospray ionization quadrupole time-of-flight tandem mass spectrometry in positive ion mode. A retro-[3+2] reaction was observed in the collision-induced dissociation spectra of protonated 2,3′-bisindolylmethanes for the first time. The mechanism of retro-[3+2] reaction was concerted or stepwise. For the concerted pathway, carbon–carbon bonds of a protonated compound simultaneously cracked and the m/z 208 ion ([C₁₅H₁₀D₂N]⁺) was observed with hydrogen–deuterium exchange labeling. The stepwise pathway goes through 1,3-hydrogen migration twice and the m/z 208 ion ([C₁₅H₁₀D₂N]⁺) and m/z 207 ion ([C₁₅H₁₁DN]⁺) were detected with deuterium labeling. In the deuterium-labeled tandem mass spectrum for one compound, only the peak at m/z 208 was present at high abundance, suggesting that the concerted pathway is more likely. In addition, the substituents have no obvious trends on the ratios of the product intensity to the base intensity, further supporting the concerted pathway.     

Theoretical observation of hexaatomic molecules containing pentacoordinate planar carbon

The smallest molecules up to date containing a D5h pentacoordinate planar carbon (PPC) atom, CBe5 and CBe54-, are presented by means of ab initio calculations. To gain a better understanding about which electronic factors contribute to their stabilization, natural bond orbital (NBO) analysis and the nucleus independent chemical shifts (NICS) were calculated. The data reported here suggest that D5h CBe5 is σ aromaticity in nature, while in D5h CBe54- π aromaticity is dominating. The classical octet rule is well satisfied in both molecules, and is one of the fundamental reasons to understand the stability of the pentagon structures. The Be5 ring serves as σ donor in D5h CBe5, and π-acceptor in D5h CBe54-. The D5h CBe54- possessing 18 valence electrons with a closed-shell electron configuration is the most plau-sible candidate for experimental detection.     

Blending Non‐Group‐3 Transition Metal and Rare‐Earth Metal into a C80 Fullerene Cage with D5h Symmetry

Rare‐earth metals have been mostly entrapped into fullerene cages to form endohedral clusterfullerenes, whereas non‐Group‐3 transition metals that can form clusterfullerenes are limited to titanium (Ti) and vanadium (V), and both are exclusively entrapped within an I_h‐C₈₀ cage. Non‐Group‐3 transition‐metal‐containing endohedral fullerenes based on a C₈₀ cage with D_5h symmetry, V_xSc_3?xN@D_5h‐C₈₀ (x=1, 2), have now been synthesized, which exhibit two variable cluster compositions. The molecular structure of VSc₂N@D_5h‐C₈₀ was unambiguously determined by X‐ray crystallography. According to a comparative study with the reported Ti‐ and V‐containing clusterfullerenes based on a I_h‐C₈₀ cage and the analogous D_5h‐C₈₀‐based metal nitride clusterfullerenes containing rare‐earth metals only, the decisive role of the non‐Group‐3 transition metal on the formation of the corresponding D_5h‐C₈₀‐based clusterfullerenes is unraveled.     

Control of Chromophore Symmetry by Positional Isomerism of Peripheral Substituents

The first example of the control of porphyrinoid chromophore symmetry based on the positional isomerism of peripheral substituents has been achieved by preparing tetraazaporphyrins (TAPs) with C_4h, D_2h, C_2v, and C_s symmetry due to the relative arrangement of peripheral tert‐butylamino and cyano groups as push and pull substituents, respectively. The four structural isomers were successfully isolated and characterized by ¹H NMR spectroscopy and X‐ray crystallography. The band morphology in the Q‐band region varies depending on the molecular symmetry due to the significant perturbation introduced into the chromophore by the push and pull substituents. The C_4h and C_2v isomers exhibit a single Q band, whereas the Q bands of the D_2h and C_s isomers show a marked splitting. The magnetic circular dichroism spectra indicate that the push–pull TAPs retain the properties of the 16‐membered 18π‐electron perimeter generally observed for porphyrinoids. Theoretical calculations have demonstrated that the perturbation introduced by the substituents lowers the D_4h symmetry of the parent TAP π‐conjugated system, and this results in significant spectral changes. A novel approach to the fine‐tuning of the spectral properties of porphyrinoids based on changes in the chromophore symmetry is described.     

A quantum chemical study of [1.1.1.1]pagodane and its related compounds

The PBE0/6-31G** quantum chemical method is used to determine the symmetry and equilibrium structural parameters of the molecules of [1.1.1.1]pagodane (C₂₀H₂₀, D _2h ), two dienes (C₂₀H₂₀, D _2h ), two diradicals (C₂₀H₂₀, C _2ν ), and two dications (C₂₀H₂₀²⁺, D _2h and C _2ν ). The energy of a highly symmetric dication with a rectangular cycle is lower by 36 kcal/mole than that of a low symmetric dication with a trapezoidal cycle. The polarization interaction with liquid methylene chloride causes its decrease by 147 kcal/mole.