Calcul des indices n et k avec Kramers Kronig
Ce calculateur estime l’indice d’extinction k à partir du coefficient d’absorption, calcule un n ponctuel à incidence normale via la réflectance, puis génère une dispersion n(λ) cohérente avec une relation de type Kramers Kronig sur une bande spectrale choisie.
Courbe spectrale n(λ) et k(λ)
Le graphe trace l’indice réel et l’indice d’extinction sur la plage choisie. La partie réelle est estimée par intégration numérique de type Kramers Kronig à partir d’un profil d’absorption lorentzien calibré sur votre valeur de α.
Guide expert du calcul des indices n et k avec Kramers Kronig
Le couple d’indices optiques n et k est au cœur de l’analyse des matériaux en photonique, en couches minces, en microélectronique et en métrologie spectroscopique. L’indice complexe s’écrit classiquement N = n + ik. La partie réelle n gouverne la vitesse de phase et la réfraction, tandis que la partie imaginaire k quantifie les pertes optiques et l’absorption. Dès que l’on cherche à décrire un matériau réel sur un intervalle spectral large, ces deux grandeurs ne peuvent pas être ajustées indépendamment. Elles sont liées par la causalité physique, et c’est précisément cette contrainte qui conduit aux relations de Kramers Kronig.
Dans la pratique industrielle, le calcul de k peut être amorcé à partir du coefficient d’absorption α, grâce à la formule bien connue k = αλ / 4π, à condition d’utiliser des unités cohérentes. Une fois k évalué à une longueur d’onde donnée, on peut combiner cette information à une mesure de réflectance normale R pour estimer n. Cette approche ponctuelle est utile pour un contrôle rapide de dépôt, un suivi de procédé, ou une première inversion de données. Toutefois, quand on veut reconstruire un spectre cohérent sur plusieurs centaines de nanomètres, il faut imposer une cohérence causale entre l’absorption et la dispersion. C’est là que les relations de Kramers Kronig deviennent essentielles.
Pourquoi Kramers Kronig est indispensable
Les relations de Kramers Kronig relient la partie réelle et la partie imaginaire de la réponse optique d’un matériau via une transformée intégrale de type valeur principale. En termes simples, si un matériau absorbe fortement autour d’une résonance, cette absorption modifie automatiquement sa dispersion. On ne peut donc pas choisir un profil de k(λ) arbitraire sans en déduire un impact sur n(λ). Cette contrainte est fondamentale en ellipsométrie spectroscopique, dans l’identification de paramètres de modèles de Lorentz, de Drude Lorentz ou de Tauc Lorentz, ainsi que dans la validation de jeux de données optiques utilisés dans la simulation d’empilements multicouches.
Le calculateur ci-dessus suit une logique utile pour l’ingénierie. D’abord, il convertit votre donnée d’absorption en indice d’extinction. Ensuite, il estime un indice réel ponctuel à partir de la réflectance mesurée à incidence normale. Enfin, il construit une dispersion simplifiée cohérente de type Kramers Kronig à partir d’un profil d’absorption lorentzien. Cette dernière étape ne remplace pas une inversion ellipsométrique complète, mais elle fournit une approximation physiquement motivée pour visualiser l’impact d’une bande d’absorption sur la dispersion de l’indice.
Les équations à connaître
- Indice complexe : N = n + ik
- Coefficient d’absorption : α = 4πk / λ
- Indice d’extinction : k = αλ / 4π
- Réflectance normale : R = ((n – 1)2 + k2) / ((n + 1)2 + k2)
- Dispersion Kramers Kronig, forme continue : n(ω) = n∞ + (2/π) P∫[Ωk(Ω)/(Ω2 – ω2)] dΩ
La présence du symbole P rappelle qu’il s’agit d’une valeur principale de Cauchy. Numériquement, cela signifie qu’il faut éviter la singularité lorsque la fréquence d’intégration coïncide avec la fréquence d’évaluation. Le calculateur contourne ce point en effectuant une intégration discrète et en excluant l’élément diagonal. Cette méthode est robuste pour de la visualisation, du prototypage et des estimations pédagogiques.
Étapes pratiques de calcul
- Choisir une longueur d’onde de référence λ0, souvent celle d’un laser de contrôle ou d’un point de mesure spectrophotométrique.
- Mesurer ou estimer le coefficient d’absorption α à cette longueur d’onde.
- Calculer k avec la relation k = αλ / 4π.
- Mesurer la réflectance normale R, puis résoudre l’expression de la réflectance pour obtenir n.
- Définir une plage spectrale et une largeur de résonance réaliste pour générer un profil d’absorption.
- Appliquer l’intégrale de Kramers Kronig pour reconstruire une dispersion cohérente de n(λ).
- Comparer le n ponctuel obtenu à partir de R avec la valeur n(λ0) prédite par le modèle KK. Un écart significatif indique souvent qu’il faut élargir la plage spectrale, raffiner le modèle d’oscillateur ou revoir la qualité des données d’entrée.
Interprétation physique des résultats
Quand k augmente, le matériau absorbe davantage la lumière. Cette absorption peut être faible dans les diélectriques transparents du visible, mais elle devient importante dans les semi conducteurs proches du gap, dans les métaux ou dans les matériaux présentant des transitions électroniques intenses. Le comportement de n est plus subtil. Loin d’une résonance, n varie souvent lentement. Au voisinage d’une bande d’absorption, on observe généralement une dispersion dite normale avant la résonance, puis une dispersion anormale près du pic, avec parfois une inflexion marquée.
Cette relation entre absorption et dispersion a des conséquences directes sur le design optique. Un ingénieur qui ajuste un empilement antireflet, un miroir de Bragg ou une couche absorbante pour cellule solaire ne peut pas modifier seulement k. Toute variation de k se répercute sur la phase optique à travers n. C’est pour cette raison que les logiciels de simulation avancés imposent presque toujours des modèles dispersifs causaux plutôt que des tableaux de valeurs arbitraires.
Exemples de valeurs usuelles
Le tableau ci-dessous donne quelques ordres de grandeur souvent cités dans la littérature pour des matériaux bien connus autour de 632,8 nm. Les valeurs peuvent varier avec la pureté, la température, la rugosité, l’orientation cristalline et la méthode de mesure, mais elles offrent un bon point de départ pour situer un résultat calculé.
| Matériau | Longueur d’onde | n approximatif | k approximatif | Commentaire technique |
|---|---|---|---|---|
| Silice fondue | 632,8 nm | 1,457 | ≈ 0 | Diélectrique très transparent dans le visible, dispersion modérée, pertes faibles. |
| Silicium cristallin | 632,8 nm | 3,88 | 0,019 | Indice élevé, absorption non nulle dans le rouge, matériau majeur en microélectronique et photovoltaïque. |
| Or | 632,8 nm | 0,18 | 3,43 | Métal noble, forte composante imaginaire, comportement typique de type Drude dans le visible proche IR. |
| Aluminium | 632,8 nm | 1,44 | 7,39 | Métal très réfléchissant avec pertes élevées, sensible à l’oxyde de surface. |
Ce que disent les statistiques de mesure en pratique
Dans un contexte de laboratoire ou de ligne pilote, l’incertitude sur n et k dépend fortement de la méthode expérimentale. Une mesure simple transmission réflexion à une seule longueur d’onde peut être suffisante pour du tri rapide, mais elle reste moins robuste qu’une ellipsométrie spectroscopique multi angle. Les plages d’erreur ci-dessous correspondent à des ordres de grandeur fréquemment rencontrés dans les rapports techniques et publications d’instrumentation.
| Méthode | Plage spectrale courante | Incertitude typique sur n | Incertitude typique sur k | Usage privilégié |
|---|---|---|---|---|
| Transmission et réflexion simples | UV à NIR | ± 0,01 à ± 0,05 | ± 5 % à ± 20 % | Contrôle rapide, échantillons homogènes, faible nombre de paramètres. |
| Ellipsométrie spectroscopique | 190 à 1700 nm selon instrument | ± 0,001 à ± 0,01 | Très faible à quelques pourcents | Couches minces, matériaux dispersifs, extraction de modèles causaux. |
| FTIR avec inversion optique | MIR à FIR | ± 0,005 à ± 0,03 | ± 2 % à ± 10 % | Phonons, polymères, céramiques, revêtements dans l’infrarouge. |
| Réflectométrie simple angle | Visible à NIR | ± 0,02 à ± 0,10 | Très dépendant du modèle | Procédés rapides et calibration initiale. |
Quand l’approximation ponctuelle est suffisante
Un calcul ponctuel de n et k est souvent suffisant dans quatre cas. Premièrement, lors d’un contrôle de routine où l’on cherche surtout à détecter une dérive de lot. Deuxièmement, pour des matériaux faiblement dispersifs sur une bande étroite. Troisièmement, dans un pré dimensionnement optique où l’on a besoin d’un ordre de grandeur crédible. Quatrièmement, quand les données disponibles sont limitées à un seul laser et à une seule mesure de réflectance.
En revanche, dès que l’on traite des semi conducteurs proches du gap, des oxydes dopés, des métaux, des couches nanostructurées ou des polymères avec bandes vibrationnelles marquées, la reconstruction spectrale par Kramers Kronig devient beaucoup plus pertinente. Elle permet de respecter la causalité et de réduire les solutions non physiques.
Bonnes pratiques pour obtenir des résultats fiables
- Vérifier l’unité de α. Beaucoup d’erreurs viennent d’une confusion entre cm-1 et m-1.
- Éviter d’utiliser une réflectance trop proche de 1 dans un modèle simple si la rugosité de surface ou les interférences de couches sont importantes.
- Choisir une plage spectrale suffisamment large. Un calcul KK tronqué sur une bande trop étroite peut déformer la ligne de base de n.
- Introduire un n∞ réaliste. Cette valeur sert d’ancrage hors résonance et influence la stabilité de l’intégration.
- Comparer toujours les résultats à des bases de données de référence quand elles existent.
- Pour les couches minces, tenir compte de l’épaisseur et des interférences. Un simple modèle de matériau massif peut être insuffisant.
Sources de référence recommandées
Pour aller plus loin, il est utile de confronter vos calculs à des ressources institutionnelles et universitaires. Les pages suivantes sont particulièrement utiles pour relier théorie, mesure et données expérimentales :
- NIST, National Institute of Standards and Technology, référence majeure pour la métrologie et les propriétés optiques de matériaux.
- MIT OpenCourseWare, cours universitaires de haut niveau sur l’électromagnétisme, l’optique et les matériaux.
- UC Berkeley Physics, ressources académiques avancées en physique de la matière et optique.
Limites à garder en tête
Le calculateur présenté ici utilise un profil lorentzien unique pour illustrer une reconstruction de type Kramers Kronig. Dans un matériau réel, la réponse optique peut résulter de plusieurs oscillateurs, d’une contribution libre de type Drude, de transitions interbandes, de phonons optiques, voire de rugosité effective. De plus, la géométrie de mesure peut introduire une anisotropie, des effets de polarisation ou des phénomènes de diffusion qui ne sont pas décrits par un modèle isotrope simple.
Malgré ces limites, cet outil répond très bien à un besoin concret : transformer des données immédiatement accessibles, comme α et R, en un premier jeu de paramètres optiques interprétables, tout en rappelant qu’une estimation crédible de la dispersion doit respecter la causalité. C’est précisément la valeur ajoutée des relations de Kramers Kronig. Elles ne sont pas seulement une élégance mathématique. Elles sont la condition qui relie un spectre d’absorption mesuré à une dispersion physiquement admissible.
En résumé
Le calcul des indices n et k avec Kramers Kronig combine trois niveaux de lecture. Le premier est local, avec le calcul direct de k à partir de l’absorption. Le second est expérimental, avec l’extraction de n depuis la réflectance. Le troisième est global, avec une reconstruction dispersive cohérente sur une bande spectrale. Si vous retenez une idée essentielle, c’est celle ci : l’absorption et la dispersion sont liées. Dès qu’un matériau absorbe, son indice réel ne peut plus être traité comme indépendant.
Cette compréhension change la manière de concevoir un revêtement, d’interpréter une mesure et d’ajuster un modèle. Elle permet aussi de détecter rapidement des résultats aberrants. Si une courbe de k propose une bande d’absorption prononcée sans variation correspondante de n, il y a de fortes chances que le modèle soit incohérent. À l’inverse, lorsqu’un profil spectral respecte qualitativement les contraintes de Kramers Kronig, l’ingénieur dispose d’une base bien plus solide pour la simulation, la fabrication et la validation.