[Click here for a PDF of this post with nicer formatting and figures if the post had any (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]

# Motivation.

Part of a classical mechanics problem set was to look at what portions of momentum and angular momentum are conserved for various fields. Since this was also a previous midterm question, I’m expecting that some intuition was expected to be used determine the form of the Lagrangians, with not much effort on finding the precise form of the potentials. Here I try to calculate one such potential explicitly.

Oddly, I can explicitly calculate the potential for the infinite homogeneous plane if I start with the force and then calculate the potential, but if I start with the potential the integral also diverges? That doesn’t make any sense, so I’m wondering if I’ve miscalculated. I’ve tried a few different ways below, but can’t get a non-divergent result if I start from the integral definition of the potential instead of deriving the potential from the force after adding up all the directional force contributions (where I find of course that all the component but the one perpendicular to the plane cancel out).

# Forces and potential for an infinite homogeneous plane.

## Calculating the potential from the force.

For the plane, with as the distance from the plane, and coordinates in the plane, as illustrated in figure (\ref{fig:infiniteSheetPotentials:infiniteSheetPotentialsFig1}). The gravitational force from an element of mass on the plane is

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{infiniteSheetPotentialsFig1}

\caption{Coordinate choice for interaction with infinite plane mass distribution.}

\end{figure}

Integrating for the total force on the test mass, noting that the sinusoidal terms vanish when integrated over a interval, we have

a substitution gives us

For

For the Lagrangian problem we want the potential

or

This is a reasonable seeming answer. Our potential is of the form

with

## Calculating the potential directly.

If we only want the potential, why start with the force? We ought to be able to work with potentials directly, and write

It’s a quick calculation to verify that is correct, and we find (provided )

To verify the signs, it’s helpful to refer to the diagram (\ref{fig:infiniteSheetPotentials:infiniteSheetPotentialsFig2}) which illustrates the position vectors. Now, suppose we calculate the potential directly

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{infiniteSheetPotentialsFig2}

\caption{Direction vectors for interaction with mass distribution.}

\end{figure}

### Naive approach, with bogus result.

Our mass density as a function of is

Again we can make a tangent substitution

for

This has the same functional form as 2.4, except with

There’s one significant and irritating difference. The integral above doesn’t have a unit value, but instead diverges (with as ).

What went wrong? Trouble can be seen right from the beginning. Consider the differential form above for

This we are integrating on , so long before we make the trig substitutions we are in trouble.

### As the limit of a finite volume.

My guess is that we have to tie the limits of the width of the plane and its diameter, decreasing the thickness in proportion to the increase in the radius. Let’s try that.

With a finite cylinder of height , radius , with a measurement of the potential directly above the cylinder at height , we have

Performing the integration and substituting

we have

For this integral we find

Taking limits and the terms

at either or tend to 1. This leaves us only with the contribution, so

With we have the structure of a differential above, but instead of expressing the derivative in the usual forward difference form

we have to use a backwards difference, which is equivalent, provided the function is continuous at

We can then form the differential

We also want to express the charge density in terms of surface charge density , and note that these are related by .

or

Looking at this result, we have the same divergent integration result as in the first attempt, and the reason for this is clear after some reflection. The limiting process for the radius and the thickness of the slice were allowed to complete independently. Before taking limits we had

Consider a similar, but slightly more general case, where we evaluate the limit

where , so that

So even without evaluating the integral we expect that we’ll have

This matches what was obtained in 2.19 by brute forcing the integral with Mathematica, and then evaluating the limit the hard way. Darn.

### As the limit of a finite volume. Take II.

Let’s try once more. We’ll consider a homogeneous cylindrical volume of radius , thickness with total mass

so that the area density is

Now we’ll reduce the thickness of the volume, keeping the total mass fixed, so that

or

We wish to evaluate

Performing the integrals, (first , then ) we find

The term is clearly killed in the limit, and we have a contribution from the first term. For the difference of logarithms above can be written as

Suppose we could also validly argue that this tends to , and the difference of square roots could also be canceled. Then we would be left with just

If we also demand that , then we have a contribution from the remaining square root and are left with

This differs from the expected result by a factor of , and we’ve had to do some very fishy root taking to even get that far. Employing l’H\^opital’s rule (or letting Mathematica attempt to evaluate the limit), we get infinities for the difference of logarithm term.

So, with a lot of cheating we get a result that’s similar to the expected, but not actually a match, and even to get that we had to take the limits in an invalid way. It looks like it’s back to the drawing board, but I’m not sure how to approach it.

After thinking about it a bit, perhaps the limiting process for the width needs to be explicitly accounted for using a delta function? Perhaps a QM like treatment where we express the integral in terms of some basis and look for the resolution of identity in that resolution?