REVIEW 2 major objections 5 minor 55 references
Correlated quantum shift vector of particle-hole excitations
T0 review · 2 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Bound excitons erase the light-polarization dependence of the quantum shift vector.
desk verdict A clean, internally consistent theory of polarization-independent excitonic shift vectors; the proof rests on exponential localization, so the strong experimental claim is a thermodynamic-limit statement rather than a universal rule. 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
The load-bearing object is the flux-threaded Bethe-Salpeter equation in relative coordinates, $H^\kappa_Q(r,r')=e^{-i\kappa\cdot(r-r')}H_Q(r,r')$, whose solutions are the exciton envelope functions $\psi_Q(r)$. For a bound exciton the envelope decays as $e^{-|r|/\xi}$, and the Wannier functions decay as $e^{-|x-R|/\xi_W}$, so the phase ansatz $\tilde\psi^\kappa_Q(r)=e^{-i\kappa\cdot r}\psi_Q(r)$ solves the flux-threaded equation up to a deviation bounded by $C e^{-L/(2\xi_M)}$ with $\xi_M=\max(\xi,\xi_W)$. This exponential control is what makes the Wilson-loop derivative $\delta R=|\nabla_\kappa\arg W_\kappa|<D e^{-L/\xi_M}$ vanish in the thermodynamic limit. The final real-space expression $R_{0\to ex}=\sum_{r,r'}\rho(r,r')[r\delta_{r,r'}+D_{cv}(r'-r)]$, with $\rho$ the exciton reduced density matrix and $D_{cv}$ the Wannier dipole difference, carries the argument into a form computable from GW+BSE data.
What would settle it
Numerically solve the flux-threaded Bethe-Salpeter equation (3) for a model with an exponentially localized exciton envelope and compute $\delta R=|\nabla_\kappa\arg W_\kappa|$ at fixed $\kappa$ for growing $L$: if $\delta R$ decays algebraically rather than exponentially, the central bound Eq. (7) fails. Experimentally, a nonzero polarization-dependent shift photocurrent at the exciton resonance of a noncentrosymmetric non-polar material such as unstrained C3z-symmetric 3R-MoS2 would contradict the vanishing result, provided free particle-hole contributions are excluded.
Extended reading notes
Core claim
On its own terms, the paper establishes that the many-body shift vector for an excitonic transition, defined through flux threading, is exponentially insensitive to the light-matter operator: $|R^{V_1}_{0\to ex}-R^{V_2}_{0\to ex}|<D e^{-L/\xi_M}$ (Eq. 7), so in the thermodynamic limit the shift vector does not depend on light polarization. Because the excitonic shift vector then transforms as a genuine vector under point-group operations (Eq. 9), a vertical transition at $Q=0$ in a noncentrosymmetric but non-polar crystal must have zero shift vector, and hence zero shift photocurrent, regardless of the transition matrix element. The same mechanism makes transitions between two excitonic states intrinsic, $R_{ex_1\to ex_2}=R_{0\to ex_2}-R_{0\to ex_1}$, and yields a compact real-space formula built from the exciton reduced density matrix and Wannier dipole moments. Delocalized free particle-hole excitations do not obey any of this: their shift vectors are finite and strongly polarization dependent, making the shift vector a diagnostic of pair localization. The authors also show the conclusions survive topological band obstructions when hybrid Wannier functions are used along the flux direction.
Load-bearing premise
The central premise is that a bound exciton's relative-coordinate envelope and the underlying Wannier functions decay exponentially, so flux insertion changes the envelope only by a phase with boundary violations suppressed as $e^{-L/(2\xi_M)}$; for weakly bound, very large excitons in finite samples that suppression is slow, and the polarization independence and vanishing results are strict only in the thermodynamic limit.
Editorial extensions
If this is right
- Vertical excitonic transitions in noncentrosymmetric but non-polar crystals have zero shift vector, so the excitonic shift photocurrent vanishes; this explains the zero shift current observed in C3z-symmetric 3R-MoS2 and its activation under strain.
- Excitonic shift vectors in polar crystals stay aligned with the polar axis and lose their light-polarization dependence, even though the absorption matrix element itself can still depend on polarization.
- Shift photocurrent measurements can serve as a geometric ruler that distinguishes bound excitonic pair waves from delocalized free particle-hole scattering states.
- Transitions between excitonic states inherit the same shift vector, making the response additive and intrinsic: $R_{ex_1\to ex_2}=R_{0\to ex_2}-R_{0\to ex_1}$.
- The real-space formula allows excitonic shift currents to be computed from ground-state-to-exciton transition matrix elements and exciton envelope functions without summing over intermediate states.
Reading between the lines
- A direct testable extension is to measure the polarization angle dependence of shift photocurrent at exciton resonances in a strained non-polar material; if the excitonic contribution dominates, the current direction should remain locked to the strain-induced polar axis rather than rotating with the light polarization.
- The same flux-threading argument should apply to any interaction-bound composite excitation with an exponentially localized internal wavefunction, such as trions or biexcitons in two-dimensional semiconductors, predicting similarly polarization-independent shift geometry.
- For weakly bound, large-radius excitons in finite samples, the exponential suppression is slow, so small residual polarization dependence proportional to $e^{-L/\xi}$ should be observable in mesoscopic samples and could be used to extract the exciton radius $\xi$.
- If confirmed, the dichotomy suggests that above-gap photocurrent signals that do depend on polarization are dominated by unbound particle-hole pairs, so excitonic and continuum contributions to shift current can be separated spectrally.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the quantum shift vector of excitonic (bound) particle-hole excitations within a real-space Bethe-Salpeter description. It argues that flux threading in the relative coordinate can be gauged away for exponentially localized exciton envelopes up to boundary corrections of order exp(-L/xi), making the many-body shift vector independent of the light-matter operator in the thermodynamic limit. Consequently the excitonic shift vector transforms as a vector under point-group symmetries, vanishing for vertical transitions in non-polar noncentrosymmetric crystals, in contrast to delocalized free particle-hole excitations. The authors support the scaling bounds with numerical solution of a flux-threaded BSE on a honeycomb lattice and provide a real-space Wannier formula for the excitonic shift vector.
Significance. If the result holds, it identifies a qualitative, non-perturbative effect of electron-hole binding on quantum geometric response and offers a sharp diagnostic of pair localization. It also rationalizes the observed absence of excitonic shift photocurrent in certain non-polar materials. The manuscript's analytical derivation is self-contained, with the exponential suppression arguments and symmetry transformation worked out in the SI, and the numerical simulations directly verify the central scaling (flat Wilson loop and exponentially decaying Thouless number for the exciton). The real-space formula for the exciton shift vector is computationally attractive for GW+BSE implementations. The main caveat is that the proof is carried out within the single-pair BSE approximation and the thermodynamic limit; this does not undermine the central claim within that stated domain but does require the conclusions to be phrased with that scope.
major comments (2)
- [Flux threading and particle-hole excitations, Eq. (5)] The step from the residual bound in Eq. (4) to the eigenstate statement in Eq. (5) omits a spectral-gap condition. The residual norm ||(H^kappa - E) exp(-i kappa.r) psi|| is exponentially small, but the distance between the ansatz and the exact eigenstate is controlled by this residual divided by the gap to the nearest other eigenstate. For a weakly bound exciton this gap is the binding energy, and for degenerate exciton manifolds the scalar-phase ansatz is not justified without an additional degeneracy-lifting argument. The authors should either prove the eigenstate bound with the gap made explicit (e.g., via a Davis-Kahan type argument) or state the non-degenerate, finite-binding-energy conditions under which Eq. (5) holds.
- [Abstract and Discussion] The central conclusion is formulated more strongly than the model assumptions warrant. The proof of Eqs. (7) and (9) is carried out in the real-space BSE single-pair description for exponentially localized Wannier functions and envelope in the thermodynamic limit; it does not cover multi-pair (higher-order BSE) contributions, weakly bound excitons at finite L, or degenerate exciton manifolds. The experimental comparison in Ref. [12] is consistent with the result but does not by itself establish the absence of all excitonic shift-current contributions in realistic materials. The authors should qualify the 'conclusively rules out' language and state the BSE/single-pair and thermodynamic-limit scope explicitly.
minor comments (5)
- [Fig. 2 and text] The text says the Wilson loop behavior is shown in Fig. 2(c), but Fig. 2(b) shows arg(W_kappa) as a function of kappa and Fig. 2(c) shows the standard deviation versus sqrt(N). Please correct the cross-references.
- [Introduction and Flux threading section] There are typos: 'posssess' in the Introduction and 'Imporantly' in the Flux threading section. These should be corrected.
- [SI Eq. (S26)] The notation R^{V12}_{ex1->ex2}, R^{V02}_{0->ex2}, and R^{V01}_{0->ex1} is introduced without defining the superscripts V01, V02, V12. Please define these consistently.
- [Table I] Table I would be more self-contained with a footnote explaining NCS, polar, and noting that there are 21 noncentrosymmetric point groups, of which 10 are polar.
- [Symmetry analysis, Eq. (9)] The notation C3z for the threefold rotation is used without a definition; please state that it denotes the C3 rotation about the z axis.
Circularity Check
No significant circularity: the derivation chain is self-contained; the V-insensitivity bound and symmetry consequences follow from explicit exponential-localization estimates, not from fitted parameters or load-bearing self-citations.
full rationale
The central derivation is self-contained. Equation (5) is not assumed; it is obtained from the deviation estimate in Eq. (4), which is explicitly evaluated in SI Sec. III using the exponential decay of the exciton envelope and the Wannier functions. Equation (7) follows by applying Eq. (5) and the flux-transformation of Wannier functions to the Wilson loop; SI Sec. IV carries out the phase-cancellation explicitly, so the polarization insensitivity is derived rather than built into the definition of the shift vector. Equation (9) is likewise derived in SI Sec. VI from the V-insensitivity already proved, plus standard point-group kinematics; it is not imported from a self-citation. The numerical model uses illustrative parameters (t1, t2, t3, Δ, Vint) and is not fitted to the target claim or to experimental data. Self-citations that appear (e.g., Refs. [11], [28], [48]) are used for background, for naming a Wilson loop, or for contrasting free-particle behavior that is also reproduced numerically in Fig. 2; none is load-bearing. The stated restrictions to a single electron-hole pair described by the Bethe-Salpeter equation and to exponentially localized envelopes/Wannier functions are genuine domain assumptions, and the large-radius-exciton case weakens the finite-size suppression, but these are model limitations rather than circular steps. No fitted parameter is renamed as a prediction, and no uniqueness claim is forced by author self-citation.
Assumptions & free parameters
free parameters (5)
- t1 =
0.8 t0
- t2 =
t0
- t3 =
t0
- Delta =
t0
- Vint =
2.5 t0
assumptions (5)
- domain assumption Exponentially localized Wannier functions exist for the valence and conduction bands.
- domain assumption The Bethe-Salpeter equation with a single electron-hole pair describes the excitations of the insulating ground state.
- domain assumption The ground state is an insulating Slater determinant with filled valence bands and empty conduction bands.
- domain assumption Bound exciton envelope functions decay exponentially in the relative coordinate.
- domain assumption The many-body shift vector definition of Resta 2024 is the correct observable for shift photocurrent.
Cite this review
Pith. "Pith review of Correlated quantum shift vector of particle-hole excitations." pith.science (2026). https://pith.science/paper/XTFBQI6B
@misc{pith2026250707182,
author = {Pith},
title = {Pith review of: Correlated quantum shift vector of particle-hole excitations},
year = {2026},
howpublished = {\url{https://pith.science/paper/XTFBQI6B}},
note = {Machine review of arXiv:2507.07182}
}
read the original abstract
Excitons are a prime example of how electron interactions affect optical response and excitation. We demonstrate that, beyond its spectra, the bound nature of an exciton's electron-hole pair produces a correlated quantum geometry: excitonic excitations possess a quantum shift vector that is independent of light polarization. We find this counterintuitive behavior has dramatic consequences for geometric response: e.g., in noncentrosymmetric but non-polar materials, vertical excitonic transitions possess vanishing shift vector zeroing their shift photocurrent; this contrasts with finite and strongly light polarization dependent shift vectors for non-interacting delocalized particle-hole excitations. This dichotomy makes shift vector a sharp diagnostic of the pair localization properties of particle-hole excitations and demonstrates the non-perturbative effects of electron interactions in excited state quantum geometric response.
Figures
Reference graph
Works this paper leans on
- [12]
- [1]
-
[2]
Ahn, G.-Y
J. Ahn, G.-Y. Guo, and N. Nagaosa, Physical Review X 10, 041041 (2020)
2020
-
[3]
“Quantum Geometric Advantage”. Author Contributions. XY and JCWS conceived and planned the project. XY performed the calculations and numerical simulations with input from JCWS; AS contributed to discussions. XY and JCWS wrote the manuscript with input from AS
- [4]
-
[5]
Onishi and L
Y. Onishi and L. Fu, Physical Review X 14, 011052 (2024)
2024
- [6]
-
[7]
Resta, Physical Review Letters 133, 206903 (2024)
R. Resta, Physical Review Letters 133, 206903 (2024)
work page 2024
Show all 55 references
-
[8]
von Baltz and W
R. von Baltz and W. Kraut, Physical Review B 23, 5590 (1981)
1981
-
[9]
Q. Ma, R. Krishna Kumar, S.-Y. Xu, F. H. Koppens, and J. C. Song, Nature Reviews Physics 5, 170 (2023)
2023
-
[10]
G. H. Wannier, Physical Review 52, 191 (1937)
1937
-
[11]
Wang and X
H. Wang and X. Qian, Science advances 5, eaav9743 (2019)
2019
-
[13]
Dong, M.-M
Y. Dong, M.-M. Yang, M. Yoshii, S. Matsuoka, S. Kita- mura, T. Hasegawa, N. Ogawa, T. Morimoto, T. Ideue, and Y. Iwasa, Nature nanotechnology 18, 36 (2023)
2023
-
[14]
Akamatsu, T
T. Akamatsu, T. Ideue, L. Zhou, Y. Dong, S. Kitamura, M. Yoshii, D. Yang, M. Onga, Y. Nakagawa, K. Watan- abe, et al. , Science 372, 68 (2021)
2021
-
[15]
Krishna Kumar, G
R. Krishna Kumar, G. Li, R. Bertini, S. Chaudhary, K. Nowakowski, J. M. Park, S. Castilla, Z. Zhan, P. A. Pantale´ on, H. Agarwal, et al. , Nature Materials , 1 (2025)
2025
-
[16]
C. Ma, S. Yuan, P. Cheung, K. Watanabe, T. Taniguchi, F. Zhang, and F. Xia, Nature 604, 266 (2022)
2022
-
[17]
R. M. Martin, Electronic structure: basic theory and practical methods (Cambridge university press, 2020)
2020
-
[18]
Marzari and D
N. Marzari and D. Vanderbilt, Physical review B 56, 12847 (1997)
1997
-
[19]
G. D. Mahan, Many-particle physics (Springer Science & Business Media, 2013)
2013
-
[20]
Rohlfing and S
M. Rohlfing and S. G. Louie, Physical Review B 62, 4927 (2000)
2000
-
[21]
Q. Niu, D. J. Thouless, and Y.-S. Wu, Physical Review B 31, 3372 (1985)
1985
-
[22]
Resta, Physical Review Letters 80, 1800 (1998)
R. Resta, Physical Review Letters 80, 1800 (1998)
1998
-
[23]
Souza, T
I. Souza, T. Wilkens, and R. M. Martin, Physical Review B 62, 1666 (2000)
2000
-
[24]
Watanabe, Physical Review B 98, 155137 (2018)
H. Watanabe, Physical Review B 98, 155137 (2018)
2018
-
[25]
Luttinger, Physical Review 84, 814 (1951)
J. Luttinger, Physical Review 84, 814 (1951)
1951
-
[26]
Kohn, Physical review 133, A171 (1964)
W. Kohn, Physical review 133, A171 (1964)
1964
-
[27]
Matsyshyn, G
O. Matsyshyn, G. Vignale, and J. C. Song, arXiv preprint arXiv:2410.21363 (2024)
2024 arXiv
-
[28]
Sipe and A
J. Sipe and A. Shkrebtii, Physical Review B 61, 5337 (2000)
2000
-
[29]
Shi and J
L.-k. Shi and J. C. Song, Physical Review B 100, 201405 (2019)
2019
-
[30]
A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Vander- bilt, and N. Marzari, Computer physics communications 178, 685 (2008)
2008
-
[31]
J. B. Haber, D. Y. Qiu, F. H. da Jornada, and J. B. Neaton, Physical Review B 108, 125118 (2023)
2023
-
[32]
Davenport, J
H. Davenport, J. Knolle, and F. Schindler, Physical Re- view Letters 133, 176601 (2024)
2024
-
[33]
W. J. Jankowski, J. J. Thompson, B. Monserrat, and R.-J. Slager, Nature Communications 16, 1 (2025)
2025
- [34]
-
[35]
It should not be confused with the center-of-mass shift vector described in Ref
R0→ex captures electric polarization changes produced by transitions (e.g., optically generated, scattering in- duced) closely tracking the properties of the relative co- ordinate important for shift responses. It should not be confused with the center-of-mass shift vector des...
-
[36]
While finite momentum transfer Q ̸= 0 can enable to access non-vanishing finite Q shift vectors [48], for far- 6 field light Q are severely constrained by the light cone via energy-momentum conservation suppressing the finite Q shift vector
-
[37]
M. I. Aroyo, ed., International Tables for Crystallogra- phy, Volume A: Space-Group Symmetry (Wiley Online Library, 2016)
2016
-
[38]
Brouder, G
C. Brouder, G. Panati, M. Calandra, C. Mourougane, and N. Marzari, Physical review letters 98, 046402 (2007)
2007
-
[39]
Sgiarovello, M
C. Sgiarovello, M. Peressi, and R. Resta, Physical Review B 64, 115202 (2001)
2001
-
[40]
Marzari, A
N. Marzari, A. A. Mostofi, J. R. Yates, I. Souza, and D. Vanderbilt, Reviews of Modern Physics 84, 1419 (2012)
2012
-
[41]
T. G. Pedersen, Physical Review B 92, 235432 (2015)
2015
-
[42]
Y.-H. Chan, D. Y. Qiu, F. H. da Jornada, and S. G. Louie, Proceedings of the National Academy of Sciences 118, e1906938118 (2021)
2021
-
[43]
Ruan, Y.-H
J. Ruan, Y.-H. Chan, and S. G. Louie, arXiv preprint arXiv:2310.09674 (2023)
2023
-
[44]
Morimoto and N
T. Morimoto and N. Nagaosa, Physical Review B 94, 035117 (2016)
2016
-
[45]
Y. Gao, Y. Zhang, and D. Xiao, Physical Review Letters 124, 077401 (2020)
2020
-
[46]
Edwards and D
J. Edwards and D. Thouless, Journal of Physics C: Solid State Physics 5, 807 (1972)
1972
-
[47]
D. J. Thouless, Physics Reports 13, 93 (1974)
1974
-
[48]
P. A. Lee and T. V. Ramakrishnan, Reviews of modern physics 57, 287 (1985)
1985
-
[49]
L.-k. Shi, D. Zhang, K. Chang, and J. C. Song, Physical Review Letters 126, 197402 (2021)
2021
-
[50]
Taherinejad, K
M. Taherinejad, K. F. Garrity, and D. Vanderbilt, Phys- ical Review B 89, 115102 (2014)
2014
-
[51]
Gresch, G
D. Gresch, G. Autes, O. V. Yazyev, M. Troyer, D. Van- derbilt, B. A. Bernevig, and A. A. Soluyanov, Physical Review B 95, 075146 (2017)
2017
-
[52]
CORRELA TED QUANTUM SHIFT VECTOR OF P AR TICLE-HOLE EXCIT A TIONS
L. Fu and C. L. Kane, Physical Review B—Condensed Matter and Materials Physics 74, 195312 (2006). 7 SUPPLEMENT AR Y INFORMA TION FOR “CORRELA TED QUANTUM SHIFT VECTOR OF P AR TICLE-HOLE EXCIT A TIONS” I. BETHE-SALPETER EQUA TION IN RELA TIVE COORDINA TES It is instructive to e...
2006
-
[53]
The BSE Hamiltonian H(rx, r′ x, p⊥, p′ ⊥) is almost diagonal in rx, r′ x, meaning that the off-diagonal matrix elements decay as e−|rx−r′ x|/ξW
-
[54]
The exciton solution ψQ(rx, p⊥) to the BSE is exponentially localized in rx. 15
-
[55]
Coupling to κx is achieved as H(rx, r′ x, p⊥, p′ ⊥) → e−iκx(rx−r′ x)H(rx, r′ x, p⊥, p′ ⊥). With the above preparation, we can now proceed to put the BSE Hamiltonian on a periodic lattice with linear size L, adopt the ansatz ψQ(rx, p⊥) → e−iκxrx ψQ(rx, p⊥) and arrive at the con...
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.