Nanoscale continuum calculation of basal dislocation core structures in graphite

In this work, we calculate the core structures of basal dislocations in graphite in a nanoscale continuum framework. The model consists of a stack of buffered Kirchhoff plates where the plates represent the covalent interactions within individual graphene sheets and the buffer layers represent the secondary interactions between them. In the mid-plane of the buffer layers, cohesive surfaces are introduced to account for the nonlinear deformations due to basal dislocations. The cohesive surface separation is governed by using an empirical 4-8 Lennard–Jones potential. Meanwhile, their relative shear sliding is governed by using a newly proposed empirical periodic stacking-fault potential. With these potentials, the core structures of full dislocations and partials are calculated and examined. It is shown that the full dislocations automatically split into partials that repel each other. The core sizes of individual partials, measured between peak stresses, are about 5?nm wide for the edge component and slightly narrower for the screw component. Since these sizes are about 10 times the lattice constant, they lend credence to our continuum model of basal dislocation cores in graphite. It is also shown that when the dislocations are densely packed on the same glide plane, i.e. in a pile-up, with spacing one to two times the core size, the split partials retain their individual identity with well-defined and well-separated stress peaks. Meanwhile, the membrane normal stresses in the graphene sheets rise considerably at the pile-up tips which, in turn, may provoke further deformation and damage modes such as kinking and delamination.