The C60 Formation Puzzle “Solved”: QM/MD Simulations Reveal the Shrinking Hot Giant Road of the Dynamic Fullerene Self-Assembly Mechanism

Formation of the icosahedral C60 fullerene via migration of single sp atoms and annihilation of sp-atom pairs

The disappearance of sp2 structural defects during abundant fullerene isomer formation is considered within the framework of the atomistic mechanism with participation of carbon atoms with sp hybridization. The study is carried out using the example of the icosahedral C60-Ih fullerene formation from the appropriate C58-C2v fullerene with a 7-ring. In this case the studied atomistic mechanism includes the following stages: (1) insertion of single carbon atoms into the fullerene from carbon vapor as an sp-atom instead of or above a bond, (2) directional migration of the sp-atom positions towards the 7-ring with decrease of energy, and (3) meeting of two sp atoms near the 7-ring with annihilation of the sp-atom pair and formation of the sp2 structure of the C60-Ih fullerene. The probabilities of all possible sp-atom positions on the appropriate C58-C2v fullerene shell are estimated as a function of temperature using the total energies of these positions obtained by spin-polarized density functional theory calculations using the PBE functional. Based on these estimations, it is shown that formation of the C60-Ih isomer is the most probable within the framework of the considered mechanism relative to other C60 isomers. The energetics of sp-atom pair annihilation in the formation of the C60-Ih isomer is also studied via DFT calculations. The advantages of the considered atomistic mechanism of the abundant fullerene isomer formation are discussed.

Detection of nitrofurazone with metal-organic frameworks and reduced graphene oxide composites: insights from molecular dynamics simulations

Two-dimensional metal-organic framework (MOF) composites were produced by incorporating Fe-MOFs into reduced graphene oxide (rGO) nanosheets to form Fe-MOF/rGO composites by hydrothermal synthesis. SEM, TEM, XRD, XPS, and measurements of contact angles were used to characterize the composites. TEM studies revealed that the rod-like-shaped Fe-MOFs were extensively dispersed on the rGO sheets. Incorporating Fe-MOF into rGO significantly improves performance due to the large surface area, chemical stability, and high electrical conductivity. The response signals for the electrochemical sensing performance of Fe-MOF/rGO-modified electrodes to nitrofurazone (NFZ) were significantly enhanced. Differential pulse voltammetry was used to detect the NFZ, and the MOF/rGO sensor possesses a lower detection limit (0.77μM) with two dynamic ranges from 0.6-60 to 128-499.3 μM and high sensitivity (1.909 μA·mM^-1·cm^-2). Moreover, the anti-interference properties of the sensor were quite reproducible and stable. To understand the mechanism responsible for the enhanced sensing performance of the composite, grand canonical Monte Carlo calculations were performed for Fe-MOF/rGO composites with five unit cells of Fe-MOF and four layers of rGO. We attributed the improvement to the fact that the interface between the Fe-MOF and rGO absorbed increased NFZ molecules. The findings reported herein confirm that such Fe-MOF/rGO composites have significantly improved electrochemical performance and practical applicability of sensing nitrofurazone.

Carbon Condensation via [4 + 2] Cycloaddition of Highly Unsaturated Carbon Chains

We present computational studies of reaction pathways for alkyne/polyyne dimerization that represent plausible early steps in mechanisms for carbon condensation. A previous computational study of the ring coalescence and annealing model of C₆₀ formation revealed that a 1,4-didehydrobenzocyclobutadiene intermediate (p-benzyne derivative) has little to no barrier to undergoing an unproductive retro-Bergman cyclization, which brings into question the relevance of that reaction pathway. The current study investigates an alternative model, which proceeds through an initial [4 + 2] cycloaddition instead of a [2 + 2] cycloaddition. In this pathway, the problematic intermediate is avoided, with the reaction proceeding via a (potentially) more kinetically stable tetradehydronaphthalene derivative. The computational studies of the [2 + 2] and [4 + 2] model systems, with increasing alkyne substitutions, reveal that the para-benzyne diradical of the [4 + 2] pathway has a significantly greater barrier to ring opening than the analogous intermediates of the [2 + 2] pathway and that alkyne substitution has little effect on this important barrier. These studies employ spin-flip, time-dependent density functional theory (SF-TDDFT) to provide suitable treatment of open-shell diradical intermediates.

Functionalized Fullerenes and Their Applications in Electrochemistry, Solar Cells, and Nanoelectronics

Carbon-based nanomaterials have rapidly advanced over the last few decades. Fullerenes, carbon nanotubes, graphene and its derivatives, graphene oxide, nanodiamonds, and carbon-based quantum dots have been developed and intensively studied. Among them, fullerenes have attracted increasing research attention due to their unique chemical and physical properties, which have great potential in a wide range of applications. In this article, we offer a comprehensive review of recent progress in the synthesis and the chemical and physical properties of fullerenes and related composites. The review begins with the introduction of various methods for the synthesis of functionalized fullerenes. A discussion then follows on their chemical and physical properties. Thereafter, various intriguing applications, such as using carbon nanotubes as nanoreactors for fullerene chemical reactions, are highlighted. Finally, this review concludes with a summary of future research, major challenges to be met, and possible solutions.

DFTB investigations of the electronic and magnetic properties of fluorographene with vacancies and with adsorbed chemical groups

The electronic and magnetic properties of fluorographene (CF) in the presence of F-vacancy defects and/or chemical groups (-OH, -CN, or -NH₂) were computationally investigated within the framework of the density functional tight-binding (DFTB) method. The current method parameterization allowed us to perform accurate electronic structure calculations (at the ab initio level of many-body methods in the particular case of CF) for hundreds of atoms in the computational cell. We show that the F-vacancy and/or chemical groups influence the magnetic structure, which depends on the number of defects and their distribution between the two sides of the graphene plane. Interestingly, we pointed out a possibility of imprinting local magnetism not only via F-vacancy and -OH combinations, but also using F-vacancies and -CN or -NH₂ groups. In such structures, the magnetic ordering and the total magnetic moments depend on their adsorption sites and their presence in the same or on opposite sides of the graphene plane. We devote particular attention to the interacting chemical group with the F-vacancies. The interaction between the adsorbed chemical group and the unpaired spins associated with the F-vacancies in CF gives rise to interesting magnetic structures. Finally, the zigzag-like direction is shown as the most preferred for the defluorination of CF. Stable ferrimagnetic zigzag chains with interesting properties are considered to be basic magnetic features in perturbed CF. Our work provides new guidelines for engineering multifunctional spintronic components using CF as a base material. We believe, in particular, that the magnetism is predominantly controlled by the F-vacancies, and the ferromagnet can ideally be regulated via the adsorption of a chemical group on a defective CF supercell.

Synthesis of Fullerenes from a Nonaromatic Chloroform through a Newly Developed Ultrahigh-Temperature Flash Vacuum Pyrolysis Apparatus

The flash vacuum pyrolysis (FVP) technique is useful for preparing curved polycyclic aromatic compounds (PAHs) and caged nanocarbon molecules, such as the well-known corannulene and fullerene C₆₀. However, the operating temperature of the traditional FVP apparatus is limited to ~1250 °C, which is not sufficient to overcome the high energy barriers of some reactions. Herein, we report an ultrahigh-temperature FVP (UT-FVP) apparatus with a controllable operating temperature of up to 2500 °C to synthesize fullerene C₆₀ from a nonaromatic single carbon reactant, i.e., chloroform, at 1350 °C or above. Fullerene C₆₀ cannot be obtained from CHCl₃ using the traditional FVP apparatus because of the limitation of the reaction temperature. The significant improvements in the UT-FVP apparatus, compared to the traditional FVP apparatus, were the replacement of the quartz tube with a graphite tube and the direct heating of the graphite tube by impedance heating instead of indirect heating of the quartz tube using an electric furnace. Because of the higher temperature range, UT-FVP can not only synthesize fullerene C₆₀ from single carbon nonaromatic reactants but sublimate some high-molecular-weight compounds to synthesize larger curved PAHs in the future.

Continuous Butadiyne Addition to Propargyl: A Radical-Efficient Pathway for Polycyclic Aromatic Hydrocarbons

Polycyclic aromatic hydrocarbons (PAHs) play a crucial role in soot inception, interstellar evolution, and nanomaterial synthesis. Although several mechanisms, such as hydrogen-abstraction acetylene/vinylacetylene addition, have previously been proposed, PAH formation and growth are not yet fully understood. We propose an alternate PAH growth mechanism wherein propargyl radical reacts with butadiyne to form larger radicals containing newly fused aromatic rings. Butadiyne is an important intermediate in hydrocarbon oxidation and carbon rich stars, while propargyl is one of the most important resonantly stabilized radicals that persists for long times. Our proposed mechanism is validated by quantum chemical calculations, elementary reaction experiments, laminar flame analysis, and kinetic modeling. Our findings challenge the conventional wisdom that radical site regeneration, being central to PAH growth, requires sequential hydrogen elimination and/or abstraction. In our proposed mechanism, PAH growth does not depend on abundant free radical consumption, and could, therefore, help explain carbonaceous nanoparticle coalescence in radical-deficient reaction environments.

High Temperature Accelerated Stone-Wales Transformation and the Threshold Temperature of IPR-C60 Formation

The Stone-Wales bond rotation isomerization of nonicosahedral C₆₀ (C_2v-C₆₀) into isolated-pentagon rule following icosahedral C₆₀ (I_h-C₆₀ or IPR-C₆₀) is a limiting step in the synthesis of I_h-C₆₀. However, extensive previous studies indicate that the potential energy barrier of the Stone-Wales bond rotation is between 6 and 8 eV, extremely high to allow for bond rotation at the temperatures used to produce fullerenes conventionally. This is also despite data indicating a possible fullerene road mechanism that necessitates low-temperature annealing. However, these previous investigations often have limiting factors, such as using the harmonic approximation to determine free energies at high temperatures or considering only the reverse I_h-C₆₀ to C_2v-C₆₀ transition as a basis. Indeed, when the difference in energy between I_h-C₆₀ and C_2v-C₆₀ is accounted for, this barrier is generally reduced by ∼1.5 eV. Thus, utilizing the recently developed density functional tight binding metadynamics (DFTB-MTD) interface, the effects of temperature on the bond rotation in the conversion of C_2v-C₆₀ to I_h-C₆₀ have been investigated. We found that Stone-Wales bond rotations are complex processes with both in-plane and out-of-plane transition states, and which transition path dominates depends on temperature. Our results clearly show that at temperatures of 2000 K, the free energy for a C_2v-C₆₀ to I_h-C₆₀ transition is only ∼4.21 eV and further reduces to ∼3.77 eV at 3000 K. This translates to transition times of ∼971 μs at 2000 K and ∼34 ns at 3000 K, indicating that defect healing is a fast process at temperatures typical of arc jet or laser ablation experiments. Conversely, below ∼2000 K, bond rotation becomes prohibitively slow, putting a lower threshold limit on the temperature of fullerene formation and subsequent annealing.

Simulating fullerene polyhedral formation from planar precursors

The synthesis path of the C₆₀-Buckyball fullerene from a planar precursor developed by Scott et al. [Science, 2002, 295, 5559] is investigated with density functional theory (DFT) methods. Various theoretically possible closing paths are analysed with respect to structural and energetic properties. The initial geometries were obtained by geometric interpolation of a cardboard-like model comprising rigid rings connected by hinges, which were then fully optimized with a selection of DFT-functionals. Analysis of the fully optimised geometries shows remarkable stability of face planarity, bond lengths and bond angles for all studied geometries, indicating soundness of the "cardboard with hinges"-model for approximating reaction paths for molecules of this type. This raises hope for development of a force field description of fullerene precursor molecules that can aid in discovery and analysis of good precursor candidates for rational synthesis of new fullerenes.

Thermochemical stabilities of giant fullerenes using density functional tight binding theory and isodesmic-type reactions

We present a systematic assessment of the density functional tight binding (DFTB) method for calculating heats of formation of fullerenes with isodesmic-type reaction schemes. We show that DFTB3-D/3ob can accurately predict Δ_f H values of the 1812 structural isomers of C₆₀ , reproduce subtle trends in Δ_f H values for 24 isolated pentagon rule (IPR) isomers of C₈₄ , and predict Δ_f H values of giant fullerenes that are in effectively exact agreement with benchmark DSD-PBEP86/def2-QZVPP calculations. For fullerenes up to C₃₂₀ , DFTB Δ_f H values are within 1.0 kJ mol^-1 of DSD-PBEP86/def2-QZVPP values per carbon atom, and on a per carbon atom basis DFTB3-D/3ob yields exactly the same numerical trend of (Δ_f H [per carbon] = 722n^-0.72  + 5.2 kJ mol^-1 ). DFTB3-D/3ob is therefore an accurate replacement for high-level DHDFT and composite thermochemical methods in predicting of thermochemical stabilities of giant fullerenes and analogous nanocarbon architectures.

Surface specificity and mechanistic pathway of de-fluorination of C60F48 on coinage metals

We provide experimental and theoretical understanding on fundamental processes taking place at room temperature when a fluorinated fullerene dopant gets close to a metal surface. By employing scanning tunneling microscopy and photoelectron spectroscopies, we demonstrate that the on-surface integrity of C60F48 depends on the interaction with the particular metal it approaches. Whereas on Au(111) the molecule preserves its chemical structure, on more reactive surfaces such as Cu(111) and Ni(111), molecules interacting with the bare metal surface lose the halogen atoms and transform to C60. Though fluorine-metal bonding can be detected depending on the molecular surface density, no ordered fluorine structures are observed. We show the implications of the metal-dependent de-fluorination in the electronic structure of the molecules and the energy alignment at the molecule-metal interface. Molecular dynamics simulations with ReaxFF reactive force field corroborate the experimental facts and provide a detailed mechanistic picture of the surface-induced de-fluorination, which involves the rotation of the molecule on the surface. Outstandingly, a thermodynamic analysis indicates that the effect of the metal surface is lowering and diminishing the energy barrier for C-F cleave, demonstrating the catalytic role of the surface. The present study contributes to in-depth knowledge of the mechanisms that affect the degree of stability of chemical species on surfaces, which is essential to advance our understanding of the chemical reactivity of metals and their role in on-surface chemical reactions.

Fullerene Thermochemical Stability: Accurate Heats of Formation for Small Fullerenes, the Importance of Structural Deformation on Reactivity, and the Special Stability of C60

We have used quantum chemistry computations, in conjunction with isodesmic-type reactions, to obtain accurate heats of formation (HoFs) for the small fullerenes C₂₀ (2358.2 ± 8.0 kJ mol^-1), C₂₄ (2566.2 ± 7.6), and the lowest-energy isomers of C₃₂ (2461.1 ± 15.4), C₄₂ (2629.0 ± 20.5), and C₅₄ (2686.2 ± 25.3). As part of this endeavor, we have also obtained accurate HoFs for several medium-sized molecules, namely 216.6 ± 1.4 for fulvene, 375.5 ± 1.5 for pentalene, 670.8 ± 2.9 for acepentalene, and 262.7 ± 2.5 for acenaphthylene. We combine the energies of the small fullerenes and previously obtained energies for larger fullerenes (from C₆₀ to C₆₀₀₀) into a full picture of fullerene thermochemical stability. In general, the per-carbon energies can be reasonably approximated by the "R+D" model that we have previously developed [Chan et al. J. Chem. Theory Comput. 2019, 15, 1255-1264], which takes into account Resonance and structural Deformation factors. In a case study on C₅₄, we find that most of the high-deformation-energy atoms correspond to the sites of the C-Cl bond in the experimentally captured C₅₄Cl₈. In another case study, we find that C₆₀ has the lowest value for the maximum local-deformation energy when compared with similar-sized fullerenes, which is consistent with its "special stability". These results are indicative of structural deformation playing an important role in the reactivity of fullerenes.

Skeletal Rearrangements of the C240 Fullerene: Efficient Topological Descriptors for Monitoring Stone–Wales Transformations

Stone–Wales rearrangements of the fullerene surface are an uncharted field in theoretical chemistry. Here, we study them on the example of the giant icosahedral fullerene C240 to demonstrate the complex chemical mechanisms emerging on its carbon skeleton. The Stone–Wales transformations of C240 can produce the defected isomers containing heptagons, extra pentagons and other unordinary rings. Their formations have been described in terms of (i) quantum-chemically calculated energetic, molecular, and geometric parameters; and (ii) topological indices. We have found the correlations between the quantities from the two sets that point out the role of long-range topological defects in governing the formation and the chemical reactivity of fullerene molecules.

Room-temperature chemical synthesis of C2

Diatomic carbon (C2) is historically an elusive chemical species. It has long been believed that the generation of C2 requires extremely high physical energy, such as an electric carbon arc or multiple photon excitation, and so it has been the general consensus that the inherent nature of C2 in the ground state is experimentally inaccessible. Here, we present the chemical synthesis of C2 from a hypervalent alkynyl-λ3-iodane in a flask at room temperature or below, providing experimental evidence to support theoretical predictions that C2 has a singlet biradical character with a quadruple bond, thus settling a long-standing controversy between experimental and theoretical chemists, and that C2 serves as a molecular element in the bottom-up chemical synthesis of nanocarbons such as graphite, carbon nanotubes, and C60.

The Effect of the Polyaromatic Hydrocarbon in the Formation of Fullerenes

Tremendous advances in nanoscience have been made since the discovery of fullerenes. However, the short timescale of the growth process and high-energy conditions of synthesis result in severe constraints to investigation of the mechanism of fullerene formation. In this work, we attempted to reveal the formation process by analyzing the variation in the yield of fullerenes under different conditions. Experiments and theoretical analysis show that the formation of fullerenes could be affected by the addition of polycyclic aromatic compounds. It is proposed that the formation of C₆₀ during arc-discharge synthesis is fragment assembling, while the yield of C_2m (m=35, 38, 39) is strongly enhanced by building-block splicing. In addition, several features of the building blocks are put forward to predict the extent of their influence to the formation of larger fullerenes C_2n (n≥42). This work not only provides essential insight into the formation process of fullerenes, but more importantly also paves the way to improving the yield of larger fullerenes selectively.

Discrete Solvent Reaction Field Calculations for One- and Two-Photon Absorptions of Solution-Phase Dimethylaminonitrostilbene Molecule

Based on the configurations generated by molecular dynamics (MD) simulations using the on-the-fly density-functional tight-bonding (DFTB) force field, we investigated performance of the discrete solvent reaction field (DRF) model coupled to time-dependent density functional theory (TD-DFT) for solvatochromic effect of one- and two-photon absorption phenomena. Dimethylaminonitrostilbene (DANS) molecule solvated in chloroform, dichloromethane, and dimethyl sulfoxide solvents was selected as a model system for our research purpose. For every selected MD/DFTB configuration, within the context of the DRF, solute molecule is represented by TD-DFT and solvent molecules are described by atomic charges and polarizabilities. The calculated one-photon absorption energies reproduce well the positive solvatochromic behavior of solvated DANS and are in good agreement with available experimental data. For the two-photon absorption cross section, even though our approach overshot the experimental data by about 20% in absolute magnitude, experimentally observed solvatochromic change was captured qualitatively in this work. At last, we examined the contributions of atomic charges and polarizabilities of solvent molecules to the solvatochromic shifts of properties of interest.

Nanostructure of Gasification Charcoal (Biochar)

In this work, we investigate the molecular composition and nanostructure of gasification charcoal (biochar) by comparing it with heat-treated fullerene arc-soot. Using ultrahigh resolution Fourier transform ion-cyclotron resonance and laser desorption ionization time-of-flight mass spectrometry, Raman spectroscopy, and high resolution transmission electron microscopy we analyzed charcoal of low tar content obtained from gasification. Mass spectrometry revealed no magic number fullerenes such as C₆₀ or C₇₀ in the charcoal. The positive molecular ion m/ z 701, previously considered a graphitic part of the nanostructure, was found to be a breakdown product of pyrolysis and not part of the nanostructure. A higher mass distribution of ions similar to that found in thermally treated fullerene soot indicates that they share a nanostructure. Recent insights into the formation of all carbon fullerenes reveal that conditions in charcoal formation are not optimal for the formation of fullerenes, but instead, curved carbon structures coalesce into fulleroid-like structures. Microscopy and spectroscopy support such a stacked, fulleroid-like nanostructure, which was explored using reactive molecular dynamics simulations.

A Simple Model for Relative Energies of All Fullerenes Reveals the Interplay between Intrinsic Resonance and Structural Deformation Effects in Medium-Sized Fullerenes

Fullerenes are sheets of sp² carbon atoms wrapped around to form spheres. With this simple consideration, we have in the present study devised and (with over 3600 DFT data points) successfully validated a simple model, termed R+D, for estimating the relative energies of fullerenes. This model contains a resonance component to account for the intrinsic differences between the π-energies of different fullerenes, and a deformation component for treating the distortions from planarity. Notably, we find that both terms (and they alone) are required to obtain good relative energies, which lends support to the formulation of the R+D model. An interesting finding is that for some medium-sized IPR fullerenes, their isomers show similar variations in the two components. We deduce that these fullerenes may represent a good opportunity for tuning molecular properties for practical applications. We hope that the promising results of the present study will encourage further investigations into fullerenes from a fundamental perspective.

Aromaticity, Coulomb repulsion, π delocalization or strain: who is who in endohedral metallofullerene stability?

Three different models for endohedral metallofullerene structure prediction are compared, revealing the physical origin of the stability of these compounds.

Performance of Density-Functional Tight-Binding in Comparison to Ab Initio and First-Principles Methods for Isomer Geometries and Energies of Glucose Epimers in Vacuo and Solution

Density functional theory (DFT) is a widely used methodology for the computation of molecular and electronic structure, and we confirm that B3LYP and the high-level ab initio G3B3 method are in excellent agreement for the lowest-energy isomers of the 16 glucose epimers. Density-functional tight-binding (DFTB) is an approximate version of DFT with typically comparable accuracy that is 2 to 3 orders of magnitude faster, therefore generally very suitable for processing large numbers of complex structures. Conformational isomerism in sugars is well known to give rise to a large number of isomer structures. On the basis of a comprehensive study of glucose epimers in vacuo and aqueous solution, we found that the performance of DFTB is on par to B3LYP in terms of geometrical parameters excluding hydrogen bonds and isomer energies. However, DFTB underestimates both hydrogen bonding interactions as well as torsional barriers associated with rotations of the hydroxy groups, resulting in a counterintuitive overemphasis of hydrogen bonding in both gas phase as well as in water. Although the associated root mean squared deviation from B3LYP within epimer isomer groups is only on the order of 1 kcal/mol, this deviation affects the correct assignment of major isomer ordering, which span less than 10 kcal/mol. Both second- as well as third-order DFTB methods are exhibiting similar deviations from B3LYP. Even after the inclusion of empirical dispersion corrections in vacuum, these deviations remain for a large majority of isomer energies and geometries when compared to dispersion-corrected B3LYP.

Simulations of the synthesis of boron-nitride nanostructures in a hot, high pressure gas volume

We performed nanosecond timescale computer simulations of clusterization and agglomeration processes of boron nitride (BN) nanostructures in hot, high pressure gas, starting from eleven different atomic and molecular precursor systems containing boron, nitrogen and hydrogen at various temperatures from 1500 to 6000 K. The synthesized BN nanostructures self-assemble in the form of cages, flakes, and tubes as well as amorphous structures. The simulations facilitate the analysis of chemical dynamics and we are able to predict the optimal conditions concerning temperature and chemical precursor composition for controlling the synthesis process in a high temperature gas volume, at high pressure. We identify the optimal precursor/temperature choices that lead to the nanostructures of highest quality with the highest rate of synthesis, using a novel parameter of the quality of the synthesis (PQS). Two distinct mechanisms of BN nanotube growth were found, neither of them based on the root-growth process. The simulations were performed using quantum-classical molecular dynamics (QCMD) based on the density-functional tight-binding (DFTB) quantum mechanics in conjunction with a divide-and-conquer (DC) linear scaling algorithm, as implemented in the DC-DFTB-K code, enabling the study of systems as large as 1300 atoms in canonical NVT ensembles for 1 ns time.

A review of theoretical study of graphene chemical vapor deposition synthesis on metals: nucleation, growth, and the role of hydrogen and oxygen

Graphene has attracted intense research interest due to its extraordinary properties and great application potential. Various methods have been proposed for the synthesis of graphene, among which chemical vapor deposition has drawn a great deal of attention for synthesizing large-area and high-quality graphene. Theoretical understanding of the synthesis mechanism is crucial for optimizing the experimental design for desired graphene production. In this review, we discuss the three fundamental steps of graphene synthesis in details, i.e. (1) decomposition of carbon feedstocks and formation of various active carbon species, (2) nucleation, and (3) attachment and extension. We provide a complete scenario of graphene synthesis on metal surfaces at atomistic level by means of density functional theory, molecular dynamics (MD), Monte Carlo (MC) and their combination and interface with other simulation methods such as quantum mechanical molecular dynamics, density functional tight binding molecular dynamics, and combination of MD and MC. We also address the latest investigation of the influences of the hydrogen and oxygen on the synthesis and the quality of the synthesized graphene.

Revisit of the Saito-Dresselhaus-Dresselhaus C2 ingestion model: on the mechanism of atomic-carbon-participated fullerene growth

We introduce a mechanistic study based on a controversial fullerene bottom-up growth model proposed by R. Saito, G. Dresselhaus, and M. S. Dresselhaus. The so-called SDD C₂ addition model has been dismissed as chemically inadmissible but here we prove that it is feasible via successive atomic-carbon-participated addition and migration reactions. Kinetic calculations on the formation of isolated pentagon rule (IPR)-obeying C₇₀ and Y₃N@C₈₀ are carried out by employing the SDD model for the first time. A stepwise mechanism is proposed with a considerably low barrier of ca. 2 eV which is about 3 eV lower than a conventional isomerization-containing fullerene growth pathway.

Transformation of doped graphite into cluster-encapsulated fullerene cages

An ultimate goal in carbon nanoscience is to decipher formation mechanisms of highly ordered systems. Here, we disclose chemical processes that result in formation of high-symmetry clusterfullerenes, which attract interest for use in applications that span biomedicine to molecular electronics. The conversion of doped graphite into a C₈₀ cage is shown to occur through bottom-up self-assembly reactions. Unlike conventional forms of fullerene, the iconic Buckminsterfullerene cage, I_h-C₆₀, is entirely avoided in the bottom-up formation mechanism to afford synthesis of group 3-based metallic nitride clusterfullerenes. The effects of structural motifs and cluster–cage interactions on formation of compounds in the solvent-extractable C₇₀–C₁₀₀ region are determined by in situ studies of defined clusterfullerenes under typical synthetic conditions. This work establishes the molecular origin and mechanism that underlie formation of unique carbon cage materials, which may be used as a benchmark to guide future nanocarbon explorations.

An understanding of how caged carbon materials self-assemble from doped graphite is a long-standing challenge. Here, the authors show that distinct bottom-up processes lead to the synthesis of high-symmetry clusterfullerenes.

Ion-induced molecular growth in clusters of small hydrocarbon chains

We report on studies of collisions between 3 keV Ar⁺ projectile ions and neutral targets of isolated 1,3-butadiene (C₄H₆) molecules and cold, loosely bound clusters of these molecules. We identify molecular growth processes within the molecular clusters that appears to be driven by knockout processes and that could result in the formation of (aromatic) ring structures. These types of reactions are not unique to specific projectile ions and target molecules, but will occur whenever atoms or ions with suitable masses and kinetic energies collide with aggregates of matter, such as carbonaceous grains in the interstellar medium or aerosol nanoparticles in the atmosphere.

Formation of endohedral metallofullerene (EMF) ions of MnC2m+ (M = La, Y, n ≤ 6, 50 ≤ 2m ≤ 194) in the laser ablation process with graphene as precursor

RATIONALE

Graphene has proved to be a versatile precursor for the generation of mono-endohedral metallofullerenes (EMFs) of alkali and alkali earth metal atoms. Thus, it is important to investigate whether graphene is also a good precursor for the generation of EMFs including multi-metal atoms or some transition metal atoms.

METHODS

Experiments were performed on a 7.0 T Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. Laser ablation mass spectra were obtained by laser irradiation on mixtures of graphene and LaCl₃ , YCl₃ or HAuCl₄ on a stainless steel plate. A 355 nm Nd:YAG laser with a typical energy of 2.5 mJ/pulse was employed in these experiments. Tandem mass spectrometry was performed by on-resonance excitation collision-induced dissociation (CID) for selected single isotopic peaks.

RESULTS

Metallofullerene La1-6C2m+/ Y1-6C2m+ ions were observed in the laser ablation mass spectra of mixtures of graphene and LaCl₃ /YCl₃ . CID and reactant gas experiments on selected ions showed that these ions had endohedral structures. Considering the experimental condition applied here, a mechanism based on a top-down process for the generation of the EMF ions is suggested. Further experiments on HAuCl₄ showed some weak signals of AuC2m+ for the first time.

CONCLUSIONS

Mono- and multi-EMF ions of MnC2m+ (M = La, Y, n = 1-6, 50 ≤ 2 m ≤ 194) were generated by laser ablation of graphene and the corresponding metal salts and detected by FT-ICR mass spectrometry. The results have shed light on the formation and generation of some new EMF species. Copyright © 2017 John Wiley & Sons, Ltd.

Formation of Nickel Clusters Wrapped in Carbon Cages: Toward New Endohedral Metallofullerene Synthesis

Despite the high potential of endohedral metallofullerenes (EMFs) for application in biology, medicine and molecular electronics, and recent efforts in EMF synthesis, the variety of EMFs accessible by conventional synthetic methods remains limited and does not include, for example, EMFs of late transition metals. We propose a method in which EMF formation is initiated by electron irradiation in aberration-corrected high-resolution transmission electron spectroscopy (AC-HRTEM) of a metal cluster surrounded by amorphous carbon inside a carbon nanotube serving as a nanoreactor and apply this method for synthesis of nickel EMFs. The use of AC-HRTEM makes it possible not only to synthesize new, previously unattainable nanoobjects but also to study in situ the mechanism of structural transformations. Molecular dynamics simulations using the state-of-the-art approach for modeling the effect of electron irradiation are performed to rationalize the experimental observations and to link the observed processes with conditions of bulk EMF synthesis.

Comprehensive theoretical study of all 1812 C60 isomers

All 1812 C₆₀ isomers are investigated with high-level quantum chemical methods to benchmark semiempirical approaches and find appropriate stability criteria.

Generalized structural motif model for studying the thermodynamic stability of fullerenes: from C60to graphene passing through giant fullerenes

A generalized motif model to describe the stability of neutral fullerenes, covering the full range of cage sizes, starting from C₆₀, going through giant fullerenes, and ultimately leading to graphene.

A theoretical study of complexes between fullerenes and concave receptors with interest in photovoltaics

The proper combination of host and guest allows controlling the stability and charge transfer capability of fullerene–concave receptor complexes.

Probing Distinct Fullerene Formation Processes from Carbon Precursors of Different Sizes and Structures

Fullerenes, cage-structured carbon allotropes, have been the subject of extensive research as new materials for diverse purposes. Yet, their formation process is still not clearly understood at the molecular level. In this study, we performed laser desorption ionization-ion mobility-mass spectrometry (LDI-IM-MS) of carbon substrates possessing different molecular sizes and structures to understand the formation process of fullerene. Our observations show that the formation process is strongly dependent on the size of the precursor used, with small precursors yielding small fullerenes and large graphitic precursors generally yielding larger fullerenes. These results clearly demonstrate that fullerene formation can proceed via both bottom-up and top-down processes, with the latter being favored for large precursors and more efficient at forming fullerenes. Furthermore, we observed that specific structures of carbon precursors could additionally affect the relative abundance of C60 fullerene. Overall, this study provides an advanced understanding of the mechanistic details underlying the formation processes of fullerene.

Self-assembly of endohedral metallofullerenes: a decisive role of cooling gas and metal-carbon bonding

The endohedral metallofullerene (EMF) self-assembly process in Sc/carbon vapor in the presence and absence of an inert cooling gas (helium) is systematically investigated using quantum chemical molecular dynamics simulations. It is revealed that the presence of He atoms accelerates the formation of pentagons and hexagons and reduces the size of the self-assembled carbon cages in comparison with analogous He-free simulations. As a result, the Sc/C/He system simulations produce a larger number of successful trajectories (i.e. leading to Sc-EMFs) with more realistic cage-size distribution than simulations of the Sc/C system. The main Sc encapsulation mechanism involves nucleation of several hexagons and pentagons with Sc atoms already at the early stages of carbon vapor condensation. In such proto-cages, both Sc-C σ-bonds and coordination bonds between Sc atoms and the π-system of the carbon network are present. Sc atoms are thus rather labile and can move along the carbon network, but the overall bonding is sufficiently strong to prevent dissociation even at temperatures around 2000 kelvin. Further growth of the fullerene cage results in the encapsulation of one or two Sc atoms within the fullerene. In agreement with experimental studies, an extension of the simulations to Fe and Ti as the metal component showed that Fe-EMFs are not formed at all, whereas Ti is prone to form Ti-EMFs with small cage sizes, including Ti@C28-Td and Ti@C30-C2v(3).

QM/MD studies on graphene growth from small islands on the Ni(111) surface

Quantum chemical molecular dynamics simulations of graphene growth from small island precursors in different carbon nucleation densities on the Ni(111) surface at high temperatures have been conducted. The results indicate that small islands are not static, i.e. lateral diffusion and vertical fluctuation are frequently observed. In the case of low carbon nucleation density, carbon atoms or small carbon patches diffuse and attach to the edge of the nuclei to expand the size of the growing carbon network. The growth of graphene precursors is accompanied by the corresponding changes in the bonding of nickel atoms with the precipitation of subsurface carbon atoms. This is because the carbon-carbon interaction is stronger than the nickel-carbon interaction. In the case of high carbon nucleation densities, the dominant ripening mechanism depends on different growth stages. In the initial stage, the coalescence of carbon islands takes place via the Smoluchowski ripening mechanism. In the later stage the Smoluchowski ripening process is damped owing to the higher diffusion barrier of larger clusters and the restriction of movement by self-assembled nickel step edges. The cross-linking mechanism eventually takes over by the coalescence of extended polyyne chains between graphene islands. In either case, the Ostwald ripening process is not found in our molecular dynamics simulations due to the stability of carbon-carbon bonds within the islands. These investigations should be instructive to the control of graphene growth in experiments.

Possible mechanism of BN fullerene formation from a boron cluster: Density-functional tight-binding molecular dynamics simulations

We simulate the formation of a BN fullerene from an amorphous B cluster at 2000 K by quantum mechanical molecular dynamics based on the density-functional tight-binding method. We run 30 trajectories 200 ps in length, where N atoms are supplied around the target cluster, which is initially an amorphous B36 cluster. Most of the incident N atoms are promptly incorporated into the target cluster to form B-N-B bridges or NB3 pyramidal local substructures. BN fullerene formation is initiated by alternating BN ring condensation. Spontaneous atomic rearrangement and N2 dissociation lead to the construction of an sp(2) single-shelled structure, during which the BN cluster undergoes a transition from a liquid-like to a solid-like state. Continual atomic rearrangement and sporadic N2 dissociation decrease the number of defective rings in the BN cluster and increase the number of six-membered rings, forming a more regular shell structure. The number of four-membered rings tends to remain constant, and contributes to more ordered isolated-tetragon-rule ring placement.

Carbon rings: a DFT study on geometry, aromaticity, intermolecular carbon–carbon interactions and stability

Strong complimentary electrostatic interactions lead to large number of C⋯C interactions in doubly antiaromatic C_4N dimers while the delocalized system of π electrons in doubly aromatic C_4N+2 dimers suggest very few intermolecular C⋯C interactions.

An adaptive quantum mechanics/molecular mechanics method for the infrared spectrum of water: incorporation of the quantum effect between solute and solvent

Quantum effects in solute–solvent interactions, such as the many-body effect and the dipole-induced dipole, are known to be critical factors influencing the infrared spectra of species in the liquid phase.