All reactive molecular dynamics simulations has carry out performed by non-equilibrium molecular dynamics method [1]}, [2]}, [3]}. The numerical atomic positions os carbon atoms are calculated by large-scale atomic/molecular massively parallel simulator (LAMMPS) code [4]}. The interatomic potential used in all calculations are a modern reactive force field ReaxFF with set parameter described by references [5]}, [6]}. The reactive force field ReaxFF is parameterized using available experimental results and first-principles calculations [5]}, [6]}, [9]}, [10]}, [11]}, [12]}, [13]}. Approximately to empirical non reactive force fields, ReaxFF is divide by partial energy contributions, as in the follow equations  [5]}:
\(E_{system}&=&E_{bond}+E_{over}+E_{under}+E_{val}\nonumber \\&&+E_{pen}+E_{tor}+E_{conj}+E_{vdW}\nonumber \\&&+E_{co} \quad ,\) 
