[Click here for a PDF of this post with nicer formatting].
In , 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.
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
Instead of guessing what Feynman means when he writes Pauli’s equation, it would be better to just check what Pauli says. In  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.
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
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.
 R.P. Feynman. Quantum Electrodynamics. Addison-Wesley Publishing Company. Reading, Massachusetts, 1961.
 W. Pauli. Wave Mechanics. Courier Dover Publications, 2000.