Introduction to Mathematical Physics/N body problem in quantum mechanics/Atoms

One nucleus, one electronEdit


This case corresponds to the study of hydrogen atom.\index{atom} It is a particular case of particle in a central potential problem, so that we apply methods presented at section ---secpotcent--- to treat this problem. Potential is here:



It can be shown that eigenvalues of hamiltonian H with central potential depend in general on two quantum numbers k and l, but that for particular potential given by equation eqpotcenhy, eigenvalues depend only on sum n=k+l.

Rotation invarianceEdit


We treat in this section the particle in a central potential problem ([#References|references]). The spectral problem to be solved is given by the following equation:


Laplacian operator can be expressed as a function of L^2 operator.

Theorem: Laplacian operator \Delta can be written as:

\Delta=-\frac{1}{r^2}L^2+\frac{\partial^2}{\partial r^2}+\frac{2}{r}\frac{\partial}{\partial r}

Proof: Here, tensorial notations are used (Einstein convention). By definition:



L_iL_i&=&\epsilon_{ijk}\epsilon_{ilm}x_{j}p_k x_l p_m\\
&=&(\delta_{jl}\delta_{km}-\delta_{jm}\delta{kl})x_{j}p_kx_l p_m\\

The writing order of the operators is very important because operator do not commute. They obey following commutation relations:

[x_i,p_j]=i \hbar \delta_{ij}



From equation eqdefmomP, we have:

p_k=-i\hbar \frac{\partial}{\partial x_k}




&=&i\hbar x_kp_k+p_kx_jx_kp_j\\
&=&i\hbar x_kp_k+p_kx_kx_jp_j

Introducing operator:

\tilde D=x_k\frac{\partial}{\partial x_k}

we get the relation:

L^2=-x^2\Delta+\tilde D^2+\tilde D

Using spherical coordinates, we get:

\tilde D=r\frac{\partial}{\partial r}


\tilde D^2=(r\frac{\partial}{\partial r})(r\frac{\partial}{\partial
r})=r^2\frac{\partial^2}{\partial r^2}+\frac{\partial}{\partial r}

So, equation eql2pri becomes:

r^2}+\frac{2}{r}\frac{\partial}{\partial r}

Let us use the problem's symmetries:


  • L_z commutes with operators acting on r
  • L_z commutes with L^2 operator L_z commutes with H
  • L^2 commutes with H

we look for a function \phi that diagonalizes simultaneously H,L^2,L_z that is such that:


Spherical harmonics Y^m_l(\theta,\phi) can be introduced now:


Spherical harmonics Y^m_l(\theta,\phi) are eigenfunctions common to operators L^2 and L_z. It can be shown that:


Looking for a solution \phi(r) that can\footnote{Group theory argument should be used to prove that solution actually are of this form.} be written (variable separation):


problem becomes one dimensional:


r})+\frac{l(l+1)}{2\mu r^2}\hbar^2+V(r)]R_{l}(r)=E_{kl}R_{l}(r)

where R(r) is indexed by l only. Using the following change of variable: R_{l}(r)=\frac{1}{r}u_{l}(r), one gets the following spectral equation:



V_e(r)=\frac{l(l+1)}{2\mu r^2}\hbar^2+V(r)

The problem is then reduced to the study of the movement of a particle in an effective potential V_e(r). To go forward in the solving of this problem, the expression of potential V(r) is needed. Particular case of hydrogen introduced at section sechydrog corresponds to a potential V(r) proportional to 1/r and leads to an accidental degeneracy.

One nucleus, N electronsEdit

This case corresponds to the study of atoms different from hydrogenoids atoms. The Hamiltonian describing the problem is:

\frac{1}{4\pi\epsilon_0}\frac{Ze^2}{r_i}+ \sum_{j\mathrel{>}

where T_2 represents a spin-orbit interaction term that will be treated later. Here are some possible approximations:

N independent electronsEdit

This approximation consists in considering each electron as moving in a mean central potential and in neglecting spin--orbit interaction. It is a ``mean field approximation. The electrostatic interaction term


is modelized by the sum \sum W(r_i), where W(r_i) is the mean potential acting on particle i. The hamiltonian can thus be written:


where h_i=-\frac{\hbar^2}{2m}\Delta_i+W(r_i).


More precisely, h_i is the linear operator acting in the tensorial product space \otimes_{i=1}^N E_i and defined by its action on function that are tensorial products:

[1_1\otimes\dots\otimes 1_{i-1}\otimes h_i \otimes 1_{i+1}\dots\otimes

It is then sufficient to solve the spectral problem in a space E_i for operator h_i. Physical kets are then constructed by anti symmetrisation (see example exmppauli of chapter chapmq) in order to satisfy Pauli principle.\index{Pauli} The problem is a central potential problem (see section

secpotcent). However, potential W(r_i) is not like 1/r as in the

hydrogen atom case and thus the accidental degeneracy is not observed here. The energy depends on two quantum numbers l (relative to kinetic moment) and n (rising from the radial equation eqaonedimrr). Eigenstates in this approximation are called electronic configurations.


For the helium atom, the fundamental level corresponds to an electronic configuration noted 1s^2. A physical ket is obtained by anti symetrisation of vector:


Spectral termsEdit

Let us write exact hamiltonian H as:


where T_1 represents a correction to H_0 due to the interactions between electrons. Solving of spectral problem associated to H_1=H_0+T_1 using perturbative method is now presented.


It is here assumed that T_2<<T_1. This assumption is called L--S coupling approximation.

To diagonalize T_1 in the space spanned by the eigenvectors of H_0, it is worth to consider problem's symmetries in order to simplify the spectral problem. It can be shown that operators L^2, L_z, S^2 and S_z form a complete set of observables that commute.


Consider again the helium atom ([ph:mecaq:Cohen73]). From the symmetries of the problem, the basis chosen is:


where L is the quantum number associated to the total kinetic moment\index{kinetic moment}:

L\in\{l_1+l_2,l_1+l_2-1, \dots,|l_1-l_2|\}

and S is the quantum number associated to total spin of the system\index{spin}:


Moreover, one has:




Table Tab. tabpauli represents in each box the value of m_Lm_S for all possible values of m_L and m_S. \begin{table}[hbt]




For an atom with two electrons, states such that L+S is odd are excluded.


We will proof this result using symmetries. We have:

|1:n,l;2:n,l';L,M_L\mathrel{>}  \\
& &=\sum_m \sum_{m'}  \mathrel{<} l,l',m,m'|L,M_L\mathrel{>}

Coefficients \mathrel{<} l,l',m,m'|L,M_L\mathrel{>} are called Glebsh-Gordan\index{Glesh-Gordan coefficients} coefficients. If l=l', it can be shown (see ([ph:mecaq:Cohen73]) that:

\mathrel{<} l,l,m,m'|L,M_L\mathrel{>} =(-1)^L \mathrel{<} l,l,m',m|L,M_L\mathrel{>}.

Action of P_{21} on |1:n,l;2:n,l';L,M_L\mathrel{>} can thus be written:

P_{21}|1:n,l;2:n,l';L,M_L\mathrel{>} =(-1)^L|1:n,l;2:n,l';L,M_L\mathrel{>}

Physical ket obtained is:

& &= \left\{
0&\mbox{ if  }L+S\mbox{ is odd }\\
|1:n,l;2:n,l';L,M_L\mathrel{>} \otimes |S,M_S\mathrel{>} &\mbox{ if  }L+S\mbox{ is even }

Fine structure levelsEdit

Finally spectral problem associated to


can be solved considering T_2 as a perturbation of H_1=H_0+T_1. It can be shown ([ph:atomi:Cagnac71]) that operator T_2 can be written T_2=\xi(r_i)\vec l_i\vec s_i. It can also be shown that operator \vec J=\vec L+\vec S commutes with T_2. Operator </math>T_2</math> will have thus to be diagonilized using eigenvectors |J,m_J> common to operators J_z and J^2. each state is labelled by:


where L,S,J are azimuthal quantum numbers associated with operators \vec L,\vec S,\vec J.

Last modified on 30 May 2009, at 17:40