Many-electron atoms

Revised notes, for the exercises class Struttura della Materia delivered for the first time on Nov. 23 & Dec. 03, 2001.

The following exercises are related to general properties of the many-electron atoms / ions, especially ground-state and magnetic properties.

They must be solved based on some (usually approximate) quantum-mechanical treatment (Hartree-Fock approximation to Schrödinger's wave equation, and consequent effective-Z approximation...).

The nuclear charge is assumed to be Z qe.

a0 = ℏ²/(me e²) is the Bohr radius and
EHa = e²/a0 is the Hartree energy.
Here me is the electron mass and e² = qe² /(4πε0) is the electromagnetic characteristic coupling constant.

For many-electron atoms reduced-mass corrections are irrelevant and conceptually meaningless.

Ground-state symmetry and multiplet excitations of partly-filled degenerate shells

When a number of electrons occupy a shell, but is insufficient to fill it, the quantum state is degenerate in terms of one-particle energetics. For example, one electron can occupy a p shell in 6 possible ways (quantum states) (including spin), while two (identical) electrons can occupy a p shell in 5×6/2=15 possible ways. This large degeneracy is (partly) lifted by the electron-electron interactions and by spin-orbit.

This leaves us with 2 problems: (i) determining which are the resulting states (ii) which is their energy ordering.

Determination of the multiplet states

To identify which states are actually present, one must account for spherical symmetry and the Pauli principle. It turns out that for not-too-large-Z atoms/ions one should first of all combine together the individual orbital and spin angular momenta of the electrons into total L and S, and then couple these together (via spin-orbit) to a total angular momentum J. Doing this in a way to respect Pauli's principle (antisymmetry of the total wavefunction with respect to exchange of any two electrons) is a tricky business.

The way to do it properly is to construct a table of all the states, labeled by their individual ml and ms. For example for two p electrons:

ml(1) ms(1)ml(2) ms(2)ML MS
|1>-1-1-20
|2>-10-1-1
|3>-10-10
|4>-10-10
|5>-10-11
|6>-1+10-1
|7>-1+100
|8>-1+100
|9>-1+101
|10>0000
|11>011-1
|12>0110
|13>0110
|14>0111
|15>1120
Note that a state like
ml(1) ms(1)ml(2) ms(2)ML MS
|1'>-1-1-20
is not listed, since state |1> is forcefully orbitally symmetric, therefore the state listed as |1> actually is a shorthand for the antisymmetric combination (|1>-|1'>)/√2. Identical considerations hold for states |10> and |15>. Reversing spin and orbital symmetry, we have analogous considerations for states |5>, |6>, |9>, |11>, and |14>.

Next, one starts with the state with maximum value of ML=2 (here, state |15>). One notes that it is orbitally symmetric under exchange of electrons 1 and 2. This state must therefore be a spin singlet (spin-antisymmetric). In addition, this state is the largest-ML component of a L=2 multiplet, made up therefore of 5 components: |L=2,ML=2>, |L=2,ML=1>, |L=2,ML=0>, |L=2,ML=-1> and |L=2,ML=-2>. The first state, as we already said, is state |15> in the table. The next one must have the same space symmetry, thus it must be an orbitally-symmetric (and spin-antisymmetric) combination of states |12> and |13>. Then |L=2,ML=0> must be again an orbitally-symmetric (and spin-antisymmetric) combination of states |7> and |8>, with some admixture of |10>, and so on. These 5 states are indicated by the symmetry label 1D2.

After locating these first 5 states, we remove them from the table (for example by drawing a bar on them), and are left with 10 states. Again we look for the largest-ML state among those left: we see that we have just one ML=1 state, the orbitally-antisymmetric (and spin-symmetric) combination of |12> and |13>.

As it is spin-symmetric it has to be a spin triplet (S=1) state, the MS=0 component of three states, two others of which are: |14> (MS=1) and |11> (MS=-1). These three states form altogether the ML=1 component of a L=1 orbital triplet. The ML=-1 component is clearly formed by: the orbitally-antisymmetric (and spin-symmetric) combination of |3> and |4> (MS=0 component), plus states |5> (MS=1) and |2> (MS=-1). The ML=0 component has a MS=0 component given by a suitable mix of |7> and |8> (no admixture of |10>... why?). The MS=1 and MS=-1 components are clearly given by |9> and |6> respectively. We have in total 3×3=9 states making up a 3PJ=0,1,2 multiplet.

15-5-9 = 1 state is left out of these cancellations: it is the remaining orbitally symmetric and spin antisymmetric combination of |7>, |8>, and |10>, the one orthogonal to the ML=0 component of the 1D2 state. This combination has only one ML=0 and MS=0 component, it is therefore a nondegenerate spin and orbital singlet 1S0 state.

In summary, we have classified the 15 states in the table into:

Ordering of the multiplet states

Now that we have organized all the states according to the total L and S (and J), the second part of the problem is to establish their energy ordering. Here Hund's rules come into play:
  1. States with higher total S "feel" a smaller electron-electron Coulomb repulsion.
  2. For a given S, states with higher total L "feel" a smaller electron-electron Coulomb repulsion.
  3. For a given S and L, spin-orbit splits states with different J. Even if the spin-orbit parameter ξ (see below) for each electron is always positive, the effective ξ acting on total L and S is positive for less than half-filled shell and negative for more than half-filled shell.
This is all we need to order up the 5 states 3P0, 3P1, 3P2, 1D2, 1S0. According to these rules, applied in sequence, the triplets 3P0/1/2 are lowest, the 2 singlets are higher up by some eV's typically. Then, between the 2 singlets, the 1D2 is below the 1S0 (typically by some eV or fraction of eV). Finally, the 3 spin-triplets, are ordered for increasing energy with increasing J (remember that L·S/ℏ2 = [J(J+1)-L(L+1)-S(S+1)]/2, as discussed below). The magnitude of these J-splittings depends on the spin-orbit ξ, which is small for low-Z atoms, and grows roughly with Z4. (For large-Z atoms the spin-orbit becomes the dominant interaction and all the L-S coupling described here loses its meaning. One must use the J-J coupling scheme...)
multiplet spectrum of two p electrons

In summary, the multiplet ordering for increasing energy results as draw in the figure at the right.

A similar procedure can be carried out, with more or less pain and trouble for any shell and any number of electrons.

When electrons belong to different shells, all the multiplet states are found by listing all those for each shell, and making the direct product, with the usual combination rules for angular momentum. For es. a p5d1 configuration yields all multiplets with total L=1,2,3 and total S=0,1 [for a total of 6×10 = 60 = (3+5+7)×(1+3) states].

Like in the single-shell case, the first two of Hund's rules provide useful indications for the ordering of the states. The J-ordering is decided case by case by which of the several spin-orbit couplings prevails.

Note: the ground-state symmetry of all the elements is reported in this colorful periodic table.

Problems:

0
  1. Redo a similar counting and sorting of states for p4.
  2. Redo a similar counting and sorting of states for d2.
  3. Redo a similar counting and sorting of states for p3.
  4. Sort out all the states for p2f13.
  5. Determine all the spherical-symmetry states in the Nb (Z=41) ground state electronic configuration 4d45s, and identify the actual ground state among them.
    RESULT: 6D½

Magnetic moments and coupling to an external static magnetic field

The orbital motion and spin of each electron carry magnetic moments
µL = -gL µB L/ℏ and µs = -gS µB L/ℏ
respectively. The universal factor µB = qe ℏ /(2 me) = 9.27402·10-24 A·m2 = 9.27402·10-24 J/T = 57.8838 µeV/T is called Bohr magneton. The "gyro-magnetic factors" for orbital and spin moments are approximately gL=1 and gS=2.

The magnetic moment of an orbitally nondegenerate state (L=0) (e.g. an half-filled atomic shell as in the ground state of the alkali elements, of N, of Mn, and analogous configurations) is purely of spin origin, thus it has a g-factor g=gS=2, irrespective of the number of electron spins cooperating to build up the total spin of that atom/ion.

Likewise, the magnetic moment of a spin-nondegenerate state (S=0) is purely of orbital origin, therefore it has a g-factor g=gL=1, irrespective of the number of electron orbital motions cooperating to build up the total L of that atom/ion.

A nondegenerate state with total L=0 and S=0 (e.g. a "closed" atomic shell) has null magnetic moment, thus no first-order coupling to an external magnetic field.

In many atomic configurations one can have nonzero both L and S. In such case, even in the absence of external magnetic field, the full degeneracy (2 L+1)·(2 S+1) is partly lifted by the spin-orbit interaction, of the form:

HSO = ξ L·S/ℏ2
[here ξ is the effective spin-orbit coupling introduced above]
Since
L·S = ½(J2-L2-S2),
in the absence of external fields, the total angular momentum J2 commutes with the Hamiltonian, has eigenvalue ℏ2J(J+1), and J is a "good quantum number". This fact reflects the invariance of the atom under global rotations, rather than under separate rotations in the spin and in the position spaces, as could occur in the absence of spin-orbit coupling.

In concrete, the eigenstates of the spin-orbit interaction are simultaneous eigenstates of the 4 commuting operators: L2, S2, J2, and Jz. The (2 L+1)·(2 S+1) states are split into (L+S-|L-S|+1) multiplets of states, each labeled |L,S,J,MJ>, with MJ (the quantum number associated to Jz).

For example, HSO splits the 6 spin-orbital states associated to the 1s²2s²2p configuration of the boron atom into the two multiplets, labeled |L,S,J,MJ>=|1,½,J,MJ>. As J=½ and J=3/2 are the two allowed values, the 6 spin-orbital states are:

|1,½,½,½> and |1,½,½,-½>
|1,½,3/2,3/2>, |1,½,3/2,½>, |1,½,3/2,-½>, and |1,½,3/2,-3/2> ,
separated in energy by 3/2 ξ = 3/2 1.2635 meV= 1.8953 meV, for the boron atom.

The coupling with an external magnetic field B leads to an additional term in the Hamiltonian:

HB = -µtot·B = -(µL+µSB = µB (gL L/ℏ + gS S/ℏ)·B = µB (gL Lz/ℏ + gS Sz/ℏ)Bz
(assuming that the z axis points in the direction of the B field). Unfortunately, the operator J2 does not commute individually with neither Lz nor Sz (it would commute with Jz=Lz+Sz, but not with the "unbalanced" combination in the above equation). For this reason, as soon as a magnetic field is turned on, the basis of spin-orbit eigenstates is no more a basis of actual eigenstates for the total Hamiltonian H=HSO+HB.

It turns out that, for very weak field the spin-orbit-split states are still fairly good approximations to the real eigenstates (people say that J and MJ are approximate quantum numbers. In the limit of weak field, one can show that the effect of the magnetic field is to split the (2 J+1) degeneracy of each J-multiplet according to:

Heff = -µeff z(L,S,J)Bz = µB gLandé(L,S,J) (Jz/ℏ) Bz
where the effective magnetic moment component of a J-multiplet is given in terms of the dimensionless Landé g-factor
gLandé(L,S,J) = [3 J(J+1) +S(S+1) -L(L+1)]/[2 J(J+1)],
and as usual it is assumed that the z axis is oriented along the direction of the external magnetic field B.

In concrete, in the weak-field limit, the main splitting between the levels is given by the spin-orbit term HSO = ξ L·S/ℏ2 = ξ[J(J+1)-L(L+1)-S(S+1)]/2. Each of these [initially (2 J+1)-degenerate] levels is then split according to ΔE(MJ)= µB gLandé(L,S,J) MJBz.

Note that the above Landé relation has nothing to tell us about the length of the total magnetic moment vector! The latter can be worked out to:

|µtot| = µB (2J(J+1) - L(L+1) + 2S(S+1))½.
The derivation of this useful relation is a recommended exercise.

In the opposite limit of strong field, the main structure of the spectrum is given by the HB. The eigenstates are approximately those of the ξ=0 no-spin-orbit limit. They are labeled by |L,S,ML,MS>. In this limit of "sufficiently large" field, spin-orbit effects can conveniently be ignored, or, if higher accuracy is needed, treated within perturbation theory.

Therefore, the magnetic energies are trivially given by Emagn(ML,MS)= µB (gLML+gSMS)·Bz

For an intermediate field, the spectrum exhibits an irregular pattern, which interpolates continuously between the two limits.

The spin-orbit parameter ξ

For each electronic state, an effective spin-orbit parameter ξ determines the zero-field splitting of the different-J states, through HSO.

If the electronic state involves a single electron (as, for example, for the optical-electron states of alkali metal atoms or for core hole states), then a modified 1-electron formula gives an useful approximation:

ξ = EHa (Zeff-SO)4 α2 [n3 l (l+½) (l+1)]-1
where the Hartree energy EHa =2 ERy = qe2/(4 πε0 a0), α is the fine-structure constant characteristic of relativistic corrections in atomic physics and l and n are the one-particle quantum numbers of the 1-electron state considered. The above formula is exact (with Zeff-SO=Z) for 1-electron atoms, and may be derived (EXERCISE!) from the relativistically correct expression of the eigenenergies of the 1-electron atom.

The effective Zeff-SO for many-electron atoms is different from the Zeff that gives approximate expressions for the level positions, and usually Zeff<Zeff-SO<Z (EXERCISE: explain the reasons for these inequalities). Both Zeff and Zeff-SO depend on n, and for small n they are close to Z, while in the limit of large n they both converge to the residual core charge (Z-1 for neutral atoms).

The universal constant α2 EHa = 1.44904 meV sets the scale of the spin-orbit interaction. Clearly, the fourth-power Zeff-SO-dependence affects the observed level splitting quite strongly. For example, the splitting of the two 2p states of H is 3/2 ξ(n=2)= 3/2 α2 EHa/48 = 0.0453 meV, while for Cs (Z=55) a separation of 44.2 meV between the J=5/2 and J=7/2 states of the lowest 4F excitation is observed, corresponding to a Zeff-SO=14.7.

For electronic configurations involving several electrons, the spin-orbit operator acts as the sum of the individual spin-orbit couplings for each electron. In LS coupling this can be seen as an effective "total" spin-orbit coupling the total angular momenta L and S. The precise value of the effective spin-orbit parameter multiplying the operator L·S can be deduced case by case. It may even become negative for a more-than-half-filled shell.

Problems:

  1. How does a 1D2 atomic state of Ca split under the action of a magnetic field of 5 T? How would this splitting change for a 1D2 of C? How would this splitting change for a magnetic field of .05 T?
  2. How can one tell if a magnetic field is weak or strong?
    RESULT: Compare the typical magnetic energy µB|B| and the typical spin-orbit energy ξ.
  3. Draw the spectrum of the 2p levels in H for |B|=0.1 T. Compute all the splittings and the absolute energies with respect to the ionization threshold. Estimate the error introduced by the weak/strong field approximation made.
  4. Draw the spectrum of the 2p levels in H for |B|=25 T. Compute all the splittings and the absolute energies with respect to dissociation. Estimate the error introduced by the weak/strong field approximation made.
  5. Remembering the Z4 dependence of ξ, draw the spectrum of the 2p levels in Be3+ for |B|=25 T. Compute all the splittings. Estimate the error introduced by the weak/strong field approximation made.
  6. The ground state of Er is [Xe]4f125d06s2, of symmetry 3H6. Knowing that the spin-orbit parameter ξ=-0.1 eV, a) for the first two excited spin-orbit states determine the symmetry and excitation energy in the absence of external fields; b) draw the spectrum in presence of a magnetic field of 3 T, and determine all the splittings within each of the three J-multiplets.
    RESULT: L=5, S=1, thus the two excitations are 3H5 and 3H4. For the ground state and the two excited states the values of L·S/ℏ2 = [J(J+1)-S(S+1)-L(L+1)]/2 are:
    3H6 --> 5
    3H5 --> -1
    3H4 --> -6
    Consequently, the excitation energies with respect to the ground state are estimated 0.6 and 1.1 eV respectively.
    The g-factors giving the splittings within these three multiplets are:
    3H6 --> 7/6
    3H5 --> 31/30
    3H4 --> 4/5
    Note: a g-factor smaller than 1, obtained from both gL and gS≥1 deserves meditation...
  7. In a Stern-Gerlach experiment an atomic Se beam is sent across a 10 cm long magnet with a field gradient of 230 T/m. How many "spots" are observed on a screen placed at 20 cm from the end of the magnet? What is the distance between the maximally-separated spots, assuming that the kinetic energy of the incoming beam is 0.040 eV?
    RESULT: About 5 mm.
  8. Same as previous problem, replacing Se with V.

X-ray, optical spectroscopies.

Specific pages are devoted to "high-energy (core)" and "low-energy (valence)" excitations, and corresponding X-ray and optical spectroscopies.
created: 26 Nov 2001
last modified: 15 Jun 2026
by Nicola Manini