The tight-binding method was employed in the calculations
to describe interactions between carbon atoms.
With this method, the molecular dynamics technique was applied
to calculate the positions and the velocities of the
atoms as function of the time. In the MD calculation,
Newton's equations of motion were solved
using the predictor-corrector algorithm. The MD time step was
0.5
10
s.
Periodic boundaries conditions
were applied to the problem.
In order to calculate the ratio of
to
bonded atoms,
as well as the radial and angle distribution functions
of the structures of amorphous carbon generated a FORTRAN
program was written.