Table des Matières

Introduction Générale 4
0.1 Note historique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
0.2 Problématique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
0.3 Contribution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

1 Rappels et Problématique 8
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2 Rappels des concepts généraux . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.1 Matériaux magnétiques . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.2 Les équations de Maxwell . . . . . . . . . . . . . . . . . . . . . . 20
1.2.3 Matériaux diélectriques . . . . . . . . . . . . . . . . . . . . . . . . 25
1.3 Problématique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.3.1 Position du problème . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29

2 Résolution du Problème et Discussion 32

2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.2 Les matériaux à Bande Interdite Photonique . . . . . . . . . . . . . . . . 33
2.2.1 Historique des matériaux BIP . . . . . . . . . . . . . . . . . . . . 33
2.2.2 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.2.3 Propagation d’une onde TEM . . . . . . . . . . . . . . . . . . . . 34

1
 2.2.4 La méthode de la matrice de transfert . . . . . . . . . . . . . . . 36
2.2.5 Les paramètres effectifs . . . . . . . . . . . . . . . . . . . . . . . . 42
2.2.6 Les structures qui utilisent des gaps photoniques . . . . . . . . . . 50
2.3 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51

3 Résolution du Problème et Discussion 53

3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.2 Résolution du problème . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.2.1 Les multicouches . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.2.2 Le tenseur de perméabilité ferromagnétique . . . . . . . . . . . . 55
3.2.3 Les équations de Maxwell . . . . . . . . . . . . . . . . . . . . . . 63
3.2.4 Matrice de transfert et l’équation de dispersion . . . . . . . . . . 67
3.2.5 Les ondes de volume et les ondes de surface . . . . . . . . . . . . 70
3.2.6 Le spectre du coefficient de réflexion . . . . . . . . . . . . . . . . 72
3.2.7 Les applications des multicouches . . . . . . . . . . . . . . . . . . 78
3.3 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80

4 Comparaisons et Résultats 82
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
4.2 Comparaison entre les deux types de propagation dans les structures mul-
ticouches . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
4.2.1 La propagation des ondes . . . . . . . . . . . . . . . . . . . . . . 82
4.2.2 Le coefficient de réflexion et transmission . . . . . . . . . . . . . . 83
4.2.3 La matrice de transfert . . . . . . . . . . . . . . . . . . . . . . . . 85
4.2.4 Les propriétés dispersives . . . . . . . . . . . . . . . . . . . . . . . 87
4.3 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89

Conclusion Générale 91

Bibliographie 93

2
Introduction Générale

3
 Introduction générale

0.1 Note historique

L’étude des propriétés et comportements des matériaux magnétiques sont toujours d’ac-
tualité à travers le monde offrant toujours de nouvelles perspectives. Cette étude a com-
mencé durant les moyens âges.. A la fin du 16ème siècle William Gilbert (1544 − 1603)
réalisait une étude expérimentale du champ magnétique d’une sphère aimantée qui l’a-
menait à assimiler la terre elle-même à un grand aimant sphérique [1] [2].
La première structure multicouche périodique [14] a été fait par l’astronome David Rit-
tenhouse en 1785. Il a eu les vis très fines coupées par horloger sur deux morceaux de
petit fil en cuivre, et dans le fil de ces vis il a mis des poils. Cette structure simple était
suffisante pour démontrer l’effet de diffraction légère. Rittenhouse a publié son travail,
mais n’a continué pas plus depuis lui a attiré peu d’attention alors.
Ce n’était pas jusqu’en 1813 que le physicien Joseph von Fraunhofer a réinventé le réseau
de diffraction. Avec ce moteur il a produit des structures avec des périodes beaucoup
plus agréables que précédemment réalisées. Les décennies suivantes, le développement
ultérieur des structures a été particulièrement permis par des spectroscopiques. Plus
tard d’autres scientifiques et ingénieurs ont été intéressés aussi, réalisant une meilleure
exécution en construisant de plus grandes structures.
Le concept des structures BIP a été développé pour la première fois par Bragg il y a moins
d’un siècle par l’intermédiaire du miroir de Bragg. Composé d’un empilement de surfaces
planes transparentes, ce dernier grâce à des phénomènes d’interférences constructives de
réfléchir 99, 5% de l’énergie incidente à condition que cette dernière attaque le miroir sous
une incidence proche de zéro degré.
Les premières structures photoniques sont apparues sous le nom de diélectriques artifi-
ciels [10]. En effet le but était de reproduire un réseau cristallin ayant des propriétés

4
particulières. C’est E. Yablonovitch qui introduisit en premier le nom de structures à
bande interdite photonique, par analogie avec les structures cristallines et leur bande
électronique interdite. Depuis 5 ans, l’activité de recherche s’est énormément consacrée
à ces structures et de nombreuses applications ont été développées.

0.2 Problématique
Notre problème consiste à l’étude de la propagation des ondes électromagnétiques dans
différents milieux tel que les milieux diélectriques alternés avec des couches d’air qui
représente des bandes interdites et les milieux diélectriques alternés avec des couches
ferromagnétiques. Dans les deux cas nous avons utilisé la méthode de la matrice de
transfert pour calculer la relation de dispersion.

0.3 Contribution
Notre travail concerne des problèmes les plus simples qui mettent en jeu l’effet des champs
électromagnétiques avec des milieux matériels. Nous sommes alors amenés à une descrip-
tion rapide du comportement des ondes électromagnétiques dans les milieux matériels.
Pour une structure multicouche nous avons étudié le tenseur de perméabilité ferromagné-
tique sous l’influence d’un champ magnétique extérieur. Ensuite, nous avons déterminé
les équations de Maxwell pour les différentes couches. A partir de la matrice de transfert
nous avons déduit la relation de dispersion qui nous a permis de préciser les ondes de
volume et de surface. Afin de calculer les coefficients de réflexion et de transmission,
nous avons fait appel à la propagation de l’onde plane.
Ce mémoire est décomposé en quatre chapitres :
1. Dans le premier chapitre, on expose quelques notions générales de notre étude. Dans
un premier temps, nous rappelons les propriétés des matériaux magnétiques et diélec-
triques. Puis, on définit les modes de propagation TE et TM. Enfin, ce chapitre se

5
termine par une position du problème.
2. Dans le deuxième chapitre, nous avons abordé l’étude de la propagation d’une onde
TEM dans une structure multicouche périodique de même nature.
3. Dans le troisième chapitre, nous somme occupée à la résolution du problème en se
basant sur la méthode de la matrice de transfert pour le calcul de la relation de dispersion
pour une multicouche périodique, ainsi que la propagation de l’onde plane pour le calcul
des coefficients de réflexion et de transmission.
4. Dans le dernier chapitre, nous avons fait une comparaison entre les deux études.
Ce mémoire se termine par une conclusion sur les travaux présentés et par les perspectives
ouvertes par ceux-ci.

6
 Premier Chapitre

Rappels et Problématique

7
Chapitre 1

Rappels et Problématique

1.1 Introduction
Afin d’étudier la propagation d’une onde électromagnétique dans un milieu périodique
ferromagnétique diélectrique sous l’effet d’un champ magnétique extérieur, on doit pren-
dre en compte les interactions entre le champ électromagnétique et le milieu. L’approche
classique que nous considérons s’appuie sur les équations de Maxwell en présence du
milieu matériel. On adjoint alors aux équations des champs les relations du milieu qui
rendent compte des différents mécanismes d’interaction.
Ce chapitre est consacré à une présentation de quelques notions de base sur les ondes TE
et TM et le milieu de propagation c’est à dire les matériaux magnétiques et les matériaux
diélectriques. Nous poserons par la suite la problématique de notre travail.

8
1.2 Rappels des concepts généraux

1.2.1 Matériaux magnétiques

Définition

Soumises à une induction magnétique, certaines substances se mettent à produire elles-

mêmes, dans le volume qu’elles occupent et à l’extérieur, une induction magnétique.
On dit qu’elle s’aimantent ou se polarisent magnétiquement. Il s’agit-là d’une propriété
générale de la matière. Toutefois, cette propriété ne se manifeste très visiblement que
dans certains matériaux appelés matériaux magnétiques [1].

Les origines électroniques du magnétisme

Les propriétés magnétiques d’un matériau sont attribuables au spin des électrons et, à
leur mouvement orbital autour du noyau (Figure (1.1)) [2].

Figure (1.1) : Structure atomique

a) Moment magnétique élémentaire, b) Mouvement des électrons.

Les électrons qui tournent sur eux-mêmes “spin” et autour du noyau “orbitale” forment
de petits dipôles magnétiques qui peuvent être simulés par des moments magnétiques

9
engendrés par des boucles élémentaires de courant (Figure (1.1(b))).
Moment magnétique orbital :
Ce moment est lié au moment cinétique (Gl ) résultant du mouvement orbital par la
relation :

e
ml = − Gl (1.1)
2m

La direction opposée des deux moments est due à la charge négative de l’électron.
Pour rendre compte du spectre de l’atome H (hydrogène), Bohr a introduit une idée,
qui consiste à imposer que le module du moment cinétique orbital de l’électron soit un
h
multiple entier de , où h désigne la constante de Planck :
2π

h
Gl = l (l entier = 0, 1, 2, .....) (1.2)
2π

Le moment magnétique ml est alors nécessairement un multiple entier d’un moment unité
qu’on appelle le magnéton de Bohr :

Gl = lµB (1.3)

Ou :

eh
µB = (1.4)
4πm

Le magnéton de Bohr est une grandeur fondamentale en magnétisme. Il vaut :

µB = 9, 273.10−24 A.m2 . C’est l’unité de moment magnétique à l’échelle atomique.

10
Moment magnétique de spin :
En tournant sur lui-même l’électron crée un moment cinétique de spin Gs qui peut pren-
1 h 1 h
dre deux valeurs distinctes Gs = + et Gs = − .
2 2π 2 2π
Proposée initialement par S.A. Goudsmit et G. Uhlenbeck pour rendre compte des spec-
tres atomiques (1925), la notion de spin a été justifiée ensuite par Dirac en application
de la mécanique quantique [1]. A ce moment cinétique de spin est associé un moment
magnétique de spin tel que :

e
ms = − Gs (1.5)
m

Il apparaît que le spin est deux fois plus efficace qu’un moment orbital dans la production
du moment magnétique.
Moment magnétique total :
Le moment magnétique atomique total est la somme des moments magnétiques orbitaux
et de ceux produits par le spin. Dans une situation fréquemment rencontrée appelée
“couplage de Russel-Saunders” ou couplage L.S, les moments de spin se couplent entre
eux pour donner le moment de spin total [2] :

S= Si (1.6)
i

Il est de même que pour les moments orbitaux, ils se couplent entre eux pour donner le
moment orbital :

L= mli (1.7)
i

11
Enfin les moments orbitaux totaux et de spin L et S se couplent entre eux pour aboutir
à un moment cinétique total J qui peut prendre toutes les valeurs comprises entre :

J = |L + S| et J = |L − S|

A partir du moment cinétique total J, on obtient facilement le moment magnétique de

l’atome :

m = JgµB (1.8)

Où g est un nombre qu’on appelle “facteur gyromagnétique de Landé” qui vaut :

J (J + 1) + S (S + 1) − L (L + 1)
g =1+
2J (J + 1)

Définition de l’aimantation

Un atome avec des moments non compensés, issus du mouvement de spin et orbital des
électrons est dit atome magnétique. Le moment magnétique totale de ces atomes est égale
à la somme vectorielle de tous les moments magnétiques proviennent du mouvement de
spin ou orbital des électrons [2].
Soit l’élément de volume (∆V ) constitué d’un nombre N d’atomes magnétiques chaque
atome est représenté par le moment magnétique mi , i = 1, ......N. La somme vectorielle
de ces moments donne le moment magnétique m comme suit :

N
m= mi (1.9)
i=1

12
L’aimantation du matériau ou le vecteur aimantation M est caractérisée par la densité
volumique des moments magnétiques et on écrit :

N
1
M = lim mi (1.10)
∆V →0 ∆V i=1

A la place de l’aimantation, les spécialistes des matériaux qui travaillent en construction

électrique préfèrent utiliser généralement la polarisation magnétique :

J = µ0 .M (1.11)

J s’exprime en Teslas, unité commode pour évaluer les performances de nombreux matéri-
aux magnétiques industriels.

Aimantation d’un milieu matériel

La réaction d’un matériau à l’application d’un champ magnétique H est caractérisée

par l’induction magnétique B, laquelle représente l’intensité du champ magnétique à
l’intérieur du matériau. A la présence de cette induction magnétique, un couple s’exerce
sur le dipôle magnétique de moment m :

Γ=m×B (1.12)

Ce qui signifie que l’induction magnétique B tend à aligner le dipôle de façon que le
moment m soit parallèle à l’induction. L’énergie potentielle du moment magnétique est
donnée par :
W = −m.B (1.13)

13
Dans le vide, cette énergie peut être exprimée par W = −µ0 mH, ou µ0 = 4π10−7 H/m
est la perméabilité du vide (l’air).
Pour évaluer la réaction d’un matériau (mesurée par B) à l’application d’un champ mag-
nétique excitateur H, on la compare à la réaction du vide à l’application du même champ.
Dans ce cas, l’induction magnétique est égale à la somme de l’induction magnétique du
vide µ0 H et de l’aimantation du matériau µ0 M :

B = µ0 (H + M) (1.14)

B
Ou H = − M est l’intensité du champ magnétique.
µ0
Pour de nombreuses substances magnétiques, M est reliée à H par l’égalité suivante :

M = χH (1.15)

Où χ est la susceptibilité magnétique de la substance. Cette susceptibilité n’est pas

nécessairement constante. Elle peut varier en fonction du champ appliqué et présente un
comportement à haute température (dans la région de la température critique Tc ). Ainsi,
dans les matériaux anisotropes la susceptibilité a la forme d’un tenseur M = [χ] H :
⎡ ⎤ ⎡ ⎤⎡ ⎤
⎢ Mx ⎥ ⎢ χxx χxy χxz ⎥ ⎢ Hx ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢ My ⎥ = ⎢ χyx χyy χyz ⎥ ⎢ Hy ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
(1.16)
⎣ ⎦ ⎣ ⎦⎣ ⎦
Mz χzx χzy χzz Hz

Dans ce cas, le champ H et le vecteur aimantation M, ont différentes directions. En

substituant l’équation (1.15) dans (1.14) on peut écrire :

B = µ0 (1 + χ)H = µ0 µr H (1.17)

14
On caractérise une substance par sa perméabilité magnétique relative par rapport au vide,
µr laquelle est déterminée par la relation :

µr = 1 + χ (1.18)

Matériaux ferromagnétiques

Le ferromagnétisme est le type de magnétisme résultant de l’alignement de moments mag-

nétiques permanents, ces moments étant orientés parallèlement les uns aux autres par
une interaction mutuelle (d’échange) appelée “couplage ferromagnétique” (Figure (1.2)).
Les matériaux ferromagnétiques présentent donc également une aimantation spontanée
importante indispensable dans les applications en génie électrique. Ils ont aussi une
température de Curie Tc , au dessus de laquelle ils deviennent paramagnétiques, leur sus-
ceptibilité suivant alors la loi de Curie-Weiss [1].

Figure (1.2) : Représentation schématique de l’alignement des moments

magnétiques dans un matériau ferromagnétique.

La perméabilité des matériaux ferromagnétiques varie de façon importante avec l’induc-

tion et dépend du passé magnétique de l’échantillon étudié, avec des valeurs différentes
selon que le champ est obtenu par valeurs croissantes ou décroissantes “phénomène d’hys-
térésis”. En particulier une substance ferromagnétique vierge soumise à l’influence d’un
champ magnétique conserve une aimantation plus ou moins importante après disparition
du champ extérieur “phénomène de rémanence” et un champ inverse et nécessaire pour

15
la faire disparaître “phénomène de coercivité”.
Le Fer (F e), le Cobalt (Co), le Nickel (Ni) et un certain nombre de leurs alliages sont fer-
romagnétiques. Quelques terres rares et certains alliages de Manganèse avec l’Aluminium
et le Cuivre le sont aussi.

Domaines magnétiques

Le phénomène d’orientation des moments magnétiques à l’intérieur des matériaux ferro-

magnétiques se produit de manière localisée. L’expérience montre que tous les matériaux
magnétiques de volume suffisant se divisent spontanément en régions plus petites qu’on
appelle “domaines magnétiques” ou “domaines de Weiss” (ils peuvent mesurer jusqu’à
1mm). A l’intérieur d’un domaine tous les moments magnétiques atomiques sont main-
tenus parallèles par les forces d’échange, de telle sorte que chaque domaine se présente
comme un petit volume aimanté à saturation (Js ) avec une orientation différente de celles
de ces voisins de telle sorte que la polarisation globale de l’échantillon peut forte et bien
être nulle. Historiquement, cette idée d’une structure en domaine a été introduite pour
la première fois par P.Weiss en 1907 pour lever l’incompatibilité entre l’existence d’un
arrangement ordonné des moments atomiques et la possibilité d’un état désaimanté. Les
interfaces entre les domaines sont appelées “parois de Bloch” (1932) à travers lesquelles
l’orientation des moments magnétiques passe progressivement d’un domaine à l’autre
(Figure (1.3)) [2].

16
 Figure (1.3) : a) Domaines magnétiques, b) Changement d’orientation des
moments magnétiques à travers une paroi de Bloch.

Processus d’aimantation et cycle d’hystérésis

Dans les paragraphes précédents, nous avons vu la structure en domaines des matériaux
ferromagnétiques. Cette structure se manifeste à une échelle qu’on appelle échelle méso-
scopique (entre l’échelle atomique et le domaine visible). Dans ce présent paragraphe nous
abordons comment une telle structure rend compte, à notre échelle (macroscopique), des
lois d’aimantation observées. Pour ce faire, considérons un monocristal ferromagnétique
dont l’orientation des domaines est globalement nulle en l’absence de champ magnétique
extérieur H. (Figure (1.4)).

17
 Figure (1.4) : En l’absence de champ magnétique extérieur, la résultante de
l’orientation des domaines magnétiques d’un monocristal ferromagnétique est nulle.

Courbe de première aimantation :

Lorsqu’on applique un champ magnétique extérieur à ce matériau, les domaines mag-
nétiques dont l’orientation est voisine de celle de H croissent aux dépens de ceux dont
l’orientation est opposée à celle de H. Les parois de Bloch se déplacent donc dans le
matériau. A la limite, lorsque le champ magnétique extérieur atteint un valeur critique
Hs , le monocristal n’est plus constitué que d’un seule domaine ferromagnétique, dont
l’orientation est la même que celle de H; l’induction atteint alors une valeur maximale,
Bs appelée induction à saturation. Ce phénomène est représenté par la courbe OA de la
Figure (1.5) [2], courbe appelée courbe de première aimantation ou courbe de magnétisa-
tion.

18
 Figure (1.5) : Courbe de première aimantation (OA) d’un matériau ferromagnétique,
et cycle d’hystérésis (ABCDEFA) de ce matériau, lorsque le champ extérieur varie
alternativement de + Hs à − Hs .

Cycle d’hystérésis :
Lorsque l’intensité du champ H diminue pour atteindre une valeur nulle, les domaines
magnétiques ont tendance à réapparaître (courbe AB, Figure (1.5)); toutefois, puisque le
déplacement des parois de Bloch n’est pas instantané à cause de l’anisotropie magnétique,
une induction non nulle Br , se manifeste dans le matériau (point B, Figure (1.5)). Cette
valeur Br est appelée induction rémanente. Il faut en fait appliquer un champ magnétique
de sens opposé à celui du champ de première aimantation pour que l’induction soit nulle
(courbe BC, Figure (1.5)). La valeur Hc du champ qui engendre cette induction nulle
correspond au champ coercitif (point C, Figure (1.5)). Lorsque l’intensité du champ H

19
croit, l’induction atteint de nouveau la valeur maximale Bs (courbe CD, Figure (1.5)).
Enfin lorsqu’on change le sens du champ et qu’on fait varier son intensité; on obtient
la courbe DEFA (Figure (1.5)). La courbe d’aimantation ainsi obtenue est une boucle
d’hystérésis ou cycle d’hystérésis (ABCDEFA), dont les paramètres caractéristiques sont
les grandeurs Bs , Br et Hc .

1.2.2 Les équations de Maxwell

Les équations de Maxwell représentent les bases de l’électromagnétisme, c’est à dire que
ces équations permettent de décrire les évolutions spatio-temporelles du champ électrique
et du champ magnétique.
Dans le cas d’un diélectrique parfait, elles s’écrivent de la manière suivante [3] :

⎧ −
→
⎪
⎪ −→− → ∂B
⎪
⎪ Rot E = −
⎪
⎪
⎪
⎪ −→∂t
⎪
⎪
⎨ −→− → ∂D
Rot H = (1.19)
⎪
⎪ ∂t
⎪
⎪ −→
⎪
⎪ div D = 0
⎪
⎪
⎪
⎪ −→
⎩ div B = 0

−
→
Où E : est le champ électrique exprimé en V.m−1 .
−
→
B : est le champ magnétique exprimé en Tesla (T ).
−
→
D : l’induction électrique (déplacement électrique) exprimé en C.m−2 .
ε : la permittivité diélectrique absolue du milieu exprimée en F.m−1 , ε = ε0 εr , εr est la
−
→ −
→
permittivité relative du milieu, dans le vide, D = ε0 E .
−→
∂D
Le terme est communément appelé densité de courant de déplacement.
∂t
t : temps exprimé en secondes (s).
µ : la perméabilité magnétique absolue du milieu exprimée en H.m−1 , µ = µ0 µr , µr est
la perméabilité magnétique relative.
−
→ −→
µr = 1 pour un milieu “non magnétique”, dans le vide B = µ0 H .

20
Développons les deux premières équations de Maxwell afin d’obtenir les relations entre
les composantes des champs.
Première équation de Maxwell :

−→
−→− → ∂H
Rot E = −µ
∂t

∂
Pour un champ électromagnétique périodique dans le temps, on peut considérer = iω.
∂t

−
→ − → −→
∇ ∧ E = −iωµ H

∂Ez ∂Ey
− = −iωµHx (1.20)
∂y ∂z
∂Ex ∂Ez
− = −iωµHy (1.21)
∂z ∂x
∂Ey ∂Ex
− = −iωµHz (1.22)
∂x ∂y

La deuxième équation donne (on pose ε = ε0 εr ) :

−→
−→− → ∂E
Rot H = ε
∂t

∂
Pour un champ électromagnétique périodique dans le temps, on peut considérer = iω.
∂t

−
→ − → −→
∇ ∧ H = iωε E

21
 ∂Hz ∂Hy
− = iωεEx (1.23)
∂y ∂z
∂Hx ∂Hz
− = iωεEy (1.24)
∂z ∂x
∂Hy ∂Hx
− = iωεEz (1.25)
∂x ∂y

En examinant ces équations, on remarque qu’il existe deux types de solutions sont pos-
sibles pour cet ensemble d’équations : des ondes avec un champ électrique transverse
“ondes TE” et des ondes avec un champ magnétique transverse “ondes TM”. Puisque
aucune variation d’indice ne survient en x, alors toutes les dérivées partielles par rapport
à x seront nulles. Les équations (1.20) à (1.25) deviennent alors séparables en 2 groupes
de 3 équations [4].

Modes Transverse Electrique (TE)

Mode TE.

Le champ électrique des ondes TE ne possède qu’une composante non nulle, il s’agit de

22
Ex . Le champ magnétique aura alors deux composantes non nulles, soit Hy et Hz [5].
Les équations (1.21), (1.22) et (1.23) s’écrivent alors :

∂Ex
= −iωµHy (1.26)
∂z
∂Ex
= iωµHz (1.27)
∂y
∂Hz ∂Hy
− = iωεEx (1.28)
∂y ∂z

Ce qui implique que l’équation (1.26) et (1.27) devient :

1 ∂Ex
Hy = − (1.29)
iωµ ∂z
β
Hz = − Ex (1.30)
ωµ

On remplace l’équation (1.29) et (1.30) dans (1.28) on aura :

∂ 2 Ex
2
+ ω 2 εµ − β 2 Ex = 0
∂z

Posons : k 2 = ω 2 εµ − β 2 .
On écrira la solution générale sous la forme :

Ex = A1 exp (ikz) + A2 exp (−ikz)

23
Modes Transverse Magnétique (TM)

Mode TM.

Le champ magnétique des ondes TM ne possède qu’une composante non nulle, il s’agit
de Hx . Le champ électrique aura alors deux composantes non nulles, soit Ey et Ez [5].
Les équations (1.20), (1.24) et (1.25) s’écrivent alors :

∂Ez ∂Ey
− = −iωµHx (1.31)
∂y ∂z
∂Hx
= iωεEy (1.32)
∂z
∂Hx
= −iωεEz (1.33)
∂y

Ce qui implique que l’équation (1.32) et (1.33) devient :

1 ∂Hx
Ey = (1.34)
iωε ∂z
β
Ez = Hx (1.35)
ωε

24
On remplace l’équation (1.34) et (1.35) dans (1.31) on aura :

∂ 2 Hx
2
+ ω 2 εµ − β 2 Hx = 0
∂z

Posons : k 2 = ω 2 εµ − β 2 .
On écrira la solution générale sous la forme :

 
Hx = A1 exp (ikz) + A2 exp (−ikz)

1.2.3 Matériaux diélectriques

Définition d’un diélectrique

Un matériau est diélectrique s’il ne contient pas de charges électriques susceptibles de

se déplacer de façon macroscopique. Autrement dit, c’est un milieu qui ne peut pas
conduire le courant électrique. A ce titre, on l’appelle parfois isolant électrique. On
compte parmi ces milieux le verre et de nombreux plastiques. Par exemple, les câbles
électriques sont souvent protégés d’un revêtement en plastique pour éviter que le courant
électrique puisse en sortir.
Malgré l’impossibilité des milieux diélectriques de conduire le courant, ils présentent de
nombreuses caractéristiques électriques. En effet les atomes qui constituent le matériau
peuvent présenter des dipôles électrostatiques qui sont susceptibles d’interagir avec un
champ électrique. Cette interaction se traduit par la création d’une polarisation reliée
à ce champ électrique, au niveau microscopique, par une polarisabilité, et au niveau
macroscopique, par la susceptibilité électrique [6].

−
→ −→
P = ε0 χ E (1.36)

25
ε0 : la permittivité diélectrique du vide.
χ : la susceptibilité électrique du matériau.

Définition de la constante diélectrique

La constante diélectrique, également nommée permittivité du vide, est une constante

physique. Elle est notée par ε0 .
La constante électrique est définie par :

1 1
ε0 = 2
= (1.37)
µ0 c 36π109

Où :
µ0 : est la constante magnétique.
c : est la vitesse de la lumière dans le vide.
Dans le système du SI, sa valeur est égale à :

ε0 = 8, 85.10−12 F/m (1.38)

Interprétation physique

La permittivité d’un matériau est définie comme le rapport entre la norme du champ
de déplacement électrique et celle du champ électrique appliqué au matériau. Pour des
champs suffisamment grands, ce rapport n’est pas constant et tend vers ε0 .
ε0 peut être vue comme la permittivité intrinsèque du vide.
Pour un matériau donné de permittivité ε, il est possible de définir la permittivité relative,
normalisée par rapport à celle du vide :

ε
εr = (1.39)
ε0

26
Cette valeur ne possède pas d’unité et est toujours supérieure à 1.

Permittivité relative de quelques isolants

Tableau de la permittivité relative de quelques isolants :

Matériaux Permittivité relative εr

Vide 1 (par définition)
Air 1,0005
Papier 2
Polyester 3,3
Téflon 2,1
PVC 5
Verre 5
Bakélite 6
Chlorure de sodium 6,12
Titanate de baryum à peu près 1500

1.3 Problématique
Notre problème se compose en deux parties :
Dans la première partie on s’intéresse à l’étude des propriétés dispersives de la propa-
gation d’une onde TEM dans une structure multicouche périodique avec deux couches
diélectriques de même nature, le matériau utilisé dans ce cas BIP (les matériaux à Bande
Interdite Photonique).
Notre but est de calculer la relation de dispersion et les paramètres effectifs du milieu
périodique : l’indice de réfraction, l’impédance d’onde, la permittivité électrique et la
perméabilité magnétique en utilisant la méthode de la matrice de transfert.
Dans la deuxième partie nous avons étudié les propriétés dispersives de la propagation des

27
ondes électromagnétiques dans une structure multicouche périodique avec deux couches
de nature différentes une couche ferromagnétique et une couche diélectrique sous l’effet
d’un champ magnétique extérieur.
Le but est de calculer la relation de dispersion en utilisant la méthode de la matrice de
transfert et le coefficient de réflexion pour une multicouche semi-infinie en utilisant la
propagation de l’onde plane.

1.3.1 Position du problème

On considère une structure multicouche finie, avec deux types de matériaux par période
(s = 2), voir la Figure (1.6). Chaque période se compose d’une couche de largeur d1 et
indice de réfraction n0 et une couche de largeur d2 et indice de réfraction n1 . De cette
façon, chaque période sera symétrique. Le nombre total de couches est l = N (2s − 1),
avec N le nombre de périodes dans la structure.

Figure (1.6) : Une structure multicouche périodique.

On définit les paramètres de la structure par :

28
 ⎧
⎪
⎨ n0 pour i = 1, 3, 5, ..., l − 1
ni = ⎪ (1.40)
⎩ n1 pour i = 2, 4, 6, ..., l
⎧
⎪
⎨ d1 pour i = 1, 3, 5, ..., l − 1
di = ⎪ (1.41)
⎩ d2 pour i = 2, 4, 6, ..., l
⎧
⎪
⎪
⎪
⎨ 0 pour i = 0
Li = ⎪ i (1.42)
⎪
⎪ dn pour i = 1, 2, 3, ..., l
⎩
n=1

On a divisé notre travail en deux parties :

- La première partie concerne la propagation de l’onde TEM dans deux couches diélec-
triques de même type.
- La deuxième partie concernant les propriétés dispersives de la propagation des ondes
électromagnétiques sous l’effet d’un champ magnétique extérieur pour les milieux péri-
odiques de nature différente, une alternance de couches ferromagnétique et diélectrique.
Notant que pour le calcul du coefficient de réflexion et transmission, on a utilisé la prop-
agation de l’onde plane qui entre avec une amplitude A est partiellement réfléchie et
partiellement transmise sur chaque interface dans la structure. L’amplitude de l’onde
réfléchie se propageant vers la gauche est dénotée par r et l’onde transmise se propageant
vers la droite est dénotée par t. Le coefficient de transmission de l’onde entrante par cette
structure est mesuré pour différentes fréquences. Pour chaque fréquence il est possible
de mesurer l’amplitude de l’onde transmise.

1.4 Conclusion
Dans ce chapitre, nous avons vu l’importance technologique des matériaux magnétiques,
qui mobilisent les chercheurs depuis des siècles. Nous avons présenté les notions de bases
du point de vue des sciences des matériaux qui nous ont permis de comprendre comment

29
le magnétisme s’insère dans l’ensemble plus vaste de la physique moderne des solides.

30
 Deuxième Chapitre

Résolution du Problème et Discussion

31
Chapitre 2

Résolution du Problème et
Discussion

2.1 Introduction
Nous débuterons ce chapitre par un historique retraçant l’origine des cristaux pho-
toniques. Après avoir définit les éléments permettant de caractériser ces cristaux, nous
pourrons inventorier les matériaux à bande interdite photonique.
Nous adresserons ensuite la question dans quelle mesure la structure périodique agit
comme un milieu homogène effectif pour l’onde optique. C’est important du point de
vue expérimental quand les paramètres matériels des différentes couches doivent être
extraite à partir des mesures conventionnelles de réflexion et de transmission. Nous
analyserons une structure simple de couche. Ceci permet une demande de règlement
unidimensionnelle qui donne la solution sous forme analytique.
Nous continuerons par étudier la propagation de l’onde TEM à la frontière entre deux
couches puis pour une période de la structure et enfin nous aurons généralisé pour N
périodes en utilisant la matrice de transfert.
Nous terminerons par une présentation des méthodes numérique utilisées pour effectuer
ce travail.

32
2.2 Les matériaux à Bande Interdite Photonique

2.2.1 Historique des matériaux BIP

La naissance des structures photoniques provient de l’optique. C’est en 1915 que le physi-
cien anglais Sir William Lawrence Bragg reçu le prix Nobel de physique en récompense
de son travail effectué sur la diffraction des rayons X par les cristaux [10].
Bragg a développé le miroir de Bragg. Ce miroir est une succession de surfaces planes
transparentes d’indices de réfraction différents. Il permet de réfléchir, grâce à des phénomènes
d’interférences constructives, 99, 5% de l’énergie incidente.
Ce n’est qu’en 1987 qu’Eli Yablonovitch de l’université de Los Angeles proposa d’étendre
le concept des miroirs de Bragg aux fréquences micro-ondes et ce pour des incidences
quelconques.
En 1991, Eli Yablonovitch réalisa le premier cristal photonique artificiel fonctionnant aux
longueurs d’onde centimétriques. Ce cristal artificiel, réalisé dans un bloc de plexiglas
usiné dans trois directions différentes de façon à reproduire la structure du diamant, est
appelé Yablonovite.

2.2.2 Définition

Les matériaux à Bande Interdite Photonique, plus simplement appelés BIP (en anglais,
Photonic Band Gap, PBG) sont composés d’un assemblage périodique de deux ou
plusieurs matériaux diélectriques ou métalliques. Ils présentent des bandes de fréquence
pour lesquelles la propagation des ondes électromagnétiques est interdite [11].

33
 Figure (2.1) : Structure périodique unidimensionnelle.

Cette structure est le dispositif périodique le plus simple. Il est connu sous le nom de
“Miroir de Bragg”. Il consiste en une alternance de couches planes de diélectrique.

2.2.3 Propagation d’une onde TEM

Considérons une onde TEM se propageant dans une direction transversale aux couches,
−
→
appelées la direction z. Pour établir l’équation de propagation du champ électrique E ,
nous pouvons écrire successivement [9] :

−→−→− → ∂ −→− →
RotRot E = −µ Rot H (2.1)
∂t

−→
−−→ − → −→ ∂2 E
graddiv E − ∆ E = −εµ 2 (2.2)
∂t

D’où :

−→
−→ ∂2 E
∆ E − εµ 2 = 0 (2.3)
∂t

34
Dans ce cas particulier, nous pouvons remplacer les dérivées par rapport au temps par
le facteur iω.

−→
d2 E −→
2
+ ω 2 εµ E = 0 (2.4)
dz

Ainsi, la solution de l’équation (2.4) est donnée par :

E (z) = A+ eikz + A− e−ikz (2.5)

Où A+ et A− sont des constantes. Nous avons introduit le nombre d’onde “k” défini
comme :

√ ω
k = ω εµ = n (2.6)
c

1
Où c = √ est la vitesse de la lumière dans le vide et n est l’indice de réfraction
µε
complexe.
Le champ magnétique H (z) lié au champ électrique E (z) dans l’équation (2.5) :

1 + ikz
H (z) = A e − A− e−ikz (2.7)
η

Où l’impédance d’onde est donnée par :

ωµ µ
η= = (2.8)
k ε

35
La Figure (2.2) montre la variation du nombre d’onde en fonction de la fréquence. Pour
le mode fondamental ou TEM : L’onde se propage parallèlement à l’axe de la structure
ω
périodique à toutes les fréquences. La constante de propagation k = n est égale à celle
c
de l’onde plane dans l’air. Ce mode peut se propager à toutes les fréquences jusqu’à la
fréquence nulle.

Figure (2.2) : La relation de dispersion.

2.2.4 La méthode de la matrice de transfert

Nous déterminons la matrice de transfert pour chaque période dans la structure avec
un nombre arbitraire de couches. Nous étudions une structure photonique finie à une
dimension constituée de plaque de diélectrique (le diélectrique utilisé dans ce cas c’est
le verre) (ε1 = 5) d’épaisseur a1 espacées par des gaps d’air (ε0 = 1) d’épaisseur a0 . La
propagation des ondes alors se fait dans cette région.

36
À la frontière entre deux couches :

Les composants tangentiels des champs doivent être continus. Ceci mène à une relation
matricielle entre l’amplitude d’onde de différents côtés :

Figure (2.3) : (a) Structure multicouche, frontière simple.

A0 = P0 (z0 ) R01 P1−1 (z0 ) A1 (2.9)

Dans ce cas précis ν (= 0, 1). Aν est une matrice colonne composée des amplitudes d’onde,
c’est à dire :
⎡ ⎤
A−
Aν = ⎢ ⎥
ν
⎣ ⎦ (2.10)
A+
ν

Pν est une matrice de propagation diagonale :

⎡ ⎤
ikν z
e 0
Pν (z) = ⎢
⎣
⎥
⎦ (2.11)
−ikν z
0 e

37
Et Rνµ est une matrice de réflexion 2 × 2 donnée par :
⎡ ⎤
1 ⎢ 1 rνµ ⎥
Rνµ = ⎣ ⎦ (2.12)
1 + rνµ rνµ 1

Avec :

1 − η ν /η µ
rνµ = (2.13)
1 + η ν /η µ

P0 R01 P1−1 dans l’équation (2.9) s’appelle la matrice de transfert. La structure entière
de couche peut être traitée maintenant par une technique directe de multiplication de
matrice.

Une seule période de la structure :

 
Les épaisseurs de couche sont a0 = a0 + a0 et a1 , et la période étant a = a0 + a1 . Pour
une seule période, la relation de matrice de transfert :

Figure (2.4) : (b) Structure multicouche, une seule période.

A0 = P0 (z0 ) × M × P0−1 (z1 ) A1 (2.14)

38
Avec :

   
M = P0 z − z0 R01 P1 z − z R10 P0 z1 − z (2.15)

−1  a0  a0
Notez que R10 = R01 . Nous définissons z − z0 = + ∆a0 et z1 − z = − ∆a0 , la
2 2
matrice M s’écrit ensuite :

M = P0 (∆a0 ) × MS × P0−1 (∆a0 ) (2.16)

Où MS est la matrice de transfert de la cellule d’unité symétrique.

⎡ ⎤
a b ⎥
MS = ⎢
⎣ ⎦ (2.17)
−b c

Avec :

1
a= 2
eik0 a0 eik1 a1 − r01
2 −ik1 a1
e
(1 − r01 )
r01
b= 2
e−ik1 a1 − eik1 a1 (2.18)
(1 − r01 )
1
c= 2
e−ik0 a0 e−ik1 a1 − r012 ik1 a1
e
(1 − r01 )

N périodes de la structure :

La matrice de transfert totale de la structure pendant N périodes est M N :

39
 Figure (2.5) : (c) Structure multicouche, N périodes.

M N = P0 (∆a0 ) × MSN × P0−1 (∆a0 ) (2.19)

La matrice M N est calculée par diagonaliser MS en employant la transformation :

⎡ ⎤
λ1 0 ⎥
X −1 MS X = Λ = ⎢
⎣ ⎦ (2.20)
0 λ2

En utilisant :

MSN = XΛN X −1 (2.21)

Avec les équations (2.17) et (2.18), les valeurs propres de λ1,2 sont :

λ1,2 = e±iθ (2.22)

40
Où θ est donné par :

1 2
cos θ = 2
cos (k0 a0 + k1 a1 ) − r01 cos (k0 a0 − k1 a1 ) (2.23)
(1 − r01 )

L’équation (2.23) est la relation de dispersion de Kronig-Penney [10], qui est familière
de la théorie de bande d’énergie de l’électron. En faisant une simple représentation
graphique, en posant y = cos θ et x = k0 a0 , il est possible de montrer que la relation
(2.23) admet des valeurs de cette fonction supérieures à +1 et inférieures à −1. Or
l’équation (2.23) n’admettra des solutions que la relation sera compris entre +1 et −1.
Ce comportement met en évidence la notion de bandes interdites et de bandes permises
décrites par la Figure (2.6).

Figure (2.6) : Bandes permises et bandes interdites.

La matrice de transformation X s’écrit :

⎡ ⎤
1 ⎢ 1 P ⎥
X= ⎣ ⎦ (2.24)
1+P P 1

41
Où :

(λ1 − a) b
P = = (2.25)
b (λ2 − 1)

En remplaçant M avec M N dans l’équation (2.14) et utilisant les équations (2.19) à

(2.22), la relation de matrice de transfert pour N périodes est :

⎡ ⎤ ⎡ ⎤ ⎡ ⎤
⎢ A− ⎥
iNθ
e 0 A−
= P0 (z0 ) P0 (∆a0 ) X ⎢ ⎥ −1 −1 ⎢ ⎥
0 −1 N
⎣ ⎦ ⎣ ⎦ X P0 (∆a0 ) P0 (zN ) ⎣ ⎦ (2.26)
A+
0 0 e−iNθ
A+
N

2.2.5 Les paramètres effectifs

La propagation de l’onde dans la structure périodique de couche est décrite par ΛN et


le passage du milieu environnant (ν = 0) par la matrice X = P0 (∆a0 ) X. Le nombre
d’onde kef f dans le matériau périodique est définit par [9] :

kef f Na = N (θ + 2πm) (2.27)

Où θ est limité à la branche principale −π ≤ Re (θ) ≤ π et m est un nombre entier lié

aux diverses zones de Brillouin.
La Figure (2.7) montre une représentation du diagramme de dispersion d’un matériau
périodique à une dimension. Ce diagramme représente le nombre d’onde effectif multi-
ωa
pliée par le pas du réseau divisé par 2π, en fonction de la fréquence . Nous pouvons
c
observer les bandes interdites et les bandes autorisées. Le nombre d’onde est purement
réel à l’intérieur des bandes.

42
 Figure (2.7) : Diagramme de dispersion.

Le milieu périodique est caractérisé par la relation de dispersion kef f = ωnef f /c de façon
analogue avec l’équation (2.6), nous trouvons l’indice de réfraction effective :

c
nef f = (θ + 2πm) (2.28)
ωa

Dans le cas symétrique ∆a0 = 0, et le passage de l’onde du milieu homogène au périodique

est décrit par la matrice X. Comme montré dans l’équation (2.24) il a exactement la
même forme que la matrice simple de réflexion de frontière Rνµ dans l’équation (2.12).
Les coefficients de réflexion ordinaires rνµ devraient être maintenant remplacés par P .

43
D’après l’équation (2.13), nous définissons l’impédance d’onde effective η ef f de la struc-
ture de couche par :

(1 + P )
η ef f = η (2.29)
(1 − P ) 0

Les paramètres effectifs nef f et η ef f peuvent être écrit sous la forme :

c √
nef f = arcsin α+ α− (2.30)
ωa

α+
η ef f = √ η0 (2.31)
α+ α−
Où :

2
[sin (k0 a0 + k1 a1 ) − r01 sin (k0 a0 − k1 a1 ) ± 2r01 sin (k1 a1 )]
α± = 2
(2.32)
(1 − r01 )

La courbe dans cette Figure a été calculée pour l’indice de réfraction donné par l’équa-
tion (2.28). L’indice de réfraction effective est purement réel à l’intérieur des bandes. À
l’intérieur des bandes permises, l’onde se propage librement. La vitesse de groupe est
plus grande dans la région de bande et approche zéro aux bords de la bande où dnef f /dω
va à l’infini. Ceci mène à une réflexion complète de l’onde aux bords de la bande.

44
 Figure (2.8) : L’indice de réfraction effective dans le cas de deux
diélectriques parfaits de la même épaisseur avec n0 = 1 et n1 = 4.

La Figure (2.9) montre l’évolution de l’impédance d’onde effective en fonction de la

fréquence. L’impédance d’onde est 1 ou infinie aux bords de la bande donc il y a une
réflexion complète de l’onde aux bords de la bande. Il devrait noter que l’impédance
d’onde est uniquement définie par l’équation (2.29).

45
 Figure (2.9) : L’impédance d’onde effective dans le cas de deux
diélectriques parfaits de la même épaisseur avec n0 = 1 et n1 = 4.

Les paramètres matériels effectifs εef f et µef f peuvent maintenant être calculés à partir
des équations [9] :

1
εef f = nef f /η ef f
c (2.33)
1
µef f = nef f η ef f
c

L’allure de la permittivité effective est décrite sur la Figure (2.10). La permittivité élec-
trique effective εef f est zéro au bord inférieur de la bande, augmente par la bande, et
approche l’infini au bord supérieur de la bande.

46
 Figure (2.10) : La permittivité électrique effective dans le cas de deux
diélectriques parfaits de la même épaisseur avec n0 = 1 et n1 = 4.

La Figure suivante montre l’allure de la perméabilité magnétique effective en fonction

de la fréquence. La perméabilité effective a un comportement opposé. Elle commence à
partir de l’infini au bord inférieur de bande et diminue à zéro en passant par la bande.
Ce comportement implique que près du bord inférieur de bande, le composant de champ
électrique de l’onde est grand tandis qu’au bord supérieur de bande le champ magnétique
est augmenté.

47
 Figure (2.11) : La perméabilité magnétique effective dans le cas de deux
diélectriques parfaits de la même épaisseur avec n0 = 1 et n1 = 4.

Dans la limite de long longueur d’onde, k0 a0 et k1 a1 être petit, et α± peut être exprimé
comme série entière en fréquence. Les expressions résultantes pour εef f et µef f sont :

1 ωa 2 (η 20 − η 21 ) n0 n1 a0 a1 (ε0 a0 + 2ε1 a1 )
εef f = (ε0 a0 + ε1 a1 ) 1 +
a c 12η 1 η 0 a2 (ε0 a0 + ε1 a1 )
2 2 2 (2.34)
1 ωa (η 0 − η 1 ) n0 n1 a0 a1 (µ0 a0 + 2µ1 a1 )
µef f = (µ0 a0 + µ1 a1 ) 1 −
a c 12η 1 η 0 a2 (µ0 a0 + µ1 a1 )

48
Les Figures suivantes montrent l’allure de la permittivité électrique effective et la per-
méabilité magnétique effective dans la limite de la longue longueur d’onde. La première
partie du côté droit de l’équation (2.34) sont les moyennes pesées familières, les cor-
rections étant proportionnelles à ω2 . Cette région est plus simple pour interpréter les
résultats données par l’équation (2.33).

Figure (2.12) : La permittivité électrique effective dans la limite de la long

longueur d’onde.

49
 Figure (2.13) : La perméabilité magnétique effective dans la limite de la long
longueur d’onde.

2.2.6 Les structures qui utilisent des gaps photoniques

Quelle que soit la structure photonique étudiée, les variations des champs dépendent des
motifs choisis. On peut classer les structures qui utilisent les gaps photoniques en deux
familles [12] :
- Celles qui utilisent la périodicité des motifs:
Les antennes : Le champ électromagnétique n’est plus piégé dans le plan de masse
arrière, on obtient ainsi moins de lobes secondaires.
Les filtres : La défection dans le plan de masse d’une ligne permet d’éliminer les ondes
de surface en introduisant une impédance de surface importante. Plus le nombre de motif
est grand, plus la bande rejetée est profonde et large.

50
Les structures cristallines artificielles : Le but est d’empêcher la propagation de
l’onde quel que soit l’angle d’incidence de l’onde excitatrice.
- Celles qui utilisent la périodicité des motifs et des ruptures de périodicité :
Les filtres très sélectifs : La rupture de la périodicité permet la transmission de l’onde
dans la bande interdite à une fréquence précise, fonction de cette rupture.
Les guides : Un couloir est crée au milieu du réseau périodique, l’onde le suit comme
dans le cas d’un guide.
Les coupleurs : En ajoutant à deux guides parallèles une jonction commune, associée
avec des éléments actifs, on peut orienter l’onde dans les guides voulus.

2.3 Conclusion
Ce chapitre rappelle brièvement l’état de l’art des matériaux à bandes interdites pho-
toniques périodiques. L’historique des BIP et leurs propriétés ont alors été présentées.
Ces matériaux artificiels ont permis aux chercheurs de remplacer des systèmes très con-
traignants par des structures simples et faciles à réaliser. De plus, ils ont apporté des
innovations dans plusieurs domaines de l’électromagnétisme.
Nous avons prouvé qu’une structure de couche symétrique peut être décrite par un
milieu effectif avec des paramètres matériels bien définis. Cependant, les paramètres
matériels montrent une dépendance très compliquée de fréquence, particulièrement quand
la longueur d’onde est comparable à la période de la structure. Dans cette gamme de
fréquence, l’extraction des paramètres matériels εef f et µef f à partir des résultats op-
tiques de réflexion et de transmission est très difficile. C’est pour cela qu’on a utilisé la
limite de la longue longueur d’onde qui nous offre une région simple pour l’interprétation
des mesures.

51
 Troisième Chapitre

Résolution du Problème et Discussion

52
Chapitre 3

Résolution du Problème et
Discussion

3.1 Introduction
Dans ce chapitre, on étudie les propriétés des ondes propres de la structure, composée
de couches alternées d’un ferromagnétique et d’un diélectrique non-magnétique et située
dans le champ magnétique de l’extérieur, et pour cela nous allons résoudre le problème en
traitant les caractéristiques dispersives des milieux magnétique et diélectrique en prenant
en compte les composantes non nulles de la perméabilité du milieu ferromagnétique sup-
posé isotrope dépendent de la fréquence. Après avoir rappelé les équations de Maxwell
pour le cas d’un diélectrique parfait nous avons traité les composantes des champs trans-
verses électriques “TE” dans les différentes couches. Ensuite nous avons utilisé la matrice
de transfert reliant les composantes tangentielles des champs électriques et magnétiques
au début et à la fin de chaque période, a fin de déterminer la relation de dispersion pour
aboutir au domaine d’existence des ondes de volume et les ondes de surface.
Puis, nous cherchons le coefficient de réflexion pour une multicouche semi-infinie dans le
cas de l’incidence normale.

53
3.2 Résolution du problème

3.2.1 Les multicouches

Définition :

Les multicouches sont des nanostructures artificielles fabriquées par des techniques de
dépôt sous ultravide afin d’obtenir des propriétés nouvelles. Un milieu multicouche péri-
odique résulte de l’assemblage d’un grand nombre de cellule motif. Chaque cellule com-
prend un nombre réduit de couche d’orientations différentes [13].
Les multicouches périodiques ferromagnétiques diélectriques sont composés d’une alter-
nance de couches ferromagnétiques et de couches diélectriques non magnétiques (Figure
(3.1)).

Figure (3.1) : Une multicouche périodique ferromagnétique diélectrique.

Notion de couche mince :

Au cours des deux dernières décennies, le développement des matériaux sous forme de
films a contribué à une explosion des performances de l’électronique professionnelle et
grand public en abaissant notamment le coût des composants par une production de
masse. Ainsi, depuis l’invention du transistor en 1947, on a assisté à l’intégration de
plusieurs milliers de composants (tels que les composants semi-conducteurs par exemple)

54
permettant de banaliser des appareils comme les calculatrices ou les micro-ordinateurs
personnels, avec des capacités de calcul et de mémoire largement supérieures aux premiers
ordinateurs.
Les films peuvent être élaborés à partir d’une gamme quasi-infinie de compositions si
on considère les éléments de base et leurs alliages tels que les matériaux conducteurs
[20], isolants, réfractaires (oxydes, nitrures, carbures) et les polymères entre autres. La
structure des films déposés peut être mono ou multicouches avec des épaisseurs pouvant
varier d’un plan atomique (quelques Angströms) à plusieurs centaines de micromètres.
Leurs propriétés physiques (électriques, magnétiques, optiques, ...) dépendent de leur
microstructure qui est étroitement liée aux paramètres de dépôt.

3.2.2 Le tenseur de perméabilité ferromagnétique

Le tenseur de Polder :

Dans cette section nous allons relier l’induction magnétique et l’intensité de magnétisa-
tion d’un milieu gyromagnétique au champ dans le milieu, on trouvera que la relation
est de la forme tensorielle. Le cas des pertes nulles a été traité par Polder ainsi que le
cas avec pertes.
L’équation de mouvement du vecteur magnétisation peut être dérivée en considérant un
dipôle magnétique élémentaire ayant un moment dipolaire M placé dans un champ mag-
nétique constant (H0 ). Dans les conditions d’équilibre, le vecteur du moment dipolaire
M se suite dans la direction H0 (la direction de x). Nous supposons maintenant que le
dipôle magnétique est incliné par une petite force externe de sorte qu’il fait un angle θ
avec H0 tel qu’il est représenté sur la Figure (3.2). Puisque le seul champ agissant sur
M est H0 , le couple exercé sur M peut également être écrit sous la forme [17] :

1 dM
C= (3.1)
γ dt

55
 Figure (3.2) : Précession du mouvement magnétique.

−
→ − → − →
Les équations C = M ∧ B donnent l’équation du mouvement pour un dipôle simple :

−
→
dM −
→ − →
= γ M ∧ H0 (3.2)
dt
Où γ désigne le rapport gyromagnétique (rapport entre le moment cinétique et le moment
magnétique).
La magnétisation totale est :

Ms = N.M (3.3)

N : est le nombre de rotations non équilibrées par unité de volume.

L’équation (3.2) devient maintenant :

56
 −→
dMs −→ − →
= γ Ms ∧ H0 (3.4)
dt

Dans le cas des micro-ondes, le champ magnétique appliqué se compose du champ mag-
−
→ −→
nétique statique H0 et du champ magnétique radio fréquence h .

−
→ − → − →
H = H0 + h (3.5)

La magnétisation totale comprend la magnétisation statique Ms et la magnétisation

micro-ondes “m”. Les équations ci-dessus s’écrivent sous d’autres formes :
⎡ ⎤ ⎡ ⎤ ⎡ ⎤ ⎡ ⎤
⎢ H0 ⎥ ⎢ hx ⎥ ⎢ Ms ⎥ ⎢ mx ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
H0 = ⎢
⎢ 0 ⎥,
⎥
h= ⎢
⎢ hy ⎥
⎥
, Ms = ⎢
⎢ 0 ⎥,
⎥
m= ⎢
⎢ my ⎥
⎥
(3.6)
⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦
0 hz 0 mz

L’équation (3.4) se décompose donc sous la forme :

dmx
= mz γhy − my γhz
dt
dmy
= −mz γ (H0 + hx ) − hz γ (Ms + mx ) (3.7)
dt
dmz
= my γ (H0 + hx ) − hy γ (Ms + mx )
dt

Dans une première approximation les termes de m et de h sont négligeables respective-

ment par rapport à M et à H0 , on peut alors écrire :

dmx
≈0
dt
dmy
= −mz γH0 − hz γMs (3.8)
dt
dmz
= my γH0 − hy γMs
dt

57
Réécrivant les dernières équations, nous avons :

..
mx ≈ 0
..
my +ω 2H my = ω M ω H hy − ωM hz (3.9)
..
mz +ω 2H mz = ωM ω H hy − ω M hz

Avec ω M = γMs .
On peut définir un tenseur de susceptibilité [χ] qui relie la magnétisation au champ
magnétiques par :

m = [χ] h (3.10)

Avec :
⎡ ⎤
⎢ 0 0 0 ⎥
⎢ ⎥
⎢ ⎥
[χ] = ⎢
⎢ 0 χyy χyz ⎥
⎥
⎣ ⎦
0 χzy χzz

Et :

ωM ωH
χyy = χzz =
ω2H − ω 2
iωM ω
χyz = χzy = 2
ωH − ω2

En remplaçant ω M par sa valeur, on obtient :

Ms ω 2H
χyy = χzz =
H0 ω2H − ω 2 (3.11)
Ms iωω H
χyz = χzy =
H0 ω2H − ω 2

58
Il y a une résonance à ω = ωH .
D’autre part :
Les composantes de l’induction magnétique sont :

by = hy + 4πmy
(3.12)
bz = hz + 4πmz

Ce qui donne, selon la relation (3.11) :

by = hy + 4πχyy + 4πχyz hz
(3.13)
bz = hz + 4πχzz + 4πχzy hy

En posant :

µ = 1 + 4πχzz
(3.14)
ik = 4πχyz

Il vient que :

by = µhy − ikhz
(3.15)
bz = ikhy + µhz

b
Le rapport correspond au tenseur de Polder [17] qui vaut :
h
⎡ ⎤
b ⎢ µ −ik ⎥
=⎣ ⎦ (3.16)
h ik µ

Si nous désirons introduire les composantes suivant la direction Ox, toujours dans le cas
limite de la Figure (3.2), nous alors trouver pour expression du tenseur de perméabilité :

59
 ⎡ ⎤ ⎡ ⎤
⎢ µxx 0 0 ⎥ ⎢ 1 0 0 ⎥
⎢ ⎥ ⎢ ⎥
⎢ ⎥ ⎢ ⎥
[µ] = ⎢
⎢ 0 µyy µzy ⎥
⎥
= ⎢
⎢ 0 µ −ik ⎥
⎥
(3.17)
⎣ ⎦ ⎣ ⎦
0 µyz µzz 0 ik µ

En remplaçant µ et k par leurs valeurs tirées des relations (3.11) et (3.14), il vient :
⎡ ⎤
⎢ 1 0 0 ⎥
⎢ 2 ⎥
⎢
⎢
4πMs ω H 4πMs ωω H ⎥
⎥
[µ] = ⎢ 0 1+ 2
−i ⎥ (3.18)
⎢ H0 ωH − ω 2 H0 ω 2H − ω2 ⎥
⎢ ⎥
⎣ 4πMs ωω H 4πMs ω2H ⎦
0 i 1 +
H0 ω 2H − ω2 H0 ω 2H − ω 2

Dans la direction Ox, la perméabilité est égale à 1, ce qui est normal puisque le matériau
est saturé.
Dans le modèle de Polder, le coefficient d’amortissement α est nul. µ et k sont donc
purement réels [17] :

4πMs ω 2H
µyy = µzz = 1 + (3.19)
H0 ω2H − ω 2
4πMs ωω H
µyz = −µzy = (3.20)
H0 ω 2H − ω2

Et dépendent de la fréquence de l’onde hyperfréquence et du champ appliqué.

Modèles de Landau-Lifshitz et de Polder :

Si en tenant compte l’amortissement, l’équation du mouvement de pression devient :

−→
dMs −→ − → −
→
= γ Ms ∧ H0 − A (3.21)
dt

60
 −
→
Où A représente l’atténuation avec :

−
→ α −→ −
→ − →
A = M ∧ M ∧ H0 (3.22)
|M|

1
Où : α = .
ωτ
De même manière que le cas de A = 0, la résolution de l’équation (3.21) donne le tenseur
de perméabilité de Polder suivant :
⎡ ⎤
   
µ − iµ −i k − ik
[µ] = ⎢
⎣    
⎥
⎦ (3.23)
i k − ik µ − iµ

Le modèle de Landau-Lifshifz tient compte des pertes magnétiques [17].

 ω M (ω + ωH ) ωH τ 2 + 1 (ω − ωH ) ω H τ 2 + 1
µ =1+ +
2ωH (ω + ωH )2 τ 2 + 1 (ω − ωH )2 τ 2 + 1
 ωM ωτ ωτ
µ = +
2ω H (ω + ω H ) τ + 1 (ω − ω H )2 τ 2 + 1
2 2
(3.24)
 ω M (ω + ω H ) ωH τ 2 + 1 (ω − ωH ) ωH τ 2 + 1
k = +
2ωH (ω + ωH )2 τ 2 + 1 (ω − ωH )2 τ 2 + 1
 ωM ωτ ωτ
k = −
2ω H (ω + ω H ) τ + 1 (ω − ω H )2 τ 2 + 1
2 2

L’évolution de ces perméabilités en fonction de la fréquence donnée par les équations

(3.19) et (3.20) est représentée sur les Figures (3.3) et (3.4).

61
 Figure (3.3) : Variation de la partie réelle de la perméabilité magnétique
en fonction de la fréquence ω.

La Figure (3.3) représente l’allure de la partie réelle de µ en fonction de la fréquence

ω. On peut vérifier que la partie réelle de la perméabilité magnétique est strictement
nulle à la résonance, maximal ou minimal environ à une largeur de la résonance. Une
telle courbe est appelée courbe de dispersion, elle décrit la variation de la perméabilité
magnétique au voisinage de la résonance.

62
 Figure (3.4) : Variation de la partie imaginaire de la perméabilité magnétique
en fonction de la fréquence ω.

La Figure (3.4) illustre le comportement de la partie imaginaire de µ en fonction de la

fréquence ω. La partie imaginaire est en revanche maximal pour ω = ω H , elle décroît
rapidement des que ω s’écarte de ω H . Cette courbe est elle est symétrique par rapport à
l’axe passant par son maximum situé à ω = ω H : c’est la courbe d’absorption.

3.2.3 Les équations de Maxwell

Dans un milieu diélectrique parfait où il n’y a ni charges, ni courants, les équations de

Maxwell se simplifient et s’écrivent dans le domaine temporel [19] :

−→
−→− → ∂H −→
Rot E = −µ = iωµ H (3.25)
∂t
63
 −→
−→− → ∂E −→
Rot H = ε = iωε E (3.26)
∂t
−→
div D = 0
−→
div B = 0

−
→ − →
Où E et H désignent respectivement les champs électriques et magnétiques.
Développons les deux premières équations de Maxwell afin d’obtenir les relations entre
les composantes des champs.
Première équation de Maxwell :

−
→ − → −→
∇ × E = iωµ H

∂Ez ∂Ey ∂Ey

− = −iωµHx ⇒ −iβEz − = −iωµHx (3.27)
∂y ∂z ∂z
∂Ex ∂Ez ∂Ex
− = −iωµHy ⇒ = −iωµHy (3.28)
∂z ∂x ∂z
∂Ey ∂Ex
− = −iωµHz ⇒ iβEx = −iωµHz (3.29)
∂x ∂y

La deuxième équation donne :

−
→ − → −→
∇ × H = iωε E

∂Hz ∂Hy ∂Hy

− = iωεEx ⇒ −iβHz − = iωεEx (3.30)
∂y ∂z ∂z
∂Hx ∂Hz ∂Hx
− = iωεEy ⇒ = iωεEy (3.31)
∂z ∂x ∂z

64
 ∂Hy ∂Hx
− = iωεEz ⇒ iβHx = iωεEz (3.32)
∂x ∂y

En examinant ces équations, on remarque qu’il existe deux types de solutions indépen-
dantes :
1- Modes Transverse Magnétique : Les équations (3.27), (3.31) et (3.32) relient
entre elles les composantes Hx , Ey et Ez .
Le champ magnétique Hx est perpendiculaire au plan d’incidence, il ne possède pas de
composante Hz dans la direction de propagation. Le champ magnétique est purement
transverse, alors que le champ électrique qui est parallèle au plan d’incidence possède
une composante Ez , pour cette raison le mode est dit transverse magnétique (ou TM).
2- Modes Transverse Electrique : Les équations (3.28), (3.29) et (3.30) relient entre
elles les composantes Ex , Hy et Hz . Le champ électrique Ex est perpendiculaire au plan
d’incidence. Le champ magnétique parallèle à ce dernier possède donc une composante
Hy et une composante Hz .
Dans le cas des modes TE. On obtient facilement les expressions suivantes pour les
champs transverses :

1 ∂Ex
(3.28) ⇒ Hy = −
iωµ ∂z
β
(3.29) ⇒ Hz = − Ex
ωµ

En substituant les expressions de Hy et Hz dans l’équation (3.30), on obtient finalement

l’équation de propagation :

∂ 2 Ex
2
+ ω 2 εµ − β 2 Ex = 0 (3.33)
∂z

65
Les composantes des champs électriques et magnétiques des modes TE dans les différentes
couches sont donnés par [14] :

La couche ferromagnétique :
⎧
⎪
⎪
⎪
⎪
⎪
Exf = A1 exp (iν f z) + A2 exp (−iν f z)
⎪
⎪
⎨ i ∂Exf kµyz
Hyf = + Exf (3.34)
⎪
⎪ k0 µ⊥ ∂z k0 µ⊥ µzz
⎪
⎪ kµyy iµzy ∂Exf
⎪
⎪
⎪
⎩ Hzf = − Exf −
k0 µzz µ⊥ k0 µzz µ⊥ ∂z

ω
Où k0 = est le module du vecteur d’onde dans le vide, ω et c sont la fréquence et la
c
vitesse de l’onde dans le vide.

La couche diélectrique :
⎧
⎪
⎪ Exd = B1 exp (iν d z) + B2 exp (−iν d z)
⎪
⎪
⎪
⎪
⎨ i ∂Exd
Hyd = (3.35)
⎪
⎪ k0 µd ∂z
⎪
⎪ k
⎪
⎪
⎩ Hzd = − Exd
k0 µd

Dans ces expressions on a introduit les composantes transversales du vecteur d’onde de

chaque couche qui sont déterminées par les relations :

ν 2f = k02 εf µ⊥ − k2
(3.36)
ν 2d = k02 εd µd − k2

Où la perméabilité effective des couches magnétiques est :

µyz µzy
µ⊥ = µyy − (3.37)
µzz
Sur la Figure suivante (Figure (3.4)), nous avons tracé la propagation du champ électrique
de l’onde électromagnétique TE dans la couche diélectrique.

66
 Figure (3.4) : Représentation tridimensionnelle de l’onde TE
dans la couche diélectrique.

3.2.4 Matrice de transfert et l’équation de dispersion

La méthode des matrices de transfert :

Elle est souvent utilisée pour déterminer les coefficients de réflexion et de transmission des
empilements unidimensionnels, mais a adaptée au cas bidimensionnel [17]. La structure
est divisée en une succession de couches et une matrice transfert permet de relier les
champs dans une couche à ceux de la couche précédente.

L’équation de dispersion :

Pour obtenir la distribution de champ dans la structure périodique multicouches on

introduit la matrice de transfert M qui relie les composantes tangentielles du champ
électromagnétique de la première interface à celles de la dernière interface.

67
Autrement dit, la relation entre les amplitudes des champs sur la frontière de la nième
couche avec des amplitudes sur la frontière de la (n − 2) ième couche est [14] :

Exn = m11 Exn−2 + m12 Hyn−2

(3.38)
Hyn = m21 Exn−2 + m22 Hyn−2

Cette équation s’écrit sous forme matricielle :

⎛ ⎞ ⎛ ⎞⎛ ⎞
⎜ Exn ⎟ ⎜ m11 m12 ⎟ ⎜ Exn−2 ⎟
⎝ ⎠=⎝ ⎠⎝ ⎠ (3.39)
Hyn m21 m22 Hyn−2

La matrice 2 × 2 est appelée matrice de transfert de la couche, elle s’écrit :

⎛ ⎞
m11 m12 ⎟
M =⎜
⎝ ⎠ (3.40)
m21 m22

Les conditions de continuité et de périodicité des composantes tangentielles sont :

Exf (d1 ) = Exd (d1 ) , Hyf (d1 ) = Hyd (d1 )

(3.41)
Exd (d) = Exf (0) exp (iν ef d) , Hyd (d) = Hyf (0) exp (iν ef d)

Pour les composantes de la matrice de transfert on obtient l’expression suivante :

kµyz µ νd
m11 = C1 C2 − i S1 C2 − ⊥ S1 S2
ν f µzz µd ν f
kµyz µd k2 µ2yz
m22 = C1 C2 − i S1 C2 + − ν f S1 S2
ν f µzz ν d µf ν f µ2zz
kµyz k0 µd (3.42)
k0 µd k0 µ⊥
m12 =− C1 S2 + S1 S2 − S1 C2
iν d ν f µzz ν d iν f
iν d kµyz ν d 1 k2 µ2yz
m21 = C1 S2 − S1 S2 − i − ν f S1 C2
k0 µd ν f µzz k0 µd k0 µ⊥ ν f µ2zz

68
Avec :

C1 = cos ν f d1 , C2 = cos ν d d2 , S1 = sin ν f d1 , S2 = sin ν d d2

En utilisant les deux premières égalités on obtient la relation de dispersion des ondes
électromagnétiques dans la multicouche infinie :

m11 + m22 µ µ2⊥ ν d ν f k2 µ2yz

cos dν ef = = C1 C2 − d + − S1 S2 (3.43)
2 2µ⊥ µ2d ν f ν d ν f ν d µ2zz

Où ν ef est la composante transversale du vecteur d’onde de l’onde qui se propage dans

la structure et joue le rôle du nombre d’onde de Bloch. La partie imaginaire du nombre
d’onde de Bloch détermine la profondeur de pénétration effective du champ de fréquence
ultra-haute dans la structure δ ef = (Im ν ef )−1 .
Les relations k (ω) sont :

ω√
k1 = εd µd
c (3.44)
ω
k2 = εf µ⊥ (ω)
c

L’ensemble des résultats des équations (3.44) est représenté sur la Figure (3.5) où nous
avons tracé la relation de dispersion. La courbe (1) représente la relation de dispersion
dans le matériau diélectrique (Polyester), les deux branches (courbe (2)) correspondent
au matériau ferromagnétique (Fer).

69
 Figure (3.5) : La relation de dispersion.

3.2.5 Les ondes de volume et les ondes de surface

Définition des ondes de surface :

Les polaritons de surface apparaissent à l’interface entre deux milieux. Ils représen-
tent des solutions particulières des équations de Maxwell qui correspondent a des ondes
se propageant parallèlement a l’interface et dont l’amplitude décroît exponentiellement
lorsque l’on s’éloigne perpendiculairement a celle-ci. C’est pourquoi ces ondes sont sou-
vent appelées ondes de surface : elles restent confinées à proximité de l’interface [20].

Définition des ondes de volume :

Les ondes de volume existent également dans les couches diélectriques et les couches fer-
romagnétiques.

70
Les domaines des ondes de volume et de surface :

La Figure (3.6) représente les dépendances de la fréquence du nombre d’onde normé

(courbes 1, 2) obtenues par les relations (3.44).

Figure (3.6) : Domaines des ondes de volume et de surface en fonction

de la fréquence et de la constante de propagation.

Les domaines où les paramètres ν d et ν f sont réels (de sorte qu’il y a des ondes de vol-
ume) sont hachurés, alors que les domaines où les paramètres ν d et ν f sont imaginaires
(de sorte que les ondes sont des ondes de surface) ne le sont pas. Les ondes de surface
se localisent sur l’interface qui sépare les couches différentes. Le domaine doublement
hachuré marque les ondes de volume appartenant à deux milieux. Le domaine simple-
ment hachuré marque les ondes de volume du côté du milieu correspondant et les ondes
de surface du côté de l’autre milieu. Les fréquences caractéristiques présentées sur la
Figure dépendent essentiellement du champ magnétisant.

71
3.2.6 Le spectre du coefficient de réflexion

Équation de Helmholtz :

Les équations de Maxwell constituent une série d’équations différentielles couplées; il

convient donc d’en déduire les équations séparées pour chacun des vecteurs du champ.
∂
Nous appliquons l’opérateur rotationnel à l’équation (3.25) et l’opérateur à (3.26); en
∂t
comparant les deux équations résultantes nous obtenons [15] :

−→−→− → ∂ −→− →
RotRot E = −µ Rot H (3.45)
∂t

−→
−−→ − → −→ ∂2 E
graddiv E − ∆ E = −εµ 2 (3.46)
∂t
Et, par conséquent, l’équation d’onde pour le champ électrique s’écrit :

−→
−→ ∂2 E
∆ E − εµ 2 = 0 (3.47)
∂t

De façon similaire on trouve pour le champ magnétique :

−→
−→ ∂2 H
∆ H − εµ 2 = 0 (3.48)
∂t

Où :

1
c= √
µε

72
Ces équations sont vérifiées par les ondes monochromatiques planes :

E (z, t) = E0 exp (kz − iωt) (3.49)

H (z, t) = H0 exp (kz − iωt) (3.50)

Dont la pulsation ω et le vecteur d’onde k sont reliés par :

ω
|k| ≡ k =
c

Le front de l’onde est perpendiculaire à k et il se déplace le long de k avec la vitesse c.

L’application des EM à (3.49) et (3.50) donne la condition de transversalité des champs
du rayonnement et une relation entre les amplitudes :

k ⊥ E0 ⊥ H0 ⊥ k (3.51)

E0 µ E0 1
= ≡ η, = =c (3.52)
H0 ε B0 µε

Le coefficient de réflexion :

On va déterminer le coefficient de réflexion pour la structure périodique examinée et

analyser sa dépendance vis à vis de la structure et de la radiation. Supposons une onde
√
plane monochromatique à la fréquence ω et au nombre d’onde ki = k0 εµ propageant
dans la direction z < 0, rencontre la frontière entre la couche 1 et la couche 2. À la
frontière de plan, une partie de l’onde est réfléchie et propage en arrière par la couche 1
et une partie de l’onde est transmis, propageant vers la droite dans la couche 2 (Figure
(3.7)) [16].

73
 Figure (3.7) : Représentation du champ électrique et magnétique pour
l’onde incidente, réfléchie et transmise.

L’onde incidente peut être caractérisée par :

µ (i) −→
EI (z, t) = Hy exp (i (ki z − ωt)) i
ε (3.53)
(i) −→
HI (z, t) = Hy exp (i (ki z − ωt)) j

Comme mentionné, il provoque une onde réfléchie, propageant vers la gauche dans la
couche 1 :

µ (r) −
→
ER (z, t) = − Hy exp (i (−ki z − ωt)) i
ε (3.54)
(r) −
→
HR (z, t) = Hy exp (i (−ki z − ωt)) j

Et une onde transmise qui continue dessus vers la droite dans la couche 2 :


µ (t) −
→
ET (z, t) =  Hy exp (i (kt z − ωt)) i
ε (3.55)
(t) −
→
HT (z, t) = Hy exp (i (kt z − ωt)) j

74
Dans cette situation le champ total d’onde dans le domaine z < 0 est la somme des
champs d’ondes incidente et réfléchie :

µ (i) −
→
Ex = EI (z, t) + ER (z, t) = Hy [exp (iki z − ωt) − r exp (−iki z − ωt)] i
ε (3.56)
−→
Hy = HI (z, t) + HR (z, t) = Hy(i) [exp (iki z − ωt) + r exp (−iki z − ωt)] j

Le coefficient complexe de réflexion r d’une structure périodique correspond au rapport

entre les amplitudes de l’onde réfléchie et de l’onde incidente :

Hy(r)
r= (i)
Hy

Le coefficient de transmission t est le rapport entre les amplitudes de l’onde transmise et

de l’onde incidente :

Hy(t)
t= (i)
Hy

La Figure (3.8) présente la propagation du champ électrique de l’onde électromagnétique

dans le domaine non magnétique.

75
 Figure (3.8) : Représentation tridimensionnelle de l’onde plane.

Pour obtenir le coefficient de réflexion, on va utiliser les expressions pour les champs dans
chacun des milieux, les conditions aux limites et les conditions de périodicité (3.42), et
en plus les conditions aux limites à la surface qui sépare le milieu homogène du ferro-
magnétique :

Ex (0) = Exf (0), Hy (0) = Hyf (0) (3.57)

A l’aide des expressions (3.41), (3.57) on obtient le système des équations dont la solution
conduit à l’expression du coefficient de réflexion R = |r|2 :

2
exp(iν ef d) − m22 − m12 ε/µ
R= (3.58)
exp(iν ef d) − m22 + m12 ε/µ

76
Il est nécessaire de poser k = 0 en conformité de l’incidence normale de l’onde sur
la structure dans les éléments m12 et m22 de la matrice de transfert de la structure
périodique étudiée.
La transmittance de la structure est égale à :

T =1−R (3.59)

La Figure (3.9) montre la variation des coefficients de réflexion et de transmission en

fonction de la fréquence ω dans la couche diélectrique. Le matériau diélectrique utilisé
dans ce cas c’est le “Polyester”.

Figure (3.9) : Le spectre du coefficient de réflexion et de transmission

pour la couche diélectrique le cas de Polyester.

Le coefficient de réflexion prend la valeur R = 0.53 dans les bandes interdites, alors que
pour les bandes permises des ondes de volume de diélectrique, le coefficient de réflexion

77
R tend vers 0 et l’énergie de l’onde incidente pénètre presque totalement dans la multi-
couche. En négligeant la relaxation, la proportion qui pénètre est égale au coefficient de
transmission T = 1 − R.

3.2.7 Les applications des multicouches

Les guides d’ondes :

Les multicouches métalliques, ces systèmes sont obtenus par dépôt alterné de couches
métalliques magnétiques et de couches isolantes. Les épaisseurs sont typiquement de 1
à 2 nm pour les couches métalliques et 3 à 4 nm pour les couches isolantes. Du fait du
mauvais mouillage du métal sur l’isolant, les couches métalliques se brisent spontanément
en îlots. La structure résultante est alors formée de plans de particules métalliques de
taille nanométrique, séparées dans les plans et entre plans par la matrice isolante. Les
multicouches est à la base de nombreux dispositifs micro-ondes.

Les miroirs ou filtres :

Les multicouches produites dans les années 1980 consistaient à obtenir une périodicité
nanométrique à une dimension de façon à réaliser des miroir à rayons X ou à neutrons.
Ces multicouches amorphes ou épitaxiées font parties de ce qu’on appelle les nouveaux
matériaux terme ayant précédé celui de nanotechnologies.
Pour s’en faire une idée il suffit de s’imaginer un millefeuille tel que chaque couche a une
épaisseur de quelques couches atomiques : on a alors une périodicité nanométrique à une
dimension. Ces matériaux existent naturellement: citons les stéarates ou les mica. Mais
alors le contraste en ce qui concerne les Rayons X.

Les super réseaux :

En 1985, c’est dans le domaine de la supraconductivité à haute température que les

multicouches épitaxiées (en anglais : superlattices) ont débouché.

78
La spintronique :

La magnétorésistance géante a été découverte dans des couches monocristallines en 1988

par deux équipes indépendantes : celle d’Albert Fert, de l’Université de Paris Sud-Orsay,
et celle menée par Peter Grünberg du Centre de recherche de Jülich (Allemagne). La
magnétorésistance géante (en anglais, Giant Magnetoresistance Effect ou GMR) est un
effet quantique observé dans les structures de films minces composées d’une alternance de
couches ferromagnétiques et de couches non magnétiques. Elle se manifeste sous forme
d’une baisse significative de la résistance observée sous l’application d’un champ mag-
nétique externe : à champ nul, les deux couches ferromagnétiques adjacentes ont une
aimantation antiparallèle car elles subissent un couplage ferromagnétique faible. Sous
l’effet d’un champ magnétique externe, les aimantations respectives des deux couches
s’alignent et la résistance de la multicouche décroît brutalement. Le spin des électrons
du métal non magnétique se répartissent équitablement de façon parallèle et antiparal-
lèle, et subissent ainsi une diffusion magnétique moins importante lorsque les couches
ferromagnétiques sont aimantées de façon parallèle.

Cristaux photoniques :

Les Cristaux photoniques sont des structures artificielles périodiques diélectriques à 1,

2 ou 3 dimensions. À une dimension, elles constituent ce que l’on appelait des filtres
interférentiels.

Traitement multicouche des lentilles optiques :

Un des problèmes dans l’utilisation des lentilles optiques, tant en lunetterie qu’en pho-
tographie, est qu’une partie de la lumière qui devrait traverser la lentille est réfléchie.
Cela a deux conséquences néfastes :
- Cette lumière est perdue (et si le nombre de lentilles à traverser est important, il peut
y avoir une quasi-extinction du faisceau lumineux).
- Cette lumière est renvoyée, ce qui est générateur de reflets parasites.

79
Une solution simple à ce problème est d’utiliser un traitement de surface sur la lentille,
constitué d’une couche transparente d’indice de réfraction différent de la lentille, et d’une
épaisseur correspondant à 1/2 longueur d’onde de la lumière incidente.
De cette manière, la lumière réfléchie arrive en opposition de phase avec la lumière inci-
dente, et les deux s’annulent.

3.3 Conclusion
Rappelons les points importants :
- Pour une structure multicouche périodique ferromagnétique diélectrique la perméabilité
magnétique de la couche ferromagnétique est une grandeur tensorielle et la permittivité
électrique de la couche diélectrique est une grandeur scalaire.
- La résolution des équations de Maxwell pour chaque couche nous conduit à deux ondes
propres, une onde de type TE et une onde de type TM, dans ce cas on a fait la restriction
de notre étude des ondes TE.
- On utilise la matrice de transfert reliant les composantes tangentielles du champ pour
déterminer la relation de dispersion, qui nous a permis de préciser le domaine d’existence
des ondes de volume et de surface.
- Pour déterminer le coefficient de réflexion et de transmission, on fait appel à la propa-
gation de l’onde plane dans le milieu périodique en utilisant les conditions aux limites et
les conditions de périodicité.

80
 Quatrième Chapitre

Comparaisons et Résultats

81
Chapitre 4

Comparaisons et Résultats

4.1 Introduction
Dans notre travail nous avons exposé les propriétés dispersives de la propagation des
ondes électromagnétiques dans les différents milieux diélectriques ferromagnétiques et les
matériaux à Bande Interdite Photonique.
Après avoir fait une comparaison entre les deux types de propagation, nous avons déduit
l’influence du changement des couches et le rôle du champ magnétique extérieur.

4.2 Comparaison entre les deux types de propaga-

tion dans les structures multicouches

4.2.1 La propagation des ondes

Multicouche ferromagnétique diélectrique :

Dans le premier cas, sous l’effet d’un champ magnétique extérieur la résolution des équa-
tions de Maxwell pour chaque couche montre qu’il existe deux types d’ondes électromag-
nétiques peuvent se propager dans le milieu périodique :
Transverse Electrique “TE” lorsque le champ électrique est parallèle à l’axe des tiges.

82
Transverse Magnétique “TM” lorsque le champ magnétique est parallèle à l’axe des tiges.

Figure (4.4) : Polarisations possible dans une structure multicouche

périodique : TE et TM.

Puisque cette étude est très compliquée nous considérons le cas d’une onde polarisée
“TE”.

Les matériaux à Bande Interdite Photonique :

Dans le deuxième cas, nous avons étudiée la propagation de l’onde TEM dans la structure
périodique BIP. Nous avons calculé les paramètres effectifs du milieu matériel.
La propagation des ondes est permises dans certaines régions et interdites dans d’autres
régions.

4.2.2 Le coefficient de réflexion et transmission

Multicouche ferromagnétique diélectrique :

L’étude d’une structure multicouche périodique formée par une alternance de couches
ferromagnétiques et diélectriques se fait en calculant le coefficient de réflexion et trans-

83
mission dans une direction éclairée par une onde plane.
La variation fréquentielle du coefficient de réflexion et de transmission obtenu présente
plusieurs zones (Figure (4.1)) illustrant les propriétés de la propagation des ondes dans
les bandes permises.
La méthode des ondes planes a été adaptée pour calculer le coefficient de réflexion et de
transmission des matériaux ferromagnétiques diélectriques périodiques.

Figure (4.1) : Le spectre du coefficient de réflexion et de transmission

pour la couche diélectrique.

Les matériaux à Bande Interdite Photonique :

Les matériaux à Bande Interdite Photonique sont des structures diélectriques ou mé-
talliques composées d’un assemblage périodique de ou plusieurs matériaux dont la péri-
odicité peut exister dans une, deux ou les trois dimensions de l’espace. Les matériaux
élémentaires composant cette structure peuvent quant à eux avoir la forme de plaque
(1D).
La propriété fondamentale d’un matériau à Bande Interdite Photonique est de réfléchir

84
les fréquences comprises dans la bande et d’être transparent aux autres fréquences.
Dans ce cas la matrice de réflexion est donnée par :
⎡ ⎤
1 ⎢ 1 rνµ ⎥
Rνµ = ⎣ ⎦ (4.1)
1 + rνµ rνµ 1

Avec :

1 − η ν /η µ
rνµ = (4.2)
1 + η ν /η µ

4.2.3 La matrice de transfert

Dans les deux cas la matrice de transfert est utilisée pour calculer la relation de dispersion.

Multicouche ferromagnétique diélectrique :

Dans le premier cas la matrice de transfert relie les composantes tangentielles du champ
électrique et magnétique au début et à la fin de chaque période.
Les composantes de la matrice de transfert sont :

kµyz µ νd
m11 = C1 C2 − i S1 C2 − ⊥ S1 S2
ν f µzz µd ν f
kµyz µd k 2 µ2yz
m22 = C1 C2 − i S1 C2 + − ν f S1 S2
ν f µzz ν d µf ν f µ2zz
kµyz k0 µd (4.3)
k0 µd k0 µ⊥
m12 =− C1 S2 + S1 S2 − S1 C2
iν d ν f µzz ν d iν f
iν d kµyz ν d 1 k2 µ2yz
m21 = C1 S2 − S1 S2 − i − ν f S1 C2
k0 µd ν f µzz k0 µd k0 µ⊥ ν f µ2zz

85
Les matériaux à Bande Interdite Photonique :

Dans le deuxième cas la matrice de transfert relie les amplitudes d’ondes et elle est donnée
par le produit de la matrice de propagation dans la première région avec la matrice de
réflexion avec la matrice de propagation dans la deuxième région :

A0 = P0 (z0 ) R01 P1−1 (z0 ) A1 (4.4)

Nous déterminons la matrice de transfert pour chaque période dans la structure avec un
nombre arbitraire de couches.

Figure (4.2) : Visualisation de structure BIP 1D.

La matrice de transfert de période pour notre cas sera alors trouvée par le choix (s = 2).
Dans la couche 1 par exemple la solution de l’équation de Helmholtz peut être écrite

86
comme somme d’une onde propageant vers la droite et d’une onde propageant vers la
gauche (Figure (4.2)).

E (z) = A+ eikz + A− e−ikz (4.5)

4.2.4 Les propriétés dispersives

Multicouche ferromagnétique diélectrique :

Dans le premier cas à partir des composantes m11 et m22 , nous avons calculé la relation
de dispersion donnée par :

m11 + m22 µ µ2⊥ ν d ν f k2 µ2yz

cos dν ef = = C1 C2 − d + − S1 S2 (4.6)
2 2µ⊥ µ2d ν f ν d ν f ν d µ2zz

Suivant les valeurs de k (ω), nous avons tracé le diagramme de dispersion illustrant la
variation de la fréquence en fonction du nombre d’onde, pour aboutir aux ondes de volume
et de surface.

Les matériaux à Bande Interdite Photonique :

Dans le deuxième cas à partir de la relation de dispersion donnée par :

1 2
cos θ = 2
cos (k0 a0 + k1 a1 ) − r01 cos (k0 a0 − k1 a1 ) (4.7)
(1 − r01 )

Avec :

(1 − η 0 /η 1 )
r01 = (4.8)
(1 + η 0 /η 1 )

87
Et :

ω
k0 = n0
c (4.9)
ω
k1 = n1
c

Nous avons des paramètres qui représentent les bandes permises k1 a1 et les paramètres
qui représente les gaps k0 a0 .
Grâce à cette étude, nous avons représenté le diagramme des bandes qui nous a permis
de mettre en évidence l’existence des bandes interdites et des bandes permises.

Figure (4.3) : Le diagramme de dispersion.

Cette équation présente la particularité de n’avoir de solution que lorsqu’elle est comprise
entre −1 et +1. Comme k0 et k1 dépendent tous deux de la pulsation ω, on parle alors de
bandes de fréquences interdites (donc aucune onde électromagnétique ne se propagera).
Le matériau périodique unidimensionnel empêche les ondes électromagnétiques de se
propager à ces fréquences considérées.

88
4.3 Conclusion
La propagation des ondes électromagnétiques dans les structures périodiques multi-
couches a fait l’objet de nombreuses études au cours des dernières années. Les matériaux
magnétiques ouvrent des perspectives particulièrement intéressantes en raison de la pos-
sibilité de transformer leurs caractéristiques en appliquant des champs électriques et
magnétiques.
Les matériaux à gaps photoniques, ou diélectriques artificiels ont été remis au goût du
jour depuis le début des années 90, on ne peut cependant pas parler de nouvelles struc-
tures mais d’un phénomène de mode dans la recherche scientifique. Si ces matériaux à
gaps photoniques semblent très intéressants puisqu’ils permettent une amélioration des
performances de nombreux systèmes existants (antennes, filtres, coupleurs).

89
Conclusion Générale

90
 Conclusion générale

Dans ce travail, nous avons étudié les propriétés dispersives de la propagation des ondes
électromagnétiques dans une structure multicouche périodique.
Dans une première étape, nous avons rappelé les concepts fondamentaux concernant les
modes de propagation des ondes Transverse Electrique “TE” et Transverse Magnétique
“TM”, ainsi que les matériaux magnétiques et diélectriques.
Dans une deuxième étape, nous avons étudiée les matériaux à Bande Interdite Photonique
(BIP) qui ont gagné beaucoup d’intérêt ces dernières années dû à leurs applications poten-
tielles dans l’électromagnétisme et l’électronique. Les structures BIP sont des systèmes
matériels périodique dont les propriétés intéressantes surgissent d’une façon compliquée
des réflexions des multiples internes provoquées par la variation matérielle. Les matéri-
aux BIP ont été les premiers à intéresser les chercheurs qui continuent à démontrer leurs
avantages pour différentes applications. Ils permettent d’introduire les mécanismes de
propagation des ondes électromagnétiques. Une succession de bandes interdites et per-
mises apparaît ainsi dans le diagramme de bande.
Enfin nous avons présentés dans la dernière partie la variation de la perméabilité mag-
nétique qui a une forme tensorielle dans le cas des corps isotropes ferromagnétiques.
L’objet de cette partie est de déterminer les éléments de tenseur en résolvant le système
d’équation de Gilbert qui tient en compte le couplage mutuel des deux domaines adja-
cents soumis à un champ magnétique statique extérieur H0 .
Pour étudier les ondes de surface, nous avons résolu le problème électromagnétique cor-
respondant à la propagation d’une onde à l’interface entre deux milieux qui resterait
confinée sur cette interface et les ondes de volume qui se trouvent dans le volume des
différentes couches. Nous en avons déduit une relation de dispersion qui relie la fréquence
avec la constante de propagation. Cette relation de dispersion a ensuite été étudiée dans
le cas d’une interface Polyester/ferromagnétique avec le Fer.
Malgré, nos efforts déployés pour la réalisation de ce modeste projet : sur l’étude des

91
phénomènes de l’interaction des ondes avec les structures multicouches, demeurent in-
suffisants. Nos ambitions c’est de courir vers l’horizon pour une recherche et une étude
plus approfondies. Il est primordial de prévoir et de penser à des nouveaux horizons pour
satisfaire les besoins de nos futurs chercheurs. Les suggestions qu’on peut apporter pour
la continuité de ce travail sont :
1- Une étude plus approfondie et plus rigoureuse sur la propagation des ondes électro-
magnétiques dans les structures multicouches.
2- Une étude similaire concernant la propagation des ondes électromagnétiques dans les
structures multicouches ferromagnétiques avec deux matériaux magnétiques de même
nature ensuite de différente nature pour voir la différence entre les deux études.
3- Pour étudier la propagation de l’onde dans les structures à Bande Interdite Photonique
on fait l’analogie avec le modèle quantique c’est à dire l’analogie entre les équations de
Maxwell et l’équation de Schrödinger.
4- Un dépistage des sources dans ce domaine.
5- Evaluation de la scientificité de ce projet, son intérêt pour les utilisateurs par ses
garanties scientifiques.

92
Bibliographie

[1] Philippe Robert. “Matériaux de l’électrotechnique”, Traité d’électricité, Presse

Polytechniques Romandes, Troisième édition, 1989.

[2] Miouat Azzouz. “Identification-Optimisation paramétrique et quantification des

pertes dans les matériaux ferromagnétiques doux”, Mémoire de Magister en élec-
trotechnique, Université de Batna, 11/02/2006.

[3] Gérard Hincelin. “Guides micro-ondes”, Electronique B8, Notes de cours, Leçon
n◦ 8 et 9.

[4] Gérard Hincelin. “Guides micro-ondes”, Electronique B8, Notes de cours, Leçon
n◦ 10.

[5] Pierre Tremblay. “Guide d’onde diélectrique à 3 couches”, Gel 21943 / Gel 66398,
Documents de référence H-2004, 2002.

[6] Jean Paul Bailon, Jean Marie Dorlot. “Des Matériaux”, Troisième édition,
Presses Internationales Polytechnique, 2001.

[7] C. Kittel. “Physique de l’état solide”, Dunod, Paris, 7 édition, 1998.

[8] Ducharne Benjamin. “Caractérisation et modélisation de circuits magnétiques :

du global vers le local”, Université Claude Bernard-Lyon, 2003.

93
 [9] H. Vainola, J. Sinkkonen, J. Aho. “Electromagnetic Effective Medium Mate-
rial Parameters of a Periodic Multilayer Structure”, Electron Physics Laboratory,
Helsinki University of Technology, 14 August 2001.

[10] Yannick Merle. “Etude de la dispersion électromagnétique dans les matériaux

périodiques diélectriques bidimensionnels”, Thèse pour obtenir le grade de Docteur
de l’université de Limoges, Soutenue le 20 Novembre 2003.

[11] Abdelwaheb Ourir. “Applications de Matériaux à Bandes Interdites Photoniques

et de métamatériaux en Télécommunications”, Thèse pour obtenir le grade de : Doc-
teur en Sciences, Spécialité Physique, Université de Paris Sud XI d’Orsay, Faculté
des Sciences, Soutenue le 5 décembre 2006.

[12] IEEE Student Branch. “Etat de l’art sur les structures à bande interdite pho-
tonique”, Réunion technique du 22 Janvier 2002.

[13] “Liste des termes, expressions et définitions adoptés et publiés au Journal officiel de
la République française : Termes généraux de la chimie”, En application du décret
n◦ 96-602 du 3 juillet 1996 relatif à l’enrichissement de la langue française.

[14] S.V. Elisséeva, D.I. Sementsov. “Dispersion des ondes électromagnétiques dans
une multicouche périodique dans un champ magnétique extérieur : théorie”, Le
département d’électronique quantique et d’optoélectronique de l’université d’état
d’Oulianovsk, Disponible sur Internet le 18 Avril 2006.

[15] P. Kusel. “Electromagnétisme des milieux continus”, Licence de Physique, Institut

Galilée, Université Paris-Nord, 2001.

[16] Imelda van de Voorde. “On the effective formulation and optimization of 1D
periodic multilayer structures”, Master’s Thesis Applied Mathematics, University of
Twente, Department of Applied Mathematics, Chair Applied Analysis and Mathe-
matical Physics (AAMP), May 2004.

94
[17] Matthieu Valetas. “Couches minces magnétiques pour applications hyper-
fréquences : étude des Samarium-Cobalt et des Néodyme-Fer-Bore par pulvérisation
radiofréquence magnétron”, Thèse pour obtenir le grade de Docteur de l’Université
de Limoges, 27 Novembre 2003.

[18] Marc Zelsmann. “Cristaux photoniques en silicium sur isolant pour le guidage,
le filtrage, l’émission et l’extraction de lumière”, Thèse pour obtenir le grade de
Docteur de l’université Joseph Fourier-Grenoble 1, Soutenue le 07 Novembre 2003.

[19] Benmansour Hayat. “Interaction d’une onde électromagnétique avec une vapeur
atomique”, Mémoire de Magister en physique électronique et modélisation, Univer-
sité de Tlemcen, Faculté des sciences, Département de physique, 2006.

[20] Jean Philippe Mulet. “Modélisation du rayonnement thermique par une approche
électromagnétique”, Thèse de doctorat en sciences, Université de Paris XI Orsay,
Soutenue le 13 Mars 2003.

[21] Richard Feynman. “Le cours de physique de Feynman”, Electromagnétisme 2,

Dunod, Paris, 1999 (Pour la version française).

95