To examine the thermal and mechanical stability, the fully atomistic classical molecular dynamics simulations were performed using the tersoff potential as implemented in LAMMPS . The fullerene structures have dimensions of 80.1*73.4 Å 2 with periodic boundary conditions and are composed of 4800 atoms. The length of the simulation box along the z-direction is 130 Å. The equations of motion are integrated us-ing the velocity-Verlet algorithm, with a time-step of 0.5 fs. Before heating and stretching, the structures were equilibrated using an NPT ensemble at a constant temperature of 300 K using a Nose-Hoover thermostat for 200 ps. The structures were continuously stretched up to their complete structural failure (fracture) by applying a maximum strain of 25%. The tensile stretching simulations were performed here by only deforming in the y directions. Their thermal stability was investigated by heating up the structures from 300 K to 10,000 K by linearly increasing the temperature during 1 ns in an NVT ensemble. The molecular dynamics snapshots and trajectories were visualized using the visualization software VMD .