Theoretical physics makes a wide use of differential equations for which only a potential solution is applied. The possibility that these equations may have a non-potential solution is ruled out and not considered. In this paper an exact non-potential solution of the continuity equation is described. The electric field of an elementary charged particle consists of two components: the known Potential Component (PC) produced by the charge and the earlier unknown Non-potential Component (NC) with a zero charge. Charged particles have both components, while a neutron has only the NC. The proton and neutron NC ensures similarity of their properties. The PC is spherically symmetric and NC is axisymmetric. Therefore, to describe an elementary particle, one should take into account both its spatial coordinates and the NC orientation. The particle interaction is determined by their NC mutual orientation. Neglecting the latter leads to indefiniteness of the interaction result. In a homogeneous electric field, the force acting on the NC is zero. Therefore, a charged particle possessing the NC will behave like a potential one. In an inhomogeneous field, the situation is principally different. Due to the NC there occurs an interaction between a neutron and a proton. The non-potential field results in the existence of two types of neutrons: a neutron and an antineutron. A neutron repels from a proton ensuring scattering of neutrons on protons. An antineutron is attracted to a proton leading to its annihilation. The NC produces the magnetic dipole moment of an elementary particle.
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 [
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 [
The electrostatic field of a point source obeys the equation
∇ ⋅ E = − e 4 π ε δ ( r ) (1)
E → 0, r → ∞ (2)
where E is the electrostatic field, e is the electron charge, ε is the dielectric constant, r is the radius-vector, and δ ( r ) is the delta (Dirac) function. It is assumed that E is a potential vector with the potential function Φ :
E = ∇ ⋅ Φ (3)
Equations (1) and (3) give the Poisson equation ( Δ is the Laplacian):
Δ ( Φ ) = − e 4 π ε δ ( r ) (4)
The solution of (4) is the following
Φ = e 4 π ε r (5)
In the Cartesian and spherical coordinates, the field vector E is written as
E = e 4 π ε r 3 r = e 4 π ε 〈 1 r 2 , 0 , 0 〉 (6)
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. 〈 E r , E ϑ , E φ 〉 . The Cartesian components will be denoted by the braces, as it is accepted.
If two potential charges are located at a distance of 2 a , then the energy of their interaction is written, according to (5), as
W i = e 2 8 π ε a (7)
The spherical and Cartesian coordinates are related as
x = r sin ϑ cos φ , y = r sin ϑ sin φ , z = r cos ϑ (8)
The flux of the vector ε E through a closed sphere centered at the origin of the coordinates
ε ∮ S E r d S = e 4 π ε ∫ 0 2π d φ ∫ 0 π E r r 2 sin ϑ d ϑ = e (9)
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 [
U + = { x r 3 , y r 3 ,0 } (10)
For the vector U +
∇ ⋅ ( U + ) = 1 − 3 cos 2 ϑ r 3 (11)
The vector U + 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
Δ Ψ = 1 r 2 ∂ ∂ r ( r 2 ∂ Ψ ∂ r ) + 1 r 2 sin ϑ ∂ ∂ ϑ ( sin ϑ ∂ Ψ ∂ ϑ ) + 1 r 2 sin 2 ϑ ∂ 2 Ψ ∂ φ 2 = 1 − 3 cos 2 ϑ r 3 (12)
the vector
U = C ( U + − ∇ Ψ ) (13)
is the solution of equation [
e 2 π ε C → C (14)
Equation (12) is solved by the variable separation method
Ψ = R ( r ) Θ ( ϑ ) Φ ( φ ) (15)
The function Ψ must depend on r , ϑ and φ , one can use R ( r ) = 1 / r and the first term in (12) vanishes. The right side of (12) is invariant with φ . Let Φ ( φ ) = const . Hence, the third term in (12) vanishes and (12) takes the form
d d ϑ ( sin ϑ d Θ d ϑ ) = ( 1 − 3 cos 2 ϑ ) sin ϑ (16)
Its integration with respect to ϑ gives:
d Θ d ϑ = − sin ϑ cos ϑ + A sin ϑ (17)
where A is an arbitrary constant, − ∞ < A < + ∞ . Further integration yields
Θ = 1 2 cos 2 ϑ + A ln tan ϑ 2 (18)
Although Θ is an unbounded function at ϑ = 0 , 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 A = 0 , and hence the second term in (18) is excluded from consideration. As a result, we have Θ = ( cos 2 ϑ ) / 2 or in the Cartesian coordinates
Θ = z 2 2 ( x 2 + y 2 + z 2 ) (19)
From the above relation and from (15) it follows that
Ψ = z 2 2 ( x 2 + y 2 + z 2 ) 3 / 2 (20)
and from (13) one obtains
U = e 4 π ε z 2 − 2 x 2 − 2 y 2 ( x 2 + y 2 + z 2 ) 5 / 2 r = e 4 π ε r 2 〈 2 − 3 cos 2 ϑ , 0 , 0 〉 (21)
Here C = 1 is selected to fulfill (22). Substitution of (21) in (1) proves that (21) is a solution of (1). The vector U is nonpotential, since ∇ × U ≡ 0 .
The vector U is axisymmetric about the z axis, since U r does not depend on φ . The flux of the vector U through a closed sphere centered at the origin of the coordinates is
ε ∮ S U r d S = e 4 π ε ∫ 0 2π d φ ∫ 0 π U r r 2 sin ϑ d ϑ = − e (22)
The flux (22) is constant and is independent of the sphere radius, which suggests that
∇ ⋅ U = − e 4 π ε δ ( r ) (23)
as is for E in (6). From (6) and (21) it follows that the vectors E and U are parallel.
The vectors U and E satisfy the boundary problem (1). Consider the linear combination of the vectors
λ 1 E + λ 2 U (24)
where λ 1 , λ 2 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 λ 1 and λ 2 are selected in this way
λ 1 = λ 2 = λ n (25)
where λ n is the arbitrary parameter, then the electric field will be represented like this
N = λ n ( U + E ) = λ n 3 e ( 2 z 2 − x 2 − y 2 ) 4 π ε r 5 r = λ n e 8 π ε 〈 1 + 3 cos 2 ϑ r 2 , 0 , 0 〉 (26)
The vector flux N in (30) through the sphere with the center at the origin of the coordinates
∬ S N r d S = λ N e 8 π ε ∫ 0 π d ϑ ∫ 0 2π ( 1 + 3 cos 2 ϑ ) d φ = 0 (27)
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 λ n describes the neutron electric field (see the details in Section 4), therefore, I will refer to the value N in (26) as a neutron component, and the parameter λ n -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
T = E + N (28)
Indeed, the flux of the vector ε T 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 [
W = ε 2 ∫ V ρ | T | 2 d V = ε 2 ∫ V ρ | E + λ n N | 2 d V (29)
Using (6) and (21), we obtain:
W = e 2 32 π 2 ε ∫ V ρ ( E r + λ n N r ) 2 d V (30)
The lower index ρ by the integration limit V ρ 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 V ρ 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
W P = e 2 32 π 2 ε ∫ V ρ E r 2 d V = e 2 32 π 2 ε ∫ ρ ∞ d r r 2 ∫ 0 2π d φ ∫ 0 π sin ϑ d ϑ = e 2 8 π ε ρ (31)
is the self-energy for the classical potential electron component. The second term
W N = λ n 2 e 2 32 π 2 ε ∫ V ρ n N r 2 d V = λ n 2 e 2 128 π 2 ε ∫ ρ n ∞ d r r 2 ∫ 0 2π d φ ∫ 0 π ( 1 + 3 cos 2 ϑ ) 2 sin ϑ d ϑ = λ n 2 e 2 10 π ε ρ n (32)
is the self-energy for the non potential component. The third term
W I = λ n e 2 16 π 2 ε ∫ V ρ E r N r d V = λ n e 2 16 π 2 ε ∫ ρ ∞ r d r ∫ 0 2π d φ ∫ 0 π 1 + 3 cos 2 ϑ r 2 sin ϑ d ϑ = 0 (33)
is the interaction energy for the potential and nonpotential components. Here W P is the potential component energy, W N is the nonpotential component energy, W I is the energy of the interaction between the potential and nonpotential components. Since W I = 0 , 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 m c 2 , where m is the particle mass and c 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 m c 2 one obtains
W P + W N = e 2 40 π ε ρ ( 5 + 4 λ n 2 ) = m c 2 (34)
The value of the multiplier λ n is unknown. To estimate the truncation radius of a charged particle, one can assume λ n = 0 , i.e. estimating the radius by the potential component one derives, as a result
ρ = e 2 8 π ε m c 2 (35)
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 m p = 1.6726218 × 10 − 27 kg [
ρ n λ n 2 = e 2 10 π ε m n c 2 = 6.13003 × 10 − 19 m (36)
Here, on the left-hand side of the equality there is unknown parameter λ n , m n 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 λ n and ρ n . 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) A = 0 . In spherical coordinates, the gradient of the second term in the right side of (18)
∇ [ A ln tan ϑ 2 ] = A r 2 〈 − ln tan ϑ 2 , 1 sin ϑ , 0 〉 (37)
Then, the interaction energy for the classical component (6) and (37) is written as
− A e 4 π ε ∫ ρ ∞ d r r 2 ∫ 0 π ln tan ϑ 2 sin ϑ d ϑ ∫ 0 2π d φ = 0 (38)
The interaction energy for the nonpotential component (21) and (37) is zero, too. Hence, the minimum energy (29) is reached with A = 0 .
The formula for the neutron electric field (26) contains the arbitrary constant λ n . There are two possible variants: λ n > 0 and λ n < 0 . 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 λ n 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 λ n < 0 and with the neutron λ n > 0 .
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 x ′ and y ′ . The second rotation is made about the y ′ axis through the angle β and the x ′ axis assumes the position x ″ . The third rotation is made about the x ″ 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 α = 0 . The matrices for the rotations about the y and x axes are written as [
A y = ( cos β 0 sin β 0 1 0 − sin β 0 cos β ) , A x = ( 1 0 0 0 cos γ − sin γ 0 sin γ cos γ ) (39)
The vectors are transformed using the matrices A = A x A y and the coordinates using the transposed matrix A ′ . Transformation of the vector N in (26) gives
N = λ n e 16 π ε q r 5 r (40)
where
q = 6 ( z 2 − x 2 ) cos 2 β cos 2 γ + ( x 2 − 2 y 2 + z 2 ) ( 3 cos γ − 1 ) + 12 x z cos 2 γ sin 2 β − 12 y z cos β sin 2 γ − 12 x y sin β sin 2 γ (41)
Because of the stress field symmetry, the angles are variable as follows
0 ≤ β < π , 0 ≤ γ < π (42)
Let us denote the external homogeneous electric field E ˜ = { E ˜ x , E ˜ y , E ˜ z } by the tilde. It has an arbitrary value and orientation, and in the spherical coordinates its radial component will be written as
E ˜ r = E ˜ x sin ϑ cos φ + E ˜ y sin ϑ sin φ + E ˜ z cos ϑ (43)
The other spherical components are of no interest, since only the radial component of the neutron component is other than zero. Here { E ˜ x , E ˜ y , E ˜ z } = const is the external field in the Cartesian coordinates. The interaction energy of the neutron component with the external field is
ε ∫ V ρ ( N ⋅ E ˜ ) d V = 0 (44)
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 [
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 x = a there is a particle with a pure potential field (6), and at the point x = − a there is a particle with an arbitrary oriented neutron field (40), then between them there is interaction with the energy
W e n = ε ∫ V ′ ρ ( E − ⋅ N + ) d V (45)
The lower indices − and + denote that the center of the given component is located on the axis of abscissa at the points − a and + a , respectively. Integration is performed over the whole space with the exception of two spherical regions with the radii ρ at the point x = + a and
k ρ (46)
at the point x = − a . 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 x = − a is
E − = e 4 π ε r − 3 r − (47)
and the neutron field at the point x = a is
N + = λ n e 16 π ε q + r + 5 r + (48)
where
q + = 6 [ z 2 − ( x − a ) 2 ] cos 2 β cos 2 γ + [ ( x − a ) 2 − 2 y 2 + z 2 ] ( 3 cos γ − 1 ) + 12 ( x − a ) z cos 2 γ sin 2 β − 12 y z cos β sin 2 γ − 12 ( x − a ) y sin β sin 2 γ (49)
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
x = a sinh τ cosh τ − cos σ , y = a sin σ cos φ cosh τ − cos σ , z = a sin σ sin φ cosh τ − cos σ (50)
The Jacobian is
Y = a 3 sin σ ( cosh τ − cos σ ) 3 (51)
The field (48) is derived from (40) by the transformation of the shift along the axis of abscissa x → x − a , and the field (47) from (6) by the shift to the opposite side x → x + a . Here r + and r − are
r + = ( x − a ) 2 + y 2 + z 2 = a 2 cosh − sinh τ cosh τ − cos σ (52)
r + = { x − a , y , z } (53)
r − = ( x + a ) 2 + y 2 + z 2 = a 2 cosh + sinh τ cosh τ − cos σ (54)
r − = { x + a , y , z } (55)
As a result, one obtains the interaction energy value (45) in the form
W ′ e n = ε ∫ − ( q + Q ) Q d τ ∫ 0 π d ϑ ∫ 0 2π ( E − ⋅ N + ) Y d φ = λ n e 2 3072 π ε a ( 1 + 3 cos 2 β − 6 cos 2 γ sin 2 β ) exp [ − 3 ( q + Q ) ] ⋅ { − 16 exp [ 3 ( q + Q ) ] ln [ 1 − exp 2 Q ] + 16 exp [ 3 ( q + Q ) ] ln ( 1 − exp [ − 2 ( q + Q ) ] ) − 5 [ 2 exp ( q + Q ) [ exp ( 2 q + 4 Q ) − 1 ] + ln coth q + Q 2 ] + exp [ 2 ( q + Q ) ] [ exp q ( 5 exp 4 Q − 15 − 6 exp 2 Q ) × ln coth Q 2 + 3 [ 2 + 5 exp 2 ( q + Q ) ] ln coth q + Q 2 ] } (56)
Here
Q = ln 2 a ρ , q = − ln k (57)
k is the coefficient from (46). The condition of ρ ≪ a is necessary to make the calculations correctly. This condition is fulfilled for real variants for which the value
Q ~ 10 ÷ 20 (58)
Using the above value Q as well as an approximate representation for small values of x
l n ( 1 + x ) ≈ x (59)
it is possible to simplify the Formula (56)
W ′ e n = λ n e 2 1536 π ε a ( 1 + 3 cos 2 β − 6 cos 2 γ sin 2 β ) ( 9 − 16 Q ) (60)
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 x = − a and a potential charge at the point x = + a , 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
W ″ e n = ε ∫ − Q Q d τ ∫ 0 π d ϑ ∫ 0 2π ( E + ⋅ N − ) Y d φ = λ n e 2 1536 π ε a ( 1 + 3 cos 2 β − − 6 sin 2 β − ) ( 9 − 16 Q ) (61)
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.
It is always possible to choose the coordinate orientation in such a way that the symmetric axis of the charge located at the point x = − a would lie in the xz plane. Due to the choice of the above orientation of the coordinate system, the angle γ − = 0 , 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 x = − a is derived from (48) by the replacement
a → − a , β + → β − , γ + → 0 (62)
As a result, one obtains the expression for the neutron at the point x = − a :
N − = λ n e 8 π ε r − 5 { ( x + a ) 2 − 2 y 2 + z 2 + 3 [ z 2 − ( x + a ) 2 ] c o s 2 β − + 6 ( x + a ) z s i n 2 β − } r − (63)
The interaction energy of the neutron-neutron will be written as
W n n = ε ∫ V ′ ρ ( N + ⋅ N − ) d V = ε ∫ − Q Q d τ ∫ 0 π d σ ∫ 0 2π ( N + ⋅ N − ) Y d φ = ± λ N 2 e 2 4056 π ε a sinh Q { 336 ( 3 cos 2 γ − 1 ) cosh Q − 48 [ ( 2 cos 2 β − 7 cos 2 ( β − β − ) ) − 39 cos 2 ( β + β − ) ] cos 2 γ + cos 2 β − ( 3 cos 2 γ − 1 ) cosh Q + 4 ln 1 + cosh Q 1 − cosh Q [ − 6 cos 2 β ( 5 + 13 cos 2 β − ) × cos 2 γ − ( 23 + 15 cos 2 β − ) ( 3 cos 2 γ − 1 ) − 3 ( − 7 − 2 cos 2 β cos 2 γ + cos 2 β − ( 1 + 46 cos 2 β cos 2 γ − 3 cos 2 γ ) + 21 cos 2 γ ) cosh 2 Q + 768 cos β cos β − cos 2 γ sin β sin β − sinh 2 Q ] } (64)
Neutrons interact with each other due to the nonpotential electric field. The energy (64) depends both on the charge signs λ n 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 Q is simplified
W n n = ± λ n 2 e 2 512 π ε a { 2 − 18 cos 2 β − 33 cos 2 ( β − β − ) + 18 cos 2 β − + 63 cos 2 ( β + β − ) − 3 cos 2 γ [ 2 + 6 cos 2 β + 2 ( 9 − 5 cos 2 β ) cos 2 β − + 32 sin 2 β sin 2 β − ] } (65)
For instance, for the parameters β = β − = γ = 0 the interaction energy value will be such
W n n = ± λ n 2 9 e 2 512 π ε a (66)
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
W = W e e + W n n + W ′ e n + W ″ e n = e 2 1536 π ε a { 192 + λ n ( 16 Q − 9 + 6 λ n ) + 6 λ n cos 2 β − ( 9 − 16 Q + 9 λ n + 30 λ n cos 2 β cos 2 γ − 27 λ n cos 2 γ ) − 3 λ n [ 2 ( 16 Q − 9 + 18 λ n ) cos 2 β cos 2 γ + ( 9 − 16 Q + 6 λ n ) cos 2 γ + 192 λ n cos 2 γ sin 2 β sin 2 β − ] } (67)
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.
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 [
μ 4 π { 0 , 0 , χ r } (68)
μ 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
M = μ 4 π ∇ × { 0 , 0 , χ r } = μ χ 4 π r 3 { − y , x , 0 } = μ χ 4 π r 2 〈 0 , 0 , sin ϑ 〉 (69)
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
μ p 4 π { 0 , 0 , z r 3 } (70)
Magnetic dipoles really exist, for instance, the magnetic electron moment. The rotor from (70) is the magnetic field induction produced by the magnetic dipole.
B = μ p 4 π ∇ × { 0 , 0 , z r 3 } = μ 3 p 4 π z r 5 { − y , x , 0 } = 3 μ p 8 π r 3 〈 0 , 0 , sin 2 ϑ 〉 (71)
In the spherical coordinates only the azimuthal component is other than zero, it decreases with a distance as 1 / r 3 . The field is symmetric relative to the plane z = 0 , 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:
W B = 1 2 μ ∭ V ρ ( B ⋅ B ) d V = μ p 2 20 π ρ 3 (72)
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
H + = z r 5 { − y ( 1 + ζ ) , x ,0 } = 1 4 r 2 〈 − 2 ζ s i n 2 φ s i n 2 ϑ , − ζ s i n 2 ϑ s i n 2 φ , 2 ( 2 + ζ − ζ c o s 2 φ ) s i n ϑ 〉 (73)
Here ζ is the constant which characterizes the deviation of the magnetic field from the potential one. At ζ = 0 the field is potential. The vector (73), i.e. the strength of the field produced by it, will be as follows
∇ × H + = 1 r 5 { 3 x z ,3 ( 1 + ζ ) y z , x 2 ( ζ − 1 ) + z 2 ( 2 + ζ ) − y 2 ( 1 + 2 ζ ) } (74)
The divergence (73) will be as follows
∇ ⋅ H + = ζ 3 ζ sin 2 ϑ sin 2 φ 2 r 3 (75)
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
∇ ⋅ Π = 3 ζ sin 2 ϑ sin 2 φ 2 r 3 (76)
Then
∇ ⋅ ( H + − Π ) = 0 (77)
Since the vector Π is potential, it can be represented as the gradient of a certain function Ψ :
Π = ∇ Ψ (78)
Substituting this relation into (76) one obtains
Δ Ψ = 3 ζ sin 2 ϑ sin 2 φ 2 r 3 (79)
In the spherical coordinate system the equation is written as [
1 r 2 ∂ ∂ r ( r 2 ∂ Ψ ∂ r ) + 1 r 2 sin ϑ ∂ ∂ ϑ ( sin ϑ ∂ Ψ ∂ ϑ ) + 1 r 2 sin 2 ϑ ∂ 2 Ψ ∂ 2 φ 2 = 3 ζ sin 2 ϑ sin 2 φ 2 r 3 (80)
As in Section 0, let us assume that the unknown function Ψ can be repre- sented in the form of the product
Ψ ( r , ϑ , φ ) = R ( r ) Θ ( ϑ ) Φ ( φ ) (81)
Assuming that R ( r ) = 1 / r one gets the equation
sin ϑ Θ ( ϑ ) ∂ ∂ ϑ ( sin ϑ ∂ Θ ( ϑ ) ∂ ϑ ) + 1 Φ ( φ ) ∂ 2 Φ ( φ ) ∂ φ 2 = ζ 3 sin 4 ϑ sin 2 φ 2 Θ ( ϑ ) Φ ( φ ) (82)
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
sin 2 φ Φ ( φ ) = C 0 (83)
Here the constant ζ is included into the arbitrary constant C 0 , then
1 Φ ( φ ) ∂ 2 Φ ∂ φ 2 = − 4 (84)
and the Equation (82) goes into
sin ϑ ∂ ∂ ϑ ( sin ϑ ∂ Θ ( ϑ ) ∂ ϑ ) − 4 Θ ( ϑ ) = 3 2 C 0 sin 4 ϑ (85)
Let us denote
t = ln tan ϑ 2 (86)
therefore,
ϑ = 2 arctan exp t (87)
Then using the differentiating rules
∂ ∂ ϑ = ∂ ∂ t ∂ t ∂ ϑ = 1 sin ϑ ∂ ∂ t (88)
so
sin ϑ ∂ ∂ ϑ = ∂ ∂ t (89)
and (85) goes into
∂ 2 Θ ∂ t 2 − 4 Θ = 3 C 0 2 cosh 4 t (90)
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 [
Integrating them one obtains
The solution of the equation (90) will be written as
Here
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
includes two arbitrary constants
produced by this vector field will be an unlimited value at
will be unlimited too. During integration (103) first, let us find a certain integral with respect to
In order to avoid the divergence with respect to
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
For the electric field of the neutron
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
Here
is the neutron dipole electric moment. Nowadays this value has not been measured yet, although there have been some attempts to do it [
Here
Hence it follows that
The coefficient of proportionality
Here it is denoted that
Modern experimental methods [
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
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 [
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
and the energy of the proton-neutron interaction (114) is written as
Using the relation (58) one derives that for
Ivanchin, A. (2017) Electrostatic Theory of Elementary Particles. World Journal of Nuclear Science and Technology, 7, 232-251. https://doi.org/10.4236/wjnst.2017.74019