REVIEW 2 major objections 2 minor 84 references
Consistent Evaluation of Operators Involving the Position Operator in the Bloch Representation: Application to the Orbital Moment
T0 review · 2 major / 2 minor · reviewed 2026-06-27 · grok-4.3
Pith's one-line read Three rules plus gauge filtration produce consistent, gauge-invariant expressions for position operators in Bloch states.
desk verdict The three rules and gauge filtration give a workable way to fix position-operator inconsistencies in Bloch calculations and reconcile the wave-packet and Wannier orbital-moment expressions, but uniqueness of the scheme is not demonstrated. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Three rules enforcing unit-cell independence, Hermitian conjugacy, and correct intraband velocity, followed by gauge filtration to remove gauge-dependent terms from position-operator expressions.
What would settle it
A tight-binding calculation of the orbital moment on a Bloch state where the wave-packet and Wannier expressions disagree before the rules are applied but agree after gauge filtration is performed.
Extended reading notes
Core claim
We propose three rules for evaluating operators involving the position operator in the Bloch representation. The rules satisfy independence from the choice of unit cell, preservation of Hermitian conjugacy for the product of operators, and recovery of the correct intraband velocity. We introduce a gauge-filtration scheme that systematically removes gauge-dependent contributions. This methodology ensures that the resulting quantities correspond to observable physical phenomena. Application of the framework reconciles the self-rotation of the wave packet with the local circulation of the Wannier function.
Load-bearing premise
The three physical conditions together with the gauge-filtration procedure are sufficient to guarantee that all resulting quantities are observable and gauge-invariant.
Editorial extensions
If this is right
- The orbital moment from wave-packet self-rotation equals the moment from Wannier-function circulation.
- Velocity, orbital moment, and electric polarization become unambiguously defined and observable.
- Products of operators that include the position operator remain Hermitian.
- All derived quantities are independent of the arbitrary choice of unit cell.
Reading between the lines
- The same rules could be applied to electric polarization to resolve analogous representation-dependent discrepancies.
- The method may clarify gauge issues that appear in time-dependent or driven systems.
- Previous numerical mismatches in the literature were likely artifacts of incomplete gauge handling rather than distinct physical mechanisms.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes three rules for evaluating operators involving the position operator in the Bloch representation. The rules are devised to satisfy three physical conditions: unit-cell independence, preservation of Hermitian conjugacy for operator products, and recovery of the correct intraband velocity. A gauge-filtration scheme is introduced to systematically remove gauge-dependent contributions, ensuring the resulting quantities are observable. Application of the framework reconciles the self-rotation of the wave packet with the local circulation of the Wannier function for the orbital moment.
Significance. If the rules and filtration procedure hold and uniquely determine gauge-invariant operators, the work would provide a systematic framework for position-containing operators in solids, standardizing calculations of orbital moments, polarization, and velocity. The explicit reconciliation of two physical pictures is a concrete strength. No machine-checked proofs or reproducible code are mentioned, but the parameter-free character of the rules (once the three conditions are accepted) is a positive feature.
major comments (2)
- [Abstract] Abstract: the central claim that the three conditions together with gauge filtration 'ensure that the quantities obtained correspond to observable physical phenomena' and produce a unique reconciliation rests on the unproven premise that these conditions select a unique operator. The manuscript demonstrates that the proposed rules satisfy the conditions but does not show that no other operator obeying the same three conditions exists or that different filtration prescriptions yield identical observables (e.g., orbital-moment matrix elements). This is load-bearing for the reconciliation result.
- [Gauge-filtration section] Gauge-filtration section (likely §III or §IV): without an explicit operator expression after filtration or a demonstration that the filtered operator is independent of the particular filtration prescription chosen, the assertion that filtration 'systematically removes gauge-dependent contributions' remains insufficiently verified for the orbital-moment application.
minor comments (2)
- The three rules are not stated explicitly in the abstract; listing them concisely would improve immediate readability.
- Notation for the position operator r and its matrix elements in the Bloch basis should be introduced with a dedicated equation early in the text to avoid ambiguity when products of operators are discussed.
Simulated Author's Rebuttal
We are grateful to the referee for the detailed review and valuable feedback on our manuscript. We respond to each major comment in turn.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that the three conditions together with gauge filtration 'ensure that the quantities obtained correspond to observable physical phenomena' and produce a unique reconciliation rests on the unproven premise that these conditions select a unique operator. The manuscript demonstrates that the proposed rules satisfy the conditions but does not show that no other operator obeying the same three conditions exists or that different filtration prescriptions yield identical observables (e.g., orbital-moment matrix elements). This is load-bearing for the reconciliation result.
Authors: We acknowledge the validity of this observation. The manuscript shows that the proposed rules satisfy the three conditions and applies them to reconcile the orbital moment expressions, but does not prove that the conditions uniquely determine the operator or that all possible filtrations give the same result. In the revised manuscript, we will modify the abstract to state that the rules provide a consistent evaluation satisfying the physical conditions, leading to the observed reconciliation, without asserting uniqueness. We will also include a short discussion on the determination of the rules by the conditions. revision: yes
-
Referee: [Gauge-filtration section] Gauge-filtration section (likely §III or §IV): without an explicit operator expression after filtration or a demonstration that the filtered operator is independent of the particular filtration prescription chosen, the assertion that filtration 'systematically removes gauge-dependent contributions' remains insufficiently verified for the orbital-moment application.
Authors: We agree that providing an explicit expression would improve clarity. We will revise the gauge-filtration section to include the explicit form of the filtered operator in the context of the orbital moment and verify that the final observable is independent of the specific filtration procedure chosen. revision: yes
Circularity Check
No significant circularity; framework adds independent content via new rules from stated physical conditions.
full rationale
The paper proposes three explicit rules and a gauge-filtration scheme, each devised to enforce independent physical requirements (unit-cell independence, Hermitian conjugacy of products, correct intraband velocity). These conditions are external to the target observables (orbital moment, wave-packet self-rotation vs. Wannier circulation) rather than tautological with them. No equations reduce a derived quantity to a fitted parameter or prior self-citation by construction; the reconciliation follows from applying the newly specified operator-evaluation procedure. The derivation is therefore self-contained and does not collapse to its inputs.
Assumptions & free parameters
assumptions (1)
- domain assumption Bloch representation applies to periodic crystals and position-operator matrix elements can be defined within it
invented entities (1)
-
gauge filtration
Cite this review
Pith. "Pith review of Consistent Evaluation of Operators Involving the Position Operator in the Bloch Representation: Application to the Orbital Moment." pith.science (2026). https://pith.science/paper/7RPJYQAP
@misc{pith2026260611679,
author = {Pith},
title = {Pith review of: Consistent Evaluation of Operators Involving the Position Operator in the Bloch Representation: Application to the Orbital Moment},
year = {2026},
howpublished = {\url{https://pith.science/paper/7RPJYQAP}},
note = {Machine review of arXiv:2606.11679}
}
read the original abstract
The position operator plays a central role in condensed-matter observables such as velocity, orbital moment, and electric polarization. In solid-state physics, the evaluation of operators incorporating the position operator has not reached a consensus, as observed in the operator-level discrepancy between the local circulation of Wannier functions and the self-rotation of wave packets. Here, to achieve a consistent evaluation of such operators, we propose three rules for evaluating operators involving the position operator in the Bloch representation. The rules are devised to satisfy physical conditions: independence from the choice of unit cell, preservation of Hermitian conjugacy for the product of operators, and recovery of the correct intraband velocity. We further address the gauge dependence of the position operator and introduce a scheme termed gauge filtration, which systematically removes gauge-dependent contributions from the operators containing the position operator. This methodology ensures that the quantities obtained from the operator evaluation correspond to observable physical phenomena. By applying our framework, we reconcile the results concerning the self-rotation of the wave packet and the local circulation of the Wannier function. We expect our proposal to establish a consistent framework for evaluating operators involving the position operator.
Reference graph
Works this paper leans on
-
[1]
First, an operator ˆOis local if its matrix ele- ments in the position basis are supported only at coincident points. More generally, inddi- mensions, it can be represented as⟨r ′| ˆO|r⟩=P n[Qd i=1 On;ni(r)(∂ i r′)ni]δ(d)(r′ −r), wheren=Pd i=1 ni; i.e., as a linear combination of derivatives of the Dirac delta function with coefficient func- tionsO n;ni(r...
-
[2]
All operators denoted by ˆOand ˆQtreated in this study, with the exception of the position operator ˆr, are cell- periodic local operators
Second, an operator ˆOis cell-periodic if it respects the discrete translational symmetry of the crystal lattice,O(r) =O(r+R) for any Bravais lattice vectorR. All operators denoted by ˆOand ˆQtreated in this study, with the exception of the position operator ˆr, are cell- periodic local operators. Typical examples, when rep- resented in the cell-periodic ...
-
[3]
V to obtain the expectation value, which may depend on the gauge choice
Evaluate the operator using the three rules estab- lished in Sec. V to obtain the expectation value, which may depend on the gauge choice
-
[4]
Apply the gauge transformation| ˜ψnk⟩= eiχn(k)|ψnk⟩for a test functionχ n(k) and systematically collect theχ n(k)-dependent terms
-
[5]
This counterterm—constructed from the Berry connec- tion and its derivatives—is defined to ensure it does not alter any gauge-invariant physical contri- butions
Remove the gauge dependence of the evaluated operator by adding a counterterm that exactly cancels the gauge-dependent contributions. This counterterm—constructed from the Berry connec- tion and its derivatives—is defined to ensure it does not alter any gauge-invariant physical contri- butions. B. Application of gauge filtration Building upon the gauge tr...
2026
-
[6]
J. M. Luttinger, Phys. Rev.135, A1505 (1964)
1964
-
[7]
E. N. Adams, Phys. Rev.107, 698 (1957)
1957
-
[8]
E. I. Blount, Phys. Rev.126, 1636 (1962)
1962
Show all 84 references
-
[9]
Karplus and J
R. Karplus and J. M. Luttinger, Phys. Rev.95, 1154 (1954)
1954
-
[10]
Zak, Phys
J. Zak, Phys. Rev.168, 686 (1968)
1968
-
[11]
Zak, Phys
J. Zak, Phys. Rev.177, 1151 (1969)
1969
-
[12]
Hasegawa and T
H. Hasegawa and T. Sasada, Phys. Rev.177, 1392 (1969)
1969
-
[13]
Kohn, Phys
W. Kohn, Phys. Rev.115, 1460 (1959)
1959
-
[14]
Adams and E
E. Adams and E. Blount, J. Phys. Chem. Solids10, 286 (1959)
1959
-
[15]
L. Roth, J. Phys. Chem. Solids23, 433 (1962)
1962
-
[16]
D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Phys. Rev. Lett.49, 405 (1982)
1982
-
[17]
Morimoto, S
T. Morimoto, S. Kitamura, and S. Okumura, Phys. Rev. B104, 075139 (2021)
2021
-
[18]
R. W. Boyd,Nonlinear Optics(Elsevier, 2020)
2020
-
[19]
J. P. Provost and G. Vallee, Commun. Math. Phys.76, 289 (1980)
1980
-
[20]
Anandan and Y
J. Anandan and Y. Aharonov, Phys. Rev. Lett.65, 1697 (1990)
1990
-
[21]
Mera and J
B. Mera and J. Mitscherling, Phys. Rev. B106, 165133 (2022)
2022
-
[22]
Avdoshkin and F
A. Avdoshkin and F. K. Popov, Phys. Rev. B107, 245136 (2023)
2023
-
[23]
Pi´ echon, A
F. Pi´ echon, A. Raoux, J.-N. Fuchs, and G. Montambaux, Phys. Rev. B94, 134423 (2016)
2016
-
[24]
Ozawa and B
T. Ozawa and B. Mera, Phys. Rev. B104, 045103 (2021)
2021
-
[25]
C. Xu, A. Haller, S. Hegde, T. Meng, and T. L. Schmidt, Quantum geometry in the dynamics of band-projected operators (2025), arXiv:2503.11425
2025
-
[26]
Vanderbilt and R
D. Vanderbilt and R. D. King-Smith, Phys. Rev. B48, 4442 (1993)
1993
-
[27]
Thonhauser, Int
T. Thonhauser, Int. J. Mod. Phys. B25, 1429 (2011)
2011
-
[28]
Souza, J
I. Souza, J. ´I˜ niguez, and D. Vanderbilt, Phys. Rev. B69, 085106 (2004)
2004
-
[29]
Resta, D
R. Resta, D. Ceresoli, T. Thonhauser, and D. Vanderbilt, Chem. Phys. Chem.6, 1815 (2005)
2005
-
[30]
Thonhauser, D
T. Thonhauser, D. Ceresoli, D. Vanderbilt, and R. Resta, Phys. Rev. Lett.95, 137205 (2005)
2005
-
[31]
Ceresoli, T
D. Ceresoli, T. Thonhauser, D. Vanderbilt, and R. Resta, Phys. Rev. B74, 024408 (2006)
2006
-
[32]
Ceresoli, U
D. Ceresoli, U. Gerstmann, A. P. Seitsonen, and F. Mauri, Phys. Rev. B81, 060409 (2010)
2010
-
[33]
Aryasetiawan and K
F. Aryasetiawan and K. Karlsson, J. Phys. Chem. Solids 128, 87 (2019)
2019
-
[34]
Bhowal and G
S. Bhowal and G. Vignale, Phys. Rev. B103, 195309 (2021)
2021
-
[35]
A. Pezo, D. Garc´ ıa Ovalle, and A. Manchon, Phys. Rev. B106, 104414 (2022)
2022
-
[36]
G. Go, D. An, H.-W. Lee, and S. K. Kim, Nano Letters 24, 5968 (2024)
2024
-
[37]
An and S
D. An and S. K. Kim, Phys. Rev. B111, 104436 (2025)
2025
-
[38]
Sundaram and Q
G. Sundaram and Q. Niu, Phys. Rev. B59, 14915 (1999)
1999
-
[39]
R. B. Atencia, A. Agarwal, and D. Culcer, Advances in Physics: X9, 2371972 (2024)
2024
-
[40]
Blount, inFormalisms of Band Theory, Solid State Physics, Vol
E. Blount, inFormalisms of Band Theory, Solid State Physics, Vol. 13, edited by F. Seitz and D. Turnbull (Aca- demic Press, 1962) pp. 305–373
1962
-
[41]
M. S. Si and G. P. Zhang, EPL149, 26001 (2025)
2025
-
[42]
Vanderbilt,Berry Phases in Electronic Structure Theory: Electric Polarization, Orbital Magnetization and Topological Insulators(Cambridge University Press, 2018)
D. Vanderbilt,Berry Phases in Electronic Structure Theory: Electric Polarization, Orbital Magnetization and Topological Insulators(Cambridge University Press, 2018)
2018
-
[43]
Resta, Phys
R. Resta, Phys. Rev. Lett.80, 1800 (1998)
1998
-
[44]
Fuchs, Proc
K. Fuchs, Proc. Phys. Soc. London, Sec. A176, 214 (1940)
1940
-
[45]
Bross and H
H. Bross and H. Weilacher, Z. Phys. A255, 325 (1972)
1972
-
[46]
D. Xiao, J. Shi, and Q. Niu, Phys. Rev. Lett.95, 137204 (2005)
2005
-
[47]
Xiao, M.-C
D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys.82, 1959 (2010)
1959
-
[48]
Moessner and J
R. Moessner and J. E. Moore,Topological Phases of Mat- ter(Cambridge University Press, 2021)
2021
-
[49]
M. L. Cohen and S. G. Louie,Fundamentals of condensed matter physics(Cambridge University Press, 2016)
2016
-
[50]
Shankar, Topology and Condensed Matter Physics , 253 (2017)
R. Shankar, Topology and Condensed Matter Physics , 253 (2017)
2017
-
[51]
Cayssol and J
J. Cayssol and J. N. Fuchs, J. Phys. Mater.4, 034007 (2021)
2021
-
[52]
Ref. [35] effectively neglects surface terms by adopting a gauge in which the Bloch states are assumed to be smooth and single-valued across the Brillouin zone (i.e., avoid- ing functional discontinuities). However, as discussed in Ref. [37], it is generally impossible to defi...
-
[53]
Juri´ c, Universe8, 129 (2022)
T. Juri´ c, Universe8, 129 (2022)
2022
-
[54]
B. Q. Song, J. D. H. Smith, and J. Wang, Position oper- ators in terms of converging finite-dimensional matrices: Exploring their interplay with geometry, transport, and gauge theory (2024), arXiv:2403.02519 [quant-ph]
2024
-
[55]
Kontani, T
H. Kontani, T. Tanaka, D. S. Hirashima, K. Yamada, and J. Inoue, Phys. Rev. Lett.102, 016601 (2009)
2009
-
[56]
H. S. Mund, S. Tiwari, J. Sahariya, M. Itou, Y. Sakurai, and B. L. Ahuja, J. Appl. Phys.110, 073914 (2011). 12
2011
-
[57]
S. Ding, A. Ross, D. Go, L. Baldrati, Z. Ren, F. Freimuth, S. Becker, F. Kammerbauer, J. Yang, G. Jakob, Y. Mokrousov, and M. Kl¨ aui, Phys. Rev. Lett.125, 177201 (2020)
2020
-
[58]
W.-Y. He, D. Goldhaber-Gordon, and K. T. Law, Nat. Commun.11, 1650 (2020)
2020
-
[59]
Y.-G. Choi, D. Jo, K.-H. Ko, D. Go, K.-H. Kim, H. G. Park, C. Kim, B.-C. Min, G.-M. Choi, and H.-W. Lee, Nature619, 52 (2023)
2023
-
[60]
J. J. Esteve-Paredes and J. J. Palacios, SciPost Phys. Core6, 002 (2023)
2023
-
[61]
The sign difference be- tween thek-derivatives acting on [O(k)] and [Q(k)] in Eq
When unbounded operators—such as the position operator—are involved in the operator product, the cyclic property of the thermodynamic trace functional may be inapplicable since the operator product does not belong to the trace class. The sign difference be- tween thek-derivati...
-
[62]
J. J. Sakurai and J. Napolitano,Modern quantum me- chanics(Cambridge University Press, 2020)
2020
-
[63]
Chang and Q
M.-C. Chang and Q. Niu, Phys. Rev. B53, 7010 (1996)
1996
-
[64]
Upon taking the thermodynamic limit, the Dirac delta functionδ(k ′ −k) emerges, provided that limN→∞ N δk′k/VB.Z. exists
-
[65]
Because the self-rotation, upon incorporating the missing terms in the Wannier-function limit, remains indepen- dent ofR, applying a normalized summation over the lattice vectors—as performed in the evaluation of the lo- cal circulation—does not alter the result
-
[66]
Iba˜ nez-Azpiroz, F
J. Iba˜ nez-Azpiroz, F. de Juan, and I. Souza, SciPost Phys.12, 070 (2022)
2022
-
[67]
Go, H.-W
D. Go, H.-W. Lee, P. M. Oppeneer, S. Bl¨ ugel, and Y. Mokrousov, Phys. Rev. B109, 174435 (2024)
2024
-
[68]
D. Go, D. Jo, C. Kim, and H.-W. Lee, Phys. Rev. Lett. 121, 086602 (2018)
2018
-
[69]
C. M. de Sterke and J. E. Sipe, Phys. Rev. A38, 5149 (1988)
1988
-
[70]
Y. Feng, J. Han, T. Shen, H. Wu, and Z. Qin, J. Opt. Soc. Am. B41, 1471 (2024)
2024
-
[71]
Shindou, R
R. Shindou, R. Matsumoto, S. Murakami, and J.-i. Ohe, Phys. Rev. B87, 174427 (2013)
2013
-
[72]
Rych ly, P
J. Rych ly, P. Gruszecki, M. Mruczkiewicz, J. W. K los, S. Mamica, and M. Krawczyk, Low Temp. Phys.41, 745 (2015)
2015
-
[73]
Mieszczak and J
S. Mieszczak and J. W. K los, Sci. Rep.12, 11335 (2022)
2022
-
[74]
Friedberg, T
R. Friedberg, T. D. Lee, and Y. Pang, J. Math. Phys. 35, 5600 (1994)
1994
-
[75]
Friedberg, T
R. Friedberg, T. D. Lee, Y. Pang, and H. C. Ren, Phys. Rev. D52, 4053 (1995). 13 Appendix A: Rationale behind the rules This section elucidates the justifications for the evaluation rules established in Sec. V, which were formulated to satisfy the physical conditions presented...
1995
-
[76]
The evaluation is conducted in the Bloch representation, utilizing the Bloch state expressed asψ nk(r) =⟨r|ψ nk⟩= eik·runk(r) =⟨r|e ik·ˆr|unk⟩, wherenis the band index
Local operator versus position operator We first establish that the proposed three rules facilitate a consistent evaluation of a cell-periodic local operator ˆO, while ensuring that the evaluation of the position operator ˆrremains independent of the choice of unit cell. The e...
-
[77]
Applying three rules to evaluate the general result When the position operator is absent from a product chain, no differentiated Dirac delta functions appear in the evaluation, rendering the procedure straightforward, as demonstrated in Eq. (A1). Thus, we focus exclusively on ...
-
[78]
[58]: ⟨W|(ˆr−r W)× ˆv|W⟩= Z dk VB.Z
Missing terms in the self-rotation In the literature incorporating self-rotation [41, 42], the associated formulations originate from the expression derived in Ref. [58]: ⟨W|(ˆr−r W)× ˆv|W⟩= Z dk VB.Z. | ˜W(k)| 2 i ℏ ⟨∂kunk| ×( ˆH(k)−E n(k))|∂kunk⟩ + Z dk VB.Z. | ˜W(k)| 2⟨i∂ku...
-
[79]
We demonstrate that the ∆ N-regularized Wannier states resolve these issues, analogous to the use of ∆ N-regularized Bloch states
Missing terms in the local circulation and the Wannier function regularization In this subsection, we analyze the issues associated with using Wannier functions to evaluate the position operator. We demonstrate that the ∆ N-regularized Wannier states resolve these issues, anal...
-
[80]
lim N→∞ ⟨WN n′R′| ˆO|WN nR⟩=⟨W n′R′| ˆO|WnR⟩
-
[81]
lim N→∞ N −1P R⟨WN n′R| ˆO|WN nR⟩= lim N→∞ N −1P R⟨Wn′R| ˆO|WnR⟩
-
[82]
lim N→∞ ⟨WN n′R′|ˆr|WN nR⟩=⟨W n′R′|ˆr|WnR⟩
-
[83]
Here, we focus on demonstrating that lim N→∞ N −1P ˜R⟨WN n′(R′+ ˜R)|ˆr× ˆr|WN n(R+ ˜R)⟩= 0 by leveraging the ∆ N- regularization alongside the established relations
lim N→∞ N −1P R⟨WN n′R|ˆr|WN nR⟩ ̸= lim N→∞ N −1P R⟨Wn′R|ˆr|WnR⟩, but lim N→∞ N −1P R⟨WN n′R|ˆr|WN nR⟩= limN→∞ N −1P R⟨WN n′(R+ ˜R)|ˆr|WN n(R+ ˜R)⟩. Here, we focus on demonstrating that lim N→∞ N −1P ˜R⟨WN n′(R′+ ˜R)|ˆr× ˆr|WN n(R+ ˜R)⟩= 0 by leveraging the ∆ N- regularization...
-
[84]
We can interpret the gauge-independent self-rotation in Eq
Different gauge-removal schemes The self-rotation and local circulation results become explicitly gauge-dependent once the previously omitted terms are incorporated. We can interpret the gauge-independent self-rotation in Eq. (B3) and local circulation in Eq. (B11) as outcomes...
Reviewed June 27, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.