Legendre polynomials, Jackson's book problem, potential

In summary, the conversation is about a problem involving two concentric spheres divided into hemispheres with different potentials. The goal is to determine the potential in the region between the two spheres using a series in Legendre polynomials. A method involving boundary conditions and integrals is discussed, with the suggestion to only calculate the coefficients up to l = 4.
  • #1
fluidistic
Gold Member
3,944
263

Homework Statement


I'm stuck in evaluating an integral in a problem. The problem can be found in Jackson's book page 135 problem 3.1 in the third edition. As I'm not sure I didn't make a mistake either, I'm asking help here.
Two concentric spheres have radii a,b (b>a) and each is divided into two hemispheres by the same horizontal plane. The upper hemisphere of the inner sphere and the lower hemisphere of the outer sphere are maintened at potential V. The other hemispheres are at zero potential.
Determine the potential in the region [itex]a \leq r \leq b[/itex] as a series in Legendre polynomials. Include terms at least up to [itex]l=4[/itex]. Check your solution against known results in the limiting case [itex]b \to \infty[/itex] and [itex]a\to \infty[/itex].


Homework Equations


Below.


The Attempt at a Solution


I know I could solve this problem via the Green function for 2 spheres, I will try once I solve it via a less tricky way to see if I can get the same result.
Attempt: I realize it's an azimuthal symmetric problem. Therefore the solution is of the form [itex]\Phi (r, \theta )= \sum _{l=0}^{\infty } [U_l^+ r^l + U_l^- r ^{-(l+1)}] P_l (\cos \theta )[/itex]. My goal is to determine these 2 "U" coefficients. In order to do this, I check out what the boundary conditions have to tell me.
1)[itex]\Phi (a, \theta )= \sum _ {l=0}^{\infty } [U_l^+ a^l +U_l ^- a^{-(l+1)}]P_l (\cos \theta )[/itex]. I know that this expression can be worth either 0 or V depending upon the value of theta, I'll use this a step later. Now I use a trick, I multiply by [itex]P_{l'} (\cos \theta ' )[/itex] and I integrate as to use the orthogonality of the Legendre polynomials.
[tex]\int _{-1}^1 \Phi (a,\theta )P_{l'} (\cos \theta ') d \cos \theta =\left ( \frac{2}{2l+1} \right ) [U_l^+ a^l +U_l ^- a^{-(l+1)}] [/tex]. But from cos theta worth -1 to 0, the integral vanishes because Phi is 0 there. Hence I'm left with [itex]\int _{-1}^0 \Phi (a,\theta )P_{l'} (\cos \theta ') d \cos \theta =V \int _{-1}^0 P_{l'} (\cos \theta ' ) d \cos \theta =-V \int _{\pi} ^{\pi/2} P_{l'} (\cos \theta ' ) \sin \theta d\theta [/itex]. I must solve this integral. I'm not sure it is "right", I probably have to consider [itex]-V \int _{\pi} ^{\pi/2} P_{l} (\cos \theta ) \sin \theta d\theta [/itex].
I've been searching lots of relations/properties of the Legendre polynomials to see if one could help me. For example I can rewrite [itex]P_l (x)[/itex] as [itex]\frac{1}{(2l+1)} \left [ \frac{d}{dx} P_{l+1 }(x)-\frac{d}{dx} P_{l-1} (x) \right ][/itex]. But it doesn't seem to help me in any way.
I'm out of ideas. If I had this integral calculated, I'd have one expression of the [itex]U_l ^+[/itex] and [itex]U_l ^-[/itex] in terms of "a".
By looking at the other boundary conditions, 2)[itex]\Phi (b, \theta )[/itex], I guess I could do a similar job and find another expression of these U's in terms of b this time. This would give me 2 equations with 2 unknown, so I could solve it... and thus the problem.
Any help is greatly appreciated.
 
Physics news on Phys.org
  • #2
I'm not sure if it is possible to find a general expression for the coefficients, but in this problem you're only asked to find the coefficients up to [itex]l = 4[/itex]. So you could simply calculate the integral by hand for the first four Legendre polynomials, no?
 
  • #3
Steely Dan said:
I'm not sure if it is possible to find a general expression for the coefficients, but in this problem you're only asked to find the coefficients up to [itex]l = 4[/itex]. So you could simply calculate the integral by hand for the first four Legendre polynomials, no?

Hey, thanks for the help!
Hmm I don't really understand even what the integral means. [itex]l'[/itex] is arbitrary? So that the integral could be for example [itex]\int _{\pi } ^{\pi/2 } P_{10239843} ( \cos \theta ) \sin \theta d \theta[/itex] or [itex]\int _{\pi } ^{\pi/2 } P_3 ( \cos \theta ) \sin \theta d \theta[/itex]. Am I mistaken?
 
  • #4
fluidistic said:
Hey, thanks for the help!
Hmm I don't really understand even what the integral means. [itex]l'[/itex] is arbitrary? So that the integral could be for example [itex]\int _{\pi } ^{\pi/2 } P_{10239843} ( \cos \theta ) \sin \theta d \theta[/itex] or [itex]\int _{\pi } ^{\pi/2 } P_3 ( \cos \theta ) \sin \theta d \theta[/itex]. Am I mistaken?

The point of the "trick" you've used is that it allows you to pick out any value of [itex]l'[/itex] that you want, and then the right hand side will be the coefficients for that choice. So yes, the idea is that you can use the equation you've derived to figure out any of the expansion coefficients that you wish -- but the catch, as you've found, is that you have to do them one at a time since there's no general answer for what the integral of an arbitrary Legendre polynomial is. That's why Jackson only asked you to compute the first four coefficients, instead of working out the entire series expansion.
 
  • #5
Steely Dan said:
The point of the "trick" you've used is that it allows you to pick out any value of [itex]l'[/itex] that you want, and then the right hand side will be the coefficients for that choice. So yes, the idea is that you can use the equation you've derived to figure out any of the expansion coefficients that you wish -- but the catch, as you've found, is that you have to do them one at a time since there's no general answer for what the integral of an arbitrary Legendre polynomial is. That's why Jackson only asked you to compute the first four coefficients, instead of working out the entire series expansion.
Ah I understand now, thanks a lot.
I'm having a problem for the coefficients [itex]U_2[/itex]'s.
From the first boundary conditions I've reached that [itex]V \int _{1}^0 P_l (x)dx=\frac{2}{2l+1} [U_l ^+a^l+U_l^-a^{-(l+1)}][/itex]. (Note that I got wrong the limits of this integral in my first post)
While for the second boundary conditions, [itex]V \int _{0}^{-1} P_l (x)dx=\frac{2}{2l+1} [U_l ^+b^l+U_l^-b^{-(l+1)}][/itex].
For [itex]l=0[/itex] this gives me [itex]U_0^+=\frac{V}{2}+\frac{b}{a-b}[/itex] while [itex]U_0^-=\frac{ab}{b-a}[/itex].
For l=1, [itex]U_1^+=\frac{3V}{4b}+\frac{V}{2b^3} \left ( 1+\frac{a}{b} \right ) \left ( \frac{3a^2b^3}{2b^3-2a^3} \right )[/itex] and [itex]U_1^-=-\frac{V}{2}\left ( 1+ \frac{a}{b} \right ) \left ( \frac{3a^2b^3}{2b^3-2a^3} \right )[/itex].
While for [itex]l=2[/itex], the system of equations to solve leads to an infinity of solutions.
I reach in this case [itex]U_2^+=U_2^- \frac{\left ( \frac{1}{a^3} -\frac{1}{b^3} \right ) }{(b^2-a )}[/itex]. Thus any value of [itex]U_2^+[/itex] determine a value for [itex]U_2^-[/itex]. I know this cannot be right...
 
  • #6
fluidistic said:
Ah I understand now, thanks a lot.
I'm having a problem for the coefficients [itex]U_2[/itex]'s.
From the first boundary conditions I've reached that [itex]V \int _{1}^0 P_l (x)dx=\frac{2}{2l+1} [U_l ^+a^l+U_l^-a^{-(l+1)}][/itex]. (Note that I got wrong the limits of this integral in my first post)
While for the second boundary conditions, [itex]V \int _{0}^{-1} P_l (x)dx=\frac{2}{2l+1} [U_l ^+b^l+U_l^-b^{-(l+1)}][/itex]

You've still got the signs wrong; the original integral is from -1 to 1, so in either case you simply cancel out the part of the integral where the potential is zero, but you don't reverse the position of the bounds.

As for the [itex]l = 2[/itex] case, it doesn't lead to infinite solutions, it actually leads to only one possible case, if I did the calculation correctly.
 
  • #7
Steely Dan said:
You've still got the signs wrong; the original integral is from -1 to 1, so in either case you simply cancel out the part of the integral where the potential is zero, but you don't reverse the position of the bounds.

As for the [itex]l = 2[/itex] case, it doesn't lead to infinite solutions, it actually leads to only one possible case, if I did the calculation correctly.

Hmm.
From theta equal 0 to pi/2, [itex]\Phi (a, \theta ) =V[/itex]. The integral goes from theta = 0 to pi. Or equivalently from [itex]\cos \theta[/itex] equal to 1 to -1... or I'm missing something? Why would it go from -1 to 1?
 
  • #8
fluidistic said:
Hmm.
From theta equal 0 to pi/2, [itex]\Phi (a, \theta ) =V[/itex]. The integral goes from theta = 0 to pi. Or equivalently from [itex]\cos \theta[/itex] equal to 1 to -1... or I'm missing something? Why would it go from -1 to 1?

Either way, this is not important, as long as you do the integral consistently. You must be making a calculational error if you think there are other solutions than the zero solution for [itex]l = 2[/itex].
 
  • #9
Steely Dan said:
Either way, this is not important, as long as you do the integral consistently. You must be making a calculational error if you think there are other solutions than the zero solution for [itex]l = 2[/itex].

I've redone entirely half of the problem (only one boundary condition) but this time with taking extreme care of the limits of integration; in Jackson's book I've found a formula for [itex]\int _0 ^1 P_l (x)dx[/itex]. If I'm not wrong, it is worth [itex](-1) \left ( \frac{-1}{2}\right ) ^{\frac{l-1}{2}} \frac{(l-2)!}{2 \left ( \frac{l+1}{2} \right ) }[/itex]. So that it seems I could get a general expression for the problem and from it I could get any term I would like.
I'll try to finish the problem when I have some time for it.
 

FAQ: Legendre polynomials, Jackson's book problem, potential

1. What are Legendre polynomials?

Legendre polynomials are a set of orthogonal polynomials that are commonly used in mathematics and physics. They are named after French mathematician Adrien-Marie Legendre and are defined by a recursive formula.

2. How are Legendre polynomials used in physics?

Legendre polynomials are used to represent solutions to Laplace's equation, which is a fundamental equation in physics that describes the potential in a region due to a given charge distribution. They are also used in the expansion of functions in spherical harmonics, which have many applications in electromagnetism and quantum mechanics.

3. What is Jackson's book problem?

Jackson's book problem is a mathematical problem named after physicist John David Jackson. It involves finding the maximum electrostatic potential on the surface of a charged conducting book, assuming the book is an infinite plane with uniform charge distribution. The solution to this problem involves the use of Legendre polynomials.

4. How are potential and potential energy related?

Potential energy is the energy associated with the position of an object in a conservative force field. The potential energy at a particular point is equal to the work done by the conservative force in moving the object from a reference point to that point. Potential, on the other hand, is the amount of potential energy per unit charge. Therefore, potential is the potential energy per unit charge at a particular point in an electric field.

5. What is the significance of potential in physics?

Potential is a crucial concept in physics as it helps us understand and analyze the behavior of electric and gravitational fields. It is also used to calculate the work done by conservative forces and the potential energy of a system. In electromagnetism, potential is related to the electric field, and in quantum mechanics, it is related to the wave function of a particle. Overall, potential plays a vital role in many areas of physics and has a wide range of applications.

Back
Top