Voir la source de Théorie de Mie

En [[optique ondulatoire]], la '''théorie de Mie''', ou '''solution de Mie''', est une solution particulière des [[équations de Maxwell]] décrivant la [[Diffusion élastique de rayonnement|diffusion élastique]] – c'est-à-dire sans changement de longueur d'onde – d'une onde électromagnétique plane par une [[Particule (physique)|particule]] [[sphère|sphérique]] caractérisée par son diamètre et son [[indice de réfraction]] complexe. Elle tire son nom du physicien allemand [[Gustav Mie]], qui la décrivit en détail en [[1908 en science|1908]]. Le travail de son prédécesseur [[Ludvig Lorenz]] est aujourd'hui reconnu comme « empiriquement équivalent »{{sfn|Hergert|Wriedt|2012|loc = Chapitre 3|p = 74}} et l'on parle parfois de la '''théorie de Lorenz-Mie'''. Elle reçut de nombreux apports du physicien [[Peter Debye]] dans les années qui suivirent.

La théorie de Mie s'applique par exemple à des gouttes d'eau de taille macroscopique impliquées dans des [[Photométéore|phénomènes optiques météorologiques]] tels que la formation des [[Arc-en-ciel|arcs-en-ciel]], des [[Couronne (phénomène optique)|couronnes]] ou des [[Gloire (phénomène optique)|gloires]]{{sfn|Hergert|Wriedt|2012|loc=Chapitre 7}}, mais on peut également l'appliquer à des particules microscopiques telles que des [[molécule]]s, en remplaçant l'indice de réfraction par la [[polarisabilité]] de la molécule grâce à l'[[équation de Lorentz et Lorenz]]. Le cas des petites particules correspond à la [[diffusion de Rayleigh]] qui donne les mêmes résultats que la théorie de Mie jusqu'à des rayons de l'ordre de {{Formule|{{Fraction|''λ''|10}}}}. Pour les plus grandes particules, la théorie exacte de Mie s'écarte notablement des résultats de Rayleigh. On parle ainsi de '''diffusion de Mie''' pour décrire les phénomènes qualitativement différents observés dans le cas de sphères de taille comparable ou supérieure à la longueur d'onde {{Formule|''λ''}}.

== Historique ==
Dans son article de 1908{{sfn|Mie1908}}, Mie étudie la diffusion de la lumière par des petites sphères d'or en solution (aujourd'hui on parlerait de nanoparticules), expliquant ainsi pourquoi des sphères de différents diamètres donnent des couleurs différentes. Étant un physicien théoricien, il s'attache cependant à trouver une solution générale à la diffusion de particules sphériques, que celles-ci soient métalliques ou diélectriques, en se plaçant dans le cadre de l'[[électromagnétisme]] précédemment formulé par [[James Clerk Maxwell|Maxwell]]. Il se place ainsi dans l'étude de la diffusion des ondes, qui « a une histoire à laquelle ont contribué certains des plus grands noms de la physique mathématique de la fin du {{s mini-|XIX|e}} et du début du {{s-|XX|e}} »<ref name=Logan1965>{{article| doi = 10.1109/PROC.1965.4055| issn = 0018-9219| volume = 53| numéro = 8| pages = 773–785| nom = Logan| prénom = N.A.| titre = Survey of some early studies of the scattering of plane waves by a sphere| journal = Proceedings of the IEEE| date = 1965-08}}</ref>.

En décrivant la distribution de charge portée par la sphère par un simple [[Dipôle électrostatique|dipôle]], [[Lord Rayleigh]] a pu expliquer simplement à partir de 1871 la diffusion par des particules de taille plus petites que la longueur d'onde. La [[diffusion de Rayleigh]] qui porte son nom peut d'ailleurs être retrouvée en ne retenant que les premiers termes dans le développement multipolaire donné par la théorie de Mie. La théorie de Mie est donc plus générale, mais il est maintenant établi que le problème de la diffusion par une sphère a été résolu plusieurs fois de manière indépendante, et ce avant même le travail de Mie{{sfn|Hergert|Wriedt|2012|p=30}}. En particulier, [[Alfred Clebsch]] l'a résolu dans un mémoire de 1863 dans le cas des {{Lien|langue=en|trad=Elastic wave|fr=Onde élastique|texte=ondes élastiques}}, et [[Ludvig Lorenz]] s'est appuyé sur les résultats de Clebsch pour les appliquer en 1890 à la diffusion des ondes électromagnétiques<ref name=Logan1965 />. Le travail de Lorenz a cependant eu peu de répercussions<ref name=Logan1965 />{{,}}{{sfn|Hergert|Wriedt|2012|p=31}}, notamment parce qu'il ne se plaçait pas dans le même formalisme que Maxwell (bien que sa théorie soit équivalente<ref name=Logan1965 />) et parce que son mémoire a été écrit en [[danois]] (bien qu'une traduction en français ait été disponible à partir de 1896<ref name=Logan1965 />). En 1909, un an après l'article de Mie, [[Peter Debye]] simplifia la dérivation du développement multipolaire du champ en introduisant des fonctions scalaires, maintenant appelées potentiels de Debye, dont dérivent les champs électriques et magnétiques. Mie ne connaissait probablement pas le travail de Lorenz et la raison pour laquelle l'histoire a retenu son nom est que son article fournissait un guide complet pour les études futures, et parce qu'il avait avancé les calculs numériques plus loin que les auteurs précédents<ref name=Logan1965 />.

En effet, pour calculer les coefficients de diffusion, il est nécessaire d'évaluer des [[Fonction spéciale|fonctions spéciales]] telles que les [[fonction de Riccati-Bessel|fonctions de Riccati-Bessel]]<ref name="HMF P445">{{Handbook of Mathematical Functions (Abramowitz et Stegun)|numéro chapitre=10.3|passage=445}}</ref> et leurs dérivées. Avant l'arrivée des ordinateurs, ces calculs représentaient un travail énorme. Aussi, et bien que ses contemporains reconnaissent le travail de Mie, son article est resté en grande partie ignoré jusque dans les années 1950{{sfn|Hergert|Wriedt|2012|loc=Section 2.3|p=55}} qui virent un intérêt grandissant dans l'étude des [[colloïde]]s. Plusieurs [[monographie]]s ayant marqué le domaine datent de cette période de réappropriation de la théorie de Lorenz-Mie, comme celles de van de Hulst en 1957, de Kerker en 1969 ou de Bohren et Huffman en 1983. Même après l'arrivée des ordinateurs, il a fallu développer des codes fiables et efficaces pour le calcul des coefficients de diffusion{{sfn|Hergert|Wriedt|2012|loc=Section 2.3|p=56}}. De tels codes sont disponibles depuis le début des années 1970 et l'on peut maintenant calculer facilement sur un ordinateur personnel la diffusion par des sphères de paramètre de taille supérieurs à {{Formule|10&thinsp;000}}. Des programmes sont disponibles en cliquant sur les liens de la section [[#Liens externes|Liens externes]] en bas de cet article.

== Géométrie de la diffusion ==
[[Fichier:Mie scattering - geometry and notations.svg|vignette|Notations utilisées dans l'article.]]

Dans la diffusion de la lumière résolue par la théorie de Mie, une [[onde plane]] monochromatique de [[vecteur d'onde]] {{Formule|§=''k'' = {{Fraction|2π|''λ''}}}} se propage dans un milieu d'[[indice de réfraction]] {{Formule|''N''}}. Cette onde incidente rencontre une sphère de rayon {{Formule|''a''}} et d'indice de réfraction {{Formule|''N<sub>1</sub>''}}. À la suite de l'interaction une onde est réfléchie, ou plus exactement diffusée par la sphère dans tout l'espace. Pour connaître son amplitude en un point {{Formule|''P''}} situé à l'extérieur de la sphère, on repère celui-ci en coordonnées sphériques {{Formule|(''r'', ''θ'', ''φ'')}}. Comme le montre le schéma ci-contre où {{Formule|§=''φ'' = 0}}, le problème est invariant par rotation selon l'angle {{Formule|''φ''}} tant que l'on ne considère pas la [[Polarisation (optique)|polarisation]] de l'onde. On peut donc choisir soit de partir d'une polarisation fixée et de laisser {{Formule|''φ''}} variable – c'est ce qui est fait dans la dérivation de la solution de Mie – soit considérer que {{Formule|''φ''}} est nul et utiliser une polarisation variable – c'est ce qui est fait dans l'exploitation des résultats.

Le '''plan de diffusion''' est le plan donné par la direction de {{Formule|'''k'''}} et par la direction d'observation {{Formule|'''e<sub>r</sub>'''}}. Il correspond donc au plan du schéma et permet de repérer la polarisation de l'onde. Une polarisation est dite perpendiculaire si elle est orientée selon {{Formule|'''e<sub>y</sub>'''}} et parallèle si elle est dans le plan du schéma. Si l'on connaît les [[Amplitude de diffusion|amplitudes de diffusion]] pour ces deux polarisations, notées respectivement {{Formule|''S''<sub>1</sub>(''θ'')}} et {{Formule|''S''<sub>2</sub>(''θ'')}}, il est possible de prévoir la diffusion pour n'importe quelle polarisation incidente en la décomposant sur ces deux polarisations de base.

Outre la polarisation de l'onde, la diffusion est entièrement déterminée par le '''paramètre de taille''' {{Formule|''x''}} et le rapport d'indice {{Formule|''m''}} définis par :
:<math>x = ka = \frac{2 \pi N a}{\lambda} \quad \text{et} \quad m = \frac{N_1}{N}.</math>

== Section efficace de diffusion ==
[[Fichier:Diffusion de Mie - Efficacité de diffusion en fonction du paramètre de taille.svg|vignette|droite|upright=1.5|Efficacité de diffusion {{Formule|''Q''<sub>diff</sub>}} en fonction du paramètre de taille {{Formule|§=''x'' = ''k·a''}} pour un rapport d'indice {{Formule|§=''m'' = 1.33}} entre les deux matériaux – correspondant à une sphère d'eau dans l'air. L'insert montre la même courbe en échelle logarithmique, mettant en évidence la loi de puissance {{Formule|''Q''<sub>diff</sub> ∝ ''x''<sup>4</sup>}} pour les sphères de petite taille.]]

La sphère a une dimension finie alors que l'onde plane incidente a une extension en théorie infinie. La section efficace de diffusion {{Formule|§=''σ''<sub>diff</sub>}} permet de caractériser l'interaction entre les deux. Elle est définie comme le rapport entre la [[Puissance (physique)|puissance]] diffusée (en [[watt]]s) et l'[[Éclairement énergétique|irradiance]] de l'onde incidente (en W⋅m{{-2}}). Intuitivement, l'ombre portée par une sphère opaque de rayon {{Formule|''a''}} est de surface {{Formule|''π a<sup>2</sup>''}}, aussi on compare souvent la section efficace de diffusion avec cette section droite en définissant '''l'efficacité de diffusion'''
:<math>Q_\text{diff} = \frac{\sigma_\text{diff}}{\pi \, a^2}.</math>

Le graphique de droite montre l'efficacité de diffusion calculée pour un rapport d'indice de 1.33 entre la particule et le milieu. Il permet de montrer plusieurs caractéristiques de la diffusion de Mie. Après une augmentation de {{Formule|''Q''<sub>diff</sub> {{=}} 0}} à {{Formule|''Q''<sub>diff</sub> {{=}} 4}}, l'efficacité de diffusion tend asymptotiquement vers 2 par des oscillations amorties. À ces oscillations de grande période se superposent des pics étroits répartis irrégulièrement. Ces caractéristiques sont présentées plus en détail ci-dessous.

=== Petites sphères – diffusion Rayleigh ===
[[Image:Sky with puffy clouds.JPG|vignette|gauche|Le ciel et les nuages sont très différents… mais leur couleur est expliquée par la même théorie !]]

Pour les petites sphères (région {{Formule|§=''x'' < 1}}), l'efficacité de diffusion augmente très rapidement de manière proportionnelle à {{Formule|''x''<sup>4</sup>}}. On retrouve là les résultats de la [[diffusion Rayleigh]], qui donnent pour cette région{{sfn|Jackson|loc = Section 10.1.2|p = 481}}

:<math>Q_\text{diff} = \frac 8 3 \left|\frac{m^2-1}{m^2+2}\right|^2 x^4.</math>

La diffusion Rayleigh est un cas limite de la diffusion de Mie. Néanmoins, elle en diffère par plusieurs aspects perceptibles lorsqu'on les compare pour des particules de tailles très différentes.On peut apprécier la différence entre la diffusion Rayleigh et la diffusion Mie en observant le ciel : pour les molécules qui constituent l'[[Atmosphère (Terre)|atmosphère]], la première explique la [[Couleur du ciel|couleur bleue du ciel]] ; pour les gouttelettes d'eau qui forment les [[nuage]]s, la seconde explique leur couleur [[blanc]]he.

En effet, pour l'air, la sphère à considérer est la sphère moyenne occupée par une molécule du gaz<ref>On peut le voir en rapprochant les équations 10.11 et 10.34 du livre de {{Harvsp|Jackson}}. Le volume correspondant est le [[volume molaire]] divisé par le [[nombre d'Avogadro]].</ref>, qui est de rayon {{Formule|§=''a'' ≈ 2 nm}}. Le spectre visible se situe donc entre {{Formule|§=''x'' ≈ 1.6 × 10<sup>-2</sup>}} pour le rouge et {{Formule|§=''x'' ≈ 3.2 × 10<sup>-2</sup>}} pour le bleu, donc précisément dans la région de la diffusion Rayleigh. À cause de la puissance 4, le bleu sera 16 fois plus diffusé que le rouge, donnant la couleur du ciel. Pour les nuages, les gouttes d'eau sont beaucoup plus grosses, de l'ordre de {{Formule|§=''a'' ≈ 10 µm}}<ref>[http://www.passion-meteo.net/decouverte/propriete_nuage.htm]</ref>, aussi le spectre visible s'étale entre {{Formule|§=''x'' ≈ 80}} et {{Formule|§=''x'' ≈ 160}}. Comme le montre le graphique, dans cette région l'intensité diffusée sera essentiellement la même pour toutes les longueurs d'onde et le nuage paraîtra blanc, ou gris dans les régions plus épaisses<ref>{{de}} « [http://www.gup.uni-linz.ac.at/thesis/diploma/paul_heinzlreiter/html/node11.html#SECTION00543000000000000000 ''Grundlagen der atmosphärischen Optik''] ».</ref>.

=== Oscillations larges – structure d'interférences ===
[[Fichier:Mie scattering - scattering efficiency as a function of refractive index and size parameter.png|vignette|Efficacité de diffusion en fonction du paramètre de taille ''x'' et du rapport d'indice ''m''.]]
Entre les très petites et les très grandes valeurs de ''x'', l'efficacité de diffusion oscille en fonction de ''x'', avec une période de l'ordre de 10 dans l'exemple précédent de la goutte d'eau. Ces oscillations peuvent également donner lieu à des phénomènes météorologiques surprenants, comme l'observation inhabituelle d'un soleil bleu dans le nord de l'[[Europe]] le {{date-|26 septembre 1950}} et les jours suivants<ref name="BlueSun1">{{article
 | issn = 0035-8711
 | volume = 111
 | pages = 478
 | nom = Wilson| prénom = R.
 | titre = The blue Sun of 1950 September
 | journal = Monthly Notices of the Royal Astronomical Society
 | consulté le = 2015-07-24
 | date = 1951
 | résumé = http://adsabs.harvard.edu/abs/1951MNRAS.111..478W
 | lire en ligne = http://adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1951MNRAS.111..478W&link_type=ARTICLE&db_key=AST&high=}}</ref>{{,}}<ref>{{article| doi = 10.1111/j.2153-3490.1951.tb00790.x
 | issn = 2153-3490
 | volume = 3| numéro = 3
 | pages = 135–140
 | nom = Ångström| prénom = Anders
 | titre = The Blue Sun of September 1950
 | journal = Tellus
 | consulté le = 2015-07-24
 | date = 1951-08-01| résumé = http://onlinelibrary.wiley.com/doi/10.1111/j.2153-3490.1951.tb00790.x/abstract}}</ref>. L'explication de ce phénomène a pu être trouvée : à la suite de grands feux ayant eu lieu en [[Alberta]] les jours précédents, des gouttes d'huile micrométrique ont été poussées dans la haute atmosphère jusqu'en Europe. Le spectre réalisé en [[Écosse]]<ref name="BlueSun1" /> s'accorde bien avec la théorie de Lorenz-Mie si la taille des particules permet au premier minimum d'oscillation de se trouver à {{Formule|§=''λ'' = 435 nm}}, dans le bleu du spectre. Au lieu d'être 16 fois plus diffusé comme dans la diffusion Rayleigh, le bleu est alors 2 fois ''moins'' diffusé et le soleil apparaît bleu.

L'image ci-contre montre comment évolue l'efficacité de diffraction, non seulement avec le paramètre de taille, mais aussi avec le rapport d'indice ''m'' entre la sphère et le milieu environnant. Un [[développement limité]] à partir de la théorie de Mie{{sfn|Bohren|Huffman|1983|loc=Section 4.4.2|p=105}}, valable pour les grands paramètres de taille, montre que les maxima ont lieu approximativement pour
:<math>2x\cdot (m-1) = (2p+1)\pi,\quad p \in \N.</math>
Les courbes correspondantes sont tracées en noir sur la figure. On voit que cette approximation donne la bonne période, mais qu'il y a un léger décalage dû aux valeurs modérées de ''x'' sur cette figure. L'interprétation physique simple est que ces oscillations correspondent aux interférences entre l'onde ayant traversé la sphère et celle étant passée à côté{{sfn|Bohren|Huffman|1983|loc=Section 4.4.2|p=106}}. En effet, le déphasage entre les deux est justement donné par {{Formule|2''x''(''m'' – 1)}}.

=== Grosses sphères – paradoxe de l'extinction ===
Dans la limite des grosses sphères, il peut paraître surprenant que la section efficace d'extinction {{Formule|§=''σ''<sub>ext</sub> = ''σ''<sub>diff</sub> + ''σ''<sub>abs</sub>}} ne tende pas simplement vers la section droite {{Formule|πa²}} de la sphère, mais vers le double de cette valeur. Intuitivement, on ne s'attend en effet pas à ce qu'un gros objet, pour lequel l'[[optique géométrique]] est valable, ne retire de l'onde 2 fois la puissance qu'il reçoit. Ce résultat, appelé '''paradoxe de l'extinction'''{{sfn|Bohren|Huffman|1983|loc=Section 4.4.3|p=107}} est général et valable que l'indice de réfraction de la sphère soit réel – auquel cas {{Formule|§=''σ''<sub>abs</sub> = 0}} – ou complexe.

Il s'explique facilement en considérant la [[diffraction de Fraunhofer]] sur un disque opaque de même diamètre que la sphère. En vertu du [[principe de Babinet]], la figure de diffraction aura la forme d'une [[tache d'Airy]]. En intégrant sur tout l'espace, on retrouve bien le résultat que{{sfn|Bohren|Huffman|1983|loc=Section 4.4.3|p=110}}

:<math>\sigma_\text{abs}\to\pi a^2\;\text{et}\;\sigma_\text{diff}\to\pi a^2, \quad \text{donc}\; \sigma_\text{ext}\to 2\pi a^2.</math>

Cependant, la puissance diffusée ne s'écartera que très peu de l'axe du faisceau. L'essentiel de la puissance se trouvera dans la première tache d'Airy, correspondant à un angle de diffraction de {{Formule|1.22π/''x''}}. Pour distinguer l'onde diffusée de l'onde transmise, il faut avoir une résolution angulaire d'autant meilleure que ''x'' est important. Ainsi, à partir de {{Formule|§=''x'' = 440}}, cet angle devient inférieur à 0.5° qui est l'angle apparent du soleil.

=== Pics fins – modes de galerie ===
[[Fichier:Whispering gallery modes sphere.png|vignette|Modes de galerie observés par fluorescence sur une sphère de [[Dioxyde de silicium|silice]] de {{unité|300|µm}} de diamètre. Un laser rouge permettant d'exciter ces modes est amené par la [[fibre optique]] visible sur la droite<ref>[http://www.techbriefs.com/component/content/article/3160 A. Matsko et al., NASA Tech Briefs NPO-44956, Sept. 1, 2008]</ref>.]]

Autour de certaines fréquences, la section efficace de diffusion a un comportement [[résonance|résonant]] ; on observe des pics qui peuvent être très fins et qui sont dus à des {{Lien|langue=en|trad=Whispering gallery wave|fr=mode de galerie|texte=modes de galerie}} dans la sphère. Ils ne sont en général pas observable sur des ensembles de sphères en raison de leur très faible largeur spectrale : la moindre dispersion en taille de la population de sphères élargira et brouillera les résonances. Ils peuvent cependant être utilisés en [[photonique]] pour créer des filtres très sélectifs ou des lignes à délai. La [[pression de radiation]] due à la diffusion au voisinage d'une résonance est également utile pour créer des [[Pince optique|pinces optiques]] à partir de sphères diélectriques{{sfn|Bohren|Huffman|1983|loc=Section 11.4}}.

=== Impact de l'absorption ===
Si la sphère n'est pas un diélectrique sans perte, on peut observer quelques changements par rapport à ce qui est décrit ci-dessus{{sfn|Bohren|Huffman|1983|loc=Section 4.4.2|p=106}}. Tout d'abord, les pics fins s'élargissent puis disparaissent rapidement car ils sont dus à des ondes qui font un grand nombre de tours dans la sphère. Pour de plus forts coefficients d'absorption, la structure d'interférence disparaît également. Il est également possible que de nouveaux pics apparaissent, liés à une résonance dans le coefficient d'absorption du matériau.

{{clr}}

== Distribution angulaire du rayonnement diffusé ==
Le cas de la diffusion par des très petites particules, telles que des [[molécule]]s, de dimensions inférieures au dixième de la longueur d'onde considérée, est un cas limite appelé [[diffusion Rayleigh]]. Pour les particules plus grosses que cette longueur d'onde, on doit prendre en compte la diffusion de Mie dans son intégralité : elle explique dans quelles directions la diffusion est la plus intense, on obtient ainsi un « patron de réémission » qui ressemble à celui des lobes d'émission d'une antenne, avec, dans le cas de grosses particules, un lobe plus intense dans le sens de l'onde incidente.

[[Image:Mie scattering.svg|centré|cadre|De gauche à droite : intensité de la diffusion Rayleigh, de la diffusion Mie pour de petites particules et de la diffusion Mie pour de grosses particules, en fonction de la direction. L'onde incidente arrive par la gauche.]]

La diffusion de Mie n'est pas toujours fortement dépendante de la longueur d'onde utilisée comme c'est le cas dans celle de Rayleigh. Elle produit donc une lumière presque [[blanc]]he lorsque le [[Soleil]] illumine de grosses particules dans l'air : c'est cette dispersion qui donne la couleur blanc laiteux à la [[brume]] et au [[brouillard]].
{{Voir aussi|Réflexion diffuse|Source lumineuse orthotrope}}

Cependant, si les solutions fournies par la diffusion de Mie sont exactes (pour des sphères), elles se présentent sous la forme d'une série infinie de fonctions spéciales, et on est souvent limité à des approches numériques.

== Résolution mathématique ==
=== Position du problème ===
[[Image:3D Miestreuung an 2um Kugel.jpg|vignette|droite|250px|Représentation de la diffusion Mie pour une particule sphérique de 2 µm de rayon, éclairée par la gauche avec de la lumière rouge (λ = 633 nm).]]

Le problème de la diffusion par une sphère est un problème d'[[électromagnétisme]] classique. On considère une [[onde plane]] monochromatique de [[Vitesse angulaire|pulsation]] ''ω'' et de [[Polarisation (optique)|polarisation]] connue se propageant dans un milieu linéaire, homogène et isotrope d'[[indice de réfraction]] complexe ''N'' = ''n'' + i ''κ''. Cette onde incidente rencontre une sphère de rayon ''a'' et d'indice de réfraction complexe ''N<sub>1</sub>'' = ''m'' x ''N'', où ''m'' est un nombre complexe caractérisant le contraste entre la sphère et le milieu environnant. Comme cela est souligné par Bohren{{sfn|Bohren|Huffman|1983|loc=Preface|p=vii}}, cette manière de décrire les milieux cache l'origine microscopique, matérielle de l'indice de réfraction pour présenter un traitement mathématique rigoureux. Dans certaines situations, il convient donc d'être prudent, par exemple lorsque le matériau est anisotrope ou lorsqu'il y a des effets optiques non-linéaires à l'intérieur de la particule ([[diffusion Raman]] et [[diffusion Brillouin|Brillouin]], [[fluorescence]] de la particule, [[Absorption saturée|saturation]] des transitions…).

Pour résoudre le problème, on commence par trouver une base d'ondes électromagnétiques de pulsation ''ω'' adaptée aux [[coordonnées sphériques]]. Elles sont généralement appelées {{Lien|langue=en|trad=Vector spherical harmonics|fr=harmoniques sphériques vectorielles}} et peuvent être dérivées des [[Harmonique sphérique|harmoniques sphériques]] scalaires habituelles de différentes manières. Comme il s'agit d'une [[base de Hilbert]], on peut alors décomposer tout champ de pulsation ω comme une somme infinie d'harmoniques sphériques vectorielles : c'est ce que l'on appelle la {{Lien|langue=en|trad=Electromagnetic_wave_equation#Multipole_expansion|fr=Electromagnetic_wave_equation#Multipole_expansion|texte=décomposition multipolaire du champ électromagnétique}}. En particulier l'onde plane incidente se décompose de la sorte. On écrit alors que le champ électromagnétique total est composé de trois parties : le champ incident ('''E<sub>i</sub>''', '''H<sub>i</sub>''') et le champ diffusé ('''E<sub>d</sub>''', '''H<sub>d</sub>''') se superposent à l'extérieur de la sphère, et le champ interne ('''E<sub>l</sub>''', '''H<sub>l</sub>''') règne à l'intérieur. En exprimant les [[condition aux limites|conditions aux limites]] à la surface de la sphère, on peut alors trouver les coefficients du développement multipolaire des champs. Cette démarche est similaire à celle suivie pour trouver les coefficients de transmission et de réflexion d'un dioptre plan, la complication mathématique venant de la géométrie sphérique à laquelle l'onde incidente n'est pas adaptée. Une fois que les coefficients du développement sont connus, on peut les utiliser pour exprimer d'autres quantités, comme la puissance diffusée dans une direction particulière ou encore la polarisation de la lumière diffusée. Les coefficients trouvés dépendent cependant du choix de la base d'harmoniques sphériques vectorielles utilisée, or différents auteurs utilisent des fonctions légèrement différentes (elles sont notamment normalisées de manière différente). C'est pourquoi les formules peuvent sembler différentes alors que ce n'est que le choix de la base qui change.

=== Développement multipolaire du champ électromagnétique ===
Le développement multipolaire d'un champ de pulsation fixée consiste à écrire ce champ non pas sur une base d'onde planes, mais sur des harmoniques sphériques vectorielles adaptées aux coordonnées sphériques. Suivant la dérivation faite par Bohren et Huffman{{sfn|Bohren|Huffman|1983|loc=Section 4.1}}, on utilise les harmoniques sphériques vectorielles <math>\mathbf M^\pm_{\ell,m}(r, \theta, \phi)</math> et <math>\mathbf N^\pm_{\ell,m}(r, \theta, \phi)</math>. On obtient alors

:<math>\begin{align}\mathbf E(r, \theta, \phi) = \sum_{\ell=0}^\infty \sum_{m=0}^\ell  & B^+_{\ell,m} \mathbf M^+_{\ell,m}(r, \theta, \phi) + B^-_{\ell,m} \mathbf M^-_{\ell,m}(r, \theta, \phi) + \ldots \\
& + A^+_{\ell,m} \mathbf N^+_{\ell,m}(r, \theta, \phi) + A^-_{\ell,m} \mathbf N^-_{\ell,m}(r, \theta, \phi).\end{align}</math>

Les modes {{Formule|'''N'''<sup>±</sup>}} apparaissant dans le développement de {{Formule|'''E'''}} sont appelés les modes électriques car ils correspondent au champ rayonné par un [[Dipôle électrostatique|dipôle]], [[quadrupôle|quadripôle]]…, multipôle électrique. Par exemple, les trois termes {{Formule|'''N'''{{indexp|1,0|+}}}}, {{Formule|'''N'''{{indexp|1,1|+}}}} et {{Formule|'''N'''{{indexp|1,1|–}}}} correspondent au [[rayonnement dipolaire électrique]] – le quatrième terme {{Formule|'''N'''{{indexp|1,0|–}}}} étant nul dans tout l'espace. Ces mêmes modes sont aussi appelés modes transverses magnétiques ou TM car le champ magnétique correspondant est en tout point perpendiculaire à {{Formule|'''e<sub>r</sub>'''}}. De la même manière, les modes où {{Formule|'''E'''}} est proportionnel à {{Formule|'''M<sup>±</sup>'''}} sont appelés les modes magnétiques ou modes transverses électriques.

{{démonstration|titre=Détails sur le développement multipolaire
|Concrètement, on peut montrer que si le champ a une pulsation {{Formule|''ω''}} unique, les [[équations de Maxwell]] se ramènent à deux [[équation de Helmholtz|équations de Helmholtz]] vectorielles indépendantes pour le champ {{Formule|'''E'''}} et pour le champ {{Formule|'''H'''}}{{sfn|Jackson|loc=Section 9.7|p=450}}{{,}}{{sfn|Bohren|Huffman|1983|loc=Section 4.1|p=83}}. Si l'on privilégie le champ {{Formule|'''E'''}}, on obtient les trois équations
:<math>\begin{cases}(\nabla^2+k^2)\mathbf E=0 , \\ \vec{\nabla}\cdot \mathbf E=0 , \\ \mathbf H = -\frac{i}{\omega \mu}\vec{\nabla}\times \mathbf E .\end{cases}</math>

La première est l'équation de Helmholtz et fait intervenir la norme du vecteur d'onde dans les deux milieux :
:<math>k^2=\omega^2\epsilon \mu = \begin{cases}\frac{\omega^2}{c^2}N^2 & \text{dans le milieu extérieur,} \\\frac{\omega^2}{c^2}m^2N^2 & \text{à l'intérieur de la sphère ;}\end{cases}</math>
la seconde exprime le fait que le champ doit être sans divergence et la troisième permet de déduire le champ {{Formule|'''H'''}} une fois les deux premières équations résolues.

Si l'on privilégie le champ {{Formule|'''H'''}}, on obtient un système de trois équations très similaire :
:<math>\begin{cases}(\nabla^2+k^2)\mathbf H=0 , \\\vec{\nabla}\cdot \mathbf H=0 , \\ \mathbf E = \frac{i}{\omega \epsilon}\vec{\nabla}\times \mathbf H .\end{cases}</math>

Vu la similitude entre les deux systèmes d'équations, toute solution pour le champ {{Formule|'''E'''}} dans le premier système fournira une solution valide pour le champ {{Formule|'''H'''}} dans le second si l'on échange le rôle de {{Formule|'''E'''}} et de {{Formule|'''H'''}} – en changeant l'amplitude d'un facteur {{Formule|''ε''/''μ''}}. Il existe plusieurs méthodes pour trouver ces solutions<ref>{{Harvsp|Jackson|loc=Section 9.7|p=450}} cite l'article suivant : {{Article| doi = 10.1016/S0031-8914(54)80068-1| issn = 0031-8914| volume = 20| numéro = 1–6| pages = 539–554| nom1 = Bouwkamp| prénom1 = C. J.| nom2 = Casimir| prénom2 = H. B. G.| titre = On multipole expansions in the theory of electromagnetic radiation| journal = Physica| consulté le = 2015-07-15| date = 1954| url = http://www.sciencedirect.com/science/article/pii/S0031891454800681}}</ref>. On peut en particulier les générer à partir de fonctions {{Formule|''Ψ''}} qui satisfassent l'[[équation de Helmholtz]] ''scalaire''{{sfn|Bohren|Huffman|1983|loc=Section 4.1|p=84}} :

:<math>(\nabla^2 + k^2) \Psi = 0,</math>

si l'on définit les champs {{Formule|'''M'''}} et {{Formule|'''N'''}} tels que :

:<math>\mathbf M = \nabla \times (\mathbf r \Psi),</math>
:<math>\mathbf N = \frac{\nabla \times \mathbf M}{k},</math>

avec {{Formule|'''r'''}} le vecteur position. Alors, en [[coordonnées sphériques]], {{Formule|'''M'''}} et {{Formule|'''N'''}} représentent deux solutions possibles pour le champ {{Formule|'''E'''}} – les champs {{Formule|'''H'''}} correspondants étant proportionnels à {{Formule|'''N'''}} et {{Formule|'''M'''}} respectivement. Par [[séparation des variables]], on trouve par exemple les deux fonctions ci-dessous :{{sfn|Bohren|Huffman|1983|loc=Section 4.1|p=87}}{{,}}<ref group="note">La notation {{Formule|''Ψ<sup>+</sup><sub>ℓ,m</sub>''}} permet de se rapprocher de celle des harmoniques sphériques et correspond à {{Formule|''Ψ<sub>emn</sub>''}} dans les notations de {{harvsp|Bohren|Huffman|1983}} (+ correspond à ''e'' et ''ℓ'' correspond à ''n'').</ref>
:<math>\Psi^+_{\ell,m} = \cos (m \phi) P_\ell^m (\cos \theta) z_\ell(k r),</math>
:<math>\Psi^-_{\ell,m} = \sin (m \phi) P_\ell^m (\cos \theta) z_\ell(k r),</math>
où ''ℓ'' et ''m'' sont deux [[entier relatif|paramètres entiers]] tels que {{Formule|0 ≤ ''m'' ≤ ''ℓ''}} pour la fonction {{Formule|''Ψ<sup>+</sup>''}} et {{Formule|0 < ''m'' ≤ ''ℓ''}} pour la fonction {{Formule|''Ψ<sup>−</sup>''}}<ref group="note">En général, le paramètre ''m'' peut être négatif et varier de {{Formule|–''ℓ''}} à {{Formule|''ℓ''}}, mais {{harvsp|Bohren|Huffman|1983}} a décidé d'utiliser les [[Harmonique sphérique#Forme « réelle » des harmoniques sphériques|harmoniques sphériques réelles]] <math>\tilde{Y}_{\ell m}</math> sans inclure la constante de normalisation et de séparer explicitement les cas {{Formule|''m'' ≥ 0}} (fonction {{Formule|''Ψ<sup>+</sup>''}}) et {{Formule|''m'' < 0}} (fonction {{Formule|''Ψ<sup>−</sup>''}}).</ref>, {{Formule|''P<sub>ℓ</sub><sup>m</sup>''}} les [[polynômes associés de Legendre|fonctions associées de Legendre]], et {{Formule|''z<sub>ℓ</sub>''}} une combinaison linéaire de deux des quatre [[fonction de Bessel sphérique|fonctions de Bessel sphériques]] : {{Formule|''j<sub>ℓ</sub>''}}, {{Formule|''y<sub>ℓ</sub>''}}, {{Formule|''h''{{indexp|''ℓ''|(1)}}}} et {{Formule|''h''{{indexp|''ℓ''|(2)}}}}<ref group="note">Il faudrait en toute rigueur multiplier par deux le nombre de fonctions de la base. Cependant, souvent une seule des fonctions de Bessel a un coefficient non nul. Les coefficients de la combinaison linéaire dépendent des conditions aux limites, en particulier en {{Formule|''r'' {{=}} 0}} et à l'infini. La fonction {{Formule|''j<sub>ℓ</sub>''}} présente l'avantage de rester réelle et finie en {{Formule|''r'' {{=}} 0}}. La fonction {{Formule|''h''{{indexp|''ℓ''|(1)}}}} quant à elle correspond, pour {{Formule|''r''}} grand, à une onde sphérique sortante : <math>h^{(1)}_\ell(kr) \equiv (-i)^{\ell+1}e^{ikr}/(kr)</math>.</ref>.

À chaque fonction {{Formule|''Ψ''<sup>±</sup>}} correspond un champ {{Formule|'''M'''}} et un champ {{Formule|'''N'''}}, soit quatre au total, formant une [[base de Hilbert]] sur laquelle on peut décomposer toutes les solutions aux équations de Maxwell :{{sfn|Bohren|Huffman|1983|loc=Section 4.1|p=90}}
:<math>\mathbf E (r, \theta, \phi) = \sum_{\ell=0}^\infty \sum_{m=0}^\ell  B^+_{\ell,m} \mathbf M^+_{\ell,m} + B^-_{\ell,m} \mathbf M^-_{\ell,m} + A^+_{\ell,m} \mathbf N^+_{\ell,m} + A^-_{\ell,m} \mathbf N^-_{\ell,m}</math>}}

On suppose que le champ électrique incident est une onde plane se propageant selon {{Formule|''z''}} (axe {{Formule|''θ'' {{=}} 0}}) et polarisée linéairement selon {{Formule|''x''}} (axe {{Formule|''φ'' {{=}} 0}}, {{Formule|''θ'' {{=}} π/2}}). Son développement multipolaire ne contient que des termes de moment {{Formule|''m'' {{=}} 1}} :{{sfn|Bohren|Huffman|1983|loc=Section 4.2|p=92}}

:<math>\mathbf {E_i} = E_0 e^{ikz} \mathbf {e_x} = E_0 \sum_{\ell=1}^\infty u_\ell(\mathbf M^-_{\ell,1}-i\mathbf N^+_{\ell,1}),\quad\text{avec}\quad u_\ell = i^\ell\frac{2\ell+1}{\ell(\ell+1)}.</math>

De la même manière, on peut décomposer le champ ('''E<sub>l</sub>''', '''H<sub>l</sub>''') à l'intérieur de la sphère. Il ne contient que les termes déjà inclus dans '''E<sub>i</sub>'''<ref>{{harvsp|Bohren|Huffman|1983|loc=Section 4.3|p=93}}</ref>

:<math>\mathbf {E_l} = E_0 \sum_{\ell=1}^\infty u_\ell(c_\ell \mathbf M^-_{\ell,1}-i d_\ell \mathbf N^+_{\ell,1}),</math>

De même pour le champ diffusé ('''E<sub>d</sub>''', '''H<sub>d</sub>''')<ref group="note">Pour être complet, il faut préciser quelle(s) fonction de Bessel sphérique est utilisée dans les expressions de '''M''' et '''N'''. Pour les champs incident et à l'intérieur de la sphère, il s'agit de la fonction ''j<sub>ℓ</sub>'' qui est finie en ''r'' = 0, point auquel ces champs existent réellement, et pour le champ diffusé il s'agit de la fonction ''h''{{indexp|''ℓ''|(1)}} qui forme une onde sphérique sortante à grande distance.</ref> :

:<math>\mathbf {E_d} = E_0 \sum_{\ell=1}^\infty u_\ell(a_\ell \mathbf M^-_{\ell,1}-i b_\ell \mathbf N^+_{\ell,1}),</math>

{{démonstration|titre=Champs '''H<sub>i</sub>''', '''H<sub>l</sub>''' et '''H<sub>d</sub>''' correspondants
|Les expressions pour les champs '''H<sub>i</sub>''', '''H<sub>l</sub>''' et '''H<sub>d</sub>''' correspondants sont les suivantes :
:<math>\mathbf {H_i} = -\frac{k}{\omega \mu}E_0 \sum_{\ell=1}^\infty u_\ell(\mathbf M^+_{\ell,1}+i\mathbf N^-_{\ell,1}),</math>
:<math>\mathbf {H_l} = -\frac{k_1}{\omega \mu_1}E_0 \sum_{\ell=1}^\infty u_\ell(d_\ell \mathbf M^+_{\ell,1}+i c_\ell \mathbf N^-_{\ell,1}),</math>
:<math>\mathbf {H_d} = -\frac{k}{\omega \mu}E_0 \sum_{\ell=1}^\infty u_\ell(b_\ell \mathbf M^+_{\ell,1}+i a_\ell \mathbf N^-_{\ell,1}).</math>

Dans l'expression de '''H<sub>l</sub>''', ''k<sub>1</sub>'' {{=}} ''m k'' et ''μ''<sub>1</sub> représentent le vecteur d'onde et la perméabilité magnétique à l'intérieur de la sphère.}}

=== Conditions limites et coefficients de diffusion ===
La valeur des coefficients {{Formule|''a<sub>ℓ</sub>''}}, {{Formule|''b<sub>ℓ</sub>''}}, {{Formule|''c<sub>ℓ</sub>''}} et {{Formule|''d<sub>ℓ</sub>''}} se trouve à partir des conditions aux limites entre le milieu extérieur et la sphère, cette dernière étant supposée ne pas avoir de charge surfacique :
:<math>(\mathbf {E_i} + \mathbf {E_d} - \mathbf {E_l}) \times \mathbf {e_r} = 0</math>
:<math>(\mathbf {H_i} + \mathbf {H_d} - \mathbf {H_l}) \times \mathbf {e_r} = 0</math>

La particule étant de rayon {{Formule|''a''}}, la lumière incidente de [[longueur d'onde]] ''λ'', on introduit pour simplifier le '''paramètre de taille''' {{Formule|''x''}} défini par :
:<math>x = ka = \frac{2 \pi N a}{\lambda}.</math>

On peut alors résoudre le système d'équations donné par les conditions aux limites pour en extraire les coefficients de diffusion{{sfn|Bohren|Huffman|1983|loc=Section 4.4|p=100}}. Si la perméabilité magnétique de la sphère est la même que celle du milieu extérieur, les coefficients de diffusion ne dépendent que du paramètre de taille {{Formule|''x''}} et du contraste d'indice {{Formule|§=''m'' = {{fraction|''N''<sub>1</sub>|''N''}}}}

:<math> a_\ell = \frac{m\, \psi_\ell(m x)\, \psi'_\ell (x) - \psi_\ell(x)\, \psi'_\ell(m x)}{m\, \psi_\ell(m x)\, \xi'_\ell (x) - \xi_\ell(x)\, \psi'_\ell(m x)},</math>

:<math>b_\ell = \frac{\psi_\ell(m x)\, \psi'_\ell (x) - m\, \psi_\ell(x)\, \psi'_\ell(m x)}{\psi_\ell(m x)\, \xi'_\ell (x) - m\, \xi_\ell(x)\, \psi'_\ell(m x)},</math>

où le prime désigne la dérivation et {{Formule|''ψ''}} et {{Formule|''ξ''}} sont les [[fonction de Riccati-Bessel|fonctions de Riccati-Bessel]]<ref name="HMF P445" /> : {{Formule|§=''ψ<sub>ℓ</sub>''(''ρ'') = ''ρ j<sub>ℓ</sub>''(''ρ'')}} et {{Formule|§=''ξ<sub>ℓ</sub>''(''ρ'') = ''ρ h''{{indexp|ℓ|(1)}}(''ρ'')}}.

=== Grandeurs déduites des coefficients ===
À partir des coefficients, on peut extraire plusieurs grandeurs physiques intéressantes. Pour une onde incidente d'amplitude {{Formule|'''E<sub>0</sub>'''}} donnée correspondant à une [[Éclairement énergétique|irradiance]] {{Formule|''I<sub>i</sub>''}}, on peut définir la puissance absorbée {{Formule|''P''<sub>abs</sub>}} comme le flux net du [[vecteur de Pointing]] entrant à travers une sphère <math>\Omega</math> de rayon {{Formule|''r'' > ''a''}} :

:<math>P_\text{abs} = -\oint_\Omega \mathbf {\Pi \cdot e_r} \,\mathrm dS,</math>

tandis que la puissance diffusée {{Formule|''P''<sub>diff</sub>}} est le flux ''sortant'' du vecteur de Pointing ''du seul champ diffusé'' :

:<math>P_\text{diff} = \oint_\Omega \mathbf {\Pi_d \cdot e_r} \,\mathrm dS.</math>

Ces puissances étant proportionnelles à l'irradiance incidente, on caractérise plutôt le rapport <math>P_\text{diff}/I_i</math> qui a la dimension d'une surface appelée [[section efficace]] de diffusion :{{sfn|Bohren|Huffman|1983|loc=Section 4.4.1|p=103}}

:<math>\sigma_\text{diff} = \frac{P_\text{diff}}{I_i} = \frac{2 \pi}{k^2} \sum_{\ell=1}^{\infty} (2\ell+1) \left( |a_\ell|^2 + |b_\ell|^2 \right).</math>

La puissance {{Formule|''P''<sub>diff</sub>}} reste sous forme électromagnétique mais elle part dans toutes les directions et de même que la puissance absorbée, elle apparaît comme soustraite à l'intensité incidente. On définit donc la section efficace d'extinction comme la puissance totale soustraite à l'onde plane incidente pour une irradiance unitaire :

:<math>\sigma_\text{ext} = \frac{P_\text{diff}+P_\text{abs}}{I_i} = \frac{2 \pi}{k^2} \sum_{\ell=1}^{\infty} (2\ell+1) \, \Re(a_\ell + b_\ell)</math>

où <math>\Re</math> désigne la partie réelle. Ces formules peuvent être écrites différemment par d'autres auteurs. Par exemple dans le livre de Jackson{{sfn|Jackson|loc=Section 10.4|p=498}}, on a

:<math>\sigma_\text{ext} = -\frac{\pi}{k^2} \sum_{\ell=1}^{\infty} (2\ell+1) \, \Re(\alpha_\ell + \beta_\ell),</math>
la différence provenant de la définition des harmoniques sphériques vectorielles : si celles-ci sont multipliées par deux, les coefficients doivent être divisés par deux pour retrouver les mêmes résultats.

== Notes et références ==
=== Notes ===
{{Références|groupe=note}}

=== Références ===
{{Références|colonnes=2}}

== Voir aussi ==
=== Articles connexes ===
* [[Diffusion Rayleigh]] / [[Effet Tyndall]] ;
* [[Diffusion Compton]]
* [[Optique ondulatoire]].

=== Liens externes ===
* {{en}} [http://www.scattport.org/diogenes.iwt.uni-bremen.de/vt/scattport/index.html ScattPort] est un portail d'information sur la diffusion de la lumière à l'intention de la communauté scientifique étudiant ce phénomène. Il regroupe en particulier les différents [http://www.scattport.org/diogenes.iwt.uni-bremen.de/vt/scattport/index.php/light-scattering-software.html programmes de calcul] de la diffusion de la lumière et les conférences dans le domaine.
* {{en}} [http://bernstein.harvard.edu/research/nearfield/fdtd/FDTD%20SERS.html Animations de la diffusion Mie] pour différentes particules ;
* {{en}} [http://omlc.ogi.edu/calc/mie_calc.html Calculateur en ligne de la diffusion Mie] ;
* {{en}} [http://www.philiplaven.com/mieplot.htm MiePlot], un programme simple (marche sans installation) permettant de nombreuses visualisations ;
* {{en}} [http://www.T-Matrix.de Implémentations de la diffusion Mie] en [[C++|C, C++]], [[Fortran]], [[Matlab]] et [[Mathematica]] ;
* [http://www.cilas.com/principe-de-granulometrie.html Une autre illustration de la théorie de Mie]

=== Bibliographie ===
* {{Article
 | lang = de
 | auteur = Gustav Mie
 | titre =''Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen''
 | périodique = ''Annalen der Physik, Leipzig''
 | volume = 25
 | année = 1908
 | pages = 377–445
 | lire en ligne = http://www.scattport.org/diogenes.iwt.uni-bremen.de/vt/laser/papers/Mie-Beitrag-zur-Optik-AnnalenPhysik1908.pdf
 | consulté le = 30 juin 2015 
}} ; {{en}} [http://www.scattport.org/diogenes.iwt.uni-bremen.de/vt/laser/papers/SAND78-6018-Mie-1908-translation.pdf traduction anglaise de 1978 par le Laboratoire Sandia] ; [http://www.scattport.org/diogenes.iwt.uni-bremen.de/vt/laser/papers/RAE-LT1873-1976-Mie-1908-translation.pdf traduction anglaise de 1976 par la NASA]  

* {{ouvrage
 | lang = en
 | auteur1= Wolfram Hergert|auteur2=Thomas Wriedt |responsabilité2=eds
 | titre = The Mie Theory
 | éditeur = Springer Berlin Heidelberg
 | lieu = Berlin, Heidelberg
 | année = 2012
 | series = Springer Series in Optical Sciences
 | url = https://link.springer.com/10.1007/978-3-642-28738-1
 | isbn = 978-3-642-28737-4
 | consulté le = 25 juin 2015
}}
* {{fr}} {{Jackson}}
* {{en}} A. Stratton, ''Electromagnetic Theory''. McGraw-Hill, New York, 1941 ;
* {{en}} H. C. van de Hulst, ''Light scattering by small particles''. Dover, New York, 1981, fac-similé de l'édition de John Wiley & Sons, New York, 1957;
* {{en}} M. Kerker, ''The scattering of light and other electromagnetic radiation''. Academic, New York, 1969 ;
* {{Ouvrage
 | lang = en
 | éditeur = Wiley
 | isbn = 978-0-471-29340-8
 | nom1 = Bohren
 | prénom1 = Craig F
 | nom2 = Huffman
 | prénom2 = Donald R
 | titre = Absorption and scattering of light by small particles
 | lieu = New York
 | année = 1983}}
* {{en}} P. W. Barber, S. S. Hill, ''Light scattering by particles: Computational Methods''. World Scientific, Singapour, 1990 ; 
* {{en}} Hong Du, ''« Mie-scattering calculation »'', ''Applied Optics'' '''43''' (9), 1951-1956 (2004).
* {{en}} Thomas Wriedt, ''« Mie theory 1908, on the mobile phone 2008 »'', ''[[J. Quant. Spectrosc. Radiat. Transf.]] '' ''' 109 ''', 1543–1548 (2008).

{{Palette|Diffusion des ondes}}
{{portail|physique|optique}}

{{DEFAULTSORT:Théorie de Mie}}
[[Catégorie:Physique quantique]]
[[Catégorie:Diffusion des ondes|Mie]]
[[Catégorie:Loi en optique|Mie]]
[[Catégorie:Propagation des ondes radioélectriques]]