pith. sign in

arxiv: 2604.22223 · v1 · submitted 2026-04-24 · ❄️ cond-mat.str-el · cond-mat.mtrl-sci

QAssemble: A Pure Python Package for Quantum Many-Body Theory

Pith reviewed 2026-05-08 10:01 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.mtrl-sci
keywords QAssemblepure Pythonquantum many-bodyGW approximationdiscrete Lehmann representationHund-Hubbard modelGreen's functionnumerical efficiency
0
0 comments X

The pith

QAssemble shows a pure-Python package can deliver practical efficiency for quantum many-body methods like GW by using modular classes and vectorized kernels with the discrete Lehmann representation.

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

QAssemble is a pure-Python package for solving the quantum many-body problem through functional approaches such as tight-binding, Hartree-Fock, and GW approximations. It organizes the calculation by representing each physical concept—crystal structure, Hamiltonian, Green's function, self-energy—as a separate class in an object-oriented design. Performance-critical steps including the polarizability bubble, Dyson equation inversion, and lattice Fourier transforms are vectorized and combined with the discrete Lehmann representation to reach usable speeds entirely in Python. Benchmarks on a five-orbital extended Hund-Hubbard model show up to 60x speedup relative to traditional loop-based Matsubara implementations, and the code is validated on the electronic structure of graphene with local and non-local interactions. The package supports both batch production runs and interactive method development.

Core claim

The paper introduces QAssemble as a unified pure-Python framework that implements functional methods including tight-binding, Hartree-Fock, and GW approximations. Physical objects are modeled as classes, and performance-critical parts like the polarizability calculation and Dyson equation solution are vectorized and accelerated using the discrete Lehmann representation. This enables practical efficiency for calculations on systems such as graphene with interactions and yields up to 60x speedup on five-orbital extended Hund-Hubbard models compared to standard Matsubara implementations.

What carries the argument

The object-oriented class structure for physical quantities combined with vectorized pure-Python kernels for the polarizability bubble, Dyson equation, and lattice Fourier transforms, paired with the discrete Lehmann representation for frequency handling.

If this is right

  • Researchers gain an accessible platform for implementing and testing new quantum many-body approximations without leaving Python.
  • Calculations involving both local and non-local interactions become feasible on systems like graphene within the same codebase.
  • The reported speedups enable more extensive parameter sweeps and production runs on multi-orbital models.
  • Interactive workflows support rapid prototyping and debugging during method development.
  • The modular class design lowers the effort required to add further approximations or observables.

Where Pith is reading between the lines

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

  • The package could lower the barrier for students and experimentalists who lack expertise in compiled languages to perform many-body calculations.
  • Similar vectorization strategies might be applied to other Python-based scientific codes that currently rely on external compiled libraries.
  • The unified architecture could simplify combining GW with additional methods such as dynamical mean-field theory in a single workflow.
  • Adoption might accelerate community-driven extensions for new physical systems beyond the models demonstrated.

Load-bearing premise

Vectorizing the pure-Python kernels and switching to the discrete Lehmann representation preserves numerical accuracy and stability equivalent to traditional loop-based Matsubara implementations.

What would settle it

A side-by-side run of the five-orbital extended Hund-Hubbard model in QAssemble versus a conventional Matsubara-frequency loop code that produces self-energies or total energies differing by more than floating-point precision.

read the original abstract

QAssemble is a pure-Python package for the quantum many-body problem. It implements various functional approaches, such as tight-binding, Hartree-Fock, and GW approximations within a unified object-oriented architecture. Each physical concept--crystal structure, Hamiltonian, Green's function, self-energy, polarizability, screened Coulomb interaction--is represented as a distinct class. The modular design prioritizes code clarity and extensibility, leveraging NumPy, SciPy, and libdlr for numerical operations. Performance-critical kernels, including the polarizability bubble, Dyson equation inversion, and lattice Fourier transforms, are systematically vectorized and combined with the discrete Lehmann representation to achieve practical efficiency within a pure-Python environment. We validate QAssemble on the electronic structure of graphene with local and non-local interactions. Furthermore, benchmarks on a five-orbital extended Hund-Hubbard model demonstrate that this strategy delivers up to a 60x speedup over traditional loop-based Matsubara implementations. QAssemble supports both batch execution for production calculations and interactive workflows for method development.

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.

Referee Report

1 major / 2 minor

Summary. The manuscript presents QAssemble, a pure-Python package for quantum many-body theory implementing tight-binding, Hartree-Fock, and GW approximations within a unified object-oriented architecture. Physical quantities such as crystal structures, Hamiltonians, Green's functions, self-energies, and polarizabilities are represented as distinct classes. Performance-critical kernels are vectorized using NumPy, SciPy, and libdlr together with the discrete Lehmann representation. Validation is reported on the electronic structure of graphene with local and non-local interactions, and benchmarks on a five-orbital extended Hund-Hubbard model claim up to 60x speedup over traditional loop-based Matsubara implementations. The package supports both batch production runs and interactive workflows.

Significance. If the numerical accuracy of the DLR-vectorized kernels is confirmed to match traditional implementations, the package offers a valuable, accessible tool for quantum many-body calculations in a pure-Python environment. The modular class-based design and emphasis on extensibility could aid method development and education, while the reported speedups would represent a practical contribution for users avoiding compiled extensions.

major comments (1)
  1. [Abstract] Abstract: The central speedup claim of up to 60x on the five-orbital extended Hund-Hubbard model requires that the vectorized pure-Python kernels with discrete Lehmann representation compute identical physical quantities (self-energy, Green's function, total energy) to the same numerical precision as loop-based Matsubara sums. No element-wise comparison of outputs between the two approaches at identical parameters is reported, leaving open the possibility that DLR truncation or reordering introduces additional error.
minor comments (2)
  1. The graphene validation is mentioned but lacks specification of the observables compared (e.g., band structure, density of states) or quantitative error metrics.
  2. Benchmarks would be strengthened by inclusion of error bars on timing data, convergence tests with respect to DLR parameters, and direct methodological comparisons beyond wall-time ratios.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive evaluation and the recommendation of minor revision. The single major comment is addressed point-by-point below; we have revised the manuscript to incorporate the requested validation.

read point-by-point responses
  1. Referee: [Abstract] Abstract: The central speedup claim of up to 60x on the five-orbital extended Hund-Hubbard model requires that the vectorized pure-Python kernels with discrete Lehmann representation compute identical physical quantities (self-energy, Green's function, total energy) to the same numerical precision as loop-based Matsubara sums. No element-wise comparison of outputs between the two approaches at identical parameters is reported, leaving open the possibility that DLR truncation or reordering introduces additional error.

    Authors: We agree that explicit numerical equivalence must be demonstrated to support the speedup claim. Although the DLR truncation error is rigorously bounded by the chosen tolerance and the graphene results are consistent with literature, we acknowledge that an element-wise comparison for the five-orbital benchmark was not included in the original submission. In the revised manuscript we have added a dedicated subsection in the benchmarks section together with a new figure that reports direct, frequency-by-frequency and orbital-by-orbital comparisons of the self-energy, Green's function, and total energy between the vectorized DLR implementation and a reference loop-based Matsubara code at identical parameters (U=2t, J=0.5t, T=0.1t). The maximum relative deviation is 8.7×10^{-13}, well below the DLR tolerance of 10^{-10}, confirming that no additional error is introduced by vectorization or the discrete Lehmann representation. The abstract has also been updated to reference this validation. revision: yes

Circularity Check

0 steps flagged

No significant circularity; performance claims rest on direct external benchmarks.

full rationale

The paper describes a software package implementing standard many-body methods (tight-binding, HF, GW) with vectorized kernels and DLR for efficiency. Central claims are validated on graphene electronic structure and timed against traditional loop-based Matsubara codes on a five-orbital Hund-Hubbard model. These are direct comparisons to independent reference implementations rather than self-referential fits, definitions, or self-citation chains. No derivation reduces to its own inputs by construction; the modular class design and numerical choices are presented as engineering decisions with empirical timing results.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The package implements established approximations using standard numerical libraries without introducing new physical axioms, free parameters, or invented entities.

pith-pipeline@v0.9.0 · 5513 in / 1136 out tokens · 86042 ms · 2026-05-08T10:01:49.910500+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

52 extracted references · 25 canonical work pages

  1. [1]

    Anderson, P. W. More Is Different.Science177, 393–396 (1972)

  2. [2]

    W.More and Different: Notes from a Thoughtful Curmudgeon (World Scientific, 2011)

    Anderson, P. W.More and Different: Notes from a Thoughtful Curmudgeon (World Scientific, 2011)

  3. [3]

    Laughlin, R. B. & Pines, D. The Theory of Everything.Proc. Natl. Acad. Sci. USA97, 28–31 (2000)

  4. [4]

    J., Zhang, X

    Rozenberg, M. J., Zhang, X. Y. & Kotliar, G. Mott-Hubbard transition in infinite dimensions.Physical Review Letters69, 1236–1239 (1992). URL https://link. aps.org/doi/10.1103/PhysRevLett.69.1236

  5. [5]

    Y., Rozenberg, M

    Zhang, X. Y., Rozenberg, M. J. & Kotliar, G. Mott transition in the d=\ensuremath{\infty}Hubbard model at zero temperature.Physical Review Letters70, 1666–1669 (1993). URL https://link.aps.org/doi/10.1103/ PhysRevLett.70.1666. 19

  6. [6]

    Janus-Faced Influence of Hund’s Rule Coupling in Strongly Correlated Materials.Physical Review Letters, 107(25), December 2011

    de’ Medici, L., Mravlje, J. & Georges, A. Janus-Faced Influence of Hund’s Rule Coupling in Strongly Correlated Materials.Physical Review Letters107, 256401 (2011). URL https://link.aps.org/doi/10.1103/PhysRevLett.107.256401

  7. [7]

    M., Kotliar, G., Lee, S.-S

    Stadler, K. M., Kotliar, G., Lee, S.-S. B., Weichselbaum, A. & von Delft, J. Differentiating Hund from Mott physics in a three-band Hubbard-Hund model: Temperature dependence of spectral, transport, and thermodynamic properties. Physical Review B104, 115107 (2021). URL https://link.aps.org/doi/10.1103/ PhysRevB.104.115107

  8. [8]

    Ryee, S., Han, M. J. & Choi, S. Hund Physics Landscape of Two-Orbital Systems. Physical Review Letters126, 206401 (2021). URL https://link.aps.org/doi/10. 1103/PhysRevLett.126.206401

  9. [9]

    URL https://www.pnas.org/doi/full/10.1073/ pnas.2001778117

    Jang, S.et al.Evolution of the Kondo lattice electronic structure above the transport coherence temperature.Proceedings of the National Academy of Sci- ences117, 23467–23476 (2020). URL https://www.pnas.org/doi/full/10.1073/ pnas.2001778117

  10. [10]

    D.et al.Magnetically mediated superconductivity in heavy fermion compounds.Nature394, 39–43 (1998)

    Mathur, N. D.et al.Magnetically mediated superconductivity in heavy fermion compounds.Nature394, 39–43 (1998). URL https://www.nature.com/articles/ 27838

  11. [11]

    URL https://www.nature.com/articles/s41586-020-2122-2

    Jiao, L.et al.Chiral superconductivity in heavy-fermion metal UTe2.Nature 579, 523–527 (2020). URL https://www.nature.com/articles/s41586-020-2122-2

  12. [12]

    Kent, P. R. C. & Kotliar, G. Toward a predictive theory of correlated materi- als.Science361, 348–354 (2018). URL https://www.science.org/doi/10.1126/ science.aat5975

  13. [13]

    Luttinger, J. M. & Ward, J. C. Ground-State Energy of a Many-Fermion System. II.Physical Review118, 1417–1427 (1960). URL https://link.aps.org/doi/10. 1103/PhysRev.118.1417

  14. [14]

    & Kadanoff, L

    Baym, G. & Kadanoff, L. P. Conservation Laws and Correlation Functions. Physical Review124, 287–299 (1961). URL https://link.aps.org/doi/10.1103/ PhysRev.124.287

  15. [15]

    Almbladh, C.-O., Barth, U. V. & Leeuwen, R. V. VARIATIONAL TOTAL ENERGIES FROMϕ- ANDψ- DERIVABLE THEORIES.International Jour- nal of Modern Physics B13, 535–541 (1999). URL https://www.worldscientific. com/doi/abs/10.1142/S0217979299000436

  16. [16]

    Hartree, D. R. & Hartree, W. Self-consistent field, with exchange, for beryllium. Proceedings of the Royal Society of London. Series A - Mathematical and Physical Sciences150, 9–33 (1997). URL https://royalsocietypublishing.org/doi/10.1098/ rspa.1935.0085. Publisher: Royal Society. 20

  17. [17]

    Slater, J. C. The Self Consistent Field and the Structure of Atoms.Physical Review32, 339–348 (1928). URL https://link.aps.org/doi/10.1103/PhysRev.32

  18. [18]

    Publisher: American Physical Society

  19. [19]

    Slater, J. C. Note on Hartree’s Method.Physical Review35, 210–211 (1930). URL https://link.aps.org/doi/10.1103/PhysRev.35.210.2. Publisher: American Physical Society

  20. [20]

    Slater, J. C. A Simplification of the Hartree-Fock Method.Physical Review 81, 385–390 (1951). URL https://link.aps.org/doi/10.1103/PhysRev.81.385. Publisher: American Physical Society

  21. [21]

    Salvat, J

    Froese Fischer, C. General Hartree-Fock program.Computer Physics Com- munications43, 355–365 (1987). URL https://www.sciencedirect.com/science/ article/pii/0010465587900531

  22. [22]

    Hedin ,\ 10.1103/PhysRev.139.A796 journal journal Phys

    Hedin, L. New Method for Calculating the One-Particle Green’s Function with Application to the Electron-Gas Problem.Physical Review139, A796–A823 (1965). URL https://link.aps.org/doi/10.1103/PhysRev.139.A796. Publisher: American Physical Society

  23. [23]

    Hybertsen, M. S. & Louie, S. G. Electron correlation in semiconductors and insu- lators: Band gaps and quasiparticle energies.Physical Review B34, 5390–5413 (1986). URL https://link.aps.org/doi/10.1103/PhysRevB.34.5390. Publisher: American Physical Society

  24. [24]

    & Sham, L

    Godby, R., Schl¨ uter, M. & Sham, L. Self-energy operators and exchange- correlation potentials in semiconductors.Physical Review B37, 10159–10175 (1988)

  25. [25]

    Quantum Computations with Cold Trapped Ions

    Massidda, S., Continenza, A., Posternak, M. & Baldereschi, A. Band-Structure Picture for MnO Reexplored: A Model GW Calculation.Physical Review Letters 74, 2323–2326 (1995). URL https://link.aps.org/doi/10.1103/PhysRevLett.74. 2323

  26. [26]

    & Baldereschi, A

    Massidda, S., Continenza, A., Posternak, M. & Baldereschi, A. Quasiparticle energy bands of transition-metal oxides within a model GW scheme.Physi- cal Review B55, 13494–13502 (1997). URL https://link.aps.org/doi/10.1103/ PhysRevB.55.13494

  27. [27]

    G., J¨ onsson, L

    Aulbur, W. G., J¨ onsson, L. & Wilkins, J. W. inQuasiparticle Calculations in Solids(eds Ehrenreich, H. & Spaepen, F.)Solid State Physics, Vol. 54 1–218 (Academic Press, 2000). URL https://www.sciencedirect.com/science/article/ pii/S0081194708602489

  28. [28]

    inSemiconductors and Insulators(ed.Fulde, P.)Electron Correlations in Molecules and Solids189–221 (Springer, Berlin, Heidelberg, 1995)

    Fulde, P. inSemiconductors and Insulators(ed.Fulde, P.)Electron Correlations in Molecules and Solids189–221 (Springer, Berlin, Heidelberg, 1995). URL https: 21 //doi.org/10.1007/978-3-642-57809-0 9

  29. [29]

    & Vollhardt, D

    Kotliar, G. & Vollhardt, D. Strongly Correlated Materials: Insights From Dynamical Mean-Field Theory.Physics Today 57, 53–59 (2004). URL https://physicstoday.aip.org/features/ strongly-correlated-materials-insights-from-dynamical-mean-field-theory

  30. [30]

    & Rozenberg, M

    Georges, A., Kotliar, G., Krauth, W. & Rozenberg, M. J. Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions. Reviews of Modern Physics68, 13–125 (1996). URL https://link.aps.org/doi/10. 1103/RevModPhys.68.13

  31. [31]

    Reviews of Modern Physics , author =

    Kotliar, G.et al.Electronic structure calculations with dynamical mean-field theory.Reviews of Modern Physics78, 865–951 (2006). URL https://link.aps. org/doi/10.1103/RevModPhys.78.865

  32. [32]

    I., Avella, A

    Anisimov, V. I., Avella, A. & Mancini, F. Electronic structure of strongly cor- related materials.AIP Conf. Proc.1297, 3–134 (2010). URL https://pubs.aip. org/aip/acp/article/1297/1/3-134/854829

  33. [33]

    & Miyake, T

    Imada, M. & Miyake, T. Electronic Structure Calculation by First Principles for Strongly Correlated Electron Systems.Journal of the Physical Society of Japan 79, 112001 (2010). URL https://journals.jps.jp/doi/10.1143/JPSJ.79.112001

  34. [34]

    & Kotliar, G

    Sun, P. & Kotliar, G. Extended dynamical mean-field theory and GW method. Physical Review B66, 085120 (2002). URL https://link.aps.org/doi/10.1103/ PhysRevB.66.085120

  35. [35]

    & Georges, A

    Biermann, S., Aryasetiawan, F. & Georges, A. First-Principles Approach to the Electronic Structure of Strongly Correlated Systems: Combining the G W Approximation and Dynamical Mean-Field Theory.Physical Review Letters90, 086402 (2003). URL https://link.aps.org/doi/10.1103/PhysRevLett.90.086402

  36. [36]

    URL https://linkinghub.elsevier.com/retrieve/ pii/S0010465524003709

    Kang, B.et al.ComDMFT v.2.0: Fully self-consistent ab initio GW+EDMFT for the electronic structure of correlated quantum materials.Computer Physics Com- munications308, 109447 (2025). URL https://linkinghub.elsevier.com/retrieve/ pii/S0010465524003709

  37. [37]

    & Coh, S

    Yusufaly, T., Vanderbilt, D. & Coh, S. Tight-Binding Formalism in the Context of the PythTB Package.arXiv preprint arXiv:1308.4040(2013)

  38. [38]

    URL http: //journals.jps.jp/doi/10.1143/JPSJS.74S.30

    Alet, F.et al.The ALPS Project: Open Source Software for Strongly Correlated Systems.Journal of the Physical Society of Japan74, 30–35 (2005). URL http: //journals.jps.jp/doi/10.1143/JPSJS.74S.30

  39. [39]

    URL https: 22 //linkinghub.elsevier.com/retrieve/pii/S0010465515001666

    Parcollet, O.et al.TRIQS: A toolbox for research on interacting quantum sys- tems.Computer Physics Communications196, 398–415 (2015). URL https: 22 //linkinghub.elsevier.com/retrieve/pii/S0010465515001666

  40. [40]

    R., et al

    Harris, C. R.et al.Array programming with NumPy.Nature585, 357–362 (2020). URL https://doi.org/10.1038/s41586-020-2649-2

  41. [41]

    Virtanen, P.et al.SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python.Nature Methods17, 261–272 (2020)

  42. [42]

    & Strand, H

    Kaye, J., Chen, K. & Strand, H. U. libdlr: Efficient imaginary time calculations using the discrete Lehmann representation.Computer Physics Communica- tions280, 108458 (2022). URL https://linkinghub.elsevier.com/retrieve/pii/ S0010465522001771

  43. [43]

    Slater, J.Quantum Theory of Atomic StructureNo. V. 1 in International series in pure and applied physics (McGraw-Hill, 1960). URL https://books.google.co. kr/books?id=dcJEAAAAIAAJ

  44. [44]

    Electron Correlation and Ferromagnetism of Transition Metals

    Kanamori, J. Electron Correlation and Ferromagnetism of Transition Metals. Progress of Theoretical Physics30, 275–289 (1963). URL https://doi.org/10. 1143/PTP.30.275

  45. [45]

    & Sawatzky, G

    Van Der Marel, D. & Sawatzky, G. A. Electron-electron interaction and localiza- tion indandftransition metals.Physical Review B37, 10674–10684 (1988). URL https://link.aps.org/doi/10.1103/PhysRevB.37.10674

  46. [46]

    Strand, H. U. R. Valence-skipping and negative- U in the d -band from repulsive local Coulomb interaction.Physical Review B90, 155108 (2014). URL https: //link.aps.org/doi/10.1103/PhysRevB.90.155108

  47. [47]

    Some remarks on the Pariser-Parr-Pople method.Theoretica chimica acta2, 219–227 (1964)

    Ohno, K. Some remarks on the Pariser-Parr-Pople method.Theoretica chimica acta2, 219–227 (1964). URL https://doi.org/10.1007/BF00528281

  48. [48]

    & Tremblay, A.-M

    Reymbaut, A., Bergeron, D. & Tremblay, A.-M. S. Maximum entropy analytic continuation for spectral functions with nonpositive spectral weight.Phys. Rev. B 92, 060509 (2015). URL https://link.aps.org/doi/10.1103/PhysRevB.92.060509

  49. [49]

    & Choi, H

    Han, M. & Choi, H. J. Parameter-free analytic continuation for quantum many- body calculations.Phys. Rev. B106, 245150 (2022). URL https://link.aps.org/ doi/10.1103/PhysRevB.106.245150

  50. [50]

    J., Thibeault, J

    Hay, P. J., Thibeault, J. C. & Hoffmann, R. Orbital interactions in metal dimer complexes.Journal of the American Chemical Society97, 4884–4899 (1975). URL https://pubs.acs.org/doi/abs/10.1021/ja00850a018

  51. [51]

    & Johnson, S

    Frigo, M. & Johnson, S. FFTW: an adaptive software architecture for the FFT. Proc. 1998 IEEE Int. Conf. Acoust. Speech Signal Process.3, 1381–1384 (1998). URL http://ieeexplore.ieee.org/document/681704/. 23

  52. [52]

    Σ# TB𝐻! GW 𝐺! 𝐺 Σ$%&𝑃 𝑊 𝑉 𝐻! Σ

    Kaye, J., Chen, K. & Parcollet, O. Discrete Lehmann representation of imaginary time Green’s functions.Physical Review B105, 235115 (2022). URL https: //link.aps.org/doi/10.1103/PhysRevB.105.235115. 24 (a) (b) Extended Data Fig. 1 Diagrammatic representation of the Hartree–Fock self-energies. (a) The Hartree self-energy is depicted as a single fermion lin...