Pith. sign in

REVIEW 2 major objections 4 minor 26 references

The paper establishes that every three-dimensional Bravais lattice has a well-defined Baldereschi point: the single wavevector at which smooth, symmetric functions of wavevector are, on average, closest to their Brillouin-zone mean, and tab

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-01 19:23 UTC pith:QSWASEKD

load-bearing objection A clean, parameter-free tabulation of Baldereschi points for all 14 Bravais lattices, with a known gap in the reducibility test that is unlikely to matter in practice. the 2 major comments →

arxiv 2607.16974 v1 pith:QSWASEKD submitted 2026-07-18 cond-mat.mtrl-sci

Baldereschi mean value points for three-dimensional Bravais lattices

N. D. Drummond This is my paper
classification cond-mat.mtrl-sci
keywords Baldereschi pointmean-value pointBravais latticeBrillouin zonequantum Monte Carlotwist averagingstar functionsfinite-size effects
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 establishes that every three-dimensional Bravais lattice has a well-defined Baldereschi point: the single wavevector at which smooth, symmetric functions of wavevector are, on average, closest to their Brillouin-zone mean. It provides formulas and tables for all fourteen lattice types, using a hierarchical procedure that identifies which star functions can be independently set to zero or minimized. This matters for quantum Monte Carlo simulations, where each calculation can use only one Bloch wavevector; the Baldereschi point is the optimal single twist. The paper also demonstrates with three test systems that the tabulated point captures twist-averaged results far better than the Gamma point.

Core claim

The central claim is that for any Bravais lattice, the Baldereschi point is determined purely by the lattice through the star functions A_n(k) = sum of exp(i k·R) over the nth shell of lattice points, and that the previously incomplete set of known points can be completed. By constructing a hierarchy of target star functions — first those whose union spans three-dimensional space, then the shortest stars that cannot be written as combinations of already-chosen targets — the paper derives the point for all fourteen Bravais lattices. In the cubic, tetragonal, orthorhombic, monoclinic, and triclinic families the point often has a simple rational form such as (1/4, 1/4, 1/4) in reciprocal lattic

What carries the argument

The mechanism is the star-function hierarchy. A star function A_n(k) is the sum of plane waves over a complete shell of symmetry-related lattice vectors; these functions form an orthogonal basis for symmetric periodic functions. The Baldereschi point is found by zeroing or minimizing the first few target star functions in order of radius, under the constraint that a star whose representation is reducible into lower-radius stars cannot be varied independently. Equation (4) — (2 r_n)^2 = sum (n_i r_i)^2 — is the test that decides whether a star carries an independent irreducible component; the algorithm in appendix A uses it to build the target list, and the numerical scheme solves the resulti

Load-bearing premise

The load-bearing premise is that equation (4) reliably detects when a star function cannot be varied independently; if an accidental degeneracy or an accidentally satisfied radius relation makes that test fail, the corresponding table entry is not the true Baldereschi point.

What would settle it

For a lattice parameter ratio near a predicted crossover (for example rhombohedral c/a = sqrt(6) or sqrt(3/2)), compute the ensemble-averaged squared deviation sum_n s_n^2 A_n(k)^2 over a fine grid of k in the Brillouin zone and compare its global minimum with the table's point; any mismatch, or any non-unique minimum, would falsify that entry. Alternatively, test a hexagonal lattice with c=a, where the paper notes equation (4) can be invalidated, and check whether the reported point still minimizes the deviation.

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

If this is right

  • For insulating crystals, a QMC calculation at the tabulated Baldereschi twist should reproduce twist-averaged energies to within a small fraction of the spread across twists, eliminating the need for multiple twist calculations in most cases.
  • The tables resolve prior inconsistencies in the literature: the orthorhombic Baldereschi point is exactly (1/4,1/4,1/4), and the rhombohedral point has different intermediate-angle behavior than an earlier report.
  • The target-star algorithm applies beyond the fourteen lattices, so Baldereschi points can be generated for arbitrary lattice geometries, supercells, or lower-dimensional systems by running the same identification procedure.
  • The free-electron-gas test shows that Baldereschi twists reduce momentum-quantization shell oscillations by about an order of magnitude relative to Gamma twists, though twist averaging remains necessary in metals.
  • Because the point depends only on the Bravais lattice, the tables apply to any crystal or supercell with a given lattice, irrespective of atomic basis.

Where Pith is reading between the lines

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

  • The same hierarchical construction could be extended to define Baldereschi sets — small sets of k-points that jointly minimize the ensemble-averaged deviation — which might outperform existing special-point grids for intermediate-size twist sets.
  • The reliance on equation (4) could be made rigorous by replacing the radius test with an explicit decomposition of each star under the lattice point group; this would remove the residual accidental-degeneracy caveat.
  • For metallic systems, the Baldereschi point's role as the single best twist for wave-function optimization suggests a combined protocol: optimize at the Baldereschi point, then twist-average the final energy — a scheme the paper motivates but does not test directly.
  • If the target-star algorithm is implemented for a specific material, the resulting point can be checked by computing the DFT band-structure energy at nearby k-points and confirming it sits at a local minimum of the mean-squared deviation, a quick validation for any code.

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

2 major / 4 minor

Summary. The paper defines the Baldereschi mean-value point of a 3D Bravais lattice as the wavevector obtained by hierarchically zeroing or minimizing the first few star functions of the lattice-point Fourier expansion, and it presents a general algorithm for identifying the target stars. Closed-form tables are given for all fourteen 3D Bravais lattices, with piecewise parameter ranges where the target-star ordering changes. The authors cross-check their analytic results with a numerical program, and they validate the practical value of the points on hBN (DFT), monolayer graphene (QMC), and the free electron gas (analytical).

Significance. If correct, the tables would provide a standard reference for choosing a single representative twist in QMC supercell calculations for any 3D lattice geometry. The derivation of the star-function minimization procedure (§3.1–3.3) is clean and parameter-free, and the tabulated points are derived from the definition rather than fitted. The external validations in §5 give concrete evidence that the points are effective in practice. However, the correctness of every table entry rests on the reducibility test of Eq. (4), which is used as a biconditional but for which only necessity is proved; this is the main correctness risk.

major comments (2)
  1. [§3.2, Eq. (4); Appendix A, step 4(b)] The algorithm treats Eq. (4) as a sufficient test for whether a star carries a reducible representation: if (2rn)^2 = Σ (nt rt)^2 for some integers, the star is assumed reducible and is not selected as a target. The derivation of Eq. (4) only gives necessity. The paper acknowledges accidental degeneracies (e.g., hexagonal c=a) and states that they do not appear to be a problem in numerical tests, but this is a spot-check, not a proof or an exhaustive scan. Because every table entry is determined by the target-star identification, a single misclassification would change the reported point. The authors should either prove sufficiency for the finite set of low-radius stars occurring in the 14 lattice types, or provide an exhaustive case analysis (e.g., by direct decomposition of the permutation representation for each star), or explicitly restrict the central claim to generic parameters.
  2. [§4, Table 1, rhombohedral rows] The rhombohedral tables are piecewise with strict inequalities, excluding c/a = √6, √(3/2), 1 and the corresponding rhombohedral angles α = π/3, π/2, α_d, 2π/3. At these values the lattice has higher symmetry (e.g., bcc, fcc, simple-cubic-like limits), and the two adjacent piecewise expressions may not agree in the limit, so the Baldereschi point at the degenerate ratio is left undefined. The abstract's claim that all fourteen Bravais lattices are tabulated is therefore not literally satisfied for these measure-zero but physically realizable parameter values. The authors should either give the points at the degenerate ratios (possibly as limits, if they are unique) or state clearly that the tables cover only nondegenerate ratios and explain how to handle the special cases.
minor comments (4)
  1. [§3.3, after Eq. (1)] The text introduces the hexagonal lattice with lattice vectors (a,0,0), (−a/2,√3a/2,0), (0,0,c), but the displayed A1(k) and A2(k) are for the alternative basis (a,0,0), (a/2,√3a/2,0), (0,0,c). This switch is easy for a reader to miss, and it is the source of the two distinct hexagonal rows in Table 1. Please label the two conventions explicitly.
  2. [Figure 1] The caption says 'Baldereschi points' (plural) but the text refers to a single Baldereschi point; the figure appears to show symmetry-equivalent copies. A short phrase stating that symmetry-equivalent points are shown would avoid confusion.
  3. [Table 1 caption, At3 column] The meaning of the At3(kb) column is not defined in the caption. It is apparently the value of the third target star function when it is minimized (or 0 when it is set to zero). Please state this explicitly, since the same column contains numerical values (e.g., 4.404, −1.608, −3) and zeros.
  4. [§3.1, Eq. (2) and following] The claim that the Baldereschi point minimizes the mean-squared deviation for sufficiently large λ is presented as intuitive; a formal statement (e.g., uniform convergence of the minimizer as λ→∞ under the stated conditions) would strengthen the definition.

Circularity Check

0 steps flagged

No significant circularity: the tabulated Baldereschi points are computed from the definition via star-function targets, with only peripheral self-citations; the Eq. (4) caveat is a correctness risk, not circularity.

full rationale

The paper's central claim — the tabulated Baldereschi points for all fourteen 3D Bravais lattices — is produced by applying its own definition: minimize the ensemble-averaged mean-squared deviation (Eq. 2) by identifying target star functions (Appendix A) and solving A_t1 = A_t2 = 0 while zeroing or minimizing A_t3 (Section 3.4). These quantities are computed from lattice geometry and symmetry alone; no fitted parameters or data-dependent inputs feed into the table values. The Section 5 validation (hBN DFT, graphene DMC, free electron gas) is external to the derivation and checks rather than defines the points. Self-citations (e.g., Refs. [12], [18], [19], [25]) appear in peripheral methodological or software contexts and are not load-bearing for the target-star identification or the tabulated wavevectors. The acknowledged caveat in Section 3.2 — that Eq. (4) may be accidentally satisfied under special lattice-parameter ratios — is a robustness/completeness limitation, not a circular step, because the derivation does not assume the target result as an input. Overall, no specific reduction of the claimed output to its inputs by construction is present.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 0 invented entities

The paper is free of fitted parameters: every tabulated point is a closed-form expression in the lattice vectors, cross-checked numerically. The load-bearing assumptions are the smoothness/symmetry framework for observables as functions of twist, the exponential-decay ensemble model that motivates the hierarchical definition, and the reducibility test (Eq. 4) whose robustness is asserted from numerical experience. No new physical entities are introduced; 'target stars' are an algorithmic construct, not an invented physical object.

axioms (5)
  • domain assumption The mean squared deviation of a symmetric periodic function from its BZ mean is exponentially dominated by the first star function for sufficiently large λ (Eq. 2), justifying the hierarchical zeroing/minimization of star functions.
    §3.1, Eq. (2). The ensemble model (zero-mean coefficients, σ_n = s0 e^−λn) motivates the definition; the paper concedes in §5 that for slowly convergent series the Baldereschi point's effectiveness is empirical, not guaranteed.
  • domain assumption Equation (4): a star carries a reducible representation of the point group iff its radius satisfies (2r_n)^2 = Σ n_i^2 r_i^2 against already-identified target stars.
    §3.2; used in appendix A step 4(b). Derived from the irrep orthogonality theorem, but the paper admits accidental degeneracies (hexagonal c = a) and possible accidental satisfaction can invalidate it; robustness is asserted from numerical tests, not proven.
  • domain assumption In insulators, QMC observables are smooth functions of the twist k_s; in metals, occupancy changes create discontinuities.
    §2 and §5.3. The optimality claim for the Baldereschi point is restricted to the smooth case; the paper shows explicitly that metals still require twist averaging (N^−1 vs N^−4/3 scaling).
  • standard math Star functions form an orthogonal basis for symmetric functions and carry point-group representations; stars spanning independent subspaces can be varied independently.
    §3.1–3.2: Fourier expansion (Eq. 1) and the orthogonality theorem for irreps; standard group-theoretic framework for BZ symmetry.
  • domain assumption Newton-Raphson/BFGS, started from Sobol'-sampled initial points, converge to the global constrained minimum rather than a spurious stationary point.
    §3.4 and appendix B. The paper reports nonglobal minima only for the rhombohedral lattice, and notes a handful of starts suffice — an empirical claim rather than a guarantee.

pith-pipeline@v1.3.0-alltime-deepseek · 16422 in / 19840 out tokens · 180506 ms · 2026-08-01T19:23:13.832751+00:00 · methodology

0 comments
read the original abstract

The Baldereschi point of a crystal is a wavevector in the Brillouin zone at which every smooth periodic function of wavevector lies close to its mean value. Although originally introduced in the context of one-electron methods, mean-value points are ideal for explicitly correlated many-electron methods such as quantum Monte Carlo simulations, which can only use a single Bloch wavevector in the Brillouin zone of a simulation supercell. We have therefore evaluated and tabulated the Baldereschi mean-value points of all fourteen three-dimensional Bravais lattices.

Figures

Figures reproduced from arXiv: 2607.16974 by N. D. Drummond.

Figure 1
Figure 1. Figure 1: Nodes of the first two nonzero star functions A1(k) and A2(k) of a hexagonal lattice with c > √ 3a within the first BZ (the cyan hexagon). The Baldereschi points are shown by the filled circles. 5 [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Difference between the DFT-PBE energy of bulk hBN with a n × n × n k-point grid and the DFT-PBE energy with a large 123 × 123 × 40 k-point grid. Results are shown for different twists ks (i.e., offsets to the grid of k points). The labels “Γ”, “(1/4, 1/4, 1/4)”, and “Baldereschi” indicate the position of ks within the Brillouin zone of the unfolded n × n × n supercell. The twist-averaged results are averag… view at source ↗
Figure 3
Figure 3. Figure 3: Kinetic energy per electron T(N) of a two-component free electron gas in a simple cubic simulation cell against the reciprocal of the system size N. Results are shown for finite simulation cells subject to pure periodic boundary conditions, twisted periodic boundary condition with the twist at the simulation cell Baldereschi point, and twist-averaging. The kinetic energies with pure periodic boundary condi… view at source ↗
Figure 4
Figure 4. Figure 4: confirms that finite-size errors in the twist-averaged kinetic energy per electron fall off as N −4/3 [6], whereas the typical finite-size errors with periodic (Γ) and twisted periodic (Baldereschi) points fall off as N −1 , despite the fact that the latter has a much smaller prefactor. Hence, although the Baldereschi point is clearly vastly preferable to the Γ point in a finite cell subject to twisted per… view at source ↗

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

26 extracted references · 12 canonical work pages

  1. [1]

    Baldereschi A 1973 Phys. Rev. B 7(12) 5212–5215 URL https://link.aps.org/doi/10.1103/PhysRevB.7.5212

    doi:10.1103/physrevb.7.5212 1973

  2. [2]

    Status Solidi B 119 9–40 URL https://onlinelibrary.wiley.com/doi/abs/10.1002/pssb.2221190102

    Evarestov R A and Smirnov V P 1983 Phys. Status Solidi B 119 9–40 URL https://onlinelibrary.wiley.com/doi/abs/10.1002/pssb.2221190102

    work page doi:10.1002/pssb.2221190102 1983

  3. [3]

    Monkhorst H J and Pack J D 1976 Phys. Rev. B 13(12) 5188–5192 URL https://link.aps.org/doi/10.1103/PhysRevB.13.5188

    doi:10.1103/physrevb.13.5188 1976

  4. [4]

    Rajagopal G, Needs R J, Kenny S, Foulkes W M C and James A 1994 Phys. Rev. Lett. 73(14) 1959–1962 URL https://link.aps.org/doi/10.1103/PhysRevLett.73.1959

    work page doi:10.1103/physrevlett.73.1959 1994

  5. [5]

    Rajagopal G, Needs R J, James A, Kenny S D and Foulkes W M C 1995 Phys. Rev. B 51(16) 10591–10600 URL https://link.aps.org/doi/10.1103/PhysRevB.51.10591

    work page doi:10.1103/physrevb.51.10591 1995

  6. [6]

    Lin C, Zong F H and Ceperley D M 2001 Phys. Rev. E 64(1) 016702 URL https://link.aps.org/doi/10.1103/PhysRevE.64.016702

    doi:10.1103/physreve.64.016702 2001

  7. [7]

    Cunningham S L 1974 Phys. Rev. B 10(12) 4988–4994 URL https://link.aps.org/doi/10.1103/PhysRevB.10.4988

    work page doi:10.1103/physrevb.10.4988 1974

  8. [8]

    Status Solidi B 80 K89–K93 URL https://onlinelibrary.wiley.com/doi/abs/10.1002/pssb.2220800239

    Bashenov V K, Bardashova M and Mutal A M 1977 Phys. Status Solidi B 80 K89–K93 URL https://onlinelibrary.wiley.com/doi/abs/10.1002/pssb.2220800239

    work page doi:10.1002/pssb.2220800239 1977

  9. [9]

    Stevanovic V 2024 A Python code for calculating the mean-value (Baldereschi’s) point for any crystal structure URL https://arxiv.org/abs/2405.00925

    Pith/arXiv arXiv 2024

  10. [10]

    Fraser L M, Foulkes W M C, Rajagopal G, Needs R J, Kenny S D and Williamson A J 1996 Phys. Rev. B 53(4) 1814–1832 URL https://link.aps.org/doi/10.1103/PhysRevB.53.1814

    doi:10.1103/physrevb.53.1814 1996

  11. [11]

    Chiesa S, Ceperley D M, Martin R M and Holzmann M 2006 Phys. Rev. Lett. 97(7) 076404 URL https://link.aps.org/doi/10.1103/PhysRevLett.97.076404

    doi:10.1103/physrevlett.97.076404 2006

  12. [12]

    Drummond N D, Needs R J, Sorouri A and Foulkes W M C 2008 Phys. Rev. B 78(12) 125106 URL https://link.aps.org/doi/10.1103/PhysRevB.78.125106

    doi:10.1103/physrevb.78.125106 2008

  13. [13]

    Azadi S and Foulkes W M C 2015 J. Chem. Phys. 143 102807 ISSN 0021-9606 URL https://doi.org/10.1063/1.4922619

    work page doi:10.1063/1.4922619 2015

  14. [14]

    Helfrich B 1985 Theor. Comp. Sci. 41 125–139 ISSN 0304-3975 URL https://www.sciencedirect.com/science/article/pii/0304397585900672

    arXiv 1985

  15. [15]

    Chadi D J and Cohen M L 1973 Phys. Rev. B 8(12) 5747–5753 URL https://link.aps.org/doi/10.1103/PhysRevB.8.5747

    work page doi:10.1103/physrevb.8.5747 1973

  16. [16]

    Sobol I 1976 USSR Comp. Math. Math. Phys. 16 236–242 ISSN 0041-5553 URL https://www.sciencedirect.com/science/article/pii/0041555376901543

    arXiv 1976

  17. [17]

    Bratley P and Fox B L 1988 ACM Trans. Math. Softw. 14 88–100 ISSN 0098-3500 URL https://doi.org/10.1145/42288.214372 13 IOP Publishing Journal vv (yyyy) aaaaaa Author et al

    arXiv 1988

  18. [18]

    Needs R J, Towler M D, Drummond N D, López Ríos P and Trail J R 2020 J. Chem. Phys. 152 154106 ISSN 0021-9606 URL https://doi.org/10.1063/1.5144288

    doi:10.1063/1.5144288 2020

  19. [19]

    The casino program is available to download for academic use, worldwide, free of cost

    Needs R J, Towler M D, Drummond N D and López Ríos P 2026 Casino quantum Monte Carlo software. The casino program is available to download for academic use, worldwide, free of cost. For commercial use of the software, please discuss with the authors. Access to casino’s Git repository is available for academic use on request to the authors. The program for...

    2026

  20. [20]

    Chulkov E V and Psakh’e S G 1979 Sov. Phys. J. 22 1237–1238 ISSN 1573-9228 URL https://doi.org/10.1007/BF00894985

    work page doi:10.1007/bf00894985 1979

  21. [21]

    Lynch R W and Drickamer H G 1966 J. Chem. Phys. 44 181–184 ISSN 0021-9606 URL https://doi.org/10.1063/1.1726442

    work page doi:10.1063/1.1726442 1966

  22. [22]

    Kristallogr

    Clark S J, Segall M D, Pickard C J, Hasnip P J, Probert M I J, Refson K and Payne M C 2005 Z. Kristallogr. 220 567–570 URL http://www.degruyter.com/view/j/zkri.2005.220. issue-5-6-2005/zkri.220.5.567.65075/zkri.220.5.567.65075.xml

    2005

  23. [23]

    Trail J R and Needs R J 2005 J. Chem. Phys. 122 014112 URL http://scitation.aip.org/content/aip/journal/jcp/122/1/10.1063/1.1829049

    work page doi:10.1063/1.1829049 2005

  24. [24]

    Trail J R and Needs R J 2005 J. Chem. Phys. 122 174109 URL http://scitation.aip.org/content/aip/journal/jcp/122/17/10.1063/1.1888569

    work page doi:10.1063/1.1888569 2005

  25. [25]

    Drummond N D, Towler M D and Needs R J 2004 Phys. Rev. B 70(23) 235119 URL https://link.aps.org/doi/10.1103/PhysRevB.70.235119

    doi:10.1103/physrevb.70.235119 2004

  26. [26]

    Umrigar C J, Toulouse J, Filippi C, Sorella S and Hennig R G 2007 Phys. Rev. Lett. 98(11) 110201 URL http://link.aps.org/doi/10.1103/PhysRevLett.98.110201 14

    work page doi:10.1103/physrevlett.98.110201 2007