Schrödinger's wave description of 1-electron atoms

Problems carried out the course Struttura della Materia I (originally for Struttura della Materia)

The following exercises focus on the 1-electron atom or ion. Their solution assumes the quantum-mechanical approach, based on Schrödinger's wave equation and its solutions.

The potential energy is assumed to be exactly Coulomb -Ze²/r (no corrections due to the finite size of the nuclear charge distribution). The nuclear charge is assumed to be +Z qe. a is the modified Bohr radius 4 πε0 ℏ²/(µ qe²) = a0 µ/me, where µ = meMn/(me+Mn) is the reduced mass of the 2-body electron-nucleus system and me is the electron mass. Fine-structure corrections (e.g. spin-orbit coupling) are to be neglected except where explicitly mentioned.

Radial part of the wavefunction:

  1. Normalize an exponentially-decaying radial wavefunction R(r)=N exp(-r/r0).
    RESULT: N=2 (r0)-3/2 [optionally multiplied by an arbitrary phase factor exp(iφ)]
    [Use has been made of the following integral: Ik= ∫0 dx xk·exp(-x) = k! ]
  2. mean value
    Find the best variational spherically-symmetric exponentially-decaying wavefunction for the 1-electron atom.
    In other words, take the normalized wavefunction of the previous exercise, and compute the expectation value of the H-like atom Hamiltonian on this state, assuming no angular dependence (so that the contributions of all the angular derivatives in the Laplacian vanish). This "mean value of the energy" depends on the parameter r0: determine for which r0 this mean energy is minimum.
    RESULT: One must find the exact ground-state wavefunction, i.e. r0=a/Z
    NOTE: this is an exceptional variational exercise where among the (infinitely) many "variational" wavefunctions one tries, there happens to be the exact ground state, which obviously turns out to be the lowest in energy. In "real life" this almost never occurs: the best variational wavefunction is usually only an approximate ground state.
  3. Compute the average (also called expectation value) of rn over the 1-electron atom ground state wavefunction, for integer n.
    RESULT: For n≥-2: <rn>=(a/2 Z)n·(2+n)!/2
  4. In a one-el. atom in its ground state, what is the probability to find the electron at a distance from the nucleus greater or equal to a given d? Evaluate this probability for d=1 m and d=1 µm.
    RESULT: Indicating with u=2 Z d/a, P(outside d)=exp(-u)(1+u+u2/2).
  5. [Written special test 2000-12-06] Indicate for which values of Z and of the quantum numbers (n,l,m) the radial distribution in a 1-electron atom is:

    P(r)=N r4 exp(-3 r/a)

    where N is the normalization constant. Compute N, <r>, <r-1> and the radius rM for which the radial distribution is maximum.
    The general structure of the radial wavefunction is:

    Rn,l(r)=Nn,l · (k r)l Ln+l2l+1(k r) exp(-k r),       with k=Z/(n a)

    where Lpq(x) is a (p-q)-degree polynomial in x with nonzero x0 term, and Nn,l is a constant.
    RESULT: l=1, n=2, Z=3, rM=4a/3

Angular part of the wavefunction:

  1. Show that the parity of the one-electron wavefunction ψ3,2,1 is (+1).
  2. Evaluate the probability to find the electron in an equatorial "conical" band of ±10° around the xy plane θ=π/2, assuming the atom is in the ψ3,2,2 orbital state.
    RESULT: P = [150 sin(10°) + 25 sin(30°) + 3 sin(50°)]/128 = 0.319104

Transition rates

(WARNING: long exercises!):
Note: in these exercises one must compute transition rates between degenerate states. To solve them, remember the thumb rule for total transition probabilities: average over the initial states and sum over the final states.
  1. Suppose you want to build a 1 W fluorescent bulb based on the 2p -> 1s transition in atomic H. Assume you would be able to place the H atoms well separated from each other so that the light emitted from one atom is unaffected by that emitted from the other atoms. Neglect all the fine and hyper-fine structure of H. How many 2p-excited atoms are needed at any given time?
    RESULTs: (In the numerics, we neglect the reduced-mass correction).
    Radial integration yields 128/81 sqrt(2/3) a/Z.
    Angular integration gives for the (1,m) -> (0,0) amplitude: 3-1/2.
    Total single-atom transition rate = (2/3)8 Z4 α3 EHa/ℏ = 6.2683×108 Z4 s-1.
    Number of atoms = 2187/32 ℏ Power/(α3 EHa² Z6) = 9.75791×108 atoms [at Power=1 W, Z=1].
  2. Compute the the spontaneous-emission rates of the np -> 1s transitions in 1-electron atoms, for n=2 (see ex. above), 3, 4 and 5.
    RESULT: In units of Z4 α3 EHa/ℏ = meqe10Z4 /[c3 (4πε0)56] = 1.6065×1010 Z4 s-1 one gets for n=2 to 5: 1/6561, 1/24576, 162/9765625, 40/4782969.
  3. Compute the ratio between the spontaneous-emission transition rates of the 4p -> 3d and the 4p -> 3s decays in H.
    RESULT: 4096/36125

Relativistic corrections

to Schrödinger's description of 1-electron atom:
  1. [Written special test 2000-12-06] In an atomic H gas, let the electron have a significant probability to occupy the levels n=1,2,3. Accounting for the fine-structure corrections: a) draw a level scheme; b) compute the splittings in eV between levels within the same-n multiplets; c) indicate, in the dipole approximation the number of emission lines present in the Balmer and Lyman series, and compute the maximum wavelength in both these series.
    [reminder: the 1-electron atom levels are expressed by

    En,j = En [1 + (Z α)2 / n (1/(j+1/2) - 3/(4n)) ]

    where En is the energy as given by the simple Z/r potential].
    RESULT: c) 121.566964 nm and 656.47268 nm
created: 14 Nov 2001
last modified: 8 Apr 2019
by Nicola Manini