The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus

Sensitive dependence of the electronic structure on the number of layers in few-layer phosphorene raises a question about the true nature of the interlayer interaction in so-called "van der Waals (vdW) solids". We performed quantum Monte Carlo calculations and found that the interlayer interaction in bulk black phosphorus and related few-layer phosphorene is associated with a significant charge redistribution that is incompatible with purely dispersive forces and not captured by density functional theory calculations with different vdW corrected functionals. These findings confirm the necessity of more sophisticated treatment of nonlocal electron correlation in total energy calculations.