In theoretical physics, not only the solution of differential equations is important, but the number of solutions as well. In quantum mechanics, as an example, boundary problems normally have several solutions and all these solutions are used to describe the behavior of an object. Discarding any of the solutions excludes one or another effect from consideration, and the description becomes inadequate and contrary to the fact. If a boundary problem has several solutions and only one of them is chosen, we should justify why we do so. For practically important boundary problems, the question of solution uniqueness as such is often not posed. It is thought that any obtained solution is automatically unique in accordance with the Cauchy-Kovalevsky theorem  ; however, the theorem considers uniqueness for a problem only at a certain point and not over the entire domain. Proofs of uniqueness are provided for particular differential equations, e.g., two-dimensional potential boundary problems. However, conditions of the uniqueness theorem, e.g. boundedness of a solution or problem domain, may contradict actual physical phenomena.
In this work, using the field of a concentrated charge (electron, proton, neutron) as an example, it is shown that the well-known equation of mathematical physics has the second, earlier unknown solution. Taking into account this solution radically alters the interpretation of some physical processes. Surely, uniqueness of the solutions of differential equations is studied by mathematicians   , but their interest (and this is quite natural) is in mathematical rather than physical problems. The solution uniqueness of the electron problem has not been studied so far.
2. Non-Potential Solution
The electrostatic field of a point source obeys the equation
where is the electrostatic field, is the electron charge, is the dielectric constant, is the radius-vector, and is the delta (Dirac) function. It is assumed that is a potential vector with the potential function :
Equations (1) and (3) give the Poisson equation ( is the Laplacian):
The solution of (4) is the following
In the Cartesian and spherical coordinates, the field vector is written as
In this work I will use pointed brackets to denote the vector component in the spherical coordinates. The sequence of the vector components in these brackets is as follows: radial, zenith, azimuth, e.g. . The Cartesian components will be denoted by the braces, as it is accepted.
If two potential charges are located at a distance of , then the energy of their interaction is written, according to (5), as
The spherical and Cartesian coordinates are related as
The flux of the vector through a closed sphere centered at the origin of the coordinates
The assumption about field potentiality (3) has no physical grounds and limits the class of the functions that can be a solution of (1). Hence, one needed to prove that (6) is a unique solution of Equation (1), but this was not done. If there is any other solution of (1), it is necessary to justify the choice of (6) as a physically feasible solution.
Let us find a different solution of (1) than the potential solution (6) by the procedure described in    . The auxiliary vector has the form
For the vector
The vector is neither purely potential nor purely vortex because its divergence and rotor are not identically equal to zero. If we manage to find a solution for the equation
is the solution of equation  , where C is an arbitrary constant. One can replace without sacrifice of generality
Equation (12) is solved by the variable separation method
The function must depend on , and , one can use and the first term in (12) vanishes. The right side of (12) is invariant with . Let . Hence, the third term in (12) vanishes and (12) takes the form
Its integration with respect to gives:
where is an arbitrary constant, . Further integration yields
Although is an unbounded function at , the field vector produced by the function is bounded at this point. It can be shown (which is done at the end of Section 2) that one of the condition of the minimum energy is , and hence the second term in (18) is excluded from consideration. As a result, we have or in the Cartesian coordinates
From the above relation and from (15) it follows that
and from (13) one obtains
Here is selected to fulfill (22). Substitution of (21) in (1) proves that (21) is a solution of (1). The vector is nonpotential, since .
The vector is axisymmetric about the z axis, since does not depend on . The flux of the vector through a closed sphere centered at the origin of the coordinates is
The flux (22) is constant and is independent of the sphere radius, which suggests that
as is for in (6). From (6) and (21) it follows that the vectors and are parallel.
The vectors and satisfy the boundary problem (1). Consider the linear combination of the vectors
where are the arbitrary constants. The choice of the constants may not affect the charge and the charge must remain constant and equal to the electron charge.
The charge value produced by the fields (6) and (21) is determined by the relations (9) and (22). If the values of the parameters and are selected in this way
where is the arbitrary parameter, then the electric field will be represented like this
The vector flux in (30) through the sphere with the center at the origin of the coordinates
is zero, which means that the electric field (26) is produced by a zero charge. This is the property of the nonpotential field, whereas it is impossible in the case of the classical potential field. The field (26) at a certain value of describes the neutron electric field (see the details in Section 4), therefore, I will refer to the value in (26) as a neutron component, and the parameter -a neutron charge. In the spherical coordinates (26), the field has only one radial component independent of the azimuthal angle , which means that the neutron field is axisymmetric, with the applicate, in this case, being the axis of symmetry. In the general case, the orientation of the neutron field will be specified if the direction of the axis of symmetry is indicated.
If the potential field (6) is added (26), one will derive a general solution of the problem of the electron (proton) in the form
Indeed, the flux of the vector through a closed surface is equal to e, i.e. the elementary charge, according to (9) and (27). The non potentiality of the electric field charge results in the appearance of a “free” component of the electric field.
The potential energy of the electrostatic field (24) is specified by the relation 
Using (6) and (21), we obtain:
The lower index by the integration limit means that the integration is performed with respect to the whole space with the exception of the sphere of the small radius with the center at the origin of the coordinates. This spherical environment has been excluded from consideration to avoid the integral divergence at the origin of the coordinates, is the radius of truncation and the environment will be referred to as a source to emphasize the fact that the electric field possessing potential and nonpotential components is its product. A particle can have both the potential and neutron component, each of them having its own radius of truncation. The integration is performed over the entire space. Relation (30) is broken down into three terms. The first term
is the self-energy for the classical potential electron component. The second term
is the self-energy for the non potential component. The third term
is the interaction energy for the potential and nonpotential components. Here is the potential component energy, is the nonpotential component energy, is the energy of the interaction between the potential and nonpotential components. Since , there is no interaction between the potential and neutron components. The presence of the neutron component only increases the intrinsic energy of the particle.
The value of can be estimated equating the electrostatic energy of the particle to the Einstein energy , where is the particle mass and is the light speed. The existence of the neutron component as well as the potential one is due to the properties of the source. The radius of truncation for different components can be different, which the lower index n by the radius does denote (32). For estimation let us assume the radius of truncation of the particle to be the same both of the potential and non potential components, since they are produced by the same source. Adding (31) and (32) and equating the sum to one obtains
The value of the multiplier is unknown. To estimate the truncation radius of a charged particle, one can assume , i.e. estimating the radius by the potential component one derives, as a result
As a matter of fact, the real radius of truncation will be slightly larger due to the energy of the neutron component. For a proton of the mass  gives , for an electron with the mass one obtains . For a neutron instead of (34) one has
Here, on the left-hand side of the equality there is unknown parameter , is the neutron mass. For a charged particle with the field (28) the parameters of the neutron component are determined by its mass, the magnetic and electric moment, each particle will have its own values of and . At present there are not enough experimental data for their unambiguous calculation.
In addition to the main component T in (28), the dipole and multipole components, whose fields are derived from (28) by differentiation with respect to the spatial coordinates, can contribute to the total field. These components can also be produced by a source just as the main components. For the multipole components the size of the sources can differ and their value will affect the total system energy.
Let explain why we put in (18) . In spherical coordinates, the gradient of the second term in the right side of (18)
Then, the interaction energy for the classical component (6) and (37) is written as
The interaction energy for the nonpotential component (21) and (37) is zero, too. Hence, the minimum energy (29) is reached with .
The formula for the neutron electric field (26) contains the arbitrary constant . There are two possible variants: and . This agrees well with the known fact that there are two types of neutrons: a neutron and an antineutron. These are two different particles. The sign for must be chosen proceeding from the interaction. A proton carries the positive charge, annihilates with an antineutron, they are attracted and, therefore, with the antineutron and with the neutron .
An arbitrary rotation can be realized by three rotations about the coordinate axes. First, a rotation about the z axis through the angle is made. In this rotation, the x and y axes take the positions and . The second rotation is made about the axis through the angle and the axis assumes the position . The third rotation is made about the axis through the angle . For the charge at the origin of the coordinates, the initial orientation is chosen such that the symmetry axis coincides with the z axis. In this case, rotation about the z axis through any angle changes nothing. Therefore, it can be assumed that . The matrices for the rotations about the y and x axes are written as 
The vectors are transformed using the matrices and the coordinates using the transposed matrix . Transformation of the vector in (26) gives
Because of the stress field symmetry, the angles are variable as follows
5. Interaction Energy
Let us denote the external homogeneous electric field by the tilde. It has an arbitrary value and orientation, and in the spherical coordinates its radial component will be written as
The other spherical components are of no interest, since only the radial component of the neutron component is other than zero. Here is the external field in the Cartesian coordinates. The interaction energy of the neutron component with the external field is
The interaction energy of a neutron and a homogeneous electric field is equal to zero. This means that the force acting on the neutron on the side of the homogeneous electric field is also equal to zero. The same equality may not be true for the nonuniform external field. That a neutron, under certain conditions, interacts with an electric field is an experimental fact  . The interaction energy of the particles with a homogeneous electric field is completely determined by the potential component. The non-potential component does not participate in the interaction. In a homogeneous electric field it is impossible to tell if an electron (a proton) possesses the non-potential component. It behaves like a particle with a purely potential field.
5.1. Interaction of the Potential Proton and Neutron
From (33) it follows that there is no interaction between the potential and neutron component of the same particle. However, if the particles are different, for instance, at the point on the axis of abscissa there is a particle with a pure potential field (6), and at the point there is a particle with an arbitrary oriented neutron field (40), then between them there is interaction with the energy
The lower indices − and + denote that the center of the given component is located on the axis of abscissa at the points and , respectively. Integration is performed over the whole space with the exception of two spherical regions with the radii at the point and
at the point . Here it is taken into account that the particles can have different radii of truncation. The potential field of the electric charge of the particle at the point is
and the neutron field at the point is
Integration (45) is performed by passing to the bipolar coordinates, there is a scalar under the integral and the Cartesian coordinates are replaced by the bipolar ones, according to the formulas
The Jacobian is
The field (48) is derived from (40) by the transformation of the shift along the axis of abscissa , and the field (47) from (6) by the shift to the opposite side . Here and are
As a result, one obtains the interaction energy value (45) in the form
k is the coefficient from (46). The condition of is necessary to make the calculations correctly. This condition is fulfilled for real variants for which the value
Using the above value as well as an approximate representation for small values of x
it is possible to simplify the Formula (56)
Thus, the energy of interaction between the neutron and the potential charge has a weak logarithmic dependence on the value of the truncation radius and is inversely proportional to the distance between the particles. It also depends on the neutron orientation and can have both positive and negative values, that is, the interaction of the potential charge with the neutron component can be both attraction and repulsion. The different value of the radii of truncation was taken into account in the initial Formula (45). However, the influence of their difference on the final result turned out negligible. Further, I am not going to take into consideration this difference.
If there is a neutron at the point and a potential charge at the point , i.e. the particles exchange their places, then in (60) it is necessary to replace the angle by . The transformation is performed according to (62), and the energy found
depends only on one angle . This energy is needed when calculating the energy of the interaction of two charged particles, each of them having its own neutron component.
5.2. Interaction of Two Arbitrary Oriented Neutrons
It is always possible to choose the coordinate orientation in such a way that the symmetric axis of the charge located at the point would lie in the xz plane. Due to the choice of the above orientation of the coordinate system, the angle , and this neutron orientation is determined only by the angle . This choice of the orientation of the coordinate system does not decrease the generality of consideration. The neutron field at the point is derived from (48) by the replacement
As a result, one obtains the expression for the neutron at the point :
The interaction energy of the neutron-neutron will be written as
Neutrons interact with each other due to the nonpotential electric field. The energy (64) depends both on the charge signs and the mutual orientation. Depending on it the same neutrons can be both attracted and repulsed. At present there are no experimental data on the neutron-neutron interaction. The formula (64) for large values of is simplified
For instance, for the parameters the interaction energy value will be such
5.3. Interaction Energy of Non-Potential Particles
The obtained results allow defining the proton-proton interaction energy. The total energy of the field of two non-potential charged particles will be written, as
Here (7), (60), (61), (65) are used as well as the assumption that the truncation radii of the potential and neutron components are the same.
6. Potential and Non-Potential Magnetic Field
A magnetic field like an electric field can have potential and non-potential components. To describe a potential magnetic field, the notion of the vector-poten- tial  is used. For a point magnetic charge (a monopole) located at the origin of the coordinates and oriented along the z-axis, the vector-potential is written as
is the magnetic constant, is the value of the monopole charge. Monopoles have not been found in reality, so it is not clear if they exist in nature. Here the notion of a monopole is used as a mathematical object to solve problems. Taking the rotor from (68), one obtains the field strength of the magnetic monopole in the form
From here it is clear that in the spherical coordinates only the azimuthal component of the monopole strength is other than zero. Differentiating (68) with respect to z and replacing by the magnetic moment p, one obtains the vector-potential of the point magnetic dipole in the form
Magnetic dipoles really exist, for instance, the magnetic electron moment. The rotor from (70) is the magnetic field induction produced by the magnetic dipole.
In the spherical coordinates only the azimuthal component is other than zero, it decreases with a distance as . The field is symmetric relative to the plane , due to it is even function of z. The lines of force are the essence of a circle. The energy of the potential component of the magnetic dipole (31) is calculated in the same way as that of the potential component of the electric charge:
In addition to the potential component of the magnetic field given by the relation (71), a magnetic field can have a non-potential component. To find it, one can use the technique described in Section 2. Let us take the auxiliary function in the form
Here is the constant which characterizes the deviation of the magnetic field from the potential one. At the field is potential. The vector (73), i.e. the strength of the field produced by it, will be as follows
The divergence (73) will be as follows
As neither (74) nor (75) are identically equal to zero, the vector (73) is neither potential nor rotational. As in the case of Section 2, it is necessary to find the potential vector so that
Since the vector is potential, it can be represented as the gradient of a certain function :
Substituting this relation into (76) one obtains
In the spherical coordinate system the equation is written as  :
As in Section 0, let us assume that the unknown function can be repre- sented in the form of the product
Assuming that one gets the equation
The left-hand side (82) is the sum of two terms, with the augend depending only on , and the addend depending only on . The right-hand side (82) is the product of two unknown functions, one of them depending only on , and the other only on . This is possible, when one of the terms is a constant value as well as the function corresponding to this term in the right-hand side (82). This assumption proves to be correct for the addend, and in this case, the value
Here the constant is included into the arbitrary constant , then
and the Equation (82) goes into
Let us denote
Then using the differentiating rules
and (85) goes into
since one has the relations
The homogeneous equation for (90) will be written as
Its linearly independent solutions are
Following the methods of solving such equations  , let us consider the values
Integrating them one obtains
The solution of the equation (90) will be written as
Here are the arbitrary constants. Substituting there in (86), one derives
Substituting (95) into (85) one can make sure that (95) is the solution of this differential equation. Using (81) one derives the solution of the equation (80) in the form
Or turning to the Cartesian coordinates
Differentiating (97) one obtains
Hence it follows (75)
Thus, one arbitrary constant in (97) has been defined. The vector
is an unlimited value at at any distances from the origin of the coordinates, although its divergence is limited. In (100) the value of the vector is given in the spherical coordinates. The total field vector
includes two arbitrary constants and. The energy density
produced by this vector field will be an unlimited value at, therefore, the total energy
will be unlimited too. During integration (103) first, let us find a certain integral with respect to, then to, after that find some integral with respect to, because a certain integral with respect to diverges at, so as a result, one derives
In order to avoid the divergence with respect to, one should use in (104):
Hence one finds
Substituting these values in (100) and (101) one obtains a non-potential vector in the form
Differentiating (105) with respect to the spatial coordinates in the Cartesian system of coordinates, one obtains multipoles of the corresponding order. For example, differentiating with respect to the coordinate z, one derives a magnetic dipole with the axis of symmetry along the z-axis in the form
This field is axially symmetric with the z-axis as the axis of symmetry, since the vector (106) is invariant with respect to the replacement. The azimuthal vector components (71), (106), (107) are proportional to one another, therefore, in the calculation of the interaction energy these components will yield a nonzero contribution to the interaction energy. This contribution, depending on the orientation of the neutron components can be both positive and negative. With the negative contribution the addition of the components (71) and (106) to (107) can decrease the total energy of the system, therefore, to correctly estimate the value of the magnetic moment of an elementary particle, it is necessary to take into consideration the influence of (71) and (106).
7. Magnetic Moment
For the electric field of the neutron (26) one has
From here it is clear that in the spherical coordinate system only the azimuthal component of the vector (107) is other than zero. This field in its structure corresponds to the field of the potential magnetic dipole (71). For the correspondence to be complete, it necessary to multiply (107) by a certain value of having the dimensions of length. Then (107) gives the electric dipole field of the kind
is the neutron dipole electric moment. Nowadays this value has not been measured yet, although there have been some attempts to do it   . According to the Maxwell equation, there is interdependence between the rotor of the electric field strength and the vector of the magnetic field in the form
Here is the source of magnetic induction (71). In modern physics the sources of the magnetic field are not taken into account in the above equation, which is explained by the fact that monopoles have not been found yet. However, the magnetic dipole sources do exist, and when necessary they should be included into the Maxwell equation. Since a stationary case is considered, i.e. all values do not depend on time, then
Hence it follows that
The coefficient of proportionality has the dimensionality s−1. If (112) is multiplied by the coefficient with the dimensions of the length, one will turn to the electric dipole (109) in the addend (112), and there will be the coefficient having the dimensions of speed before the magnetic induction. Let us choose the value of so that, where is the light speed. Then one obtains the well-known relation in electrodynamics between the electric field strength and the magnetic induction. Comparing (71) and (107) one derives the expression for the electric dipole through the magnetic moment
Here it is denoted that, which is the dipole moment. Substituting the value of the neutron magnetic moment into (113) one obtains the value of the neutron electric dipole.
Modern experimental methods   do not possess the required resolution to measure the above value of the neutron electric moment, so the question of its existence remains open. Thus, the non-potential electric field gives rise to the magnetic dipole moment of elementary particles. This magnetic dipole possesses nonzero energy of interaction with the dipole sources of the magnetic field (71) and (106). The total magnetic dipole moment will be their linear combination. Determination of the coefficients of this linear combination is outside the scope of the present paper.
As shown above, if one takes into consideration the non-potential electric field, then the known particles prove to have properties unknown earlier and the known properties find a simple explanation. It is common knowledge that scattering of elementary particles has the property of indefiniteness, which is explained by a special character of the microworld having no analogy in the macroscopic world. Taking into account the non-potential components of the electric field eliminates the “mysteriousness”. Orientation of the neutron components of the colliding particles is arbitrary, uncontrolled at the moment, and the interaction force depends on it. Therefore, the result of scattering of the particles will have a fraction of indefiniteness. Another source of indefiniteness is the fact that the neutron component is proportional to the parameter. When particles do not interact, the neutron component is in the stationary state, at which the value of is determined by the source. If the interaction is weak, then after it the orientation of the neutron components of the particles will not change. This is elastic scattering. However, if the interaction is strong, then the parameter of the neutron components may increase during collision increasing the neutron component, and the particles go into the excited state, in which the neutron components get a nonstationary addition and there appear electro- magnetic oscillations. The orientation of the neutron components after the collision may change and result in radiation. This is nonelastic scattering.
Another characteristic of nonpotentiality is the existence of the proton neu- tron, electron-neutron and neutron-neutron components of interaction. Neutrons and protons have similar properties and for this reason, they belong to the same group of particles―nucleons. In terms of the potential theory, it is not at all clear. Neutrons have no charge and they have “no” electric field, so they must not interact electrostatically with a proton. However, the interaction takes place even at small energies, for example, in the case of cold neutrons  . There is no explanation for the fact. Thus, in the potential theory the interaction of a neutron with a proton is a puzzle. For the non-potential field there is no problem like this. The energy of interaction between a neutron and a proton is derived by taking the sum of (60) and (65)
The energy of interaction of a proton with a proton is given by the relation (67). When comparing (67) and (114) one can see that the coefficient in front of the square brackets of both formulas, as can be checked by substitution of various angle values and for, is of the same order of magnitude. For instance, for one obtains for the proton-proton interaction (67)
and the energy of the proton-neutron interaction (114) is written as
Using the relation (58) one derives that for the values of (115) and (116) are of the same order of magnitude. Therefore, due to the electric field, the interaction between a neutron and a proton is approximately the same as that between a proton and proton. That explains the similarity of these particles. In terms of the potential theory, the existence of an antineutron and its annihilation with a proton are not clear either. What is the difference between a neutron and an antineutron? In the light of the non-potential theory, these are the natural consequences of the properties of the intrinsic electric fields of the particles. Choosing the corresponding sign before, one obtains the negative interaction energy (114). The electric field non-potentiality results in the elementary particles acquiring the magnetic moment. The particles also have the electric dipole moment, which has not been measured so far due to the insufficient experiment accuracy.
 Ivanchin, A.G. (1987) An Application of Potential Theory for Solving the Plane Elastic Crack Problem. In: Goerlich, P., Ed., Physica Status Solidi (a), Akademie Verlag Berlin, Berlin, 437-447.
 Wolfram Mathematica Online.
 Serebrov, A.P., et al. (1983) New Measurements of the Neutron Electric Dipole Moment with the Petersburg Nuclear Physics Institute Double Chamber Electric Dipole Moment Spectrometer. Physics of Particles and Nuclei Letters, 12, 286-296.