Module 12 : 2-D and 3-D Quantum Problems I. Particle In A 2-D Box With Infinite Walls Consider a particle trapped in a 2-D infinite box as shown below: U (x) = y=b U (x) = 0 U (x) = 0 U (x) = y=0 x=0 x=L U (x) = 0 The Time Independent Schrodinger Equation for a multi-dimensional problem is given by 2 k 2 where k 2 2 m E U . 2 Using the Cartesian form of the Laplacian, we have inside the box 2 2 2mE 2 k 2 where k 2 2 2 x y and that is zero outside the box. Since each term of the potential energy depends on at most one spatial coordinate (in this case no coordinates), we can use separation of variables to simplify our work. Thus, we choose a product solution for the wave function where X and Y are each functions of only one spatial coordinate. Xx Yy Substituting into the Time Independent Schrodinger Equation, we have 2 X 2 Y Y 2 X 2 k 2 XY x y 1 2 X 1 2 Y 2mE 2 k constant X x 2 Y y 2 2 Each term depends on different spatial coordinates which are independent. Since the equation must hold for any point inside the box, each term must be a constant! Furthermore these constants must have the dimensions of inverse square meters (i.e. wave number squared). Thus, it makes sense to write the individual constants in a manner similar to what we did with the 1-D well as follows: 1 2 X 2 kx 2 X x 1 2 Y 2 ky 2 Y y We also note the additional requirement: k x k y k 2 2mE 2 2 This is just the dot product of a vector (wave number) upon itself where k k x î k y ĵ . k ky kx From our work on the 1-D infinite square well, we know the general solutions to these two equations are given by X A Sin k x x B Cos k x x where 0 x L Y D Sin k y y E Cos k y y where 0 y b Thus, the energy wave function is A D Sin k x x Sin k y y 0 x L and 0 y b 0 otherwise The boundary condition that is zero for all values of y when x = 0 requires that B = 0. The boundary condition that is zero for all values of x when y = 0 requires that E = 0. The boundary condition that that is zero for all values of y when x = L requires that k x L n π where n = 1, 2, 3, … The boundary condition that that is zero for all values of x when y = b requires that k y b q π where q = 1, 2, 3, … By plugging the results for the components of the wave number into our energy equations, we obtain the energy eigenvalues: 2 k x k y h2 kx ky E n,q 2m 8π2 m 2 2 2 2 h2 n2π2 q2π2 E n,q 2 2 2 8π m L b h c2 n 2 q2 n 1,2,3,.. E n,q where 8 E rest L2 b 2 q 1,2,3... To find the normalization constants, we integrate over all space in the box: b L 2 2 2 1 A Sin k x x dx D Sin 2 k y y dy . 0 0 Although only the product of the constants A and D has physical meaning, it is often common to normalize each individual function X and Y independently. In this manner, you can think of the problem as simply being the product two 1-D infinite well problems from Module 10. Doing this gives us, the following: A 1 L Sin 2 2 L 2 b k x x dx 0 D 1 b Sin 0 2 k y y dy Thus, the energy eigenfunctions for this problem are 2 nπx qπy Sin Sin 0 x L and 0 y b n,q b L L b 0 otherwise with associated energy eigenvalues of h c2 n 2 q2 n 1,2,3,.. E n,q where 8 E rest L2 b 2 q 1,2,3... Since this is a 2-dimensional problem, there must be at least two quantum numbers which determine an energy state’s wave function and energy eigenvalue (one for each spatial coordinate). This is an important concept and has led us to find additional variables which effect a systems energy but which are not spatial in nature like spin!! It is also possible for two different energy wave functions (energy levels) to have the same energy eigenvalues. This is called degeneracy. One way in which this might happen is if the system possesses some symmetry. For instance if the box is a square with sides of length L then E1,3 and E3,1 state would have the same energy as there is no difference to the system is you exchange the x and y coordinates (rotate 90 degrees). This is called symmetrical degeneracy. Many of the excited energy levels in the square system would be degenerate. In other cases, degeneracy might occur just due to pure luck. This is called accidental degeneracy. For instance, our original system is degenerate for 2 2 2 2 any group of (n,q) pairs where n b q L constant . II. Hydrogen Atom A. Coulomb Potential To solve for any quantum problem, we first must find the potential energy function for the problem. In the case of the hydrogen atom or any single electron ion, the potential energy function is the Coulomb potential which can be written in Cartesian coordinates as U ke 2 Z x 2 y2 z2 where ke 1.44 eV nm . Although it is technically correct to write the Coulomb potential in Cartesian coordinates, it is a poor choice for solving the Schrodinger equation since: 1) We can’t use separation of variable when the potential has terms containing more than one variable. 2) The system and hence its boundary conditions has spherical symmetry. 2 Hence, we need to write the Coulomb potential using spherical coordinates as U ke 2 Z r where Z = 1 for Hydrogen. B. Reduced Mass - In order to reduce the hydrogen atom from a two body problem to a one body problem, we must replace the mass of the electron with the reduced mass (as previously discussed in the module on the Bohr atom) μ mM mM where M is the mass of the proton and m is the mass of the electron. C. Laplacian in Spherical Coordinates The Laplacian is much more complicated in coordinate systems other than Cartesian, but it is necessary to use these coordinate systems for reasons previously discussed. The Laplacian in various coordinate systems can be found in math handbooks like CRC or Schaum’s. The Laplacian in spherical coordinate systems is given as 1 2 1 1 2 2 r Sinθ 2 2 θ r Sin θ r r r r 2Sinθ θ 2 D. Time Independent Schrodinger Equation We can now put everything together to write the T.I.S.E. for the hydrogen atom. 1 2 Ψ 1 Ψ 1 2Ψ 2μ ke 2 Ψ 2 E r Sinθ θ r 2Sin 2θ r r 2 r r r 2Sinθ θ We now write the wavefunction as a product of three single variable functions as Ψ R r Θθ Φ . Substituting into the T.I.S.E., we have ΘΦ 2 R RΦ Θ RΘ 2 Φ 2μ ke 2 RΘΘ 2 E r Sinθ θ r 2Sin 2 θ r r 2 r r r 2Sinθ θ 1 2 R 1 1 2 2μ ke 2 2 E r Sinθ θ r 2Sin 2 θ r r 2 R r r r 2Sinθ θ We can isolate the phi dependence term by multiplying both sides by r Sin θ and rearranging the terms. 2 2 Sin 2θ 2 R Sin θ 1 2 2μ r 2 Sin 2θ ke 2 r Sinθ E R r r θ θ r 2 Sin 2 θ 2 R Sin θ 2μ r 2 Sin 2 θ ke 2 1 2 r Sin θ E R r r θ θ r 2 Because the left-hand side of the equation is independent of phi and must be equal to the right-hand side for any value of phi, the equation must be equal to a 2 constant which we will denote by m l . Thus, we now have an ordinary differential equation for phi and a partial differential equation involving theta and r. d 2Φ 2 ml d “Azimulthal Equation” Sin 2 θ 2 R Sin θ 2μ r 2 Sin 2 θ ke 2 2 E m l r Sin θ 2 R r r θ θ r “Theta-Radial Partial Differential Equation” We continue or process of isolating variables on the remaining partial differential 2 equation by dividing the equation by Sin θ to isolate theta. 1 2 R 1 2μ r 2 r Sin θ R r r Sin θ θ θ 2 1 2 R 2μ r 2 r R r r 2 2 ke 2 m l E 2 r Sin θ 2 ke 2 m l 1 E Sin θ 2 r Sin θ Sin θ θ θ Since the left-hand side must equal the right-hand side of the equation regardless of the values of r and theta, the equation must equal a constant which we will denote by l (l 1) . This gives us our two final ordinary differential equations: 1 d 2 dR 2μ r 2 r R dr dr 2 ke 2 E l l 1 r “Radial Differential Equation” ml 1 d d Sin θ l l 1 dθ Sin 2θ Sin θ dθ 2 “Theta Equation” E. Radial Equation 1. Orbital Angular Momentum and Effective Potential In order to obtain some physical understanding about the radial equation and to put it into a form where we can find its solution in math handbooks, we are going to do some algebraic rearrangement to the equation. d2R dR 2μ r 2 r 2r dr dr 2 2 2 ke 2 2 l l 1 E R 0 r 2μ r 2 d 2 R 2 dR 2μ ke 2 2 l l 1 R 0 2 E r dr 2 r dr 2μ r 2 We also note that the potential of our electron for the Coulomb potential is Ur k e2 r Thus, we can rewrite the radial equation in a general form for any problem in which the potential depends only on the radial variable and not just the Coulomb potential as d 2 R 2 dR 2μ 2 l l 1 R 0. 2 E U r dr 2 r dr 2μ r 2 By dimensional analysis and our past knowledge of the Schrodinger equation in one dimension, we know that the term in the bracket must represent a kinetic energy due to motion in the radial direction. K radial K rotation E Ur Thus, the last term in the bracket must be kinetic energy due to rotation (i.e. due to changes in theta and phi). Thus, we have the kinetic energy of rotation as 2 l l 1 K rotation 2μ r 2 We know from PHYS122 that the kinetic energy of rotation is given as L2 K rotation 2I where L is the magnitude of the object’s orbital angular momentum and I is the moment of inertia. For a point particle of mass , the moment of inertia is I μ r2 . Thus, our classical physics work gives us the rotational kinetic energy as L2 . K rotation 2μ r 2 Comparing our classical physics results with the result from the radial part of the Schrodinger equation, we now observe relationship between the constant we used to separate the theta and radial equations and the maginitude of the orbital angula momentum. L l l 1 “Orbital Angular Momentum Eigenvalues” If we again look at the radial equation at the top of the page, we see that the rotational kinetic energy reduces the available energy for radial kinetic energy just as if it was a potential energy. This is a common occurrence in a wide range of physics problems involving radial forces including planetary motion. In these cases, physicists often chose to talk about an “effective potential” which is the sum of the actual potential due to forces and the kinetic energy due to rotation which is called the centrifugal barrier potential or just the ”centrifugal barrier.” L2 2 l l 1 U centifugal 2I 2μ r 2 2 l l 1 U effective Ur 2μ r 2 Thus, we can rewrite the radial equation in a more general form that is seen in many different fields of physics as d 2 R 2 dR 2μ 2 E U effective R 0 . dr 2 r dr As an example of how our work transcends various problems, you should note that according to our formula the angular momentum of the moon can be thought of as a centrifugal barrier (repulsive potential) which prevents the moon from hitting the Earth. A similar argument explains why it is harder for a nucleus to undergo alpha decay between two nuclear states of different angular momentum. In alpha decay, the alpha particle must tunnel through the nuclear potential. In order to conserve angular momentum in decay, the difference in the angular momentum of two nuclear states must be offset by the angular momentum of the emitted alpha particle. Because the effective potential is raised when the alpha particle has angular momentum, the probability for tunneling is greatly reduced. 2. Solution for the Coulomb Potential We now wish to change our approach from looking at the general nature of the radial equation to solving it for our specific case of the Coulomb potential. d 2 R 2 dR 2μ k e 2 2 l l 1 R 0 2 E r dr 2 r dr 2μ r 2 a. Asymptotic Solution for Large Radial Distances For large r, the differential equation is approximately, d 2 R 2μ E R 0 dr 2 2 d 2R 2μ E 2 R 2 dr 2μ E d2R 2 C R where C since E 0 dr 2 This is an equation which we have solved several times previously. Since the wave function must stay finite at large values of r, the only acceptable solution is the exponential decay solution. R r e b. 2μ E r for large r Ground State Wavefunction Let us see if our asymptotic wavefunction is a possible solution to the hydrogen atom for all radial locations in some cases. We do this by using a trial wavefunction based upon our asymptotic results and plugging it into the full radial equation. Let us define the radial wave function in terms of a constant delta as R Ae r δ and find the necessary radial derivatives as dR 1 R dr d2R 1 R. dr 2 2 Plugging our results into the radial Schrodinger equation, we have r2 2r 2μ 2 2 2 2 l l 1 2 E r k e r 0 2μ 2 2 2 μ E 1 2 2μ k e r l l 1 0 2 2 r 2 In order for the polynomial to be true for all values of r, the coefficients for each term must be zero!! This gives us three conditions that must be met for the trial function to be a solution for the hydrogen atom. The 2nd order coefficient gives us 2μ E 1 2 2 0 2 E 2μ 2 The linear coefficient gives us 2μ k e 2 2 0 2 μ k e2 . Solving for delta, we see that the constant is actually the first Bohr radius. 2 a o 0.529 x10 10 m 2 μke By combining results for the linear and 2nd order coefficients, we can determine the energy eigenvalue associated with our trial wavefunction. 1 E 2μ 2 μ k e2 μ 2 k 2 e4 13.6 eV 2 2 Thus, we see that our trial wavefunction will be a solution of the hydrogen atom only if the energy is the ground state energy of the hydrogen atom. We know look at the constant term. l l 1 0 This term will only be zero if l 0 or l 1 . We will see later that the theta equation restricts l to non-negative values so we have found that l 0 . This implies that the electron in the ground state has zero angular momentum since 00 1 0 in contradiction to the Bohr atom. The ground state electron actually has some intrinsic angular momentum. The Schrodinger equation doesn’t correctly predict this. Instead, we would need to use the Dirac equation from relativistic quantum mechanics. R r A e r ao with l 0 and E 13.6 eV “Radial Wavefunction for Ground State of Hydrogen” c. Asymptotic Solution for Small Radial Distances Let us start by rewriting our radial differential equation as d 2R dR 2μ 2 2 l l 1 2 R 0 . r 2r 2 E r ke r dr 2μ dr 2 2 For small values of r, we can simplify the differential equation as d2R dR r 2r l l 1R 0 . 2 dr dr 2 One of the important things that we see when we write the equation in this form is that for small radius values the radial function is primarily effected by angular quantum number l and not the wavefunctions energy eigenvalue E!! Because the electron must have a zero probability of being found at r = 0, it is reasonable to choose a power solution for R as our test wavefunction of the form R r A r c where c is a constant. Thus, we can now compute the required derivatives as dR R c A r c-1 c dr r d 2R R c- 2 c c 1 A r c c 1 dr 2 r2 Plugging our results into the Schrodinger equation, we find that r 2 c c 1 R R 2 r c l l 1R 0 2 r r c c 1 2 c l l 1R 0 c 2 c 2 c l l 1 R 0 c 2 c l l 1 R 0 cc 1 l l 1 R 0 In order for this to be true at all values of small r, the coefficient must be zero. Thus, we have shown that c l . This tells us that the radial wavefunction near the nucleus depends on the electron’s angular momentum!! R r A r l It is important that you realize how much information can be obtained without actually solving the full radial differential equation. This approach to problem solving is an important tool in your future work!! d. Solving the Radial Schrodinger Equation in General Eisberg and Resnick shows in Appendix N how to solve the general radial equation. I will leave the details for your reading, but I do want to point out some important facts that emerge from the process. i) Each solution has a polynomial whose form depends on both the orbital angular quantum number and a second constant n which is called the “principle quantum number.” ii) The constant n must be a non-zero positive integer and must be larger than the orbital angular quantum number in order for the radial equation to have a finite solution. iii) The energy eigenvalue associated with each radial wavefunction depends only on principle quantum number and has the value E n 13.6 eV . n2 The general form of the radial wavefunctions found by actually solving the radial Schrodinger equation can written in terms of the product of our asymptotic solutions and an unknown function. The general form of the radial wavefunctions in terms of the ratio of the radial distance to the Bohr radius is given by l r r n a R r H (r/a o ) e o ao l n where H(r/ao) is a polynomial whose form depends on two constants n and l . A few of the radial wavefunctions are provided below. For more functions look at Table 8.4 on page 280 of Serway/Moses/Moyer 3rd Edition as n R nl r l 3 2 1 0 r 1 ao 2e ao 3 2 r 1 r 2ao 2 0 2 - e 2ao a o 2 1 1 2ao 3 2 r r 2ao e 3ao F. Azimuthal (Phi) Equation We now return to the azimuthal equation to obtain its solution. d 2Φ 2 ml d We have seen this equation previously. In order for this equation to stay finite as we vary phi from minus infinity to positive infinity, it is necessary that m l 0 2 (i.e. m l must be real). Thus, the solution is of the form of a complex exponential A e i ml . Since we also require that the wave function be single valued at any point in space, the argument for the complex exponent must repeat anytime we increase or decrease phi by an integer multiple of two pi. Evaluating the function at phi equal zero gives us 0 A e i 0 A 2 2 π A ei ml 2 A e i ml 2 1 m l 2 π (integer) 2π m l ...... 2, 1, 0,1, 2,..... We can determine the normalization constant by solving the equation 2π Φ 2 d 1 0 2π A 0 2 e -ml e ml d 1 2π A 2 d 1 0 A 2 1 2 A 1 2π 1 i ml e 2π “Normalized Phi Wavefunction” G. Theta Equation We know return to the theta equation which can be written as 2 1 d d ml Sin θ l l 1 0 . Sin θ dθ dθ Sin 2 θ The solutions to this equation will depend on two quantum numbers, l and m l . To find the solution to this equation, we define a new variable, x, as x Cos θ . We also find that dx Sin θ . Using these results and the chain rule of dθ Calculus, we can rewrite the theta equation as 2 1 dx d dx d ml Sin θ l l 1 0 Sin θ dθ dx dθ dx Sin 2 θ 2 - Sin θ d ml 2 d - Sin θ l l 1 0 . Sin θ dx dx Sin 2θ 2 d d ml 2 - Sin θ l l 1 dx dx Sin 2θ 0 We now use the fact that Sin θ 1 Cos θ 1 x to change the equation to 2 2 2 2 d 2 d ml x 1 l l 1 dx dx 1-x2 2 d ml 2 d 1 - x l l 1 dx dx 1- x2 0 2 d 2 d ml 1- x dx 2 2 x dx l l 1 1 - x 2 2 0 0 This equation is known as Legendre’s Associated Differential Equation (see Eq. 28.49 page 165 of the Schaum’s Math handbook ). The solutions of this equation are well know special functions of engineering and physics (See Mathematical Methods for Physicists by Arfkin or similar textbooks). Eisberg and Resnick also go through the general solution technique in Appendix N. The general results for the solution of Legendre’s Associated Differential Equation are the following: i) ii) The equation only has finite solutions if l is a non-negative integer such that l m l . The wavefunctions are normalized Associated Legendre ml Polynomials, Pl x , which are defined by the quantum number l and m l . 2l 1 l m l ! ml P x 2 l m l ! l θ lml H. Spherical Harmonics - Yl θ, m We can now combine our results for theta and phi to produce another famous set of special functions called spherical harmonics. Ylml θ, lml ml 2l 1 l m l ! ml Pl Cosθ eiml 4π l m l ! Spherical harmonics are tabulated in Table 8.3 on page 269 of Serway/Moses/Moyer 3rd Edition. You may also obtain them using relationship 28.50 for the Associated Legendre Polynomial on page 165 in the Schaum’s Math Handbook: 1 x 2 ml l P ml 2 2l l ! l d ml l 2 x 1 dx ml l 0 Example: Find Y1 Solution: 0 First find P1 1 x 2 0 0 1 P d 2 x 1 2 (1!) dx P10 1 d 2 2x x 1 2 dx 2 P10 x Cosθ Now plug your result into the relationship for the spherical harmonics. I. Y10 θ, 2(1) 1 1 0! 0 i (0) P e 4π 1 0! 1 Y10 θ, 3 1 3 Cosθ Cosθ 4π 2 π Z-Component of the Angular Momentum We know from classical physics that orbital angular momentum is defined by the equation L r p . From the definition, the x-component of the angular momentum is found to be Lz x p y - y px . If we now use the operators for position and linear momentum that we discovered when developing the Schrodinger equation several modules ago, we have L̂ z i x i y i x y . y x y x However, we can use the chain rule of Calculus and the relationships between Cartesian coordinates and spherical coordinates (see Eisberg and Resnick Appendix M) to rewrite the z-component of the angular momentum’s operator as L̂ z i Applying the operator to our spherical harmonics, we find that L̂ z Ylml θ, 2l 1 l m l ! ml iml Pl Cosθ i e 4π l m l ! L̂ z Ylml θ, i i ml 2l 1 l ml ! ml P Cosθ eiml 4π l ml ! l L̂ z Ylml θ, m l Ylml θ, . Thus, we see that the spherical harmonics are eigenfunctions for the z-component of the angular momentum operator with eigenvalues of m l . Another interesting fact that comes out of this study is that there is an uncertainty relationship between the various components of the orbital angular momentum. You can know with perfect precision the magnitude of the angular momentum and the projection of the angular momentum upon one axis, but not its other two components!! We summarize these results as follows: i) The eigenfunctions for both angular momentum ands the zcomponent of angular momentum are the spherical harmonics. ii) The z-component of the orbital angular momentum is determined by the magnetic quantum number, m l , by the relationship Lz ml iii) The location of the angular momentum vector is quantized in space due to the relationships that we have developed. It traces out a cone in the x-y plane at an angle with respect to the z-axis given by Cosθ Lz ml L l l 1 z ml l l 1 y x J. Energy Degeneracy We found in our solution of the hydrogen atom that the energy depended only upon the principle quantum number, n, and not upon the orbital angular momentum quantum number, l , or the magnetic quantum number, m l . Thus, the system has two types of degeneracy!! The first degeneracy which is related to the energy being independent of the magnetic quantum number is due to the azimuthal symmetry of the systems. Since we could have chosen any direction in space for the z-axis, our energy results are independent of any phi dependence. This degeneracy is removed if the symmetry is broken by a non-spherical potential (such as a magnetic field) being applied to the system. This removal of degeneracy can be seen by additional lines in the emission spectra which originate when an external magnetic field is applied (Zeeman Effect). In fact even when no external field is applied, very fine line splitting is observed in precise experiments because the hydrogen atom does have a small internal magnetic field being applied to the electron due to coupling between the electron’s intrinsic angular momentum (spin) and its orbital angular momentum (fine structure) as well as magnetic effect due to the spin of the proton (hyper-fine structure). We will take up the effects of magnetism and the evidence for intrinsic angular momentum in the next module. The second type of degeneracy has to due with the orbital angular momentum quantum number. This type of degeneracy is due to the 1/r nature of the Coulomb potential. For other quantum problems including multi-electron atoms where electrons make the potential non-central or for systems involving potentials which are radial but not 1/r as in the case of Van der Waals forces, this energy does depend on the angular quantum number. III. Calculating Probabilities in 3-D The basic method for calculating the probability that a particle is in some 3-D region of space is the same as in Module 3 except that the integral is now over a volume instead of a line P 2 dv . space To be able to evaluate the integral in 3-D, we must first write an infinitesimal volume of space in spherical coordinates. To determine the formula for a differential volume element, we incrementally adjusting each coordinate independently to obtain a differential length, dS. These three differential lengths form the sides of cube whose differential volume is obtained by dv dS1 dS 2 dS3 A. Cartesian Coordinates For Cartesian coordinates, we find that the three differential lengths are dS1 dx , dS 2 dy , and dS3 dz with the differential volume given by dv dS1 dS 2 dS3 dx dy dz . Thus, the probability for finding the particle in a region of space where x 1 x x 2 , y 1 y y 2 , and z 1 z z 2 is given by x 2 y2 z2 P dz dy dx . * x1 y1 z1 If the wavefunction has been written as a product of single variable wavefunctions as in the case of the particle in a box (i.e. Ψ XYZ ), the probability equation can be written in the form y2 x2 z2 P X dx Y dy Z dz 2 x1 B. 2 y1 2 z1 Spherical Coordinates In spherical coordinates, we find that the three differential lengths are dS1 dr , dS 2 r dθ , and dS3 r Sinθ d . with the differential volume given by dv dS1 dS2 dS3 r 2 Sinθ dr dθ d . You should draw a picture of a spherical coordinate system and work through this Calculus until you can quickly derive this result whenever you need it. Thus, the probability for finding the particle in a region of space where r1 r r2 , θ 1 θ θ 2 , and 1 2 is given by θ2 2 2 r2 P r R 2 l n dr 2 Ylml Sinθ dθ d θ1 1 r1 Another common thing to do with this result is to use the fact that dCosθ Sinθ and rewrite the equation as P r R r1 Cosθ1 2 2 r2 2 l n dr Ylml 2 dCosθ d . Cos θ 2 1 In many problems, we are interested in the probability that the particle is any where in some radial region of space. In this case our integrations across theta and phi are all possible values (i.e. 0 2 and 0 θ ). In this case, our previous work gives us 1 2 2 r2 P r R 2 l n dr Ylml 2 dCosθ d . -1 0 r1 Since the spherical harmonics are normalized wavefunctions, we have 1 2 Ylml 2 dCosθ d 1 . -1 0 Therefore, the radial probability for finding a particle at a distance between r1 and r2 is 2 r2 P r 2 R ln dr . r1 You should recall that in Module 3 we defined the probability density function by the equation x2 P Ρ dx . x1 By comparing our radial probability equation to our work on Module 3, we define the radial probability density function as 2 P r 2 R ln . It is important to realize that the radial probability density can go to zero when r goes to zero even though the radial wavefunction doesn’t (consider our ground state work earlier). The probability of finding the electron depends not only on its probability of occupying a specific location (determined by the square of the radial wavefunction) but also upon the number of available locations (determined by the size of the sphere through r2). To find the values for r where the probability density is a maximum or minimum (i.e. most probable and least probable radial positions for finding the particle), you can either plot the probability density function using Excel, Mathcad, Maple, etc. or you take advantage of your knowledge of Calculus I and set the derivative of the radial probability density with respect to r equal to zero and solve for the roots. d P 0 dr “Condition for Maximum or Minimum Radial Probability” You know that the probability density is zero r = 0 and r = . For other possible results, you can plug your root back in to the probability density equation to see if it is zero or a maximum.