Potential from a Quadrupole using Legendre polynml's

[Vapor88]

1. Problem Statement:
There are charges of q placed at distance +a and -a from the origin on the z-axis. There is a charge at the origin of -2q. Express the potential of this point-like linear quadrupole in Legendre polynomials

The distance between origin and point is r, the distance between (0,0,a) and r is l and the distance between (0,0,-a) is l'

The Attempt at a Solution

The potential can be expressed as
(q/(4\pi r)*(1/l + 1/l' + 2/r)

Using law of cosines

l = \sqrt{a^2 + r^2 - 2ar cos(\theta)}
l' = \sqrt{a^2 + r^2 + 2ar cos(\theta)}

Solving for 1/l and 1/l' gives
1/l = \frac{1}{r} (1 + a^2/r^2 - a/r cos (\theta))^{-1/2}
1/l' = \frac{1}{r} (1 + a^2/r^2 + a/r cos (\theta))^{-1/2}

that's supposed to be raised to the -0.5, but i don't know how to make it look right in latex

Now

V = ( \frac{q}{4 \pi \epsilon_{o}} \frac{1}{r} [(1+( \frac{a}{r})^2 - \frac{a}{r} cos (\theta))^{-1/2} + (1+ (\frac{a}{r})^2 + \frac{a}{r} cos (\theta))^{-1/2} - 2]

We had an example in class of a dipole that used the expansion

(1+x)^n = 1 + nx + ...

And only used the first two terms of it. Do I need to do the same thing to this, or can I go straight to the Legendre polynomials. I'm not sure where to go from here.

Also... If you take a^2/r^2 to go to zero, then the potential goes to zero as well. So I can't do that.

 

[Vapor88]

I just worked further into the problem, tried using that expansion the same way, don't know if I can really do that though.

But I took the a^2/r^2 + a/r cos (theta) = x

In the end, I did get something over r^3, which makes sense because it's a quadrupole, but I missed the point of the problem completely, using Legendre polynomials.

 

[vela]

The idea is to go straight to the Legendre polynomials using the well-known result

\frac{1}{\sqrt{r^2+a^2-2ar\cos\theta}}=\sum_{l=0}^\infty \frac{a^l}{r^{l+1}}P_l(\cos\theta)

when a<r.

 

[Vapor88]

Ok, so I got a bit further in the problem, he wants us to use the first three terms of binomial expansion, I only did the first two to get my results above, but they did not need Legendre polynomials. Though I do not know why my answer is incorrect. If I use the first three terms of the binomial expansion for both of the 1/(...)^1/2, then simplify, I end up with the following

\frac{q}{4 \pi \epsilon_{o} r} * [ - \frac{a^2}{r^2} + \frac{3a^4}{4r^4} + \frac{3a^2}{4r^2} cos^2( \theta)]

Now, my misunderstanding is in the Legendre polynomials. I know the definitions of Po P1 P2, etc. but I don't get how to apply them. I see that there is a cos^2 term, so I know that that's related to P2
P_{2} = \frac {(3x^{2}-1)}{2}
P_{1} = x
P_{0} = 1

/edit... P2 is supposed to be (3x^2-1)/2

And I imagine that in those terms, I have x = cos (theta)

I'm just very sketchy on the details of Legendre polynomials. Everyone in my class had my professor for quantum, I haven't taken quantum yet, so he really did not spend much time explaining it. As I understand it, the first two terms in my equation get Po applied to them, but I don't know how to do this.

 

[vela]

You made a mistake in your algebra somewhere. The cos^2 term shouldn't have a 4 in the denominator.

 

[Vapor88]

For
[1 + \frac{a^2}{r^2}- \frac{a}{r} cos( \theta)]^{-1/2}

I expand to get

[1 - \frac{1}{2} ( \frac{a^2}{r^2} - \frac{a}{r} cos( \theta)) + \frac{- \frac{1}{2} (- \frac{1}{2} - 1)}{2!} ( \frac{a^4}{r^4} - \frac{2a^3}{r^3} cos( \theta) + \frac{a^2}{r^2} cos^2( \theta))]

Which gets simplified to[1 - \frac{1}{2} ( \frac{a^2}{r^2} - \frac{a}{r} cos( \theta)) + \frac{3}{8} ( \frac{a^4}{r^4} - \frac{2a^3}{r^3} cos( \theta) + \frac{a^2}{r^2} cos^2( \theta))]

For the other term

[1 + \frac{a^2}{r^2}+ \frac{a}{r} cos( \theta)]^{-1/2}

When expanded, everything stays the same except all terms are positive

[1 - \frac{1}{2} ( \frac{a^2}{r^2} + \frac{a}{r} cos( \theta)) + \frac{3}{8} ( \frac{a^4}{r^4} + \frac{2a^3}{r^3} cos( \theta) + \frac{a^2}{r^2} cos^2( \theta))]

then, when you add em together you get

V = \frac{q}{4 \pi \epsilon_{o} r} * [ - \frac{a^2}{r^2} + \frac{3a^4}{4r^4} + \frac{3a^2}{4r^2} cos^2( \theta)]

...er there's still a 4 in the denominator... hrm

Oh! and the 1's disappeared because of the -2 seen in the first post

But yeah... I keep reworking it and get 3/8 + 3/8 = 3/4 for that last term

 

[vela]

Oh, you dropped the two way in the beginning. You wrote

1/l = [1 + \frac{a^2}{r^2} - \frac{a}{r} \cos(\theta)]^{-1/2}

but it should be

1/l = [1 + \frac{a^2}{r^2} - 2 \frac{a}{r} \cos(\theta)]^{-1/2}

 

[Vapor88]

Ok, so I believe this is correct, however, I'm still stumped as far as legendre polynomials go. I have no idea how to implement them.

<br /> V = \frac{q}{4 \pi \epsilon_{o} r} * [ - \frac{a^2}{r^2} + \frac{3a^4}{4r^4} + \frac{3a^2}{4r^2} cos^2( \theta)]<br />

Given this, I can see that it sort of looks like P2, but would that make it -1+1+3/2cos^2-1/2

This doesn't seem right at all.

 

[vela]

Ok, so I believe this is correct, however, I'm still stumped as far as legendre polynomials go. I have no idea how to implement them.

<br /> V = \frac{q}{4 \pi \epsilon_{o} r} * [ - \frac{a^2}{r^2} + \frac{3a^4}{4r^4} + \frac{3a^2}{4r^2} cos^2( \theta)]<br />

Given this, I can see that it sort of looks like P2, but would that make it -1+1+3/2cos^2-1/2

This doesn't seem right at all.

How is this different from what you wrote before? You still have the four in the denominator of the cos^2 term that doesn't belong there.

 

[Vapor88]

How is this different from what you wrote before? You still have the four in the denominator of the cos^2 term that doesn't belong there.

It's not, but I keep reworking it and get that. here is the work for how I got there

<br /> [1 - \frac{1}{2} ( \frac{a^2}{r^2} - \frac{a}{r} cos( \theta)) + \frac{3}{8} ( \frac{a^4}{r^4} - \frac{2a^3}{r^3} cos( \theta) + \frac{a^2}{r^2} cos^2( \theta))]<br />

Then you add that to

<br /> [1 - \frac{1}{2} ( \frac{a^2}{r^2} + \frac{a}{r} cos( \theta)) + \frac{3}{8} ( \frac{a^4}{r^4} + \frac{2a^3}{r^3} cos( \theta) + \frac{a^2}{r^2} cos^2( \theta))]<br />

You can see what cancels out and what gets carried through. If those two equations are correct, then the cos^2 term's coefficient is 3/4

Working back further...

For the expansion:
The 3/8 comes from n(n-1)/2! where n = -1/2

The x^2 in the expansion's third term becomes this
<br /> \frac{a^4}{r^4} + \frac {2a^3}{r^3} cos ( \theta) \frac{a^2}{r^2} cos^2( \theta)<br />

The 1's in front get canceled by the 2 in the first post. All of the cos(theta) terms cancel. The 3/8 gets distributed to the cos^2 term in both of them, then when I add them together, I get 3/4. What coefficient am I supposed to be getting, then? If I am making an algebra mistake, I am still not seeing it. Can you please point it out, I've just reworked it three times, and gotten the same thing.

 

[vela]

Read post 7 in this thread.

 

[Vapor88]

*facepalm*

thanks!

 

[Vapor88]

Alrighty...

After banging my head on the desk, I've got the following.

<br /> V = \frac{q}{4 \pi r \epsilon_{o}} [- \frac{a^2}{r^2} + \frac{3a^4}{4r^4} + \frac{3a^2}{2r^2} cos^2( \theta)]<br />

Now... to implement those Legendre polynomials.

 

[Vapor88]

Ok, here's what I know of Legendre plynml's from my notes

They are integrals such that

<br /> \int P_{l} (x) P_{m} (x)dx = 0<br />

When l =/= m

When l = m, that integral is equal to
<br /> \frac{2}{2l+1}<br />

For this problem, I want to take x to be cos(theta)? Does that turn the first two terms into terms that can be handled with Po and the third with P2?

 

[gabbagabbahey]

Alrighty...

After banging my head on the desk, I've got the following.

<br /> V = \frac{q}{4 \pi r \epsilon_{o}} [- \frac{a^2}{r^2} + \frac{3a^4}{4r^4} + \frac{3a^2}{2r^2} cos^2( \theta)]<br />

Now... to implement those Legendre polynomials.

What is P_2(\cos\theta)? How about P_0(\cos\theta)?

 

[Vapor88]

What is P_2(\cos\theta)? How about P_0(\cos\theta)?

<br /> P_{0} = 1<br />

<br /> P_{1} = cos ( \theta)<br />

<br /> P_{2} = \frac{3cos^2 ( \theta) - 1}{2}<br />

I'm fairly confident that this is correct.

Are these my P_{l} values? And would a^2/r^2 be Pm?

Does this mean that I'll have three integrals, all evaluated from -1 to 1, such that

<br /> <br /> \int P_{l} (cos \theta) P_{m} (cos \theta)d(cos \theta) = 0<br /> <br />

 

[gabbagabbahey]

<br /> P_{0} = 1<br />

<br /> P_{1} = cos ( \theta)<br />

<br /> P_{2} = \frac{3cos^2 ( \theta) - 1}{2}<br />

I'm fairly confident that this is correct.

Right. Now look at your expression for the potential..

V = \frac{q}{4 \pi r \epsilon_{o}} [- \frac{a^2}{r^2} + \frac{3a^4}{4r^4} + \frac{3a^2}{2r^2} cos^2( \theta)]

The first thing you should do is drop the \frac{a^4}{r^4}. In your Taylor expansion, you are assuming that a is very small compared to the distance from the origin at which you are measuring the potential. (If you don't make this assumption, you must keep all of the terms in the Taylor expansion) So, certainly \frac{a^4}{r^4} should be much much smaller than \frac{a^2}{r^2}, right?

After that, you should notice that

- \frac{a^2}{r^2} + \frac{3a^2}{2r^2} cos^2( \theta)=\frac{a^2}{r^2}\left(\frac{3\cos^2\theta-1}{2}\right)- \frac{a^2}{2r^2}=\frac{a^2}{r^2}P_2(\cos\theta)+\frac{a^2}{2r^2}P_0(\cos\theta)=\frac{a^2}{r^2}\left(P_2(\cos\theta)+\frac{1}{2}P_0(\cos\theta)\right)

Does this mean that I'll have three integrals, all evaluated from -1 to 1, such that

<br /> <br /> \int P_{l} (cos \theta) P_{m} (cos \theta)d(cos \theta) = 0<br /> <br />

What exactly are you planning on integrating?

Edit: You seem to have made an error in your Taylor expansion, to second order in \frac{a}{r}, I get

V\approx\frac{q}{4 \pi r \epsilon_0}\left[- \frac{a^2}{r^2} + \frac{3a^2}{r^2} cos^2( \theta)\right]=\frac{q}{2 \pi \epsilon_0}\left(\frac{a^2}{r^3}\right)P_2(\cos\theta)

 

[Vapor88]

What am I planning on integrating? lol... I'm not really sure, I'm very, very shaky on these legendre polynomials, they're completely new to me, but everyone else in my class is familiar with them.

Anyways, I see from there, when plugged into the potential equation, I get an r^-3, which is expected in quadrupoles. Is further work necessary?

 

[Vapor88]

I used the binomial expansion, not the Taylor expansion... They are different, right?

 

[gabbagabbahey]

I used the binomial expansion, not the Taylor expansion... They are different, right?

a^2+r^2\pm2ar\cos\theta is a trinomial, not a binomial. You should use a Taylor expansion instead.

 

[vela]

I used the binomial expansion, not the Taylor expansion... They are different, right?

The binomial expansion of (1+x)^1/2 is the Taylor expansion of (1+x)^1/2 about x=0.