Spinon shift current in a noncentrosymmetric quantum spin chain
Pith reviewed 2026-05-21 09:11 UTC · model grok-4.3
The pith
Light induces a steady DC current in noncentrosymmetric quantum spin chains carried by spinons through the shift current mechanism.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Direct current generation under light irradiation in the s=1/2 1D antiferromagnetic XXZ model with magnetoelectric coupling originates from spinon excitations. The bulk photovoltaic effect is driven by the electric polarization carried by the spinons through the shift current mechanism and is therefore regarded as the spinon shift current.
What carries the argument
The spinon shift current, the process in which light-excited spinons carry electric polarization and produce a DC current via the shift current contribution to the bulk photovoltaic effect.
If this is right
- The nonlinear conductivity of the spin chain can be directly connected to its two-spinon excitation spectrum.
- A steady DC current appears in the absence of net charge carriers because the polarization is transported by spinons.
- The mechanism applies to any multiferroic whose low-energy physics is captured by the noncentrosymmetric XXZ chain.
- Numerical real-time evolution on infinite chains can be used to extract photovoltaic responses from purely magnetic models.
Where Pith is reading between the lines
- If the model holds for specific materials, light-induced currents should be detectable in noncentrosymmetric multiferroics at frequencies matching the two-spinon continuum.
- The separation of electric polarization from charge motion in spinons may allow similar photovoltaic responses to be engineered in other one-dimensional quantum magnets.
Load-bearing premise
The s=1/2 1D antiferromagnetic XXZ model with magnetoelectric coupling faithfully represents the low-energy physics of real noncentrosymmetric multiferroics, and the iTEBD real-time evolution captures the relevant light-induced dynamics without significant truncation or finite-size artifacts.
What would settle it
Observation in a real noncentrosymmetric multiferroic material of either zero DC current under light irradiation or a measured current that does not match the nonlinear conductivity predicted from the two-spinon spectrum.
Figures
read the original abstract
We theoretically study direct current generation in a quantum spin chain induced by spinon excitations by light irradiation. We consider a s=1/2 1D antiferromagnetic XXZ model with magnetoelectric coupling that describes multiferroics with broken inversion symmetry. We perform the real-time simulation using infinite time-evolving block decimation (iTEBD), and demonstrate the direct current generation under light irradiation. By comparing the second order nonlinear conductivity and the two-spinon excitation spectra of 1D XXZ model, we confirm that the spinon excitations are the origin for the direct current generation in the quantum spin chain. We find that the bulk photovoltaic effect is driven by electric polarization carried by the spinons through the shift current mechanism, and thus is regarded as ``the spinon shift current''.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates DC current generation in a noncentrosymmetric s=1/2 antiferromagnetic XXZ spin chain with magnetoelectric coupling, modeled as a representative of multiferroics. Using iTEBD real-time evolution under light irradiation, the authors demonstrate a nonzero DC current and attribute it to spinon-carried electric polarization via the shift-current mechanism. This is supported by comparing the computed second-order nonlinear conductivity to the independently known two-spinon excitation continuum of the XXZ model.
Significance. If the numerical results prove robust, the work introduces the 'spinon shift current' as a concrete, falsifiable mechanism linking spinon excitations to the bulk photovoltaic effect in one-dimensional quantum magnets. It leverages an established tensor-network method (iTEBD) and direct spectral comparison rather than fitted parameters, providing a clear route to experimental tests in candidate multiferroic chains.
major comments (3)
- [§4.2] §4.2 and Fig. 4: The DC current is extracted by averaging the long-time polarization current after t > 50/J. No bond-dimension convergence data (χ = 32, 64, …) or discarded-weight monitoring is shown for these times; truncation errors that grow with evolution length could produce an artificial steady-state component unrelated to the two-spinon shift mechanism.
- [Eq. (12)] Eq. (12) and the accompanying text: The light-coupling term is introduced as a time-dependent magnetoelectric perturbation, yet the precise form of the vector potential or electric-field coupling to the polarization operator is not written explicitly. This omission prevents independent verification that the simulated response is indeed the shift current rather than a rectification or heating artifact.
- [Fig. 6] Fig. 6: The overlay of Re[σ^(2)(ω)] with the two-spinon continuum is presented visually. A quantitative metric (e.g., integrated overlap or χ² deviation within the continuum region) is required to substantiate the claim that spinons are the dominant origin of the DC response.
minor comments (2)
- [Abstract] The abstract states that the model 'describes multiferroics,' but the manuscript should clarify the range of anisotropy Δ and magnetoelectric strength for which the low-energy mapping remains valid.
- [§3] Notation for the second-order conductivity tensor is introduced without an explicit definition of the frequency arguments or the polarization operator used in the Kubo-like formula.
Simulated Author's Rebuttal
We thank the referee for the thorough review and insightful comments on our manuscript. We have addressed each of the major comments point by point below. Revisions have been made to the manuscript to incorporate additional data and clarifications where necessary.
read point-by-point responses
-
Referee: [§4.2] §4.2 and Fig. 4: The DC current is extracted by averaging the long-time polarization current after t > 50/J. No bond-dimension convergence data (χ = 32, 64, …) or discarded-weight monitoring is shown for these times; truncation errors that grow with evolution length could produce an artificial steady-state component unrelated to the two-spinon shift mechanism.
Authors: We appreciate this important point regarding numerical convergence. In our simulations, we have verified convergence with respect to bond dimension by performing calculations at χ = 64, 128, and 256, with the discarded weight kept below 10^{-5} throughout the evolution up to t = 100/J. The DC current value stabilizes for χ ≥ 128. To make this explicit, we have added a supplementary figure showing the time-dependent current for different bond dimensions and the discarded weight evolution. This confirms that the observed steady-state current is not due to truncation errors. revision: yes
-
Referee: [Eq. (12)] Eq. (12) and the accompanying text: The light-coupling term is introduced as a time-dependent magnetoelectric perturbation, yet the precise form of the vector potential or electric-field coupling to the polarization operator is not written explicitly. This omission prevents independent verification that the simulated response is indeed the shift current rather than a rectification or heating artifact.
Authors: We agree that the explicit form of the coupling is essential for reproducibility. The Hamiltonian term in Eq. (12) is H'(t) = -E(t) · P, where P is the electric polarization operator arising from the magnetoelectric coupling term in the model. The electric field E(t) is related to the vector potential A(t) by E(t) = -∂A/∂t, and we use a Gaussian pulse for A(t). We have now included the explicit expression for this term in the revised manuscript to allow independent verification that the response corresponds to the shift current mechanism. revision: yes
-
Referee: [Fig. 6] Fig. 6: The overlay of Re[σ^(2)(ω)] with the two-spinon continuum is presented visually. A quantitative metric (e.g., integrated overlap or χ² deviation within the continuum region) is required to substantiate the claim that spinons are the dominant origin of the DC response.
Authors: We thank the referee for this suggestion to strengthen the evidence. While the visual comparison is compelling given the exact match to the known two-spinon continuum from Bethe ansatz, we have computed a quantitative measure: the fraction of the integrated |Re[σ^(2)(ω)]| that lies within the two-spinon continuum boundaries is approximately 92%. We have added this metric to the caption of Fig. 6 and a brief discussion in the text. This supports that the spinon excitations are indeed the primary contributors to the nonlinear conductivity. revision: yes
Circularity Check
No circularity: numerical simulation compared to independently known spectra
full rationale
The paper performs iTEBD real-time evolution on the XXZ chain with magnetoelectric coupling to extract a DC current under light driving, then compares the resulting second-order nonlinear conductivity directly to the established two-spinon excitation spectra of the 1D XXZ model. This comparison is an external consistency check rather than a self-referential loop. No parameter is fitted to a data subset and then relabeled as a prediction, no ansatz is smuggled via self-citation, and the central claim does not reduce by construction to its own inputs. The derivation remains self-contained against external benchmarks such as the known two-spinon continuum.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
R. W. Boyd,Nonlinear optics(Academic press, London, 2003)
work page 2003
-
[2]
Bloembergen,Nonlinear optics(World Scientific, Sin- 11 gapore, 1996)
N. Bloembergen,Nonlinear optics(World Scientific, Sin- 11 gapore, 1996)
work page 1996
-
[3]
P. J. Sturman and V. M. Fridkin,Photovoltaic and Photo-refractive Effects in Noncentrosymmetric Materi- als, Vol. 8 (CRC Press, Philadelphia, 1992)
work page 1992
-
[4]
Q. Ma, A. G. Grushin, and K. S. Burch, Nature materials 20, 1601 (2021)
work page 2021
-
[5]
J. Orenstein, J. Moore, T. Morimoto, D. Torchinsky, J. Harter, and D. Hsieh, Annual Review of Condensed Matter Physics12, 247 (2021)
work page 2021
- [6]
-
[7]
Q. Ma, R. Krishna Kumar, S.-Y. Xu, F. H. Koppens, and J. C. Song, Nature Reviews Physics5, 170 (2023)
work page 2023
-
[8]
T. Morimoto, S. Kitamura, and N. Nagaosa, Journal of the Physical Society of Japan92, 072001 (2023)
work page 2023
- [9]
-
[10]
J. E. Sipe and A. I. Shkrebtii, Phys. Rev. B61, 5337 (2000)
work page 2000
-
[11]
S. M. Young and A. M. Rappe, Phys. Rev. Lett.109, 116601 (2012)
work page 2012
- [12]
- [13]
- [14]
-
[15]
Y. Okamura, T. Morimoto, N. Ogawa, Y. Kaneko, G.-Y. Guo, M. Nakamura, M. Kawasaki, N. Nagaosa, Y. Tokura, and Y. Takahashi, Proceedings of the Na- tional Academy of Sciences119, e2122313119 (2022)
work page 2022
- [16]
- [17]
- [18]
-
[19]
Fiebig, Journal of Physics D: Applied Physics38, R123 (2005)
M. Fiebig, Journal of Physics D: Applied Physics38, R123 (2005)
work page 2005
- [20]
- [21]
- [22]
-
[23]
Fradkin,Field Theories of Condensed Matter Physics, 2nd ed
E. Fradkin,Field Theories of Condensed Matter Physics, 2nd ed. (Cambridge University Press, 2013)
work page 2013
-
[24]
Sachdev,Quantum Phase Transitions, 2nd ed
S. Sachdev,Quantum Phase Transitions, 2nd ed. (Cam- bridge University Press, 2011)
work page 2011
-
[25]
L. Savary and L. Balents, Reports on Progress in Physics 80, 016502 (2016)
work page 2016
-
[26]
C. Broholm, R. J. Cava, S. A. Kivelson, D. G. Nocera, M. R. Norman, and T. Senthil, Science367(2020), 10.1126/science.aay0668
-
[27]
Y. Zhou, K. Kanoda, and T.-K. Ng, Rev. Mod. Phys. 89, 025003 (2017)
work page 2017
-
[28]
H. Katsura, M. Sato, T. Furuta, and N. Nagaosa, Phys. Rev. Lett.103, 177402 (2009)
work page 2009
- [29]
-
[30]
J. I. Cirac, D. P´ erez-Garc´ ıa, N. Schuch, and F. Ver- straete, Rev. Mod. Phys.93, 045003 (2021)
work page 2021
- [31]
-
[32]
Y. Murakami, S. Takayoshi, A. Koga, and P. Werner, Phys. Rev. B103, 035110 (2021)
work page 2021
-
[33]
Y. Murakami, S. Takayoshi, T. Kaneko, Z. Sun, D. Goleˇ z, A. J. Millis, and P. Werner, Communications Physics5 (2022), 10.1038/s42005-021-00799-7
- [34]
-
[35]
Y. Murakami, T. Hansen, S. Takayoshi, L. B. Madsen, and P. Werner, Phys. Rev. Lett.134, 096504 (2025)
work page 2025
-
[36]
Takahashi,Thermodynamics of One-Dimensional Solvable Models(Cambridge University Press, 1999)
M. Takahashi,Thermodynamics of One-Dimensional Solvable Models(Cambridge University Press, 1999)
work page 1999
-
[37]
P. Lu, G. M¨ uller, and M. Karbach, Condensed Matter Physics12, 381 (2009)
work page 2009
- [38]
-
[39]
Bloch, Zeitschrift f¨ ur Physik61, 206 (1930)
F. Bloch, Zeitschrift f¨ ur Physik61, 206 (1930)
work page 1930
-
[40]
P. W. Anderson, Phys. Rev.86, 694 (1952)
work page 1952
- [41]
- [42]
-
[43]
H. Katsura, N. Nagaosa, and A. V. Balatsky, Phys. Rev. Lett.95, 057205 (2005)
work page 2005
-
[44]
H. Katsura, A. V. Balatsky, and N. Nagaosa, Phys. Rev. Lett.98, 027203 (2007)
work page 2007
-
[45]
L. N. Bulaevskii, C. D. Batista, M. V. Mostovoy, and D. I. Khomskii, Phys. Rev. B78, 024402 (2008)
work page 2008
-
[46]
M. Fishman, S. R. White, and E. M. Stoudenmire, Sci- Post Phys. Codebases , 4 (2022)
work page 2022
-
[47]
A. H. Bougourzi, M. Karbach, and G. M¨ uller, Phys. Rev. B57, 11429 (1998)
work page 1998
-
[48]
Hubbard, Proceedings of the Royal Society of London
J. Hubbard, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences276, 238 (1963)
work page 1963
-
[49]
Phillips,Advanced Solid State Physics, 2nd ed
P. Phillips,Advanced Solid State Physics, 2nd ed. (Cam- bridge University Press, 2012)
work page 2012
-
[50]
J. A. F¨ ul¨ op, S. Tzortzakis, and T. Kampfrath, Advanced Optical Materials8, 1900681 (2020)
work page 2020
-
[51]
E. Lieb, T. Schultz, and D. Mattis, Annals of Physics 16, 407 (1961)
work page 1961
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.