The Formation and Transformation of Low-Dimensional Carbon Nanomaterials by Electron Irradiation

Low-dimensional materials based on graphene or graphite show a large variety of phenomena when they are subjected to irradiation with energetic electrons. Since the 1990s, electron microscopy studies, where a certain irradiation dose is unavoidable, have witnessed unexpected structural transformations of graphitic nanoparticles. It is recognized that electron irradiation is not only detrimental but also bears considerable potential in the formation of new graphitic structures. With the availability of aberration-corrected electron microscopes and the discovery of techniques to produce monolayers of graphene, detailed insight into the atomic processes occurring during electron irradiation became possible. Threshold energies for atom displacements are determined and models of different types of lattice vacancies are confirmed experimentally. However, experimental evidence for the configuration of interstitial atoms in graphite or adatoms on graphene remained indirect, and the understanding of defect dynamics still depends on theoretical concepts. This article reviews irradiation phenomena in graphene- or graphite-based nanomaterials from the scale of single atoms to tens of nanometers. Observations from the 1990s can now be explained on the basis of new results. The evolution of the understanding during three decades of research is presented, and the remaining problems are pointed out.

Energy Landscapes of Carbon Clusters from Tight-Binding Quantum Potentials

We calculate transformation pathways between fullerene and octahedral carbon clusters and between a buckyball and its bowl-shaped isomer. The energies and gradients are provided by efficient tight-binding potentials, which are interfaced to our Energy Landscape exploration software. From the global energy landscape, we extract the mechanistic and kinetic parameters as a function of temperature, and compare our results to selected density functional theory (DFT) (PBE/cc-pVTZ) benchmarks. Infrared spectra are calculated to provide data for experimental identification of the clusters and differentiation of their isomers. Our results suggest that the formation of buckyballs from a buckybowl will be suppressed at elevated temperatures (above around 5250 K) due to entropic effects, which may provide useful insight into the detection of cosmic fullerenes.

Bootstrap Embedding For Large Molecular Systems

Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [J. Chem. Theory Comput. 2019, 15, 4497-4506; J. Phys. Chem. Lett. 2019, 10, 6368-6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O(N³) for the integral transform and O(N) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.

Fullerene and nanotube growth: new insights using first principles and molecular dynamics

Shortly after the discovery of fullerenes, many researchers pointed out that carbon nanotubes could be considered as elongated fullerenes. However, the detailed formation mechanism for both structures has been a topic of debate for several years, and consequently it has been difficult to draw a clear connection between the two systems. While the synthesis conditions appear to be different for both fullerenes and nanotubes, here, we demonstrate that it is highly likely that, at an initial growth stage, single-walled carbon nanotubes begin to grow from a hemisphere-like fullerene cap. More importantly, by analysing the minimum-energy path, it is shown that the insertion of C2 fragments drives the transformation of this fullerene cap into an elongated structure that leads to the formation of very short carbon nanotubes.This article is part of the themed issue 'Fullerenes: past, present and future, celebrating the 30th anniversary of Buckminster Fullerene'.

STRUCTURAL, ELECTRONIC, MAGNETIC, AND OPTICAL PROPERTIES OF 3d TRANSITION METAL ENDOHEDRAL M@Ge12H12(M=Sc–Ni) CLUSTERS

The fulerine- Ni @ Ge 12 H 12 structure, which composes of four pentagons and four rhombi and is like a fullerene, has a closed-shell electronic structure, the largest HOMO–LUMO energy gap, the highest vertical ionization potential, and the lowest vertical electron affinity. All of these properties are characteristic of a magic cluster, therefore, we strongly suggest fulerine- Ni @ Ge 12 H 12 should be a magic cluster and promising as building blocks in developing cluster-assembled nanomaterials. This can be interpreted by the weak interaction between Ni and the cage together with the transference of two electrons from the 4s orbital to the 3d orbital of Ni . The magnetic moment of fulerine- M @ Ge 12 H 12( M = Sc – Ni ) varies from 0 to 3 μB, implying they have potential applications in developing new nanomaterials with tunable magnetic properties. The calculated TDDFT optical properties of fulerine- M @ Ge 12 H 12( M = Sc – Ni ) can be tuned broadly in the ultraviolet–visible region. This is very important for optoelectronic applications.

Unfolding the fullerene: nanotubes, graphene and poly-elemental varieties by simulations

Recent research progress in nanostructured carbon has built upon and yet advanced far from the studies of more conventional carbon forms such as diamond, graphite, and perhaps coals. To some extent, the great attention to nano-carbons has been ignited by the discovery of the structurally least obvious, counterintuitive, small strained fullerene cages. Carbon nanotubes, discovered soon thereafter, and recently, the great interest in graphene, ignited by its extraordinary physics, are all interconnected in a blend of cross-fertilizing fields. Here we review the theoretical and computational models development in our group at Rice University, towards understanding the key structures and behaviors in the immense diversity of carbon allotropes. Our particular emphasis is on the role of certain transcending concepts (like elastic instabilities, dislocations, edges, etc.) which serve so well across the scales and for chemically various compositions.

Ab initio calculation of carbon clusters. II. Relative stabilities of fullerene and nonfullerene C24

Chemical stabilities of six low-energy isomers of C24 derived from global-minimum search are investigated. The six isomers include one classical fullerene (isomer 1) whose cage is composed of only five- and six-membered rings (56-MRs), three nonclassical fullerene structures whose cages contain at least one four-membered ring (4-MR), one plate, and one monocyclic ring. Chemical and electronic properties of the six C24 isomers are calculated based on a density-functional theory method (hybrid PBE1PBE functional and cc-pVTZ basis set). The properties include the nucleus-independent chemical shifts (NICS), singlet-triplet splitting, electron affinity, ionization potential, and gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital (HOMO-LUMO) gap. The calculation suggests that the neutral isomer 2, a nonclassical fullerene with two 4-MRs, may be more chemically stable than the classical fullerene (isomer 1). Analyses of molecular orbital NICS show that the incorporations of 4-MRs into the cage considerably reduce paratropic contributions from HOMO, HOMO-1, and HOMO-2, which are mainly responsible for the sign change in NICS from positive for isomer 1 (42) to negative (-19) for isomer 2, although C24 clusters satisfy neither 4N+2 nor 2(N+1)2 aromaticity rule. Anion photoelectron spectra of four cage isomers, one plate, one monocyclic ring, and one tadpole isomer, as well as three bicyclic ring isomers are calculated. The simulated photoelectron spectra of mono- and bicyclic rings (with C1 symmetry) appear to match the measured HOMO-LUMO gap (between the first and second band in the experimental spectra) [S. Yang et al., Chem. Phys. Lett. 144, 431 (1988)]. Nevertheless, the nonclassical fullerene isomers 3 and 4 apparently also match the measured vertical detachment energy (2.90 eV) reasonably well. These results suggest possible coexistence of nonclassical fullerene isomers with the mono- and bicyclic ring isomers of C24(-) under the experimental conditions.

Structure and electronic properties of fullerenes C(52)q+: is C(52)2+ an exception to the pentagon adjacency penalty rule?

The structure, vibrational spectra and electronic properties of the neutral, singly and doubly charged C52 fullerenes were studied by means of the Hartree-Fock method and density functional theory. Different isomers were considered, in particular those with the lowest possible number (five or six) of adjacent pentagons, and an isomer with a four-atom ring. For neutral and singly charged species, the most stable isomer is that with the lowest number of adjacent pentagons, namely five. However, for C(52)2+, the most stable structure has six adjacent pentagons. This finding, which contradicts the pentagon adjacency penalty rule, is a consequence of complete filling of the HOMO pi shell and the near-perfect sphericity of the most stable isomer. The simulated vibrational spectra show important differences in the positions and intensities of the vibrations for the different isomers.

Bond order bond polarizability model for fullerene cages and nanotubes

It is still a challenge to accurately calculate the polarizabilities of large fullerene cages and nanotubes. In this paper, a simple bond order bond polarizability relationship for carbon was found, which allowed us to apply the bond polarizability model to any pentagon isolation rule (PIR) fullerene (cage or nanotube). Following this approach, the following simple equation, alpha=1.262n, was obtained relating the static dipole polarizability (alpha) of PIR fullerenes (cages or closed nanotubes) to their number (n) of carbon atoms. Furthermore, it was shown that the polarizabilities of C60 and C70, calculated on the basis of this model, are in excellent agreement with those obtained experimentally and by density-functional theory calculations.

Density functional theory calculations for endohedral complexes of non-πC60H60 cage with small guest molecules

Hybrid density-functional theory (B3LYP) calculations were carried out to determine the structures and energies of endohedral complexes of non-pi C(60)H(60) with H(2), CO, and LiH. It was demonstrated that the endohedral complexes of C(60)H(60) with the above three guest molecules are more stable than the corresponding complexes with C(60). Furthermore, the interaction between C(60)H(60) and the inside H(2) or CO is negligible, but the formation of the LiH-C(60)H(60) complex is exothermic with a stabilization energy of -6.0 kcal/mol. While the bond lengths of H(2) and CO changed a little when placed inside the cages, that of the LiH molecule increased and decreased inside C(60)H(60) and C(60), respectively.

Endohedral Chemistry of C60-Based Fullerene Cages

Fullerenes have unique chemistry owing to their cage structure, their richness in pi-electrons, and their large polarizabilities. They can trap atoms and small molecules to generate endohedral complexes as superconductors, drug carriers, molecular reactors, and ferroelectric materials. An important goal is to develop effective methods that can affect the behavior of the atoms and small molecules trapped inside the cage. In this paper, the quantum chemical density functional theory was employed to demonstrate that the stability and position of a guest molecule inside the C60 cage can be changed, and its orientation controlled, by modifying the C60 cage shell. The outside attachment of two hydrogen atoms to two adjacent carbon atoms located between two six-membered rings of the C60 cage affects the orientation of the LiF molecule inside and increases the stability of LiF inside the cage by 45%. In contrast, when 60 hydrogen atoms were attached to the outside surface of the C60 cage, thus transforming all C=C double bonds into single bonds, the stability of the LiF inside was reduced by 34%. If two adjacent carbon atoms were removed from C60, the stability of LiF inside this defect C60 was reduced by 41%.

Coulomb stability limit of highly charged Cq+60 fullerenes

We perform density functional theory calculations of Cq+60 and Cq+58 fullerenes as well as of transition states related to dissociation of Cq+60 into Cq-s(+)58 and C(s+)2 for charges q=0-14. For q< or =8, the lowest fission barrier corresponds to C+2 emission, whereas, for higher q values, the lowest barrier corresponds to the emission of two atomic fragments [2C]s+ with s> or =2. We conclude that the Coulomb stability limit corresponds to q=14.

Selective cap opening in carbon nanotubes driven by laser-induced coherent phonons

We demonstrate the possibility of a selective nonequilibrium cap opening of carbon nanotubes as a response to femtosecond laser excitation. By performing molecular dynamics simulations based on a microscopic electronic model we show that the laser-induced ultrafast structural changes differ dramatically from the thermally induced dimer emission. Ultrafast bond weakening and simultaneous excitation of two coherent phonon modes of different frequencies, localized in the spherical caps and cylindrical nanotube body, are responsible for the selective cap opening.

Dynamic topology of fullerene coalescence

Atomic rearrangements leading to the coalescence of fullerene cages are revealed by topological analysis. Paths found for nanotubes and C(60) consist exclusively of Stone-Wales bond rotations. Computed intermediate states show gradual evolution from separate clusters to completely fused coherent units. Molecular dynamics simulations follow the predicted routes, overcome the nucleation barrier, and reach a final annealed state. While the overall behavior resembles macroscopic sintering, the nanoscale neck undergoes quantized changes in diameter and crystalline order.

Mode-actions of the Na(+)-Ca2+ exchanger: from genes to mechanisms to a new strategy in brain disorders

Mode-actions of the Na(+)-Ca2+ exchanger from genes to mechanisms to a new strategy for brain disorders were comparatively studied in oxidative stress. In transfected Chinese hamster ovary (CHO) cells steadily expressing the Na(+)-Ca2+ exchanger's gene, Ca(2+)-efflux via an active mode of the Na(+)-Ca2+ exchanger was elicited by hydrogen peroxide (H2O2) after preincubation of the cell with a Ca(2+)-free medium, whereas Ca(2+)-influx via a reverse mode of the Na(+)-Ca2+ exchanger was dramatically evoked by H2O2 after preincubation of the cell with a Ca2+ medium, as a prelude to neuronal death. According to [45Ca2+] uptake of transfected CHO cells at given time intervals or extracellular Na+[Na+]o gradients, hyperbola, logarithmic and sigmoid curve equations of the Na(+)-Ca2+ exchanger's mode-actions were respectively defined in the absence and the presence of H2O2. The Na(+)-Ca2+ exchanger's conformational transition in oxidative stress was dominated by adenosine triphosphate (ATP)-dependent cytoskeletal redox modification, cation-pi interactions and secondary Ca2+ activation. These mechanisms were used to generate an intracellulary distributed tetra-cluster (named VISA931) for rescuing G-protein agonist-sensitive signal transduction and cortico-cerebral somatosensory evoke potential (SEP) from oxidation via activating forward operation of the Na(+)-Ca2+ exchanger, the beta-adrenergic and the P2-purinergic receptors, blocking Ca2+ influx and catalyzing the dismutation of superoxide anions (O2-.) to H2O2. In conclusion, knowledge-based drug design is a new strategy for developing promising candidates of neuroprotective agents.