Pith. sign in

REVIEW 3 major objections 5 minor 27 references

Non-linear Faraday Precession of Light Polarization in Time-Reversal Invariant Materials

T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read Light polarization rocks back and forth in nonlinear Hall metals

desk verdict Genuinely new prediction and clean numerics, but the analytic derivation of the central formula has a real second-order error that must be fixed before this is reliable. read the letter →

arxiv 2412.03656 v1 pith:DU7EPHQQ submitted 2024-12-04 cond-mat.mes-hall physics.optics

classification cond-mat.mes-hallphysics.optics PACS 78.20.Fm
keywords nonlinearHalleffectBerrycurvaturedipoleFaradayrotationpolarizationprecessionplasmaoscillationsBoltzmannequationlayeredmetalstime-reversalinvariant
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 predicts that light traveling through a time-reversal-invariant layered metal with a Berry curvature dipole will not simply rotate its polarization but will rock it back and forth as the beam passes through the material. In the weakly nonlinear regime the authors argue is relevant for experiments, the polarization angle follows $\theta(z) \approx \theta_0 \cos(\Omega z/\lambda)$, with $\Omega = (\sqrt{7}/24)(I/I_0)$ growing linearly with light intensity. The degree of polarization oscillates in step with this rocking, and the effect needs no DC current, so it can be seen in contactless optical experiments. The authors show that their analytical approximation agrees quantitatively with numerical integration of the full ODEs.

What carries the argument

The key objects are the dimensionless coupled ODEs, Eqs. (6)-(8), obtained by reducing the Maxwell-Boltzmann equations under a traveling-wave ansatz with phase velocity $v_\phi > c$. These equations describe a fictitious particle whose coordinates are the two electric-field components and an average Berry curvature, and they contain no dimensionless parameters, so all dynamics is set by initial conditions. For the weakly nonlinear regime, the paper uses the energy, angular momentum, and Laplace-Runge-Lenz vector of the two-dimensional harmonic oscillator as slowly varying constants of motion, and shows that their long-time averages close into equations whose solution gives the rocking angle and the ellipticity. In the opposite, strongly nonlinear limit the same ODEs map exactly onto a pendulum equation, which supplies the physical picture of the polarization rocking.

What would settle it

Measure the polarization rotation as a function of sample thickness in a layered material with a known Berry curvature dipole (for example WTe2 or MoTe2) at fixed intensity below the damage threshold: the claim predicts an oscillatory $\theta(z)$ with a period that grows linearly with intensity. Observing a monotonic rotation, no thickness oscillation, or a period that does not scale linearly with intensity would falsify the central claim.

Watch

Extended reading notes

Core claim

The central claim is that the nonlinear Hall effect, acting dynamically rather than as a predetermined external field, transforms the propagation of a plane wave in a layered metal into a pendulum-like rocking of the polarization state. Starting from the coupled Maxwell-Boltzmann equations for a single-band metal with a Berry curvature dipole, the authors reduce the problem to three dimensionless ODEs with no free parameters, and in the weakly nonlinear regime show that the electric field traces an ellipse whose semi-major axis direction $\theta(\xi)$ obeys $\theta = \theta_0 \cos(\Omega\xi)$ while the ellipticity oscillates as $b = (a_0 \theta_0/\sqrt{7}) \sin(\Omega\xi)$, with $\Omega = a_0^2\sqrt{7}/24$. This is the non-linear Faraday precession: a rocking, not a steady rotation, with a simultaneously oscillating degree of polarization and a light intensity that stays constant because the Berry-dipole response is nondissipative. The paper establishes this through a second-order perturbative expansion around the harmonic-oscillator constants of motion and verifies it against numerical solutions.

Load-bearing premise

The derivation requires the nonlinear Hall current to be a local, instantaneous, collisionless response of a single band; if electron scattering, interband transitions, or wavevector corrections are significant in real samples, the predicted rocking will be damped or altered.

Editorial extensions

If this is right

  • A thickness-dependent Faraday rotation measurement in a material like WTe2, MoTe2, or TaIrTe4 should reveal the polarization angle oscillating with sample thickness rather than rotating monotonically.
  • The rotation per traveled wavelength grows linearly with light intensity, giving a clean experimental knob to distinguish this effect from ordinary Faraday rotation.
  • Because the Berry-dipole contribution is nondissipative, the light intensity stays constant while the polarization state rocks, which distinguishes the effect experimentally from absorption-driven changes.
  • The time dependence of the polarization oscillations implies emission of lower-frequency radiation at a characteristic frequency $\omega_{\mathrm{BCD}} \approx \Omega \omega_p$, offering a second contactless detection channel.

Reading between the lines

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

  • If the rocking scales with intensity as predicted, the effect could serve as an all-optical, contactless probe of the Berry curvature dipole, complementing Hall-bar transport measurements.
  • Since the derivation assumes $v_\phi > c$ and a strict traveling-wave ansatz, the same pendulum structure may appear in finite-wavevector or evanescent geometries, where the generalized dispersion relation of Eq. (12) would predict modified rocking periods.
  • The pendulum mapping suggests that at high intensities the polarization can flip for a finite number of cycles when the initial state sits near the effective gravitational minimum, and detecting such flips would be a distinctive signature of the strongly nonlinear regime.
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.

Referee Report

3 major / 5 minor

Summary. The paper studies electromagnetic wave propagation along the stacking direction of a layered time-reversal-invariant metal with a Berry curvature dipole. Starting from Maxwell equations coupled to a collisionless Boltzmann equation, the authors derive an ODE system (Eqs. 6-8) for traveling-wave fields and show that in the strongly nonlinear regime the electric-field vector follows pendulum-like motion. In the weakly nonlinear regime they derive the central result, Eq. (1): the polarization direction rocks as θ(z) ≈ θ0 cos(Ω z/λ) with Ω = (√7/24)(I/I0), while the degree of polarization oscillates at the same rate. The result is supported by numerical integration of the full ODEs in Fig. 2, and an experimental detection strategy is proposed for WTe2, MoTe2, and TaIrTe4.

Significance. The central prediction is attractive and, if correct, would constitute a new nonlinear optical effect that does not require a DC current and is measurable in a contactless geometry. A clear strength is that the weak-field formula is not a fit: the dimensionless ODEs contain no free parameters, and the predicted Ω depends only on I/I0, with the numerical check in Fig. 2 providing direct verification of the reduced model. The paper also makes a falsifiable experimental prediction. However, the analytical derivation of Eq. (22) rests on a second-order expansion that is not a valid solution of the equation of motion, and the conversion from the dimensionless coordinate ξ to the physical z/λ in Eq. (1) appears to miss a factor of 2π. Both issues must be repaired before the central claim is fully supported.

major comments (3)
  1. [Supplementary Section II, Eq. (S15)] The expansion X+iY ≈ [1 + (i/3)(X0−Ẋ0)](X0+iY0) is not a valid solution of Eq. (18). For the linearly polarized initial condition θ0=0, b(0)=0, the zeroth-order solution is Z0 = a0 cos ξ, and the right-hand side of Eq. (18) reduces to −i a0² cos 2ξ. The exact second-order solution satisfying Z1(0)=Ż1(0)=0 is Z1 = (i a0²/3)(cos 2ξ − cos ξ), whereas Eq. (S15) gives Z1 = (i a0²/3)(cos² ξ + sin ξ cos ξ). This violates Y1(0)=0 and does not satisfy Z¨1+Z1 = −i a0² cos 2ξ when substituted. Consequently the averaged equations (S19)–(S27) and the resulting Eq. (22) are not justified by the stated expansion. The numerical agreement in Fig. 2 suggests the final expression may be correct, but the analytical derivation must be repaired or replaced.
  2. [Eq. (1) and Discussion] There is a missing factor of 2π in converting the dimensionless coordinate ξ to the physical distance z/λ. From Eq. (11) and the linear dispersion relation, ξ = ω_p(v_φ t − z)/√(v_φ²−c²) = q(v_φ t − z) = 2π(v_φ t − z)/λ. Since Eq. (22) gives θ(ξ) = θ0 cos(Ω_d ξ) with Ω_d = (√7/24)a0², the physical polarization angle at fixed time is θ(z) = θ0 cos(2π Ω_d z/λ). Equation (1) instead uses Ω_d z/λ, which is inconsistent with the definition of wavelength. The sentence in the Discussion stating a rotation rate of about Ω/2π ∼ 4×10⁻⁷ radians per traveled wavelength also mixes cycles and radians; the units need to be clarified and the numerical estimate corrected.
  3. [Discussion and outlook for experimental detection] The experimental projection directly extrapolates the collisionless, single-band, local nonlinear-Hall model to WTe2, MoTe2, and TaIrTe4, but no estimate is given for the electron collision rate relative to the plasma frequency, for interband absorption at the relevant frequencies, or for wavevector corrections to the local Berry-dipole response. These effects can damp or modify the predicted rocking, so the detection claim needs at least a parametric estimate or an explicit statement of the regime in which the prediction applies.
minor comments (5)
  1. [Eq. (18)] I do not find a sign inconsistency in Eq. (18): using the first integral M = Ẋ/(Y−1), which follows from d[M(Y−1)]/dξ = Ẍ together with the time-reversal-invariant initial conditions, the vector form in Eq. (18) is consistent with Eqs. (6)–(8).
  2. [Page 2, after Eq. (8)] The phrase 'consider a specific solution solution of Eqs. (6)-(8)' contains a duplicated word and should read 'consider a specific solution of Eqs. (6)-(8)'.
  3. [Fig. 1 caption] The word 'exhibtis' should be 'exhibits'.
  4. [Fig. 2 caption] The initial conditions are written as y(0) = a0 cos θ0, x(0) = a0 sin θ0, which is inconsistent with the convention in Eq. (21), where the initial semimajor axis lies at angle θ0 from the x-axis, i.e., X0(0)+iY0(0) = a0 e^{iθ0}. Please reconcile the notation.
  5. [Supplementary Section I] The statement that ω_E ≈ 0.947 is 'very close to one as expected' is not derived; since ω_E is a fitted parameter, the analytical treatment of the strong-field regime should be either completed with a derivation of ω_E or explicitly described as partly phenomenological.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the weak-field polarization rocking formula is derived from the stated Maxwell-Boltzmann equations with no fitted constants; prior self-citations supply input physics, not the predicted result.

full rationale

The central claim, Eq. (1), is derived within the paper from the Maxwell-Boltzmann equations through the traveling-wave reduction to the ODEs (6)-(8), the vector form (18), the second-order perturbative expansion (S15), and the averaging of oscillator constants of motion leading to Eqs. (S26)-(S27). The constants sqrt(7)/24 and the linear intensity dependence emerge from this algebraic expansion and are compared directly to numerical integration of the full ODEs in Fig. 2, so no fitted parameter is renamed as a prediction. The authors' own prior work is cited for the Berry curvature dipole and nonlinear Hall current formulas, which enter as assumptions or inputs, not as the target Faraday-precession result; no uniqueness theorem from prior work is invoked to forbid alternatives, and no ansatz is smuggled in via citation beyond the stated perturbative expansion. The skeptical objection that Eq. (S15) does not satisfy Eq. (18) would, if correct, be a technical derivation error, not circularity: it does not make the output equivalent to an input by construction. The paper is self-contained against its own stated model equations, and the central prediction is not controlled by fitted values or by the authors' prior claims about the final effect.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

The central weak-field prediction does not rely on any fitted constant; material inputs such as the Drude weight D and the Berry dipole d are taken from prior experiment and literature. One envelope frequency is fitted in a secondary strongly-nonlinear comparison. No new particles or forces are introduced.

free parameters (1)
  • omega_E = 0.947
    Envelope frequency of the pendulum energy in the strongly nonlinear limit, fitted to the numerical curve in Fig. 1(c). It is expected to be near 1 and does not enter the central weak-field formula Eq. (1).
assumptions (5)
  • domain assumption Semiclassical Boltzmann dynamics with Berry curvature correction.
    Eq. (5) uses the anomalous velocity and the force from the electric field; this standard model neglects interband coherence and quantum corrections.
  • domain assumption Collisionless regime with no collision integral in the Boltzmann equation.
    Eq. (4) has no collision term and the text states the focus is on the collisionless regime relevant for plasma-like oscillations; scattering would damp or modify the slow polarization dynamics.
  • domain assumption Second-order truncation of the nonlinear current and local response.
    The ODE reduction keeps terms up to second order in the electric field; higher-order nonlinearities and wavevector corrections to the Berry dipole response are neglected.
  • domain assumption Traveling-wave ansatz with phase velocity larger than the speed of light in the material.
    Fields are assumed to depend only on v_phi t - z and v_phi > c is required for the rescaling; the physical dispersion is reconstructed afterward via Eq. (12).
  • domain assumption Time-reversal-invariant initial conditions with no pre-existing Berry curvature imbalance or current.
    The ODE solutions use Xdot(0) = Ydot(0) = M(0) = 0, relevant to time-reversal-invariant systems with no injected current.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Non-linear Faraday Precession of Light Polarization in Time-Reversal Invariant Materials." pith.science (2026). https://pith.science/paper/DU7EPHQQ

@misc{pith2026241203656,
  author       = {Pith},
  title        = {Pith review of: Non-linear Faraday Precession of Light Polarization in Time-Reversal Invariant Materials},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/DU7EPHQQ}},
  note         = {Machine review of arXiv:2412.03656}
}
read the original abstract

We investigate the propagation of electromagnetic waves through materials displaying a non-linear Hall effect. The coupled Maxwell-Boltzmann equations for traveling waves can be mapped onto ordinary differential equations that resemble those for the motion of a pendulum. In the weakly non-linear regime relevant for most experiments, we find that the polarization of light displays a Faraday-like precession of its polarization direction that swings back and forth around the direction of Berry dipole vector as the light beam traverses the material. This occurs concomitantly with an oscillation of its degree of polarization, with a characteristic frequency that increases linearly with the intensity of the traveling wave. These effects could be observed by measuring thickness dependent Faraday rotations as well as the emission of lower frequency radiation associated with the polarization oscillations in materials displaying the non-linear Hall effect.

Figures

Figures reproduced from arXiv: 2412.03656 by the authors.

Figure 1
Figure 1. Numerical solution of Eqs. (6-8) with initial conditions in strongly non-linear regime: [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Numerical solution of Eqs. (6-8) with initial conditions in weakly non-linear regime: [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 20 canonical work pages

  1. [1]

    V. I. Belinicher and B. I. Sturman, Soviet Physics Uspekhi 23, 199 (1980)

  2. [2]

    Belinicher, E

    V. Belinicher, E. Ivchenko, and B. Sturman, Zh. Eksp. Teor. Fiz 83, 649 (1982)

  3. [3]

    von Baltz and W

    R. von Baltz and W. Kraut, Physical Review B23, 5590 (1981)

  4. [4]

    Aversa and J

    C. Aversa and J. E. Sipe, Physical Review B 52, 14636 (1995)

  5. [5]

    S. M. Young and A. M. Rappe, Physical review letters 109, 116601 (2012)

  6. [6]

    Morimoto and N

    T. Morimoto and N. Nagaosa, Science advances 2, e1501524 (2016)

  7. [7]

    Rangel, B

    T. Rangel, B. M. Fregoso, B. S. Mendoza, T. Morimoto, J. E. Moore, and J. B. Neaton, Physical review letters 119, 067402 (2017)

  8. [8]

    A. M. Cook, B. M. Fregoso, F. De Juan, S. Coh, and J. E. Moore, Nature communications 8, 14176 (2017)

Show all 27 references
  1. [9]

    Morimoto, M

    T. Morimoto, M. Nakamura, M. Kawasaki, and N. Nagaosa, Physical review letters 121, 267401 (2018)

  2. [10]

    Kumar, C.-H

    D. Kumar, C.-H. Hsu, R. Sharma, T.-R. Chang, P. Yu, J. Wang, G. Eda, G. Liang, and H. Yang, Nature Nanotechnology 16, 421 (2021)

  3. [11]

    L.-k. Shi, O. Matsyshyn, J. C. Song, and I. S. Villadiego, Physical Review B 107, 125151 (2023)

  4. [12]

    Onishi and L

    Y. Onishi and L. Fu, Physical Review B 110, 075122 (2024)

  5. [13]

    T. G. Rappoport, T. A. Morgado, S. Lanneb `ere, and M. G. Silveirinha, Physical Review Letters130, 076901 (2023)

  6. [14]

    Makushko, S

    P. Makushko, S. Kovalev, Y. Zabila, I. Ilyakov, A. Ponomaryov, A. Arshad, G. L. Prajapati, T. V. de Oliveira, J.-C. Deinert, P. Chekhonin, et al., Nature Electronics 7, 207 (2024)

  7. [15]

    E. Deyo, L. Golub, E. Ivchenko, and B. Spivak, arXiv preprint arXiv:0904.1917 (2009)

  8. [16]

    J. E. Moore and J. Orenstein, Physical review letters105, 026805 (2010)

  9. [17]

    Sodemann and L

    I. Sodemann and L. Fu, Physical review letters 115, 216806 (2015)

  10. [18]

    Ortix, Advanced Quantum Technologies 4, 2100056 (2021)

    C. Ortix, Advanced Quantum Technologies 4, 2100056 (2021)

  11. [19]

    Bandyopadhyay, N

    A. Bandyopadhyay, N. B. Joseph, and A. Narayan, Materials Today Electronics 8, 100101 (2024)

  12. [20]

    Vorob’ev, JETP Lett 29, 441 (1979)

    L. Vorob’ev, JETP Lett 29, 441 (1979)

  13. [21]

    Shalygin, A

    V. Shalygin, A. Sofronov, L. Vorob’ev, and Farbshtein, Physics of the Solid State 54, 2362 (2012)

  14. [22]

    S. S. Tsirkin, P. A. Puente, and I. Souza, Physical Review B97, 035158 (2018)

  15. [23]

    K ¨onig, M

    E. K ¨onig, M. Dzero, A. Levchenko, and D. Pesin, Physical Review B 99, 155404 (2019)

  16. [24]

    K. Kang, T. Li, E. Sohn, J. Shan, and K. F. Mak, Nature materials 18, 324 (2019)

  17. [25]

    Ma, S.-Y

    Q. Ma, S.-Y. Xu, H. Shen, D. MacNeill, V. Fatemi, T.-R. Chang, A. M. Mier Valdivia, S. Wu, Z. Du, C.-H. Hsu, et al., Nature 565, 337 (2019)

  18. [26]

    T. Ma, H. Chen, K. Yananose, X. Zhou, L. Wang, R. Li, Z. Zhu, Z. Wu, Q.-H. Xu, J. Yu,et al., Nature Communications13, 5465 (2022)

  19. [27]

    Kapitulnik, J

    A. Kapitulnik, J. Xia, E. Schemm, and A. Palevski, New Journal of Physics 11, 055060 (2009). 1 Non-linear Faraday Precession of Light Polarization in Time-Reversal Invariant Materials I. TRA VELING W A VES IN STRONGLY NON-LINEAR REGIME Let us first consider large electric fiel...

Pith tools

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