Block elimination method for optical near-fields

MATLAB implementation of the Block eleminination method for calculating the optical near-field properties of metal-dielectric composites.

This is a MATLAB implementation of the block elimination method from: [GENOV, Dentcho A.; SARYCHEV, Andrey K.; SHALAEV, Vladimir M. Plasmon localization and local field distribution in metal-dielectric films. Physical Review E, 2003, 67.5: 056611.].

The method is a finite element method based on RLC-circuit theory and can solve the optical near-field distributions of complex metal-dielectric composites in a very fast matrix inversion method.

Input of a structure is based on loading a .png image file where white pixels = metal, black pixels = dielectric (an example is included). The solver works in the x- and y-dimensions (z-axis is assumed infinite). All calculations are done with periodic boundary conditions. The solver will simultanously calculate the fields for both an incident x- and y-polarized optical field (field is injected along the z-axis).

The method is further described in: [FRYDENDAHL, Christian. Thesis: Experimental exploration of gold semi-continuous films in the near-and far-field. arXiv preprint arXiv:1906.01375, 2019.]