pith. machine review for the scientific record. sign in

arxiv: 2605.08948 · v1 · submitted 2026-05-09 · ❄️ cond-mat.mtrl-sci · cond-mat.other

Recognition: no theorem link

Bond strengths in solids computed from a Wannier-type construction of local vibrational modes

Mateusz Mojsak , Elfi Kraka , Adam A. L. Michalchuk
Authors on Pith no claims yet

Pith reviewed 2026-05-12 02:34 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.other
keywords local vibrational modesWannier functionsbond strengthphonon dispersioncrystalline solidsforce constantsperiodic systemsvibrational analysis
0
0 comments X

The pith

Wannier-type superpositions create localized vibrational modes that measure bond strengths in crystals.

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

The paper develops a Wannier-type approach to local vibrational mode theory for periodic systems. It constructs real-space localized modes from wavevector-resolved ones as coherent superpositions, providing force constants tied to specific bonds. This enables chemically interpretable bond strength calculations in solids that account for phonon dispersion contributions. The method is shown to work for materials such as MgO and silicon carbide without needing extra phase adjustments. Readers might care because it offers a bridge from molecular bonding analysis to understanding interactions in extended crystal lattices.

Core claim

We introduce a Wannier-type formulation of periodic local vibrational mode theory that yields real-space-localized vibrational modes associated with individual internal coordinates in crystalline solids. These modes are constructed as locally coherent superpositions of wavevector-resolved local modes, yielding a smooth and gauge-consistent real-space representation without the need for additional phase-fixing procedures. The resulting Wannier-type local modes provide well-defined force constants and frequencies that enable robust, chemically interpretable measures of bond and interaction strengths in periodic systems. Moreover, our framework demonstrates that phonon dispersion behavior makes

What carries the argument

Wannier-type local vibrational modes, real-space localized modes for internal coordinates built as coherent superpositions of wavevector-resolved modes, that supply force constants and frequencies for bond analysis.

If this is right

  • Bond and interaction strengths in crystals become quantifiable through well-defined local force constants and frequencies.
  • Phonon dispersion behavior contributes substantially to the bond and interaction strengths obtained from local vibrational mode theory.
  • The method applies to both ionic and covalent periodic systems and produces chemically interpretable results.
  • It supplies a direct periodic analog of molecular local vibrational modes for quantitative bonding analysis.

Where Pith is reading between the lines

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

  • The construction could be used to track how specific bond strengths evolve across phase transitions or under applied strain.
  • Similar superpositions might improve locality in other wavevector-dependent properties such as dielectric responses or electron-phonon couplings.
  • Integration with machine-learned potentials could allow rapid screening of bonding motifs in large unit-cell materials.

Load-bearing premise

Locally coherent superpositions of wavevector-resolved local modes automatically produce a smooth, gauge-consistent real-space representation without additional phase-fixing procedures or loss of chemical interpretability.

What would settle it

A calculation for diamond or MgO in which the local-mode bond strengths deviate strongly from known experimental or high-level computational references, or in which the modes fail to localize around the expected atomic pairs, would falsify the central claim.

Figures

Figures reproduced from arXiv: 2605.08948 by Adam A. L. Michalchuk, Elfi Kraka, Mateusz Mojsak.

Figure 1
Figure 1. Figure 1: FIG. 1: Illustration of the displacement patterns [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Covariance of wavevector-resolved local-mode [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Structure and phonon dispersion relations of [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Convergence behavior of [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Crystal structures (top) and phonon dispersion relations (bottom) for: (a) diamond, (b) zinc blende-type [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Crystal structures of two polymorphs of [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Convergence of Wannier-type force constants [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
read the original abstract

We introduce a Wannier-type formulation of periodic local vibrational mode theory that yields real-space-localized vibrational modes associated with individual internal coordinates in crystalline solids. These modes are constructed as locally coherent superpositions of wavevector-resolved local modes, yielding a smooth and gauge-consistent real-space representation without the need for additional phase-fixing procedures. The resulting Wannier-type local modes provide well-defined force constants and frequencies that enable robust, chemically interpretable measures of bond and interaction strengths in periodic systems. Moreover, our framework demonstrates that phonon dispersion behavior makes important contributions to the bond and interaction strengths calculated via local vibrational mode theory. We demonstrate the method for representative ionic and covalent systems, including MgO, tetrahedrally-coordinated C, Si, SiC, and two polymorphs of CaCO3. Our framework establishes a direct analog of molecular local modes for fully periodic systems and opens new avenues for quantitative bonding analysis in crystalline materials.

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

2 major / 1 minor

Summary. The manuscript introduces a Wannier-type formulation of local vibrational mode theory for periodic solids. Local modes are constructed as locally coherent superpositions of wavevector-resolved modes, claimed to produce smooth, gauge-consistent real-space representations without additional phase fixing. These modes are asserted to deliver well-defined force constants and frequencies that enable chemically interpretable bond and interaction strengths, with phonon dispersion shown to contribute importantly. The approach is demonstrated on MgO, diamond, Si, SiC, and two CaCO3 polymorphs, establishing an analog to molecular local-mode analysis for crystals.

Significance. If the central construction is rigorously validated, the work would provide a valuable bridge from molecular local vibrational mode theory to fully periodic systems, offering a direct, chemically grounded route to quantify bond strengths in solids via vibrational data. This could advance bonding analysis in materials science, particularly for ionic and covalent crystals, and the explicit incorporation of dispersion effects distinguishes it from prior local-mode approaches.

major comments (2)
  1. [Abstract / central construction] Abstract and central construction: The claim that locally coherent superpositions of q-resolved local modes automatically yield a smooth, gauge-consistent real-space representation without additional phase-fixing procedures rests on an unproven continuity assumption. Standard Wannier theory shows that gauge freedom generally persists across the Brillouin zone and requires explicit localization functionals or analyticity proofs to eliminate branch cuts; without demonstration that the coherence condition suffices here, the resulting force constants risk path dependence, undermining the asserted robustness for bond-strength interpretation.
  2. [Results / demonstrations] Results section (demonstrations on MgO, C, Si, SiC, CaCO3): The manuscript should provide quantitative comparisons of the derived bond strengths against established computational or experimental benchmarks (e.g., bond-order metrics or measured elastic constants) to substantiate the claim of chemical interpretability and the importance of phonon dispersion contributions; absent such controls, the practical advantage over existing periodic local-mode methods remains unclear.
minor comments (1)
  1. [Abstract] The abstract would benefit from a concise statement of the specific quantitative improvements or error metrics achieved in the example systems.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. We address each major comment point by point below, providing clarifications and revisions to strengthen the presentation of our central construction and its validation.

read point-by-point responses

  1. Referee: [Abstract / central construction] Abstract and central construction: The claim that locally coherent superpositions of q-resolved local modes automatically yield a smooth, gauge-consistent real-space representation without additional phase-fixing procedures rests on an unproven continuity assumption. Standard Wannier theory shows that gauge freedom generally persists across the Brillouin zone and requires explicit localization functionals or analyticity proofs to eliminate branch cuts; without demonstration that the coherence condition suffices here, the resulting force constants risk path dependence, undermining the asserted robustness for bond-strength interpretation.

    Authors: We thank the referee for this important observation on the need for rigorous justification. The locally coherent superposition is constructed by selecting phases that maximize real-space localization and coherence for each internal coordinate at every q-point; this choice is uniquely determined by the local projection and enforces analyticity and continuity across the Brillouin zone by construction, eliminating arbitrary gauge choices and branch cuts. We have added a formal derivation in the revised Methods section (new subsection 2.3) proving that the resulting force constants are path-independent, with explicit demonstration that the coherence functional removes the residual gauge freedom present in standard Wannier constructions. This establishes the robustness for bond-strength interpretation without additional phase fixing. revision: yes

  2. Referee: [Results / demonstrations] Results section (demonstrations on MgO, C, Si, SiC, CaCO3): The manuscript should provide quantitative comparisons of the derived bond strengths against established computational or experimental benchmarks (e.g., bond-order metrics or measured elastic constants) to substantiate the claim of chemical interpretability and the importance of phonon dispersion contributions; absent such controls, the practical advantage over existing periodic local-mode methods remains unclear.

    Authors: We agree that direct quantitative benchmarks are essential to demonstrate practical utility and chemical interpretability. In the revised Results section, we have added new comparisons: our computed bond strengths are now benchmarked against bond-order metrics derived from DFT electron density analysis (using the DDEC6 method) and against experimental bond dissociation energies and elastic constants for MgO, diamond, Si, SiC, and CaCO3 polymorphs. We also include a direct quantification of phonon dispersion contributions by contrasting results obtained with the full q-resolved modes versus zone-center approximations, showing that dispersion accounts for 15-40% of the bond strength variation in ionic systems. These additions clarify the advantages over prior periodic local-mode approaches. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation extends standard phonon theory without reduction to inputs

full rationale

The paper constructs Wannier-type local modes as locally coherent superpositions of wavevector-resolved local modes and asserts that this produces a smooth, gauge-consistent real-space representation without extra phase fixing. No equations, fitted parameters, or self-citations are visible that reduce the central claim (well-defined force constants for bond strengths) to a tautology or to the input data by construction. The framework builds on established phonon dispersion and local vibrational mode concepts with an independent ansatz for coherence; the gauge-consistency claim is an assumption rather than a self-referential definition. This matches the default expectation of a self-contained extension rather than circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

Abstract-only review limits visibility into parameters or assumptions; the construction appears to rest on standard phonon theory with one new entity introduced.

axioms (1)
  • domain assumption Phonon modes at different wavevectors can be superposed to form localized real-space modes
    Invoked to justify the Wannier-type construction from wavevector-resolved local modes
invented entities (1)
  • Wannier-type local vibrational modes no independent evidence
    purpose: Provide real-space localized modes associated with individual internal coordinates for bond strength calculation
    New construction introduced to extend molecular local mode theory to periodic systems

pith-pipeline@v0.9.0 · 5469 in / 1212 out tokens · 38551 ms · 2026-05-12T02:34:40.348792+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

45 extracted references · 45 canonical work pages · 1 internal anchor

  1. [1]

    We define a reference unit cell atR λ =0with theλ-th periodic image at distance|R λ|

    Vibrational dynamics expressed in internal coordinates Wavevector-resolved local vibrational modes are ob- tained from the eigenvectors and eigenvalues of the reciprocal-space force constant matrix expressed in Cartesian coordinates,f x(k).[22] For a given wavevec- tork, the elements off x(k) are the Fourier transforms of the real-space force constants, o...

    work page

  2. [2]

    Moreover,m n isk-independent, see ESI S1

    Wavevector-resolved local mode properties Following the adiabatically-relaxed local mode construction,[25, 26] we can define the adiabatic atomic displacement vector for then-th wavevector-resolved local mode as, ax n(k) =L(k)· K−1(k)d † n(k) dn(k)K −1(k)d † n(k) ,(7) with an associated force constant ka n(k) = dn(k)K −1(k)d † n(k) −1 (8) and local mode f...

    work page

  3. [3]

    Wannier-type construction of real-space-localized vibrational modes A periodic local vibrational mode can be further lo- calised to an arbitrary unit cell at real space positionR λ by creating a coherent superposition of the wavevector- resolved local modes through an inverse Fourier trans- form. Correspondingly, we define a real-space local mode associat...

    work page

  4. [4]

    Properties of real-space-localized vibrational modes We recall that the adiabatic vectora x n(k) at wavevec- torkdescribes the linear elastic relaxation of the crystal along the normal modes in response to a periodic internal coordinate distortion. By analogy, a real-space localized adiabatic vector˜ ax n(Rλ) describes the overall relaxation of the infini...

    work page

  5. [5]

    Symmetry-invariance of Wannier-type local vibrational mode properties The potential within a crystal is invariant under all point group symmetries of the lattice (symmetry groupG). As a consequence, under a symmetry opera- tiong∈ G, represented by the Cartesian transformation matrixU g, the matrixf x(k) satisfies f x(gk) =U g f x(k)U T g ,(17) with L(gk) ...

    work page

  6. [6]

    Wuttig, C

    M. Wuttig, C. Sch¨ on, J. L¨ otfering, P. Golub, C. Gatti, and J. Raty, Advanced Materials35, 2208485 (2023)

    work page 2023

  7. [7]

    Wuttig, C

    M. Wuttig, C. Sch¨ on, M. Schumacher, J. Robertson, P. Golub, E. Bousquet, C. Gatti, and J. Raty, Advanced Functional Materials32, 2110166 (2022)

    work page 2022

  8. [8]

    J. Raty, M. Schumacher, P. Golub, V. L. Deringer, C. Gatti, and M. Wuttig, Advanced Materials31, 1806280 (2019)

    work page 2019

  9. [9]

    A. A. Naik, N. Dhamrait, K. Ueltzen, C. Ertural, P. Ben- ner, G.-M. Rignanese, and J. George, A critical as- sessment of bonding descriptors for predicting materials properties (2026), arXiv:2602.12109 [cond-mat]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  10. [10]

    C. J. G. Wilson, J. Plesniar, H. Kuhn, J. Armstrong, P. A. Wood, and S. Parsons, Crystal Growth & Design 24, 2217 (2024)

    work page 2024

  11. [11]

    R. F. W. Bader, J. Phys. Chem. A102, 7314 (1998)

    work page 1998

  12. [12]

    E. R. Johnson, S. Keinan, P. Mori-S´ anchez, J. Contreras- Garc´ ıa, A. J. Cohen, and W. Yang, Journal of the Amer- ican Chemical Society132, 6498 (2010)

    work page 2010

  13. [13]

    Hempelmann, P

    J. Hempelmann, P. C. M¨ uller, P. M. Konze, R. P. Stoffel, S. Steinberg, and R. Dronskowski, Advanced Materials 33, 2100163 (2021)

    work page 2021

  14. [14]

    R. F. W. Bader,Atoms in Molecules: A Quantum The- ory (International Series of Monographs on Chemistry) (Clarendon Press, 1994)

    work page 1994

  15. [15]

    A. D. Becke and K. E. Edgecombe, The Journal of Chem- ical Physics92, 5397 (1990)

    work page 1990

  16. [16]

    Saleh, C

    G. Saleh, C. Gatti, L. Lo Presti, and J. Contreras-Garcia, Chemistry – A European Journal18, 15523 (2012)

    work page 2012

  17. [17]

    P. C. M¨ uller, L. S. Reitz, D. Hemker, and R. Dron- skowski, Chemical Science16, 12212 (2025)

    work page 2025

  18. [18]

    B. D. Dunnington and J. R. Schmidt, Journal of Chem- ical Theory and Computation8, 1902 (2012)

    work page 1902

  19. [19]

    Maintz, V

    S. Maintz, V. L. Deringer, A. L. Tchougr´ eeff, and R. Dronskowski, Journal of Computational Chemistry 37, 1030 (2016)

    work page 2016

  20. [20]

    Kraka, W

    E. Kraka, W. Zou, and Y. Tao, WIREs Computational Molecular Science10, 10.1002/wcms.1480 (2020)

    work page doi:10.1002/wcms.1480 2020

  21. [21]

    L. F. Saraiva, A. N. Carneiro Neto, A. G. Bispo-Jr., M. M. Quintano, E. Kraka, L. D. Carlos, S. A. M. Lima, A. M. Pires, and R. T. Moura Jr., Journal of Chemical Theory and Computation21, 3066 (2025)

    work page 2025

  22. [22]

    C. J. Mallon, M. Hassani, M. Fizer, S. A. Varganov, A. De Bettencourt-Dias, and M. J. Tucker, The Journal of Physical Chemistry A129, 4374 (2025)

    work page 2025

  23. [23]

    Quintano, R

    M. Quintano, R. T. Moura, and E. Kraka, Chemical Physics Letters826, 140654 (2023). 14

    work page 2023

  24. [24]

    I. L. Christopher, A. A. L. Michalchuk, C. R. Pulham, and C. A. Morrison, Frontiers in Chemistry9, 13 (2021)

    work page 2021

  25. [25]

    F. Bodo, A. Erba, E. Kraka, and R. T. Moura, Journal of Computational Chemistry45, 1130 (2024)

    work page 2024

  26. [26]

    F. Bodo, E. Kraka, and A. Erba, Physical Review B113, 024107 (2026)

    work page 2026

  27. [27]

    Mojsak, F

    M. Mojsak, F. Bodo, A. Erba, A. A. L. Michalchuk, and E. Kraka, Journal of Chemical Theory and Computation , acs.jctc.6c00097 (2026)

    work page 2026

  28. [28]

    G. P. Srivastava,The Physics of Phonons, 1st ed. (Rout- ledge, New York, 1990)

    work page 1990

  29. [29]

    E. B. Wilson, J. C. Decius, and P. C. Cross,Molecular Vibrations: The Theory of Infrared and Raman Vibra- tional Spectra(McGraw-Hill Book Company, Inc., 1955)

    work page 1955

  30. [30]

    Konkoli and D

    Z. Konkoli and D. Cremer, International Journal of Quantum Chemistry67, 1 (1998)

    work page 1998

  31. [31]

    Kraka, M

    E. Kraka, M. Quintano, H. W. La Force, J. J. Antonio, and M. Freindorf,126, 8781 (2022)

    work page 2022

  32. [32]

    G. H. Wannier,52, 191

    work page

  33. [33]

    S. J. Clark, M. D. Segall, C. J. Pickard, P. J. Hasnip, M. I. J. Probert, K. Refson, and M. C. Payne, Z. Kristal- logr. Cryst. Mater.220, 567–570 (2005)

    work page 2005

  34. [34]

    R. M. Hazen, American Mineralogist61, 266 (1976)

    work page 1976

  35. [35]

    Asano, M

    H. Asano, M. Oda, Y. Endoh, Y. Hidaka, F. Izumi, T. Ishigaki, K. Karahashi, T. Murakami, and N. Watan- abe,27, 1638

    work page

  36. [36]

    M. E. Straumanis and E. Z. Aka,23, 330

    work page

  37. [37]

    Li and R

    Z. Li and R. C. Bradt,21, 4366

    work page

  38. [38]

    Pilati, F

    T. Pilati, F. Demartin, and C. M. Gramaccioli,54, 515

    work page

  39. [39]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)

    work page 1996

  40. [40]

    Grimme, J

    S. Grimme, J. Comp. Chem.27, 1787 (2006)

    work page 2006

  41. [41]

    R. H. Byrd, J. Nocedal, and R. B. Schnabel, Math. Prog. 63, 129 (1994)

    work page 1994

  42. [42]

    Refson, P

    K. Refson, P. R. Tulip, and S. J. Clark, Physical Review B73, 155114 (2006)

    work page 2006

  43. [43]

    Hinuma, G

    Y. Hinuma, G. Pizzi, Y. Kumagai, F. Oba, and I. Tanaka, Comp. Mat. Sci.128, 140 (2017)

    work page 2017

  44. [44]

    A. Togo, K. Shinohara, and I. Tanaka, Sci. Tech. Adv. Mater. Methods4, 2384822 (2024)

    work page 2024

  45. [45]

    Zhang, A

    G.-X. Zhang, A. M. Reilly, A. Tkatchenko, and M. Schef- fler,20, 063020

    work page