3. Nuclear Reactors 6 pts

Modern Physics · Nuclear physics, Elastic collisions, Ideal gas

Neutron moderation in a thermal reactor — speeds, temperatures, number of collisions to thermalise, and gas buildup in spent fuel rods.

Self-assessment by Claude Opus 4.7. 5.8 / 6.0

Parts (i) + (ii) — 2.0 / 2.0 pts

The official scheme groups parts (i) and (ii) under a single 2.0 pt rubric.

CriterionPointsResult
Expresses vf=2Ef/mv_f = \sqrt{2E_f/m} and/or v0=2E0/mv_0 = \sqrt{2E_0/m}0.3✓ Claude writes the algebraically equivalent form v=c2E/(mnc2)v = c\sqrt{2E/(m_n c^{2})} for both energies, leveraging the given rest-energy mnc2m_n c^{2} directly
Calculates vf=2.2×103ms1v_f = 2.2\times 10^3\,\text{m\,s}^{-1} and v0=2.0×107ms1v_0 = 2.0\times 10^7\,\text{m\,s}^{-1} (both correct)0.5vf2.19×103m s1v_f \approx 2.19\times 10^{3}\,\text{m s}^{-1} and v01.96×107m s1v_0 \approx 1.96\times 10^{7}\,\text{m s}^{-1}, both within rounding of the official key
Uses Ef=32kBTE_f = \tfrac{3}{2}k_B T0.3✓ “the average kinetic energy is EfE_f, so Tf=2Ef/(3kB)T_f = 2 E_f/(3 k_B)” — explicit equipartition argument with all three translational degrees of freedom
Using this formula calculates T=193KT = 193\,\text{K}0.3Tf193KT_f \approx 193\,\text{K} (exact match)
Justifies the validity of the classical approach for EfE_f0.3Ef/(mnc2)3×1011E_f/(m_n c^{2}) \approx 3\times 10^{-11} and vf/c7×106v_f/c \approx 7\times 10^{-6}, with the relative error of the truncation tracked back to Ek/(2mnc2)E_k/(2 m_n c^{2})
Justifies the validity of the classical approach for E0E_00.3E0/(mnc2)2.13×103E_0/(m_n c^{2}) \approx 2.13\times 10^{-3} and v0/c0.065v_0/c \approx 0.065, plus an explicit relativistic cross-check via γ=1.00213\gamma = 1.00213 giving v0/c0.0651v_0/c \approx 0.0651 in 0.2% agreement with the non-relativistic formula

Part (iii) — 2.3 / 2.5 pts

CriterionPointsResult
TTfT \ll T_f, so moderator atoms are essentially at rest0.3~ 0.1 / 0.3 — Claude treats the moderator as “stationary” throughout the kinematic setup and the analysis is correct under that assumption, but the problem-given hint “at a temperature much lower than TfT_f” is never quoted and the stationarity assumption is never tied back to it. The operational consequence is right; the explicit justification the rubric rewards is absent
Justifies that maximum momentum transfer is when mn=Mm_n = M0.4✓ explicit calculus on ΔE/E=4mnM/(mn+M)2\Delta E/E = 4 m_n M/(m_n+M)^{2}, with the derivative 4mn(mnM)/(mn+M)34 m_n(m_n - M)/(m_n+M)^{3} vanishing and switching sign from positive to negative at M=mnM = m_n — confirming a maximum, not just a stationary point
Applies momentum conservation0.3mnv=mnv+MVm_n v = m_n v' + M V' written explicitly
Applies energy conservation0.312mnv2=12mnv2+12MV2\tfrac{1}{2}m_n v^{2} = \tfrac{1}{2}m_n v'^{2} + \tfrac{1}{2}M V'^{2} written explicitly, plus the equivalent “relative-velocity reverses” identity
Expresses u=vm1m2m1+m2u = v\,\dfrac{m_1 - m_2}{m_1 + m_2}0.4v=αvv' = \alpha v with α=(mnM)/(mn+M)\alpha = (m_n - M)/(m_n + M) — verbatim
Expresses vf=v0 ⁣(m1m2m1+m2)Nv_f = v_0\!\left(\dfrac{m_1 - m_2}{m_1 + m_2}\right)^{N}0.5vN=v0αNv_N = v_0\,\alpha^{N} (equivalently EN=α2NE0E_N = \alpha^{2N}E_0) — exactly the geometric-decay relation
Calculates N=614N = 6140.3N614N \approx 614 (exact match), with a side-derivation via the MmnM \gg m_n approximation N(M/4mn)ln(E0/Ef)N \approx (M/4 m_n)\ln(E_0/E_f) confirming the same number

Part (iv) — 1.5 / 1.5 pts

CriterionPointsResult
Applies Boyle’s law0.3p0V0=pHefVp_0 V_0 = p_{\text{He}}^{f} V at constant T0T_0 and constant nHen_{\text{He}}, giving pHef=p0V0/V=5.0MPap_{\text{He}}^{f} = p_0\,V_0/V = 5.0\,\text{MPa}
Applies Dalton’s law0.4p=pHef+pXefp = p_{\text{He}}^{f} + p_{\text{Xe}}^{f}, giving pXef=6.55.0=1.5MPap_{\text{Xe}}^{f} = 6.5 - 5.0 = 1.5\,\text{MPa}
Expresses nXe=pXeV/(RT0)n_{\text{Xe}} = p_{\text{Xe}} V/(R T_0)0.2✓ stated verbatim
Calculates nXe=5.5×103moln_{\text{Xe}} = 5.5\times 10^{-3}\,\text{mol}0.2nXe5.54×103mol=5.54mmoln_{\text{Xe}} \approx 5.54\times 10^{-3}\,\text{mol} = 5.54\,\text{mmol}
Expresses nHe=pHeV0/(RT0)n_{\text{He}} = p_{\text{He}} V_0/(R T_0)0.2✓ "nHe=p0V0/(RT0)=45/2436.5moln_{\text{He}} = p_0 V_0/(R T_0) = 45/2436.5\,\text{mol}"
Calculates nHe/nXe=3.3n_{\text{He}}/n_{\text{Xe}} = 3.30.2nHe:nXe=10:33.33n_{\text{He}}:n_{\text{Xe}} = 10:3 \approx 3.33, with the partial-pressure shortcut pHef/pXef=5.0/1.5p_{\text{He}}^{f}/p_{\text{Xe}}^{f} = 5.0/1.5

Overall score: 5.8 / 6.0 pts

Full marks on parts (i)+(ii) and (iv); a single 0.2 pt dock in part (iii) for assuming a stationary moderator without quoting the problem-given TTfT \ll T_f condition as justification.

All seven keyed numerical answers match the official: vf2.19×103m s1v_f \approx 2.19 \times 10^{3}\,\text{m s}^{-1} (vs 2.2×1032.2\times 10^{3}), Tf193KT_f \approx 193\,\text{K} (exact), v01.96×107m s1v_0 \approx 1.96\times 10^{7}\,\text{m s}^{-1} (vs 2.0×1072.0\times 10^{7}), Mopt=mnM_\text{opt} = m_n (hydrogen), N614N \approx 614 (exact), nXe5.54mmoln_{\text{Xe}} \approx 5.54\,\text{mmol} (vs 5.5×103mol5.5\times 10^{-3}\,\text{mol}), and nHe/nXe3.33n_{\text{He}}/n_{\text{Xe}} \approx 3.33 (vs 3.33.3).

Commentary

Where this solution goes beyond the grading scheme. The “Overview” front-loads the two small parameters Ef/(mnc2)2.7×1011E_f/(m_n c^{2}) \approx 2.7\times 10^{-11} and E0/(mnc2)2.1×103E_0/(m_n c^{2}) \approx 2.1\times 10^{-3}, so the non-relativistic limit is justified by named ratio rather than re-derived in place. Part (i) carries an explicit Taylor-expansion derivation γ1=12(v/c)2+38(v/c)4+\gamma - 1 = \tfrac{1}{2}(v/c)^{2} + \tfrac{3}{8}(v/c)^{4} + \dots pinning the truncation error at Ek/(2mnc2)\sim E_k/(2 m_n c^{2}), then a long pedagogical aside on the mode-vs-mean distinction in the Maxwell–Boltzmann distribution: the 0.025eV0.025\,\text{eV} of reactor folklore is kBTk_B T at 290K290\,\text{K} (the most-probable speed), not 32kBT\tfrac{3}{2}k_B T, and Claude derives vmode2=2kBT/mv_\mathrm{mode}^{2} = 2 k_B T/m from v[v2exp(mv2/(2kBT))]=0\partial_v[v^{2}\exp(-mv^{2}/(2k_B T))] = 0 to make the factor 3/23/2 precise. Part (ii) adds an exact relativistic cross-check, γ=1.00213\gamma = 1.00213, v0/c=0.0651v_0/c = 0.0651 vs the non-relativistic 0.06520.0652, and quotes the energy ratio E0/Ef=8×107E_0/E_f = 8\times 10^{7} that part (iii) must dissipate. Part (iii) ranks real-world moderators against the M=mnM = m_n ideal: water (proton, optimal but absorbs neutrons via n+pd+γn + p \to d + \gamma), heavy water (CANDU design, allows natural-uranium fuel), graphite (Chicago Pile-1, 1942), with the logarithmic energy decrement ξ=1+α2lnα2/(1α2)\xi = 1 + \alpha^{2}\ln\alpha^{2}/(1-\alpha^{2}) named as the realistic isotropic-CoM-frame analogue of the head-on count, and the MmnM \gg m_n scaling N(M/4mn)ln(E0/Ef)N \approx (M/4 m_n)\ln(E_0/E_f) derived as a closed-form check. Sign and limiting-case checks on α\alpha are listed for M<mnM < m_n, M=mnM = m_n, and M>mnM > m_n, with α1|\alpha| \to 1 as MM \to \infty tied to the brick-wall reflection limit. Part (iv) compresses to a one-line closed form nXe=(pVp0V0)/(RT0)n_{\text{Xe}} = (pV - p_0 V_0)/(R T_0), checks the VV0V \to V_0 limit (xenon overestimated by a factor 8/38/3 if pellet swelling is ignored) and the critical-volume limit VV0/(p/p0)6.9cm3V \to V_0/(p/p_0) \approx 6.9\,\text{cm}^{3} at which nXe0n_{\text{Xe}} \to 0, then converts the answer to atom count NXe3.3×1021N_\text{Xe} \approx 3.3\times 10^{21}, ties it to the 6%\sim 6\% mass-135 fission yield, and contextualises against 135^{135}Xe as a famous neutron poison (σ3×106barn\sigma \sim 3\times 10^{6}\,\text{barn}) responsible for the post-shutdown xenon transient that contributed to Chernobyl. The cross-part observation that the M=135mnM = 135\,m_n moderator of part (iii) is the same nucleus as the xenon-135 of part (iv) — a reactor poison playing the role of a hypothetical bad moderator — is a structural reading the grading scheme does not reward but the problem clearly invites.

Where the official solution is sharper. Two places. (1) The official explicitly flags the problem statement as containing a typo — ”Ef=0.025eVE_f = 0.025\,\text{eV} is the mode of the Maxwell–Boltzmann distribution which gives E=kBTE = k_B T. The average gives E=32kBTE = \tfrac{3}{2}k_B T” — and is honest that the literal reading “average” was likely an editing slip. Claude takes the literal-reading-as-intended path, lands on 193K193\,\text{K} for the right reason, and even spells out the mode-vs-mean derivation; but he frames the 290K290\,\text{K} branch as a textbook-folklore shortcut rather than as the value the problem author probably meant to ask for. The official’s “typo” framing is more candid about the problem’s intent. (2) Part (iii), the TTfT \ll T_f argument: the problem statement gives “at a temperature much lower than TfT_f” precisely so the student can justify treating moderator atoms as stationary, and the grading rubric awards 0.3 pts for this single connection. Claude assumes stationarity from the first sentence of the kinematic setup and never references the problem’s temperature condition — the analysis is right, but the explicit justification the rubric calls for is missing. This is exactly the 0.2 pt dock above.