Toward a realistic description of multilayer black phosphorus: from GW approximation to large-scale tight-binding simulations

We provide a tight-binding model parametrization for black phosphorus (BP) with an arbitrary number of layers. The model is derived from partially self-consistent $GW_0$ approach, where the screened Coulomb interaction $W_0$ is calculated within the random phase approximation on the basis of density functional theory. We thoroughly validate the model by performing a series of benchmark calculations, and determine the limits of its applicability. The application of the model to the calculations of electronic and optical properties of multilayer BP demonstrates good quantitative agreement with \emph{ab initio} results in a wide energy range. We also show that the proposed model can be easily extended for the case of external fields, yielding the results consistent with those obtained from first principles. The model is expected to be suitable for a variety of realistic problems related to the electronic properties of multilayer BP including different kinds of disorder, external fields, and many-body effects.