An arXiv preprint reports a numerical method for tracking band structures in a model with periodic potentials, including a sharply bracketed change in a band's Chern number, an integer quantity used in the study to track band topology. In the non-degenerate calculation, the lower and midpoint evaluations had Chern number 1 while the upper evaluation had 0, placing the reported transition between 3.168333 and 3.168335. The paper cautions that finite numerical resolution leaves the exact gap closure unconfirmed.
Building the numerical picture
The method is a Galerkin projection, also called Rayleigh-Ritz. It uses a finite-dimensional subspace and a basis built from Landau-operator eigenfunctions to approximate eigenpairs, the numerical energy-and-state solutions sought by the calculation. The research question is whether that spectral approach can compute band structures across weak, intermediate and strong periodic-potential regimes.
For the periodic potential, the paper uses a Fourier expansion, and its conclusions describe use of the Fast Fourier Transform. The resulting band structures are paired with Chern-number calculations made with the Fukui-Hatsugai-Suzuki method. Together, these outputs give the study a numerical view of both band structure and topology across the tested regimes.
The reported common setup used an 11 x 11 k-point grid, theta-function truncation of 10, a Fourier cutoff of 26 and 100 path points per segment. The band-structure calculations used 50 states, while the Chern-number calculations used 30 states. Those settings define the numerical resolution of the reported band and topology calculations.
A very narrow non-degenerate transition
The integer-flux tests compare one non-degenerate configuration, beta = 1 and B = 2pi, with one degenerate configuration, beta = 4 and B = 8pi. In each, the stated task is to compute band structures and associated Chern numbers for the selected configuration.
In the weak-potential calculation, the reported Chern number was 1 within machine precision. In the final non-degenerate bisection row, the lower and midpoint Chern values were 1, while the upper value was 0. The final bracket ran from 3.168333 to 3.168335.
The bracket and the gap-closure warning answer different questions. The bracket records where the calculated Chern value changed in the bisection; the paper says that numerical approximations and discretization prevent definitive confirmation of an exact gap closure at V0/|B| = 3.2. The transition is therefore a numerical estimate rather than a settled exact crossing.
Splitting and localization in the degenerate case
In the degenerate system, the calculations showed lifted degeneracy and band splitting. The third and fourth bands crossed between 3.8 and 4.0, and the reported transition was near 3.90355. Across that transition, the total Chern number changed from 1 to 0.
At high potential strength, the bands were reported to flatten again at V0/|B| = 20, while the ground-state density became highly localized at the potential wells. This extends the calculation beyond a topological label: it also reports the spatial concentration of the lowest-state density.
Rational flux changes the accounting
The rational-flux simulation fixes B = 6pi, uses a supercell flux of 3 and reports degeneracy of degree 3 for each state. In the weak regime, the per-band Chern number is 1 and the total Chern number is 3.
In the intermediate regime, C2 changed from 1 to -2 while the total Chern number stabilized at 0. At a later crossing, C2 returned to 1 and C3 changed from 1 to -2. The example shows why the reported results give both individual and total Chern numbers.
Precision without an exact closure
Because the calculation uses a finite basis, truncated theta and Fourier representations and a finite momentum grid, the reported bands and transition locations are approximate. The paper does not provide a definitive exact-gap result at V0/|B| = 3.2.
The evidence is confined to the tested weak, intermediate and strong regimes and to the selected integer- and rational-flux configurations. It therefore does not establish that the reported transition ratios would remain unchanged at a different numerical resolution.
Taken together, the reported calculations show a Galerkin method producing band structures and changing Chern-number assignments across the tested regimes, while leaving exact transition details dependent on numerical resolution. The document is an arXiv preprint dated 20 August 2026. It states that all data and code used to generate the reported results are available through its listed repository.
Publication and funding
The paper acknowledges DFG support under project 516782692 and Excellence Strategy funding EXC2075 - 390740016, as well as support from SimTech.
Paper data and sources
Original title: Diagonalization of the Landau Hamiltonian with a Periodic Potential via a Galerkin Projection Method and Applications to Topological Band Properties
Authors: Rafael Antonio Lainez Reyes, Benjamin Stamm, Hans Peter Büchler
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-20
DOI: Not available
Original paper · Full text