Mie scattering

редактировать
Mie scattering, artistic view Monostatic radar cross section (RCS) of a perfectly conducting metal sphere as a function of frequency (calculated by Mie theory). In the low-frequency Rayleigh scattering limit, where the circumference is less than the wavelength, the normalized RCS is σ/(πR) ~ 9(kR). In the high-frequency optical limit σ/(πR) ~ 1.

The Mie solutionto Maxwell's equations (also known as the Lorenz–Mie solution, the Lorenz–Mie–Debye solutionor Mie scattering) describes the scattering of an electromagnetic plane wave by a homogeneous sphere. The solution takes the form of an infinite series of spherical multipole partial waves. It is named after Gustav Mie.

The term Mie solution is also used for solutions of Maxwell's equations for scattering by stratified spheres or by infinite cylinders, or other geometries where one can write separate equations for the radial and angular dependence of solutions. The term Mie theory is sometimes used for this collection of solutions and methods; it does not refer to an independent physical theory or law. More broadly, "Mie scattering" suggests situations where the size of the scattering particles is comparable to the wavelength of the light, rather than much smaller or much larger.

Mie scattering(sometimes referred to as a non-molecular scatteringor aerosol particle scattering) takes place in the lower 4,500 metres (15,000 feet) of the atmosphere, where there may be many essentially spherical particles present with diameters approximately equal to the size of the wavelength of the incident ray. Mie scattering theory has no upper size limitation, and converges to the limit of geometric optics for large particles.

Contents

  • 1 Introduction
  • 2 Approximations
    • 2.1 Rayleigh approximation (scattering)
    • 2.2 Rayleigh–Gans approximation
    • 2.3 Anomalous diffraction approximation of van de Hulst
  • 3 Mathematics
    • 3.1 Scattering and extinction cross-sections
    • 3.2 Application to subwavelength particles
    • 3.3 Other directions of the incident plane wave
  • 4 Kerker effect
  • 5 Dyadic Green's function of a sphere
  • 6 Computational codes
  • 7 Applications
    • 7.1 Atmospheric science
    • 7.2 Cancer detection and screening
    • 7.3 Clinical laboratory analysis
    • 7.4 Magnetic particles
    • 7.5 Metamaterial
    • 7.6 Particle sizing
    • 7.7 Parasitology
  • 8 Extensions
  • 9 See also
  • 10 References
  • 11 Further reading
  • 12 External links

Introduction

Angular part of magnetic and electric vector spherical harmonics. Red and green arrows show the direction of the field. Generating scalar functions are also presented, only the first three orders are shown (dipoles, quadrupoles, octupoles).

A modern formulation of the Mie solution to the scattering problem on a sphere can be found in many books, e.g., J. A. Stratton 's Electromagnetic Theory. In this formulation, the incident plane wave, as well as the scattering field, is expanded into radiating spherical vector spherical harmonics. The internal field is expanded into regular vector spherical harmonics. By enforcing the boundary condition on the spherical surface, the expansion coefficients of the scattered field can be computed.

For particles much larger or much smaller than the wavelength of the scattered light there are simple and accurate approximations that suffice to describe the behavior of the system. But for objects whose size is similar to the wavelength, e.g., water droplets in the atmosphere, latex particles in paint, droplets in emulsions, including milk, and biological cells and cellular components, a more detailed approach is necessary.

The Mie solution is named after its developer, German physicist Gustav Mie. Danish physicist Ludvig Lorenz and others independently developed the theory of electromagnetic plane wave scattering by a dielectric sphere.

The formalism allows the calculation of the electric and magnetic fields inside and outside a spherical object and is generally used to calculate either how much light is scattered (the total optical cross section ), or where it goes (the form factor). The notable features of these results are the Mie resonances, sizes that scatter particularly strongly or weakly. This is in contrast to Rayleigh scattering for small particles and Rayleigh–Gans–Debye scattering (after Lord Rayleigh, Richard Gans and Peter Debye ) for large particles. The existence of resonances and other features of Mie scattering makes it a particularly useful formalism when using scattered light to measure particle size.

Approximations

Rayleigh approximation (scattering)

The change of sky colour at sunset (red nearest the sun, blue furthest away) is caused by Rayleigh scattering by atmospheric gas particles, which are much smaller than the wavelengths of visible light. The grey/white colour of the clouds is caused by Mie scattering by water droplets, which are of a comparable size to the wavelengths of visible light.

Rayleigh scattering describes the elastic scattering of light by spheres that are much smaller than the wavelength of light. The intensity I of the scattered radiation is given by

I = I 0 ( 1 + cos 2 ⁡ θ 2 R 2) ( 2 π λ) 4 ( n 2 − 1 n 2 + 2) 2 ( d 2) 6, {\displaystyle I=I_{0}\left({\frac {1+\cos ^{2}\theta }{2R^{2}}}\right)\left({\frac {2\pi }{\lambda }}\right)^{4}\left({\frac {n^{2}-1}{n^{2}+2}}\right)^{2}\left({\frac {d}{2}}\right)^{6},}{\displaystyle I=I_{0}\left({\frac {1+\cos ^{2}\theta }{2R^{2}}}\right)\left({\frac {2\pi }{\lambda }}\right)^{4}\left({\frac {n^{2}-1}{n^{2}+2}}\right)^{2}\left({\frac {d}{2}}\right)^{6},}

where I0is the light intensity before the interaction with the particle, R is the distance between the particle and the observer, θ is the scattering angle, n is the refractive index of the particle, and d is the diameter of the particle.

It can be seen from the above equation that Rayleigh scattering is strongly dependent upon the size of the particle and the wavelengths. The intensity of the Rayleigh scattered radiation increases rapidly as the ratio of particle size to wavelength increases. Furthermore, the intensity of Rayleigh scattered radiation is identical in the forward and reverse directions.

The Rayleigh scattering model breaks down when the particle size becomes larger than around 10% of the wavelength of the incident radiation. In the case of particles with dimensions greater than this, Mie's scattering model can be used to find the intensity of the scattered radiation. The intensity of Mie scattered radiation is given by the summation of an infinite series of terms rather than by a simple mathematical expression. It can be shown, however, that scattering in this range of particle sizes differs from Rayleigh scattering in several respects: it is roughly independent of wavelength and it is larger in the forward direction than in the reverse direction. The greater the particle size, the more of the light is scattered in the forward direction.

The blue colour of the sky results from Rayleigh scattering, as the size of the gas particles in the atmosphere is much smaller than the wavelength of visible light. Rayleigh scattering is much greater for blue light than for other colours due to its shorter wavelength. As sunlight passes through the atmosphere, its blue component is Rayleigh scattered strongly by atmospheric gases but the longer wavelength (e.g. red/yellow) components are not. The sunlight arriving directly from the Sun therefore appears to be slightly yellow, while the light scattered through rest of the sky appears blue. During sunrises and sunsets, the effect of Rayleigh scattering on the spectrum of the transmitted light is much greater due to the greater distance the light rays have to travel through the high-density air near the Earth's surface.

In contrast, the water droplets that make up clouds are of a comparable size to the wavelengths in visible light, and the scattering is described by Mie's model rather than that of Rayleigh. Here, all wavelengths of visible light are scattered approximately identically, and the clouds therefore appear to be white or grey.

Rayleigh–Gans approximation

The Rayleigh–Gans approximation is an approximate solution to light scattering when the relative refractive index of the particle is close to that of the environment, and its size is much smaller in comparison to the wavelength of light divided by |n − 1|, where n is the refractive index :

| n − 1 | ≪ 1 {\displaystyle |n-1|\ll 1}{\displaystyle |n-1|\ll 1}
k d | n − 1 | ≪ 1 {\displaystyle kd|n-1|\ll 1}{\displaystyle kd|n-1|\ll 1}

where k {\textstyle k}{\textstyle k}is the wavevector of the light (k = 2 π λ {\textstyle k={\frac {2\pi }{\lambda }}}{\textstyle k={\frac {2\pi }{\lambda }}}), and d {\displaystyle d}drefers to the linear dimension of the particle. The former condition is often referred as the "optically soft" and the approximation holds for particles of arbitrary shape.

Anomalous diffraction approximation of van de Hulst

The anomalous diffraction approximation is valid for large (compared to wavelength) and optically soft spheres; soft in the context of optics implies that the refractive index of the particle (m) differs only slightly from the refractive index of the environment, and the particle subjects the wave to only a small phase shift. The extinction efficiency in this approximation is given by

Q = 2 − 4 p sin ⁡ p + 4 p 2 ( 1 − cos ⁡ p), {\displaystyle Q=2-{\frac {4}{p}}\sin p+{\frac {4}{p^{2}}}(1-\cos p),}{\displaystyle Q=2-{\frac {4}{p}}\sin p+{\frac {4}{p^{2}}}(1-\cos p),}

where Q is the efficiency factor of scattering, which is defined as the ratio of the scattering cross-section and geometrical cross-section πa.

The term p = 4πa(n − 1)/λ has as its physical meaning the phase delay of the wave passing through the centre of the sphere, where a is the sphere radius, n is the ratio of refractive indices inside and outside of the sphere, and λ the wavelength of the light.

This set of equations was first described by van de Hulst in (1957).

Mathematics

Scattering of the plane wave, incidence direction is parallel to the z-axis, polarization is parallel to the x-axis, nanoparticle's radius is a

The scattering by a spherical nanoparticle is solved exactly regardless of the particle size. We consider scattering by a plane wave propagating along the z-axis polarized along the x-axis. Dielectric and magnetic permeabilities of a particle are ε 1 {\displaystyle \varepsilon _{1}}\varepsilon _{1}and μ 1 {\displaystyle \mu _{1}}\mu _{1}, and ε {\displaystyle \varepsilon }\varepsilon and μ {\displaystyle \mu }\mu for the environment.

In order to solve the scattering problem, we write first the solutions of the vector Helmholtz equation in spherical coordinates, since the fields inside and outside the particles must satisfy it. Helmholtz equation:

∇ 2 E + k 2 E = 0, ∇ 2 H + k 2 H = 0 {\displaystyle \nabla ^{2}\mathbf {E} +{k}^{2}\mathbf {E} =0,\ \ \ \ \nabla ^{2}\mathbf {H} +{k}^{2}\mathbf {H} =0}{\displaystyle \nabla ^{2}\mathbf {E} +{k}^{2}\mathbf {E} =0,\ \ \ \ \nabla ^{2}\mathbf {H} +{k}^{2}\mathbf {H} =0}

in addition to the Helmholtz equation, the fields must satisfy the conditions ∇ ⋅ E = ∇ ⋅ H = 0 {\displaystyle \nabla \cdot \mathbf {E} =\nabla \cdot \mathbf {H} =0}{\displaystyle \nabla \cdot \mathbf {E} =\nabla \cdot \mathbf {H} =0}and ∇ × E = i ω μ H {\displaystyle \nabla \times \mathbf {E} =i\omega \mu \mathbf {H} }{\displaystyle \nabla \times \mathbf {E} =i\omega \mu \mathbf {H} }, ∇ × H = − i ω ε E {\displaystyle \nabla \times \mathbf {H} =-i\omega \varepsilon \mathbf {E} }{\displaystyle \nabla \times \mathbf {H} =-i\omega \varepsilon \mathbf {E} }. Vector spherical harmonics possess all the necessary properties, introduced as follows:

M o e m n = ∇ × ( r ψ o e m n) {\displaystyle \mathbf {M} _{^{e}_{o}mn}=\nabla \times \left(\mathbf {r} \psi _{^{e}_{o}mn}\right)}{\displaystyle \mathbf {M} _{^{e}_{o}mn}=\nabla \times \left(\mathbf {r} \psi _{^{e}_{o}mn}\right)}— magnetic harmonics (TE)
N o e m n = ∇ × M o e m n k {\displaystyle \mathbf {N} _{^{e}_{o}mn}={\frac {\nabla \times \mathbf {M} _{^{e}_{o}mn}}{k}}}{\displaystyle \mathbf {N} _{^{e}_{o}mn}={\frac {\nabla \times \mathbf {M} _{^{e}_{o}mn}}{k}}}— electric harmonics (TM)

where

ψ e m n = cos ⁡ m φ P n m ( cos ⁡ ϑ) z n ( k r) {\displaystyle {\psi _{emn}=\cos m\varphi P_{n}^{m}(\cos \vartheta)z_{n}({k}r)}}{\displaystyle {\psi _{emn}=\cos m\varphi P_{n}^{m}(\cos \vartheta)z_{n}({k}r)}}
ψ o m n = sin ⁡ m φ P n m ( cos ⁡ ϑ) z n ( k r) {\displaystyle {\psi _{omn}=\sin m\varphi P_{n}^{m}(\cos \vartheta)z_{n}({k}r)}}{\displaystyle {\psi _{omn}=\sin m\varphi P_{n}^{m}(\cos \vartheta)z_{n}({k}r)}}

and P n m ( cos ⁡ θ) {\displaystyle P_{n}^{m}(\cos \theta)}{\displaystyle P_{n}^{m}(\cos \theta)}Associated Legendre polynomials, and z n ( k r) {\displaystyle z_{n}({k}r)}{\displaystyle z_{n}({k}r)}— any of the spherical bessel functions.

Next, we expand the incident plane wave in vector spherical harmonics:

E i n c = E 0 e i k r cos ⁡ θ e x = E 0 ∑ n = 1 ∞ i n 2 n + 1 n ( n + 1) ( M o 1 n ( 1) ( k, r) − i N e 1 n ( 1) ( k, r)) {\displaystyle \mathbf {E} _{inc}=E_{0}e^{ikr\cos \theta }\mathbf {e} _{x}=E_{0}\sum _{n=1}^{\infty }i^{n}{\frac {2n+1}{n(n+1)}}\left(\mathbf {M} _{o1n}^{(1)}(k,\mathbf {r})-i\mathbf {N} _{e1n}^{(1)}(k,\mathbf {r})\right)}{\displaystyle \mathbf {E} _{inc}=E_{0}e^{ikr\cos \theta }\mathbf {e} _{x}=E_{0}\sum _{n=1}^{\infty }i^{n}{\frac {2n+1}{n(n+1)}}\left(\mathbf {M} _{o1n}^{(1)}(k,\mathbf {r})-i\mathbf {N} _{e1n}^{(1)}(k,\mathbf {r})\right)}
H i n c = − k ω μ E 0 ∑ n = 1 ∞ i n 2 n + 1 n ( n + 1) ( M e 1 n ( 1) ( k, r) + i N o 1 n ( 1) ( k, r)) {\displaystyle \mathbf {H} _{inc}={\frac {-k}{\omega \mu }}E_{0}\sum _{n=1}^{\infty }i^{n}{\frac {2n+1}{n(n+1)}}\left(\mathbf {M} _{e1n}^{(1)}(k,\mathbf {r})+i\mathbf {N} _{o1n}^{(1)}(k,\mathbf {r})\right)}{\displaystyle \mathbf {H} _{inc}={\frac {-k}{\omega \mu }}E_{0}\sum _{n=1}^{\infty }i^{n}{\frac {2n+1}{n(n+1)}}\left(\mathbf {M} _{e1n}^{(1)}(k,\mathbf {r})+i\mathbf {N} _{o1n}^{(1)}(k,\mathbf {r})\right)}

here the superscript ( 1) {\displaystyle (1)}(1)means that in the radial part of the functions ψ o e m n {\displaystyle \psi _{^{e}_{o}mn}}{\displaystyle \psi _{^{e}_{o}mn}}are spherical Bessel functions. The expansion coefficients are obtained by taking integrals of the form

∫ 0 2 π ∫ 0 π E i n c ⋅ M o e m n ( 1) sin ⁡ θ d θ d φ ∫ 0 2 π ∫ 0 π | M o e m n ( 1) | 2 sin ⁡ θ d θ d φ {\displaystyle {\frac {\int _{0}^{2\pi }\int _{0}^{\pi }\mathbf {E} _{inc}\cdot \mathbf {M} _{^{e}_{o}mn}^{(1)}\sin \theta d\theta d\varphi }{\int _{0}^{2\pi }\int _{0}^{\pi }|\mathbf {M} _{^{e}_{o}mn}^{(1)}|^{2}\sin \theta d\theta d\varphi }}}{\displaystyle {\frac {\int _{0}^{2\pi }\int _{0}^{\pi }\mathbf {E} _{inc}\cdot \mathbf {M} _{^{e}_{o}mn}^{(1)}\sin \theta d\theta d\varphi }{\int _{0}^{2\pi }\int _{0}^{\pi }|\mathbf {M} _{^{e}_{o}mn}^{(1)}|^{2}\sin \theta d\theta d\varphi }}}

in this case, all coefficients at m ≠ 1 {\displaystyle m\neq 1}{\displaystyle m\neq 1}are zero, since the integral over the angle φ {\displaystyle \varphi }\varphi in the numerator is zero.

. Then the following conditions are imposed:

1) Interface conditions on the boundary between the sphere and the environment (which allow us to relate the expansion coefficients of the incident, internal, and scattered fields)

2) The condition that the solution is bounded at the origin (therefore, in the radial part of the generating functions ψ o e m n {\displaystyle \psi _{^{e}_{o}mn}}{\displaystyle \psi _{^{e}_{o}mn}}, Bessel spherical functions are selected for the internal field),

3) For a scattered field, the asymptotics at infinity corresponds to a diverging spherical wave (in connection with this, for the scattered field in the radial part of the generating functions ψ o e m n {\displaystyle \psi _{^{e}_{o}mn}}{\displaystyle \psi _{^{e}_{o}mn}}spherical Hankel functions of the first kind are chosen).

Scattered fields are written in terms of a vector harmonic expansion as

E s = ∑ n = 1 ∞ E n ( i a n N e 1 n ( 3) ( k, r) − b n M o 1 n ( 3) ( k, r)) {\displaystyle \mathbf {E} _{s}=\sum _{n=1}^{\infty }E_{n}\left(ia_{n}\mathbf {N} _{e1n}^{(3)}(k,\mathbf {r})-b_{n}\mathbf {M} _{o1n}^{(3)}(k,\mathbf {r})\right)}{\displaystyle \mathbf {E} _{s}=\sum _{n=1}^{\infty }E_{n}\left(ia_{n}\mathbf {N} _{e1n}^{(3)}(k,\mathbf {r})-b_{n}\mathbf {M} _{o1n}^{(3)}(k,\mathbf {r})\right)}
H s = k ω μ ∑ n = 1 ∞ E n ( a n M e 1 n ( 3) ( k, r) + i b n N o 1 n ( 3) ( k, r)) {\displaystyle \mathbf {H} _{s}={\frac {k}{\omega \mu }}\sum _{n=1}^{\infty }E_{n}\left(a_{n}\mathbf {M} _{e1n}^{(3)}(k,\mathbf {r})+ib_{n}\mathbf {N} _{o1n}^{(3)}(k,\mathbf {r})\right)}{\displaystyle \mathbf {H} _{s}={\frac {k}{\omega \mu }}\sum _{n=1}^{\infty }E_{n}\left(a_{n}\mathbf {M} _{e1n}^{(3)}(k,\mathbf {r})+ib_{n}\mathbf {N} _{o1n}^{(3)}(k,\mathbf {r})\right)}

here the superscript ( 3) {\displaystyle (3)}(3)means that in the radial part of the functions ψ o e m n {\displaystyle \psi _{^{e}_{o}mn}}{\displaystyle \psi _{^{e}_{o}mn}}are spherical Hankel functions, and E n = i n E 0 ( 2 n + 1) n ( n + 1) {\displaystyle E_{n}={\frac {i^{n}E_{0}(2n+1)}{n(n+1)}}}{\displaystyle E_{n}={\frac {i^{n}E_{0}(2n+1)}{n(n+1)}}},

Internal fields:

E 1 = ∑ n = 1 ∞ E n ( − i d n N e 1 n ( 1) ( k 1, r) + c n M o 1 n ( 1) ( k 1, r)) {\displaystyle \mathbf {E} _{1}=\sum _{n=1}^{\infty }E_{n}\left(-id_{n}\mathbf {N} _{e1n}^{(1)}(k_{1},\mathbf {r})+c_{n}\mathbf {M} _{o1n}^{(1)}(k_{1},\mathbf {r})\right)}{\displaystyle \mathbf {E} _{1}=\sum _{n=1}^{\infty }E_{n}\left(-id_{n}\mathbf {N} _{e1n}^{(1)}(k_{1},\mathbf {r})+c_{n}\mathbf {M} _{o1n}^{(1)}(k_{1},\mathbf {r})\right)}
H 1 = − k 1 ω μ 1 ∑ n = 1 ∞ E n ( d n M e 1 n ( 1) ( k 1, r) + i c n N o 1 n ( 1) ( k 1, r)) {\displaystyle \mathbf {H} _{1}={\frac {-k_{1}}{\omega \mu _{1}}}\sum _{n=1}^{\infty }E_{n}\left(d_{n}\mathbf {M} _{e1n}^{(1)}(k_{1},\mathbf {r})+ic_{n}\mathbf {N} _{o1n}^{(1)}(k_{1},\mathbf {r})\right)}{\displaystyle \mathbf {H} _{1}={\frac {-k_{1}}{\omega \mu _{1}}}\sum _{n=1}^{\infty }E_{n}\left(d_{n}\mathbf {M} _{e1n}^{(1)}(k_{1},\mathbf {r})+ic_{n}\mathbf {N} _{o1n}^{(1)}(k_{1},\mathbf {r})\right)}

k = ω c n {\displaystyle k={\frac {\omega }{c}}n}{\displaystyle k={\frac {\omega }{c}}n}is the wave vector outside the particle k 1 = ω c n 1 {\displaystyle k_{1}={\frac {\omega }{c}}{n_{1}}}{\displaystyle k_{1}={\frac {\omega }{c}}{n_{1}}}is the wave vector in the medium from the particle material, n {\displaystyle n}nand n 1 {\displaystyle n_{1}}n_{1}are the refractive indices of the medium and the particle,

After applying the interface conditions, we obtain expressions for the coefficients:

c n ( ω) = μ 1 [ ρ h n ( ρ) ] ′ j n ( ρ) − μ 1 [ ρ j n ( ρ) ] ′ h n ( ρ) μ 1 [ ρ h n ( ρ) ] ′ j n ( ρ 1) − μ [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ) {\displaystyle c_{n}(\omega)={\frac {\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho)-\mu _{1}\left[\rho j_{n}(\rho)\right]'h_{n}(\rho)}{\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}}}{\displaystyle c_{n}(\omega)={\frac {\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho)-\mu _{1}\left[\rho j_{n}(\rho)\right]'h_{n}(\rho)}{\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}}}
d n ( ω) = μ 1 n 1 n [ ρ h n ( ρ) ] ′ j n ( ρ) − μ 1 n 1 n [ ρ j n ( ρ) ] ′ h n ( ρ) μ n 1 2 [ ρ h n ( ρ) ] ′ j n ( ρ 1) − μ 1 n 2 [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ), {\displaystyle d_{n}(\omega)={\frac {\mu _{1}n_{1}n\left[\rho h_{n}(\rho)\right]'j_{n}(\rho)-\mu _{1}n_{1}n\left[\rho j_{n}(\rho)\right]'h_{n}(\rho)}{\mu n_{1}^{2}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu _{1}n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}},}{\displaystyle d_{n}(\omega)={\frac {\mu _{1}n_{1}n\left[\rho h_{n}(\rho)\right]'j_{n}(\rho)-\mu _{1}n_{1}n\left[\rho j_{n}(\rho)\right]'h_{n}(\rho)}{\mu n_{1}^{2}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu _{1}n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}},}
b n ( ω) = μ 1 [ ρ j n ( ρ) ] ′ j n ( ρ 1) − μ [ ρ 1 j n ( ρ 1) ] ′ j n ( ρ) μ 1 [ ρ h n ( ρ) ] ′ j n ( ρ 1) − μ [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ) {\displaystyle b_{n}(\omega)={\frac {\mu _{1}\left[\rho j_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)}{\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}}}{\displaystyle b_{n}(\omega)={\frac {\mu _{1}\left[\rho j_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)}{\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}}}
a n ( ω) = μ n 1 2 [ ρ j n ( ρ) ] ′ j n ( ρ 1) − μ 1 n 2 [ ρ 1 j n ( ρ 1) ] ′ j n ( ρ) μ n 1 2 [ ρ h n ( ρ) ] ′ j n ( ρ 1) − μ 1 n 2 [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ), {\displaystyle a_{n}(\omega)={\frac {\mu n_{1}^{2}\left[\rho j_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu _{1}n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)}{\mu n_{1}^{2}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu _{1}n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}},}{\displaystyle a_{n}(\omega)={\frac {\mu n_{1}^{2}\left[\rho j_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu _{1}n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)}{\mu n_{1}^{2}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu _{1}n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}},}

where

ρ = k a {\displaystyle \rho =ka}{\displaystyle \rho =ka}
ρ 1 = k 1 a {\displaystyle \rho _{1}=k_{1}a}{\displaystyle \rho _{1}=k_{1}a}with a {\displaystyle a}abeing the radius of the sphere.

j n {\displaystyle j_{n}}{\displaystyle j_{n}}and h n {\displaystyle h_{n}}h_{n}represent the spherical functions of Bessel and Hankel of the first kind, respectively.

Scattering and extinction cross-sections

Multipole decomposition spectrum of scattering cross-section by gold nanosphere with radius 100 nm Multipole decomposition spectrum of scattering cross-section by nanosphere with radius 100 nm and refractive index n=4 Multipole decomposition spectrum of scattering cross-section by silicon nanosphere with radius 100 nm

Values commonly calculated using Mie theory include efficiency coefficients for extinction Q e {\displaystyle Q_{e}}Q_e, scattering Q s {\displaystyle Q_{s}}Q_{s}, and absorption Q a {\displaystyle Q_{a}}Q_a. These efficiency coefficients are ratios of the cross section of the respective process, σ i {\displaystyle \sigma _{i}}\sigma _{i}, to the particle protected area, Q i = σ i π a 2 {\displaystyle Q_{i}={\frac {\sigma _{i}}{\pi a^{2}}}}{\displaystyle Q_{i}={\frac {\sigma _{i}}{\pi a^{2}}}}, where a is the particle radius. According to the definition of extinction,

σ e = σ s + σ a {\displaystyle \sigma _{e}=\sigma _{s}+\sigma _{a}}{\displaystyle \sigma _{e}=\sigma _{s}+\sigma _{a}}and Q e = Q s + Q a {\displaystyle Q_{e}=Q_{s}+Q_{a}}{\displaystyle Q_{e}=Q_{s}+Q_{a}}.

The scattering and extinction coefficients can be represented as the infinite series:

Q s = 2 π k 2 ∑ n = 1 ∞ ( 2 n + 1) ( | a n | 2 + | b n | 2) {\displaystyle Q_{s}={\frac {2\pi }{k^{2}}}\sum _{n=1}^{\infty }(2n+1)(|a_{n}|^{2}+|b_{n}|^{2})}{\displaystyle Q_{s}={\frac {2\pi }{k^{2}}}\sum _{n=1}^{\infty }(2n+1)(|a_{n}|^{2}+|b_{n}|^{2})}
Q e = 2 π k 2 ∑ n = 1 ∞ ( 2 n + 1) ℜ ( a n + b n) {\displaystyle Q_{e}={\frac {2\pi }{k^{2}}}\sum _{n=1}^{\infty }(2n+1)\Re (a_{n}+b_{n})}{\displaystyle Q_{e}={\frac {2\pi }{k^{2}}}\sum _{n=1}^{\infty }(2n+1)\Re (a_{n}+b_{n})}

.

.

Application to subwavelength particles

If the size of the particle is equal to several wavelengths in the material, then the scattered fields have some features. Further, we will talk about the form of the electric field since the magnetic field is obtained from it by taking the rotor.

All Mie coefficients depend on the frequency and have maximums when the denominator is close to zero (exact equality to zero is achieved for complex frequencies). In this case, it is possible, that the contribution of one specific harmonic dominates in scattering. Then at large distances from the particle, the radiation pattern of the scattered field will be similar to the corresponding radiation pattern of the angular part of vector spherical harmonics. The harmonics N o e m 1 {\displaystyle \mathbf {N} _{^{e}_{o}m1}}{\displaystyle \mathbf {N} _{^{e}_{o}m1}}correspond to electric dipoles (if the contribution of this harmonic dominates in the expansion of the electric field, then the field is similar to the electric dipole field), M o e m 1 {\displaystyle \mathbf {M} _{^{e}_{o}m1}}{\displaystyle \mathbf {M} _{^{e}_{o}m1}}correspond to the electric field of the magnetic dipole, N o e m 2 {\displaystyle \mathbf {N} _{^{e}_{o}m2}}{\displaystyle \mathbf {N} _{^{e}_{o}m2}}and M o e m 2 {\displaystyle \mathbf {M} _{^{e}_{o}m2}}{\displaystyle \mathbf {M} _{^{e}_{o}m2}}- electric and magnetic quadrupoles, N o e m 3 {\displaystyle \mathbf {N} _{^{e}_{o}m3}}{\displaystyle \mathbf {N} _{^{e}_{o}m3}}and M o e m 3 {\displaystyle \mathbf {M} _{^{e}_{o}m3}}{\displaystyle \mathbf {M} _{^{e}_{o}m3}}- octupoles, and so on. The maxima of the scattering coefficients (as well as the change of their phase to π {\displaystyle \pi }\pi ) are called multipole resonances.

The dependence of the scattering cross-section on the wavelength and the contribution of specific resonances strongly depends on the particle material. For example, for a gold particle with a radius of 100 nm, the contribution of the electric dipole to scattering predominates in the optical range, while for a silicon particle there are pronounced magnetic dipole and quadrupole resonances. For metal particles, the peak visible in the scattering cross-section is also called localized plasmon resonance.

In the limit of small particles or long wavelengths, the electric dipole contribution dominates in the scattering cross-section.

Other directions of the incident plane wave

In case of x-polarized plane wave, incident along the z-axis, decompositions of all fields contained only harmonics with m=1, but for an arbitrary incident wave this is not the case. For a rotated plane wave, the expansion coefficients can be obtained, for example, using the fact that during rotation, vector spherical harmonics are transformed through each other by Wigner D-matrixes.

In this case, the scattered field will be decomposed by all possible harmonics:

E s = ∑ n = 1 ∞ ∑ m = 0 n E 0 ( D M e m n M e m n ( 3) ( k, r) + D M o m n M o m n ( 3) ( k, r) + D N e m n N e m n ( 3) ( k, r) + D N o m n N o m n ( 3) ( k, r)) {\displaystyle \mathbf {E} _{s}=\sum _{n=1}^{\infty }\sum _{m=0}^{n}E_{0}(D_{Memn}\mathbf {M} _{emn}^{(3)}(k,\mathbf {r})+D_{Momn}\mathbf {M} _{omn}^{(3)}(k,\mathbf {r})+D_{Nemn}\mathbf {N} _{emn}^{(3)}(k,\mathbf {r})+D_{Nomn}\mathbf {N} _{omn}^{(3)}(k,\mathbf {r}))}{\displaystyle \mathbf {E} _{s}=\sum _{n=1}^{\infty }\sum _{m=0}^{n}E_{0}(D_{Memn}\mathbf {M} _{emn}^{(3)}(k,\mathbf {r})+D_{Momn}\mathbf {M} _{omn}^{(3)}(k,\mathbf {r})+D_{Nemn}\mathbf {N} _{emn}^{(3)}(k,\mathbf {r})+D_{Nomn}\mathbf {N} _{omn}^{(3)}(k,\mathbf {r}))}

Then the scattering cross section will be expressed in terms of the coefficients as follows:

C s c a = 2 π π a 2 k 2 ∑ n = 1 ∞ n ( n + 1) ( 2 n + 1) × [ ∑ m = 1 n ( n + m) ! ( n − m) ! ( | D M e m n | 2 + | D M o m n | 2 + | D N e m n | 2 + | D N o m n | 2) + 2 | D M e 0 n | 2 + 2 | D N e 0 n | 2 ]. {\displaystyle C_{sca}={\frac {2\pi }{\pi a^{2}k^{2}}}\sum _{n=1}^{\infty }{\frac {n(n+1)}{(2n+1)}}\times {\Bigl [}\sum \limits _{m=1}^{n}{\frac {(n+m)!}{(n-m)!}}(|D_{Memn}|^{2}+|D_{Momn}|^{2}+|D_{Nemn}|^{2}+|D_{Nomn}|^{2})+2|D_{Me0n}|^{2}+2|D_{Ne0n}|^{2}{\Bigr ]}.}{\displaystyle C_{sca}={\frac {2\pi }{\pi a^{2}k^{2}}}\sum _{n=1}^{\infty }{\frac {n(n+1)}{(2n+1)}}\times {\Bigl [}\sum \limits _{m=1}^{n}{\frac {(n+m)!}{(n-m)!}}(|D_{Memn}|^{2}+|D_{Momn}|^{2}+|D_{Nemn}|^{2}+|D_{Nomn}|^{2})+2|D_{Me0n}|^{2}+2|D_{Ne0n}|^{2}{\Bigr ]}.}

Kerker effect

Kerker effect is a phenomenon in scattering directionality, which occurs when different multipole responses are presented and not negligible.

Particular (dipolar) case of the Kerker effect. The total electric field of the crossed magnetic and electric dipoles radiating in phase. The radiation pattern is asymmetric, in one direction the fields are mutually destroyed, and in the other, they add up.

In 1983, in the work of Kerker, Wang and Giles the direction of scattering by particles with μ ≠ 1 {\displaystyle \mu \neq 1}{\displaystyle \mu \neq 1}was investigated. In particular, it was shown that for hypothetical particles with μ = ε {\displaystyle \mu =\varepsilon }{\displaystyle \mu =\varepsilon }backward scattering is completely suppressed. This can be seen as an extension to a spherical surface of Giles' and Wild's results for reflection at a planar surface with equal refractive indices where reflection and transmission is constant and independent of angle of incidence.

. In addition, scattering cross sections in the forward and backward directions are simply expressed in terms of Mie coefficients:

C s c a b a c k w a r d = 1 a 2 k 2 | ∑ n = 1 ∞ ( 2 n + 1) ( − 1) n ( a n − b n) | 2 {\displaystyle C_{sca}^{backward}={\frac {1}{a^{2}k^{2}}}{\bigg |}\sum _{n=1}^{\infty }{(2n+1)}(-1)^{n}(a_{n}-b_{n}){\bigg |}^{2}}{\displaystyle C_{sca}^{backward}={\frac {1}{a^{2}k^{2}}}{\bigg |}\sum _{n=1}^{\infty }{(2n+1)}(-1)^{n}(a_{n}-b_{n}){\bigg |}^{2}}
C s c a f o r w a r d = 1 a 2 k 2 | ∑ n = 1 ∞ ( 2 n + 1) ( a n + b n) | 2 {\displaystyle C_{sca}^{forward}={\frac {1}{a^{2}k^{2}}}{\bigg |}\sum _{n=1}^{\infty }{(2n+1)}(a_{n}+b_{n}){\bigg |}^{2}}{\displaystyle C_{sca}^{forward}={\frac {1}{a^{2}k^{2}}}{\bigg |}\sum _{n=1}^{\infty }{(2n+1)}(a_{n}+b_{n}){\bigg |}^{2}}

For certain combinations of coefficients, the expressions above can be minimized.

So, for example, when terms with n>1 {\displaystyle n>1}n>1can be neglected (dipole approximation), ( a 1 − b 1) = 0 {\displaystyle (a_{1}-b_{1})=0}{\displaystyle (a_{1}-b_{1})=0}, corresponds to the minimum in backscattering (magnetic and electric dipoles are equal in magnitude and are in phase, this is also called ‘first Kerker’ or ‘zero-backward intensity condition’). And ( a 1 + b 1) = 0 {\displaystyle (a_{1}+b_{1})=0}{\displaystyle (a_{1}+b_{1})=0}corresponds to minimum in forward scattering, this is also called ‘second Kerker condition’ (or ‘near-zero forward intensity condition’). For the exact solution of the problem, it is necessary to take into account the contributions of all multipoles. The sum of the electric and magnetic dipoles forms Huygens source

For dielectric particles, maximum forward scattering is observed at wavelengths longer than the wavelength of magnetic dipole resonance, and maximum backward scattering at shorter ones..

Later, other varieties of the effect were found. For example, the transverse Kerker effect, with nearly complete simultaneous suppression of both forward and backward scattered fields (side-scattering patterns), optomechanical Kerker effect, in acoustic scattering, and also found in plants.

There is also a short Video on YouTube with an explanation of the effect.

Dyadic Green's function of a sphere

Green's function is a solution to the following equation:

∇ × ∇ × G ^ ( ω, r, r ′) = ( ω c) 2 ε ( r, ω) G ^ ( ω, r, r ′) + 1 ^ δ ( r − r ′), {\displaystyle \nabla \times \nabla \times {\bf {\hat {G}}}(\omega,\mathbf {r},\mathbf {r} ')=\left({\frac {\omega }{c}}\right)^{2}\varepsilon (\mathbf {r},\omega){\bf {\hat {G}}}(\omega,\mathbf {r},\mathbf {r} ')+{\bf {\hat {1}}}\delta (\mathbf {r} -\mathbf {r} '),}{\displaystyle \nabla \times \nabla \times {\bf {\hat {G}}}(\omega,\mathbf {r},\mathbf {r} ')=\left({\frac {\omega }{c}}\right)^{2}\varepsilon (\mathbf {r},\omega){\bf {\hat {G}}}(\omega,\mathbf {r},\mathbf {r} ')+{\bf {\hat {1}}}\delta (\mathbf {r} -\mathbf {r} '),}

where 1 ^ {\displaystyle {\hat {\bf {1}}}}{\displaystyle {\hat {\bf {1}}}}— identity matrix ε ( r, ω) = ε 1 ( ω) {\displaystyle \varepsilon (\mathbf {r},\omega)=\varepsilon _{1}(\omega)}{\displaystyle \varepsilon (\mathbf {r},\omega)=\varepsilon _{1}(\omega)}для r < a {\displaystyle r{\displaystyle r<a}, and ε ( r, ω) = ε {\displaystyle \varepsilon (\mathbf {r},\omega)=\varepsilon }{\displaystyle \varepsilon (\mathbf {r},\omega)=\varepsilon }for r>a {\displaystyle r>a}{\displaystyle r>a}. Since all fields are vectorial, the Green function is a 3 by 3 matrix and is called a dyadic. If polarization P ( r) {\displaystyle \mathbf {P} (\mathbf {r})}{\displaystyle \mathbf {P} (\mathbf {r})}is induced in the system, when the fields are written as

E ω ( r) = ω 2 μ ∫ V d V ′ G ^ ( r, r ′, k) P ω ( r ′) {\displaystyle \mathbf {E} ^{\omega }({\mathbf {r} })=\omega ^{2}\mu \int \limits _{V}dV'{\hat {\bf {G}}}({\bf {r,r'}},k)\mathbf {P} ^{\omega }(\mathbf {r} ')}{\displaystyle \mathbf {E} ^{\omega }({\mathbf {r} })=\omega ^{2}\mu \int \limits _{V}dV'{\hat {\bf {G}}}({\bf {r,r'}},k)\mathbf {P} ^{\omega }(\mathbf {r} ')}

In the same way as the fields, the Green's function can be decomposed into vector spherical harmonics. Dyadic Green's function of a free space а:

G ^ 0 ( r, r ′, k) = e r ⊗ e r k 2 δ ( r − r ′) + i k 4 π ∑ n = 1 ∞ ∑ m = 0 n ( 2 − δ m, 0) 2 n + 1 n ( n + 1) ( n − m) ! ( n + m) ! ⋅ {\displaystyle {\hat {\bf {G}}}^{0}({\mathbf {r},\mathbf {r} ',k})={\frac {\mathbf {e_{r}} \otimes \mathbf {e_{r}} }{k^{2}}}\delta (\mathbf {r} -\mathbf {r} ')+{\frac {ik}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }{\displaystyle {\hat {\bf {G}}}^{0}({\mathbf {r},\mathbf {r} ',k})={\frac {\mathbf {e_{r}} \otimes \mathbf {e_{r}} }{k^{2}}}\delta (\mathbf {r} -\mathbf {r} ')+{\frac {ik}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }
{ ⋅ ( ( M e m n ( 1) [ k, r ] ⊗ M e m n ( 3) [ k, r ′ ] + M o m n ( 1) [ k, r ] ⊗ M o m n ( 3) [ k, r ′ ]) + ( N e m n ( 1) [ k, r ] ⊗ N e m n ( 3) [ k, r ′ ] + N o m n ( 1) [ k, r ] ⊗ N o m n ( 3) [ k, r ′ ])), if r < r ′ ⋅ ( ( M e m n ( 3) [ k, r ] ⊗ M e m n ( 1) [ k, r ′ ] + M o m n ( 3) [ k, r ] ⊗ M o m n ( 1) [ k, r ′ ]) + ( N e m n ( 3) [ k, r ] ⊗ N e m n ( 1) [ k, r ′ ] + N o m n ( 3) [ k, r ] ⊗ N o m n ( 1) [ k, r ′ ])), if r>r ′ {\displaystyle \left\{{\begin{array}{l}\cdot {\Bigl (}(\mathbf {M} _{emn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{emn}^{(3)}[k,\mathbf {r} ']+\mathbf {M} _{omn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{omn}^{(3)}[k,\mathbf {r} '])+({\mathbf {N} }_{emn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{emn}^{(3)}[k,\mathbf {r} ']+\mathbf {N} _{omn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{omn}^{(3)}[k,\mathbf {r} ']){\Bigr)},{\text{if }}rr'\end{array}}\right.}{\displaystyle \left\{{\begin{array}{l}\cdot {\Bigl (}(\mathbf {M} _{emn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{emn}^{(3)}[k,\mathbf {r} ']+\mathbf {M} _{omn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{omn}^{(3)}[k,\mathbf {r} '])+({\mathbf {N} }_{emn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{emn}^{(3)}[k,\mathbf {r} ']+\mathbf {N} _{omn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{omn}^{(3)}[k,\mathbf {r} ']){\Bigr)},{\text{if }}r<r'\\\cdot {\Bigl (}(\mathbf {M} _{emn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{emn}^{(1)}[k,\mathbf {r} ']+\mathbf {M} _{omn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{omn}^{(1)}[k,\mathbf {r} '])+({\mathbf {N} }_{emn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{emn}^{(1)}[k,\mathbf {r} ']+\mathbf {N} _{omn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{omn}^{(1)}[k,\mathbf {r} ']){\Bigr)},{\text{if }}r>r'\end{array}}\right.}

In the presence of a sphere, the Green's function is also decomposed into vector spherical harmonics. Its appearance depends on the environment in which the points r {\displaystyle \mathbf {r} }\mathbf {r} and r ′ {\displaystyle \mathbf {r} '}\mathbf {r} 'are located.

When both points are outside the sphere (r>a, r ′>a {\displaystyle r>a,r'>a}{\displaystyle r>a,r'>a}):

G ^ 00 ( r, r ′, k, k 1) = G ^ 0 ( r, r ′, k) + i k 4 π ∑ n = 1 ∞ ∑ m = 0 n ( 2 − δ m, 0) 2 n + 1 n ( n + 1) ( n − m) ! ( n + m) ! ⋅ {\displaystyle {\hat {\bf {G}}}^{00}({\mathbf {r},\mathbf {r} ',k,k_{1}})={\hat {\bf {G}}}^{0}({\mathbf {r},\mathbf {r} ',k})+{\frac {ik}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }{\displaystyle {\hat {\bf {G}}}^{00}({\mathbf {r},\mathbf {r} ',k,k_{1}})={\hat {\bf {G}}}^{0}({\mathbf {r},\mathbf {r} ',k})+{\frac {ik}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }
⋅ ( a n ( 0) ( ω) ( M o e m n ( 3) [ k, r ] ⊗ M o e m n ( 3) [ k, r ′ ]) + b n ( 0) ( ω) ( N o e m n ( 3) [ k, r ] ⊗ N o e m n ( 3) [ k, r ′ ])) {\displaystyle \cdot {\Bigl (}a_{n}^{(0)}(\omega)(\mathbf {M} _{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{^{e}_{o}mn}^{(3)}[k,\mathbf {r} '])+b_{n}^{(0)}(\omega)({\mathbf {N} }_{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ']){\Bigr)}}{\displaystyle \cdot {\Bigl (}a_{n}^{(0)}(\omega)(\mathbf {M} _{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{^{e}_{o}mn}^{(3)}[k,\mathbf {r} '])+b_{n}^{(0)}(\omega)({\mathbf {N} }_{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ']){\Bigr)}}

where the coefficients are :

a n ( 0) ( ω) = μ / μ 1 [ ρ 1 j n ( ρ 1) ] ′ j n ( ρ) − [ ρ j n ( ρ) ] ′ j n ( ρ 1) [ ρ h n ( ρ) ] ′ j n ( ρ 1) − μ / μ 1 [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ), {\displaystyle a_{n}^{(0)}(\omega)={\frac {\mu /\mu _{1}\left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)-\left[\rho j_{n}(\rho)\right]'j_{n}(\rho _{1})}{\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu /\mu _{1}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}},}{\displaystyle a_{n}^{(0)}(\omega)={\frac {\mu /\mu _{1}\left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)-\left[\rho j_{n}(\rho)\right]'j_{n}(\rho _{1})}{\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu /\mu _{1}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}},}
b n ( 0) ( ω) = n 2 μ 1 / μ [ ρ 1 j n ( ρ 1) ] ′ j n ( ρ) − n 1 2 [ ρ j n ( ρ) ] ′ j n ( ρ 1) n 1 2 [ ρ h n ( ρ) ] ′ j n ( ρ 1) − n 2 μ 1 / μ [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ). {\displaystyle b_{n}^{(0)}(\omega)={\frac {n^{2}\mu _{1}/\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)-n_{1}^{2}\left[\rho j_{n}(\rho)\right]'j_{n}(\rho _{1})}{n_{1}^{2}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-n^{2}\mu _{1}/\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}}.}{\displaystyle b_{n}^{(0)}(\omega)={\frac {n^{2}\mu _{1}/\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)-n_{1}^{2}\left[\rho j_{n}(\rho)\right]'j_{n}(\rho _{1})}{n_{1}^{2}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-n^{2}\mu _{1}/\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}}.}

When both points are inside the sphere (r < a, r ′ < a {\displaystyle r{\displaystyle r<a,r'<a}) :

G ^ 11 ( r, r ′, k, k 1) = G ^ 0 ( r, r ′, k 1) + i k 1 4 π ∑ n = 1 ∞ ∑ m = 0 n ( 2 − δ m, 0) 2 n + 1 n ( n + 1) ( n − m) ! ( n + m) ! ⋅ {\displaystyle {\hat {\bf {G}}}^{11}({\mathbf {r},\mathbf {r} ',k,k_{1}})={\hat {\bf {G}}}^{0}({\mathbf {r},\mathbf {r} ',k_{1}})+{\frac {ik_{1}}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }{\displaystyle {\hat {\bf {G}}}^{11}({\mathbf {r},\mathbf {r} ',k,k_{1}})={\hat {\bf {G}}}^{0}({\mathbf {r},\mathbf {r} ',k_{1}})+{\frac {ik_{1}}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }
⋅ ( c n ( 1) ( ω) ( M o e m n ( 1) [ k 1, r ] ⊗ M o e m n ( 1) [ k 1, r ′ ]) + d n ( 1) ( ω) ( N o e m n ( 1) [ k 1, r ] ⊗ N o e m n ( 1) [ k 1, r ′ ])), {\displaystyle \cdot {\Bigl (}c_{n}^{(1)}(\omega)(\mathbf {M} _{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} ]\otimes {\mathbf {M} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} '])+d_{n}^{(1)}(\omega)({\mathbf {N} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} ]\otimes {\mathbf {N} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} ']){\Bigr)},}{\displaystyle \cdot {\Bigl (}c_{n}^{(1)}(\omega)(\mathbf {M} _{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} ]\otimes {\mathbf {M} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} '])+d_{n}^{(1)}(\omega)({\mathbf {N} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} ]\otimes {\mathbf {N} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} ']){\Bigr)},}

Coefficients:

c n ( 1) ( ω) = μ 1 / μ [ ρ h n ( ρ) ] ′ h n ( ρ 1) − [ ρ 1 h n ( ρ 1) ] ′ h n ( ρ) [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ) − μ 1 / μ [ ρ h n ( ρ) ] ′ j n ( ρ 1), {\displaystyle c_{n}^{(1)}(\omega)={\frac {\mu _{1}/\mu \left[\rho h_{n}(\rho)\right]'h_{n}(\rho _{1})-\left[\rho _{1}h_{n}(\rho _{1})\right]'h_{n}(\rho)}{\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)-\mu _{1}/\mu \left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})}},}{\displaystyle c_{n}^{(1)}(\omega)={\frac {\mu _{1}/\mu \left[\rho h_{n}(\rho)\right]'h_{n}(\rho _{1})-\left[\rho _{1}h_{n}(\rho _{1})\right]'h_{n}(\rho)}{\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)-\mu _{1}/\mu \left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})}},}
d n ( 1) ( ω) = n 1 2 μ / μ 1 [ ρ h n ( ρ) ] ′ h n ( ρ 1) − n 2 [ ρ 1 h n ( ρ 1) ] ′ h n ( ρ) n 2 [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ) − n 1 2 μ / μ 1 [ ρ h n ( ρ) ] ′ j n ( ρ 1). {\displaystyle d_{n}^{(1)}(\omega)={\frac {n_{1}^{2}\mu /\mu _{1}\left[\rho h_{n}(\rho)\right]'h_{n}(\rho _{1})-n^{2}\left[\rho _{1}h_{n}(\rho _{1})\right]'h_{n}(\rho)}{n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)-n_{1}^{2}\mu /\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})}}.}{\displaystyle d_{n}^{(1)}(\omega)={\frac {n_{1}^{2}\mu /\mu _{1}\left[\rho h_{n}(\rho)\right]'h_{n}(\rho _{1})-n^{2}\left[\rho _{1}h_{n}(\rho _{1})\right]'h_{n}(\rho)}{n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)-n_{1}^{2}\mu /\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})}}.}

Source is inside the sphere and observation point is outside (r>a, r ′ < a {\displaystyle r>a,r'{\displaystyle r>a,r'<a}) :

G ^ 01 ( r, r ′, k, k 1) = i k 1 4 π ∑ n = 1 ∞ ∑ m = 0 n ( 2 − δ m, 0) 2 n + 1 n ( n + 1) ( n − m) ! ( n + m) ! ⋅ {\displaystyle {\hat {\bf {G}}}^{01}({\mathbf {r},\mathbf {r} ',k,k_{1}})={\frac {ik_{1}}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }{\displaystyle {\hat {\bf {G}}}^{01}({\mathbf {r},\mathbf {r} ',k,k_{1}})={\frac {ik_{1}}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }
⋅ ( a n ( 1) ( ω) ( M o e m n ( 3) [ k, r ] ⊗ M o e m n ( 1) [ k 1, r ′ ]) + b n ( 1) ( ω) ( N o e m n ( 3) [ k, r ] ⊗ N o e m n ( 1) [ k 1, r ′ ])) {\displaystyle \cdot {\Bigl (}a_{n}^{(1)}(\omega)(\mathbf {M} _{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} '])+b_{n}^{(1)}(\omega)({\mathbf {N} }_{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} ']){\Bigr)}}{\displaystyle \cdot {\Bigl (}a_{n}^{(1)}(\omega)(\mathbf {M} _{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} '])+b_{n}^{(1)}(\omega)({\mathbf {N} }_{^{e}_{o}mn}^{(3)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{^{e}_{o}mn}^{(1)}[k_{1},\mathbf {r} ']){\Bigr)}}

coefficients:

a n ( 1) ( ω) = [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ 1) − [ ρ 1 h n ( ρ 1) ] ′ j n ( ρ 1) [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ) − μ 1 / μ [ ρ h n ( ρ) ] ′ j n ( ρ 1), {\displaystyle a_{n}^{(1)}(\omega)={\frac {\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho _{1})-\left[\rho _{1}h_{n}(\rho _{1})\right]'j_{n}(\rho _{1})}{\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)-\mu _{1}/\mu \left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})}},}{\displaystyle a_{n}^{(1)}(\omega)={\frac {\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho _{1})-\left[\rho _{1}h_{n}(\rho _{1})\right]'j_{n}(\rho _{1})}{\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)-\mu _{1}/\mu \left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})}},}
b n ( 1) ( ω) = n n 1 [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ 1) − n n 1 [ ρ 1 h n ( ρ 1) ] ′ j n ( ρ 1) n 2 μ 1 / μ [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ) − n 1 2 [ ρ h n ( ρ) ] ′ j n ( ρ 1). {\displaystyle b_{n}^{(1)}(\omega)={\frac {nn_{1}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho _{1})-nn_{1}\left[\rho _{1}h_{n}(\rho _{1})\right]'j_{n}(\rho _{1})}{n^{2}\mu _{1}/\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)-n_{1}^{2}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})}}.}{\displaystyle b_{n}^{(1)}(\omega)={\frac {nn_{1}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho _{1})-nn_{1}\left[\rho _{1}h_{n}(\rho _{1})\right]'j_{n}(\rho _{1})}{n^{2}\mu _{1}/\mu \left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)-n_{1}^{2}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})}}.}

Source is outside the sphere and observation point is inside (r < a, r ′>a {\displaystyle ra}{\displaystyle r<a,r'>a}) :

G ^ 10 ( r, r ′, k, k 1) = i k 4 π ∑ n = 1 ∞ ∑ m = 0 n ( 2 − δ m, 0) 2 n + 1 n ( n + 1) ( n − m) ! ( n + m) ! ⋅ {\displaystyle {\hat {\bf {G}}}^{10}({\mathbf {r},\mathbf {r} ',k,k_{1}})={\frac {ik}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }{\displaystyle {\hat {\bf {G}}}^{10}({\mathbf {r},\mathbf {r} ',k,k_{1}})={\frac {ik}{4\pi }}\sum _{n=1}^{\infty }\sum _{m=0}^{n}(2-\delta _{m,0}){\frac {2n+1}{n(n+1)}}{\frac {(n-m)!}{(n+m)!}}\cdot }
⋅ ( c n ( 0) ( ω) ( M o e m n ( 1) [ k, r ] ⊗ M o e m n ( 3) [ k 1, r ′ ]) + d n ( 0) ( ω) ( N o e m n ( 1) [ k, r ] ⊗ N o e m n ( 3) [ k 1, r ′ ])) {\displaystyle \cdot {\Bigl (}c_{n}^{(0)}(\omega)(\mathbf {M} _{^{e}_{o}mn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{^{e}_{o}mn}^{(3)}[k_{1},\mathbf {r} '])+d_{n}^{(0)}(\omega)({\mathbf {N} }_{^{e}_{o}mn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{^{e}_{o}mn}^{(3)}[k_{1},\mathbf {r} ']){\Bigr)}}{\displaystyle \cdot {\Bigl (}c_{n}^{(0)}(\omega)(\mathbf {M} _{^{e}_{o}mn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {M} }_{^{e}_{o}mn}^{(3)}[k_{1},\mathbf {r} '])+d_{n}^{(0)}(\omega)({\mathbf {N} }_{^{e}_{o}mn}^{(1)}[k,\mathbf {r} ]\otimes {\mathbf {N} }_{^{e}_{o}mn}^{(3)}[k_{1},\mathbf {r} ']){\Bigr)}}

coefficients:

c n ( 0) ( ω) = [ ρ h n ( ρ) ] ′ j n ( ρ) − [ ρ j n ( ρ) ] ′ h n ( ρ) [ ρ h n ( ρ) ] ′ j n ( ρ 1) − μ / μ 1 [ ρ 1 j n ( ρ 1) ] ′ h n ( ρ), {\displaystyle c_{n}^{(0)}(\omega)={\frac {\left[\rho h_{n}(\rho)\right]'j_{n}(\rho)-\left[\rho j_{n}(\rho)\right]'h_{n}(\rho)}{\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu /\mu _{1}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}},}{\displaystyle c_{n}^{(0)}(\omega)={\frac {\left[\rho h_{n}(\rho)\right]'j_{n}(\rho)-\left[\rho j_{n}(\rho)\right]'h_{n}(\rho)}{\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-\mu /\mu _{1}\left[\rho _{1}j_{n}(\rho _{1})\right]'h_{n}(\rho)}},}
d n ( 0) ( ω) = n n 1 [ ρ h n ( ρ) ] ′ j n ( ρ) − n n 1 [ ρ j n ( ρ) ] ′ h n ( ρ) n 1 2 μ / μ 1 [ ρ h n ( ρ) ] ′ j n ( ρ 1) − n 2 [ ρ 1 j n ( ρ 1) ] ′ j n ( ρ). {\displaystyle d_{n}^{(0)}(\omega)={\frac {nn_{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho)-nn_{1}\left[\rho j_{n}(\rho)\right]'h_{n}(\rho)}{n_{1}^{2}\mu /\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)}}.}{\displaystyle d_{n}^{(0)}(\omega)={\frac {nn_{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho)-nn_{1}\left[\rho j_{n}(\rho)\right]'h_{n}(\rho)}{n_{1}^{2}\mu /\mu _{1}\left[\rho h_{n}(\rho)\right]'j_{n}(\rho _{1})-n^{2}\left[\rho _{1}j_{n}(\rho _{1})\right]'j_{n}(\rho)}}.}

Computational codes

Mie solutions are implemented in a number of programs written in different computer languages such as Fortran, MATLAB, and Mathematica. These solutions solve for an infinite series, and provide as output the calculation of the scattering phase function, extinction, scattering, and absorption efficiencies, and other parameters such as asymmetry parameters or radiation torque. Current usage of the term "Mie solution" indicates a series approximation to a solution of Maxwell's equations. There are several known objects that allow such a solution: spheres, concentric spheres, infinite cylinders, clusters of spheres and clusters of cylinders. There are also known series solutions for scattering by ellipsoidal particles. A list of codes implementing these specialized solutions is provided in the following:

A generalization that allows a treatment of more generally shaped particles is the T-matrix method, which also relies on a series approximation to solutions of Maxwell's equations.

See also external links for other codes and calculators.

Applications

Mie theory is very important in meteorological optics, where diameter-to-wavelength ratios of the order of unity and larger are characteristic for many problems regarding haze and cloud scattering. A further application is in the characterization of particles by optical scattering measurements. The Mie solution is also important for understanding the appearance of common materials like milk, biological tissue and latex paint.

Atmospheric science

Mie scattering occurs when the diameters of atmospheric particulates are similar to the wavelengths of the scattered light. Dust, pollen, smoke and microscopic water droplets that form clouds are common causes of Mie scattering. Mie scattering occurs mostly in the lower portions of the atmosphere, where larger particles are more abundant, and dominates in cloudy conditions.

Cancer detection and screening

Mie theory has been used to determine whether scattered light from tissue corresponds to healthy or cancerous cell nuclei using angle-resolved low-coherence interferometry.

Clinical laboratory analysis

Mie theory is a central principle in the application of nephelometric based assays, widely used in medicine to measure various plasma proteins. A wide array of plasma proteins can be detected and quantified by nephelometry.

Magnetic particles

A number of unusual electromagnetic scattering effects occur for magnetic spheres. When the relative permittivity equals the permeability, the back-scatter gain is zero. Also, the scattered radiation is polarized in the same sense as the incident radiation. In the small-particle (or long-wavelength) limit, conditions can occur for zero forward scatter, for complete polarization of scattered radiation in other directions, and for asymmetry of forward scatter to backscatter. The special case in the small-particle limit provides interesting special instances of complete polarization and forward-scatter-to-backscatter asymmetry.

Metamaterial

Mie theory has been used to design metamaterials. They usually consist of three-dimensional composites of metal or non-metallic inclusions periodically or randomly embedded in a low-permittivity matrix. In such a scheme, the negative constitutive parameters are designed to appear around the Mie resonances of the inclusions: the negative effective permittivity is designed around the resonance of the Mie electric dipole scattering coefficient, whereas negative effective permeability is designed around the resonance of the Mie magnetic dipole scattering coefficient, and doubly negative material (DNG) is designed around the overlap of resonances of Mie electric and magnetic dipole scattering coefficients. The particle usually have the following combinations:

  1. one set of magnetodielectric particles with values of relative permittivity and permeability much greater than one and close to each other;
  2. two different dielectric particles with equal permittivity but different size;
  3. two different dielectric particles with equal size but different permittivity.

In theory, the particles analyzed by Mie theory are commonly spherical but, in practice, particles are usually fabricated as cubes or cylinders for ease of fabrication. To meet the criteria of homogenization, which may be stated in the form that the lattice constant is much smaller than the operating wavelength, the relative permittivity of the dielectric particles should be much greater than 1, e.g. ε r>78 ( 38) {\displaystyle \varepsilon _{\text{r}}>78(38)}{\displaystyle \varepsilon _{\text{r}}>78(38)}to achieve negative effective permittivity (permeability).

Particle sizing

Mie theory is often applied in laser diffraction analysis to inspect the particle sizing effect. While early computers in the 1970s were only able to compute diffraction data with the more simple Fraunhofer approximation, Mie is widely used since the 1990s and officially recommended for particles below 50 micrometers in guideline ISO 13321:2009.

Mie theory has been used in the detection of oil concentration in polluted water.

Mie scattering is the primary method of sizing single sonoluminescing bubbles of air in water and is valid for cavities in materials, as well as particles in materials, as long as the surrounding material is essentially non-absorbing.

Parasitology

It has also been used to study the structure of Plasmodium falciparum, a particularly pathogenic form of malaria.

Extensions

In 1986, P. A. Bobbert and J. Vlieger extended the Mie model to calculate scattering by a sphere in a homogeneous medium placed on flat surface. Like Mie model, the extended model can be applied to spheres with a radius close to the wavelength of the incident light. There is a C++ code implementing Bobbert–Vlieger (BV) model. Recent developments are related to scattering by ellipsoid. The contemporary studies go to well known research of Rayleigh.

See also

References

Further reading

  • Kerker, M. (1969). The scattering of light and other electromagnetic radiation. New York: Academic.
  • Barber, P. W.; Hill, S. S. (1990). Light scattering by particles: Computational methods. Singapore: World Scientific. ISBN 978-9971-5-0813-5.
  • Mishchenko, M.; Travis, L.; Lacis, A. (2002). Scattering, Absorption, and Emission of Light by Small Particles. New York: Cambridge University Press. ISBN 978-0-521-78252-4.
  • Frisvad, J.; Christensen, N.; Jensen, H. (2007). "Computing the Scattering Properties of Participating Media using Lorenz-Mie Theory". ACM Transactions on Graphics. 26(3): 60. doi :10.1145/1276377.1276452.
  • Wriedt, Thomas (2008). "Mie theory 1908, on the mobile phone 2008". Journal of Quantitative Spectroscopy Radiative Transfer. 109(8): 1543–1548. Bibcode :2008JQSRT.109.1543W. doi :10.1016/j.jqsrt.2008.01.009.
  • Lorenz, Ludvig (1890). "Lysbevægelsen i og uden for en af plane Lysbølger belyst Kugle". Det Kongelige Danske Videnskabernes Selskabs Skrifter. 6(6): 1–62.

External links

Последняя правка сделана 2021-05-30 10:45:56
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте