[Click here for a PDF of this post with nicer formatting].

# Motivation.

In [1], Feynman writes the Pauli wave equation for a non-relativistic treatment of a mass in a scalar and vector potential electrodynamic field. That is

Is this amenable to Fourier transform solution like so many other PDEs? Let’s give it a try. It would also be interesting to attempt to apply such a computation to see if it is possible to calculate , and the first two derivatives of this expectation value. I would guess that this would produce the Lorentz force equation.

# Prep

## Fourier Notation.

Our transform pair will be written

## Interpretation of the squared momentum operator.

Feynman actually wrote

That notation I’m not familiar with, and I’ve written this as a plain old vector square. If were not an operator, then this would be a scalar, but as written this actually also includes a bivector term proportional to . To see that, lets expand this operator explicitly.

This anticommutator of the vector potential and the gradient is only a scalar has zero divergence. More generally, expanding by chain rules, and using braces to indicate the scope of the differential operations, we have

where is the spatial unit trivector, and .

This is assuming should be treated as a complex valued scalar, and not a complex-like geometric object of any sort. Does this bivector term have physical meaning? Should it be discarded or retained? If we assume discarded, then we really want to write the Pauli equation utilizing an explicit scalar selection, as in

Assuming that to be the case, our squared momentum operator takes the form

The Pauli equation, written out explicitly in terms of the gradient is then

### Confirmation.

Instead of guessing what Feynman means when he writes Pauli’s equation, it would be better to just check what Pauli says. In [2] he uses the more straightforward notation

for the vector potential dependent part of the Hamiltonian operator. This is just the scalar part as was guessed.

# Guts

Using the expansion 2.7 of the Pauli equation, and writing for the effective complex potential we have

Let’s now apply each of these derivative operations to our assumed Fourier solution from 2.2. Starting with the Laplacian we have

For the operator application we have

Putting both together we have

We can tidy this up slightly by completing the square, yielding

If this is to be zero for all , it seems clear that we need to be the solution of the first order non-linear PDE

Somewhere along the way this got a bit confused. Our Fourier transform function is somehow a function of not just wave number, but position, since by virtue of being a solution to a differential equation involving , and ? Can we pretend to not to have noticed this and continue on anyways? Let’s try the further simplification of the system by imposing a constraint of constant time potentials (). That allows for direct integration of the wave function’s Fourier transform

And inverse transforming this

By inserting the inverse Fourier transform of , we have the time evolution of the wave function as a convolution integral

Splitting out the convolution kernel, this takes a slightly tidier form

## Verification attempt.

If we apply the Pauli equation 1.1 to 3.18 does it produce the correct answer?

For the LHS we have

but for the RHS we have

So if it were not for the spatial dependence of and , we would have LHS equal to the RHS. It appears that ignoring the odd dependence in the differential equation definitely leads to trouble, and only works for constant potential distributions, a rather boring special case.

# References

[1] R.P. Feynman. *Quantum Electrodynamics*. Addison-Wesley Publishing Company. Reading, Massachusetts, 1961.

[2] W. Pauli. *Wave Mechanics*. Courier Dover Publications, 2000.