The deuteron, the nucleus of heavy hydrogen, is a very interesting, important special case. Its two nucleons, the proton and neutron, are held together by the nuclear force alone. This nuclear force is carried by the π meson. As with gravitation or the Coulomb force there is an inverse-square-of-the-distance dependence because force-carrying particles are spread over a spherical surface whose area is proportional to the square of the distance from the source.. However there is an essential difference for the nuclear force as compared to the electrostatic Coulomb force in that the particles carrying the electrostatic force, photons, do not decay whereas the π mesons of the nuclear force, a.k.a the strong force, do decay quite rapidly with distance. The proportion of the mesons which survive is an exponential function of distance. Therefore the formula for the nuclear force would be of the form

−H*e^−λ*d/d²
 

where d is the distance between the nucleons, H* and λ* are constants with dimensions of kg·m³/s² and m⁻¹, respectively. The material which follows is an attempt to work out the implications of this force formula for the quantum level phenomena of the deuteron.

At this point it is convenient to switch notation. Let r stand for the radius of the orbit of a particle with respect to the center of mass of the particle pair. For the electron in an proton-electron pair r is essentially equal to the distance d between the particles. For the proton-neutron pair of a deuteron r is essentially equal to one half of the distance between the particles. This means the formula for the force experienced by each nucleon as a function of the radius r of their orbits with respect to their center of mass is

H*e^−λ*2r/(2r)² = (H*/4)e^−(2λ*)r/r²
= He^−λr/r²
 

where H=H*/4 and λ=2λ*. It is conventient to express λ* as 1/d₀; i.e., d₀ = 1/λ*. This means that r₀ = 1/(2λ*) = ½d₀

Hideki Yukawa established that the mass of the force-carrying particles is related to the λ parameter. From the mass of the π mesons it is known that d₀ is about 1.54 fermi (1.54×10⁻¹⁵ meters). Thus r₀ is about 0.77 fermi.

Potential Functions

In mechanics the force is conveniently represented as the negative of the derivative of a potential function V(r). This means that the potential function is given by the integration of the force function with with respect to distance. For the Coulomb force of −α/r² this gives a potential function of V(r)=−α/r. Yukawa hypothesized a potential function which is just this potential function multiplied by an exponential factor; i.e., V(r)=(−α/r)e^−λr. However it is the force function that is multiplied by the exponential factor due to the decay of the force-carrying mesons; i.e.,

F = −H*e^−λd/d²
 

where d is the distance between the two nucleons and H* and λ* are constants characteristic of the nuclear force. Therefore the potential function for the nuclear force is really of the form

V(d) = ∫_d^∞(H*e^−λ*s/s²)ds
 

Quantum Mechanics

The standard procedure in quantum mechanics is to solve Schroedinger's equation for the particular potential function which applies. This works beautifully for the Coulombic potential and about a half dozen other special cases but gives no analytic results for the cases that deviate from those special cases. The Schroedinger equation cannot be solved even for the Yukawa potential function. There is no way to solve the Schroedinger equation analytically for even more complicated potential functions. Physicists then rely upon purely numerical methods of solution based upon perturbation theory.

What does work analytically is the analysis of Niels Bohr's Old Quantum Theory. Using Bohr's analysis it can be established that the angular momentum is quantized for any potential function in exactly the same way that it is for the Coulombic force. This means that angular momentum must be equal to an integral multiple of h, Planck's constant divided by 2π. This is also true whether the kinetic energy is of the Newtonian form ½mv², where m is the mass of the particle and v is its velocity or the relativistic form

m₀c²[1/(1−β²)^½−1]

where β is the velocity relative to the speed of light, v/c.

Circular Orbits

In a circular orbit the attractive force balances the centrifugal force; i.e.,

mv²/r = He^−r/r₀/r²
and hence
mv² = He^−r/r₀/r
 

This means that orbital velocity v is a function of orbit radius r; i.e.,

v = [He^−r/r₀/(mr)]^1/2
 

Angular momentum is defined as p_θ=mvr, which on the basis of the relation between v and r means

p_θ = [Hmre^−r/r₀]^1/2
 

But previous analysis found angular momentum to be quantized. That is to say,

p_θ = lh
 

where l is an integer and h is Planck's constant divided by 2π.

This means that

Hmre^−r/r₀ = l²h²
or, dividing by r₀ and solving
(r/r₀)e^−r/r₀ = l²h²/(Hmr₀)
 

Now let r/r₀ be denoted as z.

The function ze^−z rises from 0 at z=0 to a maximum of e⁻¹ at z=1 and falls asymptotically to 0 as z goes to +∞, as shown below.

For sufficiently small values of angular momentum there will be two values of r but there is a maximum angular momentum for which there are any solutions for r.

The levels shown in green above correspond to the quantity h²/(Hmr₀) multiplied by the square of an integer l. The value of this coefficient can be calculated. The value of h² in SI units is 1.1121×10⁻⁶⁸. The value of m is 1.67×10⁻²⁷ kg. As indicated above the value of r₀ is 7.7×10⁻¹⁶ m.

It has been known since 1935 that if a deuteron is hit by a sufficiently high energy photon the deuteron disintegrates into a free proton and a free neutron. The energy level required for this photodisintergration of the deuteron is 2.22457 Mev. This is the binding energy of the deuteron. This energy is 3.5638×10⁻¹³ joules. This is reasonably identified as the potential energy of each of the l=1 level nucleons. A photon that knocks one nucleon out of its orbit disassociates the deuteron.

Let r₁ be the radius of the nucleons in the deuteron. The value of r₁/r₀=2 is a significant figure. The accepted diameter of the deuteron is 4.2 fermis. The diameters of the proton and neutron are both about 1 fermi. Deducting 1 fermi, for the combined radii of the proton and neutron, from 4.2 fermis gives 3.2 fermis as the distance between the centers of mass of the two nucleons. This represents orbit radii of 1.6 fermis for both of the two nucleons. This figure divided by r₀=0.77 fermi is 2.078.

Before estimating the value of H it is advisable to consider the relativistic adjustment which should be made.

The Relativistic Case

Under the Special Theory of Relativity the kinetic energy is given by the expression

K(β) = m₀c²[(1−β²)^−½−1].
 

where β=v/c and m₀ is the rest mass of the particle.

The first two terms of this kinetic energy function are

K(v) = ½mv² + (3/8)mv⁴/c² + ...
 

Circular Orbits

The balance of the attractive nuclear strong force with the centrifugal force under Special Relativity is

mv²/r = He^−r/r₀/r²
but
m = m₀/(1-β²)^1/2
so
m₀v²/(r(1-β²)^1/2)) = He^−r/r₀/r²
 

This can be reduced to

β⁴/(1-β²)) = [He^−r/r₀/(m₀c²r)]²
 

Multiplying by the denominator on the left produces a quadratic equation in β²; i.e.,

β⁴ + β²ζ² − ζ² = 0
 

where ζ² = [He^−r/r₀/(m₀c²r)]². The solutions are

β² = ½[−ζ² ± (ζ⁴ + 4ζ²)^½]

 

The negative solution must be discarded. Thus

β² = ½[(ζ⁴ + 4ζ²)^½−ζ²]
or, equivalently
β² = ½[(ζ²(1 + 4/ζ²)^½−ζ²]
 

For large values of ζ the solution is approximately β=1.

Quantization

It is shown elsewhere that even in the relativistic case angular momentum p_θ=mvr is quantized in increments of h so

p_θ = lh
 

where l is an integer.

For a circular orbit

mv = He^−r/r₀/(vr)
 

Therefore

p_θ = [He^−r/r₀/(vr)]r
p_θ = He^−r/r₀/v
p_θ = He^−r/r₀/(cβ)
 

Since p_θ = lh it follows that

β = γe^−r/r₀/l
 

where γ=H/(hc). This is not an entirely satisfactory expression of quantization because r on the RHS is quantized in some as yet undetermined manner, yet the similarity with the non-relativistic case makes it of interest.

From the analysis of the previous section it is known that

β² = ½[(ζ²(1 + 4/ζ²)^½−ζ²]
or, equivalently
β² = ½[(1 + 4/ζ²)^½−1]ζ²
and hence
β = [½[(1 + 4/ζ²)^½−1]]^½ζ
 

where ζ = [He^−r/r₀/(m₀c²r₀)], which is the same as γ(hc/(m₀c²r₀).

In principle the two expressions for β can be equated the result solved for r as a function of l. This provides the quantization of r and subsequently that of β and the rest of the characteristics of the system.

The equating of the two expressions for β gives

[½[(1 + 4/ζ²)^½−1]]^½[He^−r/r₀/(m₀c²r)] = He^−r/r₀/(hcl)
which reduces to
[½[(1 + 4/ζ²)^½−1]]^½ = m₀cr/(hl)
 

After squaring and rearranging the above equation reduces to

(1 + 4/ζ²)^½ = 1 + 2m₀²c²r²/(h²l²)
or, after again squaring,
(1 + 4/ζ²) = 1 + 4m₀²c²r²/(h²l²) + 4[m₀²c²r²/(h²l²)]²
and after elimination of the 1's
and division by 4
1/ζ² = m₀²c²r²/(h²l²) + [m₀²c²r²/(h²l²)]²
 

Since 1/ζ=m₀c²r/(He^−r/r₀) the previous equation is equivalent to

[m₀c²r/(He^−r/r₀)]² =
m₀²c²r²/(h²l²) + [m₀²c²r²/(h²l²)]²
which upon division by (m₀c²r)² gives
[e^r/r₀/H]² = 1/(c²h²l²) + m₀²r²/(hl)⁴)
or, equivalently,
e^2r/r₀/H² = 1/(chl)² + (m₀c²)²r²/(chl)⁴
 

Previously the expression H/(hc) was defined as γ. Utilizing this term the previous equation can be put into the form

e^2r/r₀ = γ²/l² + γ⁴(m₀c²/H)²r²/l⁴
 

It is mathematically convenient to deal with z=2r/r₀ as a variable so the above equation takes the form

e^z = γ²/l² + γ⁴(m₀c²r₀/(2H))²z²/l⁴
 

To simplify matters still more for computation let (γ/l)² be denoted as ε and (m₀c²r₀/(2H)) as μ. The equation to be solved for z is then

e^z = ε + μ²ε²z²
 

Although the above equation is transcendental and cannot be solved analytically, approximations to any degree of accuracy can be readily be obtained numerically, provided of course the equation has a solution at all. Since ε and μ² are positive there will always be a negative solution, but a negative solution might be physically meaningless.

The Implication of the Analysis
for the Number of States

The equation e^z=ε+μ²ε²z² has solutions for some values of ε and μ and not for others. At the dividing point between the range for solutions and the range for no solution there would be a tangency solution. In other words there would be a value of z such that not only are the LHS and RHS equal but also the derivative of the two sides are equal. This means that

e^z = ε+μ²ε²z²
and
e^z = 2μ²ε²z
 

Equating the expressions for e^z gives a quadratic equation in z which has a solution that can be put into the form

z = 1 ± (1 − 1/(μ²ε))^½
 

Thus there will be a real valued solution for z if μ²ε≥1. The values of ε and μ reduces this condition to

μ²ε = (mc²r₀/(2hcl))² ≥ 1
which further reduces to
l ≤ mc²r₀/(2hc)
or when numerical values
are substituted into this expression
l ≤ 1.778
 

Thus l can only have the value 1 and hence there is only one angular momentum and energy state for the deuteron.

Empirical Estimates of the Model Parameters

There are several relationships involving H and the radius of the nucleon orbits r₁; i.e.,

(r₁/r₀)e^−r₀/r₀ = l²h²/(Hmr₀)
V=(H/r₀)Ei(ρ)
r₁/r₀ = z/2
where z is the solution to
z = ln(ε) + ln(1+εμ²z²)
and
ε = (H/(hcl)²
μ = m₀c²r₀/(2H)
 

The value of H which is consistent with all relationships must be found by iteration. An iteration scheme which can achieve this is as follows:

• Step 0: Set l=1.
• Step 1: Select a trial value of r₁ and compute ρ=r₁/r₀.
• Step 2: Compute H from ρe^−ρ = l²h²/(Hmr₀)
  
• Step 3: Compute Ei(ρ)=∫_ρ^∞(e^−s/s²)ds.
• Step 4: Compute γ = H/(hc).
• Step 5: Compute ε = (γ/l)².
• Step 6: Compute μ = m₀c²r₀/(2H).
• Step 7: Compute ln(ε) and εμ²
• Step 8: Solve numerically the transcendental equation
  z = ln(ε) + ln(1+εμ²z²)
• Step 9: Obtain ρ=r₁/r₀ as z/2.
• Step 10: Compute r₁=ρ*r₀. Compare this computed value of r₁ with the trial value.
• Step 11: Compute V from the equation V=(H/r₀)Ei(ρ). Compare the computed value of V with the measured value V=2.22457 Mev = 3.565×10⁻¹³ joules. If either the computed value of r₁ is not sufficiently close to the trial value of r₁ or the computed value of V is not sufficiently close to the empirical value repeat Steps 1 through 11.

For illustration, take the trial value of r₁ to be 1.54×10⁻¹⁵ so ρ=2.0. Then H is computed to be 3.24×10⁻²⁶

The values of the parameters of the transcental equation are:

γ=3.24×10⁻²⁶/(3.17×10⁻²⁶)=1.025.
Thus ε=1.051.
The value of μ is then
(1.5×10⁻¹⁰)(7.7×10⁻¹⁶)/(2*3.24×10⁻²⁶))=(1.7823)

The square of μ is then 3.01766 and εμ²=3.3386. The equation to be solved is then

z = 0.05 + ln(1+3.3386z²)
 

Its solution is 4.09025 which, since z=2r/r₀, means that r₁=(4.09025/2)×0.77=1.5756 fermis. This is not far different from the value of r=1.54 fermi upon which the value of H was based. Their potential energy is

V(r₁) = (H/r₀)(∫₂^∞(e^−s/s²)ds)
V(r₁) = [(3.24×10⁻²⁶/(7.7×10⁻¹⁶)](0.0188)
V(r₁) = 7.91065×10⁻¹³ joules = 4.9374 Mev
 

Although this is off by a factor of 2.2 from the empirical value of 2.22457 Mev found for the photodisintergration of deuterons, it is clear the trial value of H is in the ballpark.

If r₁/r₀=2.5 then H has to be equal to 4.326×10⁻²⁶. The computed value of V(r₁) is then

V(r₁) = (H/r₀)(∫_2.5^∞(e^−s/s²)ds)
V(r₁) = [4.326×10⁻²⁶/(7.7×10⁻¹⁶)](0.00792)
V(r₁) = 4.45×10⁻¹³ joules = 2.778 Mev
 

If r₁/r₀=2.7` then H has to be equal to 4.7661×10⁻²⁶. The computed value of V(r₁) is then

V(r₁) = (H/r₀)(∫_2.7^∞(e^−s/s²)ds)
V(r₁) = [4.7661×10⁻²⁶/(7.7×10⁻¹⁶)](0.005709)
V(r₁) = 3.5338×10⁻¹³ joules = 2.206 Mev
 

If r₁/r₀=2.69 then H has to be equal to 4.7320×10⁻²⁶. The computed value of V(r₁) is then

V(r₁) = (H/r₀)(∫_2.69^∞(e^−s/s²)ds)
V(r₁) = [4.732×10⁻²⁶/(7.7×10⁻¹⁶)](0.005802)
V(r₁) = 3.5656×10⁻¹³ joules = 2.2257 Mev
 

If r₁/r₀=2.6905 then H has to be equal to 2.7333×10⁻²⁶. The computed value of V(r₁) is then

V(r₁) = (H/r₀)(∫_2.691^∞(e^−s/s²)ds)
V(r₁) = [4.7333×10⁻²⁶/(7.7×10⁻¹⁶)](0.00579735)
V(r₁) = 3.5637×10⁻¹² joules = 2.22454 Mev
 

For H=4.7333×10⁻²⁶, γ=1.4961, ε=2.2383 and μ=1.2225. For these values of the parameters the solution of the transcendental equation is z=5.39404 and thus ρ=2.697. This is pretty close to the trial value of ρ=2.6904, differing from it by only about 0.25 of 1%. The value of r₁ from this computed value of ρ is thus (2.697)(0.77)=2.0767 fermis as compared with the value of (2.6904)(0.77)=2.072 fermis from the trial value. The tangential velocity of the nucleons may be obtained from the expression for β; i.e., β = γe^−r/r₀/l. Since γ/l=1.4961 and ρ=2.69, β=(1.4961)(0.0679)=0.102. The results indicate that for the deuteron that the velocity of the nucleons in the lowest energy state is 10 percent of the speed of light. Their individual kinetic energies as a ratio of their rest mass energy are given by [(1−β²)^½−1]. This means that in the deuteron the kinetic energy of the nucleons about 0.5 of 1 percent of their rest mass energy.

Since H=4.7333×10⁻²⁶ then H* is four times that value, 1.8933×10⁻²⁵. The constant for the electrostatic force is 2.31×10⁻²⁷. This means that the nuclear force is greater than the electrostatic force up to the point where e^ρ=1.8933/0.23102 or ρ= ln(8.1954)=2.10. This corresponds to a distance 2.1d₀ or 3.24 fermi. The ratio H*/(hc), which corresponds to the fine structure constant, is equal to 5.98. Since this value is greater than unity the analysis of the nuclear force cannot be achieved by a series approximation in the way that the electrostatic force can be where the fine structure constant is 1/137.036. This does not mean that the nuclear force cannot be analyzed; it just means that it has to be done by a different scheme than for the electrostatic force.

Conclusions

The results imply that the formula for nuclear force expressed with respect to particle separation d is:

F = −H*e^−d/d₀/d²
where
H* = 1.8933×10⁻²⁵ kg·m³/s²
and
d₀ = 1.54×10⁻¹⁵ m.

The radius of the single feasible orbit for the nucleons is (2.6904)(0.77)=2.072 fermis. The separation distance d of the nucleons is then 4.143 fermi. The potential energy for a nucleon in the deuteron is then 2.225 Mev. The tangential velocity of the component particles of the deuteron relative to the speed of light is 0.28. The angular velocity is then 2.04×10²³ radians per second. This corresponds to a frequency of 3.245×10²² per second.

The value of ε is 2.2383 and μ equals 1.2225. For these values the solution to the transcendental equation is 5.39404 which is 2r₁/r₀ so r₁/r₀=2.697, very close to the value derived from the photodisintegration energy of the deuteron. The two values are within 0.25 of 1% of each other.

Thus the model of the nuclear force presented here is logically consistent and empirically compatible with the data for the deuteron.

The above formula is for the force between two single nucleons. The general formula for the force between one collection of n₁ nucleons and another of n₂ nucleons is;

F = −H*(e^−d/d₀)n₁n₂/d²
where
H* = 1.8933×10⁻²⁵ kg·m³/s²
and
d₀ = 1.54×10⁻¹⁵ m.

(To be continued.)