Today, the renaissance of black phosphorus largely depends on the mechanical exfoliation method, which is accessible to produce few-layer forms from the bulk counterpart. However, the deep understanding of the exfoliation mechanism is missing. To this end, we resolve this issue by simulating the sliding processes of bilayer phosphorene based on first-principles calculations. It is found that the interlayer Coulomb interactions dictate the optimal sliding pathway, leading to the minimal energy barrier as low as , which gives a comparable surface energy of in experiment. This means that black phosphorus can be exfoliated by the sliding approach. In addition, considerable bandgap modulations along these sliding pathways are obtained. The study like ours builds up a fundamental understanding of how black phosphorus is exfoliated to few-layer forms, providing a good guide to experimental research.