Pith. sign in

REVIEW 3 major objections 5 minor 2 cited by

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.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 17:28 UTC pith:DDTHJWRA

load-bearing objection 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. the 3 major comments →

arxiv 2603.24399 v2 pith:DDTHJWRA submitted 2026-03-25 nucl-th cond-mat.str-el

Qcombo: A Python Package for Automated Commutator Calculations of Quantum Many-Body Operators

L. H. Chen , Y. Li , H. Hergert , J. M. Yao This is my paper
classification nucl-th cond-mat.str-el
keywords normal orderinggeneralized Wick theoremcommutatorIMSRGsymbolic computationmany-body operatorsflow equationsmulti-reference
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

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.

Core claim

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

What carries the argument

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.

Load-bearing premise

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.

What would settle it

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.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • 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.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • 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).

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

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.

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 (3)
  1. [Sec. 4 / Abstract] 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.
  2. [Sec. 3.2.2 / Sec. 4 benchmarks] 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
  3. [Appendix A, introductory paragraph] 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.
minor comments (5)
  1. [Sec. 3.2.2 / 3.2.3] Typos: 'atttibute' should be 'attribute'; 'tow-body' should be 'two-body'; 'ca lc ul at e' in a code block should be 'calculate'.
  2. [Sec. 3.2.5] 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.
  3. [Sec. 4, Eq. (40)] 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.
  4. [Sec. 2, Eq. (17)] 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.
  5. [References] 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.

Circularity Check

0 steps flagged

No significant circularity: automated commutator expansions are direct symbolic applications of the generalized Wick theorem, and benchmark/self-citations are external prior results used for comparison, not inputs to the derivation.

full rationale

The paper's derivation chain is a symbolic expansion: normal-ordering definitions (Eqs. 6-15), generalized Wick contraction rules (Eq. 21) whose algebraic proof is cited to the external Ref. [22], the enumeration procedure in Wick.generalizedWick (Sec. 3.2.2), and the final NO3B commutators in Appendices A-B. There are no fitted constants, no parameter-forcing of an output, no prediction-vs-fit structure, and no uniqueness theorem imported from the authors' prior work. The benchmarks are stated as 'The above flow equations are benchmarked against the SR-IMSRG(3) results [8, 10]... Furthermore, when all three-body operators are neglected, the equations reduce to the MR-IMSRG(2) expressions reported in Ref. [13].' Although Refs. [8,13] include co-author Hergert, they are previously published external derivations used as checks, not as inputs to qcombo's contraction engine. The Appendix A note that displayed expressions 'differ slightly from the direct output of the program, although they are algebraically equivalent' is a manual transcription layer and a verification risk, but not a definitional reduction. The absence of tests for the genuinely new lambda^(2)/lambda^(3) and 33X contractions is an internal validation/completeness gap, not a circularity: a missed contraction would be an implementation bug, not an input-output equivalence. The central claim remains that qcombo applies an externally proved theorem, and the paper is self-contained relative to that theorem.

Axiom & Free-Parameter Ledger

0 free parameters · 6 axioms · 0 invented entities

The paper introduces no new physical entity and no fitted free parameters. Occupation numbers n_i and irreducible density matrices lambda^(k) appear as symbolic variables. The central derivation relies on the generalized Wick theorem, antisymmetry conventions, natural-orbital diagonalization, and the NO3B truncation, as well as on the correctness of the unverified contraction enumeration in the code.

axioms (6)
  • standard math Generalized Wick theorem contraction rules (Eq. 21) and recursive normal ordering (Eq. 7)
    Foundation of the package; proof cited to Ref [22], not re-derived.
  • standard math Fermionic antisymmetry of matrix elements (Eq. 28) and permutation operators
    Used in simplification and canonicalization (Sec. 3.2.4); standard for antisymmetrized matrix elements.
  • domain assumption One-body density matrix diagonalizable in natural orbital basis with 0<=n_i<=1 (Eq. 18)
    Used throughout simplification; valid for Hermitian density matrices but restricts the reference states treated.
  • domain assumption NO3B truncation and dropping lambda^(3) from f, Gamma, W flow equations (Sec. 4)
    The paper says only zero-body term retains 3-body irreducible densities; this is a modeling approximation, not a theorem.
  • domain assumption Number conservation of operators (Sec. 5 states limitation)
    Package only handles number-conserving operators; central claim limited to this class.
  • ad hoc to paper Contraction enumeration and sign logic in Wick.generalizedWick are complete and correct
    No formal proof or independent test; bugs would invalidate all generated equations.

pith-pipeline@v1.3.0-alltime-deepseek · 23427 in / 12579 out tokens · 132928 ms · 2026-08-02T17:28:19.244758+00:00 · methodology

0 comments
read the original abstract

qcombo is a Python package for the symbolic evaluation of commutators between general quantum many-body operators expressed in normal-ordered form using the generalized Wick theorem. The package provides an automated and systematic framework for generating the corresponding algebraic expressions, significantly reducing the risk of human error in lengthy and complex analytical derivations. It is designed to assist the development and implementation of modern many-body methods in nuclear physics, quantum chemistry, and related fields. The functionality and workflow of the package are demonstrated through an application to the in-medium similarity renormalization group (IMSRG) method, which has been widely used for nuclear ab initio calculations. As a representative example, qcombo is employed to automatically generate the complete set of multi-reference IMSRG flow equations with operators truncated at the normal-ordered three-body level.

Figures

Figures reproduced from arXiv: 2603.24399 by H. Hergert, J. M. Yao, L. H. Chen, Y. Li.

Figure 1
Figure 1. Figure 1: Flowchart of the main routine easyCombo in the qcombo package. The workflow includes the construction of commutators, application of contraction rules, regularization of intermediate expressions, simplification of algebraic terms, and generation of outputs in both LATEX and amc formats. 3.2.2. Commutator After defining the operators, their commutator is evaluated using the function Wick.commutate, which co… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. High-order perturbative calculations of nuclear ground states: Automated evaluation of many-body diagrams

    nucl-th 2026-06 unverdicted novelty 7.0

    Automated MBPT up to fifth order shows convergence trends in ground-state energies of closed-shell nuclei and decomposes fourth-order terms while comparing to IMSRG.

  2. High-order perturbative calculations of nuclear ground states: Automated evaluation of many-body diagrams

    nucl-th 2026-06 conditional novelty 6.0

    Automated MBPT through fifth order yields converging ground-state energies for closed-shell nuclei up to 78Ni and exposes missing triples/quadruples in IMSRG(2).

Reference graph

Works this paper leans on

23 extracted references · 3 canonical work pages · cited by 1 Pith paper

  1. [1]

    Coester, H

    F. Coester, H. Kümmel, Short-range correlations in nu- clear wave functions, Nucl. Phys. 17 (1960) 477–485. doi:10.1016/0029-5582(60)90140-1

    doi:10.1016/0029-5582(60)90140-1 1960

  2. [2]

    R. J. Bartlett, M. Musial, Coupled-cluster theory in quan- tum chemistry, Rev. Mod. Phys. 79 (2007) 291–352.doi: 10.1103/RevModPhys.79.291

    doi:10.1103/revmodphys.79.291 2007

  3. [3]

    Hagen, T

    G. Hagen, T. Papenbrock, M. Hjorth-Jensen, D. J. Dean, Coupled-cluster computations of atomic nuclei, Rept. Prog. Phys. 77 (9) (2014) 096302.arXiv:1312.7872, doi:10.1088/0034-4885/77/9/096302

    Pith/arXiv arXiv 2014

  4. [4]

    S. D. Głazek, K. G. Wilson, Perturbative renormalization group for Hamiltonians, Phys. Rev. D 49 (8) (1994) 4214– 4218.doi:10.1103/PhysRevD.49.4214

    work page doi:10.1103/physrevd.49.4214 1994

  5. [5]

    Wegner, Flow-equations for Hamiltonians, Annalen der Physik 506 (2) (1994) 77–91.doi:10.1002/andp

    F. Wegner, Flow-equations for Hamiltonians, Annalen der Physik 506 (2) (1994) 77–91.doi:10.1002/andp. 19945060203

    doi:10.1002/andp 1994

  6. [6]

    Tsukiyama, S

    K. Tsukiyama, S. K. Bogner, A. Schwenk, In-Medium Similarity Renormalization Group for Nuclei, Phys. Rev. Lett. 106 (2011) 222502.arXiv:1006.3639,doi:10. 1103/PhysRevLett.106.222502

    Pith/arXiv arXiv 2011

  7. [7]

    F. A. Evangelista, A driven similarity renormalization group approach to quantum many-body problems, The Journal of Chemical Physics 141 (5) (2014) 054109.arXiv:https://pubs.aip.org/aip/jcp/ article-pdf/doi/10.1063/1.4890660/13374939/ 054109_1_online.pdf,doi:10.1063/1.4890660. URLhttps://doi.org/10.1063/1.4890660

    work page doi:10.1063/1.4890660/13374939/ 2014

  8. [8]

    Hergert, S

    H. Hergert, S. K. Bogner, T. D. Morris, A. Schwenk, K. Tsukiyama, The In-Medium Similarity Renormaliza- tion Group: A Novel Ab Initio Method for Nuclei, Phys. Rept. 621 (2016) 165–222.arXiv:1512.06956,doi: 10.1016/j.physrep.2015.12.007

    Pith/arXiv arXiv 2016

  9. [9]

    S. R. Stroberg, S. K. Bogner, H. Hergert, J. D. Holt, Nonempirical Interactions for the Nuclear Shell Model: An Update, Ann. Rev. Nucl. Part. Sci. 69 (2019) 307–362.arXiv:1902.06154,doi:10.1146/ annurev-nucl-101917-021120

    Pith/arXiv arXiv 2019

  10. [10]

    Heinz, A

    M. Heinz, A. Tichai, J. Hoppe, K. Hebeler, A. Schwenk, In-medium similarity renormaliza- tion group with three-body operators, Phys. Rev. C 103 (4) (2021) 044318.arXiv:2102.11172, doi:10.1103/PhysRevC.103.044318

    Pith/arXiv arXiv 2021

  11. [11]

    S. R. Stroberg, T. D. Morris, B. C. He, In-medium similarity renormalization group with flowing 3-body operators, and approximations thereof, Phys. Rev. C 110 (2024) 044316.doi:10.1103/PhysRevC.110.044316. URLhttps://link.aps.org/doi/10.1103/ PhysRevC.110.044316

    work page doi:10.1103/physrevc.110.044316 2024

  12. [12]

    B. C. He, S. R. Stroberg, Factorized approximation to the in-medium similarity renormalization group IMSRG(3), Phys. Rev. C 110 (4) (2024) 044317.arXiv:2405. 19594,doi:10.1103/PhysRevC.110.044317

    doi:10.1103/physrevc.110.044317 2024

  13. [13]

    Hergert, S

    H. Hergert, S. Binder, A. Calci, J. Langhammer, R. Roth, Ab Initio Calculations of Even Oxygen Isotopes with Chiral Two-Plus-Three-Nucleon Interactions, Phys. Rev. Lett. 110 (24) (2013) 242501.arXiv:1302.7294,doi: 10.1103/PhysRevLett.110.242501

    Pith/arXiv arXiv 2013

  14. [14]

    Gebrerufael, K

    E. Gebrerufael, K. V obig, H. Hergert, R. Roth, Ab Initio Description of Open-Shell Nuclei: Merg- ing No-Core Shell Model and In-Medium Similar- ity Renormalization Group, Phys. Rev. Lett. 118 (15) (2017) 152503.arXiv:1610.05254,doi:10.1103/ PhysRevLett.118.152503. 15

    Pith/arXiv arXiv 2017

  15. [15]

    F. A. Evangelista, Perspective: Multireference coupled cluster theories of dynamical electron correlation, The Journal of Chemical Physics 149 (3) (2018)

    2018

  16. [16]

    J. M. Yao, B. Bally, J. Engel, R. Wirth, T. R. Ro- dríguez, H. Hergert, Ab initio treatment of collective correlations and the neutrinoless double beta de- cay of 48Ca, Phys. Rev. Lett. 124 (2020) 232501. doi:10.1103/PhysRevLett.124.232501. URLhttps://link.aps.org/doi/10.1103/ PhysRevLett.124.232501

    doi:10.1103/physrevlett.124.232501 2020

  17. [17]

    E. F. Zhou, C. R. Ding, J. M. Yao, B. Bally, H. Herg- ert, C. F. Jiao, T. R. Rodriguez, Ab initio nuclear shape coexistence and emergence of island of inversion around n=20, Phys. Lett. B 865 (2025) 139464.doi:10.1016/ j.physletb.2025.139464

    arXiv 2025

  18. [18]

    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, X. Zhang, Ab initio uncertainty quantification of neutrinoless double-beta decay in 76Ge, Phys. Rev. Lett. 132 (2024) 182502. doi:10.1103/PhysRevLett.132.182502. URLhttps://link.aps.org/doi/10.1103/ PhysRevLett.132.182502

    doi:10.1103/physrevlett.132.182502 2024

  19. [19]

    Hayen, Opportunities and open questions in modernβdecay, Ann

    L. Hayen, Opportunities and open questions in modernβdecay, Ann. Rev. Nucl. Part. Sci. 74 (2024) 497–528.arXiv:2403.08485, doi:10.1146/annurev-nucl-121423-100730

    Pith/arXiv arXiv 2024

  20. [20]

    Tichai, R

    A. Tichai, R. Wirth, J. Ripoche, T. Duguet, Symme- try reduction of tensor networks in many-body theory I. Automated symbolic evaluation ofS U(2) algebra, Eur. Phys. J. A 56 (10) (2020) 272.arXiv:2002.05011, doi:10.1140/epja/s10050-020-00233-6

    Pith/arXiv arXiv 2020

  21. [21]

    Kutzelnigg, D

    W. Kutzelnigg, D. Mukherjee, Normal order and extended wick theorem for a multiconfiguration reference wave function, The Journal of Chemical Physics 107 (2) (1997) 432–449.doi:10.1063/1.474405. URLhttps://doi.org/10.1063/1.474405

    doi:10.1063/1.474405 1997

  22. [22]

    L. Kong, M. Nooijen, D. Mukherjee, An algebraic proof of generalized wick theorem, The Journal of Chemi- cal Physics 132 (23) (2010) 234107.doi:10.1063/1. 3439395. URLhttps://doi.org/10.1063/1.3439395

    doi:10.1063/1 2010

  23. [23]

    E. F. Zhou, C. R. Ding, Q. Y . Luo, J. M. Yao, H. Hergert, Ab initio mapping of the boundary of theN=20 island of inversion (3 2026).arXiv:2603.07363. 16

    arXiv 2026