How to Combine Binary Collision Approximation and Multi-Body Potential for Molecular Dynamics

Our group has been developing a hybrid simulation of the molecular dynamics (MD) and the binary collision approximation (BCA) simulation. One of the main problems of this hybridization model is that the multi-body potential suddenly appears at the moment when the simulation method switches from the BCA to the MD. This instantaneously emerged multi-body potential causes the acceleration or deceleration of atoms of the system. To solve this problem, the kinetic energy of atoms should be corrected to conserve the total energy in the system. This paper gives the solution. The hybrid simulation for hydrogen atom injection into a graphite material is executed in order to demonstrate the solution.