Marches et barrières de potentiel ; effet tunnel
Cours · physique quantique (MP), chapitre 3 · CPGE MP (2e année) · le début du chapitre, la suite sur Adloun
« Ce qui est impossible en mécanique classique n'est que peu probable en mécanique quantique. »
Le chapitre précédent a établi la fonction d'onde d'une particule libre, c'est-à-dire évoluant dans un potentiel uniforme. Nous abordons ici la situation physiquement centrale du programme MP : une particule rencontre une variation brutale du potentiel. Deux modèles élémentaires, la marche et la barrière, suffisent à faire apparaître des comportements sans aucun équivalent classique : réflexion partielle sur un obstacle franchissable, pénétration dans une région interdite (onde évanescente) et, plus spectaculaire encore, la traversée d'une barrière alors même que l'énergie de la particule est inférieure à sa hauteur — l'effet tunnel. Ces résultats gouvernent la radioactivité , le microscope à effet tunnel, les diodes tunnel et la fusion stellaire.
Conformément à l'esprit du programme, l'accent est mis sur l'interprétation des résultats et l'exploitation des formules, plus que sur la lourdeur algébrique du cas de la barrière.
3.1 Position du problème et outils
3.1.1 Potentiels constants par morceaux
On modélise l'action du milieu sur la particule par une énergie potentielle constante par morceaux : le demi-espace ou la zone étudiée est découpé en régions où prend une valeur fixe, avec des discontinuités localisées. Ce n'est qu'une idéalisation — un vrai profil varie continûment — mais elle rend les calculs analytiquement solubles tout en conservant l'essentiel de la physique.
Dans chaque région où (constante), on cherche les états stationnaires , dont la partie spatiale obéit à l'équation de Schrödinger indépendante du temps :
Selon le signe de dans une région, la solution change de nature :
- si , on pose et : les solutions sont des ondes planes (régime propagatif) ;
- si , on pose et : les solutions sont des exponentielles réelles (régime évanescent).
3.1.2 Conditions de raccordement
Aux points de discontinuité de , on ne peut pas résoudre l'équation localement ; on raccorde les solutions des différentes régions. Les conditions, admises, traduisent la régularité imposée par l'équation de Schrödinger.
En tout point où présente une discontinuité finie, la fonction d'onde spatiale et sa dérivée première sont continues :
La continuité de assure une densité de présence bien définie ; celle de découle de ce que ne peut avoir qu'une discontinuité finie (car est fini). Ces deux conditions à chaque interface fournissent exactement le nombre d'équations nécessaires pour déterminer les amplitudes.
Si devient infini (marche ou puits infini, chapitre suivant), la continuité de tombe : seule reste continue et s'annule à la paroi.
3.1.3 Courant de probabilité et coefficients et
Pour quantifier « ce qui passe » et « ce qui rebrousse chemin », on utilise le courant de probabilité. Le programme en admet l'expression pour une particule libre, où il s'interprète comme le produit d'une densité par une vitesse ; on admet ici sa forme générale, valable pour tout état stationnaire, y compris une onde évanescente. Cette sous-section, comme les coefficients et , relève du programme de MP (elle n'est pas au programme de MPI).
Pour une fonction d'onde , le courant de probabilité est
Pour une onde plane , on obtient , produit de la vitesse par la densité : c'est un flux de probabilité. En régime stationnaire à une dimension, est indépendant de (conservation du courant), ce qui traduit la conservation locale de la probabilité.
On envoie une onde incidente d'amplitude depuis . Une partie est réfléchie (amplitude , même région), une partie est transmise (amplitude , région de sortie).
On définit
La conservation du courant impose
et s'interprètent comme les probabilités que la particule soit respectivement réfléchie et transmise.
n'est pas le simple rapport des modules carrés : il faut pondérer par les vitesses (donc les vecteurs d'onde) de part et d'autre. On a (même vitesse), mais lorsque les vecteurs d'onde diffèrent ( à l'entrée, à la sortie).
3.2 Marche de potentiel
La marche de potentiel modélise de nombreuses situations réelles : un électron au voisinage d'une surface métallique (le travail de sortie crée une remontée brutale du potentiel), un neutron abordant la surface d'un noyau (chute abrupte du potentiel nucléaire), ou encore une hétérojonction entre deux semi-conducteurs (discontinuité de la bande de conduction).
On considère une marche montante située en :
3.2.1 Cas : réflexion partielle
Dans la région 1, avec . Dans la région 2, l'onde transmise se propage vers les seulement (pas de source à l'infini) : avec .
Établir et et vérifier .
Démonstration
Les conditions de raccordement en donnent :
En éliminant dans la seconde équation : , d'où
On en déduit
On vérifie .
Le résultat majeur est que dès que , c'est-à-dire dès que , même si . Classiquement, une particule d'énergie suffisante passe toujours () : elle est simplement ralentie. Quantiquement, toute variation brutale du potentiel réfléchit partiellement l'onde, exactement comme une onde lumineuse à l'interface entre deux milieux d'indices différents. Le rapport d'amplitude est d'ailleurs l'analogue du coefficient de réflexion de Fresnel en amplitude (en incidence normale, le vecteur d'onde jouant le rôle de l'indice), et celui du facteur de réflexion en énergie.
Ce phénomène subsiste pour une marche descendante ( passe de à , avec , et l'on remplace par ) : la particule accélère mais est encore partiellement réfléchie. Seule compte la discontinuité, pas son signe.
3.2.2 Cas : onde évanescente et réflexion totale
Si , la région 2 est classiquement interdite. Le vecteur d'onde devient imaginaire pur : posons avec . La solution ne diverge pas, on la retient ; l'autre, , exploserait en et est rejetée. Ainsi :
Montrer que tout en admettant une pénétration non nulle de la fonction d'onde dans la zone interdite.
Démonstration
On reprend le raccordement en remplaçant :
D'où . Ce nombre complexe a pour module : . Toute la probabilité est réfléchie. Le courant transmis est nul car est réelle à une phase globale près : , donc et est bien respecté.
Pourtant, la densité de présence dans la zone interdite n'est pas nulle :
La particule pénètre sur une profondeur caractéristique avant que l'onde ne s'éteigne.
Deux ruptures avec la mécanique classique :
- la particule est certainement réfléchie (), comme classiquement, mais avec un déphasage à la réflexion (le coefficient est complexe) ;
- il existe une probabilité non nulle de trouver la particule dans la région interdite , sur une épaisseur . Classiquement, y serait absurde. C'est cette pénétration qui, appliquée à une barrière d'épaisseur finie, engendre l'effet tunnel.
Capacité numérique : profondeur de pénétration
Estimons la profondeur de pénétration pour un électron face à une marche telle que .
import numpy as np
hbar = 1.054e-34 # J.s
m = 9.109e-31 # masse de l'electron (kg)
eV = 1.602e-19 # 1 eV en joules
def profondeur(dV_eV):
"""Profondeur de penetration delta = 1/(2 kappa) en metres."""
dV = dV_eV * eV
kappa = np.sqrt(2 * m * dV) / hbar
return 1.0 / (2 * kappa)
for dV in (0.1, 1.0, 5.0):
print(f"V0 - E = {dV:4.1f} eV -> delta = {profondeur(dV)*1e9:6.3f} nm")
# V0 - E = 1.0 eV -> delta ~ 0.098 nm : quelques dixiemes de nm, echelle atomique
La pénétration est de l'ordre de l'angström : microscopique, mais suffisante pour qu'une barrière atomiquement fine soit franchie par effet tunnel.
3.3 Barrière de potentiel et effet tunnel
On superpose maintenant deux marches : une barrière rectangulaire de hauteur et de largeur :
On s'intéresse au cas tunnel : classiquement la particule incidente rebrousse toujours chemin (). Quantiquement, l'onde évanescente n'a pas le temps de s'éteindre complètement sur la largeur et « ressort » de l'autre côté avec une amplitude réduite mais non nulle.
3.3.1 Structure des solutions
avec et . Dans la région 2, les deux exponentielles sont conservées car la région est bornée : aucune ne diverge. Il y a deux interfaces ( et ), donc quatre équations de raccordement pour quatre rapports d'amplitude — le calcul est lourd (non exigible), on en donne le résultat.
Pour , le coefficient de transmission vaut
Il décroît lorsque la hauteur ou la largeur de la barrière augmente. Quand , il ne tend pas vers mais vers (voir les exercices résolus) : une particule qui a juste l'énergie du sommet est encore partiellement réfléchie.
3.3.2 Barrière épaisse : la formule à retenir
Méthode : Approximation de barrière épaisse (ou opaque)
Lorsque la barrière est « épaisse », au sens , on a , donc . Le terme est négligeable devant le second et
Le facteur dominant est l'exponentielle : la transmission décroît exponentiellement avec
- la largeur de la barrière ;
- la racine du déficit d'énergie ;
- la racine de la masse : pour un proton (), c'est l'exposant qui est multiplié par , si bien que dans les mêmes conditions (et non ) : si un électron traverse avec , le proton traverse avec . L'effet tunnel est donc essentiellement un phénomène d'objets légers.
Cette dépendance exponentielle explique la sensibilité extrême du microscope à effet tunnel (le courant varie d'un ordre de grandeur quand la pointe s'éloigne d'un seul angström) et l'étalement gigantesque des durées de vie en radioactivité .
Pour un électron de face à une barrière de et de , , le préfacteur vaut et , alors que seul donne : le calcul est détaillé dans les exercices résolus (programme de MP). La loi en , elle, s'établit sans la formule de , par la seule décroissance de l'onde évanescente : c'est l'objet de l'exercice du TD « Hauteur, largeur, masse ».
Capacité numérique : en fonction de la largeur et de l'énergie
Traçons et comparons formule exacte et formule approchée : la décroissance exponentielle apparaît comme une droite en échelle logarithmique.
import numpy as np
import matplotlib.pyplot as plt
hbar = 1.054e-34; m = 9.109e-31; eV = 1.602e-19
def T_exact(E, V0, a):
E, V0 = E*eV, V0*eV
kappa = np.sqrt(2*m*(V0 - E))/hbar
return 1.0/(1 + V0**2/(4*E*(V0 - E)) * np.sinh(kappa*a)**2)
def T_approx(E, V0, a):
E, V0 = E*eV, V0*eV
kappa = np.sqrt(2*m*(V0 - E))/hbar
return 16*E*(V0 - E)/V0**2 * np.exp(-2*kappa*a)
a = np.linspace(0.1e-9, 1.0e-9, 200)
plt.semilogy(a*1e9, T_exact(1.0, 2.0, a), label="exact")
plt.semilogy(a*1e9, T_approx(1.0, 2.0, a), '--', label="approx epaisse")
plt.xlabel("largeur a (nm)"); plt.ylabel("T")
plt.title("Transmission tunnel : E=1 eV, V0=2 eV")
plt.legend(); plt.grid(True, which="both"); plt.show()
# La pente en echelle log vaut -2*kappa/ln(10) : signature de l'effet tunnel.
On observe une droite en échelle logarithmique (loi ) et la coïncidence de l'approximation avec l'exact dès que .
3.3.3 Applications de l'effet tunnel
- Radioactivité (Gamow, 1928). Une particule est piégée dans le noyau par la barrière coulombienne. Sa faible probabilité de traverser par effet tunnel, exponentiellement sensible à l'énergie, explique l'écart colossal des durées de vie (s à ans) : c'est la loi de Geiger--Nuttall.
- Microscope à effet tunnel (STM, 1981). Une pointe conductrice est approchée à quelques angströms d'une surface ; les électrons franchissent le vide (barrière de potentiel) par effet tunnel. Le courant, en , varie d'un facteur par angström, offrant une résolution atomique (prix Nobel 1986 à Binnig et Rohrer).
- Jonction et diode tunnel. Dans une jonction semi-conductrice mince ou une jonction Josephson supraconductrice, les porteurs traversent une fine barrière isolante par effet tunnel, à la base de composants électroniques rapides.
La fusion thermonucléaire dans les étoiles repose elle aussi sur l'effet tunnel : les protons franchissent leur répulsion coulombienne bien qu'ils n'aient pas l'énergie classique suffisante.
tunnel la barrière de potentiel que constitue le vide entre la pointe et la surface. Le courant, en , varie d'un facteur par angström :
en asservissant pendant le balayage, on cartographie la surface atome par atome.inférieure au sommet de la barrière coulombienne, est classiquement piégée entre et ; elle s'échappe par effet tunnel à travers la zone interdite (hachures). La sensibilité exponentielle de à explique des durées de vie
allant de s à ans (loi de Geiger--Nuttall).3.4 Synthèse
- Méthode générale. Dans chaque région où est constant, écrire comme onde plane (si ) ou exponentielle réelle (si ), avec
Raccorder et à chaque interface (discontinuité finie de ). Rejeter toute exponentielle divergente dans une région non bornée.
- Coefficients. et via les courants de probabilité ; toujours .
- Marche, . : réflexion partielle quantique même quand la particule a l'énergie de passer.
- Marche, . (réflexion totale) mais onde évanescente : pénétration sur dans la zone classiquement interdite.
- Barrière, = effet tunnel. : la particule traverse. Barrière épaisse () :
Décroissance exponentielle en , en et en .
- Applications. Radioactivité , microscope à effet tunnel (STM), diodes/jonctions tunnel, fusion stellaire.
3.5 Travaux dirigés
Exercice 1 — Des électrons lents qui entrent dans un métal : la marche que l'on franchit, et qui réfléchit
Dans le modèle de l'électron libre, un électron de conduction d'un métal a une énergie potentielle uniforme , plus basse que dans le vide, où on la prend nulle ; on prend . Un électron du vide, d'énergie , arrive en incidence normale sur la surface plane du métal, en : il aborde une marche descendante, pour (vide) et pour (métal). On cherche l'état stationnaire qui comporte dans le vide une onde incidente et une onde réfléchie , et dans le métal une onde transmise . On prend réel et positif.
1) Exprimer et . Pourquoi n'y a-t-il pas d'onde en dans le métal ? Établir et . Que ferait un électron classique ?
2) Montrer que dans le vide . Où sont les maxima et les minima ? Exprimer l'interfrange et le contraste en fonction de . Ces franges bougent-elles au cours du temps ? Que vaut dans le métal, et pourquoi n'y a-t-il pas de franges ?
3) Application à : calculer , , , l'interfrange et le contraste. Tracer l'allure de .
4) Cas limites. Calculer pour et pour ; montrer que quand . Pour une marche montante de hauteur abordée avec (mêmes formules, avec ), que deviennent et quand ?
5) Symétrie et paradoxe. a) Comparer pour une marche montante abordée avec et pour une marche descendante telle que . b) Dans le premier cas, calculer . Une onde transmise d'amplitude supérieure à celle de l'onde incidente viole-t-elle quelque chose ?
6) La transition réelle entre le vide et le métal s'étale sur environ . À quelle condition le modèle de la marche abrupte est-il légitime ? Le vérifier à et à .
Démonstration (Corrigé)
La stratégie. On écrit la solution dans chaque région, on élimine l'onde qui n'a pas de source, et les deux conditions de continuité en fixent les deux rapports d'amplitude. La réflexion se lit ensuite dans la densité de présence : l'onde incidente et l'onde réfléchie, qui appartiennent au même état, interfèrent.
1) La solution. Dans le vide, ; dans le métal, :
Une onde en se propagerait vers les décroissants : elle viendrait de , où il n'y a aucune source ; on la rejette. Les conditions de continuité en donnent
En reportant dans la seconde : , donc
Ce sont les formules de la marche montante du cours : elles ne supposent rien sur le signe de la marche. Ici , donc et . Un électron classique est accéléré vers le métal par la force , dirigée vers les croissants à la surface : il entre toujours, jamais il ne rebrousse chemin. Quantiquement, l'onde réfléchie a une amplitude non nulle dès que : l'électron peut être renvoyé dans le vide par une marche qui l'attire.
2) La réflexion se voit dans la densité de présence. Avec et réel :
Comme , les maxima, de valeur , sont aux points où , et les minima, de valeur , aux points où , dont : la surface est sur un minimum. La période en est l'interfrange , et le contraste vaut
L'état est stationnaire : ne dépend pas du temps, les franges sont immobiles. Ce ne sont pas deux particules qui interfèrent, mais les deux termes d'une même fonction d'onde. Dans le métal, est uniforme : une onde progressive seule a un module constant, il faut deux ondes pour faire des franges. Contrôle : en , , les deux régions se raccordent.
3) Électron de . , donc
La longueur d'onde dans le vide vaut , d'où l'interfrange ; le contraste vaut . Maxima , minima , plateau dans le métal . Une onde réfléchie d'amplitude moitié de l'onde incidente, des franges contrastées à 83 % : pour un électron lent, la réflexion est massive, alors que rien de classique ne s'oppose à son entrée.
4) Cas limites. Électron rapide, : et , contraste : on retrouve à peu près le comportement classique. En général, pour , , donc
soit ici : la marche n'est qu'une petite perturbation pour une particule très énergétique. Électron lent, : et ; quand , , la réflexion devient totale sur une marche descendante. C'est la « réflexion quantique », observée avec des atomes ultra-froids qui approchent une surface qui les attire. Marche montante, : , donc et ; pour , et déjà. Le raccord est continu avec le cas , où la réflexion est totale.
5a) Symétrie. Marche montante à : et . Marche descendante avec : . Même module, signe opposé : change de signe quand on échange et . C'est la même marche parcourue dans les deux sens (énergie cinétique d'un côté, de l'autre). La réflexion ne dépend que du rapport , pas du sens de la marche, comme la réflexion d'une onde lumineuse en incidence normale ne dépend que du rapport des indices.
5b) Le paradoxe . À sur la marche montante : , et la densité transmise vaut fois celle de l'onde incidente. Rien n'est violé : est une densité de probabilité par unité de longueur, pas un nombre de particules qui passent par seconde. Au-delà de la marche, la particule est plus lente (par de Broglie, , et ) : comme les voitures qui se resserrent sur une portion de route où l'on roule moins vite, elle passe plus de temps par unité de longueur, donc on a plus de chances de l'y trouver. Le bilan exact de ce qui passe et de ce qui revient se fait avec le courant de probabilité, qui n'est au programme qu'en MP (il y fait l'objet d'exercices résolus).
6) Validité. Le modèle abrupt est légitime si l'onde ne « voit » pas la structure de la transition, c'est-à-dire si sa longueur d'onde est grande devant l'épaisseur de celle-ci : à , et le modèle est bon. À , est de l'ordre de l'épaisseur : une transition progressive réfléchit encore moins qu'une marche abrupte, et le calculé n'est qu'un majorant d'ordre de grandeur.
Ce que l'exercice installe. Établir la solution d'une marche, montante ou descendante : une onde par région, rejet de l'onde sans source, deux conditions de continuité, . La réflexion se lit dans la densité de présence : franges immobiles d'interfrange et de contraste , plateau uniforme côté transmis. Elle est d'autant plus forte que la particule est lente devant la marche, ne dépend que de , et une amplitude transmise supérieure à n'a rien de paradoxal : une densité n'est pas un flux.
Exercices du chapitre
- Des électrons lents qui entrent dans un métal : la marche que l'on franchit, et qui réfléchit corrigé
- La bouteille à neutrons ultra-froids : réflexion totale et onde évanescente corrigé
- Le mur qui fuit : où une particule voit-elle une paroi de hauteur finie ? corrigé
- Hauteur, largeur, masse : les trois leviers de l'effet tunnel corrigé
- Lire une simulation : la densité de présence autour d'une barrière corrigé
- Garder un électron dix ans derrière 8 nm d'oxyde : la mémoire flash corrigé