4.4 Nuclear mass formula

There is more structure in Fig. 4.1 than just a simple linear dependence on A. A naive analysis suggests that the following terms should play a rôle:

  1. Bulk energy: This is the term studied above, and saturation implies that the energy is proportional to Bbulk = αA.
  2. Surface energy: Nucleons at the surface of the nuclear sphere have less neighbours, and should feel less attraction. Since the surface area goes with R2, we find Bsurface = βA.
  3. Pauli or symmetry energy: nucleons are fermions (will be discussed later). That means that they cannot occupy the same states, thus reducing the binding. This is found to be proportional to Bsymm = γ(N2 Z2)2A2.
  4. Coulomb energy: protons are charges and they repel. The average distance between is related to the radius of the nucleus, the number of interaction is roughly Z2 (or Z(Z 1)). We have to include the term BCoul = 𝜖Z2A.

Taking all this together we fit the formula

B(A,Z) = αA βA23 γ(A2 Z)2A1 𝜖Z2A13 (4.8)

to all know nuclear binding energies with A 16 (the formula is not so good for light nuclei). The fit results are given in table 4.1.


Table 4.1: Fit of masses to Eq. (4.8)
.



parameter value


α15.36 MeV
β16.32 MeV
γ90.45 MeV
𝜖0.6928 MeV




mass˙form1


Figure 4.2: Difference between fitted binding energies and experimental values, as a function of N and Z.

In Fig. 4.3 we show how well this fit works. There remains a certain amount of structure, see below, as well as a strong difference between neighbouring nuclei. This is due to the superfluid nature of nuclear material: nucleons of opposite momenta tend to anti-align their spins, thus gaining energy. The solution is to add a pairing term to the binding energy,

Bpair = A12 for N odd, Z odd A12for N even, Z even (4.9)

The results including this term are significantly better, even though all other parameters remain at the same position, see Table 4.2. Taking all this together we fit the formula

B(A,Z) = αA βA23 γ(A2 Z)2A1 δB pair(A,Z) 𝜖Z2A13 (4.10)


Table 4.2: Fit of masses to Eq. (4.10)



parameter value


α15.36 MeV
β16.32 MeV
γ90.46 MeV
δ11.32 MeV
𝜖0.6929 MeV




mass˙form2


Figure 4.3: Difference between fitted binding energies and experimental values, as a function of N and Z.


mass2


Figure 4.4: BA versus A, mass formula subtracted