schaefera said:
I know that the energy stored in a field is the same as the potential energy of the system.
dU= \varphi dq = \varphi(σ dV)
Though my book also says that dU= (1/2) ρ \varphi dV
So I'm not sure which equation to use. (I'm using Purcell.)
Your first equation should be U=\frac{1}{2}\int \varphi dq = \frac{1}{2}\int \varphi \sigma dA
The factor of 1/2 is is best seen by first looking at the discrete case of a collection of point charges. The work done when assembling the distribution can be calculated by first positioning one charge at its final location (this takes no work since there are not yet any other charges/fields present), then moving a second charge in from infinity to its final location, while holding the first charge fixed (so work is done against the electric field of the first charge), then bringing in a third charge while holding the first two fixed (so work is done against the fields of the first two charges), and so on. You end up with a sum like
W=\frac{1}{4 \pi \epsilon_0} \sum_{i=1}^{n} \sum_{j=1}_{j > i}^{n} \frac{q_i q_r}{ | \mathbf{r}_j - \mathbf{r}_i |}
where the condition j>i ensures that you don't double-count interactions. You can simplify this by
intentionally double-counting interactions and then dividing by 2 to get
W=\frac{1}{2} \sum_{i=1}^{n} q_i \left( \sum_{j=1}_{j \neq i}^{n} \frac{1}{4 \pi \epsilon_0} \frac{ q_r} {|\mathbf{r}_j - \mathbf{r}_i|} \right) = \frac{1}{2} \sum_{i=1}^{n} q_i \varphi( \mathbf{r}_i )
where \varphi( \mathbf{r}_i ) is the potential at the location \mathbf{r}_i of the charge q_i due to all the other charges.
Generalizing to a continuous distribution with \sum \to \int, q_i \to dq gives you the correct formula.
The Attempt at a Solution
If I integrate without the 1/2, I get the correct answer (8Q^2)/(3πa).
But why don't I use the 1/2?
Also, I understand that dq=σ dA, and dA=2πr^2 dr,
No, dA= r dr d\phi. Try your integration again.