arxiv: 2601.05058 · v2 · submitted 2026-01-08 · ⚛️ nucl-th

Recognition: 1 theorem link

· Lean Theorem

Benchmarking projected generator coordinate method for nuclear Gamow-Teller transitions

R. N. Chen , X. Lian , J. M. Yao , C. L. Bai

Authors on Pith no claims yet

Pith reviewed 2026-05-16 16:20 UTC · model grok-4.3

classification ⚛️ nucl-th

keywords projected generator coordinate methodGamow-Teller transitionstwo-neutrino double beta decaynuclear matrix elementsshell model48Cafp shellIMSRG

0 comments

The pith

Extended projected generator coordinate method reproduces Gamow-Teller strengths and 2νββ matrix elements in agreement with exact shell-model results for calcium isotopes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends the projected generator coordinate method to calculate Gamow-Teller transition strengths in even-even nuclei and the nuclear matrix element for two-neutrino double-beta decay. Wave functions for the odd-odd intermediate nuclei are formed as superpositions of neutron and proton quasiparticle configurations placed on vacua with constrained average odd particle numbers, after which angular momentum and particle number are restored by projection. This construction is tested on a shell-model Hamiltonian restricted to the fp shell and benchmarked directly against exact diagonalization results for calcium and titanium isotopes, including the 48Ca to 48Ti decay. The same observables are also compared with configuration-interaction calculations that use varying particle-hole truncations and with IMSRG-evolved interactions. Readers care because reliable GT matrix elements determine beta-decay rates and enter neutrino-mass searches via double-beta decay experiments.

Core claim

Within the PGCM framework the wave functions of odd-odd nuclei are constructed as superpositions of neutron and proton quasiparticle configurations built on quasiparticle vacua constrained to have, on average, odd neutron and odd proton particle numbers; angular momentum and particle numbers are restored through projection. Using a shell-model Hamiltonian defined in the fp shell, the extended method yields GT transition strengths and the 2νββ NME for 48Ca that agree well with exact shell-model solutions.

What carries the argument

Superposition of neutron and proton quasiparticle configurations on particle-number-constrained vacua, followed by angular-momentum and particle-number projection.

Load-bearing premise

The selected quasiparticle configurations on constrained vacua, once projected, already contain the essential correlations that determine the Gamow-Teller matrix elements inside the fp shell.

What would settle it

A clear mismatch between the projected-GCM GT strengths or 2νββ NME and the exact shell-model values for any fp-shell calcium or titanium isotope would show that the extension misses required correlations.

Figures

Figures reproduced from arXiv: 2601.05058 by C. L. Bai, J. M. Yao, R. N. Chen, X. Lian.

**Figure 3.** Figure 3: FIG. 3. (Color online) (a) The distribution of GT [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

**Figure 2.** Figure 2: FIG. 2. (Color online) The same as Fig. 1, but for the GT [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. (Color online) The distribution of GT transition strength for [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗

**Figure 7.** Figure 7: FIG. 7. (Color online) Cumulative NME [PITH_FULL_IMAGE:figures/full_fig_p006_7.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6. (Color online) (a) The ground-state energy of [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗

read the original abstract

In this work, we aim to achieve a minimal extension of the quantum-number projected generator coordinate method (PGCM) to describe Gamow-Teller (GT) transition strengths in even-even nuclei and to compute the NME of $2\nu\beta\beta$ decay. Within the PGCM framework, the wave functions of odd-odd nuclei are constructed as superpositions of neutron and proton quasiparticle configurations built on quasiparticle vacua constrained to have, on average, odd neutron and odd proton particle numbers. The angular momentum and particle numbers associated with the underlying mean-field states are restored through projection techniques. Using a shell-model Hamiltonian defined in the $fp$ shell, we assess the validity of this approach by benchmarking GT transitions in calcium and titanium isotopes, as well as the $2\nu\beta\beta$ decay of $^{48}$Ca to $^{48}$Ti, against exact solutions. For comparison, we also confront our results with those obtained from configuration-interaction calculations employing different particle-hole truncation schemes, both with and without in-medium similarity renormalization group (IMSRG) evolution.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

PGCM extension for GT transitions in odd-odd nuclei benchmarks cleanly against exact fp-shell diagonalization for 48Ca GT strengths and 2νββ NME.

read the letter

The paper shows that a minimal extension of the projected generator coordinate method reproduces Gamow-Teller strengths in calcium and titanium isotopes and the 2νββ matrix element for 48Ca to good accuracy against exact shell-model results in the fp shell. They build the odd-odd wave functions as superpositions of neutron-proton quasiparticle configurations on particle-number-constrained vacua, then restore angular momentum and particle number by projection. This is the concrete new piece: the quasiparticle construction for the odd-odd sector and its direct test on GT transitions and the double-beta NME. The benchmarking is done against full diagonalization plus several truncated CI schemes and IMSRG-evolved versions, all with the same Hamiltonian. That setup lets any shortfall in captured correlations appear as a numerical mismatch, which is the right way to test the method internally. The agreement looks solid enough to support the claim that the chosen superpositions pick up the essential pieces without extra collective coordinates in this space. The main limitation is that the validation stays inside the fp shell, where exact solutions exist. The paper does not yet demonstrate how the approach scales to larger model spaces or heavier nuclei where it would actually be needed. Readers will also want the precise numerical differences rather than qualitative statements of agreement. This is useful work for anyone calculating beta-decay matrix elements or double-beta decay in medium-mass nuclei. It deserves a serious referee because the benchmark is direct and the construction is reproducible.

Referee Report

0 major / 4 minor

Summary. The manuscript benchmarks a minimal extension of the projected generator coordinate method (PGCM) for Gamow-Teller (GT) transition strengths in even-even nuclei and the 2νββ nuclear matrix element (NME). Odd-odd nuclei wave functions are constructed as superpositions of neutron and proton quasiparticle configurations on particle-number-constrained vacua, followed by angular-momentum and particle-number projection. Using an fp-shell Hamiltonian, results for Ca and Ti isotopes and the 48Ca 2νββ NME are compared to exact shell-model diagonalizations, configuration-interaction calculations with varying truncations, and IMSRG-evolved interactions.

Significance. If the reported agreement holds, the work provides a direct, parameter-free validation of the extended PGCM against exact solutions in a tractable model space. This strengthens confidence in the method for GT strengths and double-beta decay NMEs, as the benchmark tests the capture of essential correlations internally via comparison to exact diagonalization of the same Hamiltonian. Such controlled benchmarks are valuable for extending the approach to heavier systems where exact methods are unavailable.

minor comments (4)

§4 (results for GT strengths): provide explicit quantitative metrics (e.g., RMS deviation or average percentage difference) between PGCM and exact shell-model values rather than qualitative statements of 'good agreement'.
§3.2 (wave-function construction): specify the selection criterion and number of quasiparticle configurations included in the superposition to allow assessment of convergence.
Table 1 or equivalent (comparison to CI truncations): clarify the precise particle-hole truncation levels used in the CI calculations and whether the same effective interaction is employed throughout.
Figure 2 (2νββ NME): add error estimates or sensitivity analysis with respect to the choice of generator coordinates.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript, the recognition of its significance as a controlled benchmark, and the recommendation for minor revision. We are pleased that the work is viewed as strengthening in the extended PGCM approach for GT strengths and 2νββ NMEs.

Circularity Check

0 steps flagged

No significant circularity: external benchmarks against exact solutions

full rationale

The paper benchmarks the minimally extended PGCM (quasiparticle configurations on particle-number-constrained vacua with J and N projection) directly against exact fp-shell shell-model diagonalizations and independent CI/IMSRG results for GT strengths and the 2νββ NME of 48Ca, using the identical Hamiltonian. No equations or self-citations reduce the reported matrix elements to fitted parameters or prior results by construction; discrepancies would appear as measurable differences in the numerical comparisons. This is a self-contained validation against external exact solutions.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach rests on standard nuclear many-body assumptions rather than new postulates; the fp-shell Hamiltonian is taken as given for benchmarking purposes, and projection operators are standard.

axioms (2)

domain assumption The fp-shell effective Hamiltonian provides a faithful representation of the low-energy structure for Ca and Ti isotopes.
Used to generate both the PGCM states and the exact reference solutions.
standard math Angular-momentum and particle-number projections applied to quasiparticle vacua restore the correct quantum numbers without significant truncation artifacts in the GT operator matrix elements.
Core technical step invoked for constructing the final wave functions.

pith-pipeline@v0.9.0 · 5498 in / 1452 out tokens · 54988 ms · 2026-05-16T16:20:46.800664+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Within the PGCM framework, the wave functions of odd-odd nuclei are constructed as superpositions of neutron and proton quasiparticle configurations built on quasiparticle vacua constrained to have, on average, odd neutron and odd proton particle numbers. The angular momentum and particle numbers ... are restored through projection techniques.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages · 4 internal anchors

[1]

Langanke and G

K. Langanke and G. Martínez-Pinedo, Nuclear weak- interaction processes in stars, Rev. Mod. Phys.75, 819 (2003)

work page 2003
[2]

Fischer, G

T. Fischer, G. Guo, K. Langanke, G. Martinez-Pinedo, Y .-Z. Qian, and M.-R. Wu, Neutrinos and nucleosynthesis of elements, Prog. Part. Nucl. Phys.137, 104107 (2024), arXiv:2308.03962 [astro-ph.HE]

work page arXiv 2024
[3]

Suzuki, Nuclear weak rates and nuclear weak pro- cesses in stars, Prog

T. Suzuki, Nuclear weak rates and nuclear weak pro- cesses in stars, Prog. Part. Nucl. Phys.126, 103974 (2022), arXiv:2205.09262 [nucl-th]. 2 Note that the value ofM 2ν effdoes not depend on the choice of thegA quench- ing factorq, while theM 2ν does. For example, theM 2ν by the shell model calculation based on Eq.(20) is 0.0539/q2 ≃0.090 MeV −1

work page arXiv 2022
[4]

Hu, L.-J

Q.-Y . Hu, L.-J. Wang, and Y . Sun, Stellar weak rates of ther p- process waiting points: Effects of strong magnetic fields, Phys. Rev. Lett.135, 042702 (2025)

work page 2025
[5]

Herczeg, Beta decay beyond the standard model, Prog

P. Herczeg, Beta decay beyond the standard model, Prog. Part. Nucl. Phys.46, 413 (2001)

work page 2001
[6]

Severijns, M

N. Severijns, M. Beck, and O. Naviliat-Cuncic, Tests of the standard electroweak model in nuclear beta decay, Rev. Mod. Phys.78, 991 (2006)

work page 2006
[7]

E. W. Otten and C. Weinheimer, Neutrino mass limit from tri- tiumβdecay, Rep. Prog. Phys.71, 086201 (2008)

work page 2008
[8]

Falkowski, M

A. Falkowski, M. González-Alonso, and O. Naviliat-Cuncic, Comprehensive analysis of beta decays within and beyond the Standard Model, JHEP04, 126, arXiv:2010.13797 [hep-ph]. 8

work page arXiv 2010
[9]

F. T. Avignone, S. R. Elliott, and J. Engel, Double beta decay, majorana neutrinos, and neutrino mass, Rev. Mod. Phys.80, 481 (2008)

work page 2008
[10]

Engel and J

J. Engel and J. Menéndez, Status and future of nuclear matrix elements for neutrinoless double-beta decay: a review, Reports on Progress in Physics80, 046301 (2017)

work page 2017
[11]

Chu and H.-S

H. Ejiri, Neutrino-mass sensitivity and nuclear matrix element for neutrinoless double beta decay, Universe6, 10.3390/uni- verse6120225 (2020)

work page doi:10.3390/uni- 2020
[12]

J. C. Hardy and I. S. Towner, Superallowed 0 + →0 + nuclearβ decays: 2020 critical survey, with implications forV ud and ckm unitarity, Phys. Rev. C102, 045501 (2020)

work page 2020
[13]

Severijns, L

N. Severijns, L. Hayen, V . De Leebeeck, S. Vanlangendonck, K. Bodek, D. Rozpedzik, and I. S. Towner,Ftvalues of the mirrorβtransitions and the weak-magnetism-induced current in allowed nuclearβdecay, Phys. Rev. C107, 015502 (2023)

work page 2023
[14]

Agostini, G

M. Agostini, G. Benato, J. A. Detwiler, J. Menéndez, and F. Vissani, Toward the discovery of matter creation with neutrinolessββdecay, Rev. Mod. Phys.95, 025002 (2023), arXiv:2202.01787 [hep-ex]

work page arXiv 2023
[15]

J. M. Yao, J. Meng, Y . F. Niu, and P. Ring, Beyond-mean- field approaches for nuclear neutrinoless double beta decay in the standard mechanism, Prog. Part. Nucl. Phys.126, 103965 (2022), arXiv:2111.15543 [nucl-th]

work page arXiv 2022
[16]

Hergert, Front

H. Hergert, A Guided Tour ofab initioNuclear Many-Body Theory, Front. in Phys.8, 379 (2020), arXiv:2008.05061 [nucl- th]

work page arXiv 2020
[17]

Weinberg, Effective chiral Lagrangians for nucleon - pion interactions and nuclear forces, Nucl

S. Weinberg, Effective chiral Lagrangians for nucleon - pion interactions and nuclear forces, Nucl. Phys. B363, 3 (1991)

work page 1991
[18]

S. R. Stroberg, S. K. Bogner, H. Hergert, and J. D. Holt, Nonempirical Interactions for the Nuclear Shell Model: An Update, Ann. Rev. Nucl. Part. Sci.69, 307 (2019), arXiv:1902.06154 [nucl-th]

work page arXiv 2019
[19]

J. M. Yao, J. Engel, L. J. Wang, C. F. Jiao, and H. Hergert, Generator-coordinate reference states for spectra and 0νββde- cay in the in-medium similarity renormalization group, Phys. Rev. C98, 054311 (2018), arXiv:1807.11053 [nucl-th]

work page internal anchor Pith review Pith/arXiv arXiv 2018
[20]

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 of48Ca, Phys. Rev. Lett.124, 232501 (2020)

work page 2020
[21]

Belley, C

A. Belley, C. G. Payne, S. R. Stroberg, T. Miyagi, and J. D. Holt,AbInitioNeutrinoless Double-Beta Decay Matrix Ele- ments for 48Ca , 76Ge , and 82Se, Phys. Rev. Lett.126, 042502 (2021), arXiv:2008.06588 [nucl-th]

work page arXiv 2021
[22]

Belley, J

A. Belley, J. M. Yao, B. Bally, J. Pitcher, J. Engel, H. Hergert, J. D. Holt, T. Miyagi, T. R. Rodríguez, A. M. Romero, S. R. Stroberg, and X. Zhang, Ab initio uncertainty quantification of neutrinoless double-beta decay in 76Ge, Phys. Rev. Lett.132, 182502 (2024)

work page 2024
[23]

Belley, T

A. Belley, T. Miyagi, S. R. Stroberg, and J. D. Holt, Ab ini- tio calculations of neutrinolessββdecay refine neutrino mass limits, (2023), arXiv:2307.15156 [nucl-th]

work page arXiv 2023
[24]

Hergert, S

H. Hergert, S. K. Bogner, T. D. Morris, S. Binder, A. Calci, J. Langhammer, and R. Roth, Ab initio multireference in- medium similarity renormalization group calculations of even calcium and nickel isotopes, Phys. Rev. C90, 041302 (2014)

work page 2014
[25]

Ring and P

P. Ring and P. Schuck,The nuclear many-body problem (Springer-Verlag, New York, 1980)

work page 1980
[26]

S. M. Lauber, H. C. Frye, and C. W. Johnson, Benchmarking angular-momentum projected Hartree–Fock as an approxima- tion, J. Phys. G48, 095107 (2021), arXiv:2104.03455 [nucl-th]

work page arXiv 2021
[27]

Gamov-Teller transitions from 14N ground to 14C ground and excited states

Y . Kanada-En’yo and T. Suhara, Gamow-Teller transitions from the 14N ground state to the 14C ground and excited states, Phys. Rev. C89, 044313 (2014), arXiv:1401.5517 [nucl-th]

work page internal anchor Pith review Pith/arXiv arXiv 2014
[28]

Konieczka, P

M. Konieczka, P. B˛ aczyk, and W. Satuła,β-decay study within multireference density functional theory and beyond, Phys. Rev. C93, 042501 (2016)

work page 2016
[29]

Morita and Y

H. Morita and Y . Kanada-En’yo, Low-energy gamow-teller transitions in deformedn=zodd-odd nuclei, Phys. Rev. C98, 034307 (2018)

work page 2018
[30]

Mi ´skiewicz, M

J. Mi ´skiewicz, M. Konieczka, and W. Satuła, Two-neutrino 0+→0+double-βdecay of Ca48 within the density-functional- theory–based no-core configuration-interaction framework, Phys. Rev. C112, 055502 (2025), arXiv:2506.13747 [nucl-th]

work page arXiv 2025
[31]

The Shell Model as Unified View of Nuclear Structure

E. Caurier, G. Martinez-Pinedo, F. Nowacki, A. Poves, and A. P. Zuker, The Shell Model as Unified View of Nuclear Structure, Rev. Mod. Phys.77, 427 (2005), arXiv:nucl-th/0402046

work page internal anchor Pith review Pith/arXiv arXiv 2005
[32]

Honma, T

M. Honma, T. Otsuka, B. A. Brown, and T. Mizusaki, New ef- fective interaction forp f-shell nuclei and its implications for the stability of then=z=28 closed core, Phys. Rev. C69, 034335 (2004)

work page 2004
[33]

The In-Medium Similarity Renormalization Group: A Novel Ab Initio Method for Nuclei

H. Hergert, S. K. Bogner, T. D. Morris, A. Schwenk, and K. Tsukiyama, The In-Medium Similarity Renormalization Group: A Novel Ab Initio Method for Nuclei, Phys. Rept.621, 165 (2016), arXiv:1512.06956 [nucl-th]

work page internal anchor Pith review Pith/arXiv arXiv 2016
[34]

Barabash, Precise Half-Life Values for Two-Neutrino Double-βDecay: 2020 Review, Universe6, 159 (2020), arXiv:2009.14451 [nucl-ex]

A. Barabash, Precise Half-Life Values for Two-Neutrino Double-βDecay: 2020 Review, Universe6, 159 (2020), arXiv:2009.14451 [nucl-ex]

work page arXiv 2020
[35]

Kotila and F

J. Kotila and F. Iachello, Phase-space factors for double-βde- cay, Phys. Rev. C85, 034316 (2012)

work page 2012
[36]

Horoi, S

M. Horoi, S. Stoica, and B. A. Brown, Shell-model calculations of two-neutrino double-βdecay rates of 48Ca with the gxpf1a interaction, Phys. Rev. C75, 034303 (2007)

work page 2007