Dipole total decay rate near a multilayer interface
dipole_direct(d = 1, wavelength, epsilon = list(incident = 1^2),
thickness = c(0, 0), Nquadrature1 = 50, Nquadrature2 = 200,
Nquadrature3 = 50, qcut = NULL, qmax = Inf, show.messages = TRUE)
Arguments
- d
- distance in nm
- wavelength
- wavelength in nm
- epsilon
- list of dielectric functions
- thickness
- list of layer thicknesses
- Nquadrature1
- quadrature points in radiative region
- Nquadrature2
- quadrature points in SPPs region
- Nquadrature3
- quadrature points in dipole image region
- qcut
- transition between regions 2 and 3
- qmax
- maximum q of region 3
- show.messages
- logical, display integration info
Details
direct application of the textbook formula using integrand_mtot; performs poorly compared to the transformed version in dipole