Nous donnons ici un très bref aperçu des méthodes utilisées dans les milieux stratifiés anisotropes. Nous nous limiterons au calcul des matrices de Jones (qui généralisent les coefficients de réflexion et transmission) pour illustrer ces méthodes.
On a vu que dans les milieux anisotropes, la relation entre les champ et déplacement électriques est de la forme $$\begin{array} {ll} \fmat{D} = \boldsymbol{\eps} \fmat{E}, & \boldsymbol{\eps}=\eps_0 \boldsymbol{\eps'} \end{array}$$ où la permittivité diélectrique relative $\boldsymbol{\eps'}$ est un tenseur réel symétrique¹. On peut procéder dans ces milieux comme on l'a fait pour les milieux isotropes à la première page du chapitre. En général, on n'a plus d'invariance sous les rotation autour de l'axe des $z$. Néanmoins, pour éviter des expressions trop compliquées, nous nous limiterons encore aux composantes de Fourier avec $k_y=0$. On obtient dans ce cas l'équation de Berreman $$\dfrac{\partial}{\partial z} \begin{bmatrix} \fff{E}_x \\ Z_0 \fff{H}_y \\ \fff{E}_y \\ Z_0 \fff{H}_x \end{bmatrix} = i \dfrac{\omega}{c} \begin{bmatrix} -\ka_x \dfrac{\eps'_{zx}}{\eps’_{zz}} & 1-\dfrac{\ka_x^2}{\eps’_{zz}} & -\ka_x \dfrac{\eps'_{zy}}{\eps’_{zz}} & 0 \\ \eps’_{xx}-\dfrac{\eps'_{xz}\eps'_{zx}}{\eps’_{zz}} & -\ka_x \dfrac{\eps'_{xz}}{\eps’_{zz}} & \eps’_{xy}-\dfrac{\eps'_{xz}\eps'_{zy}}{\eps’_{zz}} & 0 \\ 0 & 0 & 0 & -1 \\ \dfrac{\eps'_{yz}\eps'_{zx}}{\eps’_{zz}}- \eps'_{yx} & \ka_x \dfrac{\eps'_{yz}}{\eps’_{zz}} & \ka_x^2 -\eps’_{yy}+\dfrac{\eps'_{yz}\eps'_{zy}}{\eps’_{zz}} & 0 \end{bmatrix} \begin{bmatrix} \fff{E}_x \\ Z_0 \fff{H}_y \\ \fff{E}_y \\ Z_0 \fff{H}_x \end{bmatrix} \tag{1}$$ où, rappelons le, $\ka_x=ck_x/\omega$. Comme avant, nous utiliserons pour cette équation la notation $\boldsymbol\Phi’(z)=i\dfrac{\omega}{c}\boldsymbol\Delta\boldsymbol\Phi(z).$
Une fois la solution de l'équation obtenue, on peut, si on le désire, calculer les composantes $z$ des champs en utilisant les formules $$\begin{array} {ll} \fff{E}_z = \dfrac{-1}{\eps'_{zz}}\left(\ka_x Z_0 \fff{H}_y +\eps'_{zx}\fff{E}_x+\eps'_{zy}\fff{E}_y\right) , & Z_0 \fff{H}_z = \ka_x \fff{E}_y. \end{array}$$
¹ Cette propriété entraîne que le tenseur $\boldsymbol{\eps'}$ est diagonalisable. Nous n'utiliserons pas ce fait ici.
La principale différence entre l'équation (1) et son équivalent pour les milieux isotropes est que, sauf coïncidence improbable, il n'y a plus de découplage conduisant à deux polarisations indépendantes. Dès qu'il y a une couche anisotrope, on est obligé de travailler en dimension 4 plutôt qu'en dimension 2 dans l'ensemble du milieu stratifié.
L'autre différence majeure, c'est que les équations sont maintenant suffisamment compliquées pour qu'on ne puisse plus utiliser des expressions algébriques pour les valeurs et vecteurs propres. On est obligé de faire tous les calculs de manière numérique. Il est donc impossible de fournir des formules ou des analyses comme nous l'avons fait pour les milieux isotropes.
Le principe global pour calculer les coefficients de réflexion et de transmission est le même que pour les milieux isotropes, sauf qu'il faut maintenant travailler avec deux polarisations simultanément. Cela signifie que les coefficients vont devenir des matrices 2×2, étant donné qu'il y a deux types de champs incidents à considérer ainsi que deux types de champs réfléchis et transmis. Ces matrices portent le nom de « matrices de Jones ».
Pour simplifier, nous supposerons que les milieux entrant et sortant sont l'air (indice de réfraction $n=1$). On utilisera donc aussi bien dans le milieu sortant que dans le milieu entrant les vecteurs propres à quatre composantes
$$\begin{array} {l} |~p, \pm \gt = \left[\pm \cos\theta_0~~~1~~~0~~~0\right]^T \\ |~s, \pm \gt = \left[0~~~0~~~1~~~\mp \cos\theta_0\right]^T \end{array}$$ où $\theta_0$ est l'angle d'incidence. Lorsqu'on forme la matrice $\mat{M}_0$ qui a ces champs comme colonnes et qu'on l'inverse, on trouve comme lignes de la matrice $\mat{M}_0^{-1}$ les champs $$\begin{array} {l} \lt p, \pm~| = \dfrac{1}{2}~\left[\pm 1/\cos\theta_0~~~1~~~0~~~0\right] \\ \lt s, \pm~| = \dfrac{1}{2}~\left[0~~~0~~~1~~~\mp 1/\cos\theta_0 \right]. \end{array}$$
Les matrices caractéristiques sont maintenant des matrices 4×4. Pour les calculer. nous devons tout d'abord, dans chaque milieu $j$, trouver un système complet de vecteurs propres de la matrice $\boldsymbol{\Delta}_j$. Nous noterons ces vecteurs propres $|~j, k \gt$ et les valeurs propres correspondantes $\lambda_{jk}$ ($1 \le k \le 4$).
En inversant la matrice $\mat{M}_j$ qui a les $|~j, k \gt$ comme colonnes, on trouve les champs $\lt j, l~|$ avec la propriété $$\lt j, l~|~j, k \gt = \delta_{lk}$$ (ces champs sont les lignes de la matrice $\mat{M}_j^{-1}$). Une fois ces vecteurs connus, on peut calculer la matrice² $$\exp\left(i\dfrac{\omega}{c}\boldsymbol\Delta_j z\right) =\sum_{k=1}^4 |~j, k \gt \exp\left(i \frac{\omega}{c}\lambda_{jk}z\right) \lt j, k~|.$$ Le reste du calcul des matrices caractéristiques se fait comme dans le cas isotrope.
² Dans cette formule, la direction de propagation de chacun des $|~j, k \gt$ ne joue pas un rôle particulier. Ceci explique pourquoi on a éliminé les symboles $+$ et $-$ dans la notation pour les vecteurs propres et utilisé un symbole générique $\lambda$ pour les valeurs propres.
Les matrices de Jones sont définies par³ $$\begin{array} {lll} \begin{bmatrix} E_{r \mkern 1mu p} \\ E_{r \mkern 1mu s} \end{bmatrix} = \begin{bmatrix} r_{pp} & r_{ps} \\ r_{sp} & r_{ss} \end{bmatrix} \begin{bmatrix} E_{i \mkern 1mu p} \\ E_{i \mkern 1mu s} \end{bmatrix} & \mathsf{et} & \begin{bmatrix} E_{t \mkern 1mu p} \\ E_{t \mkern 1mu s} \end{bmatrix} \begin{bmatrix} t_{pp} & t_{ps} \\ t_{sp} & t_{ss} \end{bmatrix} \begin{bmatrix} E_{i \mkern 1mu p} \\ E_{i \mkern 1mu s} \end{bmatrix}, \end{array}$$ où $E_{i \mkern 1mu p}$ et $E_{i \mkern 1mu s}$, $E_{r \mkern 1mu p}$ et $E_{r \mkern 1mu s}$ sont les composantes $p$ et $s$ en $z=z_0$ des champs électriques incident et réfléchi, respectivement, alors que $E_{t \mkern 1mu p}$ et $E_{t \mkern 1mu s}$ sont les composantes $p$ et $s$ en $z=z_n$ du champ électrique transmis.
Le calcul de ces matrices peut se faire en deux étapes :
On commence avec un champ incident unitaire de type $p$ venant de la gauche : on aura à l'interface $z=z_0$ $$\boldsymbol\Phi(z_0) = |~p, + \gt + r_{pp} |~p, - \gt + r_{sp} |~s, - \gt = t_{pp} \mkern 1mu\mat{N}\mkern 1mu|~p, + \gt + t_{sp} \mkern 1mu\mat{N}\mkern 1mu|~s, + \gt. \tag{2}$$ où $\mat{N}$ est la matrice caractéristique $\mat{N}_{0\kern 2mu n}$. Ceci permet de calculer $t_{pp}$ et $t_{sp}$ en résolvant le système d'équations qui s'obtient en multipliant (2) par $\lt p, +~|$ et $\lt s, +~|$. Une fois $t_{pp}$ et $t_{sp}$ calculés, on obtient directement $r_{pp}$ et $r_{sp}$ en multipliant (2) cette fois-ci par $\lt p, -~|$ et $\lt s, -~|$.
On peut ensuite recommencer le même processus avec un champ incident de type $s$ pour calculer $r_{ps}$, $r_{ss}$, $t_{ps}$ et $t_{ss}$. Le calcul est en tout point similaire.
³ Pour alléger la notation, nous avons omis les indices $0\mkern 1mu(n+1)$ dans les éléments des matrices de Jones.
Pour le contenu de cette page, je me suis en partie inspiré du chapitre 6 du livre Hiroyuki Fujiwara. Spectroscopic Ellipsometry: Principles and Applications. John Wiley & Sons, Ltd. 2007. Notons que le sujet sera traité de manière beaucoup plus élaborée au chapitre IV.
L'étude de la propagation des ondes électromagnétiques dans les cristaux anisotropes n'est pas triviale et l'exemple suivant illustre bien la chose. Imaginons la situation la plus simple possible : un tenseur $\boldsymbol{\eps'}$ diagonal. L'équation aux valeurs propres pour la matrice $\boldsymbol{\Delta}$ se ramène alors à $$\begin{array} {lll} \lambda^2 = \eps'_{xx}(1-\kappa_x^2/\eps'_{zz}) & \mathsf{ou} & \lambda^2 = \eps'_{yy}-\kappa_x^2. \end{array}$$ Si maintenant on a $$\begin{array} {lll} \eps'_{yy} \lt \kappa_x^2 \lt \eps'_{zz} & \mathsf{ou} & \eps'_{zz} \lt \kappa_x^2 \lt \eps'_{yy}, \end{array}$$ cela signifie qu'on a simultanément deux champs qui se propagent et deux autres qui sont évanescents.