Multiscale molecular modeling of surface anisotropy and electronic properties of phosphorene