Received 31 December 2015; accepted 3 April 2016; published 6 April 2016
Graphene is an allotrope of carbon consisting of a single layer of carbon atoms connected by sp2 hybridized bonds   . It has been investigated due to its unique electrical, thermal and mechanical properties  -  . An important application of graphene is as a promising replacement for tin-doped indium oxide (ITO) as the next-generation large-area transparent conducting electrodes (TCE)   . A particular advantage of using graphene (assuming defect-free) is its unique two-dimensional electron gas properties at room temperature, leading to exceptionally high mobility  . A common method to grow large area graphene is by chemical vapor deposition (CVD)  . However, such method requires deposition temperature of approximately 1000˚C  . Furthermore, the produced graphene film needs to be transferred to the substrate of interest. Research is on-going to circumvent these limitations and a recent approach has been to investigate graphene growth from amorphous carbon (a-C)  . Understanding the graphitization mechanism of a-C is therefore essential to produce quality graphene.
In 2013, Barreiro et al. reported the growth of graphene from catalyst-free a-C by current-induced annealing  . They observed the structural evolution both with in-situ transmission electron microscopy and with molecular dynamics simulation  ; and a-C clusters were simulated as part of the work. An intriguing observation was the formation of long fibers stem from the clusters at elevated temperature. These clusters with fibers acted as carbon source and were able to heal up defects on a graphene sheet  . DFTB method has successfully been used to model and predict properties of carbon allotropes  -  . DFTB is based on Density Functional Theory but uses empirical approximations to reduce computational time  . Formalism of DFTB can be found in    . In this work, we performed molecular simulation using DFTB method on a-C cluster to further understand the transformation dynamics of sp3 to sp2, as well as the origin and processes by which the carbonous fibers emerged.
2. Simulation Details
An amorphous carbon cluster with diameter 1 nm was constructed using the Materials Visualizer within the Materials Studio environment  , as shown in Figure 1. The molecular dynamics simulation was performed using the DFTB+™ which includes the Van der Waals interaction. The dynamics were performed at temperatures of 300 K, 723 K (TFT annealing temperature), 1200 K, 1500 K and 2000 K. Duration of each dynamical simulation was 10 ps with 1 fs time-step to satisfy the Verlet assumption. Geometry optimization (structural relaxation) was performed prior to the dynamics simulation. Canonical Nose-Hoover isothermal constrained ensemble was used to supply heat to the system at constant volume and temperature. The molecular structures were optimized prior to any dynamical simulation.
3. Results and Discussion
The DFTB dynamics simulation on the a-C cluster revealed a progressive increase in the population of sp2 sites as the temperature was raised in a stepwise manner. Figure 2 shows the difference in the sp3 population between 300 K and 2000 K. There were only two sp3 sites remained in the cluster at the end of 2000 K annealing which is consistent with the work reported in Ref.  . In the absence of substrate effect the free-standing a-C cluster also expectedly evolved into a fullerene-like structure instead of graphitic layer  . Stone-Wales defect commonly found on graphene sheet  can also be seen in the outer surface of the annealed cluster. 5- and 7-member rings were observed on the outer surface of the cluster. Further study will focus on ways to minimize the defects during transformation and to promote formation of 6-member carbon rings.
This work provided detailed insight into the process by which a tetrahedral sp3 transformed into a planar sp2 at elevated temperature. The bond connecting the carbon atom on top of the tetrahedral and the center atom (indicated by an “X”) weakened and eventually dissociates as the bond length increases, as illustrated in Figure 3.
Highly reactive carbonous “fiber structures” responsible for graphene formation and defect healing has been reported  but no detail was given as to how they developed. This work shed some light on the process. Figure 4 shows two snapshots taken from the same dynamic trajectory at 1200 K. The cyan colored box highlights
Figure 1. The constructed a-C cluster consists of a mixture of sp3 and sp2 sites.
Figure 2. Population of sp3 sites (magenta colored) on the a-C cluster at (a) 300 K and (b) 2000 K.
Figure 3. The sp3-sp2 transformation at elevated temperatures: weakening and dissociation of the C-C bond as indicated by an “X”.
Figure 4. Two snapshots taken from the same trajectory at 1200 K. The enclosed cyan colored box shows (a) before and (b) after the sp3 bond was transformed into sp2, giving rise to a 7-member sp2 hybridized carbon ring and a “fiber structure”.
the dissociation of a sp3 site occurred at time equals 6.77 ps, leading to a simultaneous emergence of a 7-mem- ber ring and a “fiber structure”. Videos are available as ancillary files (1, 2, 3, and 4) showing the molecular dynamics simulations performed at temperatures of 300 K, 723 K, 1200 K, and 2000 K, respectively. At room temperature (300 K), despite there were activities of bond formations and dissociations, the overall ratio of sp3/sp2 remained almost stable with virtually no formation of fiber structure. This is consistent with real-life observation that amorphous carbon is stable at room temperature. As expected, there was notable increase in the amplitudes of vibrations with increasing temperature. The formation of carbonous fiber structures could be observed at 1200 K and this provides visual insights into the processes by which the fiber structures were being formed. At 2000 K, the carbon cluster transformed into a fullerene-like structure, consisting primarily of sp2 sites.
A molecular simulation based on DFTB was performed on a-C cluster to investigate the sp3-sp2 transformation process in detail including the formation of highly reactive carbonous fiber structure. Understanding of such processes will enable the direct deposition of graphene from a-C on a wide variety of substrates, which is essential for the mass-production of cost-effective graphene-based transparent electrode.
 Novoselov, K.S., et al. (2005) Two-Dimensional Atomic Crystals. Proceedings of the National Academy of Sciences of the United States of America, 102, 10451-10453.
 Barreiro, A., et al. (2013) Understanding the Catalyst-Free Transformation of Amorphous Carbon into Graphene by Current-Induced Annealing. Scientific Reports, 3, 1115.
 Irle, S., et al. (2012) Atomistic Mechanism of Carbon Nanostructure Self-Assembly as Predicted by Nonequilibrium QM/MD Simulations. Springer-European Academy of Sciences.
 Kotakoski, J., et al. (2011) Stone-Wales-Type Transformations in Carbon Nanostructures Driven by Electron Irradiation. Physical Review B, 83, Article ID: 245420.
 Terrones, H., et al. (2012) The Role of Defects and Doping in 2D Graphene Sheets and 1D Nanoribbons. Reports on Progress in Physics, 75, Article ID: 062501. http://dx.doi.org/10.1088/0034-4885/75/6/062501