“…However, it is exactly in these materials and atomic-scale heterostructures thereof [23][24][25][26] , that SPhPs have very recently been demonstrated to enable many novel phenomena such as hyperbolic superlensing 27,28 and negative refraction 29,30 . For instance, hexagonal boron nitride as one of the key components of van der Waals heterostructures 31 , displays many interesting nanophotonic properties due to its naturally hyperbolic character 23,28,30,32,33 . Therefore, a most general, robust, and easily implementable numera) Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6,14195 Berlin, Germany; Electronic mail: passler@fhi-berlin.mpg.de b) Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6,14195 Berlin, Germany ical formalism to analyze the optical response in arbitrarily anisotropic multilayer heterostructures is highly desirable.…”