Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I've some difficulty deriving this equation from Jackson electrodynamics (The equation after 1.30)

$$\nabla^2 \Phi_a\left({\textbf{x}}\right)=-\frac{1}{\epsilon_0}\int_{0}^{R} \frac{3a^2}{\left(r^2+a^2\right)^{\frac{5}{2}}} \left[\rho\left(\textbf{x}\right)+\frac{r^2}{6}\nabla^2\rho+...\right]r^2dr+O\left(a^2\right).$$

It follows from

$$\nabla^2 \Phi_a\left({\textbf{x}}\right)=\frac{1}{4\pi\epsilon_0}\int\rho\left(\textbf{x}'\right)\left(\frac{3a^ 2}{\left(r^2+a^2\right)^{\frac{5}{2}}}\right)d^3x'$$

by Taylor expansion of the charge density $\rho(x^{\prime})$ about the point $x$ and then doing angular integration . I've difficulty working out the details of this. Here $r=x-x^{\prime}$ and $\Phi_a$ is a function from $R^{3}$ to $R^{1}$ "the electric potential".

share|cite|improve this question
As a place to start, you might look at . It's a bit light on the details, but a lot of the related entries in the right column looked promising at first glance. – David H Mar 26 '13 at 12:59

The Taylor expansion of $\rho(\mathbf x')$ about a point $\mathbf x$ is $$ \rho(\mathbf x') = \rho(\mathbf x) + (\mathbf x'-\mathbf x)\cdot\nabla\rho(\mathbf x)+ \frac{1}{2}\sum_{ij}(x_i'-x_i)(x_j'-x_j)\frac{\partial^2\rho}{\partial x_i'\partial x_j'}(\mathbf x) + \cdots $$ The integral we are attempting to perform is $$ \int_{\mathbb R^3} d^3\mathbf x'\,\rho(\mathbf{x})f(r), \qquad f(r) = \frac{3a^2}{(r^2+a^2)^{5/2}}, \qquad r = |\mathbf x' - \mathbf x| $$ Integration is linear, so (up to some technicalities about integration of a power series term-by-term that I'll ignore) we can Taylor expand the charge density and integrate the Taylor expansion term-by-term. The first order term in the Taylor expansion gives $$ I_1(\mathbf x) = \int d^3\mathbf x'\,f(r)(\mathbf x'-\mathbf x)\cdot\nabla\rho(\mathbf x) $$ I claim that this vanishes. To see this, note first that without loss of generality (actually the fact that there is no loss of generality here must be proven, but it's relatively straightforward, and I'll leave it to you to think about) we can choose our axes so that $\nabla\rho(\mathbf x)$ points in the positive $z$-direction, namely $\nabla\rho(\mathbf x) = |\nabla\rho(\mathbf x)|\hat{\mathbf z}$. Next, make the change of variables $\mathbf r = \mathbf x - \mathbf x'$ then $d^3 \mathbf x'=d^3\mathbf r$. Putting this together gives \begin{align} I_1(\mathbf x) &= \int d^3\mathbf r \,f(r)\,\mathbf r\cdot\nabla\rho(\mathbf x) \\ &=\int dr\,d\theta\,d\phi\,r^2\sin\theta f(r) r|\nabla\rho(\mathbf x)|\cos\theta \\ &=|\nabla\rho(\mathbf x)|\int_0^\infty dr\,r^3f(r)\int_0^{2\pi} d\phi\int_0^{\pi} d\theta\,\sin^2\theta\cos\theta \end{align} but notice that the $\theta$ intgeral itself vanishes; $$ \int_0^\pi d\theta\,\sin^2\theta\cos\theta =-\int_{-\pi/2}^{\pi/2}d\eta\,\cos^2\eta\sin\eta = 0 $$ The first equality came from the change of variables $\eta = \theta-\pi/2$, and the second equality comes from integrating an odd function against an even function on a symmetric interval containing zero. It follows that $$ I_1(\mathbf x) = 0 $$ It remains to show that the integral $I_2$ that comes from integrating the second order ter in the Taylor expansion gives the $\frac{r^2}{6}\nabla^2\rho$ given in Jackson. I'll leave this to you. If you're really struggling, I'll gladly give some more hints.

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.