Ok, I'd like to take a stab at the reason for the factor of 2 error in my last problem.
When I calculated the field E (let's say on the positive plate), I did this based on the field from charge on the other plate (negative plate), ignoring the effect of the positive plate's charge upon the field E.
There may be good reasons for this. I invite comments. But I am going to use the opposite approach.
I am going to include the effect of the positive plate charge upon E (correct me if I'm wrong). Now we see that the field is E = Q/(e*A) on one side of the charge (toward the other plate) and the field is zero on the other side (toward outside of the cap).
So if Q lies on the boundary between two different fields, E= Q/(e*A) and E=0, what E do we use when calculating force F=Q*E. I propose that E=0.5Q/(e*A) is the appropriate charge. This of course would take care of the factor of 2 error in previous message.
Two justifications for this approach, a simple one and a more rigorous one.
The simple justification…. the charge lies on the boundary. If we imagine the charge to be a little sphere, it is half in the field and half out of the field. Use the average field of 0.5Q/(e*A)
A little more rigorous. Move the charge in small elements dq (perhaps stored on separate elements d_A).
That way it seems clear that dq itself has no effect on the overall field (very small element of charge), but all of the remaining charge should be considered.
When I move the first element dq from d0 to 2d0…. it immediately enters an area of zero field (the rest of the positive charge lies between it and the negative charge). The field that the first element sees is zero.
When I move the last element dq from 0 to 2d0…. it is between the rest of the charge where it sees the full field E=Q/(e*A)
So the field that is seen as each charge is moved will vary linearly from 0 (for the first differential element dq) to E=Q/(e*A) (for the last element). The average field seen is halfway between these: E=0.5Q/(e*A)
If I substitute this field E=0.5Q/(e*A) in my previous message, the mechanical energy F*d is equal to the change in potential energy.
Start here:
There is An electric field E = 0.5 Q/(e*A)
This produces a force on each plate F = Q*E = 0.5Q^2/(e*A)
By separating the plates from distance d0 to distance 2*d0, I have to move them against that force F through a distance d0. The energy W is
W = Integral (F) with respect to d
W = F* d0 = 0.5*Q^2/(e*A) * d0 = 0.5*Q^2*d0/(e*A) = 0.5*Q^2/C0 = 0.5*(C0V0)^2/C0 = 0.5*C0V0^2 = 0.5*W0.
I have to put an energy W0 into the system to move it from d0 to 2*d0. This increases stored energy from W0 to 2*W0.