^{1}

^{*}

^{1}

A formula is derived for the central nucleon-nucleon potential, based on an analysis of the physical origin of the nucleon-nucleon attraction by pion exchange. The decrease of the dynamical mass of the interaction field, exchanged pion in this case, is the principal mechanism responsible for the nuclear attraction in a similar way that the decrease of the kinetic energy of the exchange electron in the diatomic molecule is directly responsible for the covalent molecular attraction. The minimum value of this central nucleon-nucleon potential and the position of the minimum are similar with the values reported in literature for a potential calculated by lattice QCD, which shares the features of the phenomenological nucleon-nucleon potentials. The Schrodinger equation with this central nucleon-nucleon potential was solved numerically for different values of the pion mass. The binding energy increases with the decrease of the pion mass. For masses higher than the real pion mass the nucleon-nucleon system is unbound. We discuss on the two pion exchange and hard core repulsion. The minimum value of the potential for two pion exchange is comparable with the minimum value of the CD Bonn potential. For a hard core radius of 0.5 fm the binding energy is equal to the deuteron binding energy.

The effective degrees of freedom in the nuclear interaction at low energy, in particular in the nuclear bound state, are the nucleons and the pions. The pion exchange is the basic mechanism of nucleon-nucleon (NN) attraction at low energy [

As it is well known the current masses of quarks (antiquarks) u (ū) and d (đ) are very small; the nucleon and pion masses are mainly of dynamical origin (“kinetic” energy). The confinement of the quarks into nucleon and pion associate an energy, given by the Heisenberg uncertainty relation, which is just the dynamical mass.

The long range structure of the nucleon is given roughly by the pion [

At the formation of a nuclear bound state this dynamical mass decreases, which determines a mass defect and consequently a binding energy for the nucleon-nucleon state [

E Δ ≅ ℏ c ƛ π + Δ ( R ) (1)

and is lower than the initial one (that in the free nucleon). To form a bound state, the decrease of the dynamical mass of the pion degree of freedom:

Δ E = ℏ c ƛ π − ℏ c ƛ π + Δ ( R ) (2)

must be larger than the kinetic energy acquired by the system of two nucleons due to their localization at a distance R each other:

This mechanism of nuclear binding has similarities with the mechanism of molecular binding of diatomic molecules [

Let’s start with the physical interpretation of the mechanism of H 2 + ion binding presented by Feynman in [

In the H 2 + ion, since there are two protons, there is more space where the electron can have a low potential energy than in the case of hydrogen atom. The exchanged electron spreads out lowering its kinetic energy, in accord with uncertainty relation. This kinetic energy decrease is at the origin of the molecular attraction in covalent bond, in particular in the H 2 + ion [

For large distances between the two protons of the H 2 + ion the electrostatic potential energy of the exchanged electron is nearly zero over most of the space between the protons and the electron moves nearly like a free particle in empty space but with a negative energy [

p 2 2 m e = − W H (3)

where W H is the binding energy (+13.6 eV) of the hydrogen atom. This means that p is an imaginary number:

p = i 2 m e W H (4)

The probability amplitude A for a particle of definite energy to get from one place to another a distance R away is proportional to [

A ~ e ( i / ℏ ) p R R (5)

where p is the momentum corresponding to the definite energy. Replacing p one obtains that the amplitude of jumping of electron from one proton to the other, for large separation R between the two protons, will vary as [

A ~ e − ( 2 m e W H / ℏ ) R R = e − R / a 0 R (6)

where a 0 is the Bohr radius.

If the particle goes in one direction the amplitude is [

A ~ e − R / a 0 (7)

One can note that this exponential function limits drastically the amplitude of electron exchange for large separation.

The nuclear interaction which takes place between a neutron and a proton by pion exchange is described by Feynman with similar arguments [

E 2 = p 2 c 2 + m π 2 c 4 (8)

where E and p are the total energy and the momentum of the pion.

Since the exchanged pion have practically zero total energy the momentum is again imaginary [

p = i m π c (9)

This means the amplitude for the pion to jump from one nucleon to another is for large R:

A ~ e − ( m π c / ℏ ) R = e − R / ƛ π (10)

The exponential function is typical for a Yukawa potential or exponential potential, and again it limits drastically the exchange for large R.

In fact, the exponential factor is well known from the tunneling of a potential barrier of width R by a particle with an energy much lower than the barrier height. The probability of transmission is the (absolute) square of the amplitude [

P ~ e − 2 R / ƛ π (11)

The increase Δ ( R ) of the radius of localization region of the pion degree of freedom, which appears in formulae (1) and (2), is strongly limited by this exponential function (11), i.e. the probability for the exchanged pion to penetrate (tunnel) the potential barrier between the two nucleons. From a physical point of view one expects that this probability of transmission is 1 for a barrier width R → 0 . Therefore Δ ( R ) , which is proportional both to the distance R between the nucleons and to the probability of transmission of the exchanged pion, is equal to:

Δ ( R ) = R e − 2 R / ƛ π (12)

By replacing (12) in formula (2), we obtain the decrease of the dynamical mass of the pion degree of freedom, which is at the origin of nuclear attraction. In fact with sign minus this decrease is just the NN potential due to pion exchange [

V ( R ) = − Δ E = ℏ c ƛ π + Re − 2 R / ƛ π − ℏ c ƛ π = − ℏ c ƛ π * R e − 2 R / ƛ π ƛ π + R e − 2 R / ƛ π (13)

where ℏ c / ƛ π = m π c 2 .

The NN potential V(R) as a function of inter-nucleon distance R is shown in

The potential in

For a lower pion mass (360 MeV) the range of the potential gets wider also [

We solved numerically the Schrodinger equation for the potential V(R) given in relation (13), for different values of the pion mass m π :

− ℏ 2 2 μ Δ ψ ( R ) + V ( R ) ψ ( R ) = E ψ ( R ) (14)

where μ is the reduced mass of the two nucleons.

From a numerical point of view, the Schrodinger equation in central potential has been solved using the substitution: ψ ( R ) = ϕ ( R ) / R which, in turn, gives a standard Sturm-Liouville eigenvalue problem with a constant coefficient for the second order derivative:

− ℏ 2 2 μ d 2 ϕ ( R ) d R 2 + V ( R ) ϕ ( R ) = E ϕ ( R ) (15)

We look for the ground state, this means zero centrifugal energy (l = 0).

The spatial dimension has been truncated as R ∈ [ 0 , 20 ] fm with a standard discretization in equal intervals of Δ R = 10 − 3 fm . Such large radial extension is needed in order to resolve properly the states lying closely to the continuum ( E < ≅ 0 ). The eigenvalue problem is solved by means of a finite difference method with the boundary values ϕ ( 0 ) = ϕ ′ ( 0 ) = 0 and imposing an exponential decay at large radius. The resulting Hamiltonian is diagonalized and the eigenvalues (energy) are obtained.

In

The total energy of the two nucleons decreases, this means the binding energy increases with the decrease of the pion mass. For masses higher than the real pion mass the total energy gets positive (unbound state). This underline the central role of the pion as main player in the production of nuclear attraction, as largely accepted in literature [

The maximum molecular attraction is realized in the diatomic molecular bond by exchange of two electrons [_{2} in which the two hydrogen atoms put in common (exchange) their electrons.

Similarly, the maximum of NN attraction is given by the exchange of two pions. Each nucleon puts in common (exchanges) a pion degree of freedom with the other nucleon. This means two pion degrees of freedom are slightly de-localized. The NN potential in this case is practically two times the potential given by relation (13). The minimum value of the potential gets −44 MeV, which is comparable with the minimum value (−50 MeV) of the CD Bonn potential [

If one adds to this two pion exchange potential 2 V(R), where V(R) is given by relation (13), a hard core repulsion at r 0 ≤ 0.6 fm , a typical value for phenomenological nucleon-nucleon potentials [

The Schrodinger equation for this two pion exchange potential 2V(R) with hard core repulsion was solved numerically for different values of the hard core radius r 0 . In

For a value of the hard core repulsion radius equal to 0.5 fm, the binding energy is equal to the deuteron binding energy. For this hard core radius

( r 0 = 0.5 fm ) in

With pion mass increase, for pion masses higher than about 60 MeV, the total energy increases, i.e. the binding energy decreases, and gets zero at about 190 MeV, a result comparable with that obtained in [

The pion exchange is at the origin of the nuclear attraction, in a similar way that the electron exchange is at the origin of attraction in the molecular covalent bond. The decrease of the kinetic energy (in fact dynamical mass decrease from a relativistic point of view) of the exchange electron in the H 2 + ion, directly responsible for the formation of the molecular bound state, is replaced by the decrease of the dynamical mass of the pion degree of freedom in the case of the nuclear attraction by pion exchange. The slight de-localization of the pion degree of freedom, which is at the origin of this dynamical mass decrease, is drastically limited by an exponential function, which represents the probability for the pion to penetrate the potential barrier between the two nucleons. A similar exponential function exists in the case of molecular bond.

The analytical formula of the central nuclear potential (13) derived for the NN interaction by pion exchange does not contain any unknown parameter. The minimum value of the NN potential and the position of the minimum are similar with the values reported in literature for the central NN potential obtained by lattice QCD, which shares the features of the phenomenological NN potentials. A very small binding energy (0.1 MeV) was obtained by solving numerically the Schrodinger equation. The binding energy increases for pion masses lower than the real pion mass. For masses higher than the real pion mass the nucleon-nucleon system is unbound.

The fall of the potential (13) towards zero value for small R (

Since the maximum value of R e − 2 R / ƛ π is 0.26 fm (for R = ƛ π / 2 ), this means substantially smaller than ƛ π = 1.41 fm , relation (13) can be written in a good approximation as:

V ( R ) ≈ − m π c 2 R ƛ π e − 2 R / ƛ π (16)

The NN nuclear potential is proportional to the mass of the interaction field. It is also proportional to the ratio R / ƛ π , which is directly related to the slight de-localization of the pion degree of freedom and its (dynamical) mass decrease. This de-localization is drastically limited by the exponential function, which is similar to an exponential potential except the factor 2. Due to this exponential function the width of potential (16) gets larger for pion masses lower than the real pion mass and this explains the increase of the binding energy with the pion mass decrease (

In the case of two pion exchange the minimum of the potential gets comparable with the minimum value of the CD Bonn potential. The potential in this case is in a good approximation:

V 2 π ( R ) ≈ − 2 m π c 2 R ƛ π e − 2 R / ƛ π (17)

The dependence of the type x e − x of the nuclear potential, where x = 2 R / ƛ π , is similar with the attractive part of the Rydberg potential used to describe the molecular covalent bonding [

A hard core repulsion was added to this two pion exchange potential and the Schrodinger equation was solved numerically for different values of the hard core repulsion radius. For a radius of 0.5 fm the binding energy is equal to the deuteron binding energy.

Let’s compare the strength of the nuclear potential V 2 π ( R ) from relation (17)

with the coulombian potential V C = q 2 R . The ratio of the two potentials for R = ƛ π is:

V 2 π ( ƛ π ) / V C ( ƛ π ) = 2 ℏ c q 2 e − 2 = 2 1 α e − 2 (18)

where α is the e.m coupling constant, a typical result for the relative strength of the nuclear interaction to the e.m. interaction.

If we analyze the mechanism of NN interaction at quark level, we could say that by pion exchange between two nucleons some quark degrees of freedom are implicitly exchanged and in consequence are slightly de-localized. The confinement region of a quark slightly increases and accordingly its dynamical mass decreases. This suggests to interpret the nuclear interaction as a residual strong interaction.

We thank E. Dudas and V. V. Baran for helpful discussions and suggestions. Nicolae Mandache thanks to Sola Florinuta for constant support.

Mandache, N.B. and Palade, D.I. (2018) On the Mechanism of Nucleon-Nucleon Attraction by Pion Exchange. Journal of Modern Physics, 9, 1459-1468. https://doi.org/10.4236/jmp.2018.98090