Efficient structural relaxation of polycrystalline graphene models

Federico D'Ambrosio,Joris Barkema,Gerard T Barkema

DOI: https://doi.org/10.3390/nano11051242

2021-05-08

Abstract:Large samples of experimentally produced graphene are polycrystalline. For the study of this material, it helps to have realistic computer samples that are also polycrystalline. A common approach to produce such samples in computer simulations is based on the method of Wooten, Winer, and Weaire, originally introduced for the simulation of amorphous silicon. We introduce an early rejection variation of their method, applied to graphene, which exploits the local nature of the structural changes to achieve a significant speed-up in the relaxation of the material, without compromising the dynamics. We test it on a 3,200 atoms sample, obtaining a speedup between one and two orders of magnitude. We also introduce a further variation called early decision specifically for relaxing large samples even faster and we test it on two samples of 10,024 and 20,000 atoms, obtaining a further speed-up of an order of magnitude. Furthermore, we provide a graphical manipulation tool to remove unwanted artifacts in a sample, such as bond crossings.

Materials Science,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to efficiently perform structural relaxation of polycrystalline graphene models. Specifically, large - sized graphene samples produced in experiments are usually polycrystalline and contain many internal and external lattice defects. To study this material, it is necessary to generate computer - simulated samples that are also polycrystalline. Although existing methods such as the Wooten, Winer, and Weaire (WWW) algorithm can generate such samples, the computational cost is very high during the relaxation process, especially when dealing with large samples.

The paper proposes a method based on local relaxation. By using two techniques, early rejection and early decision, to accelerate the structural relaxation process, thereby significantly reducing the computational time required for each accepted bond - exchange operation. These techniques take advantage of the characteristics of local structural changes, avoid global relaxation of the entire sample, and thus achieve an acceleration effect of 1 to 2 orders of magnitude.

### Main problem summary:
1. **High computational cost**: When dealing with large - sized polycrystalline graphene samples, existing methods require frequent global relaxation, resulting in excessive computational cost.
2. **Low efficiency**: As the sample energy gradually decreases, the proportion of accepted bond exchanges becomes extremely low, and most of the computational time is wasted on rejected bond - exchange attempts.

### Solutions:
- **Early Rejection**: After bond exchange, only relax the local area, and decide whether to accept the bond exchange according to the energy estimate after local relaxation. If the energy after local relaxation exceeds the preset threshold, immediately reject the bond exchange to avoid global relaxation.
- **Early Decision**: Further simplify it to directly decide whether to accept the bond exchange after local relaxation without performing global relaxation, which is suitable for very large samples.

Through these methods, the paper shows how to significantly improve the computational efficiency of polycrystalline graphene structure relaxation while ensuring the accuracy of the results.