Integrating charge density

  • Thread starter T-chef
  • Start date
  • #1

Homework Statement

Find the electric charge centred in a sphere of radius a, centered at the origin where the electric potential is found to be (in spherical coordinates) [itex] V(r)=kr^-2 [/itex] where k is some constant.

The Attempt at a Solution

We have [itex] E=-\nabla V = -2kr^{-3} \hat{r} [/itex]
So applying Guass's law to the sphere of radius a, we get
[itex] \oint_s E \cdot da = \frac{Q}{\epsilon_0} [/itex]
And thus [itex] Q= \frac{-8\pi k \epsilon_0}{a} [/itex]

My problem is, surely the same result should be obtained by taking the triple integral of the charge density with respect to volume, but pursuing this path...
[tex] \rho=-\epsilon_0 \nabla^2 V(r) = -\epsilon_0 2kr^{-4} [/tex]
Attempting to integrate this in spherical coordinates results in,
[tex]Q= \iiint_V \rho dV = -8\pi \epsilon_0 k \int_0^a r^{-2} dr [/tex]
but due to the singularity this tends to infinity. Where did I take a wrong step?

Answers and Replies

  • #2
Science Advisor
Homework Helper
Hi T-chef! :smile:
We have [itex] E=-\nabla V = -2kr^{-3} \hat{r} [/itex]

Nooo … ∂/∂x 1/(x2 + y2 + z2) = … ? :wink:
  • #3
Hello Tiny Tim!
I was hoping by working in spherical coordinates I could make my life easier. My book assures me here grad is given by
[tex] \nabla t = \frac{\partial t}{\partial r}\hat{r} + \frac{1}{r}\frac{\partial t}{\partial \theta}\hat{\theta}+\frac{1}{rsin(\theta)}\frac{\partial t}{\partial\phi}\hat{\phi} [/tex]
So since V is independent of the angles, their partial derivatives vanish leaving me with just the first term. Have I misinterpreted working with grad in other coordinates systems?
  • #4
i think the potential for V is only the potential outside the sphere. inside the sphere the potential is affected by the charge distributions themselves
also the laplacian of the potential is NOT equivalent to taking the gradient twice, check the math and i think u will see something that looks like a dirac delta function
  • #5
I used the Laplacian as expressed in spherical coordinates, which i certainly agree, is not taking grad twice :smile:
That said, I suppose I could have used:
[tex] E=-2kr^{-3} \hat{r} [/tex]
[tex] \nabla \cdot E = \frac{\rho}{\epsilon_0} [/tex]
[tex] \rho = \epsilon_0 \nabla \cdot (-2kr^{-3} \hat{r}) [/tex]
Now in my (very limited!) knowledge of delta functions, don't I need to be applying divergence to a [itex] \frac{\bf{\hat{r}}}{r^2} [/itex] term to introduce a delta function, rather than the -3 power I've found myself with?

Related Threads on Integrating charge density