{"id":"3fb4a896-98bc-45d5-b41e-7772ff45e4d8","arxiv_id":"2411.18404","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"An iterative solver for the ATDHFB equations computes nuclear collective inertia, including time-odd mean fields, without explicit stability-matrix inversion, and matches dynamical cranking results.","lead":"This paper presents an iterative method for solving the adiabatic time-dependent Hartree-Fock-Bogoliubov equations, which compute the collective inertia of deformed atomic nuclei. The method avoids computing the full stability matrix and is tested against standard cranking calculations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fixed-point iteration in Eq. (6) lacks convergence proof; stopping on M can conceal large residuals, so the central ATDHF(B) solution claim is not yet established.","rationale":"I read the paper as claiming that the iterative solver provides a practical route to ATDHFB collective inertia without the stability matrix. The strongest independent evidence is the rotational benchmark for 20Ne and 126Ba against dynamical cranking; that is a real, nontrivial validation of the converged ATDHF results. However, that validation presupposes convergence. The reader's weakest assumption identifies exactly this: Section 2 defines the iteration and stops it on M with no proof or residual criterion. I agree this is the load-bearing point, because every number in the paper is produced by this iteration; if it converges to a non-solution or does not converge for untested soft or pairing regimes, the central claim fails regardless of benchmark agreement. I would not change the CONDITIONAL verdict: the concern is substantive but not yet shown to be realized, and the existing DC comparisons provide some evidence that for the tested rotational cases the iteration is finding the correct solution. The missing ATDHFB derivation in Ref. [5] and the 4% off-diagonal spread in Table 1 are secondary but reinforce the need for the residual-based convergence check.","tokens_in":5616,"tokens_out":10781,"duration_ms":104744,"concrete_test":"Record the residual R^{(n)}_{ph} = (ε_p − ε_h)ρ1^{(n)}_{ph} + Γ1^{(n)}_{ph} − i ˙q (∂ρ0/∂q)_{ph} after each iteration for the 20Ne and 74Ge cases, and launch the iteration from both Γ1 = 0 and Γ1 taken from the converged DC density. Accept convergence only when ‖R^{(n)}‖/‖b‖ < 10^{-6}; check that both initializations give the same M. If M stabilizes while ‖R‖ remains above threshold, the stopping criterion is invalid; if the residual converges for all tested deformation and pairing cases, the fixed-point concern is settled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the fixed-point iteration of Section 2 converges to the solution of the ATDHF(B) equation. Written as ρ1^{(n+1)} = D^{-1}(b − Lρ1^{(n)}) with D_{ph} = ε_p − ε_h and Lρ1 = Γ1(ρ1)_{ph}, Eq. (6) is a Richardson iteration for the linear system (D + L)ρ1 = b. A converged fixed point indeed solves the ATDHF equation, but the paper gives no spectral-radius bound, no contraction estimate, no residual check, and no discussion of uniqueness. The stopping criterion is the change in the scalar M (Section 2), not the residual of Eq. (6). Since M is an integrated function of ρ1, its stabilization can hide a large remaining error in ρ1; for soft or strongly paired systems the time-odd response L need not be small relative to D, so ‖D^{-1}L‖ may exceed 1 and the iteration may diverge or cycle. The ATDHFB equations themselves are deferred to Ref. [5], so the pairing case cannot be checked from this text. The agreement of converged ATDHF and DC results in Section 3.1 is evidence for the method when it converges, but it does not establish that the iteration reaches the ATDHF(B) solution in general.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes an iterative method for solving the adiabatic time-dependent Hartree-Fock(-Bogoliubov) equations, avoiding explicit construction of the stability matrix, and applies it to compute collective inertia for rotation and vibration. The time-odd density is updated via Eq. (6) in the HF single-particle basis, with the time-odd mean field recalculated at each step until the collective inertia M stops changing. Numerical tests are presented for the rotational moment of inertia of 20Ne and 126Ba compared with dynamical cranking, and for the quadrupole vibrational mass tensor of 74Ge, where the density-derivative term is evaluated by finite differences. The ATDHFB formalism itself is not given in this text: the paper states that the detailed description is deferred to a forthcoming publication [5]. The central claims are that the iterative method reproduces dynamical cranking results exactly in the full particle space and that it can be applied to superfluid nuclei.","tokens_in":5919,"tokens_out":4213,"duration_ms":38253,"significance":"If the method is correct and convergent, it offers a practical route to collective inertia for deformed and superfluid nuclei, which is important for fission and vibrational studies and would bypass the numerically prohibitive stability-matrix calculation. The ATDHF results in Section 3.1 and 3.2 agree with dynamical cranking to a high numerical precision, which is a meaningful implementation check: the comparison is not used to fit any parameter, and the two methods are distinct numerical routes within the same energy functional. The main limitation is that the paired ATDHFB case, which is the headline of the title and abstract, is not actually derived or benchmarked in this manuscript, and the convergence of the fixed-point iteration is not analyzed; the text itself acknowledges the former by deferring the equations to Ref. [5].","major_comments":[{"comment":"The manuscript does not establish that the fixed-point iteration converges to the solution of the ATDHF equation, nor that it converges at all for the cases shown. Written in the form ρ1^{(n+1)} = D^{-1}(b − Lρ1^{(n)}) with D_{ph} = ε_p − ε_h and (Lρ1)_{ph} = Γ1(ρ1)_{ph}, the scheme is a Richardson iteration for the linear system (D + L)ρ1 = b; convergence requires a bound on the iteration matrix. The stopping criterion, 'the desired precision for M' in Section 2, controls only the variation of a scalar integral of ρ1, not the residual of Eq. (6). The conclusion that the algorithm is 'rapidly converging' therefore needs support from either a contraction estimate, a residual check, or an iteration-count/residual analysis for representative deformations and pairing strengths.","section":"Section 2, Eq. (6)"},{"comment":"The only ATDHFB result, the quadrupole vibrational mass tensor of 74Ge, cannot validate the ATDHFB formalism as presented because the ATDHFB equations are not written in this manuscript, being explicitly deferred to Ref. [5]. Moreover, the comparison in Fig. 4 and Table 1 is made against the Inglis-Belyaev formula and against the numerical-differentiation spread, not against an independent solution of the paired ATDHFB equations. Therefore the abstract's claim that 'The ATDHFB equation is solved iteratively' is not substantiated within this text; the central novelty implied by the title remains unverified here.","section":"Section 3.3 and Abstract"},{"comment":"The 'exact correspondence between the ATDHF and DC moment of inertia in the full single-particle space' is demonstrated only for 20Ne with the SVT functional, and the companion case in Section 3.2 reports agreement only as a few numerical values. Because both calculations are performed in the same code (HFODD) with the same energy functional, the agreement is a self-consistency check of the implementation rather than an independent validation of the physical input; the text should state this caveat explicitly, and at least one more detailed case or a different functional would make the claim more robust.","section":"Section 3.1, Fig. 1"}],"minor_comments":[{"comment":"The notation ρ1,ph^{(n+1)} and Γ1,ph^{(n)} should specify the basis and the finite Hilbert-space truncation used in HFODD; in particular, degenerate or near-degenerate particle-hole energy denominators ε_p − ε_h require a comment on numerical stability.","section":"Section 2, Eq. (6)"},{"comment":"The statement that a 'rapidly converging iterative algorithm' was developed is not supported by any quantitative iteration counts, convergence rates, or residual norms in the manuscript.","section":"Conclusions"},{"comment":"The label 'ADB' in Fig. 4 appears to be a typo for 'ATDHFB' or 'ADB' should be defined; also, the figure caption does not explain the horizontal-axis units or the range shown.","section":"Section 3.3, Fig. 4"},{"comment":"Reference [10] lists 'Phys. Rev. C 109, L051301 (2014)'; volume 109 corresponds to 2024, so the year appears to be a typo.","section":"References"},{"comment":"The six 'different orientations' mentioned in the text are not defined in the figure caption or in the text; please specify the orientation angle or the rotation of the principal axes used.","section":"Section 3.2, Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings contribution and the companion paper [5] likely contains the missing ATDHFB derivation and convergence analysis. If so, the authors should either include a brief but self-contained statement of the paired equations and a convergence check, or explicitly scope the present text as an ATDHF-only validation. The comparison with dynamical cranking is a legitimate and useful implementation test, but the title and abstract promise more than the manuscript delivers in its current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: a sensible iterative method that avoids constructing the stability matrix, with clean ATDHF validation against dynamical cranking for rotations. The ATDHFB part of the story is mostly deferred to a companion paper, and the fixed-point iteration has no convergence analysis. Worth engaging with, but the proceedings format leaves the central paired claim under-supported.\n\nWhat is actually new: the specific combination of SVD-based adiabatic basis, self-consistent treatment of time-odd fields, and numerical differentiation for the vibrational mass tensor goes beyond earlier iterative work in Refs. [7] and the finite-amplitude implementations in Refs. [9,10]. The rotational tests are the strongest part. In particular, the 20Ne cutoff sensitivity and the 126Ba triaxial orientation invariance are clean numerical checks that the iteration, when it converges, lands on the ATDHF answer. That agreement is not circular in a harmful way; dynamical cranking is a different route through the same functional, so it is an internal consistency test rather than an external benchmark, but it is still meaningful.\n\nThe soft spots are real but proportionate. There is no contraction estimate, no spectral-radius bound, and no residual check for Eq. (6); stopping on the scalar M alone can hide large errors in rho1, especially in soft or strongly paired systems. The ATDHFB equations themselves do not appear in this text—they are deferred to Ref. [5]—so one cannot verify that the paired solver is actually solving what the authors claim. The only ATDHFB result, the 74Ge vibrational mass, has no independent benchmark and a 4% spread in the off-diagonal component. That is not fatal, but it tempers the conclusions' phrasing of 'perfect agreement' and 'rapidly converging.' Also, the conclusions say the rotational moments were obtained with the ATDHFB method, but the rotations in the paper are ATDHF without pairing; that looks like a small overstatement. No code or data are included, which is normal for a proceedings, but it limits independent checking.\n\nWho gets value from this: practitioners computing collective inertia for fission or mass tables, especially those who want to avoid the stability matrix. The paper is too short to carry the full weight of the method, but it does point to a concrete implementation in hfodd and gives a credible first validation.\n\nRecommendation: I would referee it. The desk should send it to review. The referee should ask for the ATDHFB equations or a clear pointer to the companion paper, and for a discussion of convergence and the stopping criterion. The rotational validation alone justifies the referee time.","headline":"A useful proceedings paper on avoiding the stability matrix for collective inertia, with strong ATDHF rotational tests but the ATDHFB claim deferred and convergence unproven.","tokens_in":6455,"tokens_out":2354,"would_cite":true,"duration_ms":22735,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V35","65H10"],"pacs":["21.60.Jz"],"model":"deepseek-v4-flash","headline":"The paper develops an iterative solver for the ATDHFB equations that computes nuclear collective inertia with time-odd mean fields included self-consistently, without ever constructing the two-body stability matrix.","keywords":["collective inertia","ATDHFB","time-odd mean fields","moment of inertia","Skyrme density functional theory","dynamical cranking","vibrational mass tensor","fixed-point iteration"],"falsifier":"Compute the collective inertia for a fixed deformation and pairing strength by the iterative method and by a direct solution of the ATDHFB equations (or by dynamical cranking at sufficiently small frequency) using the same interaction and basis; a disagreement beyond the numerical tolerance, or a failure of the iteration to converge, would show that the method is not delivering the exact ATDHFB inertia.","tokens_in":5440,"feed_emoji":"⚛️","tokens_out":6937,"duration_ms":53963,"temperature":0.7,"pith_summary":"The paper develops an iterative way to solve the adiabatic time-dependent Hartree-Fock-Bogoliubov (ATDHFB) equations, which are the standard microscopic route to the collective inertia of a nucleus. The iteration avoids constructing the two-body stability matrix, the step that normally makes ATDHFB calculations prohibitive for deformed or superfluid nuclei. The method includes the time-odd mean fields self-consistently, and the paper verifies it by showing that the resulting rotational inertia matches dynamical cranking values exactly when the full single-particle space is kept. It also computes the vibrational inertia mass tensor for 74Ge by numerical differentiation of the density, reaching about 1% accuracy on the diagonal components. If the method holds, collective inertia for rotation, vibration, and fission becomes practical for a much wider range of nuclei.","feed_headline":"New iteration yields exact nuclear inertia without stability matrix","feed_subtitle":"Self-consistent time-odd mean fields give rotation and vibration inertia that match dynamical cranking.","key_machinery":"The fixed-point iteration (Eq. 6) expresses the time-odd density correction $\\rho_1$ in the Hartree-Fock basis as $(\\epsilon_p-\\epsilon_h)^{-1} (i\\dot{q}\\,\\partial\\rho_0/\\partial q - \\Gamma_1)$; because the time-odd mean field $\\Gamma_1$ is a functional of $\\rho_1$, this becomes a self-consistency loop. The loop uses the singular value decomposition of $\\rho_1$ to construct the adiabatic basis of paired 'occupied' eigenstates, which lets the code evaluate time-odd densities and currents with the same machinery as the static time-even densities. The inertia is extracted from the collective kinetic energy $K = \\tfrac12 \\mathrm{Tr}(\\dot{\\rho}_0 \\chi) = \\tfrac12 M\\dot{q}^2$.","core_discovery":"The central claim is that the ATDHF/ATDHFB equations can be solved by a fixed-point iteration on the particle-hole matrix elements of the time-odd density, given explicitly by Eq. (6), without ever building the stability matrix. Starting from zero time-odd mean field (the Inglis-Belyaev approximation), each step updates the density correction, reconstructs the adiabatic basis via the singular value decomposition of the time-odd density, recomputes the time-odd mean fields, and evaluates the inertia; the iteration stops when the inertia converges. The paper shows that in the full single-particle space the ATDHF moment of inertia of 20Ne coincides exactly with dynamical cranking, that the triaxial 126Ba moments of inertia along all three axes agree with dynamical cranking to the displayed precision and are independent of the nucleus orientation, and that the quadrupole vibrational inertia of 74Ge stabilizes to within 1% as the finite-difference step for the density derivative is reduced.","pith_inferences":["If the iteration converges reliably for all deformations along a fission path, the method could replace stability-matrix-based inertia in large-scale fission calculations, where the present bottleneck is exactly the matrix inversion.","The observed sensitivity to the single-particle cutoff in 20Ne suggests that practical truncations in heavier nuclei could bias the moment of inertia, so the full-space limit is the safe operating point.","A natural testable extension is to apply this iteration to odd-mass or well-deformed rare-earth nuclei and compare with measured ground-state band moments of inertia.","The fixed-point iteration has no convergence guarantee in the paper; using a mixing or Anderson acceleration may be needed near level crossings or for strong pairing."],"forward_implications":["Rotational moments of inertia for arbitrarily deformed nuclei, including triaxial shapes, can now be computed without the numerically prohibitive two-body stability matrix.","The self-consistent time-odd mean fields raise the moment of inertia above the Inglis-Belyaev value by the expected factor of roughly 1.2 to 1.4, so the method removes the main systematic underestimate of the cranking approximation.","The same iterative solver applies to vibrational inertia, with the density-derivative term handled by numerical differentiation at about 1% accuracy for diagonal components.","Because the method is implemented in a Cartesian deformed harmonic-oscillator basis, it is directly applicable to fission paths where deformations are arbitrary and pairing is active."],"supporting_citations":[{"why":"Supplies the ATDHF decomposition of the density and the equation of motion that the iteration solves.","marker":"[6]"},{"why":"Shows the earlier iterative solution of ATDHFB equations that the present method extends.","marker":"[7]"},{"why":"Defines the Inglis cranking moment of inertia, the starting point of the iteration.","marker":"[2]"},{"why":"Extends the cranking formula to superfluid systems, giving the baseline Inglis-Belyaev inertia for the ATDHFB case.","marker":"[3]"},{"why":"Documents the 1.2 to 1.4 factor by which neglecting time-odd mean fields underestimates collective masses.","marker":"[4]"},{"why":"Provides the SVT Skyrme interaction used for the rotational inertia benchmarks.","marker":"[13]"},{"why":"Provides the SkM* interaction and pairing used for the vibrational inertia calculation.","marker":"[15]"},{"why":"Describes the solver in which the iterative method is implemented.","marker":"[11]"}],"fun_headline_variants":["Iterative method skips stability matrix for exact nuclear inertia","Self-consistent iteration nails nuclear inertia, no stability matrix","Fixed-point iteration yields collective inertia, matches cranking","Iterative ATDHFB solver: exact inertia without stability matrix","New iteration matches cranking, no stability matrix needed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fixed-point iteration is assumed to converge to the unique solution of the ATDHF/ATDHFB equations; the paper gives no proof, convergence rate, or test for multiple fixed points.","fun_headline_variants_meta":{"raw":{"variants":["Iterative method skips stability matrix for exact nuclear inertia","Self-consistent iteration nails nuclear inertia, no stability matrix","Fixed-point iteration yields collective inertia, matches cranking","Iterative ATDHFB solver: exact inertia without stability matrix","New iteration matches cranking, no stability matrix needed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000809,"raw_usage":{"total_tokens":3490,"prompt_tokens":823,"completion_tokens":2667,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":439,"completion_tokens_details":{"reasoning_tokens":2587}},"tokens_in":439,"tokens_out":2667,"duration_ms":18677,"temperature":1.0,"reasoning_tokens":2587,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:15:15.481467+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the collective inertia for a fixed deformation and pairing strength by the iterative method and by a direct solution of the ATDHFB equations (or by dynamical cranking at sufficiently small frequency) using the same interaction and basis; a disagreement beyond the numerical tolerance, or a failure of the iteration to converge, would show that the method is not delivering the exact ATDHFB inertia.","supporting_citations":[{"cited_title":"Dobaczewski and J","cited_arxiv_id":null,"evidence_quote":"Shows the earlier iterative solution of ATDHFB equations that the present method extends."},{"cited_title":"Baranger and M","cited_arxiv_id":null,"evidence_quote":"Supplies the ATDHF decomposition of the density and the equation of motion that the iteration solves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Inglis cranking moment of inertia, the starting point of the iteration."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extends the cranking formula to superfluid systems, giving the baseline Inglis-Belyaev inertia for the ATDHFB case."},{"cited_title":"Libert, M","cited_arxiv_id":null,"evidence_quote":"Documents the 1.2 to 1.4 factor by which neglecting time-odd mean fields underestimates collective masses."},{"cited_title":"Satu la et al., Phys","cited_arxiv_id":null,"evidence_quote":"Provides the SVT Skyrme interaction used for the rotational inertia benchmarks."},{"cited_title":"Bartel et al., Nucl","cited_arxiv_id":null,"evidence_quote":"Provides the SkM* interaction and pairing used for the vibrational inertia calculation."},{"cited_title":"Dobaczewski et al., J","cited_arxiv_id":null,"evidence_quote":"Describes the solver in which the iterative method is implemented."}],"review_version":1}