Why twisting angles are diverse in graphene Moiré patterns?

The interlayer energy of the twisting bilayer graphene is investigated by the molecular mechanics method using both the registry-dependent potential and the Lennard-Jones potential. Both potentials show that the interlayer energy is independent of the twisting angle θ, except in the two boundary regions θ ≈ 0° or 60^°, where the interlayer energy is proportional to the square of the twisting arc length. The calculation results are successfully interpreted by a single atom model. An important information from our findings is that, from the energy point of view, there is no preference for the twisting angle in the experimental bilayer graphene samples, which actually explains the diverse twisting angles in the experiment.