Applied Surface Science, Vol.286, 287-297, 2013
The effect of empirical potential functions on modeling of amorphous carbon using molecular dynamics method
Empirical potentials have a strong effect on the hybridization and structure of amorphous carbon and are of great importance in molecular dynamics (MD) simulations. In this work, amorphous carbon at densities ranging from 2.0 to 3.2 g/cm(3) was modeled by a liquid quenching method using Tersoff, 2nd REBO, and ReaxFF empirical potentials. The hybridization, structure and radial distribution function G(r) of carbon atoms were analyzed as a function of the three potentials mentioned above. The ReaxFF potential is capable to model the change of the structure of amorphous carbon and MD results are in a good agreement with experimental results and density function theory (DFT) at low density of 2.6 g/cm(3) and below. The 2nd REBO potential can be used when amorphous carbon has a very low density of 2.4 g/cm(3) and below. Considering the computational efficiency, the Tersoff potential is recommended to model amorphous carbon at a high density of 2.6 g/cm(3) and above. In addition, the influence of the quenching time on the hybridization content obtained with the three potentials is discussed. (C) 2013 Elsevier B.V. All rights reserved.