Ab initio study of electromagnetic modes in two-dimensional semiconductors: Application to doped phosphorene