Lecture 13-Atoms with many electrons

After our discussion of extremely simple atoms like hydrogen and helium, we will now discuss the most important properties of more complex atoms. We will see, how we can categorize them and discuss some of the general properties We started the discussion of atoms in lecture 3 by an extremely simple and powerful model, the two-level system (Jendrzejewski et al., a). We then moved on to discuss how it can emerge within the hydrogen atom and the Helium atom. For both of these we dived again into simplified schemes. Especially for Hydrogen, we saw in lecture 7 the exploding complexity of the models as we tried to describe it (Jendrzejewski et al., b). If we now want to leave these rather academic problems and turn to the other widely used atoms, molecules etc, we have towards effective models for two reasons: • Analytical solutions do not exist. • Full numerical solutions become extremely and expansive. A particularly instructive discussion of the ineffeciency of brute force numerical methods was given by Kenneth Wilson, the father of the renormalization group, in Sec 3 of (Wilson, 1975). Figure 1: A toy model, exemplifying the problem of direct numerical analysis The main idea is plotted in Fig. 1 and its description in the words of Kenny Wilson can be found in Fig. 2.

We started the discussion of atoms in lecture 3 by an extremely simple and powerful model, the two-level system (Jendrzejewski et al., a).We then moved on to discuss how it can emerge within the hydrogen atom and the Helium atom.For both of these we dived again into simplified schemes.Especially for Hydrogen, we saw in lecture 7 the exploding complexity of the models as we tried to describe it (Jendrzejewski et al., b).If we now want to leave these rather academic problems and turn to the other widely used atoms, molecules etc, we have towards effective models for two reasons: • Analytical solutions do not exist.
• Full numerical solutions become extremely and expansive.
A particularly instructive discussion of the ineffeciency of brute force numerical methods was given by Kenneth Wilson, the father of the renormalization group, in Sec 3 of (Wilson, 1975).The main idea is plotted in Fig. 1 and its description in the words of Kenny Wilson can be found in Fig. 2.
Figure 2: On the inefficiency of brute force numerics as described in Sec 3 of (Wilson, 1975) We will come back to different experimental approaches to solve problems for which we cannot devise effective theories at a later stage, when we discuss quantum simulation and quantum computation with atomic systems.But first we will try to gain a good understanding of atoms with many electrons.

Atoms with Many Electrons
As already discussed before, the N electron system cannot be solved in its full complexity, so we have to walk through the different levels of physical effects to understand what is going on.We will start out with the simplest model, which consists of

• N electrons without spin
• Bound to the point-like nucleus of charge Z, which is supposed to be not moving.
In natural units, the Hamiltonian can be written as: The system is visualized in Fig. 3.As before the Hamiltonian has two contributions: • The well understood binding of each electron to the nucleus, which allows us to make the connection with the hydrogen orbitals.
• The interaction between the electrons, which couples all the orbitals.
How can we deal with the 1/r ij terms?If many electrons are contained in the system we cannot ignore the influence of the electron-electron interaction on the system, which was already a questionable approximation in the Helium atom (Jendrzejewski et al., c).

Central Field Approximation
But we can bring back the system to the well-known problems with spherical symmetry by splitting the interaction into two parts.One will be assumed to have spherical symmetry and will be written as i S(r i ) and the rest will be treated as a perturbation.So we obtain the Hamiltonian: The perturbation is then difference between the spherically symmetric part of the interaction and the true electron-electron interaction: S(r i ) describes now the screening of the nucleus due to the other electrons.This function interpolates between (see 4): • For very small distances the electron mainly sees the strong attraction of the nucleus.
• For very large distances the other electrons mostly screen the charge of the nucleus and electron feels the 1/r potential of the hydrogen atom.
Figure 4: The screening as a function of distance.
We can now focus on the study of the part with spherical symmetry Ĥs .Given its symmetry, we can once again use the conservation of orbital angular momentum (see Lecture 5 (Jendrzejewski et al., d)).As a result, the wave function can be a product state of single-particle wave functions: where ψ i ( r i ) are solutions of the Schrödinger equation: ψ i can be split up into a radial and an angular part: For the radial part we get u(r) = r • R(r).Thus, u(r) solves the radial Schrödinger equation: The major difference to the H atom is that the degeneracy of the l levels is lifted because of V cf .

On the determination of the effective potential
We will not discuss here the different techniques to derive the central field in detail as this is quickly diving into the different numerical techniques of many-body systems.A first discussion of the different approaches can be found in chapter 10 of (Hertel and Schulz, 2015a).The different levels of sophistication are: • The Thomas-Fermi model, which assumes that the N electrons behave like a Fermi gas inside the Z/r potential.
• The Hartree method, in which the effective potential is iteratively recalculated based on the obtained solutions to the radial Schrödinger equation.
• The Hartree Fock method, in which we also take into account the proper symmetrization of the wavefunctions.

Filling up the shells
If we can ignore the spin for the determination of the energy levels, we have the following quantum numbers: • n, which is electron shell.
• l, which is the orbital angular momentum with l < n.
• s, which is the spin of the electron and it can be ±1/2 As the electrons are fermions we can fill up each of the nl states with two of them.We then write down the configuration of the electron by writing down the numbers of electrons per nl configuration.So we get the periodic table shown in Fig. 5.We will typically use the following notation: • A configuration is the distribution of the electrons over the different orbits.
• Electrons with the same n are part of the same shell.
• Electrons with the same n and l are part of the same sub-shell.
• The inner shells are typically filled and form the core.
• The outermost shell is typically named the valence shell.
Figure 5: Periodic table as taken from (com)

Alkali Atoms
Alkali atoms are the simplest to understand and widely used in the field of laser cooling.The electron configurations are Examples.
This structure can be nicely understood by the idea of screening as introduced in the central field approximation.Consider for example the Na for which one electron has to be on the third shell.Within the hydrogen atom the 3s, 3p and 3d are degenerate, however the screening will lift this degeneracy.
• Electrons in the s-shell spend a lot of time close to the nucleus and fill a strong binding potential, i.e. they have a low energy.
• On the other extreme the electrons in the d-shell have large angular momentum, which hinders them from getting close to the core.The 10 electrons can therefore efficiently screen the nucleus from the electron.The energy of the d shell is therefore quite close to the energy of the hydrogen atom.
For K the same principle applies.
• The 4s orbital is strongly shifted down to the strong nuclear potential.
• The 3d orbital is well screened by the other electrons, such that it is of higher energy than the 4s orbit.
We can describe this effect empirically through the quantum defect.We simply write down the energy levels as: where δ l is the quantum defect for a certain value of l.Some examples are summarized in Fig. 6.The effect of the quantum defect is actually so substantial that it introduces optical transitions between the different subshells of Alkali atoms.
• For Li there is a 2s → 2p transition at 671 nm.
• For Na there is a 3s → 3p transition at 589 nm.
• For K there is a 4s → 4p transition at 767 nm and 770nm.
• For Rb there is a 5s → 5p transition at 780 nm and at 795nm.
• For Cs there is a 6s → 6p transition at 852 nm and at 894nm.
What is the origin of these doublets ?The scaling of the splitting with the nucleus indicates relativistic origins and the splitting is indeed due to spin-orbit coupling, which we will discuss in the next lecture (sys).

Figure 1 :
Figure 1: A toy model, exemplifying the problem of direct numerical analysis

Figure 3 :
Figure 3: Charge distribution of three electrons bound to the nucleus.