Chapter 4 Exercises — The Semi-Empirical Mass Formula

Use the following SEMF parameters unless otherwise specified: $a_V = 15.75$ MeV, $a_S = 17.80$ MeV, $a_C = 0.711$ MeV, $a_{\text{sym}} = 23.7$ MeV, $a_P = 11.2$ MeV.

Reference experimental binding energies (AME2020) for selected nuclei:

Nucleus $Z$ $A$ $B_{\text{exp}}$ (MeV)
$^4$He 2 4 28.296
$^{12}$C 6 12 92.162
$^{16}$O 8 16 127.619
$^{40}$Ca 20 40 342.052
$^{48}$Ca 20 48 415.991
$^{56}$Fe 26 56 492.254
$^{62}$Ni 28 62 545.259
$^{90}$Zr 40 90 783.893
$^{120}$Sn 50 120 1020.543
$^{197}$Au 79 197 1559.402
$^{208}$Pb 82 208 1636.430
$^{235}$U 92 235 1783.871
$^{238}$U 92 238 1801.695

Section A: Binding Energy Calculations (Fundamental)

Problem 4.1 — Calculate each term of the SEMF separately for $^{40}$Ca ($Z = 20$, $A = 40$) and report: (a) volume, (b) surface, (c) Coulomb, (d) asymmetry, (e) pairing, and (f) total binding energy. (g) Compare to the experimental value $B_{\text{exp}} = 342.052$ MeV. (h) What is the percentage error? (i) $^{40}$Ca is doubly magic ($Z = N = 20$). Does the SEMF overestimate or underestimate the binding energy? Why?

Problem 4.2 — Repeat Problem 4.1 for $^{120}$Sn ($Z = 50$, $A = 120$). Note that $^{120}$Sn has a magic proton number ($Z = 50$). What is the sign and magnitude of the residual $B_{\text{exp}} - B_{\text{SEMF}}$?

Problem 4.3 — Calculate $B/A$ using the SEMF for $A = 12, 56, 120, 208, 238$. Plot (or tabulate) these values alongside the experimental values. At what mass number does the SEMF give the best fractional accuracy?

Problem 4.4 — For $^{56}$Fe, calculate the contribution of each SEMF term to $B/A$ (in MeV/nucleon). Which term is the largest correction to the volume term? Express each correction as a percentage of $a_V$.

Problem 4.5 — Calculate the atomic mass $M(Z,A)$ in atomic mass units (u) for $^{238}$U using your SEMF binding energy. Compare to the tabulated mass $M = 238.05079$ u. (Use $m_p = 938.272$ MeV$/c^2$, $m_n = 939.565$ MeV$/c^2$, $m_e = 0.511$ MeV$/c^2$, $1$ u $= 931.494$ MeV$/c^2$.)


Section B: The Individual Terms (Conceptual and Quantitative)

Problem 4.6 — (a) If the nuclear force did NOT saturate — that is, if every nucleon interacted equally with every other nucleon — the volume term would be replaced by a term proportional to $A(A-1)/2$. Show that the binding energy per nucleon in this case would increase linearly with $A$, rather than being approximately constant. (b) Explain why this would make heavy nuclei overwhelmingly favored over light nuclei, eliminating the peak in the $B/A$ curve. (c) Would fission still be energetically favorable in this scenario?

Problem 4.7 — The surface-to-volume ratio of a sphere of radius $R$ is $3/R$. (a) Show that the surface correction to $B/A$ scales as $A^{-1/3}$. (b) Calculate this correction for $A = 10, 50, 100, 200$. (c) At what mass number does the surface correction become less than 10% of the volume term?

Problem 4.8 — (a) Starting from the expression for the electrostatic self-energy of a uniform sphere, $U_C = (3/5)(Ze)^2/(4\pi\epsilon_0 R)$, calculate $a_C$ using $r_0 = 1.20$ fm and $e^2/(4\pi\epsilon_0) = 1.4400$ MeV$\cdot$fm. (b) How does this compare to the fitted value $a_C = 0.711$ MeV? (c) What physical effects might account for the small difference? (Hint: consider the finite size of the proton, the diffuse nuclear surface, and exchange corrections from the Pauli principle applied to protons.)

Problem 4.9 — The Coulomb energy of $^{208}$Pb is enormous. (a) Calculate the total Coulomb energy $E_C = a_C Z(Z-1)/A^{1/3}$ for $^{208}$Pb. (b) Express this as a fraction of the total binding energy. (c) Imagine you could "turn off" the Coulomb force. What would the binding energy of $^{208}$Pb be? What would $B/A$ be? (d) Would the most stable isobar at $A = 208$ still have $Z = 82$?

Problem 4.10 — Using the Fermi gas model, the asymmetry energy coefficient can be estimated as $a_{\text{sym}} \approx \epsilon_F / 3$, where $\epsilon_F \approx 38$ MeV is the nucleon Fermi energy. (a) Calculate $a_{\text{sym}}$ from this estimate and compare to the fitted value. (b) The discrepancy suggests that the potential energy also depends on the neutron-proton asymmetry. Explain physically why this should be the case.

Problem 4.11 — For the $A = 101$ isobars (odd $A$), list the known stable nuclide(s). (The answer is $^{101}$Ru, $Z = 44$.) (a) Calculate $Z_{\text{stable}}$ using the SEMF. (b) Calculate the SEMF binding energy for $Z = 42, 43, 44, 45, 46$ at $A = 101$. (c) Which of these has the highest predicted binding energy? (d) Does the SEMF correctly identify the most stable isobar?


Section C: Valley of Stability and Drip Lines

Problem 4.12 — (a) Derive the formula $Z_{\text{stable}} = A / [2 + (a_C / 4a_{\text{sym}}) A^{2/3}]$ from the condition $\partial B / \partial Z = 0$. (b) Evaluate $Z_{\text{stable}}$ for $A = 20, 40, 80, 120, 200, 250$. (c) Compare your results to the known most stable nuclei at each $A$. (d) At what mass number does the neutron-to-proton ratio $N/Z$ first exceed 1.5?

Problem 4.13 — (a) Show that the curvature of the mass parabola at fixed $A$ is $d^2 M / dZ^2 = 2(a_C A^{-1/3} + 4a_{\text{sym}} / A)$. (b) Evaluate this curvature for $A = 30$ and $A = 200$. (c) Which mass parabola is steeper? (d) What does this imply about the range of $Z$ values for which nuclei are bound at each $A$?

Problem 4.14 — For even $A$, the even-even and odd-odd mass parabolas are separated by $2\delta_0$. (a) Calculate $2\delta_0$ for $A = 100$. (b) How many stable isobars does the SEMF predict at $A = 100$? (c) The experimentally known stable isobars at $A = 100$ are $^{100}$Mo ($Z = 42$) and $^{100}$Ru ($Z = 44$). Is $^{100}$Tc ($Z = 43$, odd-odd) stable? Why or why not? (d) Can you explain why two even-even isobars can both be stable at the same $A$?

Problem 4.15Neutron separation energy. (a) Using the SEMF, calculate the neutron separation energy $S_n = B(Z,N) - B(Z,N-1)$ for $^{208}$Pb and $^{209}$Pb. (b) The experimental values are $S_n(^{208}\text{Pb}) = 7.368$ MeV and $S_n(^{209}\text{Pb}) = 3.937$ MeV. Explain the dramatic difference in terms of the $N = 126$ shell closure. (c) Does the SEMF predict this difference? Why or why not?

Problem 4.16Proton drip line estimate. The proton drip line occurs where $S_p = 0$. (a) For the tin isotopes ($Z = 50$), estimate the minimum $N$ for which $S_p > 0$ using the SEMF. (b) The lightest known bound tin isotope is $^{100}$Sn ($N = 50$). Does the SEMF prediction agree? (c) Why is $^{100}$Sn particularly interesting? (Hint: it is doubly magic.)

Problem 4.17Neutron drip line estimate. (a) For the oxygen isotopes ($Z = 8$), estimate the maximum $N$ for which $S_n > 0$ using the SEMF. (b) The experimentally determined heaviest bound oxygen isotope is $^{24}$O ($N = 16$). Does the SEMF agree? (c) $^{24}$O has $N = 16$, which is NOT a traditional magic number. Recent research suggests $N = 16$ becomes a magic number for very neutron-rich light nuclei. What does this tell you about the reliability of the SEMF near the drip lines?

Problem 4.18The end of the periodic table. (a) Using the fissility parameter $x = a_C Z^2 / (2a_S A)$, calculate $x$ for $^{238}$U, $^{294}$Og ($Z = 118$, $A = 294$), and a hypothetical nucleus with $Z = 130$, $A = 330$. (b) At what $Z$ does $x = 1$ (instantaneous fission) assuming $A = 2.5Z$ (a rough extrapolation of the stability line)? (c) What does this predict about the existence of elements beyond $Z \sim 130$?


Section D: Comparative and Analytical Problems

Problem 4.19$^{56}$Fe vs. $^{62}$Ni: Which is more tightly bound? (a) Calculate $B/A$ using the SEMF for both nuclei. (b) The experimental values are $B/A(^{56}\text{Fe}) = 8.790$ MeV and $B/A(^{62}\text{Ni}) = 8.795$ MeV. Which nucleus actually has the highest $B/A$? (c) Why is iron, not nickel, the endpoint of stellar fusion? (Hint: the relevant quantity for stellar physics is not $B/A$ but the mass per nucleon, and the proton-to-nucleon ratio matters because of electron capture in stellar cores. See Case Study 4.1 for a full discussion.)

Problem 4.20Mirror nuclei. Mirror nuclei have the proton and neutron numbers interchanged: $(Z, N) \leftrightarrow (N, Z)$. (a) Show that the SEMF predicts the binding energy difference between mirror nuclei arises solely from the Coulomb term: $\Delta B = B(Z+1, A) - B(Z, A) \approx -a_C (2Z+1)/A^{1/3} + 4a_{\text{sym}}(A-2Z-2)/A$ ... wait — actually show that for mirror nuclei with the same $A$, the asymmetry term is the same and the difference comes from the Coulomb term. (b) Calculate $\Delta B$ for the mirror pair $^{15}$O / $^{15}$N. (c) Compare to the experimental value $\Delta B = B(^{15}\text{N}) - B(^{15}\text{O}) = 3.537$ MeV. (d) The Coulomb displacement energy provides a test of charge independence of the nuclear force. Does the agreement support charge independence?

Problem 4.21Sensitivity analysis. (a) Increase $a_V$ by 5% while holding all other parameters fixed. How does $B/A$ change for $^{56}$Fe and $^{238}$U? (b) Increase $a_C$ by 5%. How does $Z_{\text{stable}}$ change for $A = 200$? (c) Which SEMF parameter has the largest effect on the position of the valley of stability? (d) Which parameter has the largest effect on the predicted mass of superheavy elements?

Problem 4.22The isobaric multiplet mass equation (IMME). For a set of isobars with the same $A$ and different $Z$, the SEMF predicts that the mass varies quadratically with $Z$ (or equivalently, with the isospin projection $T_z = (N-Z)/2$). (a) Show that $M(A, T_z) = a + b T_z + c T_z^2$ where you express $a$, $b$, and $c$ in terms of SEMF parameters. (b) This quadratic dependence on $T_z$ is called the isobaric multiplet mass equation. How well does it work experimentally? (The IMME works to better than 10 keV for most nuclear isobaric multiplets — far better than the SEMF itself.)


Section E: Computational Problems

💻 Problem 4.23SEMF implementation and visualization. Write a Python program that: (a) Implements the SEMF as a function B_SEMF(Z, A, params). (b) Plots $B/A$ vs. $A$ for the most stable isobar at each $A$ from $A = 1$ to $A = 300$. (c) Overlays experimental $B/A$ data from AME2020 (you may use the provided semf_fit.py data table or load from a file). (d) Identifies by inspection where the SEMF deviates most from experiment.

💻 Problem 4.24Fitting the SEMF. Using scipy.optimize.curve_fit, fit the five SEMF parameters to experimental binding energies for nuclei with $A \geq 20$. (a) Report the best-fit values and their uncertainties. (b) Calculate the RMS residual. (c) Plot the residuals (experiment minus SEMF) as a function of $N$ and identify the magic numbers from the residual pattern. (d) Repeat the fit excluding all nuclei within 3 units of a magic number ($N$ or $Z = 2, 8, 20, 28, 50, 82, 126$). How do the fitted parameters change? How does the RMS residual change?

💻 Problem 4.25Valley of stability visualization. Write a Python program that: (a) For each $A$ from 1 to 300, finds the $Z$ that maximizes $B_{\text{SEMF}}(Z,A)$. (b) Plots $Z_{\text{stable}}$ vs. $A$ (the beta-stability line) alongside the $N = Z$ line. (c) Overlays the experimentally known stable nuclei. (d) Shades the region between the estimated proton and neutron drip lines (where $S_p = 0$ and $S_n = 0$ respectively).

💻 Problem 4.26Mass parabolas. For $A = 100$ (even) and $A = 101$ (odd), plot the nuclear mass $M(Z,A)$ as a function of $Z$ using the SEMF. For $A = 100$, show both the even-even and odd-odd parabolas. Mark the predicted and observed stable isobars on each plot. Discuss the phenomenon of double beta decay for even-$A$ isobars.

💻 Problem 4.27Shell effects in residuals. Plot the SEMF residual ($B_{\text{exp}} - B_{\text{SEMF}}$) as a 2D color map on the $(N, Z)$ plane for all nuclei in the AME2020 database. Use a diverging colormap centered at zero. The magic numbers should appear as ridges of extra binding (positive residuals) in this plot. Annotate the magic numbers.


Section F: Challenge Problems

Problem 4.28Extended SEMF. Add a sixth term to the SEMF to approximately account for shell effects. One approach: add a correction $\Delta B_{\text{shell}} = \sum_i c_i \exp[-(N - N_i^{\text{magic}})^2 / \sigma^2] + \sum_j c_j \exp[-(Z - Z_j^{\text{magic}})^2 / \sigma^2]$ with magic numbers at 2, 8, 20, 28, 50, 82, 126. Fit the additional parameters $c_i$ (the shell correction strengths) and $\sigma$ (the width of the shell effects). How much does the RMS residual decrease? Is the improvement uniform across the chart or localized?

Problem 4.29The island of stability. Superheavy elements ($Z > 103$) exist only because of shell effects that stabilize them against the fission predicted by the liquid drop model. (a) Using the fissility parameter, estimate the half-life of an element with $Z = 114$, $A = 298$ based on the liquid drop model alone. (b) Explain qualitatively why the shell model prediction of a closed shell at $Z = 114$ (or 120) and $N = 184$ would dramatically increase the half-life. (c) The experimentally observed half-life of $^{289}$Fl ($Z = 114$) is about 2 seconds. What does this tell you about the strength of the shell stabilization?

Problem 4.30Alpha-decay Q-values from the SEMF. The Q-value for alpha decay of a nucleus $(Z, A)$ is $Q_\alpha = B(\alpha) + B(Z-2, A-4) - B(Z,A)$, where $B(\alpha) = 28.296$ MeV. (a) Derive an approximate expression for $Q_\alpha$ using the SEMF by expanding $B(Z-2, A-4)$ to first order around $(Z, A)$. (b) Show that $Q_\alpha$ generally increases with $Z$ for heavy nuclei. (c) Calculate $Q_\alpha$ for $^{238}$U using the full SEMF and compare to the experimental value $Q_\alpha = 4.270$ MeV. (d) The Geiger-Nuttall law relates the alpha-decay half-life to $Q_\alpha$. How sensitive is the predicted half-life to a 0.5 MeV uncertainty in $Q_\alpha$?


Section G: Estimation and Conceptual Problems

Problem 4.31Order-of-magnitude check. Without using a calculator, estimate $B/A$ for $^{100}$Ru ($Z = 44$, $A = 100$) by evaluating each SEMF term to the nearest 0.5 MeV/nucleon. Which terms can you neglect to within 0.5 MeV/nucleon accuracy?

Problem 4.32Design a nucleus. Suppose you could choose any combination of $Z$ and $N$ up to $A = 150$. Using the SEMF, what values of $Z$ and $N$ would give the maximum $B/A$? Is your answer a real nucleus? If so, look up its properties.

Problem 4.33The Coulomb barrier height. Using the SEMF's Coulomb term and the nuclear radius formula $R = 1.2 A^{1/3}$ fm, estimate the Coulomb barrier height (in MeV) for two identical $^{56}$Fe nuclei approaching each other. This is the minimum kinetic energy needed for iron-iron fusion. Compare this to the thermal energy $k_B T$ at the center of a massive star ($T \approx 5 \times 10^9$ K). Can stellar conditions fuse iron? Explain the implications.

Problem 4.34The proton-to-neutron ratio of the universe. In the early universe, the proton-to-neutron ratio was set by weak interaction freeze-out at about $n/p \approx 1/7$. (a) If all available neutrons were incorporated into $^4$He (as Big Bang nucleosynthesis approximately does), what fraction of the baryonic mass would be helium? (b) The SEMF shows that $^4$He is much more tightly bound per nucleon than $^2$H or $^3$He. Why does nucleosynthesis not proceed past $^4$He to build $^{12}$C and $^{56}$Fe? (Hint: consider the mass-5 and mass-8 gaps, and the low baryon density.)

Problem 4.35Neutron star crust. In the outer crust of a neutron star, nuclei exist in a lattice surrounded by a gas of electrons. The equilibrium nucleus at a given density is the one that minimizes the total energy (nuclear binding + electron kinetic energy + lattice Coulomb energy). (a) Using the SEMF, find the nucleus with the maximum binding energy per nucleon at $N/Z = 1.5$ (characteristic of neutron-rich conditions). (b) The actual outer crust composition at moderate densities is dominated by $^{56}$Fe (at low density) transitioning to increasingly neutron-rich nuclei ($^{62}$Ni, $^{64}$Ni, $^{66}$Ni, $^{86}$Kr, $^{84}$Se, ...) at higher density. Why does the equilibrium shift toward more neutron-rich nuclei as density increases? (Hint: consider the electron Fermi energy and electron capture.)


Solutions to Selected Problems

Problem 4.1 Solution:

(a–e) For $^{40}$Ca ($Z = 20$, $A = 40$, $N = 20$, even-even):

  • Volume: $15.75 \times 40 = 630.0$ MeV
  • Surface: $-17.80 \times 40^{2/3} = -17.80 \times 11.696 = -208.2$ MeV
  • Coulomb: $-0.711 \times 20 \times 19 / 40^{1/3} = -0.711 \times 380 / 3.420 = -78.9$ MeV
  • Asymmetry: $-23.7 \times (40 - 40)^2 / 40 = 0$ MeV (since $N = Z$)
  • Pairing: $+11.2 / 40^{1/2} = +11.2 / 6.325 = +1.8$ MeV

(f) Total: $B_{\text{SEMF}} = 630.0 - 208.2 - 78.9 + 0 + 1.8 = 344.7$ MeV

(g) Experimental: $B_{\text{exp}} = 342.052$ MeV

(h) Percentage error: $(344.7 - 342.1)/342.1 \times 100 = 0.8\%$

(i) The SEMF slightly overestimates the binding energy. Although $^{40}$Ca is doubly magic, which typically means extra binding beyond the SEMF, the close agreement here suggests that the SEMF parameters (which are fit to all nuclei) are partially absorbing the average shell effect. The residual pattern — comparing $^{40}$Ca to its neighbors — would reveal the shell closure more clearly than any single number.

Problem 4.12 Solution (partial):

(b) Values of $Z_{\text{stable}}$:

$A$ $Z_{\text{stable}}$ Known stable nucleus
20 10.0 $^{20}$Ne ($Z = 10$)
40 19.4 $^{40}$Ca ($Z = 20$)
80 35.5 $^{80}$Se ($Z = 34$)
120 50.5 $^{120}$Sn ($Z = 50$)
200 80.1 $^{200}$Hg ($Z = 80$)
250 96.6 (no stable nuclei at $A = 250$)

(d) $N/Z = (A - Z)/Z = A/Z - 1$. Setting $N/Z = 1.5$ gives $Z = A/2.5 = 0.4A$. From the stability formula: $Z = A / (2 + 0.015 A^{2/3})$, so $0.4A = A / (2 + 0.015 A^{2/3})$, giving $2 + 0.015 A^{2/3} = 2.5$, or $A^{2/3} = 33.3$, hence $A \approx 193$. So $N/Z > 1.5$ for $A \gtrsim 193$.

Problem 4.19 Solution (partial):

(a) For $^{56}$Fe: $B_{\text{SEMF}} = 495.7$ MeV (see worked example in Section 4.6), so $B/A = 8.85$ MeV.

For $^{62}$Ni ($Z = 28$, $A = 62$, even-even): - Volume: $15.75 \times 62 = 976.5$ MeV - Surface: $-17.80 \times 62^{2/3} = -17.80 \times 15.69 = -279.3$ MeV - Coulomb: $-0.711 \times 28 \times 27 / 62^{1/3} = -0.711 \times 756 / 3.958 = -135.8$ MeV - Asymmetry: $-23.7 \times (62 - 56)^2 / 62 = -23.7 \times 36/62 = -13.8$ MeV - Pairing: $+11.2 / 62^{1/2} = +1.4$ MeV - Total: $B_{\text{SEMF}} = 549.0$ MeV; $B/A = 8.85$ MeV

The SEMF gives essentially the same $B/A$ for both nuclei. The experimental values distinguish them: $B/A(^{62}\text{Ni}) = 8.795$ MeV vs. $B/A(^{56}\text{Fe}) = 8.790$ MeV. $^{62}$Ni is the most tightly bound nucleus per nucleon, but the difference is only 5 keV/nucleon — well within the SEMF's resolution.

Problem 4.8 Solution:

(a) $a_C = \frac{3}{5} \frac{e^2}{4\pi\epsilon_0 r_0} = \frac{3}{5} \times \frac{1.4400 \text{ MeV fm}}{1.20 \text{ fm}} = 0.600 \times 1.200 = 0.720$ MeV.

(b) The fitted value $a_C = 0.711$ MeV is about 1.3% lower than the electrostatic prediction. This is excellent agreement.

(c) The small difference can be attributed to: (i) the diffuse nuclear surface — the charge distribution is not a sharp-edged sphere but has a diffuseness $a \approx 0.5$ fm (Woods-Saxon form factor), which slightly lowers the Coulomb energy; (ii) the Coulomb exchange correction, which reduces the direct Coulomb energy by $\sim 3$–$5\%$; (iii) the finite size of the proton ($r_p \approx 0.87$ fm), which smears the charge distribution. The fact that these small corrections bring the theoretical prediction into agreement with the fit is a convincing validation of the uniform-sphere model.

Problem 4.15 Solution (partial):

(a) For $^{208}$Pb ($Z = 82$, $A = 208$, $N = 126$): $S_n = B(82, 208) - B(82, 207)$.

$B_{\text{SEMF}}(82, 208) = 1633.8$ MeV (from Section 4.6 calculation).

For $^{207}$Pb ($Z = 82$, $A = 207$, $N = 125$, odd-$A$): - Volume: $15.75 \times 207 = 3260.3$ MeV - Surface: $-17.80 \times 207^{2/3} = -17.80 \times 34.99 = -622.8$ MeV - Coulomb: $-0.711 \times 82 \times 81 / 207^{1/3} = -0.711 \times 6642 / 5.917 = -797.8$ MeV - Asymmetry: $-23.7 \times (207 - 164)^2 / 207 = -23.7 \times 1849 / 207 = -211.7$ MeV - Pairing: $0$ (odd-$A$) - Total: $B(82, 207)_{\text{SEMF}} = 1628.0$ MeV

$S_n(^{208}\text{Pb})_{\text{SEMF}} = 1633.8 - 1628.0 = 5.8$ MeV.

Experimental: $S_n(^{208}\text{Pb}) = 7.368$ MeV. The SEMF underestimates $S_n$ because it does not capture the extra binding from the $N = 126$ shell closure, which makes it especially costly to remove a neutron from the closed shell.

(b) The experimental drop from $S_n(^{208}\text{Pb}) = 7.368$ MeV to $S_n(^{209}\text{Pb}) = 3.937$ MeV — a difference of 3.4 MeV — is the direct signature of the $N = 126$ shell closure. The 127th neutron in $^{209}$Pb occupies the next shell above the $N = 126$ gap, and its separation energy is correspondingly smaller.

(c) The SEMF does NOT predict this dramatic drop. The SEMF's separation energy varies smoothly with $N$, changing by only about 0.2–0.3 MeV between consecutive isotopes. The 3.4 MeV discontinuity is a pure shell effect.

Problem 4.30 Solution (partial):

(c) For $^{238}$U ($Z = 92$, $A = 238$) decaying to $^{234}$Th ($Z = 90$, $A = 234$):

$Q_\alpha = B(^4\text{He}) + B(^{234}\text{Th}) - B(^{238}\text{U})$

$B_{\text{SEMF}}(^{238}\text{U}) = 15.75 \times 238 - 17.80 \times 238^{2/3} - 0.711 \times 92 \times 91 / 238^{1/3} - 23.7 \times (238-184)^2/238 + 11.2/\sqrt{238}$ $= 3748.5 - 683.7 - 961.5 - 290.5 + 0.7 = 1813.5$ MeV

$B_{\text{SEMF}}(^{234}\text{Th}) = 15.75 \times 234 - 17.80 \times 234^{2/3} - 0.711 \times 90 \times 89 / 234^{1/3} - 23.7 \times (234-180)^2/234 + 11.2/\sqrt{234}$ $= 3685.5 - 676.6 - 926.2 - 295.1 + 0.7 = 1788.3$ MeV

$Q_\alpha = 28.3 + 1788.3 - 1813.5 = 3.1$ MeV

The experimental value is $Q_\alpha = 4.270$ MeV. The SEMF underestimates by about 1.2 MeV — a significant discrepancy given that the Geiger-Nuttall law predicts an exponential dependence of half-life on $Q_\alpha$. A 1 MeV change in $Q_\alpha$ can change the predicted half-life by several orders of magnitude, illustrating why accurate mass predictions are essential for decay rate calculations.