Set up the Binding Energy per Nucleon Equations

u=931.494;
a=15.5;
b=16.8;
s=23.20;
d=0.72;
δ[N_,Z_]:=(Piecewise[{{34., OddQ[N] && OddQ[Z]}, {0.0, Xor[OddQ[N]&&EvenQ[Z],EvenQ[N]&&OddQ[Z]]  }, {-34., EvenQ[N] && EvenQ[Z]}}])
B[A_,Z_] := (a A - b (A^(2/3)) - s (((((A - 2 Z))^2) )/A)- ((d Z (Z-1))/(A^(1/3)))- ((δ[A-Z,Z])/(A^(3/4))))/A

Table[B[57,Z],{Z,21,29}]

Set up the Mass Parabola Equation

Mass[A_,Z_]:=((A-Z) 939.566 + Z 938.272 +Z 0.51099892- B[A,Z] A)/u

B[119,50]

Plot the Mass Parabolas

zStart=19;
array1=Table[{Z,Mass[56,Z]},{Z,zStart,31,2}];
array2=Table[{Z,Mass[56,Z]},{Z,zStart-1,31,2}];

ListLinePlot[{array1,array2}]

ListPlot[{array1,array2},PlotMarkers->{"☺","😐","☹"}]

array3a=Table[{Z,Mass[7,Z]},{Z,1,7}];
array3b=Table[{Z,Mass[23,Z]},{Z,9,14}];array3c=Table[{Z,Mass[207,Z]},{Z,79,85}];
array3d=Table[{Z,Mass[235,Z]},{Z,89,95}];
p1=ListLinePlot[array3a];
p2=ListLinePlot[array3b];
p3=ListLinePlot[array3c];
p4=ListLinePlot[array3d];

GraphicsGrid[{{p1,p2},{p3,p4}}]

Mass[3,1]
Mass[3,2]

QuantityMagnitude[IsotopeData["Sodium22","Lifetime"]]/(3.15 10^7)Log[2]

mH3=QuantityMagnitude[IsotopeData["Hydrogen3","AtomicMass"]] u
mHe4=QuantityMagnitude[IsotopeData["Helium4","AtomicMass"]] u
mCu64=QuantityMagnitude[IsotopeData["Copper64","AtomicMass"]] u
mSn119=QuantityMagnitude[IsotopeData["Tin119","AtomicMass"]] u

Binding Energy per nucleon using Semiempirical Formula

B[3,1]
B[4,2]
B[64,29]
B[119,50]

Binding energy per nucleon using known nuclear masses

me=QuantityMagnitude[ParticleData["Electron","Mass"]]
mp=QuantityMagnitude[ParticleData["Proton","Mass"]]
mn=QuantityMagnitude[ParticleData["Neutron","Mass"]]

bOvera[A_,Z_,m_] = ((mn (A-Z) + (mp+me) Z -m )/A);

bOvera[3,1,mH3]
bOvera[4,2,mHe4]
bOvera[64,29,mCu64]
bOvera[119,50,mSn119]

be[A_,Z_] := (a A - b (A^(2/3)) - s (((((A - 2 Z))^2) )/A)- ((d Z (Z-1))/(A^(1/3)))- ((δ[A-Z,Z])/(A^(3/4))))

be[13,6]
be[13,7]

mC13=QuantityMagnitude[IsotopeData["Carbon13","AtomicMass"]] u
mN13 = QuantityMagnitude[IsotopeData["Nitrogen13","AtomicMass"]] u

beAtomicMasses[A_,Z_,m_] = mn (A-Z) + (mp+me) Z -m ;

beAtomicMasses[13,6,mC13]
beAtomicMasses[13,7,mN13]

beCoulomb[A_,Z_] = ((d Z (Z-1))/(A^(1/3)));

beCoulomb[13,6]
beCoulomb[13,7]

mPa234 = QuantityMagnitude[IsotopeData["Protactinium234","AtomicMass"]] u
mPa235 = QuantityMagnitude[IsotopeData["Protactinium235","AtomicMass"]] u
mUr234 = QuantityMagnitude[IsotopeData["Uranium234","AtomicMass"]] u
mUr235 = QuantityMagnitude[IsotopeData["Uranium235","AtomicMass"]] u
mUr236 = QuantityMagnitude[IsotopeData["Uranium236","AtomicMass"]] u
mNp234 = QuantityMagnitude[IsotopeData["Neptunium234","AtomicMass"]] u
mNp235 = QuantityMagnitude[IsotopeData["Neptunium235","AtomicMass"]] u
mNp236 = QuantityMagnitude[IsotopeData["Neptunium236","AtomicMass"]] u

seN[mi_,mf_,mn_]=mf + mn -mi ;
seP[mi_,mf_,mp_] = mf + mp +me-mi ;

seP[mUr236,mPa235,mp]
seP[mUr235,mPa234,mp]
seP[mNp236,mUr235,mp]
seP[mNp235,mUr234,mp]

seN[mUr236,mUr235,mn]
seN[mUr235,mUr234,mn]
seN[mNp236,mNp235,mn]
seN[mNp235,mNp234,mn]

143/92//N