The remaining nontrivial task is the derivation of a local differential equation for and from the integral equation (37). As our essential constraint we will use, besides general principles of simplicity (like homogeneity and isotropy of space) the principle of maximal disorder, as realized by the requirement of minimal Fisher information. Using the abbreviation

| (39) |

the general solution of (37) may be written in the form

| (40) |

where the three functions have to vanish upon integration over and are otherwise arbitrary. If we restrict ourselves to an isotropic law, we may write

| (41) |

Then, our problem is to find a function which fulfills the differential equation

| (42) |

and condition (31). The method used in I for a one-dimensional situation, to determine from the requirement of minimal Fisher information, remains essentially unchanged in the present three-dimensional case. The reader is referred to the detailed explanations reported in I.

In I it has been shown that this principle of maximal disorder leads to an anomalous variational problem and to the following conditions for our unknown function :

where contains only derivatives of up to second order and does not explicitely depend on . If Eq. (43) is taken into account, the Euler-Lagrange equations of the variational problem (44) lead to the following differential equation
| (45) |

for the variable . Eq. (45) is a straightforward generalization of the corresponding one-dimensional relation [equation (68) of I] to three spatial dimensions.

Besides (45) a further (consistency) condition exists, which leads to a simplification of the problem. The function may depend on second order derivatives of but this dependence must be of a special form not leading to any terms in the Euler-Lagrange equations [according to (43) our final differential equation for and must not contain higher than second order derivatives of ]. Consequently, the first term in Eq. (45) (as well as the sum of the remaining terms) has to vanish separately and (45) can be replaced by the two equations

In I a new derivation of Fisher’s functional has been obtained, using the general solution of the one-dimensional version of (45), as well as the so-called composition law. In the present three-dimensional situation we set ourselves a less ambitious aim. We know that Fisher’s functional describes the maximal amount of disorder. If we are able to find a solution of (46), (47) that agrees with this functional (besides ’null-terms’ giving no contribution to the Euler-Lagrange equations) then we will accept it as our correct solution. It is easy to see that this solution is given by
| (48) |

where is an arbitrary constant. Eq. (48) presents again the three-dimensional (and isotropic) generalization of the one-dimensional result obtained in I. By means of the identity

| (49) |

it is easily verified that the solution (48) obeys also condition (31). Using the decomposition (16) and renaming according to , the continuity equation (3) and the second differential equation (43) respectively, take the form

The function occurring in (50), (51) is single-valued but not unique (not gauge-invariant). If now the complex-valued variable
| (52) |

is introduced, the two equations (50), (51) may be written in compact form as real and imaginary parts of the linear differential equation

| (53) |

which completes our derivation of Schrödinger’s equation in the presence of a gauge field.

Eq. (53) is in manifest gauge-invariant form. The gauge-invariant derivatives of with respect to and correspond to the two brackets in (53). In particular, the canonical momentum corresponds to the momentum operator proportional to . Very frequently, Eq. (53) is written in the form

| (54) |

with the Hamilton operator

| (55) |

Our final result, Eqs. (54), (55), agrees with the result of the conventional quantization procedure. In its simplest form, the latter starts from the classical relation , where is the Hamiltonian of a classical particle in a conservative force field, and is its energy. To perform a ”canonical quantization” means to replace and by differential expressions according to (1) and let then act both sides of the equation on states of a function space. The ’black magic’ involved in this process has been eliminated, or at least dramatically reduced, in the present approach, where Eqs. (54), (55) have been derived from a set of assumptions which can all be interpreted in physical terms.

The Hamiltonian (55) depends on the potentials and and is consequently a non-unique (not gauge-invariant) mathematical object. The same is true for the time-development operator which is an operator function of , see e.g. [34]. This non-uniqueness is a problem if is interpreted as a quantity ruling the time-evolution of a single particle. It is no problem from the point of view of the SI where and are primarily convenient mathematical objects which occur in a natural way if the time-dependence of statistically relevant (uniquely defined) quantities, like expectation values and transition probabilities, is to be calculated.