Fermi's 1934 theory of beta decay - the first quantitative field theory of the weak interaction - explicitly left open what it could not do: the relativistic (order-v/c) terms of the coupling were omitted, the covariant form of the interaction was provisional, the forbidden-transition theory was only rough, the coupling constant was fixed to order of magnitude only, and, crucially, the nuclear matrix element governing every transition rate could not be computed from first principles because nuclear structure was unknown. The interaction-form, relativistic-structure, forbidden-transition, and coupling-constant questions were decisively resolved between 1935 and 1967 (Gamow-Teller transitions; parity violation; the V-A interaction; the electroweak Standard Model; the Fermi constant to part-per-million). What survives is the matrix-element question in modern form: ab initio nuclear theory now reproduces beta-decay rates and resolves the fifty-year quenching of the axial-vector coupling constant up to tin-100, but for heavy nuclei (A above roughly 100) - the candidates for neutrinoless double-beta decay and the nuclei of the r-process - and for forbidden transitions, no controlled, uncertainty-quantified first-principles prediction of the nuclear weak matrix elements exists; model spread remains a factor of two to three, and a 2022 community roadmap explicitly identifies the remaining work. The problem posed is to compute these matrix elements from the Standard Model weak interaction and QCD-based nuclear theory with certified uncertainties, validated against measured decay data.

Background

Radioactive beta decay, in which an atomic nucleus converts a neutron into a proton (or conversely) while emitting an electron and an antineutrino (or a positron and a neutrino), was the first weak-interaction process discovered, and its continuous energy spectrum required Pauli's hypothesis of the neutral, nearly massless neutrino. In 1934 Enrico Fermi proposed the first quantitative theory of the process, built in direct analogy with the quantum theory of light emission: just as an excited atom de-excites by creating a photon, a nucleus transforms a neutron into a proton by creating an electron-neutrino pair, coupled through a point-like four-fermion (contact) interaction whose strength constant plays the role that the electric charge plays in electromagnetism. Because nucleons in nuclei move at only about a tenth of the speed of light, Fermi treated the heavy particles nonrelativistically and kept only the leading 'time-component' term of the current-current coupling; the order-v/cv/c spatial terms were dropped (this is recorded explicitly in the immediately following papers of Konopinski and Uhlenbeck (1935) and Gamow and Teller (1936), which restate Fermi's ansatz). The decay probability then factorizes into a kinematic phase-space factor and a nuclear matrix element ∫vm∗un dτ\int v_m^{*}u_n\,d\tau - the overlap of the initial and final nuclear wave functions, entering exactly as the dipole matrix element enters radiation theory. Transitions in which this leading matrix element is nonzero are called allowed; when angular-momentum and parity selection rules force the emitted lepton pair to carry orbital angular momentum, the transition is forbidden and suppressed by powers of (nuclear radius)/(lepton wavelength), grouping the known emitters onto Sargent's empirical lifetime curves - the systematic bands formed when the known emitters are plotted by decay constant against endpoint energy.

Every gap Fermi left was decisively closed in the following decades except the matrix-element problem. Konopinski and Uhlenbeck (1935) tried a derivative-coupling modification, driven largely by the radium E (bismuth-210) spectrum, which was abandoned once the excess of very slow electrons in early heavy-element data proved to be an experimental artifact (Franklin 2005 analyzes this proposal-and-refutation episode). Gamow and Teller (1936) showed that including the nucleon spin in the coupling - an axial-vector operator alongside Fermi's polar-vector one - modifies the selection rules, so that both Fermi transitions (nuclear spin unchanged) and Gamow-Teller transitions (spin change up to one unit) exist; both are realized in nature. The general relativistic interaction was classified into five independent covariant types - scalar, vector, tensor, axial-vector, and pseudoscalar couplings - and Konopinski's 1943 review completed the quantitative theory of allowed and forbidden transitions (spectrum shape factors, comparative half-lives); by 1955 the predicted forbidden spectrum shapes had been found experimentally, confirming the theory. The interaction-form question itself was settled in 1956-1958: Lee and Yang proposed that parity (mirror-reflection symmetry) might be violated in weak processes; Wu, Ambler, Hayward, Hoppes, and Hudson observed parity violation in the beta decay of polarized cobalt-60; Goldhaber, Grodzins, and Sunyar measured the helicity of the neutrino (its spin points against its momentum - it is left-handed); and Feynman and Gell-Mann showed that a unique combination of vector and axial-vector couplings with maximal parity violation - the V-A interaction - accounts for weak processes with one universal coupling. In the modern Standard Model this contact interaction is the low-energy limit of the exchange of the charged W boson of the unified electroweak theory, and the coupling is precisely known: the Fermi constant GFG_F is measured to part-per-million precision from the muon lifetime (Webber et al. 2011), and the quark-sector coupling VudV_{ud} is determined to a few parts in ten thousand from superallowed 0+→0+0^+\to0^+ nuclear beta decays - fast transitions between isobaric analogue states (states of neighboring nuclei with the same mass number, related by neutron-proton interchange) that are pure Fermi transitions - with precision tests of the conserved-vector-current hypothesis, the statement that the weak vector current is not renormalized by the strong interaction (Hardy and Towner 2009).

What survives is the matrix-element question: the nuclear matrix element must still be computed rather than fitted. In modern language these are the matrix elements of the Fermi (charge) and Gamow-Teller (spin-isospin) weak operators - and of the higher-order operators governing forbidden transitions - between exact nuclear eigenstates, including two-body weak currents, in which the weak probe couples to a correlated nucleon pair. The modern first-principles program combines chiral effective field theory (the resolution-organized low-energy expansion of QCD, the strong-interaction theory of quarks and gluons, supplying both the nuclear Hamiltonian and the weak current at a consistent order) with ab initio many-body methods that solve the nuclear many-body problem with controlled, systematically improvable approximations - the no-core shell model (exact diagonalization in a basis that treats all nucleons as active), coupled cluster (building correlations systematically on a truncated reference state), and the in-medium similarity renormalization group (a continuous transformation of the Hamiltonian into a decoupled effective form). For light and medium nuclei up to tin-100 this program now reproduces measured beta-decay rates and has resolved the fifty-year 'quenching' puzzle (measured nuclear rates systematically smaller than free-neutron-based predictions, effectively an in-medium reduction of the axial-vector coupling gAg_A to roughly 0.75 of its free value; Gysbers et al. 2019). For the heavy, open-shell nuclei (those with partly filled nucleon shells and strong collective correlations) that dominate the most demanding applications, however, no controlled first-principles prediction exists. For neutrinoless double-beta decay (0νββ0\nu\beta\beta, the hypothetical simultaneous conversion of two neutrons into two protons with emission of two electrons and no neutrinos, whose observation would show that neutrinos are their own antiparticles and allow their mass scale to be extracted), nuclear matrix elements computed in different models still differ by factors of two to three (Engel and Menéndez 2017), and ab initio calculations have reached only the lightest candidates - calcium-48, germanium-76, selenium-82 (Yao et al. 2020; Belley et al. 2021). The same limitation affects first-forbidden beta decays of fission products, which dominate reactor antineutrino spectra (Estienne et al. 2019), and the beta-decay half-lives of the thousands of neutron-rich nuclei on the rapid neutron-capture (r-)process path of heavy-element synthesis in stars (Nabi et al. 2024).

Problem Statement

The surviving problem, stated in the language of the modern theory: compute the nuclear matrix elements that govern beta decay from first principles - that is, from the Standard Model weak interaction organized at low energies as chiral effective field theory, with the nuclear Hamiltonian and the weak current (including two-body currents and the associated in-medium renormalization of the axial-vector coupling gAg_A) treated consistently at the same order, solved with controlled many-body methods - for nuclei across the nuclear chart, with quantified uncertainties, and validate the predictions against measured decay data. Concretely, the objective comprises, each with a certified uncertainty budget (effective-field-theory truncation errors, many-body and basis convergence, statistical and input low-energy-constant errors): (i) allowed-transition matrix elements (Fermi and Gamow-Teller) and the derived half-lives and comparative half-lives (log ft values) for nuclei without measured data, in particular for heavy, open-shell, often deformed nuclei with mass number A≳100A\gtrsim100; (ii) the in-medium renormalization of the axial-vector coupling beyond the light- and medium-mass regime where it is already resolved; (iii) first- and higher-forbidden transition matrix elements and spectrum shape factors; and (iv) the second-order weak matrix elements for neutrinoless double-beta decay of the heavy candidates (tellurium-130, xenon-136, neodymium-150), for which current nuclear-model spread is a factor of two to three. These four subgoals are facets of one objective - the point of Fermi's original question - namely that beta-decay observables be predicted from theory rather than fitted, for each nucleus in the set fixed below, with an uncertainty band consistent with all available measured half-lives, ft values, and spectrum shapes.

The objective is the conjunction of (i)-(iv), each accompanied by a certified uncertainty budget - that is, a band obtained by propagating the error sources named above (effective-field-theory truncation, many-body and basis convergence, statistical error, and input low-energy-constant error) to the corresponding observable and consistent with the corresponding measured data within that band. The nuclear scope is finite, not class-wide: the 'across the nuclear chart' wording above describes the reach of the method, not a per-nucleus requirement. Subgoal (iv) covers the three explicitly named candidates (tellurium-130, xenon-136, neodymium-150); subgoals (i)-(iii) cover the corresponding results for a bounded set of representative heavy open-shell nuclei with A≳100A\gtrsim100 that the solver states explicitly, together with the selection rule used to choose them, subject to a fixed coverage floor: the set must include at least one nucleus from each of the two heavy regions that dominate the applications (the medium-heavy region around A≈100A\approx100-140 and the rare-earth and actinide region above A≈140A\approx140) and at least one spherical and one deformed case. Per-nucleus coverage of the entire chart is not required.

References

  1. E. Fermi, 'Versuch einer Theorie der β-Strahlen. I', Zeitschrift für Physik 88, 161-177 (1934). DOI: 10.1007/BF01351864, https://doi.org/10.1007/BF01351864
  2. Fred L. Wilson, 'Fermi's Theory of Beta Decay' (complete English translation of Fermi 1934), American Journal of Physics 36, 1150-1160 (1968). DOI: 10.1119/1.1974382, https://doi.org/10.1119/1.1974382
  3. E. J. Konopinski and G. E. Uhlenbeck, 'On the Fermi Theory of β-Radioactivity', Physical Review 48, 7-12 (1935). DOI: 10.1103/PhysRev.48.7, https://doi.org/10.1103/PhysRev.48.7
  4. G. Gamow and E. Teller, 'Selection Rules for the β-Disintegration', Physical Review 49, 895-899 (1936). DOI: 10.1103/PhysRev.49.895, https://doi.org/10.1103/PhysRev.49.895
  5. Emil Jan Konopinski, 'Beta-Decay', Reviews of Modern Physics 15, 209-245 (1943). DOI: 10.1103/RevModPhys.15.209, https://doi.org/10.1103/RevModPhys.15.209
  6. E. J. Konopinski, 'Fermi's Theory of Beta-Decay', Reviews of Modern Physics 27, 254-257 (1955). DOI: 10.1103/RevModPhys.27.254, https://doi.org/10.1103/RevModPhys.27.254
  7. Allan Franklin, 'The Konopinski-Uhlenbeck Theory of β Decay: Its Proposal and Refutation', in 'Wrong for the Right Reasons' (Archimedes, vol. 11), pp. 209-228, Springer (2005). DOI: 10.1007/1-4020-3048-7_10, https://doi.org/10.1007/1-4020-3048-7_10
  8. T. D. Lee and C. N. Yang, 'Question of Parity Conservation in Weak Interactions', Physical Review 104, 254-258 (1956). DOI: 10.1103/PhysRev.104.254, https://doi.org/10.1103/PhysRev.104.254
  9. C. S. Wu, E. Ambler, R. W. Hayward, D. D. Hoppes, and R. P. Hudson, 'Experimental Test of Parity Conservation in Beta Decay', Physical Review 105, 1413-1415 (1957). DOI: 10.1103/PhysRev.105.1413, https://doi.org/10.1103/PhysRev.105.1413
  10. M. Goldhaber, L. Grodzins, and A. W. Sunyar, 'Helicity of Neutrinos', Physical Review 109, 1015-1017 (1958). DOI: 10.1103/PhysRev.109.1015, https://doi.org/10.1103/PhysRev.109.1015
  11. R. P. Feynman and M. Gell-Mann, 'Theory of the Fermi Interaction', Physical Review 109, 193-198 (1958). DOI: 10.1103/PhysRev.109.193, https://doi.org/10.1103/PhysRev.109.193
  12. Steven Weinberg, 'A Model of Leptons', Physical Review Letters 19, 1264-1266 (1967). DOI: 10.1103/PhysRevLett.19.1264, https://doi.org/10.1103/PhysRevLett.19.1264
  13. J. C. Hardy and I. S. Towner, 'Superallowed 0+→0+0^+\to0^+ nuclear β decays: A new survey with precision tests of the conserved vector current hypothesis and the standard model', Physical Review C 79, 055502 (2009). DOI: 10.1103/PhysRevC.79.055502, https://doi.org/10.1103/PhysRevC.79.055502
  14. D. M. Webber, V. Tishchenko, Q. Peng, S. Battu, R. M. Carey, D. B. Chitwood, J. Crnkovic, P. T. Debevec, S. Dhamija, W. Earle, A. Gafarov, K. Giovanetti, T. P. Gorringe, F. E. Gray, Z. Hartwig, D. W. Hertzog, B. Johnson, P. Kammel, B. Kiburg, S. Kizilgul, J. Kunkle, B. Lauss, I. Logashenko, K. R. Lynch, R. McNabb, J. P. Miller, F. Mulhauser, C. J. G. Onderwater, J. Phillips, S. Rath, B. L. Roberts, P. Winter, and B. Wolfe (MuLan Collaboration), 'Measurement of the Positive Muon Lifetime and Determination of the Fermi Constant to Part-per-Million Precision', Physical Review Letters 106, 041803 (2011). DOI: 10.1103/PhysRevLett.106.041803, https://doi.org/10.1103/PhysRevLett.106.041803
  15. Jonathan Engel and Javier Menéndez, 'Status and future of nuclear matrix elements for neutrinoless double-beta decay: a review', Reports on Progress in Physics 80, 046301 (2017). DOI: 10.1088/1361-6633/aa5bc5, arXiv:1610.06548, https://arxiv.org/abs/1610.06548
  16. P. Gysbers, G. Hagen, J. D. Holt, G. R. Jansen, T. D. Morris, P. Navrátil, T. Papenbrock, S. Quaglioni, A. Schwenk, S. R. Stroberg, and K. A. Wendt, 'Discrepancy between experimental and theoretical β-decay rates resolved from first principles', Nature Physics 15, 428-431 (2019). DOI: 10.1038/s41567-019-0450-7, https://doi.org/10.1038/s41567-019-0450-7
  17. M. Estienne, M. Fallot, A. Algora, J. Briz-Monago, V. M. Bui, S. Cormon, W. Gelletly, L. Giot, V. Guadilla, D. Jordan, L. Le Meur, A. Porta, S. Rice, B. Rubio, J. L. Taín, E. Valencia, and A.-A. Zakari-Issoufou, 'Updated Summation Model: An Improved Agreement with the Daya Bay Antineutrino Fluxes', Physical Review Letters 123, 022502 (2019). DOI: 10.1103/PhysRevLett.123.022502, https://doi.org/10.1103/PhysRevLett.123.022502
  18. J. M. Yao, B. Bally, J. Engel, R. Wirth, T. R. Rodríguez, and H. Hergert, 'Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of 48Ca', Physical Review Letters 124, 232501 (2020). DOI: 10.1103/PhysRevLett.124.232501, https://doi.org/10.1103/PhysRevLett.124.232501
  19. A. Belley, C. G. Payne, S. R. Stroberg, T. Miyagi, and J. D. Holt, 'Ab Initio Neutrinoless Double-Beta Decay Matrix Elements for 48Ca, 76Ge, and 82Se', Physical Review Letters 126, 042502 (2021). DOI: 10.1103/PhysRevLett.126.042502, arXiv:2008.06588, https://arxiv.org/abs/2008.06588
  20. Mannque Rho, 'Multifarious Roles of Hidden Chiral-Scale Symmetry: Quenching gAg_A in Nuclei', Symmetry 13, 1388 (2021). DOI: 10.3390/sym13081388, arXiv:2012.07418, https://arxiv.org/abs/2012.07418
  21. V. Cirigliano, Z. Davoudi, J. Engel, R. J. Furnstahl, G. Hagen, U. Heinz, H. Hergert, M. Horoi, C. W. Johnson, A. Lovato, E. Mereghetti, W. Nazarewicz, A. Nicholson, T. Papenbrock, S. Pastore, M. Plumlee, D. R. Phillips, P. E. Shanahan, S. R. Stroberg, F. Viens, A. Walker-Loud, K. A. Wendt, and S. M. Wild, 'Towards Precise and Accurate Calculations of Neutrinoless Double-Beta Decay: Project Scoping Workshop Report', Journal of Physics G: Nuclear and Particle Physics 49, 120502 (2022). DOI: 10.1088/1361-6471/aca03e, arXiv:2207.01085, https://arxiv.org/abs/2207.01085
  22. Jameel-Un Nabi, Muhammad Riaz, and Arslan Mehmood, 'Re examination of the effect of Pairing Gaps on Gamow Teller Strength Distributions and beta decay Rates', Universe 10, 128 (2024). DOI: 10.3390/universe10030128, arXiv:2409.13264, https://doi.org/10.3390/universe10030128
  23. A. Belley, J. M. Yao, B. Bally, J. Pitcher, J. Engel, H. Hergert, J. D. Holt, T. Miyagi, T. R. Rodriguez, A. M. Romero, S. R. Stroberg, and X. Zhang, 'Ab initio uncertainty quantification of neutrinoless double-beta decay in 76Ge', Physical Review Letters 132, 182502 (2024). DOI: 10.1103/PhysRevLett.132.182502, arXiv:2308.15634, https://arxiv.org/abs/2308.15634