mie_ml: mie_ml

View source: R/mie_layer.R

mie_mlR Documentation

mie_ml

Description

Far-field cross-sections

Usage

mie_ml(
  wavelength,
  epsilon,
  radii,
  n_max = 10,
  efficiency = FALSE,
  mode = c("EM", "Magnetic", "Electric"),
  order = Inf
)

Arguments

wavelength

real vector

epsilon

list of dielectric functions, from inner to outer medium

radii

concentric radii of each interface, from smaller to larger

n_max

truncation order

efficiency

logical, scale by geometrical cross-sections

mode

type of mode

order

order of multipoles

Details

Multilayered sphere illuminated by a plane wave

Value

data.frame

Author(s)

Baptiste Auguie

See Also

Other user: mie-package, mie_approximation(), mie_bh(), mie()

Examples

library(dielectric)
library(mie)
gold <- epsAg(seq(300, 800))
a <- 30
b <- 34
c <- 35

bare <- mie(gold$wavelength, gold$epsilon, radius=a, medium=1.33, efficiency=FALSE)

leps <- list(gold$epsilon, 1.33^2, 1.5^2, 1.33^2)
# leps <- list(1.5^2, gold$epsilon,  1.33^2)
la <- list(a,b,c)
coated <- mie_ml(gold$wavelength, leps, radii=la, efficiency=FALSE)

matplot(bare$wavelength, bare[, -1], type="l", lty=1,
        xlab=expression(lambda/mu*m), ylab=expression(sigma/mu*m^2))
matlines(coated$wavelength, coated[, -1],  type="l",  lty=2)

legend("topright", c(names(bare)[-1], "coated"), col=1:3, lty=c(1,1,1,2))

baptiste/mie documentation built on July 4, 2023, 1:06 p.m.