Take the 2-minute tour ×
Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. It's 100% free, no registration required.

Consider a disk with a radius $R$ (I'll use $R=1$ at various points here) that has a constant surface charge density $\sigma$. Unlike the similar problem of the field in the vicinity of a infinitely thin ring, the field directly above the disk is very well behaved. Like the case of an infinite sheet of constant surface charge, the field limits to a constant as you get closer to the surface, and a point on the surface has a defined electrical potential.

This recent paper gives equations for the specific case of points on the plane of the disk. Since we have rotational symmetry, there is only one variable, $\rho$, the distance from the axis of rotational symmetry. Graph from the paper:

Potential

As you can imagine, the point of the highest slope is the edge of the disk. Is the slope (field) at this point infinite? From what I can make out of their equations, signs seem to point to "yes". The Elliptic functions are tricky, however, and I don't trust my conclusions. Up to the point $\rho=R$, the form of the above function is:

$$ V(\rho) = V(0) \frac{2}{\pi} E(\rho)$$

Here, $E(\rho)$ is an elliptic function, and is different from the paper, because it uses the "Matlab" convention while I use the "Wikipedia" convention. The slope of E() doesn't look convincingly upright, but the derivative seems to indicate it should be. I'm also troubled by the fact that I can't seem to numerically observe a divergence at $\rho=0.9999999999$.

The physics at play are even more perplexing to me. I expected that a 2D surface, any 2D surface, would have a finite field on its face. Not so for a line, I understand that. But the ring that defines the edge of the disk only has a differential amount of charge on it.

Consider another thought experiment: start at $\rho=1$ and $z=0$, I'll denote it $(1,0)$. Now move upward, but keep $\rho$ the same. If you're at $(1,0.001)$ you should have a finite field. But what happens as $z\rightarrow 0$? Does it just jump from a defined finite value to infinity? What on Earth is going on with the electric field in the vicinity of this edge, or ANY edge of a 2D surface for that matter?

share|improve this question

1 Answer 1

up vote 6 down vote accepted

Yes, the field is infinite, but it is only log divergent near the plate, so that it is hard to see the divergence numerically. You can see this easily by solving the problem of a uniformly charged infinite plate, which is a 2d problem. Here the charges are uniform along the negative real axis, where the 2d space is imagined to be the complex plane.

This problem can be understood as follows: the 2d electrostatic field of a point charge at the origin, written as a map from C to C can be written in complex form as:

$$ E_x + i E_y = \frac{z} { 2\pi |z|^2 }= \frac{1}{ 2\pi \bar{z}}$$

it points radially outwards. This is a pure antiholomorphic function, except at the origin. It's more familiar to deal with holomorphic functions, so conjugate it!

$$ E_x - iE_y = E(z) = {1\over 2\pi z} $$

Now you want to superpose all the charges on the negative z axis. This is a simple integral:

$$ \int_{-\infty}^0 E(z-a) da = - {1\over 2\pi} \log(z) $$

where I threw away an infinite additive constant (you should think of this as calculating the potential difference between the point z=1 and any other point). This is the function with a given fixed cut discontinuity on the negative real axis.

So at the point $r,\theta$, the electric field is

$$ E_x = - {\rho\over 2\pi} \log(r) $$

$$ E_y = {\rho \theta\over 2\pi} $$

Where I have restored the $\rho$. The part in the y-direction is finite, as your intuition says--- the discontinuity is equal to the charge density (this charge density is the cut discontinuity of the electric field analytic function, which is a way of making it obvious that the electric field goes as the log--- the log function as a constant cut discontinuity). The divergent part is in the x direction, and it is only invisible in the bulk disk because when you get close to the surface, you have cancellations from the left and from the right that wash it out.

So the answer is yes, the E field is log divergent, but only the component in the plane of the disk pointing out. The solution of the disk asymptotes to the plane solution in the near disk limit.

share|improve this answer
    
I fixed about 5 TeX errors in a formula Ron. Good answer. –  Luboš Motl Oct 16 '12 at 10:01
    
What does this mean: "The divergent part is in the x direction, and it is only invisible in the bulk disk"? The divergence doesn't seem to disappear approaching from the left or right (on disk/off disk). The paper I referenced gives two equations, and in both the derivative gives an EllipticK divergence. Granted, I can't guarantee the behavior would be the same for the straight edge of an infinite sheet, but it would be truly strange if the two geometries behaved differently. –  AlanSE Oct 16 '12 at 12:56
    
The divergence is exactly the same near the disk and the plane, that's the point of the plane solution. I gave the explicit formula for the field above. There is only a divergence on the boundary circle, and it's only in the horizontal component of the field. Forget "EllipticK", it's just a stupid log divergence exactly as I gave above. –  Ron Maimon Oct 16 '12 at 15:33
    
@LubošMotl: Thanks! Although, if you notice, I try to always use Knuths "\over" rather than Lamport's "\frac". I like Knuth's way better. –  Ron Maimon Oct 16 '12 at 15:33
    
If I understand correctly, you are modeling the semi-infinite plane as a "curtain" of line charges, superimposing the fields of those line charges. However, that model breaks down just where you don't want it to, at the edge of the plane. Imagine a gaussian tube about the outermost line charge: you can shrink its radius as much as you like without decreasing the enclosed charge, so the field diverges (because the local charge density is infinite). With a surface charge density, on the other hand, the charge enclosed gets smaller as the tube shrinks, and the field stays finite (I think). –  Art Brown Oct 16 '12 at 17:15

Your Answer

 
discard

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.