Taking legendre polynomials outside the integral in a multipole expansion

mjordan2nd
Messages
173
Reaction score
1

Homework Statement



A chare +Q is distributed uniformly along the z axis from z=-a to z=+a. Find the multipole expansion.

Homework Equations



eimg182.gif


Here rho has been changed to lambda, which is just Q/2a and d^3r to dz.

The Attempt at a Solution



I have solved the problem correctly (confirmed since the answer is given in the book.) I could not figure out how to solve the problem. Just to see what happened, I took the legendre polynomial factor outside the integral, and I stumbled onto the correct answer. Was my attempt incorrect, and my correct solution a lucky coincidence, or was my method correct? If my method was correct, why are you allowed to take the legendre polynomial outside the integral? It seems to me that theta must certainly depend on z (I was convinced of this based on fig. 3.28 in Griffiths). It can also be seen from:

eimg178.gif


where r corresponds to z in my problem. Certainly in the figure above if you changed r theta would also change.
 

Attachments

  • eimg182.gif
    eimg182.gif
    1.4 KB · Views: 484
Physics news on Phys.org
You got lucky. That equation from ScienceWorld isn't even right. There should be a P_n(\cos\,\theta) outside the integral and a P_n(\cos\,\theta\,') inside the integral, where \theta\,' is one of the terms integrated over by d^3\mathbf r.
 
D H said:
You got lucky. That equation from ScienceWorld isn't even right. There should be a P_n(\cos\,\theta) outside the integral and a P_n(\cos\,\theta\,') inside the integral, where \theta\,' is one of the terms integrated over by d^3\mathbf r.

No, there's nothing wrong with his equation; it results from a straight forward expansion of |\textbf{R}-\textbf{r}|^{-1}=(R^2+r^2-2rR\cos\theta)^{-1/2} in the case R>r

mjordan2nd said:

Homework Statement


eimg182.gif


Here rho has been changed to lambda, which is just Q/2a and d^3r to dz.

It might help you to visualize what's going on here if you represent \rho in terms of delta functions and explicitly reduce the integral to a single one over z...

Certainly in the figure above if you changed r theta would also change.

Careful, since \theta is the angle between \textbf{r} and \textbf{R}, it will certainly change as \textbf{r} changes in direction, but changing the length of \textbf{r} will have no effect on the angle between it and \textbf{R}.

When you integrate over x and y, the direction of \textbf{r} will change, but the charge density is zero except for x=y=0, and hence contributes nothing to the integral. The only non-zero contribution comes when \textbf{r}=z\mathbf{\hat{z}}, and this vector points in the same direction throughout the integration over z, so \theta will be be a constant in your integral.
 
In general I disagree. In this case, it happens to work. That equation is from a wolfram article, which should have known better.

Another way to look at this is that the spherical harmonics form an orthogonal set of functions over the square integrable functions. One can represent a vector v in N-space in terms of any set of orthogonal vectors ui

\aligned<br /> \mathbf v &amp;= \sum_n a_n \mathbf u_n \\<br /> \mathbf a_n &amp;= \frac {\mathbf v \cdot \mathbf u_n}{||\mathbf u_n||}<br /> \endaligned

This concept carries over to function space spanned by a set of orthogonal polynomials, with the inner product defined in terms of an integral. The above finite sum becomes an infinite sum.

The expansion of 1/||\mathbf r - \mathbf r&#039;|| in terms of the spherical harmonics leads to the multipole expansion of the electrostatic potential. A couple of articles:
http://www.cartage.org.lb/en/themes/Sciences/Physics/Electromagnetism/Electrostatics/ElectrostaticPotential/Multipole/Multipole.htm
http://www.uni-regensburg.de/Fakultaeten/nat_Fak_IV/Theoretische_Chemie/Homeier/preprint/TC-QM-97-5/node2.html

In the special case of axial symmetry, the sectorial and tesseral harmonics terms vanish, leaving only the zonal harmonics. These can be rewritten in terms of the Legendre polynomials:

\aligned<br /> \Phi(r,\theta) &amp;= \sum_{n=0}^{\infty} \frac{a_n P_n(\cos\,\theta)}{r^{n+1}} \\<br /> a_n &amp;= \int d\mathbf r^{\prime}\, \rho(\mathbf r^{\prime}) r^{\prime n} P_n(\cos\,\theta&#039;)<br /> \endaligned

In the very special case of a charge restricted to the z axis, those P_n(\cos\,\theta&#039;) terms become 1 or -1; they nicely just drop out.

This is what I meant by he got lucky.
 
Last edited by a moderator:
I maintain my position that there is still nothing wrong with his equation

V(R)=\sum_{n=0}^{\infty} \frac{1}{R^{n+1}}\int r^n P_n(\cos\theta)\rho(\textbf{r})d^3\textbf{r}

even in the general case, so long as R&gt;r.

In the above equation, \theta is defined as the angle between the vectors \textbf{R} (the position of the field point relative to the origin) and \textbf{r} ((the position of each source point dq=\rho(\textbf{r})d^3\textbf{r} relative to the origin), not as the polar angle of \textbf{R} or \textbf{r} in general. This definition should be clear from the OP's diagram.
 
Thread 'Need help understanding this figure on energy levels'
This figure is from "Introduction to Quantum Mechanics" by Griffiths (3rd edition). It is available to download. It is from page 142. I am hoping the usual people on this site will give me a hand understanding what is going on in the figure. After the equation (4.50) it says "It is customary to introduce the principal quantum number, ##n##, which simply orders the allowed energies, starting with 1 for the ground state. (see the figure)" I still don't understand the figure :( Here is...
Thread 'Understanding how to "tack on" the time wiggle factor'
The last problem I posted on QM made it into advanced homework help, that is why I am putting it here. I am sorry for any hassle imposed on the moderators by myself. Part (a) is quite easy. We get $$\sigma_1 = 2\lambda, \mathbf{v}_1 = \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} \sigma_2 = \lambda, \mathbf{v}_2 = \begin{pmatrix} 1/\sqrt{2} \\ 1/\sqrt{2} \\ 0 \end{pmatrix} \sigma_3 = -\lambda, \mathbf{v}_3 = \begin{pmatrix} 1/\sqrt{2} \\ -1/\sqrt{2} \\ 0 \end{pmatrix} $$ There are two ways...
Back
Top