{"id":"bfabbb2f-16c0-4453-8fc1-adae0c0cd4a9","arxiv_id":"2603.24399","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new Python package automatically derives symbolic commutators of normal-ordered many-body operators and generates multi-reference IMSRG(3) flow equations.","lead":"This paper presents qcombo, a Python package that automatically expands commutators of normal-ordered many-body operators using the generalized Wick theorem, and uses it to write down MR-IMSRG(3) flow equations. It matters because long many-body commutator derivations are error-prone, and automated export to LaTeX/amc could accelerate nuclear and quantum-chemistry implementations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Contraction enumeration in Wick.generalizedWick is the load-bearing unvalidated core; the two SR/MR benchmarks leave the new λ^(2)/λ^(3) and 3-body-operator contractions in 23X/33X completely untested, so a missed contraction or sign error would silently corrupt all Appendix B equations.","rationale":"The reader's weakest assumption is correct and precisely identifies the load-bearing point: the correctness of every generated expression depends on the completeness and sign conventions of the contraction enumeration in Wick.generalizedWick, and the paper offers no formal proof or independent test of that enumeration. I considered other issues—the abstract's 'complete set' overstatement versus Sec. 4's λ^(3) truncation, and the unusual appearance of dW/ds λ^(3) in Eq. (40)—but these are secondary to the enumeration question. The density-matrix truncation is a scope/claim mismatch, not a correctness failure of the package; the dW/ds term is suspicious but could be a legitimate shorthand for a commutator contribution. The fundamental issue is that the new 23X and 33X content, which is the claimed contribution, is exactly what the two benchmark limits do not test. A randomized numerical commutator check would settle whether the enumeration is complete and correctly signed, and the same test can validate the dE/ds expression. I therefore recommend keeping the CONDITIONAL verdict: the central claim is plausible and useful, but acceptance should require the numerical cross-check to pass.","tokens_in":23647,"tokens_out":24730,"duration_ms":248876,"concrete_test":"Run a randomized numerical cross-check: in a small single-particle basis (e.g., N=6) with a correlated reference whose λ^(1), λ^(2), λ^(3) are all nonzero (e.g., a BCS or configuration-interaction state), generate random anti-symmetrized matrix elements for normal-ordered 1-, 2-, and 3-body operators. For every pair (m,n) ∈ {1,2,3}², construct the explicit operator matrices in Fock space using Eq. (7), compute the commutator matrix numerically, and compare it with the matrix reconstructed from qcombo's symbolic output for the same inputs. Repeat over many draws; any mismatch in any contraction rank isolates a missing or mis-signed contraction. As a secondary check, evaluate the right-hand side of Eq. (40) and compare with the numerically computed ⟨Φ|[η,H]|Φ⟩ to verify the dW/ds λ^(3) term.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that qcombo-automated generalized Wick contractions produce correct commutators, and that the Appendix B equations are the complete NO3B MR-IMSRG(3) set. The load-bearing condition is the completeness and sign convention of Wick.generalizedWick (Sec. 3.2.2), which enumerates all contractions and signs for products of normal-ordered m- and n-body operators. Nothing in the paper proves this enumeration: there is no formal proof, no machine-checked verification, and the only benchmarks are reductions to SR-IMSRG(3) and MR-IMSRG(2). Those limits do not exercise the new content. SR-IMSRG(3) sets λ^(k≥2)=0, so all generalized contractions involving two- and three-body irreducible densities vanish; MR-IMSRG(2) has no three-body normal-ordered operators, so the 13X, 23X, and 33X commutators — precisely the terms that make MR-IMSRG(3) new — are not tested. Appendix A states the displayed expressions 'differ slightly from the direct output of the program' after manual reformulation, adding an unverified transcription layer. A missed contraction (e.g., an overlooked pairing among the (m+n)! index permutations of two 3-body operators) or an incorrect sign for a ξ contraction in a λ^(3) term would propagate silently into every Appendix B flow equation and into the dE/ds expression, Eq. (40), where the dW/ds λ^(3) term is likewise absent in both benchmark limits. This is not a disagreement with consensus; it is an internal validation gap. The package may well be correct, but the central correctness claim is not yet secured by the presented evidence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents qcombo, a Python package for symbolic evaluation of commutators of normal-ordered many-body operators using the generalized Wick theorem. The package pipeline consists of input, commutator generation, regularization, simplification, and output in LaTeX and amc formats, with an easyCombo one-call interface. As a demonstration, the authors use qcombo to generate a set of MR-IMSRG(3) flow equations for operators truncated at the normal-ordered three-body level, collected in Appendix B. The paper claims that this is the complete NO3B MR-IMSRG(3) set, benchmarked against SR-IMSRG(3) and MR-IMSRG(2) limits.","tokens_in":24117,"tokens_out":11393,"duration_ms":115413,"significance":"If the contraction enumeration is correct, qcombo is a useful and timely tool: automated commutator evaluation is a practical bottleneck in many-body derivations, and the package is open-source, installable from PyPI, and comes with a reproducible Jupyter example. The interface with the amc package for J-scheme output is a genuine strength, as is the explicit reduction of the generated equations to known SR-IMSRG(3) and MR-IMSRG(2) limits. However, the paper’s central completeness claim is weakened by an additional truncation stated in Sec. 4, and the new content — contractions involving irreducible two- and three-body densities, and 13X, 23X, 33X commutators — is not independently validated. These gaps are fixable but currently load-bearing.","major_comments":[{"comment":"The abstract claims a 'complete set' of NO3B MR-IMSRG(3) flow equations, but Sec. 4 states that only the zero-body flow equation retains contributions up to the three-body irreducible density matrix, while the one-, two-, and three-body flow equations include only λ^(1) and λ^(2) terms. Appendix A shows that the 33X one- and two-body commutator pieces contain λ^(3) terms (e.g., λ^{cde}_{fgh} in Eq. (A.20) and λ^{efg}_{hij} in Eq. (A.21)), and these are absent from the corresponding flow equations (B.2), (B.10)-(B.11), and (B.16)-(B.17). Appendix B is therefore not the complete NO3B set claimed in the abstract; it is an additionally truncated set. The authors should either revise the claim or provide the full equations and justify that the omitted λ^(3) contributions are not needed for their stated MR-IMSRG(3) target.","section":"Sec. 4 / Abstract"},{"comment":"Correctness of every generated expression, including all of Appendix B, rests on the completeness and sign conventions of the contraction enumeration in Wick.generalizedWick. No formal proof, machine-checked verification, or independent numerical test of this enumeration is provided. The two benchmarks do not exercise the new content: SR-IMSRG(3) sets λ^(k≥2)=0, so contractions involving irreducible two- and three-body densities are untested; MR-IMSRG(2) has no three-body operators, so the 13X, 23X, and 33X commutators are untested. A missed contraction or sign error in a ξ contraction in a λ^(3) term would silently corrupt all Appendix B equations and Eq. (40). A concrete remedy is to compare qcombo output against brute-force numerical commutators in a small single-particle space (e.g., random matrices for [A^(3),B^(3)] in a 5-6 orbital basis), or to compare with an independent MR-IMSRG","section":"Sec. 3.2.2 / Sec. 4 benchmarks"},{"comment":"The statement that the displayed expressions 'differ slightly from the direct output of the program' after manual reformulation introduces an unverified transcription layer. Since the paper’s central claim is automated generation, the reader cannot check that Appendix B is exactly what easyCombo produces. Please provide a script or notebook that reproduces each Appendix B equation from the package, or specify the reformulation rules and verify them algebraically. This is especially important given the absence of independent validation noted above.","section":"Appendix A, introductory paragraph"}],"minor_comments":[{"comment":"Typos: 'atttibute' should be 'attribute'; 'tow-body' should be 'two-body'; 'ca lc ul at e' in a code block should be 'calculate'.","section":"Sec. 3.2.2 / 3.2.3"},{"comment":"The LaTeX output contains 'n^{}_c'; it should be 'n_c'. Also, the text explaining the 1/4 prefactor cancellation is terse and would benefit from a short derivation linking Eq. (36) to the final one-body equation.","section":"Sec. 3.2.5"},{"comment":"The use of dΓ/ds and dW/ds terms inside the zero-body equation is surprising at first reading. A sentence explaining that these are chain-rule reformulations of the direct zero-body commutator contributions (rather than additional dynamical terms) would improve clarity.","section":"Sec. 4, Eq. (40)"},{"comment":"The phrase 'Up to a sign factor (−1) that unifies the sign conventions' is ambiguous: the sign factor is not written explicitly. Please clarify whether ξ^i_j = λ^i_j − δ^i_j or ξ^i_j = δ^i_j − λ^i_j, and the relation to the later natural-orbital expression.","section":"Sec. 2, Eq. (17)"},{"comment":"The benchmark against MR-IMSRG(2) in Ref. [13] would be more persuasive if the conventions (natural-orbital basis, irreducible densities, operator normal ordering) used there are explicitly matched to those of this paper, since a mismatch in conventions could mask discrepancies.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful software/derivation contribution, and the open-source, reproducible workflow is a clear strength. My main concern is not that the derivation is circular or internally inconsistent, but that the new, load-bearing parts — generalized contraction enumeration for higher-rank densities and the full 33X commutator — are not independently tested. In addition, the 'complete set' claim in the abstract is stronger than what Sec. 4 actually delivers, since λ^(3) contributions are dropped from the one-, two-, and three-body flow equations. Both issues are fixable with additional verification and a textual revision, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper introduces qcombo, a Sympy-based tool that computes commutators of normal-ordered many-body operators via the generalized Wick theorem, and uses it to generate the NO3B MR-IMSRG(3) flow equations in Appendix B. The package is real, installable, and clearly described; the appendix equations are a substantial resource for the IMSRG community. That is the main value.\n\nWhat the paper does well: it automates a genuinely tedious step, handles the natural-orbital simplification, exports to the amc format, and the worked one-body/two-body example is transparent. The reduction to SR-IMSRG(3) and MR-IMSRG(2) limits is a useful sanity check.\n\nThe soft spot is the one I would flag before building on anything in Appendix B. The correctness of every generated expression rests on Wick.generalizedWick's enumeration of all contractions and signs. The paper offers no proof, no machine-checked verification, and no independent test of that enumeration. Worse, the two benchmark limits they compare against exactly remove the terms that are new here: SR-IMSRG(3) kills all two- and three-body irreducible densities, and MR-IMSRG(2) has no three-body operators. So the 13X, 23X, and 33X commutators — where the combinatorial complexity actually lives — are effectively unexercised. A missed pairing or sign error would silently flow into every equation in Appendix B, including the energy derivative. The statement in Sec. 4 that the printed expressions \"differ slightly from the direct output of the program\" adds a small transcription risk on top of that.\n\nI also think the abstract's \"complete set\" overstates things. Section 4 explicitly keeps only one- and two-body irreducible densities in the one-, two-, and three-body flow equations; only the zero-body equation keeps the three-body density. That is a defensible truncation, but the abstract should say so.\n\nNone of this makes the paper unserious. It is internally consistent, the authors clearly know the formalism, and the package is a real contribution. But the central correctness claim is not yet secured by the presented evidence. If I were the editor, I would send it to peer review and ask for independent numerical checks of the 33X terms — e.g., comparing a few generated commutators against a brute-force antisymmetrized evaluation in a small model space, or a machine-checked proof of the contraction enumeration. Until then, I would be hesitant to cite or build on the full expression set.","headline":"A genuinely useful automated Wick-contraction package and a large MR-IMSRG(3) expression set, but the central correctness claim rests on an unvalidated contraction enumerator that the benchmark limits do not actually exercise.","tokens_in":24564,"tokens_out":2513,"would_cite":false,"duration_ms":25146,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"qcombo automates the symbolic derivation of commutators between normal-ordered many-body operators, and the paper uses it to generate the complete set of MR-IMSRG(3) flow equations, replacing error-prone hand calculations.","keywords":["normal ordering","generalized Wick theorem","commutator","IMSRG","symbolic computation","many-body operators","flow equations","multi-reference"],"falsifier":"Take a commutator short enough to verify independently—for example [1B,2B] contracting to 1B—and compare qcombo's symbolic output with the known analytic result; then test a rank-3 term by evaluating both sides of the commutator identity with random numerical density matrices and fractional occupation numbers. Any discrepancy in terms or signs would indicate a defect in the contraction enumeration.","tokens_in":23578,"feed_emoji":"⚛️","tokens_out":7256,"duration_ms":62395,"temperature":0.7,"pith_summary":"The paper introduces qcombo, a Python package that symbolically evaluates commutators between normal-ordered many-body operators using the generalized Wick theorem. It sets out to show that such commutators—traditionally derived by hand, tediously and with frequent errors—can be generated automatically, simplified, and exported in symbolic form. The demonstration: qcombo generates the complete set of multi-reference IMSRG flow equations with operators truncated at the normal-ordered three-body level (MR-IMSRG(3)). If correct, this removes a major source of human error in many-body derivations and makes higher-rank truncations practical for nuclear structure and quantum chemistry.","feed_headline":"Automated tool derives full MR-IMSRG(3) flow equations","feed_subtitle":"Hand-deriving three-body commutators invites errors; qcombo generates, simplifies, and exports them symbolically.","key_machinery":"The load-bearing mechanism is the contraction enumeration in the generalized-Wick routine: every pairing between upper and lower indices of the two normal-ordered operators is generated, each with a sign fixed by index ordering, and the commutator is obtained by subtracting the two operator orderings. The package then filters terms by the body rank of the resulting normal-ordered operator and simplifies using the antisymmetry of two- and higher-body matrix elements, renaming of dummy indices, and the natural-orbital diagonal form λ^i_j = n_i δ^i_j (with ξ = λ − 1). The output step restores the combinatorial prefactor (k!)²/(m!)²(n!)² needed for the matrix-element equation.","core_discovery":"The central claim is that qcombo produces correct symbolic commutators for general normal-ordered many-body operators, replacing hand derivations even at three-body rank. The package enumerates all contractions prescribed by the generalized Wick theorem, computes the commutator as product minus swapped product, simplifies using antisymmetry of matrix elements and the natural-orbital diagonalization of the one-body density matrix, and outputs LaTeX or input ready for an angular-momentum-coupled code. The representative result is the entire MR-IMSRG(3) flow-equation set—zero-, one-, two-, and three-body flow equations—with each contribution labeled by operator rank and density-matrix rank. The","pith_inferences":["A formal proof or a battery of random numerical tests of the contraction enumeration would make the package trustworthy beyond its own benchmark cases; the present paper does not provide one, so confidence currently rests on reduction to prior special cases.","The same strategy should extend to number-breaking operators and possibly to bosonic or non-scalar operators, though the paper explicitly limits its current scope to number-conserving operators and does not claim these extensions.","Publishing the full MR-IMSRG(3) equations in appendices effectively turns the paper into a reference dataset that other groups can use to validate their own symbolic or numerical implementations of IMSRG(3)."],"forward_implications":["The full MR-IMSRG(3) flow equations become available as a starting point for numerical implementations, enabling systematic study of three-body truncation effects in open-shell nuclei.","The same automation extends to other many-body methods that require normal-ordered commutators, such as multi-reference coupled-cluster with full triplets.","Users can regenerate the flow equations quickly under different generator choices or truncation schemes, lowering the barrier to exploring variants of IMSRG.","Because the resulting expressions are symbolic and machine-readable, they can be inspected, checked, and exported directly to the angular-momentum-coupled format needed for practical nuclear-structure calculations."],"fun_headline_variants":["Qcombo automates full MR-IMSRG(3) flow equation derivation","Python package derives three-body commutators symbolically","Automated symbolic commutators for nuclear many-body theory","Qcombo replaces hand derivations in IMSRG flow equations","From Wick theorem to flow equations automatically"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The correctness of every generated expression rests on the completeness and sign conventions of the contraction enumeration in the generalized-Wick routine; a missed contraction or wrong sign in a high-rank term would silently propagate into all published flow equations, and the paper gives no formal proof or independent numerical test of that enumeration.","fun_headline_variants_meta":{"raw":{"variants":["Qcombo automates full MR-IMSRG(3) flow equation derivation","Python package derives three-body commutators symbolically","Automated symbolic commutators for nuclear many-body theory","Qcombo replaces hand derivations in IMSRG flow equations","From Wick theorem to flow equations automatically"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000278,"raw_usage":{"total_tokens":1452,"prompt_tokens":670,"completion_tokens":782,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":414,"completion_tokens_details":{"reasoning_tokens":700}},"tokens_in":414,"tokens_out":782,"duration_ms":7258,"temperature":1.0,"reasoning_tokens":700,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T17:28:19.244758+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a commutator short enough to verify independently—for example [1B,2B] contracting to 1B—and compare qcombo's symbolic output with the known analytic result; then test a rank-3 term by evaluating both sides of the commutator identity with random numerical density matrices and fractional occupation numbers. Any discrepancy in terms or signs would indicate a defect in the contraction enumeration.","supporting_citations":[],"review_version":1}