  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.
