R.P. Feynman

(Received June 8, 1950)

The validity of the rules given in previous papers for the solution of problems in quantum electrodynamics is established. Starting with Fermi’s formulation of the field as a set of harmonic oscillators, the effect of the oscillators is integrated out in the Lagrangian form of quantum mechanics. There results an expression for the effect of all virtual photons valid to all orders in ${e}^{2}/\mathrm{\hslash c}$. It is shown that evaluation of this expression as a power series in ${e}^{2}/\mathrm{\hslash c}$ gives just the terms expected by the aforementioned rules.

In addition, a relation is established between the amplitude for a given process in an arbitrary unquantized potential and in a quantum electrodynamical field. This relation permits a simple general statement of the laws of quantum electrodynamics.

A description, in Lagrangian quantum-mechanical form, of particles satisfying the Klein-Gordon equation is given in an Appendix. It involves the use of an extra parameter analogous to proper time to describe the trajectory of the particle in four dimensions.

A second Appendix discusses in the special case of photons, the problem of finding what real processes are implied by the formula for virtual processes.

Problems of the divergences of electrodynamics are not discussed.

In two previous papers^{1} rules were given for the calculation of the matrix element for any process in electrodynamics, to each order in ${e}^{2}/\mathrm{\hslash c}$. No complete proof of the equivalence of these rules to the conventional electrodynamics was given in these papers. Secondly, no closed expression was given valid to all orders in ${e}^{2}/\mathrm{\hslash c}$. In this paper these formal omissions will be remedied.^{2}

In paper II it was pointed out that for many problems in electrodynamics the Hamiltonian method is not advantageous, and might be replaced by the over-all space-time point of view of a direct particle interaction. It was also mentioned that the Lagrangian form of quantum mechanics^{3} was useful in this connection. The rules given in paper II were, in fact, first deduced in this form of quantum mechanics. We shall give this derivation here.

The advantage of a Lagrangian form of quantum mechanics is that in a system with interacting parts it permits a separation of the problem such that the motion of any part can be analyzed or solved first, and the results of this solution may then be used in the solution of the motion of the other parts. This separation is especially useful in quantum electrodynamics which represents the interaction of matter with the electromagnetic field. The electromagnetic field is an especially simple system and its behavior can be analyzed completely. What we shall show is that the net effect of the field is a delayed interaction of the particles. It is possible to do this easily only if it is not necessary at the same time to analyze completely the motion of the particles. The only advantage in our problems of the form of quantum mechanics in C is to permit one to separate these aspects of the problem. There are a number of disadvantages, however, such as a lack of familiarity, the apparent (but not real) necessity for dealing with matter in non-relativistic approximation, and at times a cumbersome mathematical notation and method, as well as the fact that a great deal of useful information that is known about operators cannot be directly applied.

It is also possible to separate the field and particle aspects of a problem in a manner which uses operators and Hamiltonians in a way that is much more familiar. One abandons the notation that the order of action of operators depends on their written position on the paper and substitutes some other convention (such that the order of operators is that of the time to which they refer). The increase in manipulative facility which accompanies this change in notation makes it easier to represent and to analyze the formal problems in electrodynamics. The method requires some discussion, however, and will be described in a succeeding paper. In this paper we shall give the derivations of the formulas of II by means of the form of quantum mechanics given in C.

The problem of interaction of matter and field will be analyzed by first solving for the behavior of the field in terms of the coordinates of the matter, and finally discussing the behavior of the matter (by matter is actually meant the electrons and positrons). That is to say, we shall first eliminate the field variables from the equations of motion of the electrons and then discuss the behavior of the electrons. In this way all of the rules given in the paper II will be derived.

Actually, the straightforward elimination of the field variables will lead at first to an expression for the behavior of an arbitrary number of Dirac electrons. Since the number of electrons might be infinite, this can be used directly to find the behavior of the electrons according to hole theory by imagining that nearly all the negative energy states are occupied by electrons. But, at least in the case of motion in a fixed potential, it has been shown that this hole theory picture is equivalent to one in which a positron is represented as an electron whose space-time trajectory has had its time direction reversed. To show that this same picture may be used in quantum electrodynamics when the potentials are not fixed, a special argument is made based on a study of the relationship of quantum electrodynamics to motion in a fixed potential. Finally, it is pointed out that this relationship is quite general and might be used for a general statement of the laws of quantum electrodynamics.

Charges obeying the Klein-Gordon equation can be analyzed by a special formalism given in Appendix A. A fifth parameter is used to specify the four-dimensional trajectory so that the Lagrangian form of quantum mechanics can be used. Appendix B discusses in more detail the relation of real and virtual photon emission. An equation for the propagation of a self-interacting electron is given in Appendix C.

In the demonstration which follows we shall restrict ourselves temporarily to cases in which the particle’s motion is non-relativistic, but the transition of the final formulas to the relativistic case is direct, and the proof could have been kept relativistic throughout.

The transverse part of the electromagnetic field will be represented as an assemblage of independent harmonic oscillators each interacting with the particles, as suggested by Fermi.^{4} We use the notation of Heitler.^{5}

The Hamiltonian for a set of non-relativistic particles interacting with radiation is, classically, $H={H}_{p}+{H}_{I}+{H}_{c}+{H}_{t\text{}r}$, where ${H}_{p}+{H}_{I}={\sum}_{n}\frac{1}{2}{{m}_{n}^{-1}\left({p}_{n}-{e}_{n}{A}^{t\text{}r}\left({x}_{n}\right)\right)}^{2}$ is the Hamiltonian of the particles of mass ${m}_{n}$, charge ${e}_{n}$, coordinate ${x}_{n}$, and momentum ${p}_{n}$, and their interaction with the transverse part of the electromagnetic field. This field can be expanded into plane waves

$${A}^{tr}\left(x\right)={\left(8\pi \right)}^{\frac{1}{2}}{\sum}_{K}\left[{e}_{1}\left({q}_{K}^{\left(1\right)}\mathrm{cos}\left(K\xb7x\right)+{q}_{K}^{\left(3\right)}\mathrm{sin}\left(K\xb7x\right)\right)+{e}_{2}\left({q}_{K}^{\left(2\right)}\mathrm{cos}\left(K\xb7x\right)+{q}_{K}^{\left(4\right)}\mathrm{sin}\left(K\xb7x\right)\right)\right]$$ | (1) |

where ${e}_{1}$ and ${e}_{2}$ are two orthogonal polarization vectors at right angles to the propagation vector $K$, magnitude $k$. The sum over $K$ means, if normalized to unit volume, $\frac{1}{2}\int {d}^{3}K/8{\pi}^{3}$, and each ${q}_{K}^{\left(r\right)}$, can be considered as the coordinate of a harmonic oscillator. (The factor $\frac{1}{2}$ arises for the mode corresponding to $K$ and to $-K$ is the same.) The Hamiltonian of the transverse field represented as oscillators is

$${H}_{t\text{}r}=\frac{1}{2}\sum _{K}^{}\sum _{r=1}^{4}\left({\left({\U0001d4c5}_{K}^{\left(r\right)}\right)}^{2}+{{k}^{2}\left({q}_{K}^{\left(r\right)}\right)}^{2}\right)$$

where ${\U0001d4c5}_{K}^{\left(r\right)}$ is the momentum conjugate to ${q}_{K}^{\left(r\right)}$. The longitudinal part of the field has been replaced by the Coulomb interaction^{6} ,

$${H}_{c}=\frac{1}{2}{\sum}_{n}{\sum}_{m}{e}_{n}{e}_{m}/{r}_{n\text{}m}$$

where ${r}_{n\text{}m}^{2}={\left({x}_{n}-{x}_{m}\right)}^{2}$. As is well known, when this Hamiltonian is quantized one arrives at the usual theory of quantum electrodynamics. To express these laws of quantum electrodynamics one can equally well use the Lagrangian form of quantum mechanics to describe this set of oscillators and particles. The classical Lagrangian equivalent to this Hamiltonian is $L={L}_{p}+{L}_{I}+{L}_{c}+{L}_{tr}$ where

$$\begin{array}{rl}{L}_{p}& =\frac{1}{2}{\sum}_{n}{m}_{n}{x}_{n}^{2}\\ {L}_{I}\text{}& ={\sum}_{n}{e}_{n}{{x}^{.}}_{n}\xb7{A}^{tr}\left({x}_{n}\right)\\ {L}_{t\text{}r}\text{}& =\frac{1}{2}{\sum}_{K}{\sum}_{r}\left({\left({\stackrel{.}{q}}_{K}^{\left(r\right)}\right)}^{2}-{{k}^{2}\left({q}_{K}^{\left(r\right)}\right)}^{2}\right)\\ {L}_{c}\text{}& =-\frac{1}{2}{\sum}_{n}{\sum}_{m}{e}_{n}{e}_{m}/{r}_{m\text{}n}.\end{array}$$ | (2) |

When this Lagrangian is used in the Lagrangian forms of quantum mechanics of C, what it leads to is, of course, mathematically equivalent to the result of using the Hamiltonian $H$ in the ordinary way, and is therefore equivalent to the more usual forms of quantum electrodynamics (at least for non-relativistic particles). We may, therefore, proceed by using this Lagrangian form of quantum electrodynamics, with the assurance that the results obtained must agree with those obtained from the more usual Hamiltonian form.

The Lagrangian enters through the statement that the functional which carries the system from one state to another is $\mathrm{exp}\left(iS\right)$ where

$$S=\int Ldt={S}_{p}+{S}_{I}+{S}_{c}+{S}_{t\text{}r}.$$ | (3) |

The time integrals must be written as Riemann sums with some care; for example,

$${S}_{I}={\sum}_{n}\int {e}_{n}{{x}^{.}}_{n}\left(t\right)\xb7{A}^{tr}\left({x}_{n}\left(t\right)\right)d\text{}t$$ | (4) |

becomes according to C, Eq. (19)

$${S}_{I}={\sum}_{n}{\sum}_{i}\frac{1}{2}{e}_{n}\left({x}_{n,i+1}-{x}_{n,i}\right)\xb7\left({A}^{tr}\left({x}_{n,i+1}\right)+{A}^{tr}\left({x}_{n,i}\right)\right)$$ | (5) |

so that the velocity ${{x}^{.}}_{n,i}$ which multiplies ${A}^{tr}\left({x}_{n,i}\right)$ is

$${{x}^{.}}_{n,i}=\frac{1}{2}{\epsilon}^{-1}\left({x}_{n,i+1}-{x}_{n,i}\right)+\frac{1}{2}{\epsilon}^{-1}\left({x}_{n,i}-{x}_{n,i-1}\right).$$ | (6) |

In the Lagrangian form it is possible to eliminate the transverse oscillators as is discussed in C, Section 13. One must specify, however, the initial and final state of all oscillators. We shall first choose the special, simple case that all oscillators are in their ground states initially and finally, so that all photons are virtual. Later we do the more general case in which real quanta are present initially or finally. We ask, then, for the amplitude for finding no quanta present and the particles in state ${\chi}_{{t}^{\u2033}}$, at time ${t}^{\u2033}$, if at time ${t}^{\prime}$ the particles were in state ${\psi}_{{t}^{\prime}}$ and no quanta were present.

The method of eliminating field oscillators is described in Section 13 of C. We shall simply carry out the elimination here using the notation and equations of C. To do this, for simplicity, we first consider in the next section the case of a particle or a system of particles interacting with a single oscillator, rather than the entire assemblage of the electromagnetic field.

We consider a harmonic oscillator, coordinate $q$, Lagrangian $L=\frac{1}{2}\left({\stackrel{.}{q}}^{2}-{\omega}^{2}{q}^{2}\right)$ interacting with a particle or system of particles, action ${S}_{p}$, through a term in the Lagrangian $q\left(t\right)\gamma \left(t\right)$ where $\gamma \left(t\right)$ is a function of the coordinates (symbolized as $x$) of the particle. The precise form of $\gamma \left(t\right)$ for each oscillator of the electromagnetic field is given in the next section. We ask for the amplitude that at some time ${t}^{\u2033}$ the particles are in state ${\chi}_{{t}^{\u2033}}$, and the oscillator is in, say, an eigenstate $m$ of energy $\omega \left(m+\frac{1}{2}\right)$ (units are chosen such that $\hslash =c=1$) when it is given that at a previous time ${t}^{\prime}$ the particles were in state ${\psi}_{{t}^{\prime}}$. and the oscillator in $n$. The amplitude for this is the transition amplitude [see C, Eq. (61)]

$${\langle {\chi}_{{t}^{\u2033}}{\phi}_{m}\left|1\right|{\psi}_{{t}^{\prime}}{\phi}_{n}\rangle}_{{S}_{p}+{S}_{0}+{S}_{I}}=\int \int {{\chi}_{{t}^{\u2033}}}^{*}\left({x}_{{t}^{\u2033}}\right){{\phi}_{m}}^{*}\left({q}_{{t}^{\u2033}}\right)\mathrm{exp}i\left({S}_{p}+{S}_{0}+{S}_{I}\right)$$ | (7) |

$$\xb7{\phi}_{n}\left({q}_{{t}^{\prime}}\right){\psi}_{{t}^{\prime}}\left({x}_{{t}^{\prime}}\right)d{x}_{{t}^{\u2033}}d{x}_{{t}^{\prime}}d{q}_{{t}^{\u2033}}d{q}_{{t}^{\prime}}\mathcal{D}x\left(t\right)\mathcal{D}q\left(t\right)$$ | (7) |

where $x$ represents the variables describing the particle, ${S}_{p}$ is the action calculated classically for the particles for a given path going from coordinate ${x}_{{t}^{\prime}}$ at ${t}^{\prime}$ to ${x}_{{t}^{\u2033}}$ at ${t}^{\u2033}$, ${S}_{0}$ is the action $\int \frac{1}{2}\left({\stackrel{.}{q}}^{2}-{\omega}^{2}{q}^{2}\right)dt$ for any path of the oscillator going from ${q}_{{t}^{\prime}}$ at ${t}^{\prime}$ to ${q}_{{t}^{\u2033}}$ at ${t}^{\u2033}$, while

$${S}_{I}=\int q\left(t\right)\gamma \left(t\right)dt,$$ | (8) |

the action of interaction, is a functional of both $q\left(t\right)$ and $x\left(t\right)$, the paths of oscillator and particles. The symbols $\mathcal{D}x\left(t\right)$ and $\mathcal{D}q\left(t\right)$ represent a summation over all possible paths of particles and oscillator which go between the given end points in the sense defined in C, Eq. (9). (That is, assuming time to proceed in infinitesimal steps, $\epsilon $, an integral over all values of the coordinates $x$ and $q$ corresponding to each instant in time, suitably normalized.)

The problem may be broken in two. The result can be written as an integral over all paths of the particles only, of $\left(\mathrm{exp}i{S}_{p}\right)\xb7{G}_{m\text{}n}$ :

$${\langle {\chi}_{{t}^{\u2033}}{\phi}_{m}\left|1\right|{\psi}_{{t}^{\prime}}{\phi}_{n}\rangle}_{{S}_{p}+{S}_{0}+{S}_{I}}={\langle {\chi}_{{t}^{\u2033}}\left|{G}_{m\text{}n}\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}$$ | (9) |

where ${G}_{m\text{}n}$ is a functional of the path of the particles alone (since it depends on $\gamma \left(t\right)$) given by

$${G}_{m\text{}n}={\langle {\phi}_{m}\left|\mathrm{exp}i\int q\left(t\right)\gamma \left(t\right)dt\right|{\phi}_{n}\rangle}_{{S}_{0}}=\int {{\phi}_{m}}^{*}\left({q}_{{t}^{\u2033}}\right)\mathrm{exp}i\left({S}_{0}+{S}_{I}\right){\phi}_{n}\left({q}_{{t}^{\prime}}\right)d{q}_{{t}^{\prime}}d{q}_{{t}^{\u2033}}\mathcal{D}q\left(t\right)\text{}=\int {{\phi}_{m}}^{*}\left({q}_{j}\right)\mathrm{exp}i\epsilon \sum _{i=0}^{j-1}\left[\frac{1}{2}{{\epsilon}^{-2}\left({q}_{i+1}-{q}_{i}\right)}^{2}-\frac{1}{2}{\omega}^{2}{{q}_{i}}^{2}+{q}_{i}{\gamma}_{i}\right]\text{}\xb7{\phi}_{n}\left({q}_{0}\right)d{q}_{0}{a}^{-1}d{q}_{1}{a}^{-1}d{q}_{2}\cdots {a}^{-1}d{q}_{j}$$ | (10) |

where we have written the $\mathcal{D}q\left(t\right)$ out explicitly (and have set $a={\left(2\pi i\epsilon \right)}^{\frac{1}{2}}$, ${t}^{\u2033}-{t}^{\prime}=j\epsilon $, ${q}_{{t}^{\prime}}={q}_{0}$, ${q}_{{t}^{\u2033}}={q}_{j}$). The last form can be written as

$${G}_{m\text{}n}=\int {{\phi}_{m}}^{*}\left({q}_{j}\right)k\left({q}_{j},{t}^{\u2033};{q}_{0},{t}^{\prime}\right){\phi}_{n}\left({q}_{0}\right)d{q}_{0}d{q}_{j}$$ | (11) |

where $k\left({q}_{j},{t}^{\u2033};{q}_{0},{t}^{\prime}\right)$ is the kernel [as in I, Eq. (2)] for a forced harmonic oscillator giving the amplitude for arrival at ${q}_{j}$ at time ${t}^{\u2033}$ if at time ${t}^{\prime}$ it was known to be at ${q}_{0}$. According to C it is given by

$$k\left({q}_{j},{t}^{\u2033};{q}_{0},{t}^{\prime}\right)={\left(2\pi i{\omega}^{-1}\mathrm{sin}\omega \left({t}^{\u2033}-{t}^{\prime}\right)\right)}^{-\frac{1}{2}}\mathrm{exp}iQ\left({q}_{j},{t}^{\u2033};{q}_{0},{t}^{\prime}\right)$$ | (12) |

where $Q\left({q}_{j},{t}^{\u2033};{q}_{0},{t}^{\prime}\right)$ is the action calculated along the classical path between the end points ${q}_{j},{t}^{\u2033};{q}_{0},{t}^{\prime}$, and is given explicitly in C.^{7} It is

$$Q=\frac{\omega}{2\mathrm{sin}\omega \left({t}^{\u2033}-{t}^{\prime}\right)}\left[\left({{q}_{j}}^{2}+{{q}_{0}}^{2}\right)\mathrm{cos}\omega \left({t}^{\u2033}-{t}^{\prime}\right)-2{q}_{j}{q}_{0}+\frac{2{q}_{j}}{\omega}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}\gamma \left(t\right)\mathrm{sin}\omega \left(t-{t}^{\prime}\right)dt\phantom{\rule{0ex}{0ex}}+\frac{2{q}_{0}}{\omega}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}\gamma \left(t\right)\mathrm{sin}\omega \left({t}^{\u2033}-t\right)dt-\frac{2}{{\omega}^{2}}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}{\int}_{{t}^{\prime}}^{t}\gamma \left(t\right)\gamma \left(s\right)\mathrm{sin}\omega \left({t}^{\u2033}-t\right)\mathrm{sin}\omega \left(s-{t}^{\prime}\right)dsdt\right]$$ | (13) |

The solution of the motion of the oscillator can now be completed by substituting (12) and (13) into (11) and performing the integrals. The simplest case is for $m,n=0$for which case^{8}

$${\phi}_{0}\left({q}_{0}\right)={\left(\omega /\pi \right)}^{\frac{1}{2}}\mathrm{exp}\left(-\frac{1}{2}\omega {{q}_{0}}^{2}\right)\mathrm{exp}\left(-\frac{1}{2}i\omega {t}^{\prime}\right)$$

so that the integrals on ${q}_{0}$, ${q}_{j}$ are just Gaussian integrals. There results

$${G}_{0\text{}0}=\mathrm{exp}\left(-\frac{1}{2\omega}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}{\int}_{{t}^{\prime}}^{t}\mathrm{exp}\left(-i\omega \left(t-s\right)\right)\gamma \left(t\right)\gamma \left(s\right)dtds\right)$$

a result of fundamental importance in the succeeding developments. By replacing $t-s$ by its absolute value $\left|t-s\right|$ we may integrate both variables over the entire range and divide by 2. We will henceforth make the results more general by extending the limits on the integrals from $-\infty $ to $+\infty $. Thus if one wishes to study the effect on a particle of interaction with an oscillator for just the period ${t}^{\prime}$ to ${t}^{\u2033}$ one may use

$${G}_{0\text{}0}=\mathrm{exp}\left(-\frac{1}{4\omega}{\int}_{-\infty}^{\infty}{\int}_{-\infty}^{\infty}\mathrm{exp}\left(-i\omega \left|t-s\right|\right)\gamma \left(t\right)\gamma \left(s\right)dtds\right)$$ | (14) |

imagining in this case that the interaction $\gamma \left(t\right)$ is zero outside these limits. We defer to a later section the discussion of other values of $m$, $n$.

Since ${G}_{0\text{}0}$ is simply an exponential, we can write it as $\mathrm{exp}\left(iI\right)$, consider that the complete “action” for the system of particles is $S={S}_{p}+I$ and that one computes transition elements with this “action” instead of ${S}_{p}$ (see C, Sec. 12). The functional $I$ which is given by

$$I=\frac{1}{4}i{\omega}^{-1}\int \int \mathrm{exp}\left(-i\omega \left|t-s\right|\right)\gamma \left(s\right)\gamma \left(t\right)dsdt$$ | (15) |

is complex, however; we shall speak of it as the complex action. It describes the fact that the system at one time can affect itself at a different time by means of a temporary storage of energy in the oscillator. When there are several independent oscillators with different interactions, the effect if they arc all in the lowest state at ${t}^{\prime}$ and ${t}^{\u2033}$, is the product of their separate ${G}_{0\text{}0}$ contributions. Thus the complex action is additive, being the sum of contributions like (15) for each of the several oscillators.

We can now apply these results to eliminate the transverse field oscillators of the Lagrangian (2). At first we can limit ourselves to the case of purely virtual transitions in the electromagnetic field, so that there is no photon in the field at ${t}^{\prime}$ and ${t}^{\u2033}$. That is, all of the field oscillators are making transitions from ground state to ground state.

The ${{\gamma}_{K}}^{\left(r\right)}$, corresponding to each oscillator ${{q}_{K}}^{\left(r\right)}$, is found from the interaction term ${L}_{I}$ [Eq. (2b)], substituting the value of ${A}^{tr}\left(x\right)$ given in (1). There results, for example,

$$\begin{array}{cc}{{\gamma}_{K}}^{\left(1\right)}& ={\left(8\pi \right)}^{\frac{1}{2}}\sum _{n}{e}_{n}\left({e}_{1}\xb7{x}_{n}^{.}\right)\mathrm{cos}\left(K\xb7{x}_{n}\right)\\ {{\gamma}_{K}}^{\left(3\right)}& ={\left(8\pi \right)}^{\frac{1}{2}}\sum _{n}{e}_{n}\left({e}_{1}\xb7{x}_{n}^{.}\right)\mathrm{sin}\left(K\xb7{x}_{n}\right)\end{array}$$ | (16) |

the corresponding results for ${{\gamma}_{K}}^{\left(2\right)}$, ${{\gamma}_{K}}^{\left(4\right)}$, replace ${e}_{1}$, by ${e}_{2}$.

The complex action resulting from oscillator of coordinate ${{q}_{K}}^{\left(1\right)}$, is therefore

$${{I}_{K}}^{\left(1\right)}=\frac{8\pi i}{4k}\sum _{n}\sum _{m}\int \int {e}_{n}{e}_{m}\mathrm{exp}\left(-ik\left|t-s\right|\right)\times \left({e}_{1}\xb7{x}_{n}^{.}\left(t\right)\right)\left({e}_{1}\xb7{x}_{m}^{.}\left(s\right)\right)\mathrm{cos}\left(K\xb7{x}_{n}\left(t\right)\right)\mathrm{cos}\left(K\xb7{x}_{m}\left(s\right)\right)dsdt$$

The term ${{I}_{K}}^{\left(3\right)}$ exchanges the cosines for sines, so in the sum ${{I}_{K}}^{\left(1\right)}+$${{I}_{K}}^{\left(3\right)}$, the product of the two cosines, $\mathrm{cos}A\xb7\mathrm{cos}B$ is replaced by $\left(\mathrm{cos}A\mathrm{cos}B+\mathrm{sin}A\mathrm{sin}B\right)$or $\mathrm{cos}\left(A-B\right)$. The terms ${{I}_{K}}^{\left(2\right)}+$${{I}_{K}}^{\left(4\right)}$, give the same result with ${e}_{2}$ replacing ${e}_{1}$. The sum $\left({e}_{1}\xb7V\right)\left({e}_{1}\xb7{V}^{\prime}\right)+\left({e}_{2}\xb7V\right)\left({e}_{2}\xb7{V}^{\prime}\right)$ is $V\xb7{V}^{\prime}-{k}^{-2}\left(K\xb7V\right)\left(K\xb7{V}^{\prime}\right)$ since it is the sum of the products of vector components in two orthogonal directions, so that if we add the product in the third direction (that of K) we construct the complete scalar product. Summing over all K then, since $\sum K=\frac{1}{2}\int {d}^{3}K/8{\pi}^{3}$ we find for the total complex action of all of the transverse oscillators,

$${I}_{tr}=i\sum _{n}\sum _{m}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}dt{\int}_{{t}^{\prime}}^{{t}^{\u2033}}ds\int {e}_{n}{e}_{m}\mathrm{exp}\left(-ik\left|t-s\right|\right)\times \left[{x}_{n}^{.}\left(t\right)\xb7{x}_{m}^{.}\left(s\right)-{k}^{-2}\left(K\xb7{x}_{n}^{.}\left(t\right)\right)\left(K\xb7{x}_{m}^{.}\left(s\right)\right)\right]\text{}\xb7\mathrm{cos}\left(K\xb7\left({x}_{n}^{.}\left(t\right)-{x}_{m}^{.}\left(s\right)\right)\right)\frac{{d}^{3}K}{8{\pi}^{2}k}.$$ | (17) |

This is to be added to ${S}_{p}+{S}_{c}$ to obtain the complete action of the system with the oscillators removed.

The term in $\left(K\xb7{x}_{n}^{.}\left(t\right)\right)\left(K\xb7{x}_{m}^{.}\left(s\right)\right)$ can be simplified by integration by parts with respect to $t$ and with respect to $s$ [note that $\mathrm{exp}\left(-ik\left|t-s\right|\right)$ has a discontinuous slope at $t=s$, or break the integration up into two regions]. One finds

$${I}_{tr}=R-{I}_{c}+{I}_{\mathrm{transient}}$$ | (18) |

where

$$R=-i\sum _{n}\sum _{m}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}dt{\int}_{{t}^{\prime}}^{{t}^{\u2033}}ds\int {e}_{n}{e}_{m}\times \mathrm{exp}\left(-ik\left|t-s\right|\right)\left(1-{{x}^{.}}_{n}\left(t\right)\xb7{{x}^{.}}_{m}\left(s\right)\right)\phantom{\rule{0ex}{0ex}}\xb7\mathrm{cos}K\xb7\left({x}_{n}\left(t\right)-{x}_{m}\left(s\right)\right){d}^{3}K/8{\pi}^{2}k$$ | (19) |

and

$${I}_{c}=-\sum _{n}\sum _{m}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}dt\int {e}_{n}{e}_{m}\mathrm{cos}K\xb7\left({x}_{n}\left(t\right)-{x}_{m}\left(t\right)\right){d}^{3}K/4{\pi}^{2}{k}^{2}$$ | (20) |

comes from the discontinuity in slope of $\mathrm{exp}\left(-ik\left|t-s\right|\right)$ at $t=s$. Since

$$\int \mathrm{cos}\left(K\xb7R\right){d}^{3}K/4{\pi}^{2}{k}^{2}={\int}_{0}^{\infty}{\left(k\text{}r\right)}^{-1}\mathrm{sin}\left(k\text{}r\right)dk/\pi ={\left(2r\right)}^{-1}$$

this term ${I}_{c}$ just cancels the Coulomb interaction term ${S}_{c}=\int {L}_{c}dt$. The term

$${I}_{\mathrm{transient}}=-\sum _{n}\sum _{m}{e}_{n}{e}_{m}\int \frac{{d}^{3}K}{4{\pi}^{2}{k}^{2}}\times \left\{{\int}_{{t}^{\prime}}^{{t}^{\u2033}}\left[\mathrm{exp}\left(-ik\left({t}^{\u2033}-t\right)\right)\mathrm{cos}K\xb7\left({x}_{n}\left({t}^{\u2033}\right)-{x}_{m}\left(t\right)\right)\phantom{\rule{0ex}{0ex}}+\mathrm{exp}\left(-ik\left(t-{t}^{\prime}\right)\right)\mathrm{cos}K\xb7\left({x}_{n}\left(t\right)-{x}_{m}\left({t}^{\prime}\right)\right)\right]dt+{\left(2k\right)}^{-1}i\left[\mathrm{cos}K\xb7\left({x}_{n}\left({t}^{\u2033}\right)-{x}_{m}\left({t}^{\u2033}\right)\right)+\mathrm{cos}K\xb7\left({x}_{n}\left({t}^{\prime}\right)-{x}_{m}\left({t}^{\prime}\right)\right)-2\mathrm{exp}\left(-ik\left({t}^{\u2033}-{t}^{\prime}\right)\right)\mathrm{cos}K\xb7\left({x}_{n}\left({t}^{\prime}\right)-{x}_{m}\left({t}^{\u2033}\right)\right)\right]\right\}$$ | (21) |

is one which comes from the limits of integration at ${t}^{\prime}$ and ${t}^{\u2033}$, and involves the coordinates of the particle at either one of these times or the other. If ${t}^{\prime}$ and ${t}^{\u2033}$ are considered to be exceedingly far in the past and future, there is no correlation to be expected between these temporally distant coordinates and the present ones, so the effects of ${I}_{\mathrm{transient}}$ will cancel out quantum mechanically by interference. This transient was produced by the sudden turning on of the interaction of field and particles at ${t}^{\prime}$ and its sudden removal at ${t}^{\u2033}$. Alternatively we can imagine the charges to be turned on after ${t}^{\prime}$ adiabatically and turned off slowly before ${t}^{\u2033}$ (in this case, in the term ${L}_{c}$, the charges should also be considered as varying with time). In this case, in the limit, ${I}_{\mathrm{transient}}$ is zero.^{9} Hereafter we shall drop the transient term and consider the range of integration of $t$ to be from $-\infty $ to $+\infty $, imagining, if one needs a definition, that the charges vary with time and vanish in the direction of either limit.

To simplify $R$ we need the integral

$$J=\int \mathrm{exp}\left(-ik\left|t\right|\right)\mathrm{cos}\left(K\xb7R\right){d}^{3}K/8{\pi}^{2}k={\int}_{0}^{\infty}\mathrm{exp}\left(-ik\left|t\right|\right)\mathrm{sin}\left(k\text{}r\right)dk/2\pi r$$ | (22) |

where $r$ is the length of the vector $R$. Now

$${\int}_{0}^{\infty}\mathrm{exp}\left(-ikx\right)dk=\underset{\epsilon \to 0}{\mathrm{lim}}\left(-{i\left(x-i\epsilon \right)}^{-1}\right)=-i{x}^{-1}+\pi \delta \left(x\right)=\pi {\delta}_{+}\left(x\right)$$

where the equation serves to define ${\delta}_{+}\left(x\right)$ [as in II, Eq. (3)]. Hence, expanding $\mathrm{sin}\left(kr\right)$ in exponentials find

$$J=-{\left(4\pi r\right)}^{-1}\left({\left(\left|t\right|-r\right)}^{-1}-{\left(\left|t\right|+r\right)}^{-1}\right)+{\left(4ir\right)}^{-1}\left(\delta \left(\left|t\right|-r\right)-\delta \left(\left|t\right|+r\right)\right)$$ | (23) |

$$=-{\left(2\pi \right)}^{-1}{\left({t}^{2}-{r}^{2}\right)}^{-1}+{\left(2i\right)}^{-1}\delta \left({t}^{2}-{r}^{2}\right)$$ | (23) |

$$=-\frac{1}{2}i{\delta}_{+}\left({t}^{2}-{r}^{2}\right)$$ | (23) |

where we have used the fact that

$$\delta \left({t}^{2}-{r}^{2}\right)={\left(2r\right)}^{-1}\left(\delta \left(\left|t\right|-r\right)+\delta \left(\left|t\right|+r\right)\right)$$

and that $\delta \left(\left|t\right|+r\right)=0$ since both $\left|t\right|$ and $r$ are necessarily positive.

Substitution of these results into (19) gives finally,

$$R=-\frac{1}{2}\sum _{n}\sum _{m}{\int}_{-\infty}^{+\infty}{\int}_{-\infty}^{+\infty}{e}_{n}{e}_{m}\left(1-{{x}^{.}}_{n}\left(t\right)\xb7{{x}^{.}}_{m}\left(s\right)\right)\times {\delta}_{+}\left({\left(t-s\right)}^{2}-{\left({x}_{n}\left(t\right)-{x}_{m}\left(s\right)\right)}^{2}\right)dtds.$$ | (24) |

The total complex action of the system is then^{10} ${S}_{p}+R$. Or, what amounts to the same thing; to obtain transition amplitudes including the effects of the field we must calculate the transition element of $\mathrm{exp}\left(iR\right)$:

$${\langle {\chi}_{{t}^{\u2033}}\left|\mathrm{exp}iR\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}$$ | (25) |

under the action ${S}_{p}$ of the particles, excluding interaction. Expression (24) for $R$ must be considered to be written in the usual manner as a Riemann sum and the expression (25) interpreted as defined in C [Eq. (39)]. Expression (6) must be used for ${x}_{n}^{.}$ at time $t$.

Expression (25), with (24), then contains all the effects of virtual quanta on a (at least non-relativistic) system according to quantum electrodynamics. It contains the effects to all orders in ${e}^{2}/\hslash c$ in a single expression. If expanded in a power series in ${e}^{2}/\hslash c$, the various terms give the expressions to the corresponding order obtained by the diagrams and methods of II. We illustrate this by an example in the next section.

We shall not be much concerned with the non-relativistic case here, as the relativistic case given below is as simple and more interesting. It is, however, very similar and at this stage it is worth giving an example to show how expressions resulting from (25) are to be interpreted according to the rules of C. For example, consider the case of a single electron, coordinate $x$, either free or in an external given potential (contained for simplicity in^{11} ${S}_{p}$, not in $R$). Its interaction with the field produces a reaction back on itself given by $R$ as in (24) but in which we keep only a single term corresponding to $m=n$. Assume the effect of $R$ to be small and expand $\mathrm{exp}\left(iR\right)$ as $1+iR$. Let us find the amplitude at time ${t}^{\u2033}$ of finding the electron in a state $\psi $ with no quanta emitted, if at time ${t}^{\prime}$ it was in the same state. It is

$${\langle {\psi}_{{t}^{\u2033}}\left|1+iR\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}={\langle {\psi}_{{t}^{\u2033}}\left|1\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}+i{\langle {\psi}_{{t}^{\u2033}}\left|R\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}$$

where ${\langle {\psi}_{{t}^{\u2033}}\left|1\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}=\mathrm{exp}\left[-iE\left({t}^{\u2033}-{t}^{\prime}\right)\right]$ if $E$ is the energy of the state, and

$${\langle {\psi}_{{t}^{\u2033}}\left|R\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}=-\frac{1}{2}{e}^{2}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}dt{\int}_{{t}^{\prime}}^{{t}^{\u2033}}ds{\langle {\psi}_{{t}^{\u2033}}\left|\left(1-{{x}^{.}}_{t}\xb7{{x}^{.}}_{s}\right){\delta}_{+}\left({\left(t-s\right)}^{2}-{\left({x}_{t}-{x}_{s}\right)}^{2}\right)\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}.$$ | (26) |

Here ${x}_{s}=x\left(s\right)$, etc. In (26) we shall limit the range of integrations by assuming $s<t$, and double the result.

The expression within the brackets ${\langle \text{}\rangle}_{{S}_{p}}$on the right-hand side of (26) can be evaluated by the methods described in C [Eq. (29)]. An expression such as (26) can also be evaluated directly in terms of the propagation kernel $K\left(2,1\right)$ [see I, Eq. (2)] for an electron moving in the given potential.

The term ${x}_{t}^{.}\xb7{x}_{s}^{.}$ in the non-relativistic case produces an interesting complication which does not have an analog for the relativistic case with the Dirac equation. We discuss it below, but for a moment consider, in further detail expression (26) but with the factor $\left(1-{x}_{s}^{.}\xb7{x}_{t}^{.}\right)$ replaced simply by unity.

The kernel $K\left(2,1\right)$ is defined and discussed in I. From its definition as the amplitude that the electron be found at ${x}_{2}$ at time ${t}_{2}$, if at ${t}_{1}$ it was at ${x}_{1}$, we have

$$K\left({x}_{2},{t}_{2};{x}_{1},{t}_{1}\right)={\langle \delta \left(x-{x}_{2}\right){t}_{2}\left|1\right|\delta \left(x-{x}_{1}\right){t}_{1}\rangle}_{{S}_{p}}$$ | (27) |

that is, more simply $K\left(2,1\right)$ is the sum of $\mathrm{exp}\left(i{S}_{p}\right)$ over all paths which go from space time point 1 to 2.

In the integrations over all paths implied by the symbol in (26) we can first integrate over all the ${x}_{i}$ variables corresponding to times ${t}_{i}$ from ${t}^{\prime}$ to $s$, not inclusive, the result being a factor $K\left({x}_{s},s;{x}_{{t}^{\prime}},{t}^{\prime}\right)$ according to (27). Next we integrate on the variables between $s$ and $t$ not inclusive, giving a factor $K\left({x}_{t},t;{x}_{s},s\right)$ and finally on those between $t$ and ${t}^{\u2033}$ giving $K\left({x}_{{t}^{\u2033}},{t}^{\u2033};{x}_{t},t\right)$. Hence the left-hand term in (26) excluding the ${x}_{t}^{.}\xb7{x}_{s}^{.}$ factor is

$$-{e}^{2}\int dt\int ds\int {\psi}^{*}\left({x}_{{t}^{\u2033}},{t}^{\u2033}\right)K\left({x}_{{t}^{\u2033}},{t}^{\u2033};{x}_{t},t\right){\delta}_{+}\left({\left(t-s\right)}^{2}-{\left({x}_{t}-{x}_{s}\right)}^{2}\right)\xb7K\left({x}_{t},t;{x}_{s},s\right)K\left({x}_{s},s;{x}_{{t}^{\prime}},{t}^{\prime}\right)$$ | (28) |

$$\times \psi \left({x}_{{t}^{\prime}},{t}^{\prime}\right){d}^{3}{x}_{{t}^{\u2033}}{d}^{3}{x}_{t}{d}^{3}{x}_{s}{d}^{3}{x}_{{t}^{\prime}}$$ | (28) |

which in improved notation and in the relativistic case is essentially the result given in II.

We have made use of a special case of a principle which may be stated more generally as

$${\langle {\chi}_{{t}^{\u2033}}\left|F\left({x}_{1},{t}_{1};{x}_{2},{t}_{2};\cdots {x}_{k},{t}_{k}\right)\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}=\int {\chi}^{*}\left({x}_{{t}^{\u2033}}\right)K\left({x}_{{t}^{\u2033}},{t}^{\u2033};{x}_{1},{t}_{1}\right)\xb7K\left({x}_{1},{t}_{1};{x}_{2},{t}_{2}\right)\cdots \text{}\times K\left({x}_{k-1},{t}_{k-1};{x}_{k},{t}_{k}\right)K\left({x}_{k},{t}_{k};{x}_{{t}^{\prime}},{t}^{\prime}\right)\xb7F\left({x}_{1},{t}_{1};{x}_{2},{t}_{2};\cdots {x}_{k},{t}_{k}\right)\text{}\times \psi \left({x}_{{t}^{\prime}}\right){d}^{3}{x}_{{t}^{\u2033}}{d}^{3}{x}_{1}{d}^{3}{x}_{2}\cdots {d}^{3}{x}_{k}{d}^{3}{x}_{{t}^{\prime}}$$ | (29) |

where $F$ is any function of the coordinate ${x}_{1}$, at time ${t}_{1}$, ${x}_{2}$ at ${t}_{2}$ up to ${x}_{k}$, ${t}_{k}$, and, it is important to notice, we have assumed ${t}^{\u2033}>{t}_{1}>{t}_{2}>\cdots {t}_{k}>{t}^{\prime}$.

Expressions of higher order arising for example from ${R}^{2}$ are more complicated as there are quantities referring to several different times mixed up, but they all can be interpreted readily. One simply breaks up the ranges of integrations of the time variables into parts such that in each the order of time of each variable is definite. One then interprets each part by formula (29).

As a simple example we may refer to the problem of element of the transition element

$$\langle {\chi}_{{t}^{\u2033}}\left|\int U\left(x\left(t\right),t\right)dt\int V\left(x\left(s\right),s\right)ds\right|{\psi}_{{t}^{\prime}}\rangle $$

arising, say, in the cross term in $U$ and $V$ in an ordinary second order perturbation problem (disregarding radiation) with perturbation potential $U\left(x,t\right)+V\left(x,t\right)$. In the integration on $s$ and $t$ which should include the entire range of time for each, we can split the range of $s$ into two parts, $s<t$ and $s>t$. In the first case, $s<t$, the potential $V$ acts earlier than $U$, and in the other range, vice versa, so that

$$\langle {\chi}_{{t}^{\u2033}}\left|\int U\left({x}_{t},t\right)dt\int V\left({x}_{s},s\right)ds\right|{\psi}_{{t}^{\prime}}\rangle ={\int}_{{t}^{\prime}}^{{t}^{\u2033}}dt{\int}_{{t}^{\prime}}^{t}ds\int {\chi}^{*}\left({x}_{{t}^{\u2033}}\right)K\left({x}_{{t}^{\u2033}},{t}^{\u2033};{x}_{t},t\right)U\left({x}_{t},t\right)K\left({x}_{t},t;{x}_{s},s\right)$$ | (30) |

$$\times V\left({x}_{s},s\right)K\left({x}_{s},s;{x}_{{t}^{\prime}},{t}^{\prime}\right)\psi \left({x}_{{t}^{\prime}}\right){d}^{3}{x}_{{t}^{\u2033}}{d}^{3}{x}_{t}{d}^{3}{x}_{s}{d}^{3}{x}_{{t}^{\prime}}+{\int}_{{t}^{\prime}}^{{t}^{\u2033}}dt{\int}_{t}^{{t}^{\u2033}}ds\int {\chi}^{*}\left({x}_{{t}^{\u2033}}\right)K\left({x}_{{t}^{\u2033}},{t}^{\u2033};{x}_{s},s\right)\text{}\times V\left({x}_{s},s\right)K\left({x}_{s},s;{x}_{t},t\right)U\left({x}_{t},t\right)K\left({x}_{t},t;{x}_{{t}^{\prime}},{t}^{\prime}\right)\psi \left({x}_{{t}^{\prime}}\right){d}^{3}{x}_{{t}^{\u2033}}{d}^{3}{x}_{s}{d}^{3}{x}_{t}{d}^{3}{x}_{{t}^{\prime}}$$ | (30) |

so that the single expression on the left is represented by two terms analogous to the two terms required in analyzing the Compton effect. It is in this way that the several terms and their corresponding diagrams corresponding to each process arise when an attempt is made to represent the transition elements of single expressions involving time integrals in terms of the propagation kernels $K$.

It remains to study in more detail the term in (26) arising from ${x}^{.}\left(t\right)\xb7{x}^{.}\left(s\right)$ in the interaction. The interpretation of such expressions is considered in detail in C, and we must refer to Eqs. (39) through (50) of that paper for a more thorough analysis. A similar type of term also arises in the Lagrangian formulation in simpler problems, for example the transition element

$$\langle {\chi}_{{t}^{\u2033}}\left|\int {x}^{.}\left(t\right)\xb7A\left(x\left(t\right),t\right)dt\int {x}^{.}\left(s\right)\xb7B\left(x\left(s\right),s\right)ds\right|{\psi}_{{t}^{\prime}}\rangle $$

arising say, in the cross term in A and B in a second-order perturbation problem for a particle in a perturbing vector potential $A\left(x,t\right)+B\left(x,t\right)$. The time integrals must first be written as Riemannian sums, the velocity (see (6)) being replaced by ${x}^{.}=\frac{1}{2}{\epsilon}^{-1}\left({x}_{i+1}-{x}_{i}\right)+\frac{1}{2}{\epsilon}^{-1}\left({x}_{i}-{x}_{i-1}\right)$ so that we ask for the transition element of

$$\sum _{i}\sum _{j}\left[\frac{1}{2}\left({x}_{i+1}-{x}_{i}\right)+\frac{1}{2}\left({x}_{i}-{x}_{i-1}\right)\right]\xb7A\left({x}_{i},{t}_{i}\right)\times \left[\frac{1}{2}\left({x}_{j+1}-{x}_{j}\right)+\frac{1}{2}\left({x}_{j}-{x}_{j-1}\right)\right]\xb7B\left({x}_{j},{t}_{j}\right)$$ | (31) |

In C it is shown that when converted to operator notation the quantity $\left({x}_{i-1}-{x}_{i}\right)/\epsilon $ is equivalent (nearly, see below) to an operator,

$$\left({x}_{i+1}-{x}_{i}\right)/\epsilon \u27f6i\left(H\text{}x-x\text{}H\right)$$ | (32) |

operating in order indicated by the time index $i$ (that is after ${x}_{l}$’s for $l\le i$ and before all ${x}_{l}$’s for $l>i$). In nonrelativistic mechanics $i\left(H\text{}x-x\text{}H\right)$ is the momentum operator ${p}_{x}$ divided by the mass $m$. Thus in (31) the expression $\frac{1}{2}\left({x}_{i+1}-{x}_{i}\right)+\frac{1}{2}\left({x}_{i}-{x}_{i-1}\right)\xb7A\left({x}_{i},{t}_{i}\right)$ becomes $\epsilon \left(p\xb7A+A\xb7p\right)/2m$. Here again we must split the sum into two regions $j<i$ and $j>i$ so the quantities in the usual notation will operate in the right order such that eventually (31) becomes identical with the right-hand side of Eq. (30) but with $U\left({x}_{t},t\right)$ replaced by the operator

$$\frac{1}{2m}\left(\frac{1}{i}\frac{\partial}{\partial {x}_{t}}\xb7A\left({x}_{t},t\right)+A\left({x}_{t},t\right)\xb7\frac{1}{i}\frac{\partial}{\partial {x}_{t}}\right)$$

standing in the same place, and with the operator

$$\frac{1}{2m}\left(\frac{1}{i}\frac{\partial}{\partial {x}_{s}}\xb7B\left({x}_{s},s\right)+\frac{1}{i}B\left({x}_{s},s\right)\xb7\frac{\partial}{\partial {x}_{s}}\right)$$

standing in the place of $V\left({x}_{s},s\right)$. The sums and factors ϵ have now become $\int dt\int ds$.

This is nearly but not quite correct, however, as there is an additional term coming from the terms in the sum corresponding to the special values, $j=i$, $j=i+1$ and $j=i-1$. We have tacitly assumed from the appearance of the expression (31) that, for a given $i$, the contribution from just three such special terms is of order ${\epsilon}^{2}$. But this is not true. Although the expected contribution of a term like $\left({x}_{i+1}-{x}_{i}\right)\left({x}_{j+1}-{x}_{j}\right)$ for $j\ne i$ is indeed of order ${\epsilon}^{2}$, the expected contribution of ${\left({x}_{i+1}-{x}_{i}\right)}^{2}$ is $+i\text{}\epsilon \text{}{m}^{-1}$[C, Eq. (50)], that is, of order ϵ. In nonrelativistic mechanics the velocities are unlimited and in very short times ϵ the amplitude diffuses a distance proportional to the square root of the time. Making use of this equation then we see that the additional contribution from these terms is essentially

$$i\text{}{m}^{-1}\epsilon \sum _{i}A\left({x}_{i},{t}_{i}\right)\xb7B\left({x}_{i},{t}_{i}\right)=i\text{}{m}^{-1}\int A\left(x\left(t\right),t\right)\xb7B\left(x\left(t\right),t\right)dt$$

when summed on all $i$. This has the same effect as a first-order perturbation due to a potential $A\xb7B/m$. Added to the term involving the momentum operators we therefore have an additional term^{12}

$$\frac{i}{m}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}dt\int {\chi}^{*}\left({x}_{{t}^{\u2033}}\right)K\left({x}_{{t}^{\u2033}},{t}^{\u2033};{x}_{t},t\right)A\left({x}_{t},t\right)\xb7B\left({x}_{t},t\right)\xb7K\left({x}_{t},t;{x}_{{t}^{\prime}},{t}^{\prime}\right)\psi \left({x}_{{t}^{\prime}}\right){d}^{3}{x}_{{t}^{\u2033}}{d}^{3}{x}_{t}{d}^{3}{x}_{{t}^{\prime}}.$$ | (33) |

In the usual Hamiltonian theory this term arises, of course, from the term ${A}^{2}/2m$ in the expansion of the Hamiltonian

$$H={\left(2m\right)}^{-1}{\left(p-A\right)}^{2}={\left(2m\right)}^{-1}\left({p}^{2}-p\xb7A-A\xb7p+{A}^{2}\right)$$

while the other term arises from the second-order action of $p\xb7A+A\xb7p$. We shall not be interested in non-relativistic quantum electrodynamics in detail. The situation is simpler for Dirac electrons. For particles satisfying the Klein-Gordon equation (discussed in Appendix A) the situation is very similar to a four-dimensional analog of the non-relativistic case given here.

Expressions (24) and (25) and their proof can be readily generalized to the relativistic case according to the one electron theory of Dirac. We shall discuss the hole theory later. In the non-relativistic case we began with the proposition that the amplitude for a particle to proceed from one point to another is the sum over paths of $\mathrm{exp}\left(i\text{}{S}_{p}\right)$, that is, we have for example for a transition element

$$\langle \chi \left|1\right|\psi \rangle =\underset{\epsilon \to 0}{\mathrm{lim}}\int \cdots \int {\chi}^{*}\left({x}_{N}\right){\Phi}_{p}\left({x}_{N},{x}_{N-1},\cdots {x}_{0}\right)\xb7\psi \left({x}_{0}\right){d}^{3}{x}_{0}{d}^{3}{x}_{1}\cdots {d}^{3}{x}_{N}$$ | (34) |

where for $\mathrm{exp}\left(i\text{}{S}_{p}\right)$ we have written ${\Phi}_{p}$, that is more precisely,

$${\Phi}_{p}={\Pi}_{i}{A}^{-1}\mathrm{exp}i\text{}S\left({x}_{i+1},{x}_{i}\right).$$

As discussed in C this form is related to the usual form of quantum mechanics through the observation that

$${\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}={A}^{-1}\mathrm{exp}\left[i\text{}S\left({x}_{i+1},{x}_{i}\right)\right]$$ | (35) |

where ${\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}$ is the transformation matrix from a representation in which $x$ is diagonal at time ${t}_{i}$ to one in which $x$ is diagonal at time ${t}_{i+1}={t}_{i}+\epsilon $ (so that it is identical to ${K}_{0}\left({x}_{i+1},{t}_{i+1};{x}_{i},{t}_{i}\right)$ for the small time interval $\epsilon $). Hence the amplitude for a given path can also be written

$${\Phi}_{p}={\Pi}_{i}{\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}$$ | (36) |

for which form, of course, (34) is exact irrespective of whether ${\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}$ can be expressed in the simple form (35).

For a Dirac electron the ${\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}$ is a $4\times 4$ matrix (or ${4}^{N}\times {4}^{N}$ if we deal with $N$ electrons) but the expression (34) with (36) is still correct (as it is in fact for any quantum-mechanical system with a sufficiently general definition of the coordinate ${x}_{i}$). The product (36) now involves operators, the order in which the factors are to be taken is the order in which the terms appear in time.

For a Dirac particle in a vector and scalar potential (times the electron charge $e$) $A\left(x,t\right)$, ${A}_{4}\left(x,t\right)$, the quantity ${\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}^{\left(A\right)}$ is related to that of a free particle to the first order in $\epsilon $ as

$${\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}^{\left(A\right)}={\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}^{\left(0\right)}\mathrm{exp}\left[-i\left(\epsilon {A}_{4}\left({x}_{i},{t}_{i}\right)-\left({x}_{i+1}-{x}_{i}\right)\xb7A\left({x}_{i},{t}_{i}\right)\right)\right].$$ | (37) |

This can be verified directly by substitution into the Dirac equation.^{13} It neglects the variation of $A$ and ${A}_{4}$ with time and space during the short interval $\epsilon $. This produces errors only of order ${\epsilon}^{2}$ in the Dirac case for the expected square velocity ${\left({x}_{i+1}-{x}_{i}\right)}^{2}/{\epsilon}^{2}$ during the interval ϵ is finite (equaling the square of the velocity of light) rather than being of order $1/\epsilon $ as in the non-relativistic case. [This makes the relativistic case somewhat simpler in that it is not necessary to define the velocity as carefully as in (6); $\left({x}_{i+1}-{x}_{i}\right)/\epsilon $ is sufficiently exact, and no term analogous to (33) arises.]

Thus ${\Phi}_{p}^{\left(A\right)}$ differs from that for a free particle, ${\Phi}_{p}^{\left(0\right)}$, by a factor ${\Pi}_{i}\mathrm{exp}-i\left(\epsilon {A}_{4}\left({x}_{i},{t}_{i}\right)-\left({x}_{i+1}-{x}_{i}\right)\xb7A\left({x}_{i},{t}_{i}\right)\right)$ which in the limit can be written as

$$\mathrm{exp}\left\{-i\int \left[{A}_{4}\left(x\left(t\right),t\right)-{x}^{.}\left(t\right)\xb7A\left(x\left(t\right),t\right)\right]dt\right\}$$ | (38) |

exactly as in the non-relativistic case.

The case of a Dirac particle interacting with the quantum-mechanical oscillators representing the field may now be studied. Since the dependence of ${\Phi}_{p}^{\left(A\right)}$ on $A$, ${A}_{4}$ is through the same factor as in the non-relativistic case, when $A$, ${A}_{4}$ are expressed in terms of the oscillator coordinates $q$, the dependence of $\Phi $ on the oscillator coordinates $q$ is unchanged. Hence the entire analysis of the preceding sections which concern the results of the integration over oscillator coordinates can be carried through unchanged and the results will be expression (25) with formula (24) for $R$. Expression (25) is now interpreted as

$$\langle {\chi}_{{t}^{\u2033}}\left|\mathrm{exp}i\text{}R\right|{\psi}_{{t}^{\prime}}\rangle =\underset{\epsilon \to 0}{\mathrm{lim}}\int {\chi}^{*}\left({x}_{{t}^{\u2033}}^{\left(1\right)},{x}_{{t}^{\u2033}}^{\left(2\right)}\cdots \right)\times \prod _{n}\left({\Phi}_{p,n}^{\left(0\right)}{d}^{3}{x}_{{t}^{\u2033}}^{\left(n\right)}{d}^{3}{x}_{{t}^{\u2033}-\epsilon}^{\left(n\right)}\cdots {d}^{3}{x}_{{t}^{\prime}}^{\left(n\right)}\right)\mathrm{exp}\left(i\text{}R\right)\psi \left({x}_{{t}^{\prime}}^{\left(1\right)},{x}_{{t}^{\prime}}^{\left(2\right)}\cdots \right)$$ | (39) |

where ${{\Phi}_{p,n}}^{\left(0\right)},$ the amplitude for a particular path for particle $n$ is simply the expression (36) where ${\left({x}_{i+1}|{x}_{i}\right)}_{e}$ is the kernel ${K}_{0,n}\left({{x}_{i+1}}^{\left(n\right)},{t}_{i+1};{{x}_{i}}^{\left(n\right)},{t}_{i}\right)$ for a free electron according to the one electron Dirac theory, with the matrices which appear operating on the spinor indices corresponding to particle ($n$) and the order of all operations being determined by the time indices.

For calculational purposes we can, as before, expand $R$ as a power series and evaluate the various terms in the same manner as for the non-relativistic case. In such an expansion the quantity ${x}^{.}$($t$) is replaced, as we have seen in (32), by the operator $i\left(H\text{}x-x\text{}H\right)$, that is, in this case by $\alpha $ operating at the corresponding time. There is no further complicated term analogous to (33) arising in this case, for the expected value of ${\left({x}_{i+1}-{x}_{i}\right)}^{2}$ is now of order ${\epsilon}^{2}$ rather than $\epsilon $.

For example, for self-energy one sees that expression (28) will be (with other terms coming from those with ${x}^{.}$($t$) replaced by $\alpha $ and with the usual $\beta $ in back of each ${K}_{0}$ because of the definition of ${K}_{0}$ in relativity theory)

$${\langle {\psi}_{{t}^{\u2033}}\left|R\right|{\psi}_{{t}^{\prime}}\rangle}_{{S}_{p}}=-{e}^{2}\int {\psi}^{*}\left({x}_{{t}^{\u2033}}\right){K}_{0}\left({x}_{{t}^{\u2033}},{t}^{\u2033};{x}_{t},t\right){\beta a}_{\mu}\xb7{\delta}_{+}\left({\left(t-s\right)}^{2}\u2013{\left({x}_{t}\u2013{x}_{s}\right)}^{2}\right){K}_{0}\left({x}_{t},t;{x}_{s},s\right){\mathrm{\beta \alpha}}_{\mu}$$ | (40) |

$$\xb7{K}_{0}\left({x}_{s},s;{x}_{{t}^{\prime}},{t}^{\prime}\right)\mathrm{\beta \psi}\left({x}_{{t}^{\prime}}\right){d}^{3}{x}_{{t}^{\u2033}}{d}^{3}{x}_{t}{d}^{3}{x}_{s}{d}^{3}{x}_{{t}^{\prime}}dtds,$$ | (40) |

where ${a}_{4}=1,{a}_{1,2,3}={a}_{x,y,z}$ and a sum on the repeated index μ is implied in the usual way; ${a}_{\mu}{b}_{\mu}={a}_{4}{b}_{4}-{a}_{1}{b}_{1}-{a}_{2}{b}_{2}-{a}_{3}{b}_{3}$. One can change ${\beta a}_{\mu}$ to ${\gamma}_{\mu}$ and ${\psi}^{*}$ to $\stackrel{\_}{\psi}\beta .$ In this manner all of the rules referring to virtual photons discussed in **II** are deduced; but with the difference that ${K}_{0}$ is used instead of ${K}_{+}$ and we have the Dirac one electron theory with negative energy states (although we may have any number of such electrons).

Since in (39) we have an arbitrary number of electrons, we can deal with the hole theory in the usual manner by imagining that we have an infinite number of electrons in negative energy states.

On the other hand, in paper **I** on the theory of positrons, it was shown that the results of the hole theory in a system with a given external potential ${A}_{\mu}$ were equivalent to those of the Dirac one electron theory if one replaced the propagation kernel, ${K}_{0}$, by a different one, ${K}_{+}$, and multiplied the resultant amplitude by factor $C$_{ν} involving ${A}_{\mu}$. We must now see how this relation, derived in the case of external potentials, can also be carried over in electrodynamics to be useful in simplifying expressions involving the infinite sea of electrons.

To do this we study in greater detail the relation between a problem involving virtual photons and one involving purely external potentials. In using (25) we shall assume in accordance with the hole theory that the number of electrons is infinite, but that they all have the same charge, $e$. Let the states ${\psi}_{{t}^{\prime}},{\chi}_{{t}^{\u2033}},$ represent the vacuum plus perhaps a number of real electrons in positive energy states and perhaps also some empty negative energy states. Let us call the amplitude for the transition in an external potential ${B}_{\mu}$, but excluding virtual photons, ${T}_{0}\left[B\right]$, a functional of ${B}_{\mu}$(1). We have seen (38)

$${T}_{0}\left[B\right]=\langle {\chi}_{{t}^{\u2033}}\left|\mathrm{exp}\right|{\psi}_{{t}^{\prime}}\rangle $$ | (41) |

where

$$P=-\sum _{n}\int \left[{B}_{4}\left({x}^{\left(n\right)}\left(t\right),t\right)\u2013{x}^{.\left(n\right)}\left(t\right)\xb7B\left({x}^{\left(n\right)}\left(t\right),t\right)\right]dt$$

by (38). We can write this as

$$P=-\sum _{n}\int {B}_{\mu}\left({{x}_{v}}^{\left(n\right)}\left(t\right)\right){{\stackrel{.}{x}}_{\mu}}^{\left(n\right)}\left(t\right)dt$$

where ${x}_{4}\left(t\right)=t$ and ${\stackrel{.}{x}}_{4}=1$, the other values of μ corresponding to space variables. The corresponding amplitude for the same process in the same potential, but *including* all the virtual photons we may call,

$${T}_{{e}^{2}}\left[B\right]=\langle {\chi}_{{t}^{\u2033}}\left|\mathrm{exp}\left(\mathrm{iR}\right)\mathrm{exp}\left(\mathrm{iP}\right)\right|{\psi}_{{t}^{\prime}}\rangle .$$ | (42) |

Now let us consider the effect on ${T}_{{e}^{2}}\left[B\right]$ of changing the coupling ${e}^{2}$ of the virtual photons. Differentiating (42) with respect to ${e}^{2}$ which appears only^{14} in $R$ we find

$$d{T}_{{e}^{2}}\left[B\right]/d\left({e}^{2}\right)=\langle {\chi}_{{t}^{\u2033}}\left|\u2014\frac{i}{2}\sum _{n}\sum _{m}\int \int dtds{{\stackrel{.}{x}}_{\mu}}^{\left(n\right)}\left(t\right){{\stackrel{.}{x}}_{\mu}}^{\left(m\right)}\left(s\right)\xb7{\delta}_{+}\left({\left({{x}_{v}}^{\left(n\right)}\left(t\right)-{{x}_{v}}^{\left(m\right)}\left(s\right)\right)}^{2}\right)\mathrm{exp}\left(R+P\right)\right|{\psi}_{{t}^{\prime}}\rangle .$$ | (43) |

We can also study the first-order effect of a change of ${B}_{\mu}$:

$$d{T}_{{e}^{2}}\left[B\right]/\delta {B}_{\mu}\left(1\right)=-i\langle {\chi}_{{t}^{\u2033}}\left|\sum _{n}\int dt{\stackrel{.}{x}}_{\mu}^{\left(n\right)}{\delta}^{4}\left({x}_{\alpha}^{\left(n\right)}\left(t\right)-{x}_{\alpha ,1}\right)\xb7\mathrm{exp}i\left(R+P\right)\right|{\psi}_{{t}^{\prime}}\rangle $$ | (44) |

where ${x}_{a,1}$ is the field point at which the derivative with respect to ${B}_{\mu}$, is taken^{15} and the term (current density) $\u2013{\sum}_{n}\int d\text{}t{{\stackrel{.}{x}}_{\mu}}^{\left(n\right)}\left(t\right){\delta}^{4}\left({{x}_{\alpha}}^{\left(n\right)}\left(t\right)\u2013{x}_{a,1}\right)$ is just $\delta \text{}P/\delta \text{}{B}_{\mu}\left(1\right)$. The function ${\delta}^{4}\left({{x}_{\alpha}}^{\left(n\right)}{\mathrm{\u2013x}}_{\alpha ,1}\right)$ means $\delta \left({{x}_{4}}^{\left(n\right)}\u2013{x}_{4,1}\right)\xb7\delta \left({x}_{3}^{\left(n\right)}-{x}_{3,1}\right)\xb7\delta \left({x}_{2}^{\left(n\right)}-{x}_{2,1}\right)\xb7\delta \left({x}_{1}^{\left(n\right)}-{x}_{1,1}\right)$ that is, $\delta \left(2,1\right)$ with ${x}_{\alpha ,2}={x}_{\alpha}^{\left(n\right)}\left(t\right)$. A second variation of $T$ gives, by differentiation of (44) with respect to ${B}_{\nu}\left(2\right)$,

$${\delta}^{2}{T}_{{e}^{2}}\left[B\right]/\delta {B}_{\mu}\left(1\right)\delta {B}_{\nu}\left(2\right)=-\langle {\chi}_{{t}^{\u2033}}\left|\sum _{n}\sum _{m}\int dt\int ds{\stackrel{.}{x}}_{\mu}^{\left(n\right)}\left(t\right){\stackrel{.}{x}}_{\nu}^{\left(m\right)}\left(s\right)\text{}\xb7{\delta}^{4}\left({x}_{\alpha}^{\left(n\right)}\left(t\right)-{x}_{\alpha ,1}\right){\delta}^{4}\left({x}_{\beta}^{\left(n\right)}\left(s\right)-{x}_{\beta ,2}\right)\mathrm{exp}i\left(R+P\right)\right|{\psi}_{{t}^{\prime}}\rangle .$$

Comparison of this with (43) shows that

$$d{T}_{{e}^{2}}\left[B\right]/d\left({e}^{2}\right)=\frac{1}{2}i\int \int \left({\delta}^{2}{T}_{{e}^{2}}\left[B\right]/\delta {B}_{\mu}\left(1\right)\delta {B}_{\mu}\left(2\right)\right)\times {\delta}_{+}\left({s}_{12}^{2}\right)d{\tau}_{1}d{\tau}_{2}$$ | (45) |

where ${s}_{12}^{2}=\left({x}_{\mu ,1}-{x}_{\mu ,2}\right)\left({x}_{\mu ,1}-{x}_{\mu ,2}\right)$.

We now proceed to use this equation to prove the validity of the rules given in II for electrodynamics. This we do by the following argument. The equation can be looked upon as a differential equation for ${T}_{{e}^{2}}\left[B\right]$. It determines ${T}_{{e}^{2}}\left[B\right]$ uniquely if ${T}_{o}\left[B\right]$ is known. We have shown it is valid for the hole theory of positrons. But in I we have given formulas for calculating ${T}_{o}\left[B\right]$ whose correctness relative to the hole theory we have there demonstrated. Hence we have shown that the ${T}_{{e}^{2}}\left[B\right]$ obtained by solving (45) with the initial condition ${T}_{o}\left[B\right]$ as given by the rules in I will be equal to that given for the same problem by the second quantization theory of the Dirac matter field coupled with the quantized electromagnetic field. But it is evident (the argument is given in the next paragraph) that the rules^{16} given in II constitute a solution in power series in ${e}^{2}$ of the Eq. (45) [which for ${e}^{2}=0$ reduce to the ${T}_{0}\left[B\right]$ given in I]. Hence the rules in II must give, to each order in ${e}^{2}$, the matrix element for any process that would be calculated by the usual theory of second quantization of the matter and electromagnetic fields. This is what we aimed to prove.

That the rules of II represent, in a power series expansion, a solution of (45) is clear. For the rules there given may be stated as follows: Suppose that we have a process to order $k$ in ${e}^{2}$ (i.e., having $k$ virtual photons) and order $n$ in the external potential ${B}_{\mu}$. Then, *the matrix element for the process with one more virtual photon and two less potentials is that obtained from the previous matrix by choosing from the *$n$* potentials a pair, say *${B}_{\mu}\left(1\right)$* acting at 1 and* ${B}_{\nu}\left(2\right)$* acting at 2, replacing them by *$i{e}^{2}{\delta}_{\mathrm{\mu \nu}}{\delta}_{+}\left({s}_{12}^{2}\right)$*, adding the results for each way of choosing the pair, and dividing by *$k+1$*, the present number of photons*. The matrix with no virtual photons $\left(k=0\right)$ being given to any $n$ by the rules of I, this permits terms to all orders in ${e}^{2}$ to be derived by recursion. It is evident that the rule in italics is that of II, and equally evident that it is a word expression of Eq. (45). [The factor $\frac{1}{2}$ in (45) arises since in integrating over all $d{\tau}_{1}$and $d{\tau}_{2}$ we count each pair twice.The division by $k+1$ is required by the rules of II for, there, each diagram is to be taken only once, while in the rule given above we say what to do to add one extra virtual photon to $k$ others. But which one of the $k+1$ is to be identified at the last photon added is irrelevant. It agrees with (45) of course for it is canceled on differentiating with respect to ${e}^{2}$ the factor ${\left({e}^{2}\right)}^{k+1}$ for the $\left(k+1\right)$ photons.]

The relation implied by (45) between the formal solution for the amplitude for a process in an arbitrary unquantized external potential to that in a quantized field appears to be of much wider generality. We shall discuss the relation from a more general point of view here (still limiting ourselves to the case of no photons in initial or final state).

In earlier sections we pointed out that as a consequence of the Lagrangian form of quantum mechanics the aspects of the particles’ motions and the behavior of the field could be analyzed separately. What we did was to integrate over the field oscillator coordinates first. We could, in principle, have integrated over the particle variables first. That is, we first solve the problem with the action of the particles and their interaction with the field and then multiply by the exponential of the action of the field and integrate over all the field oscillator coordinates. (For simplicity of discussion let us put aside from detailed special consideration the questions involving the separation of the longitudinal and transverse parts of the field.) Now the integral over the particle coordinates for a given process is precisely the integral required for the analysis of the motion of the particles in an unquantized potential. With this observation we may suggest a generalization to all types of systems.

Let us suppose the formal solution for the amplitude for some given process with matter in an external potential ${B}_{\mu}\left(1\right)$ is some numerical quantity ${T}_{0}$. We mean matter in a more general sense now, for the motion of the matter may be described by the Dirac equation, or by the Klein-Gordon equation, or may involve charged or neutral particles other than electrons and positrons in any manner whatsoever. The quantity ${T}_{0}$ depends of course on the potential function ${B}_{\mu}\left(1\right)$; that is, it is a functional ${T}_{0}\left[B\right]$ of this potential. We assume we have some expression for it in terms of ${B}_{\mu}$ (exact, or to some desired degree of approximation in the strength of the potential).

Then the answer ${T}_{{e}^{2}}\left[B\right]$ to the corresponding problem in quantum electrodynamics is ${T}_{0}\left[{A}_{\mu}\left(1\right)+{B}_{\mu}\left(1\right)\right]\mathrm{exp}\left(i{S}_{0}\right)$ summed over all possible distributions of field ${A}_{\mu}\left(1\right)$, wherein ${S}_{0}$ is the action for the field ${S}_{0}=-{\left(8\pi {e}^{2}\right)}^{-1}{\sum}_{\mu}\int \left({\left(\partial {A}_{\mu}/\partial t\right)}^{2}-{\left(\nabla {A}_{\mu}\right)}^{2}\right){d}^{3}xdt$ the sum on $\mu $ carrying the usual minus sign for space components.

If $F\left[A\right]$ is any functional of ${A}_{\mu}\left(1\right)$ we shall represent by ${\u200a}_{0}{\left|F\left[A\right]\right|}_{0}$ this superposition of $F\left[A\right]\mathrm{exp}\left(i{S}_{0}\right)$ over distributions of ${A}_{\mu}$ for the case in which there are no photons in initial or final state. That is, we have

$${T}_{{e}^{2}}\left[B\right]={\u200a}_{0}{\left|{T}_{0}\left[A+B\right]\right|}_{0}.$$ | (46) |

The evaluation of ${\u200a}_{0}{\left|F\left[A\right]\right|}_{0}$ directly from the definition of the operation ${\u200a}_{0}{\left|\right|}_{0}$ is not necessary. We can give the result in another way. We first note that the operation is linear,

$${\u200a}_{0}{\left|{F}_{1}\left[A\right]+{F}_{2}\left[A\right]\right|}_{0}={\u200a}_{0}{\left|{F}_{1}\left[A\right]\right|}_{0}+{\u200a}_{0}{\left|{F}_{2}\left[A\right]\right|}_{0}$$ | (47) |

so that if $F$ is represented as a sum of terms each term can be analyzed separately. We have studied essentially the case in which $F\left[A\right]$ is an exponential function. In fact, what we have done in Section 4 may be repeated with slight modification to show that

$${\u200a}_{0}{\left|\mathrm{exp}\left(-i\int {j}_{\mu}\left(1\right){A}_{\mu}\left(1\right)d{\tau}_{1}\right)\right|}_{0}=\mathrm{exp}\left(-\frac{1}{2}i\text{}{e}^{2}\int \int {j}_{\mu}\left(1\right){j}_{\mu}\left(2\right){\delta}_{+}\left({s}_{12}^{2}\right)d{\tau}_{1}d{\tau}_{2}\right)$$ | (48) |

where ${j}_{\mu}\left(1\right)$ is an arbitrary function of position and time for each value of $\mu $.

Although this gives the evaluation of ${\u200a}_{0}{\left|\text{}\right|}_{0}$ for only a particular functional of ${A}_{\mu}$ the appearance of the arbitrary function ${j}_{\mu}\left(1\right)$ makes it sufficiently general to permit the evaluation for any other functional. For it is to be expected that any functional can be represented as a superposition of exponentials with different functions ${j}_{\mu}\left(1\right)$ (by analogy with the principle of Fourier integrals for ordinary functions). Then, by (47), the result of the operation is the corresponding superposition of expressions equal to the right-hand side of (48) with the various $j$’s substituted for ${j}_{\mu}$.

In many applications $F\left[A\right]$ can be given as a power series in ${A}_{\mu}$:

$$F\left[A\right]={f}_{0}+\int {f}_{\mu}\left(1\right){A}_{\mu}\left(1\right)d{\tau}_{1}+\int \int {f}_{\mu \nu}\left(1,2\right){A}_{\mu}\left(1\right){A}_{\nu}\left(2\right)d{\tau}_{1}d{\tau}_{2}+\cdots $$ | (49) |

where ${f}_{0},{f}_{\mu}\left(1\right),{f}_{\mu \nu}\left(1,2\right)\cdots $ are known numerical functions independent of ${A}_{\mu}$. Then by (47)

$${\u200a}_{0}{\left|F\left[A\right]\right|}_{0}={f}_{0}+\int {f}_{\mu}\left(1\right){\u200a}_{0}{\left|{A}_{\mu}\left(1\right)\right|}_{0}d{\tau}_{1}+\int \int {f}_{\mu \nu}\left(1,2\right){\u200a}_{0}{\left|{A}_{\mu}\left(1\right){A}_{\nu}\left(2\right)\right|}_{0}d{\tau}_{1}d{\tau}_{2}+\cdots $$ | (50) |

where we set ${\u200a}_{0}{\left|1\right|}_{0}=1$ (from (48) with ${j}_{\mu}=0$). We can work out expressions for the successive powers of ${A}_{\mu}$ by differentiating both sides of (48) successively with respect to ${j}_{\mu}$ and setting ${j}_{\mu}=0$ in each derivative. For example, the first variation (derivative) of (48) with respect to ${j}_{\mu}\left(3\right)$ gives

$${\u200a}_{0}{\left|-i{A}_{\mu}\left(3\right)\mathrm{exp}\left(-i\int {j}_{\nu}\left(1\right){A}_{\nu}\left(1\right)d{\tau}_{1}\right)\right|}_{0}=-i\text{}{e}^{2}\int {\delta}_{+}\left({s}_{34}^{2}\right){j}_{\mu}\left(4\right)d{\tau}_{4}\phantom{\rule{0ex}{0ex}}\times \mathrm{exp}\left(-\frac{1}{2}i\text{}{e}^{2}\int \int {j}_{\nu}\left(1\right){j}_{\nu}\left(2\right){\delta}_{+}\left({s}_{12}^{2}\right)d{\tau}_{1}d{\tau}_{2}\right).$$ | (51) |

Setting ${j}_{\mu}=0$ gives

$${\u200a}_{0}{\left|{A}_{\mu}\left(3\right)\right|}_{0}=0.$$

Differentiating (5l) again with respect to ${j}_{\nu}\left(4\right)$ and setting ${j}_{\nu}=0$ shows

$${\u200a}_{0}{\left|{A}_{\mu}\left(3\right){A}_{\nu}\left(4\right)\right|}_{0}=i\text{}{e}^{2}{\delta}_{\mu \nu}{\delta}_{+}\left({s}_{34}^{2}\right)$$ | (52) |

and so on for higher powers. These results may be substituted into (50). Clearly therefore when ${T}_{0}\left[B+A\right]$ in (46) is expanded in a power series and the successive terms are computed in this way, we obtain the results given in **II**.

It is evident that (46), (47), (48) imply that ${T}_{{e}^{2}}\left[B\right]$ satisfies the differential equation (45) and conversely (45) with the definition (46) implies (47) and (48). For if ${T}_{0}\left[B\right]$ is an exponential

$${T}_{0}\left[B\right]=\mathrm{exp}\left(-i\int {j}_{\mu}\left(1\right){B}_{\mu}\left(1\right)d{\tau}_{1}\right)$$ | (53) |

$${T}_{{e}^{2}}\left[B\right]=\mathrm{exp}\left[-\frac{1}{2}i\text{}{e}^{2}\int \int {j}_{\mu}\left(1\right){j}_{\mu}\left(2\right){\delta}_{+}\left({s}_{12}^{2}\right)d{\tau}_{1}d{\tau}_{2}\right]\xb7\mathrm{exp}\left[-i\int {j}_{\nu}\left(1\right){B}_{\nu}\left(1\right)d{\tau}_{1}\right].$$ | (54) |

Direct substitution of this into Eq. (45) shows it to be a solution satisfying the boundary condition (53). Since the differential equation (45) is linear, if ${T}_{0}\left[B\right]$ is a superposition of exponentials, the corresponding superposition of solutions (54) is also a solution.

Many of the formal representations of the matter system (such as that of second quantization of Dirac electrons) represent the interaction with a fixed potential in a formal exponential form such as the left-[hand] side of (48), except that ${j}_{\mu}\left(1\right)$ is an operator instead of a numerical function, Equation (48) may still be used if care is exercised in defining the order of the operators on the right-hand side. The succeeding paper will discuss this in more detail.

Equation (45) or its solution (46), (47), (48) constitutes a very general and convenient formulation of the laws of quantum electrodynamics for virtual processes. Its relativistic invariance is evident if it is assumed that the unquantized theory giving ${T}_{0}\left[B\right]$ is invariant. It has been proved to be equivalent to the usual formulation for Dirac electrons and positrons (for Klein-Gordon particles see Appendix A). It is suggested that it is of wide generality. It is expressed in a form which has meaning even if it is impossible to express the matter system in Hamiltonian form; in fact, it only requires the existence of an amplitude for fixed potentials which obeys the principle of superposition of amplitudes. If ${T}_{0}\left[B\right]$ is known in power series in $B$, calculations of ${{T}_{e}}^{2}\left[B\right]$ in a power series of ${e}^{2}$ can be made directly using the italicized rule of Sec. 7. The limitation to virtual quanta is removed in the next section.

On the other hand, the formulation is unsatisfactory because for situations of importance it gives divergent results, even if ${T}_{0}\left[B\right]$ is finite. The modification proposed in II of replacing ${\delta}_{+}\left({{s}_{12}}^{2}\right)$ in (45), (48) by ${f}_{+}\left({{s}_{12}}^{2}\right)$ is not satisfactory owing to the loss of the theorems of conservation of energy or probability discussed in II at the end of Sec. 6. There is the additional difficulty in positron theory that even ${T}_{0}\left[B\right]$ is infinite to begin with (vacuum polarization). Computational ways of avoiding these troubles are given in II and in the references of footnote 2.

The case in which there are real photons in the initial or the final state can be worked out from the beginning in the same manner.^{17} We first consider the case of a system interacting with a single oscillator. From this result the generalization will be evident. This time we shall calculate the transition element between an initial state in which the particle is in state ${\psi}_{{t}^{\prime}}$ and the oscillator is in its $n$th eigenstate (i.e., there are $n$ photons in the field) to a final state with particle in ${\chi}_{{t}^{\u2033}}$, oscillator in $m$th level. As we have already discussed, when the coordinates of the oscillator are eliminated the result is the transition element $\langle {\chi}_{{t}^{\u2033}}\left|{G}_{mn}\right|{\psi}_{{t}^{\prime}}\rangle $ where

$${G}_{m\text{}n}=\int {{\phi}_{m}}^{*}\left({q}_{j}\right)k\left({q}_{j},{t}^{\u2033};{q}_{0},{t}^{\prime}\right){\phi}_{n}\left({q}_{0}\right)d{q}_{0}d{q}_{j}$$

where ${\phi}_{m}$, ${\phi}_{n}$$$ are the wave functions [8] for the oscillator in state $m$, $n$ and $k$ is given in (12). The ${G}_{mn}$ can be evaluated most easily by calculating the generating function

$$g\left(X,Y\right)=\sum _{m}\sum _{n}{G}_{mn}{X}^{m}{{Y}^{n}\left(m!n!\right)}^{-\frac{1}{2}}$$ | (55) |

for arbitrary $X$, $Y$. If expression (11) is substituted in the left-hand side of (55), the expression can be simplified by use of the generating function relation for the eigenfunctions [8] of the harmonic oscillator

$$\sum _{n}{\phi}_{n}\left({q}_{0}\right){Y}^{n}{\left(n!\right)}^{-\frac{1}{2}}={\left(\omega /\pi \right)}^{\frac{1}{2}}\mathrm{exp}\left(-\frac{1}{2}i\text{}\omega \text{}{t}^{\prime}\right)\times \mathrm{exp}\frac{1}{2}\left[\omega {{q}_{0}}^{2}-{\left(Y\mathrm{exp}\left[-i\text{}\omega \text{}{t}^{\prime}\right]-{\left(2\omega \right)}^{\frac{1}{2}}{q}_{0}\right)}^{2}\right]$$

Using a similar expansion for the ${{\phi}_{m}}^{*}$ one is left with the exponential of a quadratic function of ${q}_{0}$ and ${q}_{j}$. The integration on ${q}_{0}$ and ${q}_{j}$ is then easily performed to give

$$g\left(X,Y\right)={G}_{00}\mathrm{exp}\left(XY+i\text{}{\beta}^{*}X+i\text{}\beta \text{}Y\right)$$ | (56) |

from which expansion in powers of $X$ and $Y$ and comparison to (3.11) gives the final result

$${G}_{mn}={G}_{00}{\left(m!n!\right)}^{-\frac{1}{2}}\sum _{r}\frac{m!}{\left(m-r\right)!r!}\frac{n!}{\left(n-r\right)!r!}\times r!{\left(i\text{}{\beta}^{*}\right)}^{m-r}{\left(i\text{}\beta \right)}^{n-r}$$ | (57) |

where ${G}_{00}$ is given in (14) and

$$\begin{array}{cc}\beta & ={\left(2\omega \right)}^{-\frac{1}{2}}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}\gamma \left(t\right)\mathrm{exp}\left(-i\text{}\omega \text{}t\right)dt,\\ {\beta}^{*}& ={\left(2\omega \right)}^{-\frac{1}{2}}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}\gamma \left(t\right)\mathrm{exp}\left(+i\text{}\omega \text{}t\right)dt,\end{array}$$ | (58) |

and the sum on $r$ is to go from 0 to $m$ or to $n$ whichever is the smaller. (The sum can be expressed as a Laguerre polynomial but there is no advantage in this.)

Formula (57) is readily understandable. Consider first a simple case of absorption of one photon. Initially we have one photon and finally none. The amplitude for this is the transition element of ${G}_{01}=i\text{}\beta \text{}{G}_{00}$ or $\langle {\chi}_{{t}^{\u2033}}\left|i\text{}\beta \text{}{G}_{00}\right|{\psi}_{{t}^{\prime}}\rangle $. This is the same as would result if we asked for the transition element for a problem in which all photons are virtual but there was present a perturbing potential $-{\left(2\omega \right)}^{-\frac{1}{2}}\gamma \left(t\right)\mathrm{exp}\left(i\text{}\omega \text{}t\right)$ and we required the first-order effect of this potential. Hence photon absorption is like the first order action of a potential varying in time as $\gamma \left(t\right)\mathrm{exp}\left(i\text{}\omega \text{}t\right)$ that is with a positive frequency (i.e., the sign of the coefficient of $t$ in the exponential corresponds to positive energy). The amplitude for emission of one photon involves ${G}_{10}=i\text{}{\beta}^{*}{G}_{00}$, which is the same result except that the potential has negative frequency. Thus we begin by interpreting $i\text{}{\beta}^{*}$ as the amplitude for emission of one photon $i\text{}\beta $ as the amplitude for absorption of one.

Next for the general case of $n$ photons initially and finally we may understand (57) as follows. We first neglect Bose statistics and imagine the photons as individual distinct particles. If we start with $n$ and end with $m$ this process may occur in several different ways. The particle may absorb in total $n-r$ of the photons and the final $m$ photons will represent $r$ of the photons which were present originally plus $m-r$ new photons emitted by the particle. In this case the $n-r$ which are to be absorbed may be chosen from among the original $n$ in $n!/\left(n-r\right)!r!$ different ways, and each contributes a factor $i\text{}\beta $, the amplitude for absorption of a photon. Which of the $m-r$ photons from among the $m$ arc emitted can be chosen in $m!/\left(m-r\right)!r!$ different ways and each photon contributes a factor $i{\beta}^{*}$ in amplitude. The initial $r$ photons which do not interact with the particle can be re-arranged among the final $r$ in $r!$ ways. We must sum over the alternatives corresponding to different values of $r$. Thus the form of ${G}_{mn}$ can be understood. The remaining factor ${\left(m!\right)}^{-\frac{1}{2}}{\left(n!\right)}^{-\frac{1}{2}}$ may be interpreted as saying that in computing probabilities (which therefore involves the square of ${G}_{mn}$) the photons may be considered as independent but that if $m$ are actually equal the statistical weight of each of the states which can be made by rearranging the m equal photons is only $1/m!$. This is the content of Bose statistics; that $m$ equal particles in a given state represents just one state, i.e., has statistical weight unity, rather than the $m!$ statistical weight which would result if it is imagined that the particles and states can be identified and rearranged in $m!$ different ways. This holds for both the initial and final states of course. From this rule about the statistical weights of states the derivation of the blackbody distribution law follows.

The actual electromagnetic field is represented as a host of oscillators each of which behaves independently and produces its own factor such as ${G}_{mn}$. Initial or final states may also be linear combinations of states in which one or another oscillator is excited. The results for this case are of course the corresponding linear combination of transition elements.

For photons of a given direction of polarization and for sin or cos waves the explicit expression for β can be obtained directly from (58) by substituting the formulas (16) for the γ’s for the corresponding oscillator. It is more convenient to use the linear combination corresponding to running waves. Thus we find the amplitude for absorption of a photon of momentum **K**, frequency $k={\left(K\xb7K\right)}^{\frac{1}{2}}$ polarized in direction **e** is given by including a factor $i$ times

$${\beta}_{K,e}={\left(4\pi \right)}^{\frac{1}{2}}{\left(2k\right)}^{-\frac{1}{2}}\sum _{n}{e}_{n}{\int}_{{t}^{\prime}}^{{t}^{\u2033}}\mathrm{exp}\left(-i\text{}k\text{}t\right)\times \mathrm{exp}\left(iK\xb7{x}_{n}\left(t\right)\right)e\xb7{x}_{n}^{.}\left(t\right)dt$$ | (59) |

in the transition element (25). The density of states in momentum space is now ${\left(2\pi \right)}^{-3}{d}^{3}K$. The amplitude for emission is just $i$times the complex conjugate of this expression, or what amounts to the same thing, the same expression with the sign of the four vector $k,\text{'}$ (?) reversed. Since the factor (59) is exactly the first-order effect of a vector potential

$${A}^{PH}={\left(2\pi /k\right)}^{\frac{1}{2}}e\mathrm{exp}\left(-i\left(k\text{}t-K\xb7x\right)\right)$$

of the corresponding classical wave, we have derived the rules for handling real photons discussed in II.

We can express this directly in terms of the quantity ${T}_{{e}^{2}}\left[B\right]$, the amplitude for a given transition without emission of a photon. What we have said is that the amplitude for absorption of just one photon whose classical wave form is ${{A}_{\mu}}^{PH}\left(1\right)$ (time variation $\mathrm{exp}\left(-i\text{}k\text{}{t}_{1}\right)$ corresponding to positive energy $k$) is proportional to the first order (in ϵ) change produced in ${T}_{{e}^{2}}\left[B\right]$ on changing $B$ to $B+\epsilon \text{}{A}^{PH}$. That is, more exactly,

$$\int \left({\mathrm{\delta T}}_{{e}^{2}}\left[B\right]/\delta {B}_{\mu}\left(1\right)\right){{A}_{\mu}}^{PH}\left(1\right)d{\tau}_{1}$$ | (60) |

is the amplitude for absorption by the particle system of one photon, ${A}^{PH}$. (A superposition argument shows the expression to be valid not only for plane waves, but for spherical waves, etc., as given by the form of ${A}^{PH}$.) The amplitude for emission is the same expression but with the sign of the frequency reversed in ${A}^{PH}$. The amplitude that the system absorbs two photons with waves ${A}_{\mu}^{P{H}_{1}}$ and ${A}_{\nu}^{P{H}_{2}}$ is obtained from the next derivative,

$$\int \int \left({\delta}^{2}{T}_{{e}^{2}}\left[B\right]/\delta {B}_{\mu}\left(1\right)\delta {B}_{\nu}\left(2\right)\right){{A}_{\mu}}^{P{H}_{1}}\left(1\right){{A}_{\nu}}^{P{H}_{2}}\left(2\right)d{\tau}_{1}d{\tau}_{2}$$

the same expression holding for the absorption of one and emission of the other, or emission of both depending on the sign of the time dependence of ${A}^{P{H}_{1}}$ and ${A}^{P{H}_{2}}$. Larger photon numbers correspond to higher derivatives, absorption of ${l}_{1}$ emission of ${l}_{2}$ requiring the $\left({l}_{1}+{l}_{2}\right)$ derivaties. When two or more of the photons are exactly the same (e.g., ${A}^{P{H}_{1}}={A}^{P{H}_{2}}$) the same expression holds for the amplitude that ${l}_{1}$ are absorbed by the system while ${l}_{2}$ are emitted. However, the statement that initially $n$ of a kind are present and $m$ of this kind are present finally, does not imply ${l}_{1}=n$ and ${l}_{2}=m$. It is possible that only $n-r={l}_{1}$ were absorbed by the system and $m-r={l}_{2}$ emitted, and that $r$ remained from initial to final state without interaction. This term is weighed by the combinatorial coefficient ${\left(m!n!\right)}^{-\frac{1}{2}}\left(\begin{array}{c}m\\ r\end{array}\right)\left(\begin{array}{c}n\\ r\end{array}\right)r!$ and summed over the possibilities for $r$ as explained in connection with (57). Thus once the amplitude for virtual processes is known, that for real photon processes can be obtained by differentiation.

It is possible, of course, to deal with situations in which the electromagnetic field is not in a definite state after the interaction. For example, we might ask for the total probability of a given process, such as scattering, without regard for the number of photons emitted. This is done of course by squaring the amplitude for the emission of $m$ photons of a given kind and summing on all $m$. Actually the sums and integrations over the oscillator momenta can usually easily be performed analytically. For example, the amplitude, starting from vacuum and ending with $m$ photons of a given kind, is by (56) just

$${G}_{m0}={\left(m!\right)}^{-\frac{1}{2}}{{G}_{00}\left(i\text{}{\beta}^{*}\right)}^{m}.$$ | (61) |

The square of the amplitude summed on $m$ requires the product of two such expressions (the $\gamma \left(t\right)$ in the β of one and in the other will have to be kept separately) summed on $m$:

$${\sum}_{m}{{G}_{m0}}^{*}{{G}_{m0}}^{\prime}={\sum}_{m}{{G}_{00}}^{*}{{{G}_{00}}^{\prime}\left(m!\right)}^{-1}{{\beta}^{m}\left({\beta}^{\prime *}\right)}^{m}={{G}_{00}}^{*}{{G}_{00}}^{\prime}\mathrm{exp}\left({\mathrm{\beta \beta}}^{\prime *}\right)$$

In the resulting expression the sum over all oscillators is easily done. Such expressions can be of use in the analysis in a direct manner of problems of line width, of the Bloch-Nordsieck infra-red problem, and of statistical mechanical problems, but no such applications will be made here.

The author appreciates his opportunities to discuss these matters with Professor H. A. Bethe and Professor J. Ashkin, and the help of Mr. M. Baranger with the manuscript.

In this Appendix we describe a formulation of the equations for a particle of spin zero which was first used to obtain the rules given in II for such particles. The complete physical significance of the equations has not been analyzed thoroughly so that it may be preferable to derive the rules directly from the second quantization formulation of Pauli and Weisskopf. This can be done in a manner analogous to the derivation of the rules for the Dirac equation given in I or from the Schwinger-Tomonaga formulation [2] in a manner described, for example, by Rohrlich.^{18} The formulation given here is therefore not necessary for a description of spin zero particles but is given only for its own interest as an alternative to the formulation of second quantization.

We start with the Klein-Gordon equation

$${\left(i\partial /\partial {x}_{\mu}-{A}_{\mu}\right)}^{2}\psi ={m}^{2}\psi $$ | (A1) |

for the wave function ψ of a particle of mass $m$ in a given external potential ${A}_{\mu}$. We shall try to represent this in a manner analogous to the formulation of quantum mechanics in C. That is, we try to represent the amplitude for a particle to get from one point to another as a sum over all trajectories of an amplitude $\mathrm{exp}\left(i\text{}S\right)$ where $S$ is the classical action for a given trajectory. To maintain the relativistic invariance in evidence the idea suggests itself of describing a trajectory in space-time by giving thc four variables ${x}_{\mu}\left(u\right)$ as functions of some fifth parameter $u$ (rather than expressing ${x}_{1}$, ${x}_{2}$, ${x}_{3}$ in terms of ${x}_{4}$). As we expect to represent paths which may reverse themselves in time (to represent pair production, etc., as in I) this is certainly a more convenient representation, for all four functions ${x}_{\mu}\left(u\right)$ may be considered as functions of a parameter $u$ (somewhat analogous to proper time) which increase as we go along the trajectory, whether the trajectory is proceeding forward $\left(d{x}_{4}/du>0\right)$ or backward $\left(d{x}_{4}/du<0\right)$ in time.^{19} We shall then have a new type of wave function $\phi \left(x,u\right)$ a function of five variables, $x$ standing for the four ${x}_{\mu}$. It gives the amplitude for arrival at point ${x}_{\mu}$ with a certain value of the parameter $u$. We shall suppose that this wave function satisfies the equation

$$i\partial \phi /\partial u=-\frac{1}{2}{\left(i\partial /\partial {x}_{\mu}-{A}_{\mu}\right)}^{2}\phi $$ | (A2) |

which is seen to be analogous to the time-dependent Schrödinger equation, $u$ replacing the time and the four coordinates of space-time ${x}_{\mu}$ replacing the usual three coordinates of space.

Since the potentials ${A}_{\mu}\left(x\right)$ are functions only of coordinates ${x}_{\mu}$, and are independent of $u$, the equation is separable in $u$ and we can write a special solution in the form $\phi =\mathrm{exp}\left(\frac{1}{2}i\text{}{m}^{2}u\right)\psi \left(x\right)$ where $\psi \left(x\right)$, a function of the coordinates ${x}_{\mu}$ only, satisfies (A1) and the eigenvalue $\frac{1}{2}{m}^{2}$ conjugate to $u$ is related to the mass $m$ of the particle. Equation (A2) is therefore equivalent to the Klein-Gordon Eq. (A1) provided we ask in the end only for the solution of (A1) corresponding to the eigenvalue $\frac{1}{2}{m}^{2}$ for the quantity conjugate to $u$.

We may now proceed to represent Eq. (A2) in Lagrangian form in general and without regard to this eigenvalue condition. Only in the final solutions need we apply the eigenvalue condition. That is, if we have some special solution $\phi \left(x,u\right)$ of (A2) we can select that part corresponding to the eigenvalue $\frac{1}{2}{m}^{2}$ by calculating

$$\psi \left(x\right)={\int}_{-\infty}^{\infty}\mathrm{exp}\left(-\frac{1}{2}i\text{}{m}^{2}u\right)\phi \left(x,u\right)du$$ | (A3) |

and thereby obtain a solution ψ of Eq. (A1).

Since (A2) is so closely analogous to the Schrödinger equation, it is easily written in the Lagrangian form described in C, simply by working by analogy. For example if $\phi \left(x,u\right)$ is known at one value of $u$ its value at a slightly larger $u+\epsilon $ is given by

$$\phi \left(x,u+\epsilon \right)=\int \mathrm{exp}i\text{}\epsilon \left[-\frac{{\left({x}_{\mu}-{{x}_{\mu}}^{\prime}\right)}^{2}}{2{\epsilon}^{2}}-\frac{1}{2}\left(\frac{\left({x}_{\mu}-{{x}_{\mu}}^{\prime}\right)}{\epsilon}\right)\left({A}_{\mu}\left(x\right)+{A}_{\mu}\left({x}^{\prime}\right)\right)\right]\xb7\phi \left({x}^{\prime},u\right){d}^{4}{\tau}_{{x}^{\prime}}2{\left(\pi \text{}i\text{}\epsilon \right)}^{\frac{3}{2}}{\left(-2\pi \text{}i\text{}\epsilon \right)}^{-\frac{1}{2}}$$ | (A4) |

where ${\left({x}_{\mu}-{{x}_{\mu}}^{\prime}\right)}^{2}$ means $\left({x}_{\mu}-{{x}_{\mu}}^{\prime}\right)$$\left({x}_{\mu}-{{x}_{\mu}}^{\prime}\right)$, ${d}^{4}{\tau}_{{x}^{\prime}}=d{{x}_{1}}^{\prime}d{{x}_{2}}^{\prime}d{{x}_{3}}^{\prime}d{{x}_{4}}^{\prime}$ and the sign of the normalizing factor is changed for the ${x}_{4}$ component since the component has the reversed sign in its quadratic coefficient in the exponential, in accordance with our summation convention ${a}_{\mu}{b}_{\mu}={a}_{4}{b}_{4}-{a}_{1}{b}_{1}-{a}_{2}{b}_{2}-{a}_{3}{b}_{3}$. Equation (A4), as can be verified readily as described in C, Sec. 6, is equivalent to first order in ϵ to Eq. (A2). Hence, by repeated use of this equation the wave function at ${u}_{0}={n}_{\epsilon}$ can be represented in terms of that at $u=0$ by:

$$\phi \left({x}_{\nu ,n},{u}_{0}\right)=\int \mathrm{exp}-\frac{i\text{}\epsilon}{2}\sum _{i=1}^{n}\left[{\left(\frac{{x}_{\mu ,i}-{x}_{\mu ,i-1}}{\epsilon}\right)}^{2}+{\epsilon}^{-1}\left({x}_{\mu ,i}-{x}_{\mu ,i-1}\right)\left({A}_{\mu}\left({x}_{i}\right)+{A}_{\mu}\left({x}_{i-1}\right)\right)\right]$$ | (A5) |

$$\xb7\phi \left({x}_{\nu ,0},0\right)\prod _{i=0}^{n-1}\left({d}^{4}{\tau}_{i}/4{\pi}^{2}{\epsilon}^{2}i\right).$$ | (A5) |

That is, roughly, the amplitude for getting from one point to another with a given value of ${u}_{0}$ is the sum over all trajectories of $\mathrm{exp}\left(i\text{}S\right)$ where

$$S=-\underset{0}{\overset{{u}_{0}}{\int}}\left[\frac{1}{2}{\left(d{x}_{\mu}/du\right)}^{2}+\left(d{x}_{\mu}/du\right){A}_{\mu}\left(x\right)\right]du,$$ | (A6) |

when sufficient care is taken to define the quantities, as in C. This completes the formulation for particles in a fixed potential but a few words of description may be in order.

In the first place in the special case of a free particle we can define a kernel ${k}^{\left(0\right)}\left(x,{u}_{0};{x}^{\prime},0\right)$ for arrival from ${{x}_{\mu}}^{\prime},0$ to ${x}_{\mu}$ at ${u}_{0}$ as the sum over all trajectories between these points of $\mathrm{exp}-i{\int}_{0}^{{u}_{0}}\frac{1}{2}{\left(d{x}_{\mu}/du\right)}^{2}du$. Then for this case we have

$$\phi \left(x,{u}_{0}\right)=\int {k}^{\left(0\right)}\left(x,{u}_{0};{x}^{\prime},0\right)\phi \left({x}^{\prime},0\right){d}^{4}{\tau}_{{x}^{\prime}}$$ | (A7) |

and it is easily verified that ${k}_{0}$ is given by

$${k}^{\left(0\right)}\left(x,{u}_{0};{x}^{\prime},0\right)={\left(4{\pi}^{2}{{u}_{0}}^{2}i\right)}^{-1}\mathrm{exp}-{i\left({x}_{\mu}-{{x}_{\mu}}^{\prime}\right)}^{2}/2{u}_{0}$$ | (A8) |

for ${u}_{0}>0$ and by $0$, by definition, for ${u}_{0}<0$. The corresponding kernel of importance when we select the eigenvalue $\frac{1}{2}{m}^{2}$ is^{20}

$$\begin{array}{cc}2i\text{}{I}_{+}\left(x,{x}^{\prime}\right)& ={\int}_{\mathrm{\u2013\infty}}^{\infty}{k}^{\left(0\right)}\left(x,{u}_{0};{x}^{\prime},0\right)\mathrm{exp}\left(-\frac{1}{2}i\text{}{m}^{2}{u}_{0}\right)d{u}_{0}\\ & ={\int}_{0}^{\infty}d{{u}_{0}\left(4{\pi}^{2}{{u}_{0}}^{2}i\right)}^{-1}\mathrm{exp}-\frac{1}{2}i\left({m}^{2}{u}_{0}+{{{u}_{0}}^{-1}\left({x}_{\mu}-{x}_{{\mu}^{\prime}}\right)}^{2}\right)\end{array}$$ | (A9) |

(the last extends only from ${u}_{0}=0$ since ${k}_{0}$ is zero for negative ${u}_{o}$) which is identical to the ${I}_{+}$ defined in II.^{21} This may be seen readily by studying the Fourier transform, for the transform of the integrand on the right-hand side is

$$\int {\left(4{\pi}^{2}{u}_{0}^{2}i\right)}^{-1}\mathrm{exp}\left(i\text{}p\xb7x\right)\mathrm{exp}-\frac{1}{2}i\left({m}^{2}{u}_{0}+{{x}_{\mu}}^{2}/{u}_{0}\right){d}^{4}{\tau}_{x}=\mathrm{exp}-\frac{1}{2}i{u}_{0}\left({m}^{2}-{{p}_{\mu}}^{2}\right)$$

so that the ${u}_{0}$ integration gives for the transform of ${I}_{+}$ just $1/\left({{p}_{\mu}}^{2}-{m}^{2}\right)$ with the pole defined exactly as in II. [21] Thus we are automatically representing the positrons as trajectories with the time sense reversed.

If ${\Phi}^{\left(0\right)}\left[x\left(\mu \right)\right]=\mathrm{exp}-i{\int}_{0}^{{u}_{0}}\frac{1}{2}{\left(d{x}_{\mu}/du\right)}^{2}du$ is the amplitude for a given trajectory ${x}_{\nu}\left(u\right)$ for a free particle, then the amplitude in a potential is

$${\Phi}^{\left(A\right)}\left[x\left(\mu \right)\right]={\Phi}^{\left(0\right)}\left[x\left(\mu \right)\right]\mathrm{exp}-i{\int}_{0}^{{u}_{0}}\left(d{x}_{\mu}/du\right){A}_{\mu}\left(x\right)du.$$ | (A10) |

If desired this may be studied by perturbation methods by expanding the exponential in powers of ${A}_{\mu}$.

For interpretation, the integral in (A10) must be written as a Riemann sum, and if a perturbation expansion is made, care must be taken with the terms quadratic in the velocity, for the effect of $\left({x}_{\mu ,i+1}-{x}_{\mu ,i}\right)\left({x}_{\nu ,i+1}-{x}_{\nu ,i}\right)$ is not of order ${\epsilon}^{2}$ but is $i\text{}{\delta}_{\mu \nu}\epsilon $. The “velocity” $d{x}_{\mu}/du$ becomes the momentum operator ${p}_{\mu}=+i\partial /\partial {x}_{\mu}$ operating half before and half after ${A}_{\mu}$, just as in the non-relativistic Schrödinger equation discussed in Sec. 5. Furthermore, in exactly the same manner as in that case, but here in four dimensions, a term quadratic in ${A}_{\mu}$ arises in the second-order perturbation terms from the coincidence of two velocities for the same value of $u$.

As an example, the kernel ${k}^{\left(A\right)}\left(x,{u}_{0};{x}^{\prime},0\right)$ for proceeding from ${x}_{\mu}^{\prime},0$ to ${x}_{\mu},{u}_{0}$ in a potential ${A}_{\mu}$ differs from ${k}^{\left(0\right)}$ to first order in ${A}_{\mu}$ by a term

$$-i{\int}_{0}^{{u}_{0}}d\text{}u{k}^{\left(0\right)}\left(x,{u}_{0};y,u\right)\frac{1}{2}\left({p}_{\mu}{A}_{\mu}\left(y\right)+{A}_{\mu}\left(y\right){p}_{\mu}\right){k}^{\left(0\right)}\left(y,{u}_{0};{x}^{\prime},0\right)d{\tau}_{y}$$ | (A11) |

the ${p}_{\mu}$ here meaning $+i\partial /\partial {y}_{\mu}$. The kernel of importance on selecting the eigenvalue $\frac{1}{2}{m}^{2}$ is obtained by multiplying this by $\mathrm{exp}\left(-\frac{1}{2}i\text{}{m}^{2}{u}_{0}\right)$ and integrating ${u}_{0}$ from 0 to $\infty $. The kernel ${k}^{\left(0\right)}\left(x,{u}_{0};y,u\right)$ depends only on ${u}^{\prime}={u}_{0}-u$ and in the integrals on $u$ and ${u}_{0}$; ${\int}_{0}^{\infty}d{u}_{0}{\int}_{0}^{\infty}du\mathrm{exp}\left(-\frac{1}{2}i{m}^{2}{u}_{0}\right)\dots $, can be written, on interchanging the order of integration and changing variables to $u$ and ${u}^{\prime}$, ${\int}_{0}^{\infty}du{\int}_{0}^{\infty}d{u}^{\prime}\mathrm{exp}\left(-\frac{1}{2}i\text{}{m}^{2}\left(u+{u}^{\prime}\right)\right)\dots $. Now the integral on ${u}^{\prime}$ converts ${k}^{\left(0\right)}\left(x,{u}_{0};y,u\right)$ to $2i\text{}{I}_{+}\left(x,y\right)$ by (A9), while that on $u$ converts ${k}^{\left(0\right)}\left(y,u;{x}^{\prime},0\right)$ to $2i\text{}{I}_{+}\left(y,{x}^{\prime}\right)$, so the result becomes

$$\int 2i{I}_{+}\left(x,y\right)\left({p}_{\mu}{A}_{\mu}+{A}_{\mu}{p}_{\mu}\right){I}_{+}\left(y,{x}^{\prime}\right){d}^{4}{\tau}_{y}$$ | (A12) |

as expected. The same principle works to any order so that the rules for a single Klein-Gordon particle in external potentials given in II, Section 9, are deduced.

The transition to quantum electrodynamics is simple for in (A6) we already have a transition amplitude represented as a sum (over trajectories, and eventually ${u}_{0}$) of terms, in each of which the potential appears in exponential form. We may make use of the general relation (54). Hence, for example, one finds for the case of no photons in the initial and final states, in the presence of an external potential ${B}_{\mu}$, the amplitude that a particle proceeds from $\left({x}_{\mu}^{\prime},0\right)$ to $\left({x}_{\mu},{u}_{0}\right)$ is the sum over all trajectories of the quantity

$$\mathrm{exp}-i\left[\frac{1}{2}{\int}_{0}^{{u}_{0}}{\left(\frac{d{x}_{\mu}}{du}\right)}^{2}du+{\int}_{0}^{{u}_{0}}\frac{d{x}_{\mu}}{du}{B}_{\mu}\left(x\left(u\right)\right)du+\frac{{e}^{2}}{2}{\int}_{0}^{{u}_{0}}{\int}_{0}^{{u}_{0}}\frac{d{x}_{\mu}\left(u\right)}{du}\frac{d{x}_{\nu}\left({u}^{\prime}\right)}{d{u}^{\prime}}\phantom{\rule{0ex}{0ex}}\times {\delta}_{+}\left({\left({x}_{\mu}\left(u\right)-{x}_{\mu}\left({u}^{\prime}\right)\right)}^{2}\right)dud{u}^{\prime}\right].$$ | (A13) |

This result must be multiplied by $\mathrm{exp}\left(-\frac{1}{2}i\text{}{m}^{2}{u}_{0}\right)$ and integrated on ${u}_{0}$ from zero to infinity to express the action of a Klein-Gordon particle acting on itself through virtual photons. The integrals are interpreted as Riemann sums, and if perturbation expansions arc made, the necessary care is taken with the terms quadratic in velocity. When there are several particles (other than the virtual pairs already included) one use a separate $u$ for each, and writes the amplitude for each set of trajectories as the exponential of $-i$ times

$$\frac{1}{2}\sum _{n}{\int}_{0}^{{{u}_{0}}^{\left(n\right)}}{\left(\frac{d{{x}_{\mu}}^{n)}}{du}\right)}^{2}du+\sum _{n}{\int}_{0}^{{{u}_{0}}^{\left(n\right)}}\frac{d{{x}_{\mu}}^{\left(n\right)}}{du}{B}_{\mu}\left({x}^{\left(n\right)}\left(u\right)\right)du$$ | (A14) |

$$+\frac{{e}^{2}}{2}\sum _{n}\sum _{m}{\int}_{0}^{{{u}_{0}}^{\left(n\right)}}{\int}_{0}^{{{u}_{0}}^{\left(m\right)}}\frac{d{{x}_{\nu}}^{\left(n\right)}\left(u\right)}{du}\frac{d{{x}_{\nu}}^{\left(m\right)}\left({u}^{\prime}\right)}{d{u}^{\prime}}$$ | (A14) |

$$\times {\delta}_{+}\left({\left({x}_{\mu}^{\left(n\right)}\left(u\right)-{x}_{\mu}^{\left(m\right)}\left({u}^{\prime}\right)\right)}^{2}\right)dud{u}^{\prime}$$ | (A14) |

where ${{x}_{\mu}}^{\left(n\right)}\left(u\right)$ are the coordinates of the trajectory of the nth particle.^{22} The solution should depend on the ${{u}_{0}}^{\left(n\right)}$ as $\mathrm{exp}\left(-\frac{1}{2}i{m}^{2}\sum _{n}{{u}_{0}}^{\left(n\right)}\right)$.

Actually, knowledge of the motion of a single charge implies a great deal about the behavior of several charges. For a pair which eventually may turn out to be a virtual pair may appear in the short run as two “other particles.” As a virtual pair, that is, as the reverse section of a very long and complicated single track we know its behavior by (A13). We can assume that such a section can be looked at equally well, for a limited duration at least, as being due to other unconnected particles. This then implies a definite law of interaction of particles if the self-action (A13) of a single particle is known. (This is similar to the relation of real and virtual photon processes discussed in detail in Appendix B.) It is possible that a detailed analysis of this could show that (A13) implied that (A14) was correct for many particles. There is even reason to believe that the law of Bose-Einstein statistics and the expression for contributions from closed loops could be deduced by following this argument. This has not yet been analyzed completely, however, so we must leave this formulation in an incomplete form. The expression for closed loops should come out to be ${C}_{\nu}=\mathrm{exp}+L$ where $L$, the contribution from a single loop, is

$$L=2{\int}_{0}^{\infty}l\left({u}_{0}\right)\mathrm{exp}\left(-\frac{1}{2}i{m}^{2}{u}_{0}\right)d{u}_{0}/{u}_{0}$$ | (A15) |

where $l\left({u}_{0}\right)$ is the sum over all trajectories which close on themselves $\left({x}_{\mu}\left({u}_{0}\right)={x}_{\mu}\left(0\right)\right)$ of $\mathrm{exp}\left(iS\right)$ with $S$ given in (A6), and a final integration $d{\tau}_{x\left(0\right)}$ on ${x}_{\mu}\left(0\right)$ is made. This is equivalent to putting

$$l\left({u}_{0}\right)=\int \left({k}^{\left(A\right)}\left(x,{u}_{0};x,0\right)-{k}^{\left(0\right)}\left(x,{u}_{0};x,0\right)\right)d{\tau}_{x}.$$ | (A16) |

The term ${k}^{\left(0\right)}$ is subtracted only to simplify convergence problems (as adding a constant independent of ${A}_{\mu}$, to $L$ has no effect).

If one has a general formula for all virtual processes he should be able to find the formulas and states involved in real processes. That is to say, we should be able to deduce the formulas of Section 9 directly from the formulation (24), (25) (or its generalized equivalent such as (46), (48)) without having to go all the way back to the more usual formulation. We discuss this problem here.

That this possibility exists can be seen from the consideration that what looks like a real process from one point of view may appear as a virtual process occurring over a more extended time.

For example, if we wish to study a given real process, such as the scattering of light, we can, if we wish, include in principle the source, scatterer, and eventual absorber of the scattered light in our analysis. We may imagine that no photon is present initially, and that the source then emits light (the energy coming say from kinetic energy in the source). The light is then scattered and eventually absorbed (becoming kinetic energy in the absorber). From this point of view the process is virtual; that is, we start with no photons and end with none. Thus we can analyze the process by means of our formula for virtual processes, and obtain the formulas for real processes by attempting to break the analysis into parts corresponding to emission, scattering, and absorption.^{23}

To put the problem in a more general way, consider the amplitude for some transition from a state empty of photons far in the past (time ${t}^{\prime}$) to a similar one far in the future ($t={t}^{\u2033}$). Suppose the time interval to be split into three regions $a$, $b$, $c$ in some convenient manner, so that region $b$ is an interval ${t}_{2}>t>{t}_{1}$ around the present time that we wish to study. Region $a$, (${t}_{1}>t>{t}^{\prime}$), precedes $b$, and $c$, (${t}^{\u2033}>t>{t}_{2}$), follows $b$. We want to see how it comes about that the phenomena during $b$ can be analyzed by a study of transitions ${g}_{ji}\left(b\right)$ between some initial state $i$ at time ${t}_{1}$, (which no longer need be photon-free), to some other final state $j$ at time ${t}_{2}$. The states $i$ and $j$ are members of a large class which we will have to find out how to specify. (The single index $i$ is used to represent a large number of quantum numbers, so that different values of $i$ will correspond to having various numbers of various kinds of photons in the field, etc.) Our problem is to represent the over-all transition amplitude, $g\left(a,b,c\right)$, as a sum over various values of $i$, $j$ of a product of three amplitudes,

$$g\left(a,b,c\right)={\sum}_{i}{\sum}_{j}{g}_{0j}\left(c\right){g}_{ji}\left(b\right){g}_{i0}\left(a\right);$$ | (B1) |

first the amplitude that during the interval $a$ the vacuum state makes transition to some state $i$, then the amplitude that during $b$ the transition to $j$ is made, and finally in $c$ the amplitude that the transition from $j$ to some photon-free state 0 is completed. The mathematical problem of splitting $g\left(a,b,c\right)$ is made definite by the further condition that ${g}_{ji}\left(b\right)$ for given $i$, $j$ must not involve the coordinates of the particles for times corresponding to regions $a$ or $c$, ${g}_{i0}\left(a\right)$ must involve those only in region $a$, and ${g}_{0j}\left(c\right)$ only in $c$.

To become acquainted with what is involved, suppose first that we do not have a problem involving virtual photons, but just the transition of a one-dimensional Schrödinger particle going in a long time interval from, say, the origin $o$ to the origin $o$, and ask what states $i$ we shall need for intermediary time intervals. We must solve the problem (B1) where $g\left(a,b,c\right)$ is the sum over all trajectories going from $o$ at ${t}^{\prime}$ to $o$ at ${t}^{\u2033}$ of $\mathrm{exp}iS$ where $S=\int Ldt$. The integral may be split into three parts $S={S}_{a}+{S}_{b}+{S}_{c}$ corresponding to the three ranges of time. Then $\mathrm{exp}\left(iS\right)=\mathrm{exp}\left(i{S}_{a}\right)\mathrm{exp}\left(i{S}_{b}\right)\mathrm{exp}\left(i{S}_{c}\right)$ and the separation (B1) is accomplished by taking for ${g}_{i0}\left(a\right)$ the sum over all trajectories lying in $a$ from $o$ to some end point ${x}_{{t}_{1}}$ of $\mathrm{exp}\left(i{S}_{a}\right)$, for ${g}_{ji}\left(b\right)$ the sum over trajectories in $b$ of $\mathrm{exp}\left(i{S}_{b}\right)$ between end points ${x}_{{t}_{1}}$ and ${x}_{{t}_{2}}$ and for ${g}_{0j}\left(c\right)$ the sum of $\mathrm{exp}\left(i{S}_{c}\right)$ over the section of the trajectory lying in $c$ and going from ${x}_{{t}_{2}}$ to $o$. Then the sum on $i$ and $j$ can be taken to be the integrals on ${x}_{{t}_{1}}$, ${x}_{{t}_{2}}$, respectively. Hence the various states $i$ can be taken to correspond to particles being at various coordinates $x$. (Of course any other representation of the states in the sense of Dirac’s transformation theory could be used equally well. Which one, whether coordinate, momentum, or energy level representation, ;is of course just a matter of convenience and wee cannot determine that simply from (B1).)

We can consider next the problem including virtual photons. That is, $g\left(a,b,c\right)$ now contains an additional factor $\mathrm{exp}\left(iR\right)$ where $R$ involves a double integral $\int \int $ over all time. Those parts of the index $i$ which correspond to the particle states can be taken in the same way as though $R$ were absent. We study now the extra complexities in the states produced by splitting the $R$. Let us first (solely for simplicity of the argument) take the case that there are only two regions $a$, $c$ separated by time ${t}_{0}$, and try to expand

$$g\left(a,c\right)={\sum}_{i}{g}_{0i}\left(c\right){g}_{i0}\left(a\right).$$ | (B2) |

The factor $\mathrm{exp}\left(i\text{}R\right)$ involves $R$ as a double integral which can be split into three parts ${\int}_{a}{\int}_{a}+{\int}_{c}{\int}_{c}+{\int}_{a}{\int}_{c}$ for the first of which both $t$, $s$ are in $a$, for the second both are in $c$, for the third one is in $a$ the other in $c$. Writing $\mathrm{exp}\left(i\text{}R\right)$ as $\mathrm{exp}\left(i\text{}{R}_{cc}\right)\mathrm{exp}\left(i\text{}{R}_{aa}\right)\mathrm{exp}\left(i\text{}{R}_{ac}\right)$ shows that the factors ${R}_{cc}$ and ${R}_{aa}$ produce no new problems for they can be taken bodily into ${g}_{0i}\left(c\right)$ and ${g}_{i0}\left(a\right)$ respectively. However, we must disentangle the variables which are mixed up in $\mathrm{exp}\left(i\text{}{R}_{ac}\right)$.

The expression for ${R}_{ac}$, is just twice (24) but with the integral on $s$ extending over the range $a$ and that for $t$ extending over $c$. Thus $\mathrm{exp}\left(i\text{}{R}_{ac}\right)$ contains the variables for times in $a$ and in $c$ in a quite complicated mixture. Our problem is to write $\mathrm{exp}\left(i\text{}{R}_{ac}\right)$ as a sum over possibly a vast class of states $i$ of the product of two parts, like${{h}_{i}}^{\prime}\left(c\right){h}_{i}\left(a\right)$, each of which involves the coordinates in one interval alone.

This separation may be made in many different ways, corresponding to various possible representations of the state of the electromagnetic field. We choose a particular one. First we can expand the exponential, $\mathrm{exp}\left(i\text{}{R}_{ac}\right)$, in a power series, as $\sum _{n}{{i}^{n}\left(n!\right)}^{-1}{\left({R}_{ac}\right)}^{n}$. The states $i$ can therefore be subdivided into subclasses corresponding to an integer $n$ which we can interpret as the number of quanta in the field at time ${t}_{0}$. The amplitude for the case $n=0$ clearly just involves $\mathrm{exp}\left(i\text{}{R}_{aa}\right)$ and $\mathrm{exp}\left(i\text{}{R}_{cc}\right)$ in the way that it should if we interpret these as the amplitudes for regions $a$ and $c$, respectively, of making a transition between a state of zero photons and another state of zero photons.

Next consider the case $n=1$. This implies an additional factor in the transitional element; the factor ${R}_{ac}$. The variables are still mixed up. But an easy way to perform the separation suggests itself. Namely, expand the ${\delta}_{+}\left({\left(t-s\right)}^{2}-{\left({x}_{n}\left(t\right)-{x}_{m}\left(s\right)\right)}^{2}\right)$ in ${R}_{ac}$, as a Fourier integral as

$$i\int \mathrm{exp}\left(-i\text{}k\left|t-s\right|\right)\mathrm{exp}\left(-i\xb7\left({x}_{n}\left(t\right)-{x}_{m}\left(s\right)\right)\right){d}^{3}K/4{\pi}^{2}k.$$ | (B3) |

For the exponential can be written immediately as a product of $\mathrm{exp}+i\left(K\xb7{x}_{m}\left(s\right)\right)$, a function only of coordinates for times $s$ in $a$ (suppose $s<t$), and $\mathrm{exp}-i\left(K\xb7{x}_{m}\left(s\right)\right)$ (a function only of coordinates during interval $c$). The integral on ${d}^{3}K$ can be symbolized as sum over states $i$ characterized by the value of K. Thus the state with $n=1$ must be further characterized by specifying a vector **K**, interpreted as the momentum of the photon. Finally the factor $\left(1-{x}_{n}^{.}\left(t\right)\xb7{x}_{m}^{.}\left(s\right)\right)$ in ${R}_{ac}$ is simply the sum of four parts each of which is already split (namely 1, and each of the three components in the vector scalar product). Hence each photon of momentum K must still be characterized by specifying it as one of four varieties; that is, there are four polarizations.^{24} Thus in trying to represent the effect of the past $a$ on the future $c$ we are lead to invent photons of four polarizations and characterized by a propagation vector K.

The term for a given polarization and value of K (for $n=1$) is clearly just $-{\beta}_{c}{{\beta}_{a}}^{*}$ where the ${\beta}_{a}$ is defined in (59) but with the time integral extending just over region $a$, while ${\beta}_{c}$ is the same expression with the integration over region $c$. Hence the amplitude for transition during interval a from a state with no quanta to a state with one in a given state of polarization and momentum is calculated by inclusion of an extra factor $i\text{}{{\beta}_{a}}^{*}$ in the transition element. Absorption in region $c$ corresponds to a factor $i\text{}{\beta}_{c}$.

We next turn to the case $n=2$. This requires analysis of ${R}_{a{c}^{2}}$. The ${\delta}_{+}$ can be expanded again as a Fourier integral, but for each of the two ${\delta}_{+}$ in $\frac{1}{2}{R}_{a{c}^{2}}$ we have a value of **K** which may be different. Thus we say we have two photons, one of momentum **K** and one momentum ${K}^{\prime}$ and we sum over all values of $K$ and ${K}^{\prime}$. (Similarly each photon is characterized by its own independent polarization index.) The factor $\frac{1}{2}$ can be taken into account neatly by asserting that we count each possible pair of photons as constituting just one state at time ${t}_{0}$. Then the $\frac{1}{2}$ arises for the sum over all $K$, ${K}^{\prime}$ (and polarizations) counts each pair twice. On the other hand, for the terms representing two identical photons ($K={K}^{\prime}$) of like polarization, the $\frac{1}{2}$ cannot be so interpreted. Instead we invent the rule that a state of two like photons has statistical weight $\frac{1}{2}$ as great as that calculated as though the photons were different. This, generalized to $n$ identical photons, is the rule of Bose statistics.

The higher values of $n$ offer no problem. The $1/n!$ is interpreted combinatorially for different photons, and as a statistical factor when some are identical. For example, for all $n$ identical one obtains a factor ${\left(n!\right)}^{-1}\left(-{\beta}_{c}{{\beta}_{a}}^{*}\right)$ so that ${\left(n!\right)}^{-\frac{1}{2}}{\left(i\text{}{\beta}_{a}^{*}\right)}^{n}$ can be interpreted as the amplitude for emission (from no initial photons) of $n$ identical photons, in complete agreement with (61) for ${G}_{m0}$.

To obtain the amplitude for transitions in which neither the initial nor the final state is empty of photons we must consider the more general case of the division into three time regions (B1). This time we see that the factor which involves the coordinates in an entangled manner is $\mathrm{exp}i\left({R}_{ab}+{R}_{bc}+{R}_{ac}\right)$. It is to be expanded in the form $\sum _{i}$$\sum _{j}$${h}_{i}^{\u2033}\left(c\right){h}_{ij}^{\prime}\left(b\right){h}_{j}\left(a\right)$. Again the expansion in power series and development in Fourier series with a polarization sum will solve the problem. Thus the exponential is ${\sum}_{r}$${\sum}_{{l}_{1}}$${\sum}_{{l}_{2}}$${\left(i\text{}{R}_{ac}\right)}^{r}{\left(i\text{}{R}_{ab}\right)}^{{l}_{1}}{\left(i\text{}{R}_{bc}\right)}^{{l}_{2}}{\left({l}_{1}!\right)}^{-1}{\left({l}_{2}!\right)}^{-1}{\left(r!\right)}^{-1}$. Now the $R$ are written as Fourier series, one of the terms containing ${l}_{1}+{l}_{2}+r$ variables **K**. Since ${l}_{1}+r$ involve $a$, ${l}_{2}+r$ involve $c$ and ${l}_{1}+{l}_{2}$ involve $b$, this term will give the amplitude that ${l}_{1}+r$ photons are emitted during the interval $a$, of those ${l}_{1}$ are absorbed during $b$ but the remaining $r$, along with ${l}_{2}$ new ones emitted during $b$ go on to be absorbed during the interval $c$. We have therefore $n={l}_{1}+r$ photons in the state at time ${t}_{1}$ when $b$ begins, and $m={l}_{2}+r$ at ${t}_{2}$ when $b$ is over. They each are characterized by momentum vectors and polarizations. When these are different the factors ${\left({l}_{1}!\right)}^{-1}{\left({l}_{2}!\right)}^{-1}{\left(r!\right)}^{-1}$ are absorbed combinatorially. When some are equal we must invoke the rule of the statistical weights. For example, suppose all ${l}_{1}+{l}_{2}+r$ photons are identical. Then ${R}_{ab}=i\text{}{\beta}_{b}{{\beta}_{a}}^{*}$, ${R}_{bc}=i\text{}{\beta}_{c}{{\beta}_{b}}^{*}$, ${R}_{ac}=i\text{}{\beta}_{c}{{\beta}_{a}}^{*}$ so that our sum is

$${\sum}_{{l}_{1}}{\sum}_{{l}_{2}}{\sum}_{r}{\left({l}_{1}!{l}_{2}!r!\right)}^{-1}{\left(i\text{}{\beta}_{c}\right)}^{{l}_{2}+r}{\left(i\text{}{\beta}_{b}\right)}^{{l}_{1}}{\left(i\text{}{\beta}_{b}^{*}\right)}^{{l}_{2}}{\left(i\text{}{\beta}_{a}^{*}\right)}^{{l}_{1}+r}.$$ | (B4) |

Putting $m={l}_{2}+r$, $n={l}_{1}+r$, this is the sum on $n$ and $m$ of

$${\left(i\text{}{\beta}_{c}\right)}^{m}{\left(m!\right)}^{-\frac{1}{2}}\left[{\sum}_{r}{\left(m!n!\right)}^{\frac{1}{2}}{\left(\left(m-r\right)!\left(n-r\right)!r!\right)}^{-1}\times {\left(i\text{}{\beta}_{b}^{*}\right)}^{m-r}{\left(i\text{}{\beta}_{b}\right)}^{n-r}\right]{\left(n!\right)}^{-\frac{1}{2}}{\left(i\text{}{\beta}_{a}^{*}\right)}^{n}.$$ | (B5) |

The last factor we have seen is the amplitude for emission of $n$ photons during interval $a$, while the first factor is the amplitude for absorption of $m$ during $c$. The sum is therefore the factor for transition from $n$ to $m$ identical photons, in accordance with (57). We see the significance of the simple generating function (56).

We have therefore found rules for real photons in terms of those for virtual. The real photons are a way of representing and keeping track of those aspects of the past behavior which may influence the future.

If one starts from a theory involving an arbitrary modification of the direct interaction ${\delta}_{+}$ (or in more general situations) it is possible in this way to discover what kinds of states and physical entities will be involved if one tries to represent in the present all the information needed to predict the future. With the Hamiltonian method, which begins by assuming such a representation, it is difficult to suggest modifications of a general kind, for one cannot formulate the problem without having a complete representation of the characteristics of the intermediate states, the particles involved in interaction, etc. It is quite possible (in the author’s opinion, it is very likely) that we may discover that in nature the relation of past and future is so intimate for short durations that no simple representation of a present may exist. In such a case a theory could not find expression in Hamiltonian form.

An exactly similar analysis can be made just as easily starting with the general forms (46), (48). Also a coordinate representation of the photons could have been used instead of the familiar momentum one. One can deduce the rules (60), (61). Nothing essentially different is involved physically, however, so we shall not pursue the subject further here. Since they imply [23] all the rules for real photons, Eqs. (46), (47), (48) constitute a compact statement of all the laws of quantum electrodynamics. But they give divergent results. Can the result *after* charge and mass renormalization also be expressed to all orders in ${e}^{2}/\hslash \text{}c$ in a simple way?

An attempt has been made to find a differential wave equation for the propagation of an electron interacting with itself, analogous to the Dirac equation, but containing terms representing the self-action. Neglecting all effects of closed loops, one such equation has been found, but not much has been done with it. It is reported here for whatever value it may have.

An electron acting upon itself is, from one point of view, a complex system of a particle and a field of an indefinite number of photons. To find a differential law of propagation of such a system we must ask first what quantities known at one instant will permit the calculation of these same quantities an instant later. Clearly, a knowledge of the position of the particle is not enough. We should need to specify: (1) the amplitude that the electron is at $x$ and there are no photons in the field, (2) the amplitude the electron is at $x$ and there is one photon of such and such a kind in the field, (3) the amplitude there are two photons, etc. That is, a series of functions of ever increasing numbers of variables. Following this view, we shall be led to the wave equation of the theory of second quantization.

We may also take a different view. Suppose we know a quantity ${\Phi}_{{e}^{2}}\left[B,x\right]$, a spinor function of ${x}_{\mu}$, and functional of ${B}_{\mu}\left(1\right)$, defined as the amplitude that an electron arrives at ${x}_{\mu}$, with no photon in the field when it moves in an arbitrary external unquantized potential ${B}_{\mu}\left(1\right)$. We allow the electron also to interact with itself, but ${\Phi}_{{e}^{2}}$ is the amplitude at a given instant that there happens to be no photons present. As we have seen, a complete knowledge of this functional will also tell us the amplitude that the electron arrives at $x$ and there is just one photon, of form ${{A}_{\mu}}^{PH}\left(1\right)$ present. It is, from (60), $\int \left({\mathrm{\delta \Phi}}_{{e}^{2}}\left[B,x\right]/\delta \text{}{B}_{\mu}\left(1\right)\right){{A}_{\mu}}^{PH}\left(1\right)d{\tau}_{1}$.

Higher numbers of photons correspond to higher functional derivatives of ${\Phi}_{{e}^{2}}$. Therefore, ${\Phi}_{{e}^{2}}\left[B,x\right]$ contains all the information requisite for describing the state of the electron-photon system, and we may expect to find a differential equation for it. Actually it satisfies ($\nabla ={\gamma}_{\mu}\partial /\partial \text{}{x}_{\mu}$, $B={\gamma}_{\mu}{B}_{\mu}$),

$$\left(i\nabla -m\right){\Phi}_{{e}^{2}}\left[B,x\right]=B\left(x\right){\Phi}_{{e}^{2}}\left[B,x\right]+i\text{}{e}^{2}{\gamma}_{\mu}\int {\delta}_{+}\left({s}_{x1}^{2}\right)\left({\mathrm{\delta \Phi}}_{{e}^{2}}\left[B,x\right]/\delta {B}_{\mu}\left(1\right)\right)d{\tau}_{1}$$ | (C1) |

as may be seen from a physical argument.^{25} The operator $\left(i\nabla -m\right)$ operating on the $x$ coordinate of ${\Phi}_{{e}^{2}}$ should equal, from Dirac’s equation, the changes in ${\Phi}_{{e}^{2}}$ as we go from one position $x$ to a neighboring position due to the action of vector potentials. The term $B\left(x\right){\Phi}_{{e}^{2}}$ is the effect of the external potential. But ${\Phi}_{{e}^{2}}$ may also change for at the first position $x$ we may have had a photon present (amplitude that it was emitted at another point l is ${\mathrm{\delta \Phi}}_{{e}^{2}}/\delta {B}_{\mu}\left(1\right)$) which was absorbed at $x$ (amplitude photon released at 1 gets to $x$ is ${\delta}_{+}\left({{s}_{x1}}^{2}\right)$ where ${{s}_{x1}}^{2}$ is the squared invariant distance from 1 to $x$) acting as a vector potential there (factor ${\gamma}_{\mu}$). Effects of vacuum polarization are left out.

Expansion of the solution of (C1) in a power series in $B$ and ${e}^{2}$ starting from a free particle solution for a single electron, produces a series of terms which agree with the rules of II for action of potentials and virtual photons to various orders. It is another matter to use such an equation for the practical solution of a problem to all orders in ${e}^{2}$. It might be possible to represent the self-energy problem as the variational problem for $m$, stemming from (C1). The ${\delta}_{+}$ will first have to be modified to obtain a convergent result.

We are not in need of the general solution of (C1). (In fact, we have it in (46), (48) in terms of the solution ${T}_{0}\left[B\right]={\Phi}_{0}\left[B,x\right]$ of the ordinary Dirac equation $\left(i\nabla -m\right){\Phi}_{0}\left[B,x\right]=B{\Phi}_{0}\left[B,x\right]$. The general solution is too complicated, for complete knowledge of the motion of a self-acting electron in an arbitrary potential is essentially all of electrodynamics (because of the kind of relation of real and virtual processes discussed for photons in Appendix B, extended also to real and virtual pairs). Furthermore, it is easy to see that other quantities also satisfy (C1). Consider a system of many electrons, and single out some one for consideration, supposing all the others go from some definite initial state $i$ to some definite final state $f$. Let ${\Phi}_{{e}^{2}}\left[B,x\right]$ be the amplitude that the special electron arrives at $x$, there are no photons present, and the other electrons go from $i$ to $f$ when there is an external potential ${B}_{\mu}$ present (which ${B}_{\mu}$ also acts on the other electrons). Then ${\Phi}_{{e}^{2}}$ also satisfies (C1). Likewise the amplitude with closed loops (all other electrons go vacuum to vacuum) also satisfies (C1) including all vacuum polarization effects. The various problems correspond to different assumptions as to the dependence of ${\Phi}_{{e}^{2}}\left[B,x\right]$ on ${B}_{\mu}$ in the limit of zero ${e}^{2}$. The Eq. (C1) without further boundary conditions is probably too general to be useful.

^{1} | R P. Feynman, Phys. Rev. 76, 749 (1949), hereafter called I, and Phys. Rev. 76, 769 (1949), hereafter called II. |

^{2} | See in this connection also the papers of S. Tomonaga, Phys. Rev. 74, 224 (l948); S. Kanesawa and S. Tomonaga, Prog. Theoret. Phys. 3, 101 (1948); J. Schwinger, Phys. Rev. 76, 790 (1949); F. Dyson, Phys. Rev. 75, 1736 (1949); W. Pauli and F. Villars, Rev. Mod. Phys. 21, 434 (1949). The papers cited give references to previous work. |

^{3} | R. P. Feynman, Rev. Mod. Phys. 20, 367 (1948), hereafter called C. |

^{4} | E. Fermi, Rev. Mod. Phys. 4, 87 (1932). |

^{5} | W. Heitler, The Quantum Theory of Radiation, second edition (Oxford University Press, London, 1944). |

^{6} | The term in the sum for $n=m$ is obviously infinite but must be included for relativistic invariance. Our problem here is to re-express the usual (and divergent) form of electrodynamics in the form given in II. Modifications for dealing with the divergences are discussed in II and we shall not discuss them further here. |

^{7} | That (12) is correct, at least insofar as it depends on ${q}_{0}$, can be seen directly as follows. Let $\stackrel{\_}{q}\left(t\right)$ be the classical path which satisfies the boundary condition $\stackrel{\_}{q}\left({t}^{\prime}\right)={q}_{0}$, $\stackrel{\_}{q}\left({t}^{\u2033}\right)={q}_{j}$. Then in the integral defining $k$ replace each of the variables ${q}_{i}$ by ${q}_{i}={\stackrel{\_}{q}}_{i}+{y}_{i}$, (${\stackrel{\_}{q}}_{i}=\stackrel{\_}{q}\left({t}_{i}\right)$) that is, use the displacement ${y}_{i}$ from the classical path ${\stackrel{\_}{q}}_{i}$ as the coordinate rather than the absolute position. With the substitution ${q}_{i}={\stackrel{\_}{q}}_{i}+{y}_{i}$ in the action ${S}_{0}+{S}_{I}=\int \left(\frac{1}{2}{\stackrel{.}{q}}^{2}-\frac{1}{2}{\omega}^{2}{q}^{2}+\gamma q\right)dt=$$\int \left(\frac{1}{2}{\stackrel{.}{\stackrel{\_}{q}}}^{2}-\frac{1}{2}{\omega}^{2}{\stackrel{\_}{q}}^{2}+\gamma \stackrel{\_}{q}\right)dt+\int \left(\frac{1}{2}{\stackrel{.}{y}}^{2}-\frac{1}{2}{\omega}^{2}{y}^{2}\right)dt$ the terms linear in $y$ drop out by integrations by parts using the equation of motion $\stackrel{\_}{q}=-{\omega}^{2}\stackrel{\_}{q}+\gamma \left(t\right)$ for the classical path, and the boundary conditions $y\left({t}^{\prime}\right)=y\left({t}^{\u2033}\right)=0$. That this should occur should occasion no surprise, for the action functional is an extremum at $q\left(t\right)=\stackrel{\_}{q}\left(t\right)$ so that it will only depend to second order in the displacements $y$ from this extremal orbit $\stackrel{\_}{q}\left(t\right)$. Further, since the action functional is quadratic to begin with, it cannot depend on $y$ more than quadratically. Hence ${S}_{0}+{S}_{I}=Q+\int \left(\frac{1}{2}{\stackrel{.}{y}}^{2}-\frac{1}{2}{\omega}^{2}{y}^{2}\right)dt$ so that since $d{q}_{i}=d{y}_{i}$, $k\left({q}_{j},{t}^{\u2033};{q}_{0},{t}^{\prime}\right)=\mathrm{exp}\left(iQ\right)\int \mathrm{exp}\left(i\int \frac{1}{2}\left({\stackrel{.}{y}}^{2}-{\omega}^{2}{y}^{2}\right)dt\right)\mathcal{D}y\left(t\right)$. The factor following the $\mathrm{exp}i\text{}Q$ is the amplitude for a free oscillator to proceed from $y=0$at $t={t}^{\prime}$ to $y=0$ at ${t}^{\prime}={t}^{\u2033}$ and does not therefore depend on ${q}_{0}$, ${q}_{j}$, or $\gamma \left(t\right)$, being a function only of ${t}^{\u2033}-{t}^{\prime}$. [That it is actually ${\left(2\pi \text{}i\text{}{\omega}^{-1}\mathrm{sin}\omega \left({t}^{\u2033}-{t}^{\prime}\right)\right)}^{-\frac{1}{2}}$ can be demonstrated either by direct integration of the $y$ variables or by using some normalizing property of the kernels $k$, for example that ${G}_{00}$ for the case $\gamma =0$ must equal unity.] The expression for $Q$ given in C on page 386 is in error, the quantities ${q}_{0}$ and ${q}_{j}$ should be interchanged. |

^{8} | It is most convenient to define the state ${\phi}_{n}$ with the phase factor $\mathrm{exp}\left[-i\text{}\omega \left(n+\frac{1}{2}\omega \right){t}^{\prime}\right]$ and the final state with the factor $\mathrm{exp}\left[-i\text{}\omega \left(m+\frac{1}{2}\omega \right){t}^{\u2033}\right]$ so that the results will not depend on the particular times ${t}^{\prime}$, ${t}^{\u2033}$ chosen. |

^{9} | One can obtain the final result, that the total interaction is just $R$, in a formal manner starting from the Hamiltonian from which the longitudinal oscillators have not yet been eliminated. There are for each $K$ and $\mathrm{cos}$ or $\mathrm{sin}$, four oscillators ${q}_{\mu K}$ corresponding to the three components of the vector potential ($\mu =1,2,3$) and the scalar potential ($\mu =4$). It must then be assumed that the wave functions of the initial and final state of the $K$ oscillators is the function $\left(k/\pi \right)\mathrm{exp}\left[-\frac{1}{2}k\left({q}_{1K}^{2}+{q}_{2K}^{2}+{q}_{3K}^{2}-{q}_{4K}^{2}\right)\right]$. The wave function suggested here has only formal significance, of course, because the dependence on ${q}_{4K}$ is not square integrable, and cannot be normalized. If each oscillator were assumed actually in the ground state, the sign of the ${q}_{4K}$ term would be changed to positive, and the sign of the frequency in the contribution of these oscillators would be reversed (they would have negative energy). |

^{10} | The classical action for this problem is just $Sp+{R}^{\prime}$ where ${R}^{\prime}$ is the real part of the expression (24). In view of the generalization of the Lagrangian formulation of quantum mechanics suggested in Section 12 of C, one might have anticipated that $R$ would have been simply ${R}^{\prime}$. This corresponds, however, to boundary conditions other than no quanta present in past and future. It is harder to interpret physically. For a system enclosed in a light tight box, however, it appears likely that both $R$ and ${R}^{\prime}$ lead to the same results. |

^{11} | One can show from (25) how the correlated effect of many atoms at a distance produces on a given system the effects of an external potential. Formula (24) yields the result that this potential is that obtained from Liénard and Wiechert by retarded waves arising from the charges and currents resulting from the distant atoms making transitions. Assume the wave functions χ and ψ can be split into products of wave functions for system and distant atoms and expand $\mathrm{exp}\left(i\text{}R\right)$ assuming the effect of any individual distant atom is small. Coulomb potentials arise even from nearby particles it they are moving slowly. |

^{12} | The term corresponding to this for the self-energy expression (26) would give an integral over ${\delta}_{+}\left({\left(t-{t}^{\prime}\right)}^{2}-{\left({x}_{t}-{x}_{{t}^{\prime}}\right)}^{2}\right)$ which is evidently infinite and leads to the quadratically divergent self-energy. There is no such term for the Dirac electron, but there is for Klein-Gordon particles. We shall not discuss the infinities in this paper as they have already been discussed in II. |

^{13} | Alternatively, note that Eq. (37) is exact for arbitrarily large $\epsilon $ if the potential ${A}_{\mu}$ is constant. For if the potential in the Dirac equation is the gradient of a scalar function ${A}_{\mu}=\partial \chi /\partial {x}_{\mu}$ the potential may be removed by replacing the wave function by $\psi =\mathrm{exp}\left(-i\text{}\chi \right){\psi}^{\prime}$ (gauge transformation). This alters the kernel by a factor $\mathrm{exp}\left[-i\left(\chi \left(2\right)-\chi \left(1\right)\right)\right]$ owing to the change in the initial and final wave functions. A constant potential ${A}_{\mu}$ is the gradient of $\chi ={A}_{\mu}{x}_{\mu}$ and can be completely removed by this gauge transformation, so that the kernel differs from that of a free particle by the factor $\mathrm{exp}\left[-i\left({A}_{\mu}{x}_{\mu 2}-{A}_{\mu}{x}_{\mu 1}\right)\right]$ as in (37). |

^{14} | In changing the charge ${e}^{2}$ we mean to vary only the degree to which virtual photons are important. We do not contemplate changes in the influence of the external potentials. If one wishes, as $e$ is raised the strength of the potential is decreased proportionally so that $B$_{μ}, the potential times the charge $e$, is held constant. |

^{15} | The functional derivative is defined such that if $T$[$B$] is a number depending on the functions ${B}_{\mu}$(l), the first order variation in $T$ produced by a change from ${B}_{\mu}$ to ${B}_{\mu}+{\mathrm{\u25b3B}}_{\mu}$, is given by $T\left[B+\u25b3\text{}B\right]\mathrm{\u2013T}\left[B\right]=\int \left(\mathrm{\delta T}\left[B\right]/\delta \text{}{B}_{\mu}\left(1\right)\right){\mathrm{\u25b3B}}_{\mu}\left(1\right)d{\tau}_{1}$ the integral extending over all four-space ${x}_{\alpha ,1}.$ |

^{16} | That is, of course, those rules of II which apply to the unmodified electrodynamics of Dirac electrons. (The limitation excluding real photons in the initial and final states is removed in Sec. 8.) The same arguments clearly apply to nucleons interacting via neutral vector mesons, vector coupling. Other couplings require a minor extension of the argument. The modification to the ${\left({x}_{i+1}|{x}_{i}\right)}_{\epsilon}$,as in (37), produced by some couplings cannot very easily be written without using operators in the exponents. These operators can be treated as numbers if their order of operation is maintained to be always their order in time. This idea will be discussed and applied more generally in a succeeding paper. |

^{17} | For an alternative method starting directly from the formula (24) for virtual photons, see Appendix B. |

^{18} | F. Rohrlich (to be published). |

^{19} | The physical ideas involved in such a description are discussed in detail by Y. Nambu, Prog. Theor. Phys. 5, 82 (1950). An equation of type (A2) extended to the case of Dirac electrons has been studied by V. Fock, Physik Zeits. Sowjetunion 12, 404 (1937). |

^{20} | The factor $2i$ in front of ${I}_{+}$ is simply to make the definition of ${I}_{+}$ here agree with that in I and II. In II it operates with $p\xb7A+A\xb7p$ as a perturbation. But the perturbation coming from (A3) in a natural way by expansion of the exponential is $-\frac{1}{2}i\left(p\xb7A+A\xb7p\right)$. |

^{21} | Expression (A8) is closely related to Schwinger's parametric integral representation of these functions. For example, (A8) becomes formula (45) of F. Dyson, Phys. Rev. 75, 486 (1949)) for ${\Delta}_{F}\equiv {\Delta}^{\left(1\right)}-2\text{}i\text{}\stackrel{\_}{\Delta}\equiv 2\text{}i\text{}{I}_{+}$ if ${\left(2\alpha \right)}^{-1}$ is substituted for ${u}_{0}$. |

^{22} | The form (A10) suggests another interesting possibility for avoiding the divergences of quantum electrodynamics in this case. The divergences arise from the ${\delta}_{+}$ function when $u={u}^{\prime}$. We might restrict the integration in the double integral such that $\left|u-{u}^{\prime}\right|>\delta $ where δ is some finite quantity, very small compared with ${m}^{-2}$. More generally, we could keep the region $u={u}^{\prime}$ from contributing by including in the integrand a factor $F\left(u-{u}^{\prime}\right)$ where $F\left(x\right)\to 1$ for $x$ large compared to some δ, and $F\left(0\right)=0$ (e.g., $F\left(x\right)$ acts qualitatively like $1-\mathrm{exp}\left(-{x}^{2}{\delta}^{-2}\right)$. (Another way might be to replace $u$ by a discontinuous variable, that is, we do not use the limit in (A4) as $\epsilon \to 0$ but set $\epsilon =\delta $.) The idea is that two interactions would contribute very little in amplitude if they followed one another too rapidly in $u$. It is easily verified that this makes the otherwise divergent integrals finite. But whether the resulting formulas make good physical sense is hard to see. The action of a potential would now depend on the value of $u$ so that Eq. (A2), or its equivalent, would not be separable in $u$ so that $\frac{1}{2}{m}^{2}$ would no longer be a strict eigenvalue for all disturbances. High energy potentials could excite states corresponding to other eigenvalues, possibly thereby corresponding to other masses. This note is meant only as a speculation, for not enough work has been done in this direction to make sure that a reasonable physical theory can be developed along these lines. (What little work has been done was not promising.) Analogous modifications can also be made for Dirac electrons. |

^{23} | The formulas for real processes deduced in this way are strictly limited to the case in which the light comes from sources which are originally dark, and that eventually all light emitted is absorbed again. We can only extend it to the case for which these restrictions do not hold by hypothesis, namely, that the details of the scattering process are independent of these characteristics of the light source and of the eventual disposition of the scattered light. The argument of the text gives a method for discovering formulas for real processes when no more than the formula for virtual processes is at hand. But with this method belief in the general validity of the resulting formulas must rest on the physical reasonableness of the above-mentioned hypothesis. |

^{24} | Usually only two polarizations transverse to the propagation vector $K$ are used. This can be accomplished by a further rearrangement of terms corresponding to the reverse of the steps leading from (17) to (19). We omit the details here as it is well known that either formulation gives the same results. See II, Section 8. |

^{25} | Its general validity can also be demonstrated mathematically from (45). The amplitude for arriving at $x$ with no photons in the field with virtual photon coupling ${e}^{2}$ is a transition amplitude. It must, therefore, satisfy (45) with ${T}_{{e}^{2}}\left[B\right]={\Phi}_{{e}^{2}}\left[B,x\right]$ for any $x$. Hence show that the quantity ${C}_{{e}^{2}}\left[B,x\right]=\left(i\text{}\nabla -m-B\left(x\right)\right){\Phi}_{{e}^{2}}\left[B,x\right]-i\text{}{e}^{2}{\gamma}_{\mu}\int {\delta}_{+}\left({s}_{x1}^{2}\right)\left({\mathrm{\delta \Phi}}_{{e}^{2}}\left[B,x\right]/\delta {B}_{\mu}\left(1\right)\right)d{\tau}_{1}$ also satisfies Eq. (45) by substituting ${C}_{{e}^{2}}\left[B,x\right]$ for ${T}_{{e}^{2}}\left[B\right]$ in (45) and using the fact that ${\Phi}_{{e}^{2}}\left[B,x\right]$ satisfies (45). Hence if ${C}_{0}\left[B,x\right]=0$ then ${C}_{{e}^{2}}\left[B,x\right]=0$ for all ${e}^{2}$. But ${C}_{{e}^{2}}\left[B,x\right]=0$ means that ${\Phi}_{{e}^{2}}\left[B,x\right]$ satisfies (C1). Therefore, that solution ${\Phi}_{{e}^{2}}\left[B,x\right]$ of (45) which also satisfies $\left(i\text{}\nabla -m-B\left(x\right)\right){\Phi}_{0}\left[B,x\right]=0$ (the propagation of a free electron without virtual photons) is a solution of (C1) as we wished to show. Equation (C1) may be more convenient than (45) for some purposes for it does not involve differentiation with respect to the coupling constant, and is more analogous to a wave equation. |