pith. sign in

arxiv: 2410.09587 · v2 · submitted 2024-10-12 · ❄️ cond-mat.mtrl-sci · math-ph· math.MP

A polynomial GCD certificate for exact flat bands in finite-range Bloch Hamiltonians

Ivan Damnjanovi\'c , Milan Damnjanovi\'c , Ivanka Milo\v{s}evi\'c , Dragan Stevanovi\'c This is my paper

Pith reviewed 2026-05-23 18:40 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci math-phmath.MP
keywords flat bandsBloch Hamiltoniantight-binding modelspolynomial GCDcharacteristic polynomialperiodic latticesdispersionless bands
0
0 comments X

The pith

The monic GCD of the coefficient polynomials in the Bloch Hamiltonian's characteristic polynomial identifies the exact flat band energies.

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

This paper formulates a polynomial GCD certificate for exact flat bands in finite-range periodic tight-binding Hamiltonians. By writing the characteristic polynomial of the Bloch Hamiltonian as a Laurent polynomial in the wavevector variables, the monic GCD of its energy-dependent coefficients is shown to isolate the largest factor that depends only on energy. Its roots mark the flat band energies and their multiplicities give the common algebraic multiplicities across the entire Brillouin zone. A sympathetic reader would care because the approach supplies a symbolic computation tool for tuning hopping parameters to produce desired flat bands without pointwise diagonalization. The certificate is invariant under unit cell choice and Bloch gauge.

Core claim

Writing the characteristic polynomial of the Bloch Hamiltonian as a Laurent polynomial P_L(z, λ) = sum_t c_t(λ) z^t, the monic greatest common divisor G_L(λ) = gcd_t c_t(λ) is precisely the maximum factor of P_L that depends only on the energy variable. Its roots are exactly the exact flat-band energies, and their multiplicities give common algebraic multiplicities of these flat bands throughout the Brillouin zone.

What carries the argument

The monic greatest common divisor G_L(λ) of the coefficients c_t(λ) in the Laurent polynomial form of the characteristic polynomial.

If this is right

  • The roots of G_L(λ) are exactly the exact flat-band energies.
  • Multiplicities in G_L(λ) give the common algebraic multiplicities of the flat bands throughout the Brillouin zone.
  • The certificate is invariant under choice of unit cell and Bloch gauge.
  • It functions as a symbolic tool for engineering hopping parameters in lattice models.
  • It applies to examples including kagome, dice, and octahedron-chain lattices.

Where Pith is reading between the lines

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

  • The GCD step can be performed once symbolically before any numerical diagonalization of the Hamiltonian.
  • The method isolates the detection of dispersionless eigenvalues from the separate eigenvector analysis required for compact localized states or topology.

Load-bearing premise

The characteristic polynomial of the Bloch Hamiltonian admits a Laurent polynomial representation whose coefficients share common factors independent of the wavevector precisely when exact flat bands exist.

What would settle it

For the kagome lattice model known to possess an exact flat band at a specific energy E, compute G_L(λ) and verify whether (λ - E) divides it with the stated multiplicity.

read the original abstract

We formulate a polynomial GCD certificate for exact flat bands in finite-range periodic tight-binding Hamiltonians. Writing the characteristic polynomial of the Bloch Hamiltonian as a Laurent polynomial \( P_L(\mathbf{z},\lambda)=\det(\lambda I-H_B(\mathbf{z}))=\sum_{\mathbf{t}}c_{\mathbf{t}}(\lambda)\mathbf{z}^{\mathbf{t}}, \) we show that the monic greatest common divisor \(G_L(\lambda)=\gcd_{\mathbf{t}}c_{\mathbf{t}}(\lambda)\) is precisely the maximum factor of \(P_L\) that depends only on the energy variable. Its roots are exactly the exact flat-band energies, and their multiplicities give common algebraic multiplicities of these flat bands throughout the Brillouin zone. The coefficient-vanishing criterion underlying this statement is known in the flat-band and periodic-graph literature; the contribution emphasized here is the compact GCD formulation, its unit cell and Bloch-gauge invariance, and its use as a symbolic computation tool for hopping parameter engineering. The method is illustrated on kagome, dice and octahedron-chain examples, including weighted kagome and dice lattices. The certificate detects exact dispersionless eigenvalues; compact localized states, band touching and topological character must be analyzed in a subsequent eigenvector or projector calculation.

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

0 major / 3 minor

Summary. The manuscript formulates a polynomial GCD certificate for exact flat bands in finite-range periodic tight-binding Hamiltonians. It expresses the characteristic polynomial of the Bloch Hamiltonian as the Laurent polynomial P_L(z, λ) = det(λI − H_B(z)) = ∑_t c_t(λ) z^t and shows that the monic GCD G_L(λ) = gcd_t c_t(λ) is the largest factor of P_L independent of z; its roots are the exact flat-band energies and its multiplicities give the common algebraic multiplicities of those eigenvalues over the Brillouin zone. The construction is shown to be invariant under unit-cell redefinitions and Bloch-gauge choices, and is illustrated on kagome, dice, and octahedron-chain lattices (including weighted variants) as a symbolic tool for hopping-parameter engineering. The underlying coefficient-vanishing criterion is noted as previously known; the contribution is the compact GCD formulation.

Significance. If the algebraic claim holds, the method supplies a parameter-free, machine-checkable certificate that directly extracts flat-band energies and their minimal multiplicities from the characteristic polynomial without requiring explicit diagonalization or eigenvector analysis at each k-point. It is directly implementable in computer-algebra systems, separates flat-band detection from subsequent compact-localized-state or topological analysis, and applies uniformly to any finite-range periodic graph whose Bloch Hamiltonian yields a Laurent polynomial. These features make it a practical addition to the flat-band literature for both theoretical classification and material-design workflows.

minor comments (3)
  1. [Method section (near the definition of G_L)] The manuscript states that the GCD construction is unit-cell and Bloch-gauge invariant but does not display the explicit transformation law for the coefficients c_t under a change of unit cell; adding one short paragraph or appendix equation would make the invariance claim self-contained.
  2. [Examples section] In the weighted-kagome and weighted-dice examples, the symbolic GCD output is reported but the intermediate coefficient polynomials c_t(λ) are not listed; including them (even in an appendix) would allow readers to reproduce the GCD step by hand.
  3. [Introduction] The abstract and introduction cite the coefficient-vanishing criterion as known in the flat-band literature; adding one or two explicit references to the original sources of that criterion would strengthen the novelty statement.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading and positive evaluation of the manuscript. The recommendation to accept is appreciated. No major comments were raised that require point-by-point rebuttal or manuscript changes.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper's central claim is that G_L(λ) = gcd_t c_t(λ) extracts the flat-band energies from the characteristic polynomial P_L(z, λ) = det(λI - H_B(z)) written in Laurent form. This follows directly from the algebraic fact that a nonzero Laurent polynomial vanishes identically on the torus if and only if every coefficient vanishes, so the common roots of the c_t(λ) are exactly the energies where P_L(z, E) ≡ 0 for all z. The paper explicitly states that the underlying coefficient-vanishing criterion is already known in the literature and presents the GCD only as a compact, invariant reformulation for symbolic computation. No fitted parameters, self-referential predictions, uniqueness theorems imported from the authors' prior work, or ansatzes appear; all steps are standard facts about polynomial rings and determinants with no reduction to the paper's own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The method rests on standard facts of polynomial algebra over Laurent polynomial rings; no free parameters, ad-hoc axioms, or new entities are introduced in the abstract.

axioms (2)
  • standard math The characteristic polynomial of a finite-range Bloch Hamiltonian can be expressed as a Laurent polynomial in the momentum variables z with coefficients that are polynomials in the energy variable λ.
    Invoked when writing P_L(z, λ) = sum_t c_t(λ) z^t; this is a standard property of finite-range tight-binding models.
  • standard math The monic GCD of the coefficient polynomials c_t(λ) extracts exactly the common energy-dependent factors independent of z.
    Core algebraic statement underlying the certificate; follows from properties of polynomial rings.

pith-pipeline@v0.9.0 · 5783 in / 1560 out tokens · 21760 ms · 2026-05-23T18:40:08.062570+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

30 extracted references · 30 canonical work pages

  1. [1]

    Zhang, Y

    Z. Zhang, Y. Wang, K. Watanabe, T. Taniguchi, K. Ueno, E. Tutu c, B. J. LeRoy, Nature Phys. 16 (2020) 1093

    work page 2020

  2. [2]

    S. Lisi, X. Lu, T. Benschop, T. A. de Jong, P. Stepanov, J. R. Du ran, F. Margot, I. Cucchi, E. Cappelli, A. Hunter, A. Tamai, V. Kandyba, A. Giampietr i, A. Barinov, J. Jobst, V. Stalman, M. Leeuwenhoek, K.Watanabe, T. T aniguchi, L. Rademaker et al., Nature Phys. 17 (2021) 189

    work page 2021

  3. [3]

    Danieli, A

    C. Danieli, A. Andreanov, and S. Flach, Phys. Rev. B 102 (2020) 041116(R)

    work page 2020

  4. [4]

    M. Kang, L. Ye, S. Fang, J. S. You, A. Levitan, M. Han, J. I. Fac io, C. Jozwiak, A. Bostwick, E. Rotenberg, M. K. Chan, R. D. McDonald, D. Graf, K. K aznatcheev, E. Vescovo, D. C. Bell, E. Kaxiras, J. van den Brink, M. Richter, M. P rasad Ghimire et al., Nature Mater. 19 (2020) 163 . 14

    work page 2020

  5. [5]

    Z. Liu, H. Mao, Y. Nakajima, J. Paglione, Z. Yin, S. Borisenko, Phy s. Rev. B 109 (2024) 075103

    work page 2024

  6. [6]

    D.-S.Ma, Y. Xu, C. S. Chiu, N. Regnault, A. A. Houck, Z. Song, B. A . Bernevig, Phys. Rev. Lett. 125 (2020) 266403

    work page 2020

  7. [7]

    Calugaru, A

    D. Calugaru, A. Chew, L. Elcoro, Y. Xu, N. Regnault, Z.-D. Song, B. A. Bernevig, Nat. Phys. 18 (2022) 185

    work page 2022

  8. [8]

    Y.-F. Wang, H. Yao, C.-D. Gong, D. N. Sheng, Phys. Rev. B 86 (2012) 201101(R)

    work page 2012

  9. [9]

    Y. Yang, C. Roques-Carmes, S. E. Kooi, H. Tang, J. Beroz, E. M azur, I. Kaminer, J. D. Joannopoulos, M. Soljacic, Nature 613 (2023) 42

    work page 2023

  10. [10]

    F. Xie, Z. Song, B. Lian, B. A. Bernevig, Phys. Rev. Lett. 124 (2020)167002

    work page 2020

  11. [11]

    J. M. Lee, C. Geng, J. W. Park, M. Oshikawa, S.-S. Lee, H.W. Yeo m,G. Y. Cho, Phys. Rev. Lett. 124 (2020) 137002

    work page 2020

  12. [12]

    Balents, C

    L. Balents, C. R. Dean, D. K. Efetov, A. F. Young, Nat. Phys. 16 (2020) 725–733

    work page 2020

  13. [13]

    M. M. Al Ezzi, J. Hu, A. Ariando, F. Guinea, S. Adam, Phys. Rev. Lett. 132 (2024) 126401

    work page 2024

  14. [14]

    Lieb, Phys

    E.H. Lieb, Phys. Rev. Lett. 62 (1989) 1201

    work page 1989

  15. [15]

    Calugaru, A

    D. Calugaru, A. Chew, L. Elcoro, Y. Xu, N. Regnault, Z.-D. Song , B.A. Bernevig, Nat. Phys. 18 (2022) 185–189

    work page 2022

  16. [16]

    A. J. Kollar, M. Fitzpatrick, P. Sarnak, A. A. Houck, Commun. M ath. Phys. 376 (2020) 1909

    work page 2020

  17. [17]

    C. S. Chiu, A. N. Carroll, N. Regnault, A. A. Houck, Phys. Rev. R es. 4 (2022) 023063

    work page 2022

  18. [18]

    Morfonios, M

    C.V. Morfonios, M. Rontgen, M. Pyzh, P. Schmelcher, Phys. Re v. B 104 (2021) 035105

    work page 2021

  19. [19]

    Bradlyn, L

    B. Bradlyn, L. Elcoro, J. Cano, M. G. Vergniory, Z. Wang, C. Fe lser, M. I. Aroyo, B. A. Bernevig, Nature 547 (2017) 298

    work page 2017

  20. [20]

    H. C. Po, A. Vishwanath and H. Watanabe, Nat. Commun. 8 (2017) 50

    work page 2017

  21. [21]

    L. Chen, T. Mazaheri, A. Seidel, X. Tang, J. Phys. A: Math. The or. 47 (2014) 152001

    work page 2014

  22. [22]

    C. S. Chiu, D.-S. Ma, Z.-D. Song, B. A. Bernevig, A. A. Houck, Ph ys. Rev. Res. 2 (2020) 043414. 15

    work page 2020

  23. [23]

    Sutherland, Phys

    B. Sutherland, Phys. Rev. B 34 (1986) 5208

    work page 1986

  24. [24]

    B. Sturmfels, Solving systems of polynomial equations, CBMS Re gional Confer- ences in Mathematics 97, American Mathematical Society, Providen ce, RI, USA, 2002, available at https://math.berkeley.edu/ bernd/cbms.pdf

    work page 2002

  25. [25]

    Ramachandran, A

    A. Ramachandran, A. Andreanov and S. Flach, Phys. Rev. B 96 (2017) 161104(R)

    work page 2017

  26. [26]

    Leykam, A

    D. Leykam, A. Andreanov and S. Flach, Advances in Physics: X 3 (2018) 1473052

    work page 2018

  27. [27]

    T. X. Hoang, D. Leykam and Y. Kivshar, Phys. Rev. Lett. 132 (2024) 043803

    work page 2024

  28. [28]

    Miloˇ sevi´ c, S

    I. Miloˇ sevi´ c, S. Dmitrovi´ c, T. Vukovi´ c, A. Dimi´ c, M. Damnjanovi´ c, J. Phys. A: Math. Theor. 53 (2020) 455204

    work page 2020

  29. [29]

    Dmitrovi´ c, T

    S. Dmitrovi´ c, T. Vukovi´ c, I. Miloˇ sevi´ c, M. Damnjanovi´c, J. Phys. A: Math. Theor. 55 (2022) 385201

    work page 2022

  30. [30]

    Damnjanovi´ c, Topological phases of layers: elementary band representat ions, NANOTEXNOLOGY 2020 (Thessaloniki, Greece 4 – 11 July 2020)

    M. Damnjanovi´ c, Topological phases of layers: elementary band representat ions, NANOTEXNOLOGY 2020 (Thessaloniki, Greece 4 – 11 July 2020). 16

    work page 2020