{"id":"0496136b-3b0b-4d85-ad9c-136d1987e902","arxiv_id":"1908.01873","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A vibration-driven generator coordinate method, using quasiparticle Tamm-Dancoff modes as coordinates, brings GCM 0νββ matrix elements and 4+ spectra closer to shell-model benchmarks.","lead":"This paper introduces a generator-coordinate method that builds its basis from the lowest nuclear vibration modes computed at the mean-field level, then projects out good quantum numbers. Tested on six double-beta decay nuclei against exact shell model results, it reduces the overestimate of decay matrix elements seen in standard GCM while improving excited-state spectra.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed NME improvement is not yet shown robust: no λ or mode-number sensitivity study, and the remaining QTDA-to-SM gap is comparable to the claimed CHFB-to-QTDA improvement.","rationale":"The paper is a credible proof-of-principle: it gives a concrete construction of QTDA-generated reference states via the generalized Thouless transformation (Eqs. 8–15), benchmarks against exact shell-model diagonalization in the same model space, and reports improvement in 4+ energies and total 0νββ NMEs for all three tested nuclei. Those are real, reproducible-sounding results. My stress-test pass does not find an internal inconsistency in the formalism; the generalized Thouless step is standard and the projected Hill-Wheeler diagonalization is well posed. The load-bearing weakness is exactly the one the reader identified: the method's success is demonstrated at a single, non-physically-calibrated amplitude λ=1 and a fixed 15-mode truncation, with no sensitivity or convergence evidence. This matters more than usual because the claimed NME improvement (0.13–0.39) is comparable to the remaining disagreement with the shell model (0.20–0.45), and because the QTDA-GCM ground-state energies are no better than the CHFB-GCM ones. The Fermi matrix elements even worsen, so the total-NME improvement is partly a cancellation between GT and F errors. The paper's own closing remarks concede that adding more reference states or combining with CHFB states may be necessary. None of this falsifies the method; it means the central claim is conditional on the reference-space construction being in a stable regime, which has not been shown. A λ-sweep and mode-number test is cheap relative to the 3D projection cost already paid and would settle whether the improvement is systematic. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":10736,"tokens_out":11084,"duration_ms":141824,"concrete_test":"Rerun the QTDA-GCM calculation for 130Te (or 136Xe) with the same SVD Hamiltonian and projection machinery, varying the amplitude over λ = 0.25, 0.5, 1.0, 1.5, 2.0, 3.0 and the number of QTDA modes over 5, 10, 15, 20, 25, 30, recording total M0ν and the 4+ excitation energy for each run. If M0ν changes by more than 0.10–0.15 (the magnitude of the claimed CHFB-to-QTDA improvement in 130Te) or the 4+ energy shift reverses sign across the λ range, the headline improvement is not robust; if M0ν is flat to within ~0.05 across λ ≥ 1 and mode counts ≥ 15, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—QTDA-driven GCM brings 0νββ NMEs closer to shell-model values—depends on the stability of a reference space fixed by two undocumented choices: the single amplitude λ=1 in Eq. (8) for every QTDA mode, and the truncation to the 15 lowest QTDA phonons, described only as 'a computationally tractable choice.' The QTDA eigenvectors carry no intrinsic physical scale; their normalization is set by the diagonalization, so λ=1 is a hand-set displacement with no physical criterion. In Table III the remaining QTDA-GCM-to-SM gap in M0ν is 0.38 (124Sn), 0.45 (130Te), and 0.20 (136Xe), while the claimed improvement over CHFB-GCM is 0.23, 0.13, and 0.39; the residual error is therefore comparable to or larger than the effect being claimed. Moreover, the Fermi contributions move away from SM for all three nuclei, and Table I shows QTDA-GCM ground-state energies are tens to a few hundred keV worse than CHFB-GCM. The paper itself flags that combining with CHFB reference states or adding more states 'may be a useful or even necessary strategy.' Without a λ sweep and a mode-number growth test, the improved NMEs and 4+ energies could be coincidental features of a small, heuristically chosen subspace rather than a systematic property of the method.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new variant of the generator-coordinate method (GCM) called QTDA-driven GCM. Starting from a Hartree-Fock-Bogoliubov vacuum, the authors solve the quasiparticle Tamm-Dancoff approximation (QTDA) to obtain low-lying two-quasiparticle vibrational modes, then use Thouless' theorem with a fixed amplitude lambda=1 to generate non-orthogonal reference states. These states are projected onto good angular momentum and particle number and mixed via the Hill-Wheeler equation. The method is benchmarked against standard constrained HFB (CHFB) GCM and full shell-model (SM) diagonalization for the 0nubb candidate nuclei 124Sn/Te, 130Te/Xe, and 136Xe/Ba, all using the same SVD Hamiltonian. The authors report that QTDA-driven GCM lowers the excitation energies of 4+ states relative to CHFB-GCM and brings the 0nubb Gamow-Teller and total matrix elements closer to the SM values, although ground-state energies and Fermi matrix elements are not improved.","tokens_in":11069,"tokens_out":3837,"duration_ms":41331,"significance":"If the central claim is robust, the method is significant: it offers a deterministic, Hamiltonian-driven way to incorporate vibrational and broken-pair correlations into GCM without hand-selected constraint fields, addressing a known limitation of standard GCM for spherical and weakly deformed nuclei. The benchmarking is honest and well controlled: the same SVD Hamiltonian is used for CHFB-GCM, QTDA-GCM, and exact SM, and the comparison against numerically exact SM results provides a strong test bed. The paper also clearly identifies the new ingredient relative to prior two-quasiparticle GCM work (collective QTDA superpositions rather than individual two-quasiparticle excitations). However, the central quantitative improvement is not yet shown to be systematic: the key choices of lambda=1 and 15 QTDA modes are heuristic, no sensitivity study is provided, and the residual discrepancies with the SM are comparable to the claimed improvements. The method is promising but the evidence presented is preliminary.","major_comments":[{"comment":"The choice lambda=1 for every QTDA mode is not justified by any physical criterion. Because the QTDA eigenvectors obtained from Eq. (7) determine only a ray and carry no intrinsic scale, lambda=1 is a hand-set displacement in the Thouless evolution. The paper states only that 'a single value of lambda = 1 worked sufficiently well.' No lambda-sensitivity study is shown, so the claimed improvements over CHFB-GCM in Tables I and III could be an artifact of this tuning. A sweep over lambda and a demonstration that the projected observables are stable (or converge) is needed before the improvement can be attributed to the method rather than to the particular amplitude.","section":"Eq. (8) and surrounding text"},{"comment":"The truncation to the 15 lowest-energy QTDA phonons is described only as 'a computationally tractable choice.' No test is given of convergence with respect to the number of QTDA modes. In Table III the remaining QTDA-GCM-to-SM gap in M0nu is 0.38 (124Sn), 0.45 (130Te), and 0.20 (136Xe), while the claimed improvement over CHFB-GCM is 0.23, 0.13, and 0.39; the residual error is comparable to or larger than the effect being claimed. Without a mode-number growth test, the improved NMEs and 4+ energies could be coincidental features of a small, heuristically chosen subspace rather than a systematic property of the method.","section":"Section on reference-state selection"},{"comment":"The Fermi contribution M_F moves away from the SM value for all three nuclei under QTDA-GCM: for 124Sn, CHFB-GCM gives -0.51, QTDA-GCM -0.73, and SM -0.47; for 130Te, -0.47, -0.69, -0.44; for 136Xe, -0.32, -0.50, -0.40. Thus the improvement in the total M0nu is entirely due to the Gamow-Teller component, and the Fermi part is made worse. The abstract's claim of 'significant improvement in 0νββ nuclear matrix elements' should be qualified to specify that the improvement is not uniform across the GT, Fermi, and tensor components.","section":"Table III"},{"comment":"The QTDA-GCM ground-state energies are 49-441 keV above the CHFB-GCM values and up to 1.8 MeV above the SM values (e.g., 124Te: QTDA-GCM -22.641 MeV, CHFB-GCM -23.082 MeV, SM -24.446 MeV). Since the Hill-Wheeler diagonalization of Eq. (17) is variational with respect to the chosen reference space, this indicates that the 15 QTDA-generated states miss some correlations that are captured by the CHFB basis. The paper acknowledges this and suggests that combining the two bases 'may be a useful or even necessary strategy,' but this undercuts the present claim that QTDA-driven GCM is competitive with or superior to CHFB-GCM for ground-state properties and reinforces the need for a systematic study of the reference space.","section":"Table I"}],"minor_comments":[{"comment":"The text contains several typographical errors: 'Tamm-Damcoff' should be 'Tamm-Dancoff'; 'we need to file new quasiparticle states' should be 'fill'; 'as these nuclides have axially symmetric HFB minimua' should be 'minima'; and 'the jj55-shell configuration space that compromises' should be 'comprises.'","section":"Section II"},{"comment":"The caption lists comparisons only with CHFB-GCM and SM, but the text states that the results are also in 'reasonable agreement with ... experimental spectra.' If experimental levels are shown or intended, they should be identified in the caption; otherwise the sentence should be revised.","section":"Figure 1 caption"},{"comment":"The sentence 'While the underestimation, also mentioned in Ref. [44], suggests the effective charges should be adjusted' is slightly unclear: it would be useful to state explicitly that the comparisons in this paper are intended as method benchmarks, not as fits to experiment, so the effective-charge issue does not affect the conclusions.","section":"Table II and text"}],"recommendation":"major_revision","confidential_remarks":"This is a promising methods paper with an honest and well-controlled benchmark, but the central quantitative claim is not yet robust to the two main heuristic choices (lambda=1 and 15 modes). The lack of any sensitivity or convergence study is the main obstacle; the paper is suitable for this journal if that gap is addressed. I do not see grounds for rejection, but the current version overstates the significance of the NME improvement given that the Fermi component worsens and residual errors are comparable to the claimed effect."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper proposes a genuinely new GCM construction—use collective QTDA (quasiparticle Tamm-Dancoff) modes, rather than hand-picked constrained fields or individual two-quasiparticle states, to generate reference states via Thouless' theorem, then project onto good angular momentum and particle number and mix. Applied to 0νββ decay candidates, it moves the 4+ states and the 0νββ matrix elements toward exact shell-model results compared with the earlier CHFB-GCM. That is a real, useful step, and the benchmarking is honest: same SVD Hamiltonian, same projection machinery, and direct comparison to numerically exact shell-model diagonalization. The tables are clear and the paper does not oversell what it cannot show.\n\nWhat is actually new: prior work used only the lowest individual two-quasiparticle excitations; here the collective QTDA superpositions are the generator coordinates, and full three-dimensional angular momentum projection is included. That combination is not in the cited literature. The connection to RPA is acknowledged, and the credit to prior MCSM and Thouless-based thinking is fair.\n\nSoft spots, in proportion. The central weakness is the unexamined parameter choice. λ=1 in Eq. (8) is justified only as \"worked sufficiently well,\" and the 15 lowest QTDA modes are selected as \"computationally tractable.\" There is no λ sweep and no mode-number growth test. Since the QTDA eigenvectors carry no intrinsic scale, λ=1 is a hand-set displacement, and the 15-mode cutoff is a heuristic. The stress-test concern is on target: for 124Sn and 130Te the remaining QTDA-to-SM gap is larger than the claimed CHFB-to-QTDA improvement, and the Fermi contributions move further from SM in all three nuclei. The ground-state energies are also somewhat worse than CHFB-GCM (50–350 keV). None of this sinks the paper—the direction of the NME shift is consistent, and the authors explicitly flag that combining with CHFB states or adding more reference states \"may be a useful or even necessary strategy.\" But it means the systematicity of the improvement is unproven, not established.\n\nWho this is for: people working on generator-coordinate methods, nuclear structure near closed shells, and 0νββ matrix elements. It deserves a serious referee, and the right referee will ask for a sensitivity study of λ and mode number. I would engage with it, and I would cite it as a method contribution, but I would not yet treat the NME improvement as a robust quantitative result. Recommend send to review, with the caveat that a revision should include those sensitivity tests.","headline":"A promising GCM variant that moves 0νββ NMEs toward shell-model values, but the heuristic λ=1 and 15-mode truncation need sensitivity tests before the improvement can be called systematic.","tokens_in":11577,"tokens_out":1535,"would_cite":true,"duration_ms":18093,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["21.60.Jz","23.40.Hc","27.60.+j"],"model":"deepseek-v4-flash","headline":"A generator-coordinate method driven by low-lying quasiparticle vibrational modes brings neutrinoless double-beta decay matrix elements closer to shell-model values.","keywords":["generator coordinate method","quasiparticle Tamm-Dancoff approximation","Thouless theorem","neutrinoless double-beta decay","angular momentum projection","particle number projection","shell-model benchmark","Hartree-Fock-Bogoliubov"],"falsifier":"Run the same calculation with several λ values (say 0.5, 1, and 2) and with 5, 15, and 30 QTDA modes and compare against the shell-model results for 124Sn, 130Te, and 136Xe; if the 4+ energies or 0νββ matrix elements move by more than the CHFB-to-shell-model gap across those choices, the reported improvement is not robust.","tokens_in":10536,"feed_emoji":"⚛️","tokens_out":10248,"duration_ms":98045,"temperature":0.7,"pith_summary":"This paper proposes a new way to build generator-coordinate bases for nuclear many-body calculations: instead of hand-picking external constraint fields, solve the quasiparticle Tamm-Dancoff (QTDA) equations on a Hartree-Fock-Bogoliubov ground state and use the low-lying vibrational modes as generator coordinates via Thouless' theorem. The reference states are then projected onto good angular momentum and particle number and mixed through the Hill-Wheeler equation. For six nuclei in neutrinoless double-beta decay candidate chains, benchmarked against full shell-model diagonalization with the same interaction, the method lowers 4+ excitation energies that standard constrained HFB-GCM overestimates and brings the 0νββ matrix elements of 124Sn, 130Te, and 136Xe closer to the shell-model values. If the pattern holds, this offers a systematic route to the correlations needed for reliable 0νββ matrix elements without relying on intuition about which collective coordinates matter.","feed_headline":"Vibrational modes narrow the double-beta matrix-element gap","feed_subtitle":"Low-lying QTDA phonons as generator coordinates fix 4+ states and bring 0νββ matrix elements closer to shell model.","key_machinery":"The engine is the QTDA phonon, a collective superposition of two-quasiparticle excitations obtained by diagonalizing the Hamiltonian in the two-quasiparticle space built on the HFB vacuum, interpreted as a small-amplitude vibrational mode of the energy landscape. Thouless' theorem, exp(λ Z_r)|Φ0>, turns such a one-body operator into a new non-orthogonal quasiparticle vacuum, with a Cholesky step restoring the fermion anticommutation relations. Fifteen low-lying QTDA phonons plus the HFB vacuum, after projection to good angular momentum and particle number and solution of the Hill-Wheeler equation, form the GCM basis. The fixed amplitude λ=1 appears in Eq. (8) and is the one tuning parameter used in lieu of sampling the amplitude.","core_discovery":"The central claim is that the correlations missing from a standard constrained Hartree-Fock-Bogoliubov generator-coordinate calculation—vibrational motion and broken-pair excitations—can be supplied by taking the lowest quasiparticle Tamm-Dancoff modes as generator coordinates. Exponentiating each QTDA operator on the HFB vacuum produces a new non-orthogonal vacuum; after angular-momentum and particle-number projection, the Hill-Wheeler diagonalization mixes these states. Using 15 low-lying QTDA modes at a single amplitude λ=1, the authors find that for the weakly deformed nuclei 124Sn, 124Te, 130Te, 130Xe, 136Xe, and 136Ba the 4+ states come down toward the shell model, and for the three 0νββ emitters the total matrix element shrinks from CHFB-GCM values of 2.76, 2.52, and 2.35 to 2.53, 2.39, and 1.96, against shell-model values of 2.15, 1.94, and 1.76. The improvement shows up mainly in the Gamow-Teller part, while ground-state energies stay 50 to 350 keV above the CHFB-GCM values.","pith_inferences":["The fixed amplitude λ=1 and the 15-mode cutoff are the two tuning dials; a convergence scan over λ and mode number would show whether the shell-model agreement is systematic or accidental.","The pattern in Table III—Gamow-Teller moves toward the shell model while the small Fermi component moves slightly away—suggests the method is capturing spin-isospin-like correlations, not all missing physics; a QTDA-plus-CHFB combined basis could target both.","For candidates with triaxial or strongly deformed minima, QTDA modes with K≠0 will mix angular momentum projections, so the method's full 3D projection machinery becomes essential; 76Ge is a natural test case.","The same basis-generation idea could be carried from QTDA to QRPA operators, which would put two-particle-two-hole correlations into the reference states and might cure the remaining ground-state-energy deficit; the paper notes this as future work."],"forward_implications":["Near-spherical and weakly deformed nuclei get lower 4+ excitation energies than CHFB-GCM provides, so the method captures the vibrational and broken-pair correlations that push those states down.","The 0νββ matrix elements of 124Sn, 130Te, and 136Xe are reduced relative to CHFB-GCM and move toward full shell-model values, mainly through the Gamow-Teller contribution, which is the part a neutrino experiment's interpretation depends on.","2+ energies and B(E2; 0+→2+) strengths stay close to shell-model and adopted experimental values, so adding QTDA reference states does not spoil the quadrupole collectivity the standard GCM already describes.","Because the QTDA modes are generated from the Hamiltonian itself, the method gives a systematic way to include non-collective and pairing correlations in GCM without guessing constraint fields.","Ground-state energies are slightly worse than CHFB-GCM, so combining QTDA-generated states with the usual constrained HFB states is a natural next step."],"supporting_citations":[{"why":"Supplies the Thouless theorem, the QTDA equations, projection operators, and the Hill-Wheeler framework on which the whole method rests.","marker":"[8]"},{"why":"Provides the CHFB-GCM baseline results this work compares against and documents the roughly 30% overestimate of 0νββ matrix elements that motivates the new correlations.","marker":"[35]"},{"why":"Provide the full shell-model energies, B(E2) values, and 0νββ matrix elements used as the exact benchmarks.","marker":"[44, 45]"},{"why":"Defines the SVD shell-model Hamiltonian that fixes the model space and interaction shared by all three calculations.","marker":"[46]"},{"why":"Gives the quasiparticle transformation and Cholesky construction that turns a Thouless-evolved state into a proper quasiparticle vacuum.","marker":"[43]"},{"why":"Defines the 0νββ transition operator and the CD-Bonn short-range correlation parametrization used in the matrix-element calculation.","marker":"[31]"},{"why":"The Monte Carlo shell model's use of Thouless' theorem to explore the energy landscape inspired the deterministic QTDA-mode generation.","marker":"[17, 18]"},{"why":"Supply the angular-momentum projection and natural-state diagonalization procedures used to build and solve the projected GCM kernels.","marker":"[22, 23]"}],"fun_headline_variants":["QTDA phonons as generator coordinates sharpen 0νββ matrix elements","Vibrational generator coordinates bring double-beta matrix elements in line","Low-lying phonon modes improve double-beta decay predictions","Mixing rotational and vibrational modes tightens neutrinoless double-beta numbers","Vibrational modes added to GCM shrink 0νββ matrix-element gap"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the 15 lowest QTDA modes, each evolved by the same fixed amplitude λ=1, span the correlations missing from the CHFB-GCM basis; if that subspace is too small or the amplitude is mis-scaled, the agreement with the shell model could be coincidental.","fun_headline_variants_meta":{"raw":{"variants":["QTDA phonons as generator coordinates sharpen 0νββ matrix elements","Vibrational generator coordinates bring double-beta matrix elements in line","Low-lying phonon modes improve double-beta decay predictions","Mixing rotational and vibrational modes tightens neutrinoless double-beta numbers","Vibrational modes added to GCM shrink 0νββ matrix-element gap"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000815,"raw_usage":{"total_tokens":3571,"prompt_tokens":946,"completion_tokens":2625,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":2543}},"tokens_in":562,"tokens_out":2625,"duration_ms":19057,"temperature":1.0,"reasoning_tokens":2543,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:00:47.842247+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same calculation with several λ values (say 0.5, 1, and 2) and with 5, 15, and 30 QTDA modes and compare against the shell-model results for 124Sn, 130Te, and 136Xe; if the 4+ energies or 0νββ matrix elements move by more than the CHFB-to-shell-model gap across those choices, the reported improvement is not robust.","supporting_citations":[{"cited_title":"Neacsu and M","cited_arxiv_id":null,"evidence_quote":"Defines the SVD shell-model Hamiltonian that fixes the model space and interaction shared by all three calculations."},{"cited_title":"Chen and J","cited_arxiv_id":null,"evidence_quote":"Gives the quasiparticle transformation and Cholesky construction that turns a Thouless-evolved state into a proper quasiparticle vacuum."},{"cited_title":"Vogel and M","cited_arxiv_id":null,"evidence_quote":"Defines the 0νββ transition operator and the CD-Bonn short-range correlation parametrization used in the matrix-element calculation."}],"review_version":1}