{"id":"a2515874-8e11-40eb-bcc6-a3a1b8939aa2","arxiv_id":"2412.03656","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"In weak fields, the polarization plane of light in a Berry-dipole metal rocks as θ(z) ≈ θ0 cos(Ω z/λ), with Ω proportional to the light intensity.","lead":"This paper predicts that light passing through a layered metal with a non-linear Hall effect rocks its polarization direction back and forth, like a pendulum, instead of rotating steadily. The effect could give a contactless optical way to measure the Berry curvature dipole in time-reversal-invariant materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The analytical derivation of Eq. (1) relies on a second-order expansion (Supp. Eq. S15) that does not satisfy the equation of motion (18); the final formula may be correct, but the derivation needs correction.","rationale":"The paper's central claim is Eq. (1): the polarization direction rocks with Ω=(√7/24)(I/I0). This claim is derived from the Maxwell-Boltzmann reduction to Eqs. (6)–(8), then from a weak-field perturbative expansion in the supplementary. The numerical solution of the ODEs in Fig. 2 agrees with the analytical formula for small amplitudes, so the ODE-level result appears empirically sound. However, the analytical derivation as written contains a concrete internal inconsistency: Eq. (S15), which is used to compute the averaged time derivatives of the oscillator constants, does not actually solve Eq. (18) even to second order. I verified this explicitly for the clean case of linear polarization along the Berry dipole, where the correct second-order solution is Z1=(i a0^2/3)(cos2ξ−cosξ), whereas Eq. (S15) gives a different function that neither satisfies the equation nor the initial conditions. Because Eqs. (S19)–(S27) and ultimately the coefficient √7/24 are derived from Eq. (S15), the analytical route to the headline result is not self-consistent. The numerical check in Fig. 2 suggests that the final formula may nevertheless be correct, and a corrected derivation might well recover √7/24. I therefore do not reject the paper, but I would require the authors to fix the expansion or provide an alternative derivation before accepting. The reader's identified weakest assumption was the collisionless, single-band material model; that is a valid concern about applicability to real WTe2/MoTe2/TaIrTe4 samples, but it is not the same as the internal inconsistency I found. The reader's additional comment about an apparent sign inconsistency in Eq. (18) is not correct: the relation M=Ẋ/(Y−1) follows from the first integral Ẋ−M(Y−1)=0 under the stated time-reversal-invariant initial conditions, and Eq. (18) is consistent with that relation. Thus the stress-test identifies a different, more concrete load-bearing concern, and the verdict remains CONDITIONAL pending correction of the analytical derivation.","tokens_in":11501,"tokens_out":39235,"duration_ms":355080,"concrete_test":"Re-derive the second-order solution of Eq. (18) for θ0=0 using variation of parameters, and compare it with Eq. (S15). Concretely, substitute Z1=(i a0^2/3)(cos2ξ−cosξ) into Eq. (18) and confirm it satisfies the equation and initial conditions, whereas Eq. (S15) does not. Then recompute the averaged angular-momentum derivative ⟨Ldot⟩ and ⟨ε_y dot⟩ using the correct expansion to see whether Eqs. (S19)–(S27) and the coefficient √7/24 in Eq. (22) survive. If the coefficient changes, the central claim Eq. (1) must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central weak-field result is derived in Supplementary Section II from the expansion X+iY ≈ [1 + (i/3)(X0−Ẋ0)](X0+iY0) (Eq. S15). This expansion is not a valid solution of Eq. (18). For a beam initially polarized along the Berry dipole (θ0=0, b(0)=0), the leading solution is Z0=a0 cosξ. The exact second-order solution with Z1(0)=Ż1(0)=0 is Z1=(i a0^2/3)(cos2ξ−cosξ), which satisfies Z¨1+Z1=−i a0^2 cos2ξ. Eq. (S15) instead gives Z1=(i a0^2/3)(cos^2ξ+sinξ cosξ); substituting this yields Z¨1+Z1=i a0^2/6−i a0^2(cos2ξ+sin2ξ)/2, which does not equal the required forcing. The error is also visible at ξ=0, where Eq. (S15) gives Z1(0)=i a0^2/3 instead of zero. The averaged equations (S19)–(S27) are therefore not justified by the stated expansion. The numerical agreement in Fig. 2 suggests the final Ω may be correct for the ODE system, but the analytical derivation must be repaired or replaced. Note: the reader's concern about a sign inconsistency in Eq. (18) is not supported; M=Ẋ/(Y−1) follows from the first integral Ẋ−M(Y−1)=0 under time-reversal-invariant initial conditions, and Eq. (18) is consistent with it.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":11826,"tokens_out":16959,"duration_ms":160894,"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":[{"comment":"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.","section":"Supplementary Section II, Eq. (S15)"},{"comment":"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.","section":"Eq. (1) and Discussion"},{"comment":"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.","section":"Discussion and outlook for experimental detection"}],"minor_comments":[{"comment":"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).","section":"Eq. (18)"},{"comment":"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)'.","section":"Page 2, after Eq. (8)"},{"comment":"The word 'exhibtis' should be 'exhibits'.","section":"Fig. 1 caption"},{"comment":"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.","section":"Fig. 2 caption"},{"comment":"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.","section":"Supplementary Section I"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and the central idea is interesting. The two load-bearing issues—the invalid second-order expansion in the supplementary material and the missing factor of 2π in Eq. (1)—are repairable within the manuscript's framework, so rejection is not warranted. I do not see a circularity problem in the use of Refs. 11 and 17; they are the standard sources for the Berry curvature dipole and the nonlinear Hall effect, and they are used as definitions rather than as support for the new result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's what I make of it. The paper makes a genuinely new prediction: in a time-reversal-invariant layered metal with a Berry curvature dipole, a traveling wave's polarization rocks back and forth with propagation distance, with a rate linear in intensity. That's not in the earlier nonlinear Hall work, which treated the drive as fixed, and it does not need a DC current the way current-induced optical activity does. The reduction of Maxwell-Boltzmann to the three ODEs (6)-(8) is clean, and the numerical check in Fig. 2 shows the rocking behavior and the frequency ∝ a0^2/√7 in the full ODE system. That part is solid and worth taking seriously.\n\nThe problem is in the analytical derivation of the central formula. The expansion in Supp. Eq. (S15), Z = [1 + (i/3)(X0 - X0dot)]Z0, is not a solution of Eq. (18) at second order. Take the simple case of initial linear polarization along the Berry dipole, θ0=0. The exact second-order correction with Z1(0)=Z1dot(0)=0 is Z1 = (i a0^2/3)(cos2ξ - cosξ). Equation (S15) instead gives Z1 = (i a0^2/3)(cos^2ξ + sinξ cosξ), which has Z1(0)=i a0^2/3 and fails the forced equation. So the averaged equations (S19)-(S27) are not justified by the stated expansion. The numerical agreement in Fig. 2 suggests the final Ω may be the right answer for the ODE system, but the derivation as written has to be repaired or replaced. This is a load-bearing flaw, not a typo.\n\nTwo smaller points. The reader flagged a sign inconsistency in Eq. (18), but that's a false alarm: M = Xdot/(Y-1) follows from the first integral and Eq. (18) is consistent with it. The strongly nonlinear section uses one fitted envelope frequency ω_E ≈ 0.947; that's minor and clearly labeled. The collisionless single-band model is a stated idealization, so I would not treat the numbers for WTe2 and MoTe2 as quantitative, but that's a scope limitation rather than an error.\n\nWho is this for? People in nonlinear Hall and Berry-phase optics. It deserves a serious referee: the physics is novel, the numerics support the phenomenon, and the flaw is fixable. I would not cite Eq. (1) in my own work until the derivation is corrected, but I would send this to review and ask for the weak-field calculation to be redone properly.","headline":"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.","tokens_in":12343,"tokens_out":5748,"would_cite":false,"duration_ms":51830,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["78.20.Fm"],"model":"deepseek-v4-flash","headline":"Light polarization rocks back and forth in nonlinear Hall metals","keywords":["nonlinear Hall effect","Berry curvature dipole","Faraday rotation","polarization precession","plasma oscillations","Boltzmann equation","layered metals","time-reversal invariant"],"falsifier":"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.","tokens_in":11287,"feed_emoji":"🔄","tokens_out":5926,"duration_ms":48092,"temperature":0.7,"pith_summary":"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.","feed_headline":"Light's polarization rocks back and forth in nonlinear metals","feed_subtitle":"Theory maps Maxwell-Boltzmann equations to a pendulum and predicts a measurable rocking of the polarization plane.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Defines the Berry curvature dipole tensor that controls the nonlinear Hall effect used throughout.","marker":"[17]"},{"why":"Reports the nonlinear Hall effect in WTe2, supplying a candidate material and parameter scale.","marker":"[24]"},{"why":"Provides a second experimental observation of the nonlinear Hall effect in WTe2, supporting the material choice.","marker":"[25]"},{"why":"Measures the nonlinear Hall conductivity in MoTe2, from which the paper extracts $E_0 \\approx 10^8$ V/m and $I_0 \\approx 3 \\times 10^{13}$ W/m$^2$.","marker":"[26]"},{"why":"Documents interferometric rotation sensitivities near $10^{-9}$ radians, the detection threshold for the predicted effect.","marker":"[27]"}],"fun_headline_variants":["Nonlinear Hall effect swings light's polarization","Light's polarization rocks in nonlinear metals","Pendulum-like swing predicted for light in nonlinear Hall","Berry curvature makes polarization precess back and forth","Nonlinear Faraday precession: light's polarization oscillates"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Nonlinear Hall effect swings light's polarization","Light's polarization rocks in nonlinear metals","Pendulum-like swing predicted for light in nonlinear Hall","Berry curvature makes polarization precess back and forth","Nonlinear Faraday precession: light's polarization oscillates"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00031,"raw_usage":{"total_tokens":1745,"prompt_tokens":899,"completion_tokens":846,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":774}},"tokens_in":515,"tokens_out":846,"duration_ms":8604,"temperature":1.0,"reasoning_tokens":774,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:15:07.493073+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the nonlinear Hall effect in WTe2, supplying a candidate material and parameter scale."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Measures the nonlinear Hall conductivity in MoTe2, from which the paper extracts $E_0 \\approx 10^8$ V/m and $I_0 \\approx 3 \\times 10^{13}$ W/m$^2$."},{"cited_title":"Kapitulnik, J","cited_arxiv_id":null,"evidence_quote":"Documents interferometric rotation sensitivities near $10^{-9}$ radians, the detection threshold for the predicted effect."}],"review_version":1}