Voir la source de Atome d'hydrogène

{{2autres|la théorie quantique de l'atome d'hydrogène|l'élément chimique|Hydrogène|le corps simple H{{ind|2}}|Dihydrogène}}
{{article général|atome}}

L''''atome d'hydrogène''' est le plus simple de tous les [[atome]]s du [[tableau périodique des éléments|tableau périodique]], étant composé d'un [[proton]] et d'un [[électron]]<ref>On se limite ici à l'[[isotope]] le plus abondant naturellement, qui sera seul considéré sauf exception. Pour mémoire le [[deutérium]] à un noyau constitué d'un proton et d'un neutron, et le [[tritium]] (qui est radioactif) d'un proton et de deux neutrons. L'influence de ces nucléons supplémentaires est relativement faible.</ref>. Il correspond au premier élément de la [[Tableau périodique des éléments|classification périodique]].

La compréhension des interactions au sein de cet atome au moyen de la théorie quantique fut une étape importante qui a notamment permis de développer la théorie des [[atome à N électrons|atomes à N électrons]]. C'est pour comprendre la nature de son [[spectre d'émission]], discret, alors que la théorie classique prévoyait un spectre continu, que [[Niels Bohr]] a introduit en 1913 un premier modèle quantique de l'atome (cf. [[Modèle de Bohr|atome de Bohr]]). L'approfondissement de son étude a permis de valider les théories de la [[physique quantique]] au fur et à mesure des progrès accomplis : d'abord l'ancienne [[théorie des quanta]]<ref>par l'intermédiaire de l'[[Modèle de Bohr|atome de Bohr]]</ref>, ensuite la [[mécanique quantique]] non relativiste de [[Erwin Schrödinger|Schrödinger]], la mécanique quantique relativiste de [[Paul Dirac|Dirac]], et enfin la [[théorie quantique des champs]].

Dans le cadre de la [[mécanique quantique]], l'atome d'hydrogène est un [[problème à deux corps]] [[analyse fonctionnelle (mathématiques)|analytiquement]] soluble, du moins dans le cas d'un modèle non-relativiste d'un hamiltonien où est seule prise en compte l'interaction coulombienne entre le proton et l'électron, considérés comme ponctuels. Il est ainsi possible d'en déduire les [[niveau d'énergie|niveaux d'énergie]], et de les comparer aux mesures des [[raie spectrale|raies spectrales]].

L'étude théorique de l'atome d'hydrogène n'a pas qu'un intérêt purement académique, et limité à ce seul atome: en fait, elle est d'une importance considérable en [[physique atomique|physique atomique et moléculaire]]. Tout d'abord, elle permet de comprendre directement les spectres des [[ion]]s dit ''[[hydrogénoïde]]s'', c'est-à-dire qui ont perdu tous leurs électrons sauf un (par exemple He{{exp|+}}, Li{{exp|2+}}). De façon plus générale les concepts dégagés par l'étude de ce modèle permettent d'expliquer la structure des niveaux d'énergie et les spectres des atomes à plusieurs électrons, qui peuvent être compris dans le cas de modèles à électrons indépendants (approche de champ moyen).

== Données expérimentales et position du problème en mécanique quantique ==
=== Aspects expérimentaux et historique ===
Au cours du {{s-|XIX}}, l'amélioration des techniques de l'optique (notamment avec l'usage des [[réseau de diffraction|réseaux plans]]) permet le développement de la [[spectroscopie]]. Dès 1859 [[Gustav Kirchhoff|Kirchhoff]] et [[Robert Wilhelm Bunsen|Bunsen]] découvrent que les raies spectrales des [[élément chimique|éléments chimiques]] sont caractéristiques de ceux-ci. Cette technique permettra la découverte de nombreux nouveaux éléments parmi lesquels le [[césium]] (en 1860), le [[rubidium]] (en 1861) et surtout l'[[hélium]] en 1868 par [[Jules Janssen|Janssen]] et [[Joseph Norman Lockyer|Lockyer]].

En 1885, [[Johann Jakob Balmer|Balmer]] découvre que les longueurs d'onde des 4 raies d'émission dans le visible de l'hydrogène, identifiées peu avant par [[Anders Jonas Ångström|Angström]] comme étant situées respectivement à {{Unité|656.3|nm}}, {{Unité|486.1|nm}}, {{Unité|434.0|nm}} et {{Unité|410.2|nm}}, sont données par la formule empirique suivante, dite de Balmer :
:<math>\bar{\nu}_n=\frac{1}{\lambda_n}=R_\mathrm{H}\left(\frac{1}{4}-\frac{1}{n^2}\right)</math> avec <math>n>2</math>,
avec <math>R_\mathrm{H}=1,09677 \cdot 10^7 \, \text{m}^{-1}</math> (valeur moderne), dite ''constante de Rydberg'' de l'atome d'hydrogène. Cette formule, dite de Balmer, est rapidement généralisée par [[Johannes Rydberg|Rydberg]] et [[Walther Ritz|Ritz]] à la suite de la découverte de nouvelles « séries » de raies dans d'autres domaines du spectre, sous la forme suivante, dite de Rydberg-Ritz :
:<math>\bar{\nu}_n=\frac{1}{\lambda_n}=R_\mathrm{H}\left(\frac{1}{p^2}-\frac{1}{n^2}\right)</math> avec <math>n > p</math>.
Dans cette formule, dite de Rydberg-Ritz, ''p'' est l'indice de la ''série'', et ''n'' l'indice de la ''raie''. La [[formule de Balmer]] correspond à la série ''p=2'' (cette série est ainsi nommée « [[série de Balmer]] »). Différentes séries convergentes de raies seront progressivement mises en évidence dans le [[spectre de l'atome d'hydrogène]] :
* <math>p=1</math> : [[série de Lyman]], dans l'ultraviolet ;
* <math>p=3</math> : [[série de Paschen]], dans l'infrarouge proche ;
* <math>p=4</math> et <math>p=5</math>, séries de [[série de Brackett|Brackett]] et de [[série de Pfund|Pfund]] dans l'infrarouge lointain.

Parallèlement, la [[formule de Rydberg|formule de Rydberg-Ritz]] est généralisée à d'autres éléments que l'hydrogène, notamment les ions dit ''hydrogénoïdes'' et certains alcalins, au prix d'une modification de la [[constante de Rydberg]], et l'utilisation de nombres non-entiers pour l'indice de raies (« correction de Rydberg<ref>Cf. notamment Herzberg, ''{{Langue|en|Atomic Spectra and Atomic Structure}}'', Dover, 1944.</ref> »).

Dans le même temps, l'existence même de spectres de raies plutôt que de spectres continus n'est pas alors expliquée par la théorie classique, et posera un problème dès le début de l'essor des théories sur la structure de l'[[atome]]. En effet, l'[[expérience de Rutherford]] montre en 1911 que l'atome est constitué d'un noyau positivement chargé, concentrant l'essentiel de la masse dans un rayon environ {{Unité|100000}} fois plus petit que l'atome lui-même, et d'électrons négativement chargés « en orbite » autour de ce noyau sous l'action de l'interaction électrostatique. Si l'origine de l'émission ou l'absorption de radiations électromagnétiques par les électrons peuvent facilement se comprendre dans le cadre de la théorie classique de l'électromagnétisme, une difficulté de taille apparaît. En effet, la théorie prévoit naturellement que les spectres doivent être ''continus'' et non pas de raies. Elle est bien sûr par suite impuissante à expliquer les régularités dans certains spectres mis en évidence par les formules de Rydberg ou de Rydberg-Ritz.

Au demeurant, l'existence même de l'atome ne peut être expliquée par la théorie classique : l'électron étant accéléré dans le champ du noyau, les [[équations de Maxwell]] prédisent qu'il doit rayonner de l'énergie et finir par « s'écraser » sur le noyau.

En 1913, [[Niels Bohr]] propose un modèle empirique permettant de rendre compte de la stabilité de l'atome d'hydrogène et de l'existence des raies de son spectre. Dans ce modèle (cf. [[Modèle de Bohr|atome de Bohr]] pour plus de détails), l'électron est en interaction électrostatique avec le noyau, supposé infiniment lourd, et se déplace sur des orbites circulaires de rayon ''r'' similaires à celles des planètes autour du Soleil (cf. [[problème à deux corps]]). Toutefois Bohr, s'inspirant de la théorie naissante dites des [[théorie des quanta|quanta]], postule que les seules orbites permises sont celles pour lesquelles le [[moment cinétique]] de l'[[électron]] sont quantifiées, avec la condition suivante : <math>m_e r v=n\hbar</math> (''m{{ind|e}}'' étant la masse de l'électron) avec ''n'' entier strictement positif. Sur de telles orbites il est postulé que l'électron ne rayonne pas d'énergie. Les différentes orbites admissibles ont donc des rayons ''quantifiées'' donnés par <math>r_n = n^2a_0</math> avec <math>n=1,2,\ldots,+\infty</math> et <math>a_0=\frac{4\pi\varepsilon_0\hbar^2}{m_e e^2}</math> ''rayon de Bohr'' de l'atome d'hydrogène, de valeur environ 53 pm.

Cette condition implique une ''quantification de l'énergie'' de l'électron qui s'écrit alors sous la forme <math>E_n=-\frac{E_1}{n^2}</math> avec <math>E_1 = \frac{m_e e^4}{32\pi^2\varepsilon_0^2\hbar^2}</math> énergie d'ionisation de l'atome (''e'' est la [[charge élémentaire]]).

Dans ce contexte, l'émission ou l'absorption de lumière est expliquée par Bohr comme le passage d'une orbite d'énergie supérieure (resp. inférieure) ''E{{ind|p}}'' à une autre d'énergie inférieure ''E{{ind|n}}'' avec ''n>p''(resp. supérieure), la [[longueur d'onde]] de la lumière émise étant alors donnée par la formule d'Einstein comme liée à la différence d'énergie entre les deux orbites :
:<math>\frac{hc}{\lambda}=E_1\left(\frac{1}{p^2}-\frac{1}{n^2}\right)</math>, d'où pour le [[nombre d'onde]] de la transition <math>\frac{1}{\lambda}=\frac{E_1}{hc}\left(\frac{1}{p^2}-\frac{1}{n^2}\right)</math>.

Cette théorie permet ainsi de démontrer la formule de Rydberg-Ritz précédente, en identifiant <math>R_\mathrm{H}</math> et <math>\frac{E_1}{hc}</math>, ce qui est vérifié par le calcul qui donne en effet <math>E_1=13,6 \,\mathrm{eV}</math> avec la formule précédente et bien sûr <math>E_1=hcR_\mathrm{H}</math>.

Si la théorie de Bohr permet d'expliquer la formule de Rydberg-Ritz, et l'existence de spectres de raies pour les atomes, elle le fait au prix d'une hypothèse ''ad hoc'' difficilement conciliable avec la théorie classique. Si l'[[expérience de Franck et Hertz]] apporte dès 1914 une confirmation expérimentale du modèle de Bohr, et malgré l'amélioration de celui-ci par [[Arnold Sommerfeld|Sommerfeld]] pour tenir compte d'orbites elliptiques, et expliquer en partie l'existence de la [[structure fine]] du spectre de l'atome d'hydrogène mise en évidence expérimentalement, les difficultés conceptuelles demeurent et ne seront résolues qu'avec le développement de la [[mécanique quantique]] dans les années qui suivront. L'atome d'hydrogène est aujourd'hui celui pour lequel la théorie quantique arrive à décrire le spectre avec la plus grande précision, en plein accord avec l’expérience. Son étude est indispensable pour aborder la théorie des atomes à plusieurs électrons ainsi que celles des [[molécule]]s, en permettant d'introduire de nombreux concepts fondamentaux en [[physique atomique]] et en [[chimie]] (cf [[atomistique]]), en particulier la notion centrale d'[[orbitale atomique]].

=== Description générale de l'atome d'hydrogène ===
L'atome d'hydrogène est constitué d'un noyau comprenant un seul proton de masse <math>m_\mathrm{p}=1,67265
 \cdot 10^{-27}\,\mathrm{kg} </math> et de charge <math>q=+e</math>, et d'un électron de masse <math>m_\mathrm{e}=</math> {{unité|9,10953 e−31 kg}} et de charge <math>q_\mathrm{e}=-e</math>, <math>e</math> étant la charge élémentaire (<math>e=</math> {{unité|1,602189 e−19 C}}). Le proton comme l'électron sont des [[fermion]]s de [[spin]]s respectifs notés <math>I</math> (spin nucléaire) et <math>s</math> (spin électronique), de même valeur <math>I=s=\tfrac12</math>.

=== Position du problème ===
Selon la première approximation, considérée par Schrödinger, l'atome d'hydrogène est étudié dans le cadre d'un modèle quantique exactement soluble, dans lequel le proton et l'électron, considérés comme ponctuels, interagissent uniquement par l'intermédiaire de l'interaction électrostatique de Coulomb, l'énergie potentielle du système étant : <math>V(r)=\frac{q_p q_e}{4\pi\varepsilon_0 \; r}=-\frac{e^2}{4\pi\varepsilon_0 \; r}</math>, avec <math>r\equiv\left|\vec{r}_1-\vec{r}_2\right|</math>, où <math>\vec{r}_1</math> et <math>\vec{r}_2</math> désignent les positions respectives du noyau et de l'électron par rapport à un repère d'origine ''O'' arbitraire<ref>En mécanique quantique les positions respectives du noyau et de l'électron sont données par les opérateurs dits « vectoriels » <math>\widehat{\vec{r}}_1</math> et <math>\widehat{\vec{r}}_2</math>, chacun d'eux correspondant en fait à l'ensemble de trois « composantes » que sont les opérateurs « positions » liés à chacune des coordonnées, <math>\widehat{\mathbf{\mathit{x}}}_1\text{ , }\widehat{\mathbf{\mathit{y}}}_2\text{ , etc.}</math>. Toutefois en représentation position, qui sera seule utilisée, ces opérateurs se ramènent à la multiplication de la fonction d'onde <math>\Psi(\vec{r},t)</math> par la coordonnée correspondante, il est alors possible d'utiliser les notations vectorielles ordinaires <math>\vec{r}_1</math> et <math>\vec{r}_2</math> pour désigner les positions respectives du noyau et de l'électron.</ref>.

Dans un tel modèle les effets relativistes de tout ordre, ainsi que l'influence du spin tant électronique que nucléaire, et les diverses corrections prenant en compte la singularité de la distribution de charge de l'électron comme du noyau, ne sont pas pris en compte. En toute rigueur, tenir compte des effets relativistes nécessiterait d'utiliser l'[[équation de Dirac]] pour l'électron dans un champ extérieur<ref>Cf. {{landau|tome 4}}, chapitre IV.</ref>. Toutefois, l'atome d’hydrogène étant un système où les effets relativistes demeurent faibles, tout comme ceux liés aux interactions entre les spins nucléaire et électronique, il est possible de considérer les différents effets, très faibles devant l'interaction électrostatique largement dominante, comme des termes de perturbations<ref>Ainsi, les termes dits de « structure fine » sont en valeur moyenne {{Unité|1000}} fois plus faibles que ceux du hamiltonien non perturbé, et ceux dits de « structure hyperfine » sont eux-mêmes près de {{Unité|2000}} fois plus faibles que ceux de structure fine, cf. {{Cohen}}, tome II, chapitre XII.</ref>. Ils sont néanmoins très importants pour expliquer les structures dites [[structure fine|fine]] et [[structure hyperfine|hyperfine]] des [[raie spectrale|raies spectrales]]. Négliger ces différents termes présente l'avantage de pouvoir considérer une situation dans laquelle l'équation de Schrödinger peut être résolue exactement, tout en reproduisant avec une bonne précision les différents résultats expérimentaux sur le spectre de l'atome d'hydrogène notamment.

À cette approximation, l'étude de l'atome d'hydrogène, qui conduit à la détermination de ses états électroniques, est celle du mouvement d'un électron dans un champ central (cf. l'article [[mouvement dans un champ central symétrique]]), et correspond en fait à un « problème à deux corps quantique ». Il est alors possible d'utiliser les résultats généraux de ce type de mouvement : dans le référentiel du centre de masse du système l'étude du mouvement se ramène à celui d'une ''particule fictive'' de masse <math>\mu=\frac{m_e m_p}{m_e + m_p} \approx m_e</math> dont la position est donnée par <math>\vec{r}=\vec{r}_1-\vec{r}_2</math>, se déplaçant avec une énergie potentielle électrique <math>V(r)=-\frac{e^2}{4\pi \varepsilon_0 r}</math>, le hamiltonien du système se réduisant à :
:<math>\widehat{\mathbf{\mathit{H}}}_0=\frac{\widehat{\mathbf{\mathit{p}}}^2}{2\mu}-\frac{e^2}{4\pi \varepsilon_0 r}</math><ref>L'indice « 0 » est là pour rappeler qu'il s'agit du hamiltonien non perturbé.</ref>{{,}}<ref>Bien entendu le mouvement libre du centre de masse est négligé.</ref>.

En pratique, comme ''m{{ind|p}} >> m{{ind|e}}'', le centre de masse est pratiquement confondu avec le centre du noyau, et il est fréquent d'effectuer l'approximation du noyau ''infiniment lourd'', pour laquelle <math>\mu \approx m_e</math> dans les formules : la « particule fictive » se confond alors avec l'[[électron]]. Dans la suite, aucune distinction ne sera faite entre l'électron et la particule fictive.

== Détermination des états propres électroniques ==
La fonction d'onde des états propres électroniques s'écrit sous la forme générale (en représentation position) : <math>\psi(\vec{r},\sigma)=\psi(\vec{r})\chi(\sigma)</math>, où <math>\psi(\vec{r})</math> est la fonction d'onde ''spatiale'' et <math>\chi(\sigma)</math> la fonction d'onde de ''[[spin]]'' (''s'' est la « variable de spin » indiquant la valeur de la projection du spin sur un axe donné). Ceci étant, le hamiltonien <math>\widehat{\mathbf{\mathit{H}}}_0</math> n'agit pas sur la variable de spin de l'électron, et donc les états électroniques sont uniquement déterminés par la donnée de <math>\psi(\vec{r})</math>, la prise en compte du spin ''s'' de l'électron résulte uniquement en une dégénérescence double de chaque état propre électronique puisque ''s=1/2'' et donc il existe deux valeurs possibles de la projection du spin de l'électron <math>m_s=\pm\tfrac{1}{2}</math><ref>Cette prise en compte du spin comme simple « cause de multiplicité » de l'état électronique est bien sûr directement liée au caractère non-relativiste de la théorie adoptée : le développement de l'équation de Dirac à l'ordre <math>1/c^2</math> ferait apparaître un terme d'interaction ''spin-orbite'' de la forme <math>\xi(r)\widehat{\mathbf{\mathit{l}}}\cdot\widehat{\mathbf{\mathit{s}}}</math>, responsable d'une levée partielle de la dégénérescence liée au spin, cf. {{Cohen}}, tome II, chapitre XII.</ref>.

=== Équation de Schrödinger – Séparation radiale-angulaire ===
Le hamiltonien étant indépendant du temps, l'étude quantique de l'atome d'hydrogène implique de résoudre l'équation de Schrödinger stationnaire suivante pour la fonction d'onde <math>\psi(\vec{r})</math>, en représentation position :
:<math>\widehat{\mathbf{\mathit{H}}}_0\psi=E\psi</math> soit <math>\left[\frac{\hbar^2\Delta}{2\mu}+\left(E+\frac{e^2}{4\pi\varepsilon_0 \; r}\right)\right]\psi(\vec{r})=0</math>.

Dans la théorie du [[mouvement dans un champ central symétrique]] il apparaît qu'en utilisant les coordonnées sphériques ''(r,θ,ϕ)'' comme l'implique la symétrie du hamiltonien, la fonction d'onde correspond aux états propres communs de <math>\widehat{\mathbf{\mathit{H}}}_0</math> et des opérateurs liés au [[moment cinétique orbital]] de l'électron <math>\widehat{\mathbf{\mathit{L}}}^2</math> et <math>\widehat{\mathbf{\mathit{L}}}_z</math>. Elle se décompose donc en un terme ''radial'' et un terme angulaire, ce dernier correspondant aux [[harmonique sphérique|harmoniques sphériques]] (''k'' est le « [[nombre quantique]] » associé à l'énergie, ''a priori'' réel et continu) :
:<math>\psi(\vec{r})=R_{k,\ell}(r)Y_{\ell,m}(\theta,\phi)=\frac{u_{k,\ell}(r)}{r}Y_{\ell,m}(\theta,\phi)</math>,
la fonction <math>u_{k,\ell}(r)</math> étant solution de l'équation de Schrödinger radiale obtenue par substitution de l'expression précédente de la fonction d'onde dans l'équation de Schrödinger :
:<math>\left[\frac{\hbar^2}{2\mu}\frac{d^2}{dr^2}+\left(E_{k,\ell}+\frac{e^2}{4\pi\varepsilon_0 \; r}-\frac{\hbar^2 \ell(\ell+1)}{2\mu r^2}\right)\right] u_{k,\ell}(r)=0</math><ref>La notation <math>E_{k,\ell}</math> indique qu'a priori l'énergie de l'état propre correspondant à la fonction d'onde dépend des nombres quantiques ''k'' et <math>\ell</math>.</ref>.

Du point de vue physique, les solutions acceptables de l'équation de Schrödinger précédente sont celles qui ont un comportement régulier sur le domaine <math>[0,+\infty[</math>, et normalisées selon la condition :
:<math>\int_{0}^{+\infty}\left|R_{k,\ell}(r)\right|^2 r^2 \mathrm dr= \int_{0}^{+\infty} \left|u_{k,\ell}(r)\right|^2 \mathrm dr = 1</math>.

Ceci montre que physiquement <math>u_{k,\ell}(r)</math> est la fonction d'onde correspondant à un problème unidimensionnel, celui d'une particule astreinte à évoluer dans le potentiel dit ''effectif'' <math>U_{\mathrm{eff}}(r)=-\frac{e^2}{4\pi\varepsilon_0 \; r}+\frac{\hbar^2 \ell(\ell+1)}{2\mu r^2}</math>, et infini pour <math>r<0</math>.

''Remarque :'' pour les états continus la condition de normalisation précédente devra être prise « au sens des distributions » : <math>\int_{0}^{+\infty} u_{k',\ell}^*(r)u_{k,\ell}(r) \mathrm dr = \delta(k'-k)</math>, ''k'' et ''k’'' étant réels.

=== Résolution de l'équation radiale ===
La résolution mathématique de l'équation radiale nécessite de suivre plusieurs étapes. En premier lieu, il convient de mettre l'équation radiale précédente sous une forme adimensionnelle, puis de rechercher des solutions à l'équation ainsi obtenue qui soient '''régulières''' à l'origine et tendent vers 0 à l'infini (sinon elle ne serait pas de [[carré sommable]]). Un changement de fonction est alors nécessaire pour achever la résolution.

Tous calculs faits, les résultats sont les suivants :
* les états liés ont des valeurs propres de l'énergie quantifiées et données par <math>E_n = -\frac{E_\mathrm{I}}{n^2}</math>,
:avec <math>E_\mathrm{I} = \frac{\mu e^4}{32\pi^2 \varepsilon_0^2 \hbar^2}</math>, et <math>n=1,2,\ldots,+\infty</math>, qui correspond à l'énergie d'ionisation de l'atome d'hydrogène et est identique à celle obtenue dans le [[modèle de Bohr]] ;
* les fonctions d'onde normalisées correspondantes ont pour expressions : <math>\psi_{n,\ell,m}(r,\theta,\phi)=\sqrt {{\left ( \frac{2}{n r_0} \right )}^3\frac{(n-\ell-1)!}{2n(n+\ell)!} }L_{n-\ell-1}^{2\ell+1}\left(\tfrac{2r}{nr_0}\right)\left(\dfrac{2r}{nr_0}\right)^{\ell}e^{-\tfrac{r}{nr_0}}Y_{\ell,m}(\theta,\phi)</math>,
:avec <math>r_0 \equiv \frac{4\pi\varepsilon_0\hbar^2}{\mu e^2}</math>, qui correspond dans le cas où l'on assimile [[masse réduite]] et masse de l'électron au [[rayon de Bohr]] de l'atome d’hydrogène, <math>L_{n-\ell-1}^{2\ell+1}</math> [[polynôme de Laguerre|polynômes généralisés de Laguerre]], et <math>n=1,2,...,+\infty \text{ , } \ell=0,1,...(n-1) \text{ , } m=-\ell,-\ell+1,...,\ell-1,\ell</math>.

{{démonstration|contenu=
==== Passage à une équation sans dimension ====
L'équation radiale se met aisément sous la forme :
:<math>\left[\frac{d^2}{dr^2}+\left(\frac{2\mu E}{\hbar^2}+\frac{2 \mu e^2}{4\pi\varepsilon_0\hbar^2}\frac{1}{r}-\frac{\ell(\ell+1)}{r^2}\right)\right]u_{k,\ell}(r) = 0</math>.

L'analyse dimensionnelle montre sans difficulté que la quantité <math>\frac{\mu e^2}{4\pi\varepsilon_0\hbar^2} \equiv \frac{1}{r_0}</math> a la dimension de l'inverse d'une longueur. Cela donne l'échelle « naturelle » du problème et il est possible de poser <math>r_0 \equiv \frac{4\pi\varepsilon_0\hbar^2}{\mu e^2}</math><ref>Dans le cas où l'on assimile masse réduite et masse de l'électron cette quantité n'est autre que le [[rayon de Bohr]] <math>a_0 = \frac{4\pi\varepsilon_0\hbar^2}{m_e e^2}\approx\text{53 pm}</math>, c'est-à-dire le rayon de la première orbite dans le [[modèle de Bohr|modèle semi-classique de Bohr]].</ref>.

Il est dès lors possible d'introduire naturellement une variable sans dimension <math>x \equiv \frac{r}{r_0}</math>, l'équation précédente peut alors être mise sous une forme adimensionnelle :
:<math>\frac{d^2}{dx^2}u_{k,\ell}+\left(\frac{2\mu r_0^2 E_{k,\ell}}{\hbar^2}+\frac{2}{x}-\frac{\ell(\ell+1)}{x^2}\right)u_{k,\ell}(x)=0</math>.
Sans perte de généralité, il est possible de poser :
:<math>\frac{2\mu r_0^2 E_{k,\ell}}{\hbar^2}=\pm\frac{1}{\beta^2}</math>,
<math>\beta=\beta(k,\ell)</math> étant réel, le signe (-) étant pour le cas où ''E<0''. Avec cette notation l'équation à résoudre est de la forme :
:<math>\frac{d^2}{dx^2}u_{k,\ell}+\left(\frac{2}{x}-\frac{\ell(\ell+1)}{x^2}\pm \frac{1}{\beta^2}\right)u_{k,\ell}(x)=0</math>.

==== Comportement asymptotique des solutions de l'équation radiale ====
Si <math>x\rightarrow + \infty</math> l'équation adimensionnelle précédente devient :
:<math>\frac{d^2}{dx^2}u_{k,\ell}\approx\mp \frac{u_{k,\ell}(x)}{\beta^2}</math>,
équation dont les solutions générales dépendent du signe de ''E'' :
* ''E>0'': dans ce cas <math>u_{k,\ell}(x) \sim A e^{-\tfrac{ix}{\beta}}+B e^{\tfrac{ix}{\beta}}</math> quand <math>x\rightarrow + \infty</math> (''A'' et ''B'' étant deux constantes d'intégration). Une telle solution correspond à une particule ''libre'', soit un état non-lié, continu: il s'agit là de l'atome ionisé.
* ''E<0'': en tenant compte de la nécessité d'avoir <math>u_{k,\ell}(x)</math> bornée le comportement asymptotique est de la forme <math>u_{k,\ell}(x) \sim Ae^{-\tfrac{x}{\beta}}</math>. Dans ce cas il est clair que la probabilité de présence de l'électron à grande distance du noyau tend vers 0, par suite cela correspond à un état ''lié'', ce qui était attendu physiquement.

Dans la suite il est uniquement utile de considérer le cas ''E<0'', et donc de rechercher les solutions de :
:<math>\frac{d^2}{dx^2}u_{k,\ell}+\left(\frac{2}{x}-\frac{\ell(\ell+1)}{x^2}- \frac{1}{\beta^2}\right)u_{k,\ell}(x)=0</math>,
ayant pour comportement asymptotique <math>u_{k,\ell}(x) \approx Ae^{-\tfrac{x}{\beta}}</math>.

==== Comportement à l'origine des solutions de l'équation radiale ====
De façon générale, pour un mouvement dans un champ central symétrique, les solutions régulières de l'équation radiale sont de la forme <math>Ar^{\ell+1}</math> au voisinage de l'origine (cf. [[Mouvement dans un champ central symétrique|article détaillé]]), par suite <math>u_{k,\ell}(x) \underset{r\rightarrow 0^+}{\sim} Ax^{\ell+1}</math>.

==== Détermination des états propres ====
Du fait des comportements asymptotique et au voisinage de l'origine des solutions physiquement acceptable de l'équation pour les fonctions d'onde radiales des états liés de l'atome d'hydrogène, il est utile d'introduire la fonction auxiliaire <math>\chi_{k,\ell}(y)</math> telle que :
:<math>u_{k,\ell}(y)=\chi_{k,\ell}(y) y^{\ell+1}e^{-\tfrac{y}{2}}</math>,
avec pour nouvelle variable <math>y=\tfrac{2x}{\beta}</math>.

Avec ce changement de fonction et de variable, l'équation radiale adimensionnelle précédente devient :
:<math>y\chi''_{k,\ell}+\left((2\ell+1)+1-y\right)\chi'_{k,\ell}+(\beta-\ell-1)\chi_{k,\ell}(y)=0</math>.

Cette équation différentielle a pour solution la [[fonction hypergéométrique confluente]] <math>_1F_1(-(\beta-\ell-1);(2\ell+1);y)</math>. Or cette fonction n'a de comportement régulier à l'infini que si <math>\beta-\ell-1</math> est un ''entier positif ou nul'', donc puisque <math>\ell</math> est un entier <math>\beta \equiv n =1,2,\ldots,+\infty</math>. La nécessité de régularité à l'infini des solutions physiquement acceptables de l'équation (17) implique que les états liés soient '''quantifiés''', les valeurs de <math>\ell</math> admissibles étant telles que :<math>\ell=0,\ldots,(n-1)</math>, (17bis).

Dans ce cas, la fonction <math>_1F_1(-(n-\ell-1);(2\ell+1);y)</math> se réduit alors (par troncation de la série hypergéométrique) aux [[polynôme de Laguerre|polynômes généralisés de Laguerre]] <math>L_{n-\ell-1}^{2\ell+1}(x)</math>.

L'état propre de fonction d'onde <math>\psi_{n,\ell,m}(\vec{r})</math> a pour énergie propre :
:<math>E_n = -\frac{E_\mathrm{I}}{n^2}\text{, avec }E_\mathrm{I} = \frac{\mu e^4}{32\pi^2 \varepsilon_0^2 \hbar^2}\text{, et }n=1,2,\ldots,+\infty</math>.

La quantité <math>E_\mathrm{I}</math> correspond à l'énergie d'ionisation de l'atome d'hydrogène et est identique à celle obtenue dans le [[modèle de Bohr]].

Finalement les fonctions d'onde radiales solutions de l'équation de Schrödinger radiale sont de la forme :
:<math>u_{n,\ell}=A \; L_{n-\ell-1}^{2\ell+1}\left(\tfrac{2r}{nr_0}\right)\left(\tfrac{2r}{nr_0}\right)^{\ell+1}e^{-\tfrac{r}{nr_0}}</math>,
''A'' étant une constante de normalisation.

Regroupant l'ensemble des résultats précédents les fonctions d'onde normalisées pour un état propre <math>|n,\ell,m\rangle</math> ont pour expression :
:<math>\psi_{n,\ell,m}(r,\theta,\phi)=\sqrt {{\left ( \frac{2}{n r_0} \right )}^3\frac{(n-\ell-1)!}{2n(n+\ell)!} }L_{n-\ell-1}^{2\ell+1}\left(\tfrac{2r}{nr_0}\right)\left(\dfrac{2r}{nr_0}\right)^{\ell}e^{-\tfrac{r}{nr_0}}Y_{\ell,m}(\theta,\phi)</math>,

avec <math>\begin{cases} n=1,2,...,+\infty \\ \ell=0,1,...(n-1) \\ m=-\ell,-\ell+1,...,\ell-1,\ell \end{cases}</math>.
}}

L'énergie de l'électron dans un état propre <math>\psi_{n,\ell,m}(\vec{r})</math> ne dépend donc ''pas'' du nombre quantique <math>\ell</math>. À un même état propre de nombre quantique {{mvar|n}} correspondent alors :
* {{mvar|n}} valeurs distinctes de <math>\ell</math>, de {{math|0}} à <math>n-1</math> ;
* pour chaque valeur de <math>\ell</math>, <math>(2\ell+1)</math> valeurs de {{mvar|m}} liées à l'existence de la dégénérescence dite essentielle<ref>La non-dégénérescence par rapport à {{mvar|m}} provient du fait de l'isotropie du hamiltonien, aucune direction de l'espace n'étant privilégiée, et est donc dite pour cette raison ''essentielle''. Cette dégénérescence peut être levée dans un champ extérieur, par exemple magnétique : on aura une séparation de la raie spectrale en plusieurs raies en présence du champ, c'est l'[[effet Zeeman]].</ref> des niveaux d'énergie par rapport à ce nombre quantique ;
* les deux valeurs possibles <math>m_s=\pm \tfrac{1}{2}</math> de la projection du spin de l'électron selon O''z''.
Finalement, chaque état de nombre quantique {{mvar|n}} est donc <math>2\sum_{k=0}^{n-1} (2k+1) = 2n^2</math> fois dégénéré.

La dégénérescence par rapport au nombre quantique <math>\ell</math> (dit aussi nombre quantique azimutal) est un cas de ''dégénérescence accidentelle'', lié en fait à l'existence d'une symétrie additionnelle pour le champ coulombien (cf. article [[vecteur de Runge-Lenz]]){{,}}<ref>La dégénérescence « accidentelle » est fondamentalement liée à la nature coulombienne du champ subi par l'électron. Son équivalent classique est l'existence de trajectoires fermées (ellipses, à la limite cercle) dans le cas d'un champ coulombien attractif (comme le champ gravitationnel : cf. [[problème à deux corps]]). Si le champ réel n'est pas rigoureusement coulombien, par exemple du fait de l'existence de perturbations, cette dégénérescence peut être levée, et on aura <math>E_{n,\ell}</math>.</ref>.

== Propriétés des états propres électroniques - orbitales atomiques ==
L'étude des propriétés des états propres électroniques, c'est-à-dire des fonctions d'onde radiale <math>R_{n,\ell}(r)</math> et angulaire <math>Y_{\ell,m}(\theta,\phi)</math> permet de déterminer la répartition spatiale des densités de probabilités de présence de l’électron. De façon générale, la probabilité de trouver l'électron dans l'[[Volume#Volume élémentaire|élément de volume]] <math>\mathrm dV=r^2\mathrm dr \mathrm d\Omega</math> (<math>\mathrm d\Omega</math> désignant l'[[angle solide]] élémentaire autour de la direction de <math>\vec{r}</math>) situé autour du point de coordonnées <math>(r,\theta,\phi)</math> est donnée par:
:<math>\mathrm dP=\left|\psi_{n,\ell,m}\left(\vec{r}\right)\right|^2\mathrm dV=r^2\left|R_{n,\ell}(r)\right|^2\mathrm dr \left|Y_{\ell,m}(\theta,\phi)\right|^2\mathrm d\Omega</math>,

par suite la densité de probabilité correspondante est le produit d'une densité ''radiale'' donnée par <math>r^2\left|R_{n,\ell}(r)\right|^2</math>, et d'une densité de probabilité ''angulaire'', soit <math>\left|Y_{\ell,m}(\theta,\phi)\right|^2</math>.

Par suite, il est possible de distinguer diverses « catégories » d'états d'énergie selon la valeur du nombre quantique (dit parfois ''secondaire'') <math>\ell</math>, qui détermine largement la répartition angulaire de la probabilité de présence de l'électron. À chacune des valeurs de <math>\ell</math> est associée une fonction d'onde de type particulier appelé [[orbitale atomique]]:
* si <math>\ell = 0</math> la densité de probabilité angulaire ne dépend pas de <math>(\theta,\phi)</math>, puisque <math>Y_{00}=1/\sqrt{4\pi}</math>. Dans ce cas, et à un facteur de normalisation près, la densité de probabilité de présence de l'électron est purement radiale et possède la symétrie ''sphérique''. Ces orbitales sont dites de type '''''s''''' (de l'anglais "sharp", terminologie historique provenant de la forme des raies spectrales);
* si <math>\ell = 1</math> il est possible de montrer que la densité de probabilité de présence possède une symétrie ''axiale''. Ces orbitales sont dites de type '''''p''''' (de l'anglais "principal");
* si <math>\ell = 2</math> il est possible de montrer que la densité de probabilité de présence possède des plans de symétrie. Ces orbitales sont dites de type '''''d''''' (de l'anglais "diffuse");
* si <math>\ell = 3</math> la densité de probabilité de présence possède des propriétés de symétrie plus complexe. Ces orbitales sont dites de type '''''f''''' (de l'anglais "fundamental").

Il convient dans un premier temps de s'intéresser au cas le plus simple, celui de l'état fondamental de l'atome d'hydrogène, qui correspond à <math>n=1</math> et donc à <math>\ell=m=0</math>.

=== État fondamental de l'atome d'hydrogène ===
En fait, dans le cas de l'atome d'hydrogène, on peut trouver la solution de l'état fondamental (c'est-à-dire de plus basse énergie) '''rigoureusement''', en s'aidant uniquement du [[Principe d'incertitude|principe d'incertitude d'Heisenberg]]. C'est une façon très élégante de procéder, sans beaucoup de mathématiques.

En effet, très vite, (en 1929), [[Werner Heisenberg]] fait comprendre un des points-clefs de la mécanique quantique : Les [[Grandeur physique|grandeurs physiques]] ne sont plus des fonctions <math>f \left (\vec r , \vec p \right )</math> de l'espace de la position et de la vitesse (appelé en mécanique classique hamiltonienne, l'espace des phases) : '''cet espace n'est pas pertinent en mécanique quantique'''. Les grandeurs physiques <math>f \left (\vec r , \vec p \right )</math> doivent être remplacées par des [[opérateur linéaire|opérateurs linéaires]] [[observable]]s sur un [[espace de Hilbert|espace vectoriel de Hilbert]] et les valeurs propres, réelles, de ces matrices seront les valeurs expérimentalement mesurées. Comme l'opérateur position <math>\hat x </math> et l'opérateur impulsion <math>\hat p_x</math> '''[[Loi commutative|ne commutent pas]]''', il en résulte le théorème d'[[Principe d'incertitude|inégalité de Heisenberg]]) :

<center>
{| border=0
|-----
|
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-----
| <math>[\hat p,\hat x] = i\hbar \Rightarrow \Delta p \cdot \Delta x \ge {\hbar \over 2}</math>.
|}
|
| |
|}</center>


Alors, dans le cas d'égalité stricte — on dit que l'inégalité a été saturée à sa limite — la [[saturation des inégalités d'Heisenberg]] donne un moyen '''rigoureux''' de calculer la fonction d'onde,<math>\Psi_{1s}(x,y,z)</math>, de l'état fondamental de l'atome d'hydrogène.

Ce problème de valeur propre et de vecteur propre est donc résolu, dans l'article [[Saturation des inégalités d'Heisenberg]], pour la plus basse énergie (cf. [[Modèle de Bohr|atome de Bohr]]) ; et cela donne :

<center>
{| border=0
|-----
|
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-----
| <math>E_1= -{me^4 \over 2\hbar^2}</math>
|}
|
| |
|}</center>

<center>
{| border=0
|-----
|
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-----
| <math>\psi_{1s}(r) = N \cdot e^{-r/a_0}</math>
|}
|
| |
|}</center>

''N'' étant la constante, réelle, dite de normalisation de la probabilité.

=== Vérification ===
On va se contenter, ici, de vérifier que ceci est vrai en insérant directement cette solution dans l'équation de Schrödinger.

Premièrement, dans cette équation, la variable temps se sépare immédiatement :
: <math>\Psi(x,y,z,t) = \psi(x,y,z) e^{-i{Et \over \hbar}}</math>
dans ce cas dit stationnaire, cela amène à trouver les valeurs propres de l'opérateur linéaire H dans l'espace ''L''{{exp|2}} des fonctions des trois variables ƒ(''x'', ''y'', ''z'') à valeur complexe, de carré sommable :
: <math>\mathbf{H}\psi(x,y,z) = - { \hbar^2 \over 2m } \nabla^2 \psi(x, y, z) - { e^2 \over r} \psi(x,y,z) = E \psi(x,y,z) </math>.

Or, dans ce cas, cette fonction uniquement de ''r'' a pour Laplacien, la valeur usuelle <math> \Delta f(r) = f'' + {2 \over r} f'</math>.

De plus, on se sert évidemment des [[système d'unités atomiques|unités atomiques]], qui a été introduit à cet usage. Cela revient à faire <math>[ \hbar = m = e^2 = 1]</math> dans les calculs ; Landau (p142) appelle ce [[système d'unités]] coulombiennes :
: ƒ" = ƒ, ƒ' = &minus;ƒ ;
donc il s'agit de vérifier si :
: &minus;1/2 &centerdot; (ƒ + 2/''r'' &centerdot; (&minus;ƒ)) + 1/''r'' &centerdot; ƒ = &minus;1/2 &centerdot;ƒ
qui est vrai.

=== Densité de probabilité de présence ===
On en déduit aussitôt la probabilité d''p'' de trouver l'électron à une distance du noyau comprise entre ''r'' et ''r'' + d''r'' : elle est donnée par d''p'' = ''P''(''r'')• d''r'' :
: <math>P(r) = 4 \left(\frac{1}{a_0}\right)^3 r^2 e^\frac{-2r}{a_0}</math>.

Sur le graphique de la densité de probabilité, la distance au noyau est donnée en multiple du premier rayon de Bohr, on voit immédiatement que la probabilité est maximale au premier rayon de Bohr :

<!-- Le graphique va bientôt être rectifié. - (en particulier, il faut enlever ce malencontreux 4 qui perturbe la lecture, merci d'avance). (Il est mieux représenté dans [[atome]])-->

=== Orbitale 1s ===
Cette solution s'appelle en chimie l'orbitale 1s.

On pourra vérifier le [[théorème du viriel]] :
: moyenne de 1/''r '' = <1/''r''> = 1/''a{{ind|0}}''

et le théorème d'Ehrenfest :
: moyenne de 1/''r{{exp|2}}'' = <1/''r{{exp|2}}''> = 2/''a{{exp|2}}''

La moyenne de ''r'' n'est pas ''a'', mais (3/2)·''a'' ; [''de manière générale, l'inverse de la moyenne n'est pas la moyenne de l'inverse'' ].

Et la moyenne de ''r{{exp|2}}'' vaut : 3''a{{exp|2}}'',
donc la variance de ''r'' (<''r''{{exp|2}}>-<''r''>{{exp|2}} )vaut (3-9/4)·''a{{exp|2}}'' = 0,75·''a{{exp|2}}'' ; soit un écart-type de 0,866·''a'', ce qui est très grand.

L'électron est dit délocalisé dans l'espace. L'espace occupé reste malgré tout de volume fini : au-delà d'une distance de 3''a'', la probabilité de détecter l'électron devient très faible. Cet espace occupé décrit l'orbitale de l'électron : typiquement en [[chimie quantique]], on convient formellement de tracer la méridienne de la surface qui englobe en gros 98 % de chance d'y trouver l'électron :
: ici ''r ''= 3/2 + 1,732 ~ 3,2·''a''. Ceci est très conventionnel.

=== Note : espace des impulsions ===
L'opérateur impulsion a évidemment une moyenne nulle (symétrie sphérique), mais l'opérateur ''P{{exp|2}}'' vaut 2''m''·''E{{ind|c}}'', dont la valeur moyenne est par le théorème du viriel
: <''P{{exp|2}}''> = -2''m''·''E{{ind|c}}'', soit en unités atomiques +2·1/2 = 1.
Donc la variance de ''P'' vaut <math> <\psi|P^2|\psi> =1 \cdot {me^2 \over \hbar}</math>.

On retrouve bien (heureusement!) ce dont on était parti dans l'article ''[[Saturation des inégalités d'Heisenberg]]''.

Mais, on peut aller un peu plus loin [ ''ne jamais perdre de vue que l'espace des impulsions joue un rôle équivalent à celui des positions, bien qu'il soit moins étudié en chimie''] :

Remarque : représentation dite des impulsions
: La transformée de Fourier de <math>\psi(\vec r)</math> est <math>\phi(\vec p)</math>, avec la même règle de Born bien sûr : <math>|\phi(\vec p)|^2</math> donne la densité de probabilité dans l'espace des impulsions. Le calcul de la transformée de Fourier de exp-r donne 1/(1+''p''²)² et donc on peut calculer de même la distribution des impulsions et retrouver la variance de ''p'', et la valeur moyenne de l'[[énergie cinétique]] : il est très important de faire ces calculs pour bien comprendre que l'électron, bien que dans un état stationnaire, ne cesse de « bouger » : il est aussi délocalisé en impulsion.
: En fait, ce n'est pas une particule, ce n'est pas une onde, c'est une entité nouvelle, la « particlonde », qui ne satisfait plus les équations de la mécanique classique (dans la version dite d'Hamilton), ni l'équation des ondes de l'optique, ni l'équation de la diffusion, mais cette drôle d'équation, l'équation de Schrödinger qui ne se laisse appréhender '''concrètement''' que dans la vision lagrangienne de [[Paul Dirac|Dirac]] et [[Richard Feynman|Feynman]] (on parle d'[[intégrale de chemin]]s (l'équation de Schrödinger est ainsi parfois appelée équation de cheminement)). Si bien que l'énergie cinétique n'est pas du tout négligeable, puisque égale à 50 % de l'énergie potentielle en module.

Remarque : [[Hagen Kleinert|Kleinert]], élève de Feynman, a réussi à donner l'interprétation du « cheminement » dans le cas de l'atome d'hydrogène, ce qui reste une prouesse. En ce sens, pour les chimistes, le seul vrai progrès notable depuis Hartree-Fock et Clementi a été (pour l'atome à N électrons), la notion de densité fonctionnelle de Kohn (Nobel de chimie 1998).

=== Conclusion ===
Il faut garder en mémoire toujours ces deux aspects, '''le couple''' [<math>\Psi(r), \Phi(p)</math>], pour bien comprendre l'aspect non statique, mais stationnaire de cette délocalisation de l'électron.

Beaucoup de livres proposent comme règle empirique : si l'électron est localisé dans une région de l'ordre de ''r'' = ''a'', lui donner une énergie cinétique de l'ordre de <math>{\hbar^2 \over 2ma^2}</math>. Dans le cas présent, cela, donne une énergie totale <math>E ={\hbar^2 \over 2ma^2} - {e^2 \over a}</math> dont le minimum est bien <math>- {e^2 \over 2a}</math>, où ''a'' est le rayon de Bohr :
: <math>a = {\hbar^2 \over me^2}</math>.

C'est une façon simple et élégante d'introduire les OdG ([[ordre de grandeur]]) de l'atome, souvent reprise dans les bons ouvrages.

== Orbitales ==
=== Résonance acoustique ===
Nous allons obtenir la structure de l'atome d'hydrogène en dénombrant les modes de vibration de la sphère. Le mode de vibration fondamental, dont la fréquence est la plus élevée, est toujours sphérique. Lorsque la vitesse des ondes est constante, il correspond à une longueur d'onde égale à deux fois le diamètre, comme pour une corde vibrante où le fondamental a une longueur d'onde double de la longueur de la corde lorsqu'elle est encastrée ou libre à chacune de ses extrémités.

Dans un cube, les résonances apparaissent lorsque les demi-longueurs d'onde sont des fractions entières du côté du cube avec un ventre au centre du cube. Le premier harmonique du cube présente un nœud au milieu, comme pour la corde vibrante. Son [[nombre quantique principal]] est n=2. Il y a trois façons de le placer, selon les trois directions de l'espace, ce qui donne trois modes de vibration de même énergie, autrement dit dégénérés. Il peut aussi y avoir un nœud au centre du cube. Il y a donc 4 possibilités.
=== Résonance de l'atome ===
Il en est de même pour la sphère qui aura un nœud au centre soit un nombre quantique secondaire ''l''=1. Comme pour le cube, il y a trois façons de le placer, selon les trois directions de l'espace, ce qui donne trois modes de vibration de même énergie, autrement dit dégénérés. Il peut aussi y avoir un nœud de symétrie sphérique, en son centre. Il y a donc 4 possibilités comme pour le cube.

Lorsqu'on augmente la fréquence de vibration, c'est-à-dire l'énergie de vibration, le nombre de nœuds augmente d'une unité à chaque fois. On obtient ainsi les harmoniques successifs qui ne sont généralement pas des harmoniques au sens musical du terme. Les harmoniques, au sens mathématique, du tambour, par exemple, ne sont pas des harmoniques au sens musical car ils ne sont pas des multiples entiers du fondamental. Il en est de même dans l'atome.

Dans l'atome d'hydrogène où la vitesse des ondes de de Broglie est fonction du potentiel électrostatique du noyau, le mode fondamental correspond au nombre quantique principal n=1. La théorie de Schrödinger fait apparaître deux nombres quantiques supplémentaires, le nombre quantique secondaire ''l'' et le nombre quantique magnétique m, nuls pour le mode fondamental n=1.

On n'utilise pas les coordonnées cartésiennes pour la sphère, mais la colatitude <math>\theta</math> et la longitude <math>\varphi</math>. On peut se passer de la variable r grâce à Bohr car la théorie de Schrödinger prévoit les mêmes niveaux d'énergie. L'axe principal est vertical, celui pour lequel <math>\theta=0</math>. Pour le premier harmonique, n = 2, on a trois orientations possibles pour le nœud, un nœud selon l'équateur, les deux autres des méridiens. On pourrait prendre des méridiens perpendiculaires, mais il revient au même de n'en prendre qu'un seul et de le faire tourner dans un sens ou dans l'autre, ce qui correspond au nombre quantique magnétique ''m'' = ± 1. On fait donc varier ''m'' entre - 1 et + 1.

En résumé, le nombre quantique principal ''n'' donne le nombre de nœuds. Le nombre quantique secondaire ''l'' < ''n'' donne le nombre de configurations possibles pour les nœuds et ''m'' les numérote de - (''l'' - 1) à ''l'' - 1.

Les orbitales sont d'abord représentées de façon simplifiée avec les nœuds des harmoniques sphériques sans les ventres (les « larmes »). La représentation est polaire, comme la Terre, en coordonnées sphériques.

On donne ensuite une représentation sous forme de « larmes » puis un tableau récapitulatif.

Il y a deux fois plus d'états quantiques que de modes de vibration en vertu du principe d'exclusion de Pauli.

=== Couche K ===
==== 1 orbitale sphérique (1s) ====
C'est l'état fondamental, de symétrie sphérique 1s, un nœud de vibration sphérique, qu'on peut placer soit sur la périphérie, à l'infini, soit sur le noyau. Les nombres quantiques correspondants sont :

:[[Image:Circle - black simple.svg|50px]] <math>n=1, l=0, m=0</math>

Il n'y a qu'un mode de vibration car les valeurs ''m'' = + 0 = ''m'' = - 0. En vertu du principe d'exclusion de Pauli, la couche K n'a qu'une seule orbitale et ne peut contenir que deux électrons au maximum. Avec un électron, on a l'hydrogène. Avec deux, on a l'hélium.

=== Couche L ===
==== 1 orbitale sphérique (2s) ====
Elle comprend une orbitale sphérique 2s, soit un nœud de vibration sphérique et deux états quantiques, donc deux éléments (Li et Be) :
:[[Image:L2s.svg|50px]] <math> n = 2 ; l = 0 ; m = 0</math>

==== 3 orbitales (2p) ====
Une orbitale de symétrie de révolution et deux autres avec un méridien, de symétrie axiale. Le méridien pouvant tourner dans un sens ou dans l'autre, il y a deux valeurs du nombre quantique magnétique ''m'' :

:[[Image:L2p0.svg|50px]] <math>n = 2 ; l = 1 ; m = 0</math>
:[[Image:L2p1.svg|50px]] <math>n = 2 ; l = 1 ; m = \pm 1</math>
En additionnant les orbitales des couches K et L, on a 5 orbitales soit, en vertu du principe de Pauli, 10 électrons et un numéro atomique N = 10 correspondant au néon. Cela permet non seulement de comprendre l'atome d'hydrogène mais aussi de construire la [[tableau périodique des éléments|table de Mendeleiev]].

=== Couche M ===
==== 1 orbitale (3s) sphérique ====
:[[Image:M3s.svg|50px]] <math>n = 3 ; l = 0 ; m = 0</math>

==== 3 orbitales (3p) ====
Comme pour la couche L sauf qu'il y a un nœud sphérique en plus.
:[[Image:M3p0.svg|50px]] <math>n = 3 ; l = 1 ; m = 0</math>
:[[Image:M3p1.svg|50px]] <math>n = 3 ; l = 1 ; m = \pm 1</math>

==== 5 orbitales (3d) ====
Les orbitales m ? 0 sont doubles.
:[[Image:M3d0.svg|50px]] <math>n = 3 ; l = 2 ; m = 0 </math> ( <math>d_{z^2}</math> symétrie de révolution)
:[[Image:M3d1.svg|50px]] <math>n = 3 ; l = 2 ; m = \pm 1</math> (trèfle à quatre lobes <math>d_{xy}</math> et <math>d_{x^2-y^2}</math>)
:[[Image:M3d2.svg|50px]] <math>n = 3 ; l = 2 ; m = \pm 2</math> (trèfle à quatre lobes <math>d_{xz}</math> et <math>d_{yz}</math>)

On peut aussi représenter les orbitales d de la couche M, en larmes d'eau :
:[[Image:Orbitales d.jpg]]

Pour plus de détails voir [[harmonique sphérique|harmoniques sphériques]]. On remarquera l'anneau sur l'orbitale <math>d_{z^2}</math> de révolution qu'on retrouve sur les atomes hydrogénoïdes ou les [[atome de Rydberg|atomes de Rydberg]].

=== Récapitulatif des couches K, L, M ===
:[[Image:HAtomOrbitals.png]]

Cette figure résume les modes de vibration qu'on rencontre dans les trois couches K, L, M. Chaque couche reprend la couche précédente avec un nœud de plus.

Les orbitales de gauche, sphériques, sont simples. Les orbitales p sont triples avec un nœud plan et les orbitales d quintuples avec deux nœuds plans.

Chaque couche contient les couches inférieures, par exemple en dessous de la sous-couche 3p, on a les sous-couches 1s, 2s, 2p et 3s. D'après le principe d'exclusion de Pauli, le nombre maximal d'électrons dans une sous-couche doit être pair. Considérons le cas du sodium Na. Toutes les couches seront remplies jusqu'à la sous-couche 3s avec un électron célibataire. Sa structure électronique s'écrit <math>1s^22s^22p^63s^1</math> ou <math>[Ne]s^1</math>, ce qui fait un total de Z = 11 électrons, numéro atomique du sodium Na qui suit immédiatement Ne. On marie l'électron célibataire du sodium pour obtenir le magnésium. On remplit la couche 3p avec six électrons pour obtenir l'argon Ar. Tous les gaz rares ont une sous-couche externe de type p sauf l'hélium avec une sous-couche s.

Au-delà de l'argon, il y a une anomalie due à la répulsion électrostatique entre les électrons : le mode 4s se met à la place de 3d. Le modèle de l'atome d'hydrogène ne s'applique donc que jusqu'à l'argon. Au-delà, l'ordre des sous-couches diffère de celui prévu par le modèle de la cavité sphérique comme de celui de Schrödinger. Pour prévoir l'ordre réel on a des règles empiriques (Aufbau, [[Friedrich Hund|Hund]], [[Règle de Klechkowski|Klechkowski]] ou Madelung). Connaissant ces anomalies, on peut construire la [[Tableau périodique des éléments|Table périodique des éléments]].

== Notes et références ==
{{Références}}

== Voir aussi ==
* [[Orbitale de l'atome d'hydrogène]]

=== Articles connexes ===
* [[Atome]]
* [[Modèle de Bohr|Atome de Bohr]]
* [[Saturation des inégalités d'Heisenberg]]
* [[Orbitale atomique]]
* [[Orbitale de l'atome d'hydrogène]]
* [[Harmonique sphérique|Harmoniques sphériques]]
* [[Mouvement dans un champ central symétrique]]
* [[vecteur de Runge-Lenz]]
* [[Atome à N électrons]]
* [[Tableau périodique des éléments|Classification périodique]]
* [[Dihydrogène]]
* [[Raie à 21 centimètres]]

=== Bibliographie ===
* {{Basdevant}}
* {{Cohen}}
* {{Messiah}}
* {{Landau|tome 3}}
* S.G. Karshenboim ''et al.'' (éditeurs) ; ''The hydrogen atom - Precision physics of simple atomic systems'', Lecture Notes in Physics 570, Springer-Verlag (2001). Recueil d'articles de revue sur l'état de l'art en [[spectroscopie atomique]], mesures de fréquences, et mesures de constantes fondamentales. Niveau troisième cycle universitaire.
* Victor Guillemin & Shlomo Sternberg ; ''Variations on a Theme by Kepler'', Providence RI, American Mathematical Society (1990), ASIN 0821810421. Très beau livre.
* Bruno Cordani ; ''The Kepler Problem - Group Theoretical Aspects, Regularization and Quantization, with Application to the Study of Perturbations'', Progress in Mathematical Physics '''29''', Birkhäuser (2003), {{ISBN|3-7643-6902-7}}.
* Stephanie F. Singer, ''Linear Symmetry and Predictions in Hydrogen Atom'', Undergraduate Text in Mathematics, Springer-Verlag (2005), {{ISBN|0-387-24637-1}}. De niveau plus simple, mais très moderne et beaucoup de références web.
* E. Hansch, A. Schawlaw & G. Series ; ''Le spectre de l'hydrogène atomique'', Pour La Science 19 ({{date-|mai 1979}}) 46.

{{portail|physique|chimie}}

[[Catégorie:Modèle quantique]]
[[Catégorie:Atome]]
[[Catégorie:Physique de l'hydrogène]]

[[pl:Wodór atomowy]]