In this work we re-examine a model of the nucleons that involve the weak interaction which was once considered by Heisenberg; that is a neutron may have the structure of a dwarf hydrogen-like atom. We formulate a quantum dynamics for the Heisenberg model of the neutron associated with interaction that involves the beta decay in terms of a mixed Coulomb-Yukawa potential and the More General Exponential Screened Coulomb Potential (MGESCP), which has been studied and applied to various fields of physics. We show that all the components that form the MGESCP potential can be derived from a general system of linear first order partial differential equations similar to Dirac relativistic equation in quantum mechanics. There are many interesting features that emerge from the MGESCP potential, such as the MGESCP potential can be reduced to the potential that has been proposed to describe the interaction between the quarks for strong force in particle physics, and the energy spectrum of the bound states of the dwarf hydrogen-like atom is continuous with respect to distance. This result leads to an unexpected implication that a proton and an electron may also interact strongly at short distances. We also show that the Yukawa potential when restrained can generate and determine the mathematical structures of fundamental particles associated with the strong and weak fields.
Despite that the mathematical formulation of quantum mechanics has been highly developed and the theory has been successfully applied into all domains of applied sciences with the most accuracies that can be achieved by experiments, many fundamental physical processes at the quantum level that involve quantum mechanics still remain a mystery. In particular, one of the profound epistemological problems that continue to exist is the question of whether microscopic phenomena are in fact continuous or progressing in quantum jumps. In an article entitled ARE THERE QUANTUM JUMPS? Schrödinger wrote: “… A great many of our educated contemporaries, not equipped with the mathematical apparatus to follow our more technical deliveries, are yet deeply concerned with many general questions; one of the most stirring among them certainly is whether actually natura facit saltus or no …” [
Since the energy spectrum of the emitted electron in the beta minus decay is continuous therefore Heisenberg’s model of the neutron as a dwarf hydrogen-like atom cannot be realised if we only apply the Coulomb potential to describe the system. Instead, in this work we will show that a continuous spectrum of energy can be obtained by applying a mixed Coulomb-Yukawa potential of the form
V ( r ) = − α e − β r r + Q r , (1)
where Q , α and β are physical parameters that will need to be determined [
V ( r ) = − V 0 e − 2 α r r − V 0 r − V 0 α e − 2 α r , (2)
where V 0 is the potential depth and the parameter α is the strength coupling constant [
In the present state of physics there are four confirmed types of interactions between physical objects, which are the gravitational interaction, the electromagnetic interaction, the strong interaction, and the weak interaction. Except for the gravitational interaction, the other three types of interactions can be mathematically formulated so that they can comply with quantum mechanics, especially in the so-called standard model of particle physics [
∑ i = 1 n ∑ j = 1 n a i j r ∂ ψ i ∂ x j = ∑ i = 1 n ( ∑ j = 1 n b i j r V j + c i r ) ψ i + d r , r = 1 , 2 , ⋯ , n (3)
The system of equations given in Equation (3) can be rewritten in a matrix form as
( ∑ i = 1 n A i ∂ ∂ x i ) ψ = − i ( ∑ i = 1 n q B i V i + m σ ) ψ + i k J , (4)
where ψ = ( ψ 1 , ψ 2 , ⋯ , ψ n ) T , ∂ ψ / ∂ x i = ( ∂ ψ 1 / ∂ x i , ∂ ψ 2 / ∂ x i , ⋯ , ∂ ψ n / ∂ x i ) T with A i , B i and σ are matrices representing the quantities a i j r , b i j r , c j r , which are assumed to be constant in this work, and k J is a matrix that represents the quantity d r , where k is a dimensional constant. As normal convention, the letter i in front of the last two terms in Equation (4) is the imaginary number i. While the quantities q, m and k J represent physical entities related directly to the physical properties of the particle under consideration, the quantities V i represent the potentials of an external field, such as an electromagnetic field or the matter field of a quantum particle. In order to formulate a physical theory from the system of equations given in Equation (4), it is necessary to determine the unknown quantities A i , B i and σ , as well as the mathematical conditions that they must satisfy, such as commutation relations between them. The commutation relations between the matrices can be determined if we apply the operator ∑ i = 1 n A i ∂ / ∂ x i on the left on both sides of Equation (4) as follows
Since the quantities A i , B i , σ , q and m are assumed to be constant, Equation (5) becomes
( ∑ i = 1 n A i 2 ∂ 2 ∂ x i 2 + ∑ i = 1 n ∑ j > i n ( A i A j + A j A i ) ∂ 2 ∂ x i ∂ x j ) ψ = ( − i ( ∑ i = 1 n A i ∂ ∂ x i ) ( ∑ j = 1 n q B j V j + m σ ) ) ψ − i ( ∑ i = 1 n q B i V i + m σ ) ( ( ∑ j = 1 n A j ∂ ∂ x j ) ψ ) + ∑ i = 1 n A i ∂ ( i k J ) ∂ x i (6)
= − i ( ∑ i = 1 n ∑ j = 1 n q A i B j ∂ V j ∂ x i ) ψ − ( ∑ i = 1 n ∑ j > i n q 2 ( B i B j + B j B i ) V i V j − 2 i ∑ i = 1 n q m B i V i σ − m 2 σ 2 ) ψ − i ( ∑ i = 1 n q B i V i + m σ ) ( i k J ) + ∑ i = 1 n A i ∂ ( i k J ) ∂ x i
In order to describe the dynamics of a particular physical system, undetermined parameters given in Equation (4) must be specified accordingly. To obtain a system of partial differential equations similar to Dirac equation for an arbitrary field ( γ μ ( i ∂ μ − q V μ ) − m ) ψ = 0 [
( ∑ i = 1 4 γ i ∂ ∂ x i ) ψ = − i ( ∑ i = 1 4 q γ i V i + m ) ψ + i k J , (7)
where the matrices A i and B i have been identified with Dirac matrices γ i as follows
γ 1 = ( 1 0 0 0 0 1 0 0 0 0 − 1 0 0 0 0 − 1 ) , γ 2 = ( 0 0 0 1 0 0 1 0 0 − 1 0 0 − 1 0 0 0 ) , γ 3 = ( 0 0 0 − i 0 0 i 0 0 i 0 0 − i 0 0 0 ) , γ 4 = ( 0 0 1 0 0 0 0 − 1 − 1 0 0 0 0 1 0 0 ) (8)
For the case of Dirac equation, the matrices γ i and σ satisfy the following conditions
γ i 2 = ± 1 (9)
γ i γ j + γ j γ i = 0 for i ≠ j (10)
It is seen from Equation (4) that the quantity k J represents an internal source which is associated with a dynamical process of a quantum particle. For example, a quantum particle is undergoing some form of physical evolution that changes its physical structure, such as an accumulation of mass during its formation. In fact we will show that this physical process can be formulated using the MGESCP potential in which the energy spectrum depends continuously on distance. In terms of the gamma matrices γ i then Equation (7) can be rewritten in a covariant form similar to Dirac equation for an arbitrary field with an internal source k J as
( γ μ ( i ∂ μ − q V μ ) − m ) ψ = i k J . (11)
In this case Equation (6) also reduces to the following equation
If the quantities V i are the four-potential of an electromagnetic field given by the identification ( V 1 , V 2 , V 3 , V 4 ) = ( V , A x , A y , A z ) then Equation (12) can be used to determine the dynamics of the components of the wavefunction ψ = ( ψ 1 , ψ 2 , ψ 3 , ψ 3 ) T , where the term ∂ V j / ∂ x i − ∂ V i / ∂ x j are the components of the electric field E and the magnetic field B .
Now we will discuss how free quantum particles can create their own physical fields in which wavefunctions can be identified as potentials. Therefore we set ∑ i = 1 n q B i V i = 0 . Equations (7) and (12) for free particles reduce to equations
γ μ ∂ μ ψ = − i m ψ + i k J , (13)
γ i 2 ∂ 2 ψ ∂ x i 2 = − m 2 ψ + m k J + γ i ∂ ( i k J ) ∂ x i . (14)
In the following we will consider two cases depending on the conditions that are applied to the internal source k J in which either k J = 0 or ≠ 0 . For the case k J = 0 , Equation (13) reduces to a form similar to Dirac equation for a free particle
γ μ ∂ μ ψ = − i m ψ . (15)
For massive particles in which m ≠ 0 , all components of the wavefunction ψ μ satisfy the Klein-Gordon equation
∂ 2 ψ μ ∂ t 2 − ∂ 2 ψ μ ∂ x 2 − ∂ 2 ψ μ ∂ y 2 − ∂ 2 ψ μ ∂ z 2 = − m 2 ψ μ . (16)
And, in particular, if the wavefunctions are time-independent then we obtain
∂ 2 ψ μ ∂ x 2 + ∂ 2 ψ μ ∂ y 2 + ∂ 2 ψ μ ∂ z 2 = m 2 ψ μ . (17)
In this case the wavefunctions ψ μ can be viewed as static Yukawa potential
ψ μ ( r ) = − α e − β r r , (18)
where α and β are undetermined dimensional constants [
For massless time-independent particles, the Klein-Gordon equation given in Equation (17) reduces to Laplace equation
∂ 2 ψ μ ∂ x 2 + ∂ 2 ψ μ ∂ y 2 + ∂ 2 ψ μ ∂ z 2 = 0. (19)
Solutions to Laplace equation can be written in the form
ψ μ ( x , y , z ) = Q r . (20)
In this case the wavefunctions ψ μ can be viewed as static Coulomb potential, where Q is an undetermined dimensional constant, which is associated with the electromagnetic interaction between elementary particles.
As mentioned in the introduction, we will discuss possible bound states of a dwarf hydrogen-like atom which can be identified with a neutron therefore we will need to use a more general form of Yukawa potential, which is the MGESCP potential given as in Equation (2). Since the MGESCP potential has an extra term of the form ψ μ = a e − m r , therefore we now need to show how to derive this form of potential from the Dirac-like equation with an internal source given in Equation (13). Now, Dirac-like wavefunctions ψ μ satisfy the following Klein-Gordon-like equation with a source
∂ 2 ψ μ ∂ t 2 − ∂ 2 ψ μ ∂ x 2 − ∂ 2 ψ μ ∂ y 2 − ∂ 2 ψ μ ∂ z 2 = − m 2 ψ μ + m k J + γ i ∂ ( i k J ) ∂ x i . (21)
In particular, if the wavefunctions are time-independent then we obtain
∂ 2 ψ μ ∂ x 2 + ∂ 2 ψ μ ∂ y 2 + ∂ 2 ψ μ ∂ z 2 = m 2 ψ μ − m k J − γ i ∂ ( i k J ) ∂ x i . (22)
It can be verified that a solution of the form ψ μ = a e − m r , where a and m are constants, satisfies the following equation
∂ 2 ψ μ ∂ x 2 + ∂ 2 ψ μ ∂ y 2 + ∂ 2 ψ μ ∂ z 2 = m 2 ψ μ − 2 m a e − m r r . (23)
By comparing Equation (23) to Equation (22), we obtain the following equation for the internal quantity k J
m k J + γ i ∂ ( i k J ) ∂ x i = 2 a m e − m r r . (24)
A differential equation for the quantity k J can be determined by using the matrices γ i given in Equation (8), which is written in an explicit form as
( 0 0 0 1 0 0 1 0 0 − 1 0 0 − 1 0 0 0 ) ∂ ( i k J ) ∂ x + ( 0 0 0 − i 0 0 i 0 0 i 0 0 − i 0 0 0 ) ∂ ( i k J ) ∂ y + ( 0 0 1 0 0 0 0 − 1 − 1 0 0 0 0 1 0 0 ) ∂ ( i k J ) ∂ z = ( 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ) ( − m k J + 2 a m e − m r r ) . (25)
From Equation (25) we obtain the following equations for the quantity k J
− m k J + 2 a m e − m r r = 0 , (26)
∂ ( k J ) ∂ z = 0 , ∂ ( k J ) ∂ x − i ∂ ( k J ) ∂ y = 0 , ∂ ( k J ) ∂ x + i ∂ ( k J ) ∂ y = 0. (27)
The equations given in Equation (27) show that the source k J is constant and from Equation (26) this also results in the constancy of the Yukawa potential which can be written as
e − m r r = k J 2 a . (28)
Now we apply the results that have been obtained into the MGESCP potential given in Equation (2), one of whose components has the form ψ μ = V 0 α e − 2 α r . By comparing this potential to ψ μ = a e − m r we have a = V 0 α and m = 2 α . With the Yukawa potential that is restrained by the condition given in Equation (28) the MGESP potential given in Equation (2) is reduced to
V ( r ) = − V 0 r ( 1 + ( 1 + α r ) e − 2 α r ) = − V 0 r + k J 2 α ( 1 + α r ) = − V 0 r + k J 2 r + k J 2 α . (29)
Except for the constant k J / 2 α , the potential given in Equation (29) has the form that is similar to the potential that describes an interaction between two fundamental quarks as proposed in the theory of quantum chromodynamics, namely, V ( r ) = A / r + B r . This type of potential describes interactions between two quarks that can be represented in the following picture shown in
For small values of the distance r the potential manifests as a Coulomb potential V = A / r , however, for large values of distance the potential acts as a linear potential with respect to the distance V = B r . The linear potential shown in
In this section we discuss further the restraint to the Yukawa potential given in Equation (28) which has been shown to reduce the MGESCP potential to the potential that is proposed for the interaction between the quarks for strong force in particle physics. We now show that in fact the restrained Yukawa potential actually generates and determines mathematical structures of physical objects that may be identified with quantum mediators of the weak and strong interactions. Instead of giving a mathematical analysis of the restrained Yukawa potential given in Equation (28), as an illustration, we simply plot possible shapes that can be generated and determined by a restrained Yukawa potential from the relation given in Equation (28), namely, e − m r = ( k J / 2 a ) r , with different values given to the parameters m and k J / 2 a . Together, they possess a remarkable difference in their topological structures that may underlie physical effects that are associated with elementary quantum particles [
For the case m = 0.001 with k J / 2 a = γ = 1 , 2 , 3 , 4 we have the following possible shapes for elementary quantum particles as shown in Figures 3-6.
For the case m = 0.01 with k J / 2 a = γ = 1 , 2 , 3 , 4 we have the following possible shapes for elementary quantum particles as shown in Figures 7-10.
For the case m = 0.1 with k J / 2 a = γ = 0.1 , 1 , 2 , 3 we have the following possible shapes for elementary quantum particles as shown in Figures 11-14.
For the case m = 1 with k J / 2 a = γ = 0.001 , 0.01 , 0.1 , 1 we have the following possible shapes for elementary quantum particles as shown in Figures 15-18.
For the case m = 5 with k J / 2 a = γ = 10 − 11 , 10 − 3 , 0.01 , 0.1 we have the following possible shapes for elementary quantum particles as shown in Figures 19-22.
For the case m = 10 with k J / 2 a = γ = 10 − 26 , 10 − 21 , 10 − 18 , 10 − 13 we have the following possible shapes for elementary quantum particles as shown in Figures 23-26.
In this section we will discuss whether a neutron can be modelled as a dwarf hydrogen-like atom with the two different mixed potentials given in Equations (1) and (2) so that it can be used to explain the physical processes associated with the beta minus decay. As shown in Section 2, these potentials can be formed from the three forms of potentials that have been derived from the Dirac equations. First we consider the mixed potential that is formed from the Coulomb potential and the Yukawa potential as given in Equation (1). As shown below, this form of potential can be used to explain how an electron can be repelled from a dwarf hydrogen-like atom composed of a proton and an electron. By differentiating Equation (1), we can obtain the following equations
d V d r = α e − β r [ β r + 1 r 2 ] − Q r 2 , d 2 V d r 2 = − α e − β r [ β 2 r + 2 β r 2 + 2 r 3 ] + 2 Q r 3 . (30)
From Equation (30), the corresponding force of interaction F ( r ) = − d V / d r is obtained as
F ( r ) = − α e − β r [ β r + 1 r 2 ] + Q r 2 . (31)
According to classical electrodynamics, the net force acting on the electron must be zero when it is circulating in stable orbits. If we assume the net force acting on the electron to vanish when it moves in a stationary orbit of finite radius r = R = 1 / β , i.e. F ( 1 / β ) = 0 , then from Equation (31) we obtain the relation
β 2 ( Q − 2 α e ) = 0. (32)
Since β ≠ 0 , we require Q = 2 α / e . The mixed potential given in Equation (1) now takes the form
V ( r ) = − e Q 2 e − r R r + Q r . (33)
And the corresponding force of interaction F ( r ) = − d V / d r is
F ( r ) = − e Q 2 e − r R [ 1 R r + 1 r 2 ] + Q r 2 . (34)
In order to investigate further we need to know the nature of the stationary point at r = R . From Equation (30), the second derivative of V ( r ) at r = R is found as
d 2 V d r 2 = − Q 2 R 3 . (35)
Since we are considering the case for which the mixed potential given in Equation (33) is applied to the bound system of two charges of opposite signs, namely a proton and an electron, therefore the quantity Q can be written as Q = − K q 1 q 2 , where K > 0 is a coupling constant. If q 1 is the charge of the proton, q 1 = q , and q 2 is the charge of the electron, q 2 = − q , then we have Q = K q 2 > 0 . Then from Equation (35) we obtain the condition d 2 V / d r 2 < 0 , therefore V ( r ) has a local maximum at r = R . Since F ( r ) = − d V / d r , the force is attractive for r < R and repulsive for r > R . This situation is seen similar to the case of weak interaction of beta minus decay in which a neutron turns into a proton and emits an electron and an anti-neutrino. First the neutron turns into a dwarf hydrogen-like atom, whose bound states will be described below using the MGESCP potential, then the electron moves in an orbit of zero force and then it is repelled from the dwarf hydrogen atom by a repulsive force. The force given in Equation (34) is rewritten as follows
F ( r ) = − e K q 2 2 e − r R [ 1 R r + 1 r 2 ] + K q 2 r 2 . (36)
The force given in Equation (36) is assumed to be a weak force. Since R is the radius of the stationary orbit therefore we can assume that the electron is ejected from the stationary orbit because the equilibrium of this system at r = R is unstable.
It can be seen that the whole process of beta decay is a complicated physical process that actually undergoes through many different physical configurations of the system, therefore it can only be described approximately by many different dynamics, only if we can formulate the whole physical process under a mathematical formulation that can give rise to each physical configuration by some form of limit associated with mathematical parameters that are used to describe the whole system. With this in mind in the following we will discuss in terms of Schrödinger wave mechanics a more complete structure of a neutron as a dwarf hydrogen-like atom using a more complete MGESCP potential given in Equation (2). Consider the time-independent Schrödinger equation
− ℏ 2 2 μ ∇ 2 ψ ( r ) + V ( r ) ψ ( r ) = E ψ ( r ) , (37)
in which V ( r ) is the More General Exponential Screened Coulomb Potential (MGESCP) given in Equation (2). Since the MGESCP potential is spherically symmetric, Equation (37) can be written in the spherical polar coordinates as
− ℏ 2 2 μ ( 1 r 2 ∂ ∂ r ( r 2 ∂ ∂ r ) − L 2 ℏ 2 r 2 ) ψ ( r ) + ( − V 0 r ( 1 + ( 1 + α r ) e − 2 α r ) ) ψ ( r ) = E ψ ( r ) , (38)
where the orbital angular momentum operator L 2 is given by
L 2 = − ℏ 2 ( 1 sin θ ∂ ∂ θ ( sin θ ∂ ∂ θ ) + 1 sin 2 θ ∂ 2 ∂ 2 ϕ ) . (39)
Solutions of Equation (38) can be found using the separable form
ψ n l ( r ) = R n l ( r ) Y l m ( θ , ϕ ) , (40)
where R n l is a radial function and Y l m is the spherical harmonic. Applying Equation (40), Equation (38) is reduced to the system of equations
L 2 Y l m ( θ , ϕ ) = l ( l + 1 ) ℏ 2 Y l m ( θ , ϕ ) , (41)
( − ℏ 2 2 μ ( d 2 d r 2 + 2 r d d r ) + l ( l + 1 ) ℏ 2 2 μ r 2 − V 0 r ( 1 + ( 1 + α r ) e − 2 α r ) ) R n l ( r ) = E R n l ( r ) . (42)
It has been shown that the radial solution R n l ( r ) to Equation (42) can be obtained as
R n l ( r ) = N n l r ( − 1 + 4 l ( l + 1 ) + 1 ) / 2 e − β r L n 1 + 4 l ( l + 1 ) + 1 ( 2 β r ) , (43)
and the corresponding energy spectrum E n l is given by
E n l = − V 0 α e − 2 α r − μ 2 ℏ 2 ( ( V 0 + V 0 e − 2 α r ) n + l + 1 ) 2 , (44)
where β 2 = ( 2 μ ( V 0 + V 0 e − 2 α r ) / ℏ 2 ( 2 n + 1 + 4 l ( l + 1 ) + 1 ) ) 2 [
E n l = − k J 2 r − μ 2 ℏ 2 ( ( V 0 + ( k J / 2 α ) r ) n + l + 1 ) 2 . (45)
Now we may interpret this continuous spectrum of energy with respect to distance as the energy spectrum of massive mediators associated with strong force described by the potential given in Equation (29). When a physical particle is created it is being created continuously until it reaches the size that is required for the system. This process happens in a very short time therefore it seems like an instantaneous creation. In particle physics, the parameter α of the exponential term is expressed in terms of the mass m of a force carrier as α = m c / 2 ℏ . Therefore when the mass of the force carrier is being continuously created the parameter α is being getting larger, at the same time the radius r is also getting bigger, therefore the term e − 2 α r → 0 and also the term α e − 2 α r → 0 . The mass that is accumulated by the force carrier must be supplied by the neutron. When the force carrier with required mass hits the electron, the latter will move further from the proton. On the other hand the MGESCP potential is reduced to the mixed Coulomb-Yukawa potential when the process of creation of the force carrier is complete. This form of potential provides a repulsive force to move the electron away.
As a further remark on the forms of potentials given in Equations (1) and (2), it is seen that the potential given in Equation (2) cannot be reduced to the potential given in Equation (1). And as a consequence the solutions given in Equation (43) cannot be reduced to solutions to a Schrödinger wave equation that uses a potential of the form given in Equation (1). In order to reduce to the potential given in Equation (1) we would need to consider a more general form of potential which can be written as
V ( r ) = − α e − β r r + Q r + K e − γ r . (46)
where K , Q , α , β and γ are physical parameters that will need to be determined. This potential would describe a more complete physical process of beta decay when it is applied to the Schrödinger wave equation given in Equation (37). However, whether the Schrödinger wave equation with this form of potential could be solved to obtain exact solutions similar to solutions given in Equation (43) requires more rigorous mathematical investigations.
In this work we have discussed a quantum dynamics of Heisenberg’s model of the neutron associated with the beta minus decay through the weak and strong interactions in which a neutron may have the structure of a dwarf hydrogen-like atom. It has been shown that the whole process of beta decay is a complicated physical process that actually undergoes many different physical states of configuration of the system. Therefore, the physical process of beta decay should be described by many different dynamics rather than a single one, only if we can formulate the whole physical process under a mathematical formulation that can give rise to each state by some form of limit associated with mathematical parameters that are used to describe the whole system. For example, a more complete mathematical formulation would require a more general form of potential as given in Equation (46) in Section 4. However, in this work we have still been able to discuss in terms of Schrödinger wave mechanics a quantum dynamics of the neutron as a dwarf hydrogen-like atom using a more complete potential which has been studied and applied to various fields of physics, the so-called the More General Exponential Screened Coulomb Potential (MGESCP). We have shown that the MGESCP potential can be derived from a Dirac-like system of equations which can be reduced from a general system of linear first-order partial differential equations. There is a particular advantage for our approach to deriving a Dirac-like system of equations from a general system of differential equations, that are Dirac-like wavefunctions which can be interpreted according to the purpose of the mathematical investigation of a physical system. They can be used to represent different type of physical objects rather than exclusive mathematical wavefunctions that are used to calculate the probability of the outcome of an experimental result. From the MGESCP potential we have also been able to derive a form of potential that is used to describe interaction between quarks in strong interaction and, in particular, the energy spectrum of the bound state of the dwarf hydrogen-like atom that is continuous with respect to space. This may be a manifestation of a continuous creation of an elementary particle in space. We have also shown that the Yukawa potential can generate and determine the physical shapes of fundamental particles associated with the strong and weak fields. As a consequence, it seems reasonable to suggest that the functional potential in physics may have physical mechanisms to generate mediators of associated physical fields, and these mechanisms can be formulated in terms of differentiable manifolds and their corresponding direct sums of prime manifolds as discussed in our works on the possibility to formulate physics in terms of differential geometry and topology.
We would like to thank the reviewers for their critical and constructive comments and we would also like to thank Jane Gao of the administration of JMP for her editorial advice during the preparation of this work.
The author declares no conflicts of interest regarding the publication of this paper.
Ho, V.B. (2019) A Quantum Dynamics of Heisenberg Model of the Neutron Associated with Beta Decay. Journal of Modern Physics, 10, 1065-1082. https://doi.org/10.4236/jmp.2019.109069