{"id":"f9e98592-6c32-4304-a49f-622ab9e2fba7","arxiv_id":"2411.11156","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A new disk-geometry analysis unambiguously extracts all six nonlinear conductivity components and reveals a giant, symmetry-broken nonlinear Hall effect in Bernal bilayer graphene.","lead":"Researchers derived the analytical potential distribution for nonlinear current flow through a disk-shaped sample and used it to extract the full nonlinear conductivity tensor from angle-resolved measurements. They report a giant nonlinear Hall response in Bernal bilayer graphene whose magnitude exceeds known microscopic mechanisms by orders of magnitude.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The extraction is only as strong as the assumed uniform, local, diffusive transport; nonlinear locality is asserted, not independently tested, so the reported C could be an effective fitting artifact.","rationale":"The paper is technically impressive: the analytic disk solution is explicit, the fit is overdetermined (48 measurements for 6 parameters, 72 for sample 2), the second sample reproduces the main observations, and the linear-RRMSE analysis does support diffusive transport at the selected operating points. The reader's conditional verdict is therefore reasonable. However, the reader's weakest-assumption identification focuses on the Joule-Lenz decomposition and the lack of consensus on nonlinear Onsager relations. That issue is less load-bearing than stated: the extracted C is defined as the transverse, non-dissipative (j_nd dot E = 0) part of the local nonlinear current, and this is a direct algebraic decomposition of the fitted tensor, not an Onsager-dependent choice. The more fragile premise is the assumed spatial uniformity and locality of the nonlinear conductivity itself. The linear fits show that the linear response is described by a single local tensor, but the nonlinear term is a new object; a spatially textured nonlinear response could produce angular signatures that the uniform model would misattribute to A, B, and C. The paper's 'remarkable agreement' is the main evidence against this, but a statistical cross-validation of the nonlinear fit and subset stability checks would make the uniformity claim conclusive. Since the requested concrete test is feasible from existing data, the appropriate verdict remains CONDITIONAL rather than ACCEPT or REJECT, and no verdict change is needed from the reader's assessment.","tokens_in":91228,"tokens_out":20420,"duration_ms":211587,"concrete_test":"Refit the six nonlinear parameters independently on disjoint subsets of the 48/72 measurements: (a) configurations with source-drain opening angle near 180 degrees versus near 90 degrees, and (b) the first four versus last four rotation angles within each configuration. If the extracted A, B, C disagree beyond propagated fitting error, or if the full-model residual significantly exceeds the 2-omega noise floor, the uniform-local model fails for the nonlinear response. As a calibration, also generate synthetic V2-omega data from the extracted tensor with the measured noise level, refit, and confirm that C is recovered within error; this checks whether the specific wiring set actually identifies all six parameters rather than relying on the claimed overdetermination.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim has two steps: solve the disk electrostatics with uniform, local, diffusive conductivities (SI Sec. III states these are the only two fundamental assumptions), then fit the six nonlinear parameters to 48/72 angle-resolved points. The excellent visual agreement of the fits is real evidence, but it does not by itself validate the model at the nonlinear level. The paper supports the linear conductivity with linear-RRMSE fits and uses RRMSE to map the diffusive regime (SI Fig. S9), yet the nonlinear extraction is not cross-validated in the same way: no nonlinear RRMSE is reported against the 2-omega noise floor, no subset stability of the extracted A, B, C vectors is shown, and the two samples have the same disk diameter, so finite-size or ballistic corrections cannot be distinguished. If the nonlinear conductivity is spatially nonuniform (for example due to isospin domain textures in PIP2 or Sym4), or if near-ballistic corrections in Sym4 modify the angular pattern, the analytic basis functions (Eqs. 4-6) are the wrong model; the least-squares inversion will nevertheless return six parameters, and the resulting C would be an effective artifact rather than a local nonlinear Hall conductivity. The Onsager/decomposition issue flagged by the reader is less decisive here: the C vector is, by construction, the transverse j dot E = 0 part of the nonlinear tensor, so the load-bearing assumption is locality and diffusivity, not the Joule-Lenz labeling. For this reason the paper remains CONDITIONAL, not ACCEPT.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytical solution for the second-order (nonlinear) electrostatic potential in a disk-shaped two-dimensional conductor with arbitrary uniform linear and nonlinear conductivity tensors, expressed as a linear combination of three basis functions (Eqs. 4-6). The authors apply this to angle-resolved second-harmonic measurements on Bernal bilayer graphene disks with eight leads, fitting the six components of the nonlinear tensor simultaneously from 48 (Sample 1) or 72 (Sample 2) independent wiring installations. They report a large non-dissipative nonlinear Hall-like component C (e.g., 6.88±2.05 µm/(Ω·V) in PIP2 and 668±339 µm/(Ω·V) in Sym4), a small threefold component, a reorientation of B and C across the PIP2-Sym4 transition, and invariance of the C-vector under D reversal, which they argue rules out the Berry-curvature-dipole mechanism. They conclude that none of the known microscopic mechanisms can account for the magnitude and structure.","tokens_in":91601,"tokens_out":5277,"duration_ms":57971,"significance":"If the extraction is valid, the paper provides a valuable general framework for nonlinear transport measurements, solving a real ambiguity in Hall-bar geometries (the Methods section shows that a Hall bar cannot disentangle C from the 3-fold terms). The analytical basis functions and the overdetermined fitting procedure are strong assets; the two-sample reproducibility, the quadratic-regime check, the quantum oscillation phase identification, and the D-reversal test are genuine independent checks. However, the central quantitative claim of a giant nonlinear Hall conductivity is conditional on the assumed uniform, local, diffusive transport at the nonlinear level, which is not independently validated, and on the Joule-Lenz criterion for separating dissipative from non-dissipative currents.","major_comments":[{"comment":"The extraction rests on the assumption of uniform, local, diffusive conductivities; SI Sec. III identifies these as the only fundamental assumptions. The linear full fits and linear RRMSE (SI Figs. S3, S4, S9) support the linear part of that assumption, but no nonlinear analogue is reported: there is no nonlinear RRMSE against the 2ω noise floor, no stability test of the extracted A, B, C when subsets of the 48/72 wiring installations are removed, and the theoretical curves in Figs. 3a, 4a-b, S2c-d are generated by the same fitted tensor, so the close agreement is an in-sample fit-quality measure rather than a prediction. This leaves open the possibility that a spatially nonuniform or nonlocal nonlinear conductivity (e.g., isospin domain textures, or near-ballistic corrections) would produce an equally good six-parameter fit with different effective A, B, C. I would like to see a nonlinear RRMSE, a leave-one-configuration-out or bootstrap stability analysis, and, ideally, a prediction for configurations not used in the fit.","section":"SI Sec. III; Eqs. (4)-(6); Fig. 3a"},{"comment":"Both samples have the same 2 µm disk diameter, so the model cannot be tested against finite-size or ballistic corrections by comparing samples. The paper itself notes that the larger Sym4 uncertainties come from proximity to the ballistic regime, and SI Fig. S9 shows a sharp RRMSE onset as the sample leaves the diffusive regime. If in the Sym4 phase the mean free path (~1 µm from SI Fig. S4) is not negligible compared with the disk diameter, the analytic basis functions (Eqs. 4-6) may be the wrong model, and the least-squares inversion would still return six parameters, turning C into an effective fitting artifact rather than a local nonlinear Hall conductivity. A sample with a different diameter, or an explicit estimate of finite-size corrections, is needed to support the reported magnitudes, especially the Sym4 value.","section":"SI Sec. I.B; SI Sec. I.C; main text around Sym4 uncertainties"},{"comment":"The identification of the extracted C as a 'nonlinear Hall conductivity' relies on the Joule-Lenz criterion j_nd·E = 0 to separate non-dissipative from dissipative parts of the nonlinear tensor. The authors correctly note that there is no consensus on Onsager relations for the nonlinear conductivity tensor. Since the magnitude, orientation, and D-reversal behavior of C are all defined with respect to this criterion, the central claim is not model-independent. This is not a circularity, but it is a correctness risk: if a different decomposition is physically mandated, the quoted C values and the conclusion that the response is non-dissipative would change. A concrete test would be to derive or adopt Onsager-type constraints for the second-order response and check whether the reported C remains non-dissipative under that alternative.","section":"SI Sec. II.D"}],"minor_comments":[{"comment":"The text says 'Eq. (656) is an exact analog of the linear least squares fit'; this should be Eq. (13).","section":"Methods, 'Details of the nonlinear conductivity extraction'"},{"comment":"The phrase 'unambiguously extract' should be qualified as 'within the uniform, local, diffusive model', since the extraction is unambiguous only conditional on those assumptions.","section":"Abstract and main text"},{"comment":"There is a typo in the first paragraph: 'exmaple' should be 'example'.","section":"SI Sec. I.A"},{"comment":"The phrase 'with 2 large Fermi surfaces (FSs)' is awkward; consider 'with two large Fermi surfaces' for consistency with the later '4 large Fermi surfaces'.","section":"Main text, phase identification paragraph"},{"comment":"No data or code availability statement is provided; given the methodological claims, depositing the fitting code and the processed angle-resolved datasets would substantially strengthen the reproducibility of the extraction.","section":"Reproducibility"}],"recommendation":"major_revision","confidential_remarks":"This is a technically strong paper with a genuinely useful methodology, but the central experimental claim is less secure than the presentation suggests because the nonlinear model is not independently validated. The stress-test concern about nonlinear locality is the decisive one; the Onsager/decomposition issue is secondary but should be clarified. I would be willing to reconsider after the authors provide nonlinear cross-validation (e.g., nonlinear RRMSE, subset stability, withheld-configuration predictions) and explicitly address finite-size corrections, especially for the Sym4 phase."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is a methods paper that earns its keep. The authors analytically solve the nonlinear potential problem on a disk with an arbitrary uniform conductivity tensor and give closed-form basis functions (Eqs. 4-6). That is a real step beyond the Hall bar, which only sees two of six tensor components and cannot separate the Hall piece from the 3-fold dissipative piece. The measurement protocol—six configurations by eight wirings, 48 points fitted with six parameters—is well overdetermined, and the linear fits put the diffusive regime on reasonable footing.\n\nThe extraction itself is convincing. The in-sample agreement is visually excellent, and the D-reversal invariance plus the roughly 180-degree reorientation across PIP2/Sym4 are checks that do not depend on the fit. Sample 2 reproduces the structure. I believe the tensor extraction works within the stated model.\n\nSoft spots, in proportion. The model assumes uniform, local, diffusive transport; the authors say that themselves in SI Sec. III. The problem is that the nonlinear regime is never cross-validated. There is no nonlinear RRMSE against the 2-omega noise floor, no subset-stability of A/B/C, and both samples have the same 2-micrometer disk, so ballistic or finite-size corrections are not distinguished. The angular agreement is fit quality, not prediction. So 'unambiguously extract' is too strong; 'extract within the uniform local model' is accurate.\n\nSecond, the dissipative/non-dissipative split follows the Joule-Lenz criterion j dot E = 0. The authors note there is no consensus on nonlinear Onsager relations. That does not invalidate the math—C is defined as the part orthogonal to E—but the label 'nonlinear Hall conductivity' carries physical weight the criterion alone doesn't supply. The magnitudes are three to seven orders above known mechanisms; that is a genuine mystery, and the paper says so.\n\nThird, no data or code are released. For a measurement-plus-analysis framework, reviewers should ask for them.\n\nBottom line: the paper deserves a serious referee. The method is citable and will be reused; the BBG observation is interesting and possibly important. Value comes to experimentalists in nonlinear transport and theorists who analyze such data. I'd send it to review, with a request for data/code and one nonlinear validation check. The claim should be softened in the abstract.","headline":"A genuinely useful disk-geometry toolkit for extracting the nonlinear conductivity tensor; the giant nonlinear Hall claim is credible but 'unambiguous' outruns the uniformity assumption.","tokens_in":92145,"tokens_out":2672,"would_cite":true,"duration_ms":73782,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A disk-geometry analytical solution extracts the full second-order conductivity tensor from angle-resolved measurements in Bernal bilayer graphene, revealing a giant non-dissipative nonlinear Hall conductivity up to $668\\…","keywords":["nonlinear Hall effect","Bernal bilayer graphene","nonlinear conductivity tensor","disk geometry","second-harmonic transport","angle-resolved transport","Berry curvature dipole","rotational symmetry breaking"],"falsifier":"A direct calorimetric experiment on the same disk—measuring the local heat produced by the second-order current and checking whether the component identified as $C$ generates zero heat for all orientations of the electric field—would settle whether the Joule–Lenz split is correct.","tokens_in":91070,"feed_emoji":"⚡","tokens_out":8564,"duration_ms":76754,"temperature":0.7,"pith_summary":"In two-dimensional electron systems, the second-order (I²) current is described by a rank-3 conductivity tensor with six independent components, and measuring them requires more than a Hall-bar geometry. This paper solves the nonlinear electrostatics problem for a disk-shaped sample with arbitrary linear and nonlinear tensors, giving the full spatial distribution of the I² potential in closed form. Applying the solution to angle-resolved second-harmonic measurements on Bernal bilayer graphene, the authors extract all six components and identify a nonlinear Hall vector $C$ with magnitude up to $668\\ \\mu$m/($\\Omega\\cdot$V) in the Sym4 phase and $6.88\\ \\mu$m/($\\Omega\\cdot$V) in the PIP2 phase—values orders of magnitude beyond known mechanisms. They also find that the dissipative and Hall parts share a single mirror axis, rotate almost 180° across the PIP2–Sym4 transition, and are insensitive to reversal of the displacement field, ruling out the Berry-curvature-dipole mechanism. The broader claim is that disk-based angle-resolved nonlinear transport is a systematic tool for identifying the symmetry and magnitude of nonlinear response in quasi-2D materials.","feed_headline":"Giant nonlinear Hall effect found in bilayer graphene via disk probe","feed_subtitle":"An analytic disk solution extracts all six nonlinear tensor components; the Hall part is orders of magnitude beyond known mechanisms.","key_machinery":"The load-bearing object is the closed-form solution for the $I^2$ electrostatic potential $\\Phi^{(2)}$ on a disk with arbitrary uniform linear and nonlinear conductivity tensors. For isotropic linear response, the potential at the boundary is a linear combination of three basis functions, each proportional to one of the complex parameters $\\Xi^{(2)}_-$, $\\Xi^{(2)}_+$, and $\\Xi^{(2)}_0$, which depend on the source-drain angle with distinct phase windings. Fitting these basis functions to 48 independent wiring configurations yields the six real components of the nonlinear tensor. The non-dissipative Hall part is isolated by the Joule–Lenz criterion $j_{\\rm nd}\\cdot E = 0$, which selects the component $C$ with $j^{\\rm nl}_{\\rm Hall} = \\hat{z} \\times E\\,(C\\cdot E)$; the dissipative longitudinal part $B$ and the threefold part $A$ are defined analogously. The derivation uses the Neumann Green's function of the disk (and of the ellipse for anisotropic linear response) to convert the Poisson equation for $\\Phi^{(2)}$ into a quadrature.","core_discovery":"The paper's central discovery is that the full DC nonlinear conductivity tensor of a two-dimensional electron system can be extracted unambiguously from second-harmonic measurements on a disk, and that in Bernal bilayer graphene this extraction reveals a giant nonlinear Hall effect. The nonlinear current is decomposed into three contributions: a threefold-symmetric term $A$, a purely longitudinal dissipative term $B$, and a non-dissipative Hall term $C$, the last of which is the nonlinear Hall vector. In the measured phases, $A$ is vanishingly small while $B$ and $C$ are large—$|C| = 6.88 \\pm 2.05\\ \\mu$m/($\\Omega\\cdot$V) in the partially isospin-polarized phase (PIP2) and $|C| = 668 \\pm 339\\ \\mu$m/($\\Omega\\cdot$V) in the symmetric four-Fermi-surface phase (Sym4)—and the two vectors are orthogonal, indicating a single mirror plane. The authors show that neither the Berry curvature dipole mechanism nor disorder, interaction, or quantum-geometric mechanisms can account for the structure or magnitude of the response, and that the $C$ and $B$ vectors rotate by roughly 180° across the PIP2–Sym4 transition while remaining robust within each phase.","pith_inferences":["A direct test of the Joule–Lenz split would be to measure local heat generation at second harmonic: if the component identified as $C$ is truly non-dissipative, it should produce zero heat for every orientation of the electric field.","The same disk formalism could be extended to third- and higher-harmonic responses, where the tensor structure has more independent components; analogous basis-function decompositions would provide a systematic route to extracting those higher-order tensors.","If the giant $C$ reflects a spontaneous order parameter that breaks rotational symmetry, temperature- and field-dependent measurements of the $C$-vector orientation could map the symmetry axes of that order parameter without requiring a Hall bar.","The extraction assumes local, uniform, diffusive transport; extending the basis-function decomposition to nonlocal or ballistic corrections would test whether the reported magnitude of $C$ survives beyond the diffusive regime."],"forward_implications":["Hall-bar geometries can never uniquely determine the nonlinear Hall effect: only two of the six nonlinear tensor components are accessible, and even with mirror symmetry the equations remain underdetermined, so future claims of nonlinear Hall response need disk-type multi-contact measurements.","The same disk solution provides a ready-made analysis tool for any quasi-2D material: after measuring a small number of angle-resolved configurations, the full tensor can be predicted for all 840 possible wiring installations of an eight-contact disk.","The observed giant $C$ values and their insensitivity to displacement-field reversal rule out the Berry curvature dipole as the origin of the nonlinear Hall effect in Bernal bilayer graphene.","The roughly 180° rotation of $B$ and $C$ across the PIP2–Sym4 transition shows that the second-harmonic response is a sensitive probe of isospin ordering, able to detect symmetry changes that do not show up in the linear conductivity.","Because the threefold component $A$ is negligibly small, skew scattering and side-jump mechanisms—which necessarily produce a large threefold contribution—are excluded, putting a strong constraint on microscopic theories."],"supporting_citations":[{"why":"Supplies the exact solution for the linear electrostatic potential on a disk, which is the starting point for the nonlinear extraction and is used to fit the linear conductivity tensor.","marker":"[30]"},{"why":"Introduces the Berry curvature dipole mechanism and the coordinate-free form $C = D$ that the authors compare against and ultimately exclude as the origin of the observed signal.","marker":"[6]"},{"why":"Provides the angle-resolved transport and relative-root-mean-square-error (RRMSE) fitting methodology, and reports a similar apparent decoupling between linear and nonlinear rotational symmetry.","marker":"[12]"},{"why":"Reports the previously observed giant nonlinear conductivity in graphene moiré superstructures, which the PIP2 value here exceeds by about a factor of two.","marker":"[29]"},{"why":"Provides the compressibility phase diagram used to identify the PIP2 and Sym4 phases and the transition line across which the $C$-vector rotates.","marker":"[32]"},{"why":"Early experimental report of a nonlinear Hall effect in WTe2, used as an example where a transverse nonlinear voltage was interpreted as nonlinear Hall.","marker":"[7]"},{"why":"Hall-bar observation of the time-reversal-symmetric nonlinear Hall effect; the paper argues that such data cannot disentangle the nonlinear Hall component from the threefold component.","marker":"[8]"},{"why":"Provides an estimate of the skew-scattering contribution in graphene-based systems (about $0.01\\ \\mu$m/($\\Omega\\cdot$V)), used to show that known mechanisms fall orders of magnitude short.","marker":"[19]"}],"fun_headline_variants":["Giant nonlinear Hall effect in bilayer graphene from disk measurements","Analytic disk solution extracts full tensor, finds giant Hall effect","Unexpectedly large nonlinear Hall response in Bernal bilayer graphene","Disk probe unlocks nonlinear tensor, revealing giant Hall effect in graphene"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that the extracted component $C$ is a genuine nonlinear Hall conductivity rests on the assumption that the heat-free part of the nonlinear current—defined by $j\\cdot E = 0$—is the correct definition of the non-dissipative nonlinear response; if a different separation is physical, the magnitude, orientation, and displacement-field-reversal invariance of $C$ would no longer uniquely determine a nonlinear Hall effect.","fun_headline_variants_meta":{"raw":{"variants":["Giant nonlinear Hall effect in bilayer graphene from disk measurements","Analytic disk solution extracts full tensor, finds giant Hall effect","Unexpectedly large nonlinear Hall response in Bernal bilayer graphene","Disk probe unlocks nonlinear tensor, revealing giant Hall effect in graphene"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000675,"raw_usage":{"total_tokens":3084,"prompt_tokens":971,"completion_tokens":2113,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":587,"completion_tokens_details":{"reasoning_tokens":2043}},"tokens_in":587,"tokens_out":2113,"duration_ms":92353,"temperature":1.0,"reasoning_tokens":2043,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:52:03.245998+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct calorimetric experiment on the same disk—measuring the local heat produced by the second-order current and checking whether the component identified as $C$ generates zero heat for all orientations of the electric field—would settle whether the Joule–Lenz split is correct.","supporting_citations":[{"cited_title":"He , author G","cited_arxiv_id":null,"evidence_quote":"Reports the previously observed giant nonlinear conductivity in graphene moiré superstructures, which the PIP2 value here exceeds by about a factor of two."}],"review_version":1}