There are shell structures of protons and neutrons in nuclei. These shells are manifested in terms of the stability of the nuclides. There are nuclides which are more stable such that when a shell is filled. The shell structure is also manifested in terms of the incremental binding energies as additional nucleons are added. For example consider the incremental binding energies for the isotopes of tin (Sn):

The breakpoints come at certain numbers, called magic numbers, that represent filled shells. The sawtooth patterns comes from the formation of neutron pairs. Protons also form pairs. The declining slope of the incremental binding energy pattern indicates that the binding energy is a quadratic function of the number of nucleons.

The slopes and curvature of the relationships differ for different shells. Instead of trying to allow for such differences in a statistical analysis of all 2931 nuclides this analysis looks at the cases in which the number of protons is between 29 and 50 and the number of neutrons is also between 29 and 50. The neutron and proton shells have capacity of 22 nucleons each.

There 285 nuclides satisfying those conditions. These nuclides have an average binding energy of 666.36 MeV.

The explanatory variables for the binding energies are the numbers of protons p and the number of neutrons n in the shell. In order to capture the effect of the pairing of nucleons the numbers of the nucleons are expressed as the number of pairs of each of the nucleons and whether there are singleton (unpaired) nucleons. The number of pairs of protons and neutrons is denoted as #pp and #nn, respectively. From these variables three additonal explanatory variables are created (#pp)², (#pp)*(#nn) and (#nn)². Additionally there are the variables, s_p and s_n, which are equal to 1 if a singleton proton or singleton neutron is present in a nuclide.

The results of the regression are

BE = 823.41387 + 2.04835#pp + 23.80268#nn
- 0.932325(#pp)² + 1.15657(#pp)*(#nn) - 0.38569(#nn)²
+ 2.61761s_p + 8.733316s_n

[545.1] [11.9] [193.3]
[-66.2] [69.9] [-53.8]
[15.9] [53.8]

R² = 0.999695

Standard error = 1.51064 Mev

Coefficient of Variation = (1.51064 MeV)/(1074.06MeV) = 0.0014065

The numbers in the square brackets below the coefficients are the t-ratios for the coefficients. For the regression coefficient to be statistically signficant at the 95 percent level of confidence its t-ratio must be roughly 2.0 or larger.

The binding energy of the nuclide with 50 protons and 50 neutrons (Silicon) is 847.2MeV. This is significantly different from the regression constant of 823.41387 MeV.

There was a possibility that if a singleton proton and a singleton neutron were present in a nuclide they would form a pair which would enhance the binding energy. When such a variable was included in the regression its coefficient was not statisically significant. (The t-ratio was 0.5.)

Conclusions

The numbers of protons and neutrons in the nuclides having 51 to 82 of each explain all but 0.03 of 1 percent in the variation in binding energies of these nuclides. The regression equation gives estimates of binding energies that are accurate to roughly ± 0.14 of 1 percent.