Layer-Resolved Resonance Intensity of Evanescent Polariton Modes in Anisotropic Multilayers (Matlab Files)
Abstract
The attached Matlab files allow the calculation of layer-resolved mode intensity after evanescent excitation for any number of arbitrarily anisotropic or isotropic, absorbing and non-absorbing multilayer systems. The formalism builds on the electric field distributions obtained from our 4x4 transfer matrix formalism (Passler and Paarmann, JOSA B 34, 2128 (2017)) in its corrected form (Passler & Paarmann, Erratum, JOSA B 36, 3246 (2019)) with the respective Matlab (zenodo link) and Python (zenodo link) implementations. The code further extends a previous implementation of layer-resolved absorption,(Passler et al., PRB B 101, 165425 (2020), Zenodo code: matlab, python). Here, now also evanescent excitation is treated, allowing for a layer- or depth-resolved analysis of multilayer systems. We develop a empirical model based on Poynting vector analysis that is valid also for evanescent excitation.The code provided here includes the transfer matrix implementation. The respective publication was recently submitted, with its preprint available on the arXiv. We also note a small improvement of the TMM formalism over the previous versions(TMM, Layer resolved absorption). We noticed that in the case of birefringent layers, in some very special cases the mode ordering and eigenmode calculations according to Xu et al Phys. Rev. B 61, 1740–1743 (2000) used in our previous code produced faulty results, likely linked to the ambuigity of mode sorting in the case of birefringence. Therefore, we now instead directly use the Berreman eigenmodes for birefringent layers, which removes the necessity of mode sorting and thereby also circumvents the rare issues encountered with the previous implementations. All calculations preformed with this new implementation did not show any problems. Apart from the core code file (passler_layer_resolved_absorption_core.m) and an auxilliary function for generation of the dielectric tensors for all layers (passler_epsTarray_generator.m), we also provide two example scripts which demonstrate the capabilities of the formalism. These sample scripts also generate the figures shown in the associated publication (arXiv). For any questions please contact Alex Paarmann,
- Publication:
-
Zenodo
- Pub Date:
- DOI:
- 10.5281/zenodo.7034721
- Bibcode:
- 2022zndo...7034721P