REVIEW 3 major objections 5 minor 29 references
Relaxation operator for quasiparticles in a solid
T0 review · 3 major / 5 minor · reviewed 2026-08-27 · deepseek-v4-flash
Pith's one-line read A relaxation operator pairing each coherence with its time-reversed partner and weighted by Bloch overlaps preserves charge conservation for arbitrary dispersion, replaces the standard low-frequency model, and yields correct static and…
desk verdict A genuinely useful generalization of the charge-conserving relaxation operator to Bloch states, but the claimed universality rests on an energy-only relaxation-rate assumption that is stated rather than justified. 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 time-reversal-paired relaxation operator itself: $R_{ckdq}=-\gamma_{ckdq}(\rho_{ckdq}-G_{ckdq}\rho_{d-q,c-k})$, where $G_{ckdq}$ is the Bloch-overlap ratio defined by Eq. (36). The operator is engineered so that when it is inserted into the continuity-equation condition, Eq. (25), the paired terms cancel exactly after lattice averaging; $G_{ckdq}$ is precisely the coefficient that makes this cancellation work. This coefficient carries the band-eigenvector information, so the operator adapts to any dispersion rather than assuming plane waves or parabolic bands. It also produces a damped oscillator denominator in the susceptibility, which is the structural reason the new response satisfies the fundamental relation between conductivity and polarizability.
What would settle it
Measure the low-frequency linear response of high-mobility monolayer graphene at finite wavevector $\kappa$ and $\omega\ll\gamma$: the new model predicts a real, relaxation-independent susceptibility with vanishing imaginary part, whereas the standard model Eq. (2) produces a spurious relaxation-dependent imaginary part. Seeing the spurious part, or finding a microscopic scattering calculation in which $\gamma_{ckdq}\neq\gamma_{d-q,c-k}$ breaks the cancellation behind Eq. (15), would settle the central claim.
Extended reading notes
Core claim
The central claim is that coherence relaxation in a solid should be written as $R_{ckdq}=-\gamma_{ckdq}(\rho_{ckdq}-G_{ckdq}\rho_{d-q,c-k})$, with $G_{ckdq}=(u^*_{c,-k}\cdot u_{d,-q})/(u^*_{d,q}\cdot u_{c,k})$ and $\gamma_{ckdq}=\gamma_{d-q,c-k}$, rather than as the diagonal relaxation of Eq. (2). This antisymmetric-in-quasimomentum structure makes the off-diagonal contribution to the particle density vanish, so the continuity equation is preserved for arbitrary dispersion and both intraband and interband transitions are included. Solving the coupled equations for the two paired coherences gives a generalized Lindhard susceptibility, Eq. (59), whose denominator is the damped-oscillator form $(E_{ck}-E_{dq})^2-\hbar^2\omega^2-2i\hbar^2\omega\gamma_{ckdq}$ instead of the shifted pole of the standard model. In the static finite-$\kappa$ limit this susceptibility is real and independent of relaxation constants, describing a current-free equilibrium; in the $\kappa\to 0$ limit it yields the interband response plus a Drude intraband term, reproducing the Ohmic DC conductivity.
Load-bearing premise
The load-bearing premise is that the coherence relaxation rate depends only on transition energy, so $\gamma_{ckdq}=\gamma_{d-q,c-k}$; real scattering rates generally depend on band indices and momenta, and if this symmetry fails the cancellation that preserves the continuity equation and the final response formulas would need correction.
Editorial extensions
If this is right
- Response functions computed with the new operator satisfy $\sigma(\omega,\kappa)=-i\omega\chi(\omega,\kappa)$ exactly, removing the order-$\gamma/\omega$ ambiguity between independently calculated conductivity and polarizability.
- In the static limit at finite $\kappa$, the susceptibility is real and independent of the relaxation rate, describing the equilibrium current-free state of a closed system; in the uniform-field limit $\kappa\to 0$, the same expression gives interband response plus a Drude intraband term with the Ohmic DC conductivity.
- The generalized dissipative Lindhard formula, Eq. (59), replaces the substitution $\omega\to\omega+i\gamma$ for screening and plasmon problems in two-dimensional layers, including graphene.
- For monolayer graphene the low-frequency ($\omega\lesssim\gamma$) susceptibility from the new model differs sharply from the standard model; the standard model is invalid in this regime, while the new model gives finite behavior consistent with the static and Drude limits.
- Because the construction only needs the Bloch overlap ratio $G_{ckdq}$, it extends to arbitrary band dispersion, including the massless Dirac spectrum of graphene and the broken-time-reversal Weyl semimetal model treated in the paper.
Reading between the lines
- The paper does not spell this out, but because $G_{ckdq}$ is built purely from eigenvector overlaps, the operator can be evaluated directly from any tight-binding or first-principles band calculation, giving a parameter-free form for coherence relaxation once an energy-dependent rate is supplied.
- A natural testable extension is to compare low-frequency absorption or screening of high-quality graphene at $\omega\sim\gamma$ with Eq. (59): the standard model predicts an extra relaxation-dependent contribution at finite wavevector that the new model forbids.
- The limit-ordering subtlety suggests that for confined or finite samples the crossover from the equilibrium response to Ohmic conduction is controlled by boundary conditions rather than by the relaxation model alone, which links this construction to mesoscopic transport calculations in small graphene devices.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a phenomenological relaxation operator for the single-particle density matrix of quasiparticles in a solid, of the form R_ckdq = -gamma_ckdq (rho_ckdq - G_ckdq rho_{d-q c-k}), with G_ckdq fixed by Bloch-wave overlaps so that the continuity equation is preserved. The authors claim this operator is universal for arbitrary quasiparticle dispersion, preserves the continuity equation, reproduces the correct static current-free limit for finite wave vector and the Ohmic DC limit in a uniform field, and includes both intra- and interband transitions. They apply the construction to graphene and to a TRS-breaking Weyl model, derive a generalized Lindhard-type susceptibility in Eq. (59), and compare with the standard relaxation operator Eq. (2), finding large differences at frequencies comparable to or below the relaxation rate.
Significance. If the construction is correct, it would provide a relatively simple phenomenological fix to a known deficiency of the standard relaxation operator, which violates the continuity equation and the relation sigma = -i omega chi at low frequencies. The paper is explicit and self-contained, and it delivers concrete analytic formulas for graphene, including a comparison with the standard model. The main strengths are the transparent derivation of the operator from the continuity requirement and the demonstration that the static and DC limits are recovered depending on the order of limits. However, the claimed universality is narrower than the abstract suggests, because the construction relies on a symmetry assumption for the relaxation rate, and the paper leaves open the positivity question for the proposed non-Lindblad operator. The factor-of-two discrepancy in the Drude limit also needs clarification. With these points addressed, the work would be a useful contribution to the phenomenological master-equation literature.
major comments (3)
- [§II.B, Eq. (27)] The continuity-preserving property of the relaxation operator depends critically on the assumption gamma_ckdq = gamma_{d-q c-k}, stated before Eq. (22) as the condition that the relaxation rate depends only on the transition energy. In the rearrangement leading to Eq. (27), this symmetry is needed to factor a common gamma_ckdq out of the two terms; without it, the cancellation that fixes G via Eq. (28) no longer occurs. For realistic momentum- and band-dependent scattering rates (e.g., from electron-phonon or impurity scattering), this condition fails, so the operator does not preserve the continuity equation in general. The abstract and conclusions describe the operator as 'universal' and 'valid for arbitrary dispersion,' which is too strong; the universality claim should be explicitly qualified to hold only for relaxation rates with the required reversal symmetry.
- [§III.C, Eq. (60) and following] The intraband contribution to the susceptibility in the uniform-field limit, second term of Eq. (60), gives a Drude-type denominator with 2i gamma, whereas the standard Drude model and the standard relaxation operator Eq. (2) give i gamma. The authors acknowledge this 'up to a factor of 2 in the definition of a relaxation constant' but do not reconcile it. Since gamma is a free parameter in both phenomenological models, the factor can be absorbed by redefinition, but then the same gamma also controls the interband transitions in Eq. (59), so the numerical comparison in Section IV and the physical interpretation of gamma are affected. The paper should either fix the convention for gamma, discuss the relation to the standard Drude scattering time, or state clearly that gamma has a different meaning in the two models and adjust the comparison accordingly.
- [§I.C and §II.B] The paper does not address whether the proposed relaxation operator Eq. (26) preserves the positive semidefiniteness of the density matrix. The discussion in §I.C concerning Eq. (5) is for a real-basis operator and relies on specific arguments from Refs. [6,19]; the new operator contains generally complex coefficients G_ckdq and couples density-matrix elements in a different way, so those arguments do not directly transfer. Since Eq. (1) is a master equation and the proposed operator is not of Lindblad form, the authors should at least discuss the domain of validity or provide a positivity check, or else explicitly state that positivity is an open issue for this phenomenological model.
minor comments (5)
- [§II.B, before Eq. (22)] The assumption that the relaxation rate 'depends only on its energy' is stronger than the condition actually used, gamma_ckdq = gamma_{d-q c-k}. Please state the precise symmetry condition, since the cancellation in Eq. (27) uses only the reversal symmetry and not the full energy-only dependence.
- [§III.C and Appendix C/D] The text in §III.C says 'it will be shown in Appendix C' regarding the independence of the uniform-field limit from |G|, but the relevant calculation appears in Appendix D. Please correct the cross-reference.
- [Abstract and Conclusions] The phrase 'universal' in the abstract and 'generalize the Lindhard formula' in Section V overstate what is demonstrated. The continuity preservation and the sigma = -i omega chi relation are enforced by the construction of G, so they are consistency checks rather than independent predictions. I recommend rewording these passages to emphasize that the operator is constructed to satisfy the required properties, not that it is derived from first principles or that the properties are empirically verified.
- [Eq. (65)] In the zero-frequency limit, the expression contains delta((E_c,q+kappa - E_dq)^2) and is then converted to (1/2) delta(E_c,q+kappa - E_dq). The factor of 1/2 is correct, but it would be helpful to include the identity explicitly for readers.
- [Throughout] There are a few typographical errors, e.g., 'describe the the procedure' before Appendix B and an unclosed parenthesis in the sentence following Eq. (60). A careful proofreading pass is recommended.
Circularity Check
No significant circularity: the operator is constructed from the continuity constraint, and the graphene response formulas are algebraically derived, not fitted.
full rationale
No step in the paper reduces a claimed prediction to an input by construction or to a self-citation chain. The core object is an explicit ansatz R_ckdq=-gamma_ckdq(rho_ckdq-G_ckdq rho_d-q c-k) in Eq. (26); the matrix G is then solved from the requirement that the continuity sum in Eq. (25) vanish, yielding Eq. (36). Thus preserving the continuity equation is a design constraint used to determine G, not a prediction that G later verifies. The subsequent results (sigma=-i omega chi, the static limit, and the Drude limit) are consistency checks of that construction; the paper does not claim them as independent empirical predictions. The graphene susceptibility Eqs. (63), (67), (68) and the comparison with the standard model are derived algebraically from the stated model, with no parameter fitted to the compared quantities. The assumption gamma_ckdq=gamma_d-q c-k is stated explicitly before Eq. (22) as an input needed for the cancellation; it limits applicability to energy-only relaxation rates, but this is a stated limitation rather than a circular reduction, since the final formulas do not presuppose the conclusion. Self-citations [6,7,9] motivate the antisymmetric form and document known failures of Eq. (2), but the derivation of G and of the susceptibility is self-contained and does not rest on any unverified uniqueness theorem from those works.
Assumptions & free parameters
free parameters (1)
- Coherence relaxation rate gamma_ckdq =
Not fitted; set to 10^14 s^-1 for graphene plots
assumptions (7)
- domain assumption A quantum master equation with a phenomenological relaxation operator R(rho) is a valid starting point for open condensed-matter systems.
- domain assumption Time-reversal symmetry of the isolated system, or an averaged formulation when TRS is broken, allows pairing transitions with time-reversed partners.
- ad hoc to paper The relaxation rate for a given transition depends only on its energy, so gamma_ckdq = gamma_{d-q c-k}.
- domain assumption Diagonal elements of the density matrix are not perturbed in linear response, so only off-diagonal relaxation matters.
- domain assumption The non-Lindblad relaxation operator preserves positivity of the density matrix at times exceeding the Markov averaging time.
- domain assumption For the averaged description, the quasimomentum interval satisfies delta k times lattice constant much less than 1, and the Hamiltonian can be truncated near an extremum.
- ad hoc to paper The condition |G_ckdq|^2 = 1 holds for the systems treated, including the TRS-breaking Weyl model in Appendix C.
Cite this review
Pith. "Pith review of Relaxation operator for quasiparticles in a solid." pith.science (2026). https://pith.science/paper/STVHOY4K
@misc{pith2026200909021,
author = {Pith},
title = {Pith review of: Relaxation operator for quasiparticles in a solid},
year = {2026},
howpublished = {\url{https://pith.science/paper/STVHOY4K}},
note = {Machine review of arXiv:2009.09021}
}
read the original abstract
Popular models of the phenomenological relaxation operators that are widely used in the master equation formalism for open condensed-matter systems have significant flaws ranging from limited applicability to violation of fundamental physical principles. We propose a relatively simple universal model of the relaxation operator which is free from these flaws, has a correct static limit, correct direct-current limit in a uniform electric field, includes both interband and intraband transitions, and is valid for an arbitrary dispersion of quasiparticles in a solid. We use the proposed operator to generalize the Lindhard formula and derive explicit expressions for the relaxation operator for Dirac materials with an unconventional energy spectrum of quasiparticles, such as graphene and Weyl semimetals. We compare the linear susceptibility spectra for graphene obtained with different relaxation models and show that the proposed relaxation operator leads to physically meaningful behavior of the susceptibility at low frequencies whereas the existing models become completely invalid.
Figures
Reference graph
Works this paper leans on
-
[1]
H. P. Breuer and F. Petruccione, The theory of open quantum systems (Oxford University Press, Great Clarendon Street, 2002)
2002
-
[2]
Blum, Density matrix theory and applications; 3rd ed
K. Blum, Density matrix theory and applications; 3rd ed. , Springer Series on Atomic Optical and Plasma Physics (Springer, Berlin, 2012)
work page 2012
-
[3]
Y. Wang, M. Tokman, and A. Belyanin, Phys. Rev. A 91, 033821 (2015)
work page 2015
- [4]
-
[5]
N. D. Mermin, Phys. Rev. B 1, 2362 (1970). 35
work page 1970
- [6]
- [7]
-
[8]
J. Salmilehto, P. Solinas, and M. M¨ ott¨ onen, Phys. Rev. A85, 032110 (2012)
work page 2012
Show all 29 references
-
[9]
M. D. Tokman, Phys. Rev. A 79, 053415 (2009)
2009
-
[10]
G. S. Atwal and N. W. Ashcroft, Phys. Rev. B 65, 115109 (2002)
2002
-
[11]
Haug and S
H. Haug and S. W. Koch, Quantum Theory of the Optical and Elec- tronic Properties of Semiconductors , 5th ed. (WORLD SCIENTIFIC, 2009) https://www.worldscientific.com/doi/pdf/10.1142/7184
2009 doi
-
[12]
Fain and Y
V. Fain and Y. Khanin, Quantum Electronics: Basic Theory , International series of mono- graphs in natural philosophy (MIT Press (MA), 1969)
1969
-
[13]
Lindblad, Communications in Mathematical Physics 40, 147 (1975)
G. Lindblad, Communications in Mathematical Physics 40, 147 (1975)
1975
-
[14]
Scully and M
M. Scully and M. Zubairy, Quantum Optics (Cambridge University Press, 1997)
1997
-
[15]
L. D. Landau and E. M. Lifshitz, Phys. Z. Sowjetunion 8, 153 (1935)
1935
-
[16]
L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory , Teoretich- eskaia fizika (Elsevier Science, 2013)
2013
-
[17]
Fleming, N
C. Fleming, N. I. Cummings, C. Anastopoulos, and B. L. Hu, Journal of Physics A: Mathe- matical and Theoretical 43, 405304 (2010)
2010
-
[18]
Stokes, A
A. Stokes, A. Kurcz, T. P. Spiller, and A. Beige, Phys. Rev. A 85, 053805 (2012)
2012
-
[19]
W. J. Munro and C. W. Gardiner, Phys. Rev. A 53, 2633 (1996)
1996
-
[20]
Gantmakher and I
V. Gantmakher and I. Levinson, Carrier Scattering in Metals and Semiconductors , Modern problems in condensed matter sciences (North-Holland, 1987)
1987
-
[21]
Berestetskii, L
V. Berestetskii, L. Pitaevskii, and E. Lifshitz, Quantum Electrodynamics , v. 4 (Elsevier Science, 1982)
1982
-
[22]
Y. Wang, M. Tokman, and A. Belyanin, Phys. Rev. B 94, 195442 (2016)
2016
-
[23]
A. R. Kutayiah, M. Tokman, Y. Wang, and A. Belyanin, Phys. Rev. B 98, 115410 (2018)
2018
-
[24]
Vafek and A
O. Vafek and A. Vishwanath, Annual Review of Condensed Matter Physics 5, 83 (2014)
2014
-
[25]
Orlita, D
M. Orlita, D. Basko, M. Zholudev, F. Teppe, W. Knap, V. Gavrilenko, N. Mikhailov, S. Dvoretskii, P. Neugebauer, C. Faugeras, et al. , Nature Physics 10, 233 (2014)
2014
-
[26]
Katsnelson, Graphene: Carbon in Two Dimensions (Cambridge University Press, 2012)
M. Katsnelson, Graphene: Carbon in Two Dimensions (Cambridge University Press, 2012)
2012
-
[27]
X. Yao, M. Tokman, and A. Belyanin, Phys. Rev. Lett. 112, 055501 (2014). 36
2014
-
[28]
Q. Chen, A. R. Kutayiah, I. Oladyshkin, M. Tokman, and A. Belyanin, Phys. Rev. B 99, 075137 (2019)
2019
-
[29]
Q. Chen, M. Erukhimova, M. Tokman, and A. Belyanin, Phys. Rev. B 100, 235451 (2019). 37
2019
Reviewed August 27, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.