Size effects and strain localization in atomic-scale cleavage modeling

In this work, we study the adhesion and decohesion of Cu(1 0 0) surfaces using density functional theory (DFT) calculations. An upper stress to surface decohesion is obtained via the universal binding energy relation (UBER), but the model is limited to rigid separation of bulk-terminated surfaces. When structural relaxations are included, an unphysical size effect arises if decohesion is considered to occur as soon as the strain energy equals the energy of the newly formed surfaces. We employ the nudged elastic band (NEB) method to show that this size effect is opposed by a size-dependency of the energy barriers involved in the transition. Further, we find that the transition occurs via a localization of bond strain in the vicinity of the cleavage plane, which resembles the strain localization at the tip of a sharp crack that is predicted by linear elastic fracture mechanics.