===== Nov 06 (Fri) Legendre polynomials, Radial wave equation ===== ** Responsible party: liux0756, Dagny ** **To go back to the lecture note list, click [[lec_notes]]**\\ **previous lecture note: [[lec_notes_1104]]**\\ **next lecture note: [[lec_notes_1109]]**\\ **Main class wiki page: [[home]]** Please try to include the following * main points understood, and expand them - what is your understanding of what the points were. * expand these points by including many of the details the class discussed. * main points which are not clear. - describe what you have understood and what the remain questions surrounding the point(s). * Other classmates can step in and clarify the points, and expand them. * How the main points fit with the big picture of QM. Or what is not clear about how today's points fit in in a big picture. * wonderful tricks which were used in the lecture.\\ \\ Today we focus on the following two points. * rest of Legendre polynomial * Radial part of Schrodinger equation ==== Part I: Legendre polynomial ==== In the last class the power series method is used: U(\xi)=\sum_{n=0}^\infty a_n \xi^2 Substitute it into differential equation derives the recursive relation: a_k=\frac{(\nu+1)(\nu+2)...(\nu+k)(-\nu)(-\nu+1)...(-\nu+k-1)}{(k!)^2} a_0 a_0 is determined by normalization. Now consider the convergence requirement. Just as the simple harmonic oscillator, for normalizable solutions the power series must terminate. And if \nu is an integer, a_k will be zero if k is large enough. If \nu is negative, there exists a corresponding positive \nu that leads to the same recursive relation. So \nu can be limited to non-negative numbers: \nu=0,1,2,... \alpha=\nu(\nu+1)=0,2,6,12,... Generally P_\nu (z) is called Legendre function, \nu is any real number. If convergence is important, \nu=l is integer, we deal with the Legendre polynomial P_l (z) ==== Part II: Radial part of Schrodinger equation ==== The equation is written as: \frac{d}{dr} (r^2 \frac{dR}{dr})- \frac{2mr^2}{\hbar^2}[V(r)-E]R=l(l+1)R l(l+1) is related with the square of angular momentum: L^2=l(l+1)\hbar^2 Define R(r)=\frac{u(r)}{r}, then the differential equation is transformed from R(r) to u(r): -\frac{\hbar^2}{2m} \frac{d^2u}{dr^2}+[V+\frac{\hbar^2}{2m} \frac{l(l+1)}{r^2}]u=Eu This equation is similar to 1D Schrodinger equations discussed in Chapter 2. The equation above cannot be solved further before one knows the potential distribution in the system. === Example === Since we cannot go much further without specifying potential energy, now consider a 3D infinite square well. When r>a, the potential goes to infinity. When r, the potential is 0. The radial differential equation is: \frac{d^2u}{dr^2}-[\frac{l(l+1)}{r^2}-k^2]u=0 where k^2=\frac{2mE}{\hbar^2} The solution is Bessel function j_l(kr) Now take a look at the 0th order solution. When l=0, the differential equation is \frac{d^2u}{dr^2}=-k^2u The solution is: u(r)=Ae^{ikr}+Be^{-ikr}=A'sinkr+B'coskr The boundary condition is u(a)=0 At r=0, if u is not 0, then R=\frac{u}{r} will go to infinity. So u(0)=0 should be satisfied. (However, this is not quite correct because even if the value of wave function is infinite, the wave function may still be able to be normalized: \int_0^a \frac{1}{r^2}dV = \int_0^a \frac{1}{r^2} r^2 sin\theta dr d\theta d\phi=\int_0^a sin\theta dr d\theta d\phi is normalizable - finite) So the B'coskr term should be 0. (In fact this term is allowed) ------------------------------------------- **To go back to the lecture note list, click [[lec_notes]]**\\ **previous lecture note: [[lec_notes_1104]]**\\ **next lecture note: [[lec_notes_1109]]**\\