{"id":"48f0a5e7-0a87-42a7-b7be-e7bc8fc7fb31","arxiv_id":"2506.21464","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A kinetic model shows that an anisotropic Fermi surface in 2D materials drives a confined, purely growing TE electromagnetic instability that generates out-of-plane magnetic fields.","lead":"This paper predicts that 2D materials with an anisotropic Fermi surface can host a purely growing electromagnetic instability that generates out-of-plane magnetic fields without an oscillating drive. The finding could offer a new way to create tunable magnetic textures in atomically thin materials, though only transient structures are demonstrated so far.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Collisionless Vlasov assumption is the load-bearing gap; no estimate of collisional damping vs growth rate is given, so the material-level claim remains untested.","rationale":"The reader's weakest assumption identifies the collisionless Vlasov equation as the main vulnerability, and I agree that this is the most load-bearing concern. However, the mechanism is slightly different from the reader's phrasing: for an anisotropic Fermi surface that is an equilibrium property of the band structure, collisions do not relax the anisotropy itself; they damp the collective mode. The paper's own Section V limitation statement is therefore in-scope evidence and must be weighed, as the instructions require. The manuscript offers no quantitative comparison between the growth rate and realistic collision rates, so the step from the idealized kinetic model to 'bidimensional materials' is unsupported. I do not see an internal contradiction in the linear analysis: the secular equation (7) follows from the linearized Vlasov-Maxwell system, the anisotropic case yields a purely growing mode with Re ω=0, and the sign of the anisotropy term is consistent with the instability band q<q_crit. The paper's bespoke simulations appear to confirm the linear growth rate, which is independent supporting evidence. A secondary issue is a likely typographical error in Eq. (12), where the Fermi-Dirac exponent contains μ_x/μ_y rather than μ_x/T; the subsequent use of polylogarithms with α=μ_x/T indicates the intended distribution uses μ_x/T. This is a correctable typo, not a fundamental flaw, and it does not change the verdict. Since the reader already assigned CONDITIONAL based on the collisionless assumption, and my analysis supports that assessment without moving it to a different verdict, the appropriate outcome is UNCHANGED: CONDITIONAL remains the right verdict.","tokens_in":6256,"tokens_out":21408,"duration_ms":239136,"concrete_test":"Add a Bhatnagar-Gross-Krook collision term −ν(f−f0) to Eq. (3), re-derive the linearized dispersion relation for the quadratic-band case, and compute the threshold collision frequency ν_c at which the maximum growth rate first vanishes. Then compare ν_c with the electron-electron or electron-phonon scattering rate in a concrete candidate, e.g., trilayer graphene with α=10 and γ=5 as in Fig. 1. If ν_scattering > ν_c, the instability would not grow; if ν_scattering < ν_c, the collisionless prediction is robust. The same check can be run in the Quetzal solver by adding the BGK term and observing whether the magnetic-field growth is suppressed for ν = 0.1 Γ_max and ν = Γ_max.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires the anisotropic Fermi surface to persist long enough for the purely growing TE instability to amplify. Section V explicitly states that in dense, cold 2D conductors 'it is unlikely that such an anisotropy can be maintained, due to the presence of collisions and their relaxation of the distribution function towards equilibrium.' For a non-equilibrium temperature anisotropy that statement is exactly right, but for an anisotropic Fermi surface arising from band structure the equilibrium distribution is itself anisotropic, and collisions drive the system toward that equilibrium rather than away from it. The actual gap is that Eq. (3) is a collisionless Vlasov equation, and the growth rates derived in Eqs. (13) and (15) therefore describe a collisionless system. If the momentum-relaxation rate ν exceeds the maximum growth rate Γ_max, the mode is overdamped and no magnetic field is generated. The manuscript does not provide an estimate of ν for any candidate material, nor a value of Γ_max in physical units, so the realizability of the central claim in a real bidimensional material is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a collisionless kinetic model of electrons in a two-dimensional material with an anisotropic Fermi surface, considering both quadratic and linear low-energy band dispersions. Linearizing the Vlasov equation (3) about the anisotropic Fermi-Dirac equilibria (11) and (12) and coupling the response to the modified Maxwell system with sheet-current boundary conditions (4)-(5), the authors derive the secular equation (7) and the approximate dispersion relations (13) and (15). For both band types the unstable mode is purely growing (Re omega = 0), unstable wavenumbers lie in a finite band bounded by qcrit in Eq. (17), and the fields are confined along z because kappa_z = sqrt(q^2 + Gamma^2/c^2) becomes real for growing modes. Fully kinetic simulations reproduce the theoretical growth rates in the linear phase and, in the nonlinear phase, show the formation and turbulent breakdown of out-of-plane magnetic field structures.","tokens_in":6371,"tokens_out":15661,"duration_ms":191502,"significance":"If the results hold, they identify a new route to spontaneous out-of-plane magnetic field generation in 2D materials, with explicit, analytically derived growth rates and cutoff wavenumbers that are not fitted to the simulations. The contrast between unconfined stable TE waves and confined growing modes is a sharp, testable prediction. The numerical code validates the linear algebra within the same model, although it does not independently benchmark the model assumptions. The main caveat is that the entire analysis is collisionless, and the material-level realizability of the instability depends on growth being faster than momentum relaxation, which the manuscript does not quantify.","major_comments":[{"comment":"The Vlasov equation (3) has no collision operator, so the growth rates in Eqs. (13) and (15) are those of an ideal collisionless system. Section V correctly notes that in dense, cold 2D conductors a temperature anisotropy would be relaxed by collisions, and it then shifts to an anisotropic Fermi surface; for that equilibrium the collision operator does not erase the equilibrium anisotropy, but it does damp the perturbation. The manuscript never estimates the momentum relaxation rate nu for the cited materials ([16]-[19]) nor gives Gamma_max in physical units, so the condition Gamma_max >> nu is left unexamined. Without such an estimate, the claim that the instability generates magnetic fields 'in bidimensional materials' is not established; please add a relaxation-time estimate for at least one candidate material, or explicitly restrict the conclusions to the collisionless model.","section":"Section V, Eq. (3)"},{"comment":"The linear-band equilibrium distribution is printed as f0 = [1 + exp((v_F/T) sqrt(p_x^2 + (mu_x^2/mu_y^2) p_y^2) - mu_x/mu_y)]^{-1}; the subtracted quantity is the anisotropy ratio mu_x/mu_y, not mu_x/T. Eq. (15) and Fig. 1 use alpha = mu_x/T through Li_n(-e^alpha), so the stated equilibrium is not the Fermi-Dirac distribution on which the linear-band dispersion relation is based. Please correct Eq. (12) to subtract mu_x/T and re-verify the integrals leading to Eq. (15).","section":"Eq. (12)"}],"minor_comments":[{"comment":"The text in Section V.A says the figure shows the linear-band case, while the Fig. 2 caption says 'quadratic dispersion charge carriers'; please reconcile this discrepancy.","section":"Section V.A, Fig. 2"},{"comment":"The sentence introducing the boundary conditions contains a duplicated 'and' ('continuity conditions ... and and discontinuity condition'); please proofread.","section":"Section IV, Eqs. (4a)-(4b)"},{"comment":"The text describes Eq. (11) as an anisotropy in the chemical potential, but the expression implements an anisotropic effective mass through the factor mu_x/mu_y multiplying p_y^2 with a single chemical potential mu_x; please clarify the intended physical meaning and notation.","section":"Eq. (11)"},{"comment":"The integral defining the dispersion function ZFD_n does not show its integration limits explicitly; please state them for completeness.","section":"Eq. (9)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript sits between plasma physics and condensed-matter theory. The analytic derivation and numerical checks are solid within the stated model, but the materials-science framing makes the missing collisionality estimate important; an editorial decision may benefit from a referee with condensed-matter transport expertise to judge whether the growth rate can plausibly exceed momentum relaxation in the cited anisotropic materials."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news here is that the 2D Weibel instability, driven by an anisotropic Fermi surface rather than a temperature anisotropy, works analytically for both linear and quadratic bands. The growth rate vanishes as q→0, there is a finite q_crit, and the unstable modes have real kappa_z—so they stay confined to the material plane, unlike stable TE waves. The linear-stage simulations match the theory, and the nonlinear runs show filamentary magnetic fields that break up turbulently. That is a solid, honest piece of kinetic theory, and I agree with the reader that no circularity is evident: the growth rates are derived, not fitted.\n\nWhere I would push back is on the severity of the collisionless assumption. The stress-test note gets this exactly right: for a Fermi surface that is anisotropic in equilibrium, collisions do not relax the distribution to isotropic—they relax to the anisotropic equilibrium. So the problem is not that the anisotropy cannot be maintained. The real gap is quantitative: the paper never compares the collision/damping rate ν to the maximum growth rate Γ_max for any candidate material. If ν > Γ_max, the mode is overdamped. The authors mention collisions but do not give numbers, and they do not provide code or data for independent checks. That is a substantive gap for the material-level claim, though not a fatal one for the idealized model.\n\nThe derivation of the secular equation (Eq. 7) is also compressed; a referee will want to see the intermediate steps, especially for the linear-band case with the elliptic integrals. The nonlinear saturation is simulation-only, and the paper itself concedes that stable magnetic domains were not achieved.\n\nOverall: this is a serious paper, worth a real referee. The central mechanism is plausible and the analytic framework is reusable, even if the experimental relevance remains unquantified. The right request is for the authors to add a collision-damping estimate for a concrete material (e.g., trilayer graphene), expand the secular-equation derivation, and consider releasing the Quetzal code or at least a detailed benchmark. I would not cite it in the next year for any quantitative claim, but I would cite it for the analytical dispersion relations once the damping question is addressed.","headline":"A credible analytic derivation of a 2D solid-state Weibel-type instability from Fermi-surface anisotropy, with the main open question being whether real materials can grow the mode before collisions damp it.","tokens_in":6943,"tokens_out":1578,"would_cite":false,"duration_ms":21538,"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":"An anisotropic Fermi surface in a two-dimensional conductor makes transverse-electric waves grow as a purely magnetic instability, generating confined out-of-plane magnetic fields.","keywords":["two-dimensional materials","Fermi surface anisotropy","Weibel instability","transverse electric mode","Vlasov equation","kinetic simulation","magnetic field generation"],"falsifier":"Measure the out-of-plane magnetic-field spectrum in a clean 2D conductor with an anisotropic Fermi surface while keeping the electron distribution anisotropic: the claim predicts purely growing (zero real frequency) magnetic fluctuations for wavenumbers below $q_{\\mathrm{crit}}$, confined to the plane. Observing only damped or oscillating modes, or growth rates far below the collisionless prediction, would falsify the central claim.","tokens_in":5977,"feed_emoji":"🧲","tokens_out":5942,"duration_ms":65477,"temperature":0.7,"pith_summary":"This paper claims that an anisotropy of the Fermi surface in a two-dimensional conductor turns the normally unsupported transverse-electric electromagnetic mode into a purely growing instability. Using a collisionless kinetic (Vlasov) description coupled to Maxwell's equations in a 2D geometry, the authors derive implicit dispersion relations for the growth rate in both quadratic and linear low-energy bands. The instability, identified as the two-dimensional version of the Weibel instability, grows with zero real frequency, is confined to the material plane, and produces structured out-of-plane magnetic fields. Fully kinetic nonlinear simulations show these magnetic structures form and then break down turbulently at saturation. If correct, this gives a mechanism for spontaneously generating tunable magnetic domains in materials with anisotropic Fermi surfaces.","feed_headline":"Anisotropic Fermi surfaces trigger magnetic fields in 2D materials","feed_subtitle":"Purely growing, confined transverse-electric waves appear when unequal chemical potentials distort the Fermi surface.","key_machinery":"The load-bearing object is an anisotropic Fermi-Dirac equilibrium with an elliptical Fermi surface, set by a chemical-potential anisotropy ratio $\\gamma$, combined with the collisionless Vlasov equation. The argument linearizes Vlasov around this equilibrium and couples the perturbed distribution to a modified Maxwell system obtained by Fourier transforming in the plane and imposing boundary conditions at the material plane; the central secular equation (7) then produces the implicit dispersion relations (13) for quadratic bands and (15) for linear bands. Quadratic bands involve the two-dimensional Fermi-statistics plasma dispersion function $Z^{\\mathrm{FD}}_1(\\zeta)$, while linear bands require an elliptic integral $E_0(\\gamma)$ because the velocity is independent of momentum magnitude. The same machinery yields a critical wavenumber $q_{\\mathrm{crit}}$ (Eq. 17) that bounds the unstable band.","core_discovery":"The paper's central claim is that an elliptical distortion of the Fermi surface, modeled by an anisotropy ratio $\\gamma = \\mu_y/\\mu_x$ between chemical potentials along the two in-plane directions, destabilizes transverse-electric waves in a 2D material. For both quadratic bands and linear bands, the linearized Vlasov-Maxwell system yields modes whose frequency is purely imaginary, $\\mathrm{Re}\\,\\omega = 0$, meaning they grow without oscillating, with growth rates given by implicit relations (13) and (15). Unstable modes exist only for wavenumbers $0 < q < q_{\\mathrm{crit}}$, and because $\\kappa_z = \\sqrt{q^2 + \\mathrm{Im}(\\omega)^2/c^2}$ remains real, the fields are evanescent and confined along the perpendicular direction, in contrast to stable TE waves in these systems, which radiate away. Nonlinear simulations confirm the linear growth and show that the instability saturates through a secondary longitudinal instability that disrupts the coherent magnetic filaments, generating an incoherent low-amplitude out-of-plane magnetic field.","pith_inferences":["If a short pump could transiently impose momentum-space anisotropy faster than scattering relaxes it, the same instability could act as an ultrafast source of confined magnetic-field textures even in materials whose equilibrium Fermi surface is isotropic.","Applying the same kinetic formalism to tilted or warped Fermi surfaces would likely replace the elliptic integral $E_0(\\gamma)$ with a band-geometry-dependent factor, offering a way to engineer the unstable wavenumber band through band-structure design.","Adding a collision operator to the Vlasov equation would set a concrete threshold: the instability should survive only when the growth rate exceeds the momentum relaxation rate, giving an experimental criterion in terms of carrier mobility and sample temperature.","The turbulent breakdown of magnetic filaments at saturation resembles magnetic turbulence seen in collisionless plasma settings, suggesting that 2D materials with driven anisotropy could serve as a compact laboratory analog of that dynamics."],"forward_implications":["In any 2D conductor with a persistent elliptical Fermi-surface anisotropy, transverse-electric waves should spontaneously grow rather than radiate, generating out-of-plane magnetic-field filaments aligned with the low-chemical-potential direction.","The unstable band is finite: only wavenumbers below $q_{\\mathrm{crit}}$ grow, with the growth-rate maximum near the middle of that band, so the emergent magnetic structure has a characteristic length scale set by density, band parameters, and the anisotropy ratio.","The instability saturates through a secondary longitudinal instability, so coherent magnetic domains are transient; stabilizing them into permanent magnetic domains remains an open problem.","Because the chemical potential and material parameters set the growth rate and spectral range, the effect may be gate-tunable in suitably anisotropic 2D conductors."],"supporting_citations":[{"why":"Supplies the original Weibel instability mechanism that the paper adapts to a two-dimensional solid-state setting.","marker":"[14]"},{"why":"Provides the foundational analysis of Weibel-mode growth that the dispersion-relation treatment extends.","marker":"[15]"},{"why":"Gives trilayer graphene with a triangular Fermi surface as a concrete material realizing the anisotropic Fermi-surface condition.","marker":"[16]"},{"why":"Establishes the baseline result that transverse-electric waves are not supported in 2D materials, which the instability overturns.","marker":"[10]"},{"why":"Supplies the Wigner/Moyal kinetic equation from which the collisionless Vlasov model is obtained.","marker":"[13]"}],"fun_headline_variants":["Fermi anisotropy yields purely growing magnetic modes in 2D","Elliptical Fermi surface seeds magnetic instability in 2D","2D materials: Fermi anisotropy drives magnetic fields","Fermi surface anisotropy drives confined magnetic growth"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation assumes the electron distribution stays anisotropic with no collisions; if ordinary scattering relaxes the Fermi-surface anisotropy faster than the instability can grow, the predicted magnetic fields will not appear in a real material.","fun_headline_variants_meta":{"raw":{"variants":["Fermi anisotropy yields purely growing magnetic modes in 2D","Elliptical Fermi surface seeds magnetic instability in 2D","2D materials: Fermi anisotropy drives magnetic fields","Fermi surface anisotropy drives confined magnetic growth"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000637,"raw_usage":{"total_tokens":2885,"prompt_tokens":844,"completion_tokens":2041,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":460,"completion_tokens_details":{"reasoning_tokens":1987}},"tokens_in":460,"tokens_out":2041,"duration_ms":17908,"temperature":1.0,"reasoning_tokens":1987,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:24:40.854639+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the out-of-plane magnetic-field spectrum in a clean 2D conductor with an anisotropic Fermi surface while keeping the electron distribution anisotropic: the claim predicts purely growing (zero real frequency) magnetic fluctuations for wavenumbers below $q_{\\mathrm{crit}}$, confined to the plane. Observing only damped or oscillating modes, or growth rates far below the collisionless prediction, would falsify the central claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the foundational analysis of Weibel-mode growth that the dispersion-relation treatment extends."},{"cited_title":"Qi and A","cited_arxiv_id":null,"evidence_quote":"Gives trilayer graphene with a triangular Fermi surface as a concrete material realizing the anisotropic Fermi-surface condition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the baseline result that transverse-electric waves are not supported in 2D materials, which the instability overturns."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Wigner/Moyal kinetic equation from which the collisionless Vlasov model is obtained."}],"review_version":1}