Different descriptions of electrostatic energy

Mexicorn
Messages
2
Reaction score
0
This is actually a question pertaining to a paper I'm trying to understand (PRB 73, 115407 (2006)), but I decided to put it here just to be safe.

Homework Statement


The paper I'm reading involves starting with an electrostatic energy contribution, and rewriting it with a green's function solution to the Poisson equation rather than the standard 1/r form to incorporate a dielectric. The step I'm missing, though, is in how the author formulates the energy contribution at first.

Homework Equations


The electrostatic energy contribution to the total energy is described like so:
-\int \frac{\epsilon(r)}{8\pi}\left|\nabla V(r)\right|^{2}d^{3}r +\int \rho(r) V(r)d^{3}r.
After incorporating the green's function solution (which is not relevant to my holdup), the term becomes the more familiar (to me at least):
\frac{1}{2}\int \rho(r) V(r)d^{3}r. (the sign on this term is ambiguous as it is defined as positive in the paper and negative in an online formulation written by the same author)


The Attempt at a Solution


Starting from the top equation..
-\int \frac{\epsilon(r)}{8\pi}\left|\nabla V(r)\right|^{2}d^{3}r +\int \rho(r) V(r)d^{3}r
and plugging in the identity: \textbf{E} = - \nabla V(r), I am left with
\int \frac{\epsilon(r)}{8\pi}\nabla V(r) \cdot \textbf{E} d^{3}r+\int \rho(r) V(r)d^{3}r <br />
Then using Gauss' theorem yields:
\int \frac{\epsilon(r)}{8\pi}\nabla \cdot (V(r)\textbf{E}) d^{3}r-\int \frac{\epsilon(r)}{8\pi} V(r) (\nabla \cdot \textbf{E}) d^{3}r+\int \rho(r) V(r)d^{3}r <br />
By the Divergence theorem, the first term is equivalent to \oint (V(r) \textbf{E}) \cdot d\textbf{A} which is zero under the assumption that the potential vanishes at infinity. The last step is where I am struggling. At first I thought I would merely plug in the identity \nabla \cdot \textbf{E} = 4\pi \rho, but this still leaves me with a dangling permittivity within the integral. I tried going through everything again substituting in the electric displacement instead of the field, but then my last step will have a \rho_f instead of the total charge.

Anyone know where I'm going wrong with this?
 
Physics news on Phys.org
Mexicorn said:
At first I thought I would merely plug in the identity \nabla \cdot \textbf{E} = 4\pi \rho,

The correct formula is \nabla \cdot \textbf{D} = 4\pi \rho

ehild
 
Shouldn't the form be \nabla \cdot \textbf{D} = \rho_{free} which would have me missing a part of the total charge?

Going through the steps using \textbf{D}=-\epsilon(r)\nabla V(r) gives me the following result:
-\int \frac{1}{2} V(r) \rho_{free}(r) d^{3}r + \int \rho_{total}(r) V(r) d^{3}r

I don't see how I can combine a total charge with just a free charge (coming from the dielectric displacement). This model is used where there will be dielectric areas which should accumulate bound charges...

Thank you, ehild, for taking a look at this.
 
Are you sure the charge in the second integral is not free charge?

ehild
 
Prove $$\int\limits_0^{\sqrt2/4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx = \frac{\pi^2}{8}.$$ Let $$I = \int\limits_0^{\sqrt 2 / 4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx. \tag{1}$$ The representation integral of ##\arcsin## is $$\arcsin u = \int\limits_{0}^{1} \frac{\mathrm dt}{\sqrt{1-t^2}}, \qquad 0 \leqslant u \leqslant 1.$$ Plugging identity above into ##(1)## with ##u...
Back
Top