“…In addition to reflection and transmission coefficients, we also calculate the full electric field distributions throughout the heterostructure, which is particularly useful for analysis of polariton modes and their associated local field enhancements. To demonstrate the capabilities of our algorithm, we present simulation results for several model systems where SPhPs can be excited using prism coupling in the Otto geometry 20 : 6H-SiC, GaN/SiC, SiC/GaN/SiC, and α-quartz, covering a number of phenomena such as critical coupling to SPhPs 21,22 , index-sensing 35 , modesplitting 36 , wave-guiding 37 , and polaritons in hyperbolic media 29,32 .…”