The Hydrogen Atom Analyzed Using Commutator Algebra I

©Fernando Caracena 2013

On a previous blog on the one-dimensional harmonic oscillator, we showed that the algebra of non-commuting operators in quantum mechanics could be exploited systematically to solve for the entire spectrum of state vectors and energy eigenvalues. Algebra sufficed for most of the calculation. In the following series we will develop the techniques of Lie algebra using non-commuting operators to solve for the state vectors and corresponding eigenvalues for the electrons in a hydrogen atom. In doing so, we will cover some techniques for systematically solving problems that have several, simultaneous eigenvalues. In three dimensional problems such a that of atomic structure, angular momentum plays a role in organizing the motion of electrons. We will find that angular momentum is the source of two additional eigenvalues in addition to that of energy.

The Problem of the Hydrogen Atom

The electron moving in a stationary state in a hydrogen atom satisfies the time independent Schrödinger Equation,
H|Φ> =En |Φ>,                                                           (1a)
where the Hamiltonian operator is the following
= p2/2m – k e2/R2,                                                       (1b)
and En
are its eigenvalues.

Preliminaries Old Quantum Theory Solution

Bohr solved for the energy eigenvalues as follows by using the old quantum theory, which are negative because they represent bound states of the electron:
En = –½ (k e2)2me2 (1/n2).                                        (1c)
[Free electron states are represented by positive values.]
The constant k is given by
k=1/(4 πε0).
The expression for energy eigenvalues can be simplified by writing it as
En = -½ k e2/a0 (1/n2),                                               (1d)
where  a0 is known as the Bohr Radius
a0 = ħ/(k e2 me)
a0 = 4πε0ħ/(e2 me).                                                     (1e)
Modern quantum theory finds the same energy eigenvalues, but finds that there are two more eigenvalues associated with the angular momentum of the orbital electron.

Eigenvalues and Casimir Operators

Casimir operators are those of the group of operators that commute among themselves, that is their order of multiplication does not matter. Commutation relations play an important role in modern quantum theory. They are a calculational tools that reduce much of the  mathematics of quantum mechanics to algebra.

Commutation relations are written for operators, some of which do not necessarily commute, that is, their product depends on the order that it is taken, such as in the following, in which the commutator is defined as
[A, B]= A B – B A.                                                      (2a)
If order does not matter, such as between X and Y, then
their commutator is zero.

Note that we write quantum operators in upper case, bold letters.
The following is a very useful theorem that we will use repeatedly
[A B, C]= [A, C] B + A [B, C].                                  (2b)
Note that the proof falls out just by writing out the operators
A B C A B = ACB  C A B + A B A C B.
Another almost trivial theorem is left to the reader to verify
[A, B]=[B, A].                                                          (2c)

Heisenberg's Uncertainty relations apply to expectation values of noncommuting opertors.

If operators, such as X and P, fail to commute, their values cannot be defined simultaneously with arbitrary precision. For example, Heisenberg's Uncertainty relations apply to components of position and corresponding components of momentum

Δx Δpxħ/2 ,                                                           (3a)
Δy Δpy ≥ ħ/2 ,                                                          (3b)
Δz Δpz ≥ ħ/2 ,                                                          (3c)

where the 'Δ's represent the root mean square (RMS) expectation values defined for x, y and z as

Δx2 = <Φ| (X <Φ|X|Φ>)2|Φ>,                            (4a)
Δy2 = <Φ| (Y – <Φ||Φ>)2|Φ>,                            (4b)
Δz2 = <Φ| (Z <Φ|Z|Φ>)2|Φ>.                            (4c)

Each pair of conjugate variables, such as coordinates and momentum components, are represented by operators that have the commutation relations:

[X, Px] = i ħ,                            (5a)
[Y, Py] = i ħ,                            (5b)
[Z, Pz] = i ħ.                            (5c)

The Heisenberg Uncertainty relations between them result mathematically from the failure of their representative operators to commute, which is not shown here, but perhaps in a later blog.

All of the quantities of a quantum mechanical system that represent variables that can be measured simultaneously with arbitrary precision are represented by operators that commute among themselves, i. e., they are Casimir operators.

Some further properties involving commutators

Commutators can be used instead of derivatives, thereby reducing a lot of calculation in quantum physics to algebra. For example, (5a) is just the beginning of a series of derivatives that can be taken of the various powers and functions of X:

[X2 , Px] = i ħ(2X) = i ħ d(X2)/dX                                 (6a)
[Xn, Px] = i ħ nXn-1  =i ħ dXn/dX                                    (6b)
[f(X), Px] = i ħ df(X)/dX.                                              (6c)                                            
Note that (6a) follows directly by applying (2b) then (5a). The reader can prove (6b) by mathematical induction, and (6c) by expanding f(X) in a Taylor Series in powers of X.

Identifying the Casimir operators for the hydrogen atom

The energy operator, or Hamiltonian, is a candidate for one of the Casimir operators for the hydrogen atom. What are some other possibilities? A good place to look is at angular momentum, which as an operator is defined as follows
L = R x P.
The components of the angular momentum vector components (Lx, Ly, Lz) are each the following operators:
Lx= YPzZPy,                                                          (7a)
Ly= ZPxXPz,                                                          (7b)
Lz= XPyYPx .                                                         (7c)

Note that the components defined in (7a), (7b) and (7c) do ont commute among themselves, which we will show for LxLy:
[Lx, Ly] = Y[Pz, Z] Px+ PyX[Z, Pz]   ;
[Lx, Ly] = -i ħ YPx+ i ħ PyX              ;
[Lx, Ly] = i ħ [XPy-YPx]                    ;
[Lx, Ly] = i ħ Lz  .                                                      (8a)

 It is left as an exercise for the reader to derive the following two commutation relations:
[Ly, Lz] = i ħ Lx ,                                                     (8b)
[Lz, Lx] = i ħ Ly .                                                     (8c)

Note that each component commutes with itself,
[Lx, Lx] =0  ,                                                            (8d)
[Ly, Ly]= 0  ,                                                            (8e)
[Lz, Lz]= 0  .                                                            (8f)
It is left as an exercise for the reader to show why.

Nest we show that Lz commutes with L2, where
L2 = Lx2+Ly2+ Lz2.                                               (9a)
This is easy, by showing that Lz  commutes with the square of each component:
z, Lx2] = [Lz, Lx] Lx + Lx [Lz, Lx] ;                    (9b)
[Lz, Lx2] = i ħ (Ly Lx + Lx Ly)  ;                            (9c)
[Lz, Ly2] = [Lz, Ly] Ly + Ly [Ly, Lx]    ;                 (9d)
[Lz, Ly2] = -i ħ ( Lx Ly + Ly Lx );
[Lz, Lz2] = 0 .                                                         (9e)
It follows that
[Lz, L2] =0 .                                                            (9f)

Finally, we will show that the radius of the orbit (R) commutes with all angular momentum components:
R= √[X2 + Y2 +Z2] ;                                                               (10a)
R/∂X = X/R ;                                                                          (10b)
R/∂Y = Y/R ;                                                                          (10c)
R/∂Z = Z/R ;                                                                          (10d)
[Lz,f(R)] = X [Py,f(R)] Y [Px, f(R) ] ;                                    (10e)
[Lz,f(R)] = – i ħ X Y / R df(R)/dR + iħ Y X / R df(R)/dR ;       (10f)
[Lz,f(R)] = 0 ;                                                                          (10g)

It is left as an exercise for the reader to show that
[Lx,f(R)] = 0 ;                                                                          (10h)
[Ly,f(R)] = 0 ;                                                                          (10i)
It follows that
[f(R), L2] =0.                                                                           (10j)
Finally, we show that P2 commutes with every component of angular momentum
[Lz, Px2] = 2 [X,Px]Px Py                                                        (11a)
[Lz, Pz] = 0                                                                              (11b)
[Lz, Pz2] = 0                                                                            (11c)
[Lz, Px2] = 2 iħ Px Py                                                              (11d)
[Lz, Py2] = [XPy,Py2] – [Y,Py2]Px                                           (11e)
[Lz, Py2] = 0 -2 iħ Py Px                                                          (11f)
[Lz, Px2 + Py2] = 2 iħ [Px ,Py].                                              (11g)
It follows that
[Lz, P2] = 0.                                                                             (11h)
 It is left as an exercise for the reader t show that
[Lx, P2] = 0                                                                              (11j)
Ly, P2] = 0.                                                                             (11k)

Using the form of the Hamiltonian (1b) it follows that
[Lx, H] = 0 ;                                                                             (12a)
[Ly, H] = 0 ;                                                                             (12b)
[Lz, H] = 0 ;                                                                             (12c)
and therefore
[L2, H] = 0.                                                                              (12d)

The Casimir operators for the hydrogen atom are therefore  H, L2, and Lz . The choice of the latter component of angular momentum is arbitrary, but by convention it is taken as Lz .

In the next blog in this series on the hydrogen atom, we will derive solutions for three simultaneous eigenvalues of the Casimir operators:
H|Φ> =En |Φ>,                                                                       (12e)
L2|Φ> =α  |Φ>
Lz|Φ> =β  |Φ>.                                                                      (12g)






This entry was posted in algebra, mathematics, quantum physics, State Vectors. Bookmark the permalink.

2 Responses to The Hydrogen Atom Analyzed Using Commutator Algebra I

Leave a Reply

Your email address will not be published. Required fields are marked *