Pour interpréter correctement les spectres d'émission des atomes, plusieurs physiciens ont été amenés à proposer une théorie nouvelle reposant sur une notion de fonction d'onde. Nous allons voir avec Maple quelques-unes des propriétés de cette théorie qui est à la base de la Chimie moderne.

Atome d'hydrogène

En mécanique classique, l'équation du mouvement d'un électron soumis au champ électrique d'un proton supposé immobile s'écrit
E = p^2/(2*m)-e^2/(4*pi*epsilon[0]*r) où p=mv est la quantité de mouvement
Cette énergie E est une constante du mouvement. Si E est négative, les trajectoires décrites par l'électron sont des ellipses.
En mécanique quantique, un électron soumis au champ électrique du proton est décrit par l'équation de Schrödinger suivante

E*Psi(r) = -h^2*Delta*Psi(r)/(2*m*(2*pi)^2)-e^2*Psi... Psi(r) est appelée la fonction d'onde. abs(Psi(r))^2 est la probabilité de présence de l'électron au point r . La normalisation de la probabilité impose que int(int(int(abs(Psi(r))^2,r)*sin(theta),theta),phi)... On recherche une solution en coordonnées sphériques de la forme
Psi(r) = u[n](r)*Y[lm](theta,phi)/r Y[lm](theta,phi) sont les harmoniques sphériques, l un entier positif ou nul, et m un entier compris entre -l et +l.
On pose :

E[i] = -m*e^4/(2*(2*epsilon[0]*h)^2) et a[0] = h^2*epsilon[0]/(pi*m*e^2) . En effectuant le changement de variables suivant E1=E/Ei et r=r/a0, on obtient une équation pour la fonction d'onde
diff(u[lm](rho),`$`(rho,2))+(-l*(l+1)/(rho^2)+2/rho...


On s'intéresse aux solutions pour l=0 (orbitales s). Vérifier que si l'on prend E1=1/n^2 où n est un nombre entier strictement positif, on peut obtenir une solution explicite. Déterminer pour n=1 et n=2, la constante _C1 pour que la probabilité soit normée à l'unité.
Tracer les fonctions pour les valeurs de n allant de 1 à 5 et pour les valeurs de l permises correspondantes.

> eq:=diff(u(r),r$2)+(-l*(l+1)/r^2+2/r-1/n^2)*u(r);

eq := diff(u(r),`$`(r,2))+(-l*(l+1)/(r^2)+2/r-1/(n^...

> res1:=subs(_C2=0,dsolve({subs(n=1,l=0,eq)},u(r)));

res1 := u(r) = _C1*exp(-r)*r

> res11:=solve(int(rhs(res1)^2,r=0..infinity)=1);

res11 := 2, -2

> res111:=subs(_C1=res11[1],rhs(res1));

res111 := 2*exp(-r)*r

> res2:=subs(_C2=0,dsolve({subs(n=2,l=0,eq)},u(r)));

res2 := u(r) = _C1*exp(-1/2*r)*r*(-2+r)

> res21:=solve(int(rhs(res2)^2,r=0..infinity)=1);

res21 := 1/4*sqrt(2), -1/4*sqrt(2)

> res211:=subs(_C1=res21[2],rhs(res2));

res211 := -1/4*sqrt(2)*exp(-1/2*r)*r*(-2+r)

> res3:=subs(_C2=0,dsolve({subs(n=3,l=0,eq)},u(r)));

res3 := u(r) = _C1*exp(-1/3*r)*r*(27-18*r+2*r^2)

> res31:=solve(int(rhs(res3)^2,r=0..infinity)=1);

res31 := 2/243*sqrt(3), -2/243*sqrt(3)

> res311:=subs(_C1=res31[1],rhs(res3));

res311 := 2/243*sqrt(3)*exp(-1/3*r)*r*(27-18*r+2*r^...

> res4:=subs(_C2=0,dsolve({subs(n=4,l=0,eq)},u(r)));

res4 := u(r) = _C1*exp(-1/4*r)*r*(-192+144*r-24*r^2...

> res41:=solve(int(rhs(res4)^2,r=0..infinity)=1);

res41 := 1/768, -1/768

> res411:=subs(_C1=res41[1],rhs(res4));

res411 := 1/768*exp(-1/4*r)*r*(-192+144*r-24*r^2+r^...

> res5:=subs(_C2=0,dsolve({subs(n=5,l=0,eq)},u(r)));

res5 := u(r) = _C1*exp(-1/5*r)*r*(9375-7500*r+1500*...

> res51:=solve(int(rhs(res5)^2,r=0..infinity)=1);

res51 := 2/234375*sqrt(5), -2/234375*sqrt(5)

> res511:=subs(_C1=res51[1],rhs(res5));

res511 := 2/234375*sqrt(5)*exp(-1/5*r)*r*(9375-7500...

> plot({seq((res||i||11)^2,i=1..5)},r=0..50);