ReaxFF Molecular Dynamics Simulation for the Graphitization of Amorphous Carbon: A Parametric Study

A parametric study of ReaxFF for molecular dynamics simulation of graphitization of amorphous carbon was conducted. The responses to different initial amorphous carbon configurations, simulation time steps, simulated temperatures, and ReaxFF parameter sets were investigated. The results showed that a time step shorter than 0.2 fs is sufficient for the ReaxFF simulation of carbon using both Chenoweth 2008 and Srinivasan 2015 parameter sets. The amorphous carbon networks produced using both parameter sets at 300 K are similar to each other, with the first peak positions of pair distribution function curves located between the graphite sp2 bond peak position and the diamond sp3 bond peak position. In the graphitization process, the graphene fragment size increases and the orientation of graphene layers transforms to be parallel with each other with the increase of temperature and annealing time. This parallel graphene structure is close to the crystalline graphite. Associated with this graphitization is the ...