  We present the molecular hyperdynamics algorithm and its implementation to
the nonorthogonal tight-binding model NTBM and the corresponding software. Due
to its multiscale structure, the proposed approach provides the long time scale
simulations (more than 1 s), unavailable for conventional molecular dynamics.
No preliminary information about the system potential landscape is needed for
the use of this technique. The optimal interatomic potential modification is
automatically derived from the previous simulation steps. The average time
between adjusted potential energy fluctuations provides an accurate evaluation
of physical time during the hyperdynamics simulation. The main application of
the presented hyperdynamics method is the study of thermal-induced defects
arising in the middle-sized or relatively large atomic systems at low
temperatures. To validate the presented method, we apply it to the C$_{60}$
cage and its derivative C$_{60}$NH$_{2}$. Hyperdynamics leads to the same
results as a conventional molecular dynamics, but the former possesses much
higher performance and accuracy due to the wider temperature region. The
coefficient of acceleration achieves 10$^{7}$ and more.
