
Of course, these quotes mean nothing else than that the quantities proposed as a norm do not actu-
ally constitute one. They cannot have the physical meaning of a probability integral. Considering that
quantum physics is based on the notion of probability, this seems rather problematic for the physical
justification of the Bethe-Salpeter equation.
The motivation for integral equations for a multi-time wave function can also be approached from a
second angle – one that was dear to Detlef Dürr: the quantization of Wheeler-Feynman (WF) electro-
dynamics [7–11]; see [12–14] for some of Detlef Dürr’s works on the topic. This theory pursues the idea
that the ultraviolet divergence problem of classical Maxwell-Lorenz electrodynamics can be avoided by
“integrating out” the fields. The result is a dynamics where interactions between particles occur directly
and exactly when particle world-lines are light-like separated.
The discussions with Detlef Dürr about finding a suitable quantum version of that theory sparked
my personal curiosity about the subject. After studying previous proposals for quantizations of WF
electrodynamics [15–18], [19, chap. 8], which all encounter their own difficulties, it seemed to me that
integral equations for a multi-time wave function might be a more promising way forward. In [20], I laid
out how these types of integral equations make it possible to transfer the principle of direct interactions
along light cones, that forms the core of WF electrodynamics, to the quantum level. This was done
in a way that retains the Dirac-Schrödinger equation with a spin-dependent Coulomb potential as the
non-relativistic limit, thus staying close to empirically successful models. In a series of papers [21–24],
my co-authors and me were able to prove that multi-time integral equations provide a fully relativistic
and interacting quantum dynamics which does not suffer from the ultraviolet divergence problem, even
for singular light-cone interactions [24].
However, the question of probability conservation was left open and, as we have seen for the BS
equation, there is reason for concern. Equations with interaction terms which are non-local in the
configuration space of quantum mechanics do typically not imply local conservation laws. It is then at
best unclear whether global probability conservation holds true. Historically, Feynman himself saw the
problem of probability conservation as one of the central obstacles to quantizing WF electrodynamics, as
he reports in his Nobel lecture:
I found that if one generalized the action from the nice Lagrangian forms [...] to these forms
[...] then the quantities which I defined as energy, and so on, would be complex. The energy
values of stationary states wouldn’t be real and probabilities of events wouldn’t add up to
100%. – Feynman, 1965 [25]
I don’t think we have a completely satisfactory relativistic quantum-mechanical model, even
one that doesn’t agree with nature, but, at least, agrees with the logic that the sum of
probability of all alternatives has to be 100%. – Feynman, 1965 [25]
In view of these difficulties, if one is not ready to dismiss multi-time integral equations altogether, one
may conclude that the equations are not exactly the right ones and that some modification is in order.
From the point of view that the interaction term in the BS equation which quantum electrodynamics
suggests (an infinite series of Feynman diagrams in need of renormalization) is not the most simple and
natural, such a modification seems easy to accept. In addition, the argument that non-local interaction
terms usually preclude local conservation laws does not apply to global conservation laws, leaving room
for logical possibilities which may not have been sufficiently explored. I am going to adopt these positions
here. This makes it possible to prove that for certain classes of integral equations the global probability
integral is, in fact, conserved.
2 The integral equation
For simplicity, we focus on the case of N= 2 Dirac particles. Moreover, we set c=1=~. Then the class
of integral equations we shall study reads:
ψ(x1, x2) = ψfree(x1, x2) + iZd4x0
1d4x0
2G1(x1−x0
1)G2(x2−x0
2)K(x0
1, x0
2)ψ(x0
1, x0
2).(1)
Here ψ:R4×R4→C4⊗C4is a multi-time wave function with 16 complex components for two particles.1
ψfree(x1, x2)is a solution of the free Dirac equation in x1, x2, i.e.:
(−iγµ
k∂xµ
k+mk)ψfree(x1, x2)=0, k = 1,2.(2)
1One can also study (1) on a sub-domain of R4×R4, e.g., on the set Sof space-like configurations.
2