“…Here we have assumed an exponential form for the contribution to the energy as a function of depth. Hung et al [2] have shown that the electron density in slabs of increasing widths rapidly approaches its bulk value .The total energy E tot can be calculated as:Where N A is the number of atoms in-plane, N D is the number of layers (although the sum goes to infinity here, in practice N D is very large but not infinite).Calling N = N A .N D , the total number of atoms, we obtain:Hung et al [3] have done a first-principles density functional theory calculation where they obtained the energy as a function of separation of two surface-terminated bulks. Call that energy curve E(γd 0 ), where d 0 = ca 0 is the equilibrium interplanar distance, a 0 is the equilibrium lattice parameter, and c is a geometric factor that depends on the orientation and the crystal structure.…”