The mechanical behavior of discrete atoms has been studied with molecular dynamics whose computational time is proportional to the square of the number of atoms, O(N2). Recently, a faster algorithm, Atomic-scale Finite Element Method (AFEM) with computational time proportional to the number of atoms, O(N), had been developed. The main idea of AFEM, compared with conventional finite element method is to replace nodes with atoms and elements with electric forces between atoms. When interpreting a non-linear system, it is necessary to use an iteration scheme.
A simulation of molecular dynamics based on the Verlet's method was conducted in order to validate AFEM in one dimension. The speed of AFEM was investigated in one and two dimensional atomic systems. The results showed that the computational time of AFEM is approximately proportional to the number of atoms, and the absolute computation time appears to be small. The frameworks of AFEM not only for multi-body potential but also pair potential are presented. Finally, AFEM was applied to analyze and interpret the mechanical behavior of a carbon nanotube and a quartz. The buckling behavior of carbon nanotube showed a good agreement with the results illustrated in the original literature.
