The Central Force Problem: Hydrogen Atom B. Ramachandran Separation of Variables The Schrödinger equation for an atomic system with Z protons in the nucleus and one electron outside is h µ Ze ψ = Eψ, r 4πε where ε is the permittivity of vacuum, which is a constant having the value 8.85487 J C m 4πε =.65 J C m, and µ is the reduced mass of the nucleus-electron system. The quantity appearing above is the Laplacian which, in Cartesian coordinates, is given by = = x + y + z. Transforming the Laplacian to spherical polar coordinates, we get We recognize that h µ = h µ r r r r + sin θ θ sin θ θ + sin θ φ. L = h sin θ θ sin θ θ + sin θ With this substitution, and making the well-justified assumption that µ m e, the mass of the electron, Eq. becomes [ h m e r r r r + L m e r Ze 4πε r E ] φ. ψr, θ, φ =. 3 Thepresenceofthe L operator suggests that the angular part of the solution, i.e., the part that depends on the angles, are the spherical harmonics Yl m θ, φ. Therefore, we choose the following function to attempt a separation of variables: ψr, θ, φ =RrYl m θ, φ. 4 Substituting Eq. 4 into Eq. 3 and simplifying, we get R r + r [ R r + m e h E + Ze 4πε r ll + r ] R =. 5 The notation can be further simplified by defining the following dimensionless parameters: α = m ee h,β = Ze m e 4πε h α.
Note that α will be positive only for negative values of E. This means that the parameter α will become complex for positive values of E, thereby restricting real solutions to those with negative energy. Since zero potential energy corresponds to infinite separation of the charges Ze and e, the negative energies correspond to the bound states of the atom. Using these definitions, Eq. 5 can be simplified to R r + r R [ α αβ + r r ] ll + R =. 6 r Solving the Radial Equation We could now attempt a series solution of Eq. 6 as was done in the case of the harmonic oscillator. However, this would lead to a three-term recursion relationship which is much harder to deal with than a two-term one. Therefore, we define a new variable ρ =αr and express Eq. 6 in terms of this variable as T ρ + ρ where T ρ is simply Rr expressed in terms of the new variable. T [ β ] ll + + ρ 4 ρ ρ T =, 7 Let us examine the asymptotic behavior of the solution. As ρ, Eq. 7 reduces to T ρ 4 T =, which has the solution T ρ =exp±ρ/. The solution corresponding to the positive exponent grows exponentially as ρ increases. This makes it impossible for the solution to be square-integrable. Therefore, we choose the solution corresponding to the negative exponent. Let us now assume that the general solution will be of the form T ρ =Gρexp ρ/. 8 Substituting Eq. 8 into Eq. 7 and simplifying yields ρ G ρ+ρ ρ G ρ+{β ρ ll +} Gρ =. 9 This equation must be valid for all values of ρ. Therefore, at ρ =,wehave ll +Gρ =; l =,,,... Since this equation has to be satisfied for nonzero values of l, this suggests that Gρ must have the form ρ s Lρ, withs>. We now expand Lρ in a power series as Lρ = a j ρ j. j= Substituting this expression for Gρ into Eq. 9 and collecting terms with the same power of ρ, weget the following recursion relationship for the coefficients a j : a j+ a j = l + j + β j +j +l +.
Now, we must verify whether the series will converge as j. From Eq., it is clear that a j+ Lim = j a j j. To understand what such a series might resemble, consider a Taylor expansion of the function e ρ : e ρ ρ j = j! = b j ρ j, j= j= b j+ j! Lim = j b j j +! = for large j. j + j 3 Truncation of the Series: Quantized Energies This suggests that Lρ, and therefore Gρ, behaves like e ρ asymptotically. This, of course, leads to a non-square-integrable solution, and is also in contradiction to the conclusion from Eq. 8 that T ρ behaves like e ρ/ as ρ. The way out of this difficulty, of course, is to truncate the series at some value of j as in the case of the harmonic oscillator, say j = k, sothat a k+ a k = i.e., β = l + k +. l + k + β k +k +l + =, Since l + k + must be an integer, having possible values,, 3,..., wedefine n = l + k +so β = n; n =,,.... Also, since the truncation relationship must be valid for all values of k including k =, we see that n l + or l n. Thus, the series truncation also imposes a limit on the values l can assume. Now, substituting this into the definition of β, weget n Ze m e = 4πε h α, or n = Ze m e 4πε h Rearranging, and recognizing that the energies depend on the integer n, weget We now define a parameter a = 4πε h, so that m e e h m e E. E n = Z e 4 m e 3π ε ; n =,,.... 3 h n E n = e 4πε a Z n. 4 3
4 Atomic Units The parameter a has units of length verify this! and is equal to.5977 m. This parameter has a special place in quantum chemistry because it is the unit of length called the Bohr in atomic units au, which are a set of units commonly used in quantum chemistry calculations. In these units, the unit of mass is the mass of the electron i.e., m e =au, and h =so that for the case of the Hydrogen atom Z =, the ground state energy is E = /. The unit of energy is the quantity e /4πε a called the Hartree and so, we see that the ground state of the hydrogen atom is half of a Hartree below the ionization limit i.e., E =. From this information, we calculate that Hartree =7.5 ev =4.35979 8 J, where ev stands for the electron volt, a commonly used unit of energy in quantum chemistry, defined as ev = ec. Volt =.69 9 J. 5 The Rydberg Constant The energy expression of Eqs. 3 or 4 can be used to calculate the value of the Rydberg constant appearing in Bohr s original expression for the frequencies of the lines appearing in the emission spectrum of the hydrogen atom. Bohr s expression was λ = R n n ; n >n. Recognizing that E = hν = hc/λ, we get from Eq. 4, E = h Z m e a n n ; n >n, λ = h Z 4π cm e a n n. The value of the Rydberg constant R thus calculated is 9,677 cm. Modern experiments have established the correctness of this value to about 8 significant figures. 6 Atomic Orbitals It is obvious from Eq. 5 that the solutions of the radial equation depend on the value of l. Now, the series truncation condition above also brings in a quantum number n. The radial solutions are, therefore, written as Rnlr, and including normalization, are given by Rnlr = NnlGnl ρ e ρ/, { Z 3 } n l! Nnl = na n [n + l!] 3, G nl ρ = ρ l L l+ n+ ρ, 4
where the polynomials denoted by L s ρ are called the associated Laguerre polynomials of degree r s, r and are defined as L s r ρ = ds dρ s L rρ, where the L r ρ are the Laguerre polynomials given by L r ρ =e ρ dr dρ r ρ r e ρ. Substituting for α in the defintion of ρ, wegetρ = Zr/na. We now write the solutions to the hydrogen atom problem as ψ nlm r, θ, φ = R nl ry m l θ, φ, 5 n =,,..., l =,,..., n, m =, ±, ±,...,±l. The convention followed in labeling atomic states is to use a letter symbol to denote various values of l, so that we use Value of l : 3 4 5 Letter symbol: s p d f g h... Therefore, the state nlm = is referred to as the p state and so on. Each nlm state defines what is known as an atomic orbital. Let us now examine the nature of the wave functions. We define two types of probability functions for this purpose. The quantity ψ nlm r, θ, φ r dr sin θdθdφ gives the probability of finding the electron in a volume element enclosed by the limits r, r + dr; θ, θ + dθ; φ, φ + dφ. However, when dealing with a spherical species like an atom, it is more useful to average over the angles and get the radial distribution function, S nl r, which is obtained as S nl r = r R nl r dr = r R nl r dr. π π sin θdθ dφ Yl m θ, φ The final result comes about because the integrals over the angles yield unity because the spherical harmonics are normalized: π π sin θdθ dφ Yl m θ, φ =. 6 The physical significance of this is that r R nl r dr gives the probability for finding the electron on a spherical shell of radius r and thickness dr. We will examine the shapes of the radial functions and the two types of probability distribution functions using class handouts. We will also combine the radial functions with the spherical harmonics i.e., generate atomic orbitals and examine the nodal structure of these three-dimensional functions. Below, we tabulate a few solutions of the Schrödinger equation Eq.. For convenience, we define 5
σ = nρ/ =Zr/a. n l m ψ nlm r Common name 3/ Z e σ π a 4 π ± 3 Z 8 3π 3 3 ± 3/ Z a 3/ Z σe σ/ cos θ 4 π a 3/ Z σe σ/ sin θ cos φ 4 π a 3/ Z σe σ/ sin θ sin φ 4 π a a Z 8 π 3/ Z 8 π a Z 8 π 3 3 ± 3 ± s σe σ/ s p z p x p y 3/ 7 8σ +σ e σ/3 3s a 3/ σ6 σe σ/3 cos θ 3p z σ 6 σ e σ/3 sin θ cos φ 3px 3/ σ6 σe σ/3 sin θ sin φ 3py a 3/ Z σ e σ/3 3 cos θ 3d 8 6π z a 3/ Z σ e σ/3 sin θ cos θ cos φ 3d xz a 3/ Z σ e σ/3 sin θ cos θ sin φ 3d 8 yz π a 3/ Z σ e σ/3 sin θ cos φ 8 π a 3d x y 8 π 8 π Z a 3/ σ e σ/3 sin θ sin φ 3d xy 7 Most Probable and Average Values of r The most probable and average distance of the electron from the nucleus is of considerable interest in chemistry and physics. The most probable distance of the electron from the nucleus is given by the position of the highest maximum in S nl r. Recall that the derivative of a function goes to zero at its maximum or minimum. For example, in the case of the s orbital, we set the derivative of Sr = and get d dr σ e σ =σ e σ =, 7 where we have divided both sides by the multiplicative constants present in S r to eliminate them. Since the exponential term in the derivative is non-negative for all real values of σ, we conclude that the most probable value of r corresponds to σ =,orrmp = a /Z. For the hydrogen atom, this is exactly equal to the Bohr radius, which is the radius of the ground state orbit in Bohr s planetary model of the 6
hydrogen atom.the average value of r is given by the expectation value π π r nl = ψ nlm r, θ, φ rψ nlm r, θ, φr sin θdrdθdφ π π = R nl r rr nlrr dr Yl m θ, φ sin θdθdφ = r 3 Rnlr dr { [ = n a + ll +. Z n ]} The final result is obtained by a somewhat lengthy and clever use of the method of integration by parts. 8 Hydrogenic atoms in a magnetic field Note that the energy of a one-electron atom, given in Eq. 4, is independent of quantum numbers l and m. This means that states with the same value of n but different values of l and m have the same energy. Let us examine what happens to these states when a magnetic field is applied along the Z-axis. Consider a circular loop of conducting wire of radius r through which a charge Q is moving with velocity v. Since current is the charge flow per unit time, we get I = Qv πr. The value of the magnetic dipole µ associated with this current is given by IA, wherea is the area enclosed by the loop, i.e., A = πr. Therefore, Therefore, the magnetic dipole vector is given by For an electron, and µ = Qv πr = πr Qrp Qvr = m. µ = Q m r p = Q m L. µ = e m e L, 8 µ = e h m e ll +=µb where the quantity µb is a constant called the Bohr magneton. ll +, 9 Now let us consider a magnetic field B applied to the atom. between the field B and the magnetic dipole µ is given by The potential energy of interaction VB = µ B Yep. You guessed it. Homework assignment! 7
From Eq. 8, we get VB = e/mel B. If we assume that the field is oriented along the Z-direction, i.e., the components of the field along the x and y directions are zero, then, e VB = îlx + me ĵl y + kl kb z e = BL z. m e In operator terminology, therefore, V B = eb Lz. m e Therefore, the Schrödinger equation for a hydrogen atom in a magnetic field is Ĥ + VB R nl ryl m θ, φ =E n + µ B BmR nl ryl m θ, φ, where the Hamiltonian operator Ĥ corresponds to the Schrodinger equation we started with in Eq., E n is given by Eq. 4, and µ B is the Bohr magneton defined in Eq. 9. Therefore, in a magnetic field, states with nonzero values of m will aquire different energies. 8