Complex vector potential of solid sphere and Heaviside function

[Shinobii]

Homework Statement

I have done this problem for the case of a spherical shell, however, I am not understanding how to go about this for a solid sphere.

Homework Equations

[tex] \vec{A} = \frac{1}{4 \pi} \int_{\phi' = 0}^{2 \pi} \int_{-1}^1 \int_0^R \rho_o \Theta(R-r) \sum_{l=0}^\infty \sum_{m=-l}^l \frac{4 \pi}{2l +1} \frac{r_<^l}{r_>^{l+1}}Y_{l,m}^*(\theta',\phi')Y_{l,m}(\theta,\phi) r'^2 d(\cos(\theta'))d\phi' [/tex]

The Attempt at a Solution

For the case of a shell, there is a delta function which makes life easy.

My question is, what do I do with this Heaviside function? Do I treat it differently for r < R and r > R? This is my first time encountering this function.

Any hints would be greatly appreciated.

 

[vela]

You're missing a dr in your integral. Given the limits on your integral, you always have ##r\le R##, so the step function is equal to 1.

 

[Shinobii]

Hi, sorry I meant to put in the dr'.

I guess when I am integrating from 0 to r > R, I would split the integral up and only the integral from 0 to R would contribute while the other integral (R to r) would be zero.

 

[vela]

What are the limits on your r' integral supposed to be? Are they 0 and R as you've written? Should the r inside the step function be r'?

 

[Shinobii]

Yes sorry, the r is r' inside the Heaviside function.

Do we split up the integral like the case of the electric potential? If so,

For r<R:

The integral, like solving for the electric potential, would have to broken up into two components, that is [itex] \int_{\infty}^r \rightarrow \int_{\infty}^R + \int_R^r. [/itex] Would the first integral = 0 by definition of the Heaviside function?

For r>R:

The integral is from [itex] \int_{\infty}^r [/itex] and therefore 0?

Again solving it for the shell case none of this was an issue...

 

[vela]

I'm still unclear on what you're calculating here. The righthand side doesn't appear to be a vector. Where are you getting ##r## and ##\infty## for your limits?

 

[Shinobii]

Yeah, I messed up real good. I forgot the cross product [itex](\omega \times r')[/itex] from the velocity.

Using the complex current density,

[tex] J_x + iJ_y = -i \sqrt{\frac{3}{2\pi}}\frac{q \omega r'}{ R^3} \Theta(R-r')Y_{1,1}. [/tex]

And using this to find the complex vector potential I get,

[tex] A_x + iA_y = -i\sqrt{\frac{3}{2\pi}}\frac{q \omega}{ R^3} \int_{\phi' = 0}^{2 \pi} \int_{-1}^1 \int_0^r \Theta(R-r')Y_{1,1}(\theta',\phi') \sum_{l=0}^\infty \sum_{m=-l}^l \frac{4 \pi}{2l +1} \frac{r_<^l}{r_>^{l+1}}Y_{l,m}^*(\theta',\phi')Y_{l,m}(\theta,\phi) r'^3 dr'd(\cos(\theta'))d\phi' [/tex]

From here I believe I can use orthogonality to kill off some terms setting l=m=1, and I get

[tex]A_x + iA_y = \sqrt{\frac{3}{2\pi}}\frac{q \omega}{ R^3} \int_0^r \Theta(R-r')Y_{1,1}(\theta',\phi') \frac{4 \pi}{2(1) +1} \frac{r_<^{(1)}}{r_>^{(1)+1}} r'^3 dr'[/tex]

I hope this is right thus far. . . And yes I am clearly not understanding what I am integrating over for the r' integration.

 

[vela]

You want to integrate over all space, so the limits should be 0 to infinity. Since the step function vanishes for r'>R, you can simply integrate from 0 to R.

 

[Shinobii]

You are correct, thanks! Sorry again for my sloppy mathematics, it's my pet peeve yet I am sloppy myself...

 

[Shinobii]

One more quick question on notation, would this be correct?

[tex] \int_0^{\infty} r' \, \Theta(R-r') dr' \rightarrow \int_0^{R} r' \, dr'. [/tex]

Would I simply change the limits of integration and remove the Heaviside function from the equation?

 

[vela]

Sorry for the late response. Yes, that's correct.