The semi-empirical mass formula
In nuclear physics, the semi-empirical mass formula is a formula for
calculating the approximate nuclear binding energy \(B\) of an atomic nucleus
with atomic number \(Z\) and mass number \(A\) :
$$
B=a_{1} A-a_{2} A^{2 / 3}-a_{3} \frac{Z^{2}}{A^{1 / 3}}-a_{4} \frac{(A-2
Z)^{2}}{A}+\frac{a_{5}}{A^{1 / 2}}
$$
where, in units of millions of electron volts, the constants are \(a_{1}=15.8,
a_{2}=18.3\), \(a_{3}=0.714, a_{4}=23.2\), and
$$
a_{5}= \begin{cases}0 & \text { if } A \text { is odd, } \\ 12.0 & \text { if
} A \text { and } Z \text { are both even, } \\ -12.0 & \text { if } A \text
{ is even and } Z \text { is odd. }\end{cases}
$$
a) Write a program that takes as its input the values of \(A\) and \(Z\), and
prints out the binding energy for the corresponding atom. Use your program to
find the binding energy of an atom with \(A=58\) and \(Z=28\). (Hint: The correct
answer is around \(500 \mathrm{MeV}\) )
b) Modify your program to print out not the total binding energy B, but the
binding energy per nucleon, which is \(B / A\).
c) Now modify your program so that it takes as input just a single value of
the atomic number \(Z\) and then goes through all values of \(A\) from \(A=Z\) to
\(A=3 Z\), to find the one that has the largest binding energy per nucleon. This
is the most stable nucleus with the given atomic number. Have your program
print out the value of \(A\) for this most stable nucleus and the value of the
binding energy per nucleon
d) Modify your program again so that, instead of taking \(Z\) as input, it runs
through all values of \(Z\) from 1 to 100 and prints out the most stable value
of \(A\) for each one. At what value of \(Z\) does the maximum binding energy per
nucleon occur? (The true answer, in real life, is \(Z=28\), which is nickel.)