Kinetic coefficients of two-dimensional electrons with strong Zeeman splitting

A. P. Dmitriev; P. S. Alekseev; Yu. O. Alekseev

REVIEW 4 major objections 4 minor 2 references

Kinetic coefficients of two-dimensional electrons with strong Zeeman splitting

T0 review · 4 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper derives the microscopic relaxation rates that govern hydrodynamic flow in a Zeeman-split two-component 2D electron fluid, showing that only the relative velocity between the two subbands decays and that shear stresses relax indepe

desk verdict Solid strong-splitting rate calculation with two clean structural results, but it does not cover the small-splitting crossover it advertises; worth refereeing after a cleanup. read the letter →

arxiv 2603.03105 v2 pith:JD774C6L submitted 2026-03-03 cond-mat.mes-hall

Yu. O. Alekseev , P. S. Alekseev , A. P. Dmitriev This is my paper · ORCID

classification cond-mat.mes-hall MSC 82C4082C7082D20 PACS 72.10.-d72.20.My

keywords Zeemansplittingtwo-componentelectronfluidrelaxationrateshydrodynamictransportshearviscositymagnetoresistancekineticequationelectron-electroncollisions

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tries to establish the quantitative kinetic coefficients—electron-electron relaxation-rate matrices for the first and second angular harmonics—for a two-dimensional electron gas split by a strong Zeeman field into two subbands. It shows that the first-harmonic matrix has eigenvalues {0, λ(1)}: a state with equal subband velocities does not relax, while the relative velocity relaxes at a specific rate. The second-harmonic matrix is diagonal, so shear stresses in the two subbands relax independently. These rates feed into linear hydrodynamic equations with per-subband shear and Hall viscosities, which the authors argue can quantitatively explain the magnitude of magnetoresistance in recent experiments.

What carries the argument

The central object is the relaxation-rate matrix Γ(m) for the m-th angular harmonic of the two-component distribution function, whose eigenvectors encode which combinations of the two subbands' perturbations are conserved or damped. For m=1 the matrix is a rank-one structure (37) with the null eigenvector corresponding to a Galilean boost (equal velocities) and the other eigenvector to zero total current, relaxing at λ(1) = α(1)11(1+a²). For m=2 the off-diagonal rates vanish identically, a result derived from the integral identity I = ∫ ... = 0 for a ≠ ±1 (49), which reduces to the vanishing of a total derivative of a function G; pairwise cancellation of momentum flux (51) gives the physical

What would settle it

Measure the magnetoresistance of a high-mobility 2D electron gas as a function of in-plane magnetic field from B|| = 0 through the onset of Zeeman splitting; if the amplitude of the positive saturating magnetoresistance does not match the curves predicted by these rates, or if spin-drag experiments detect a nonzero cross-subband second-moment relaxation at moderate splitting, the central claim would be falsified.

Watch

Extended reading notes

Core claim

The central claim is that in the strongly-split, degenerate regime the electron-electron collision integrals produce relaxation matrices with a precise structure: Γ(1) has eigenvalues {0, α(1)11(1+a²)} for a = pF1/pF2, meaning inter-subband momentum exchange dampens only the relative velocity, and Γ(2) is diagonal with zero off-diagonal elements, so the second-moment (shear-stress) perturbations of the two subbands do not entrain each other. The off-diagonal cancellation is traced to a kinematic integral identity (49) and a geometric argument that the momentum-flux change cancels pairwise. Explicit integrals (35)–(52) give the rates, which then enter the hydrodynamic balance equations (56) w

Load-bearing premise

The calculation assumes the two subbands are strongly split and degenerate, T ≪ εF1, εF2 and |εF1 − εF2| ≫ T, and the paper itself notes that the final expressions for the rates cannot be continued to a = 1 (equal Fermi momenta), so the results do not cover the small-splitting regime where the two-component transition actually occurs.

Editorial extensions

If this is right

The derived rates give a parameter-free (up to interaction potential) replacement for the phenomenological relaxation parameters in the two-component viscous fluid model.
The hydrodynamic equations predict that friction between subbands is proportional to relative velocity, so equal-velocity flow is dissipation-free; this is a direct consequence of momentum conservation.
Shear viscosity of each subband is determined by its own second-harmonic relaxation time; no cross-subband shear entrainment, simplifying the transport equations.
If correct, the model can explain the amplitude of the positive saturating magnetoresistance observed in experiments, not just its sign.
The analytic asymptotics for b ~ 1, rs ≪ 1 give scalings α(2)11 ~ (1-b)^{-1}(1-b²)^{-1/2} r_s² ln(1/r_s), λ(1) ~ (1-b²)^{-3/2} r_s² ln(1/r_s), β(2)11 ~ (1-b)^{-2} r_s² ln(1/r_s), which can be tested.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The paper leaves open the small-splitting limit a→1, b→0, where the two-component transition occurs; a separate treatment of near-degenerate subbands would be needed to describe the full negative-to-positive magnetoresistance crossover.
The diagonal structure of Γ(2) suggests a testable prediction: spin-drag or shear-entrainment measurements between the two spin subbands should show no cross-relaxation of the second moment, a feature that could distinguish Zeeman-split systems from valley-split or subband-split systems with different interaction matrix elements.
The kinematic cancellation (49) may generalize to higher even harmonics, implying that the 'no entrainment' property holds for all even moments; odd harmonics beyond the first remain parametrically small.
If the matrix element had a finite exchange contribution between different Zeeman subbands (e.g., in systems with spin-orbit coupling), the diagonal structure of Γ(2) would break down; this could serve as a probe of spin-orbit effects.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The strong-splitting part of this paper is solid and worth having. The two central structural results are genuinely clean: the zero eigenvalue of Γ(1) follows from momentum conservation (the equal-velocity mode does not relax), and the remaining rate λ(1) is identified. The diagonality of Γ(2) — no cross-subband relaxation of the second harmonics — is supported by a total-derivative identity for a≠1, and the physical picture in Fig. 3 is plausible. The calculation is parameter-free: given the Rytova–Keldysh potential and standard material constants, the rates are explicit integrals. That is reproducible and a real contribution.

The soft spots are real but not fatal to that core. The paper cannot cover the regime where the motivating crossover happens. It assumes |εF1−εF2| ≫ T, and the authors themselves state after Eq. (34) that the a→1 limit of the final rate expressions is incorrect, even when finite. Since a→1 exactly as the Zeeman splitting goes to zero, and the experiments show the negative-to-positive magnetoresistance transition while the splitting grows from zero, the advertised quantitative explanation of those experiments is not delivered. The calculation applies to the strongly split regime only. That is a scope gap, not a contradiction in the math.

There are also editing problems: the text after Eq. (50) says α(2)_11 = 0 and α(2)_22 = 0, contradicting Eq. (46) and Fig. 4(c); Eq. (29) carries an extra T; the second line of Eq. (56) uses Δj1 where Δj2 is needed; and the placeholder '[the original text]' survives in the Results section. None of these break the main argument, but they will confuse a referee.

The hydrodynamic equations are assembled rather than derived from the kinetic equation. For a paper whose purpose is to supply kinetic coefficients, that is acceptable if the starting equations are clear; here they are asserted, so a referee should ask for a derivation or a statement that they follow from the standard moment expansion.

The central structural results hold up, so the paper deserves a serious referee. Send it out, but ask the authors to either sharpen the claim about experiments or move it to a stated limitation, and to fix the typos.

Referee Report

4 major / 4 minor

Summary. The paper derives microscopic kinetic coefficients for a two-dimensional electron gas split into two Zeeman subbands by a strong in-plane magnetic field. Starting from the linearized kinetic equation with Rytova-Keldysh electron-electron interactions, the authors compute the relaxation matrices Γ^(1) and Γ^(2) for the first and second angular harmonics of the two-component distribution function. The central structural results are: (i) Γ^(1) has a zero eigenvalue corresponding to equal hydrodynamic velocities of the two subbands (Galilean invariance), with the only nonzero eigenvalue λ^(1)=α^(1)_11(1+a²) governing relaxation of the relative velocity; (ii) Γ^(2) is diagonal, α^(2)_12=α^(2)_21=0, so shear-stress modes in the two subbands relax independently. These rates are used to write two-component Navier-Stokes equations (56) with per-subband shear and Hall viscosities. The paper argues that these equations provide a quantitative microscopic basis for explaining puzzling magnetotransport experiments in tilted magnetic fields.

Significance. If correct, the calculation supplies parameter-free microscopic relaxation rates for a two-component hydrodynamic electron fluid, replacing phenomenological parameters used in prior work (ref. 22). The structural results are supported by elegant arguments: the zero eigenmode follows from momentum conservation, and the vanishing of α^(2)_12 is traced to an exact total-derivative integral identity (Eqs. (49)-(50)) for a≠1. The paper contains no fitted parameters and gives explicit integral expressions for all rates. The significance is therefore high within the strong-splitting regime. However, as detailed below, the claimed connection to the motivating experiments is weakened by the paper's own restriction away from a→1, the regime in which the negative-to-positive magnetoresistance crossover occurs.

major comments (4)

[Eq. (50) and following paragraph] The paper explicitly states that taking the limit a→1 in the final rate expressions is incorrect, even when the a=1 values are finite. Since a=pF1/pF2→1 as the Zeeman splitting goes to zero, the entire small-splitting regime is outside the validity of Eqs. (35)-(55). The singular behavior is visible in Eq. (24): at a=1 and θ2=π the denominator |sin θ2| vanishes and the delta-function argument becomes identically zero, so the factorization used to obtain (26)-(29) is non-uniform. The motivating experiments (refs. 19-21) show the evolution from giant negative to positive saturating magnetoresistance as the splitting grows from zero, i.e., precisely through the a→1, b→0 region. Thus the hydrodynamic equations (56), with the rates computed here, cannot be integrated through the crossover, and the paper's stated goal of quantitatively explaining those experiments is not delivered. The authors
[Physical interpretation, Eq. (51)] The text states 'α^(2)_11=0, α^(2)_22=0' but the context and Eq. (47) clearly show that the off-diagonal elements vanish: α^(2)_12=α^(2)_21=0. This typo is potentially confusing because the diagonal elements α^(2)_11 and α^(2)_22 are computed and displayed in Eq. (46) and in Fig. 4(c). Please correct to α^(2)_12=α^(2)_21=0.
[Eqs. (56)-(57)] The physical argument around Eq. (51) claims that the absence of off-diagonal second-harmonic relaxation 'does not require the a≠1 condition to be fulfilled.' However, the only rigorous proof given is the integral identity (49)-(50), which holds for |a|≠1. At a=1, the derivation of the delta-function form (24) and the subsequent exchange of energy and angular integrations break down. The paper therefore makes two statements that are in tension: one asserting a≠1 is needed for the formulas, and another asserting the conclusion is independent of a. If the authors believe the off-diagonal vanishing is exact for all a, they should provide a careful derivation valid at a=1; otherwise, they should state the result is proven for a≠1 and discuss any continuity assumptions needed for the hydrodynamic equations.
[Conclusion] The hydrodynamic equations are written using the rates derived in the strong-splitting limit. The viscosity coefficients η_xx,i and η_xy,i in Eq. (57) depend on τ_2,i=(α^(2)_ii+β^(2)_ii)^{-1}, which are only computed under the assumption T≪|εF1−εF2|. The equations are therefore not valid in the small-splitting regime where the two-component crossover occurs. This is a scope limitation, not a logical error, but it undermines the concluding claim that solving these equations 'will make it possible to explain the magnitude of the magnetoresistance observed in experiments.' The authors should temper this claim or provide a small-splitting analysis.

minor comments (4)

[Eq. (35)] The manuscript contains several typographical errors and awkward phrasings (e.g., 'the lastleads', 'in the works19–21 the evolution', 'is taken into account'). A careful proofreading is recommended.
[Figure 4] The definition of C in Eq. (36) appears before the integral in Eq. (35) is fully introduced; the notation would be clearer if the constant were defined after the integral expression. Also, check the dimensions: the prefactors in Eqs. (26)-(29) and (35)-(36) should be verified to yield relaxation rates with units of inverse time.
[References] The horizontal axis is described as 'splitting in Fermi energy units' but the parameter b=μB/εF is used later in the text. Please define b explicitly near the figure and specify its sign convention consistently with the positive/negative semi-axis description.
[Eq. (49)] Reference 26 (Rytova) is listed with an unusual format; please provide the full journal reference. Also, some references have inconsistent formatting (e.g., 'Phys. Rev. Lett. B' in ref. 20).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: kinetic coefficients are derived from explicit collision integrals; self-citations are peripheral.

full rationale

The paper's central deliverables—the first-harmonic relaxation matrix eigenstructure (Eqs. 37–43) and the diagonal second-harmonic matrix (Eqs. 46–50)—are derived within the paper from the linearized Boltzmann collision integrals (Eqs. 8–9), the Rytova–Keldysh matrix elements (Eqs. 11–17), and the stated low-temperature, large-splitting kinematics (Eqs. 22–29). No parameter is fitted to experimental data: the material constants (κ=12.9, m=0.067m0, n=2.9×10^11 cm^-2) are fixed inputs, and the plotted relaxation rates are evaluations of the explicit integrals (35)–(52). The zero eigenvalue for common-velocity motion follows from momentum conservation and is verified as an eigenvector identity (38)–(40); the vanishing of the off-diagonal second-harmonic rates is established by the integral identity (49)–(50), not assumed. Self-citations (refs 7, 11, 13–17, 22, 25) are contextual or methodological; none is load-bearing in the sense of importing an unverified uniqueness claim, and the key results are not obtained by renaming a cited result. The paper itself flags a genuine limitation: "in the limit of a→1 in the final expressions for the rates is incorrect, even if at a=1 they turn out to be finite" (paragraph after Eq. 34), and the promised explanation of the small-splitting negative-to-positive magnetoresistance crossover in experiments 19–21 is therefore not fully delivered; but that is a scope gap, not a circular reduction. There is also an evident typo after Eq. (50), where "α(2)_11 = 0, α(2)_22 = 0" should read α(2)_12 = α(2)_21 = 0; this does not alter the derived Eq. (47). No load-bearing step reduces to its own input.

Assumptions & free parameters 0 free parameters · 8 assumptions · 0 invented entities

No free parameters are fitted to data; the relaxation-rate calculation is parameter-free given the Rytova-Keldysh screening model and material constants taken from the target experiment (m = 0.067m0, κ = 12.9, T = 0.3 K, n = 2.9×10^11 cm^-2 used for Fig. 4). No new particles, forces, or conserved quantities are introduced: the two components are the physical spin subbands of the Zeeman-split 2DEG. The central results rest on the domain assumptions listed, especially the strong-splitting condition and the absence of exchange between opposite-spin electrons.

assumptions (8)

domain assumption Degenerate Fermi-liquid limit with large Zeeman splitting: T ≪ εF1,2 and |εF1 − εF2| ≫ T, so |p| ≈ pF and scattering angle and energy transfer are independent.
Invoked in §Relaxation of angular harmonics; the integral reduction (22)-(23) and the δ-function mapping (24)-(25) depend on it. Breaks down at a = 1 as the authors note after Eq. (34).
domain assumption Rytova-Keldysh screened Coulomb potential as the electron-electron interaction, with no self-consistent screening correction.
Eq. (11); the numerical rates (Fig. 4) and the analytic asymptotics (53)-(55) depend on this potential choice. A point-interaction model is used only as a cross-check.
domain assumption Impurity scattering is neglected entirely; only electron-electron collisions are considered.
Stated in §Basic equations: 'impurity density is very low, therefore, the scattering of electrons by impurities is not taken into account.'
domain assumption No spin-orbit interaction and no transitions between Zeeman subbands.
Stated in the Conclusion: the calculation ignores inter-subband transitions, 'which is realized if the spin-orbit interaction is not taken into account.' Spin-orbit coupling in GaAs is weak but nonzero.
domain assumption Opposite-spin electrons are treated as distinguishable: the antisymmetrized intersubband matrix elements are M' = U(P) and M'(Q) = U(Q), with no exchange difference term.
Eqs. (12)-(14). This removes the exchange contribution that exists for intrasubband scattering and is stated in the Conclusion to be the reason the second-moment relaxation matrix is diagonal.
domain assumption Intrasubband collisions are dominated by frontal collisions θ2 ≈ π, reducing the angular δ-function to δ(θ2 − π).
Eq. (30); used to compute the β(m) rates. Standard in the single-component 2D treatment (ref 25), but an approximation in the present two-subband setting.
standard math Standard linearized Boltzmann collision-integral form and angular-harmonic decomposition for a two-component Fermi liquid.
Eqs. (5)-(9); the rate-matrix definition (6) and the harmonic expansion are textbook kinetic theory.
domain assumption The hydrodynamic balance equations (56) are obtained from the kinetic equation by moment expansion with the computed relaxation-rate structure.
The text writes (56) directly ('the hydrodynamic equations for the flow take the form') without showing the moment expansion; the step is standard but asserted rather than derived in the paper.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Kinetic coefficients of two-dimensional electrons with strong Zeeman splitting." pith.science (2026). https://pith.science/paper/JD774C6L

@misc{pith2026260303105,
  author       = {Pith},
  title        = {Pith review of: Kinetic coefficients of two-dimensional electrons with strong Zeeman splitting},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JD774C6L}},
  note         = {Machine review of arXiv:2603.03105}
}

read the original abstract

In modern nanostructures with very low defect densities, has recently been realized a hydrodynamic regime of electric transport, in which two-dimensional (2D) electrons form a viscous fluid due to frequent electron-electron collisions. Many bright transport phenomena have been observed in these systems. Of particular interest are two-component hydrodynamic electron systems, where a richer variety of phenomena becomes possible, than in one-component systems. A simplest way to implement and control a two-component 2D electron system is to place a structure with 2D electrons in a magnetic field with a large component in the 2D plane, that leads to a Zeeman splitting of the electron energy spectrum into two subbands. Here we develop a microscopic model of hydrodynamic transport in such system. By solving the kinetic equation, we calculate the electron-electron relaxation rates of the first and second angular harmonics of the two-component distribution function. Then we derive the hydrodynamic balance equations with the kinetic coefficient containing these rates. Namely, are taken into account the shear viscosity in each fluid component and the effect of the friction between the two components. The last leads to equalization of the hydrodynamic velocities in the two subbands. The obtained equations can be used to explain the results of puzzling magnetotransport experiments in ultra-pure nanostructures in a strong oblique magnetic field.

Figures

Figures reproduced from arXiv: 2603.03105 by the authors.

**Figure 1.** (a): Equilibrium electron distribution of the two components and the corresponding parameters of the two components. (b): Notation for the scattering angles of two electrons. B|| simply as B below. In the described system, due to Zeeman splitting, two subbands with different energies are formed. ε1,2; p = p 2 2m ± µB, µB ≫ kBT . (1) Total particle density n consists of the particle concentrations in these subbands … view at source ↗

**Figure 2.** Kinematics of inter-zone (a) and intra-zone collisions (b). The initial and final pulses correspond to the conservation of the total momentum during collisions and small energy changes. These values are calculated based on the obtained matrix elements, taking into account the assumptions made. Let us first describe the calculation of the contribution of intersubband collisions α (m) . The equation for harmonic ampli… view at source ↗

**Figure 3.** Visual explanations of kinematic the reasons for the [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

Figures from the paper (1 more)

**Figure 4.** Figure 4: Different contributions to the relaxation rates of the [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

2 extracted references · 2 linked inside Pith

[1]

Kinetic coefficients of two-dimensional electrons with strong Zeeman splitting Yu. O. Alekseev, P. S. Alekseev, A. P. Dmitriev Ioffe Institute, Saint Petersburg, Russia In nanostructures with two-dimensional (2D) electrons and very low defect densities, a hydro- dynamic transport regime has recently been realized. In this regime, 2D electrons form a visco...

arXiv 2015
[2020]

S Rytova

26 N. S Rytova. Screened potential of a point charge in a thin film.Moscow University Physics Bulletin, 56:30; see also the on–line version of this article: arXiv:1806.00976, 1967

arXiv 1967

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.

[1] [1]

Kinetic coefficients of two-dimensional electrons with strong Zeeman splitting Yu. O. Alekseev, P. S. Alekseev, A. P. Dmitriev Ioffe Institute, Saint Petersburg, Russia In nanostructures with two-dimensional (2D) electrons and very low defect densities, a hydro- dynamic transport regime has recently been realized. In this regime, 2D electrons form a visco...

arXiv 2015

[2] [2020]

S Rytova

26 N. S Rytova. Screened potential of a point charge in a thin film.Moscow University Physics Bulletin, 56:30; see also the on–line version of this article: arXiv:1806.00976, 1967

arXiv 1967