Résonance ferromagnétique

Modèle:Leçon du jour Cette leçon nécessite une bonne connaissance préalable de la leçon : Introduction à l'électromagnétisme des milieux matériels

Préliminaires

Précession d'un dipôle magnétique

On sait qu'un dipôle magnétique de moment $\overrightarrow{𝔪}$ plongé dans un champ magnétique $\overrightarrow{B}=B{\overrightarrow{u}}_z$ est soumis :

• à la force $\overrightarrow{F}=\overrightarrow{\mathrm{\nabla }}(\overrightarrow{𝔪}\cdot \overrightarrow{B})$
• au moment $\overrightarrow{\mathrm{\Gamma }}=\overrightarrow{𝔪}\wedge \overrightarrow{B}$

On applique le théorème du moment cinétique au dipôle :

$\frac{\mathrm{d}\overrightarrow{L}}{\mathrm{d}t}=\overrightarrow{\mathrm{\Gamma }}$, soit $\frac{\mathrm{d}\overrightarrow{𝔪}}{\mathrm{d}t}=\gamma \overrightarrow{𝔪}\wedge \overrightarrow{B}$

Modèle:Propriété

Comportement de l'aimantation

On peut relier l'aimantation d'un milieu magnétique en un point à la densité du moment magnétique moyen dans un volume mésoscopique ΔV autour du point :

$\overrightarrow{M}=n_v⟨\overrightarrow{𝔪}⟩$

• n_v est la densité volumique de dipôles magnétiques dans ΔV
• $⟨\overrightarrow{𝔪}⟩$ est le moment magnétique moyen des dipôles de ΔV

On peut donc écrire que $\overrightarrow{M}$ suit la même loi d'évolution que $\overrightarrow{𝔪}$ :

$\frac{\mathrm{d}\overrightarrow{M}}{\mathrm{d}t}=\gamma \overrightarrow{M}\wedge \overrightarrow{B}$

Dans le cas où il y a dissipation d'énergie, un terme de relaxation apparaît

$\frac{\mathrm{d}\overrightarrow{M}}{\mathrm{d}t}=\gamma (\overrightarrow{M}\wedge \overrightarrow{B})+\overrightarrow{R}$ avec $\overrightarrow{R}=\frac{\delta }{||\overrightarrow{M}||}(\overrightarrow{M}\wedge \frac{\partial \overrightarrow{M}}{\partial t})$

Gyrotropie

Dans cette partie, on va s'attacher à déterminer le tenseur de susceptibilité magnétique ${\overline{\overline{\chi }}}_m$, défini par :

$\overrightarrow{M}={\overline{\overline{\chi }}}_m\overrightarrow{H}$

Pour ce faire, on remplace $\overrightarrow{B}$ par $\overrightarrow{H}$ dans l'équation de précession en sachant que $\overrightarrow{B}={\mu }_0(\overrightarrow{M}+\overrightarrow{H})$:

$\begin{array}{cc}\frac{\mathrm{d}\overrightarrow{M}}{\mathrm{d}t}& =\gamma {\mu }_0(\overrightarrow{M}\wedge (\overrightarrow{M}+\overrightarrow{H}))+\overrightarrow{R}\\ & =-\gamma (\overrightarrow{M}\wedge \overrightarrow{H})+\overrightarrow{R}\end{array}$

On pose ${\gamma }_m=-\gamma {\mu }_0$ le coefficient magnétomécanique.

Décomposons les champs $\overrightarrow{M}$ et $\overrightarrow{H}$ en deux composantes :

• Une composante d'équilibre (on a ${\overrightarrow{M}}_0$ et ${\overrightarrow{H}}_0$ colinéaires)
• Une composante dépendante du temps

$\{\begin{array}{c}\overrightarrow{M}={\overrightarrow{M}}_0+\overrightarrow{M}(t)\\ \overrightarrow{H}={\overrightarrow{H}}_0+\overrightarrow{H}(t)\end{array}$

Pour trouver ${\overline{\overline{\chi }}}_m$ grâce à l'équation $\overrightarrow{M}(t)={\overline{\overline{\chi }}}_m\overrightarrow{H}(t)$, il faut linéariser le système, ce qui est possible lorsque

$\{\begin{array}{c}||\overrightarrow{M}(t)||\ll ||{\overrightarrow{M}}_0||\\ ||\overrightarrow{H}(t)||\ll ||{\overrightarrow{H}}_0||\end{array}$

$\begin{array}{cc}\frac{\mathrm{d}\overrightarrow{M}}{\mathrm{d}t}& =-{\gamma }_m(\overrightarrow{M}\wedge \overrightarrow{H})+\overrightarrow{R}\\ & =-{\gamma }_m(({\overrightarrow{M}}_0+\overrightarrow{M}(t))\wedge ({\overrightarrow{H}}_0+\overrightarrow{H}(t)))+\frac{\delta }{||{\overrightarrow{M}}_0+\overrightarrow{M}(t)||}(({\overrightarrow{M}}_0+\overrightarrow{M}(t))\wedge \frac{\partial ({\overrightarrow{M}}_0+\overrightarrow{M}(t))}{\partial t})\\ & =-{\gamma }_m({\overrightarrow{M}}_0\wedge \overrightarrow{H}(t)+\overrightarrow{M}(t)\wedge {\overrightarrow{H}}_0+\overrightarrow{M}(t)\wedge \overrightarrow{H}(t))+\frac{\delta }{||{\overrightarrow{M}}_0+\overrightarrow{M}(t)||}({\overrightarrow{M}}_0\wedge \frac{\partial \overrightarrow{M}(t)}{\partial t}+\overrightarrow{M}(t)\wedge \frac{\partial \overrightarrow{M}(t)}{\partial t})\\ & \approx -{\gamma }_m({\overrightarrow{M}}_0\wedge \overrightarrow{H}(t)+\overrightarrow{M}(t)\wedge {\overrightarrow{H}}_0)+\frac{\delta }{||{\overrightarrow{M}}_0||}({\overrightarrow{M}}_0\wedge \frac{\partial \overrightarrow{M}(t)}{\partial t}+\overrightarrow{M}(t)\wedge \frac{\partial \overrightarrow{M}(t)}{\partial t})\\ \end{array}$

Les grandeurs étant sinusoïdales, $\frac{\partial }{\partial t}=j\omega$, donc :

$\begin{array}{cc}j\omega \overrightarrow{M}(t)& =-{\gamma }_m{\overrightarrow{M}}_0\wedge \overrightarrow{H}(t)-{\gamma }_m\overrightarrow{M}(t)\wedge {\overrightarrow{H}}_0+\frac{\delta }{M_0}(j\omega M_0{\overrightarrow{u}}_z\wedge \overrightarrow{M}(t)+j\omega \overrightarrow{M}(t)\wedge \overrightarrow{M}(t))\\ & \approx -{\gamma }_m{\overrightarrow{M}}_0\wedge \overrightarrow{H}(t)-{\gamma }_m\overrightarrow{M}(t)\wedge {\overrightarrow{H}}_0+\frac{\delta }{M_0}j\omega M_0{\overrightarrow{u}}_z\wedge \overrightarrow{M}(t)\\ & =-{\gamma }_m{M}_0{\overrightarrow{u}}_z\wedge \overrightarrow{H}(t)+{\gamma }_m{H}_0{\overrightarrow{u}}_z\wedge \overrightarrow{M}(t)+j\delta \omega {\overrightarrow{u}}_z\wedge \overrightarrow{M}(t)\\ & ={\overrightarrow{u}}_z\wedge (-{\gamma }_m{M}_0\overrightarrow{H}(t)+{\gamma }_m{H}_0\overrightarrow{M}(t)+j\delta \omega \overrightarrow{M}(t))\\ \end{array}$

On pose $\{\begin{array}{c}{\underset{\_}{\omega }}_H={\gamma }_m{H}_0+j\delta \omega \\ {\omega }_H=\mathrm{\Re }({\underset{\_}{\omega }}_H)={\gamma }_m{H}_0\\ {\omega }_M={\gamma }_m{M}_0\end{array}$

On aboutit à $j\omega \overrightarrow{M}(t)={\overrightarrow{u}}_z\wedge ({\underset{\_}{\omega }}_H\overrightarrow{M}-{\omega }_M\overrightarrow{H})$

On projette suivant ${\overrightarrow{u}}_x$ et ${\overrightarrow{u}}_y$ :

$\{\begin{array}{c}j\omega M_x=-{\underset{\_}{\omega }}_H{M}_y+{\omega }_M{H}_y\\ j\omega M_y={\underset{\_}{\omega }}_H{M}_x-{\omega }_M{H}_x\end{array}\mathrm{d}\mathrm{o}\mathrm{n}\mathrm{c}\{\begin{array}{c}j\omega M_x+{\underset{\_}{\omega }}_H{M}_y={\omega }_M{H}_y\\ -{\underset{\_}{\omega }}_H{M}_x+j\omega M_y=-{\omega }_M{H}_x\end{array}$

Le but du calcul étant de trouver ${\overline{\overline{\chi }}}_m$, on inverse ce système :

$\{\begin{array}{c}{M}_x=\frac{{\omega }_M({\underset{\_}{\omega }}_H{H}_x+j\omega H_y)}{{\underset{\_}{\omega }}_H^2-{\omega }^2}\\ M_y=\frac{{\omega }_M(-j\omega H_x+{\underset{\_}{\omega }}_H{H}_y)}{{\underset{\_}{\omega }}_H^2-{\omega }^2}\end{array}\mathrm{s}\mathrm{o}\mathrm{i}\mathrm{t}\{\begin{array}{c}{M}_x=\frac{{\omega }_M{\underset{\_}{\omega }}_H}{{\underset{\_}{\omega }}_H^2-{\omega }^2}{H}_x+j\frac{{\omega }_M\omega }{{\underset{\_}{\omega }}_H^2-{\omega }^2}{H}_y\\ M_y=-j\frac{{\omega }_M\omega }{{\underset{\_}{\omega }}_H^2-{\omega }^2}{H}_x+\frac{{\omega }_M{\underset{\_}{\omega }}_H}{{\underset{\_}{\omega }}_H^2-{\omega }^2}{H}_y\end{array}$

On pose

$\chi =\frac{{\omega }_M{\underset{\_}{\omega }}_H}{{\underset{\_}{\omega }}_H^2-{\omega }^2}$

${\chi }_a=\frac{{\omega }_M\omega }{{\underset{\_}{\omega }}_H^2-{\omega }^2}$

Modèle:Encadre

Le tenseur ${\overline{\overline{\chi }}}_m$ est hermitien. Cette propriété est à la source du phénomène de gyrotropie.

Résonance ferromagnétique

Composante résonante

Pour étudier la résonance, on va s'intéresser aux composantes circulaires de $\overrightarrow{M}$ et $\overrightarrow{H}$ dans la base $({\overrightarrow{u}}_x,{\overrightarrow{u}}_y)$.

On introduit ainsi les composantes polarisées à droite (indice +) et à gauche (indice -) de ces deux champs :

$\{\begin{array}{c}{M}_{+}=M_x+j{M}_y\\ H_{+}=H_x+j{H}_y\end{array}\mathrm{e}\mathrm{t}\{\begin{array}{c}{M}_{-}=M_x-j{M}_y\\ H_{-}=H_x-j{H}_y\end{array}$

On introduit de même les susceptibilités magnétiques à droite et à gauche par les relations suivantes :

$\{\begin{array}{c}{M}_{+}={\chi }_{+}{H}_{+}\\ M_{-}={\chi }_{-}{H}_{-}\end{array}$

La relation matricielle $\left[\begin{array}{c}{M}_x\\ M_y\end{array}\right]=\left[\begin{array}{cc}\chi & j{\chi }_a\\ -j{\chi }_a& \chi \end{array}\right]\left[\begin{array}{c}{H}_x\\ H_y\end{array}\right]$ donne après développement les expressions de ces susceptibilités :

$\{\begin{array}{c}{\chi }_{+}=\chi +{\chi }_a=\frac{{\omega }_M({\underset{\_}{\omega }}_H+\omega )}{{\underset{\_}{\omega }}_H^2-{\omega }^2}=\frac{{\omega }_M}{{\underset{\_}{\omega }}_H-\omega }\\ {\chi }_{-}=\chi -{\chi }_a=\frac{{\omega }_M({\underset{\_}{\omega }}_H-\omega )}{{\underset{\_}{\omega }}_H^2-{\omega }^2}=\frac{{\omega }_M}{{\underset{\_}{\omega }}_H+\omega }\end{array}$

Modèle:Propriété

Étude de la résonance

Séparons partie réelle et partie imaginaire. On pose ${\chi }_{+}={{\chi }_{+}}^{\prime }-j{{\chi }_{+}}^{″}$

• La partie réelle influe sur la réfraction du milieu
• La partie imaginaire représente l'atténuation due aux pertes dans le milieu.

$\begin{array}{cc}{\chi }_{+}& =\frac{{\omega }_M}{{\underset{\_}{\omega }}_H-\omega }\\ & =\frac{{\omega }_M({\underset{\_}{\omega }}_H^{*}-\omega )}{({\underset{\_}{\omega }}_H-\omega )({\underset{\_}{\omega }}_H^{*}-\omega )}\\ & =\frac{{\omega }_M({\underset{\_}{\omega }}_H^{*}-\omega )}{{\underset{\_}{\omega }}_H{\underset{\_}{\omega }}_H^{*}-\omega ({\underset{\_}{\omega }}_H^{*}+{\underset{\_}{\omega }}_H)+{\omega }^2}\\ & =\frac{{\omega }_M({\underset{\_}{\omega }}_H^{*}-\omega )}{{\omega }_H^2+(\delta \omega {)}^2-2\omega {\omega }_H+{\omega }^2}\\ & =\frac{{\omega }_M({\omega }_H-\omega )}{({\omega }_H-\omega {)}^2+(\delta \omega {)}^2}-j\frac{\delta {\omega }_M\omega }{({\omega }_H-\omega {)}^2+(\delta \omega {)}^2}\\ \end{array}$

La partie réelle de la susceptibilité magnétique change brusquement de signe au voisinage de la fréquence de résonance. C'est également (et surtout) l'endroit où l'atténuation est la plus forte. En pratique, on repère la résonance ferromagnétique par le fait que l'intensité reçue par les détecteurs chute brusquement.