{"id":"b4bdd3fe-8dfb-4b3a-92f0-7ab5d1a4196d","arxiv_id":"2508.09067","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"A new gyrokinetic scheme, GYRAZE, shows the critical magnetic field angle for a monotonic sheath rises with the electron gyroradius to Debye length ratio, but stays below typical divertor angles.","lead":"This paper models the thin electric layer where a magnetised plasma meets an absorbing wall at a very shallow magnetic field angle, and derives when a simple monotonic voltage profile can exist. Its code GYRAZE predicts the range of angles and the ion and electron distributions arriving at the wall, which matters for fusion divertor design.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Critical-angle result depends on an assumed monotonic potential; without evidence that physical sheaths are monotonic below divertor angles, the claim's applicability is unverified.","rationale":"The reader's weakest_assumption identifies the same concern: the monotonic potential profile is assumed, and if the true physical sheath is non-monotonic at grazing incidence, the critical-angle result does not apply. This is exactly the load-bearing vulnerability of the paper's central claim. The abstract itself acknowledges the assumption, so this is not a hidden internal contradiction, but a limitation that affects the interpretation and practical relevance. The specific worry is that the numerical scheme may restrict the solution space a priori, making alpha_c an artifact rather than a physical threshold. The concrete test I propose directly probes whether non-monotonic solutions exist below alpha_c; if they do, the paper's comparison to divertor angles is misleading. I also considered the asymptotic ordering lambda_D/rho_S -> 0, but that is a standard boundary-layer limit and less likely to change the qualitative behavior; the monotonicity constraint is more central because it changes the solution space. Since the full manuscript is not available for inspection, the reader's UNVERDICTED verdict remains appropriate, and my concern does not move it to a different category; it strengthens the need to see the derivation and any validation. Thus the verdict should remain UNCHANGED. Agreement with the reader is 'agree' because we both focus on the monotonic assumption as the weakest point.","tokens_in":748,"tokens_out":5173,"duration_ms":57369,"concrete_test":"Run a fully kinetic particle-in-cell (PIC) simulation in 1D-2V for the same parameters as a GYRAZE case (grazing incidence, finite gamma = rho_e/lambda_D, small but finite lambda_D/rho_S) at a magnetic angle below the predicted alpha_c, without imposing any monotonicity constraint. If the PIC simulation converges to a steady-state non-monotonic potential profile, the critical angle is not a physical existence threshold and the monotonic assumption is the cause. If no non-monotonic steady solution exists and all converged profiles are monotonic, the assumption is empirically validated for that parameter set.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that, in the grazing-incidence limit lambda_D/rho_S -> 0 with finite gamma = rho_e/lambda_D, a monotonic electrostatic potential profile exists only above a critical angle alpha_c(gamma), and alpha_c remains below typical divertor magnetic field angles. The abstract explicitly states: 'A monotonic electrostatic potential profile, assumed in this work.' This assumption is load-bearing because the entire existence threshold is derived within the restricted class of monotonic profiles. If the GYRAZE iterative scheme enforces monotonicity, it cannot represent non-monotonic or oscillatory potential structures that may be physical in the grazing-incidence regime. Such non-monotonic sheaths are known in other kinetic sheath treatments and could exist for alpha < alpha_c. If so, alpha_c is not a physical boundary for sheath existence but an artifact of the constraint, and the comparison to divertor angles ('still typically smaller than the magnetic field angle') does not support the implied relevance to fusion devices. No evidence is provided in the abstract that the physical sheath in this regime is indeed monotonic; the assumption is simply stated. This is the weakest point in the logical chain from assumptions to the practical conclusion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a magnetized plasma in contact with an absorbing planar wall at grazing magnetic incidence (α ≪ 1), using a gyrokinetic treatment that retains the finite ratio γ = ρ_e/λ_D. Building on the authors' earlier iterative scheme [2,3], it proposes a code called GYRAZE to solve simultaneously for the quasineutral magnetic presheath and the non-neutral Debye sheath in the asymptotic limit λ_D/ρ_S → 0. The central claim is that, within the class of monotonic electrostatic potential profiles, a steady solution exists only for magnetic field angles above a critical angle, and that this critical angle increases with γ while remaining typically smaller than the magnetic field angle at divertor targets of a fusion device.","tokens_in":1066,"tokens_out":2858,"duration_ms":36025,"significance":"If fully substantiated, the result would be relevant to modeling magnetic presheaths and Debye sheaths at grazing incidence, an important regime for divertor physics. The explicit focus on finite γ and the simultaneous treatment of both sheath scales is timely. The paper introduces a named code and a testable prediction (a critical-angle threshold), which are strengths. However, at the level of the abstract alone, the claim is conditional on an assumed monotonic profile and on an asymptotic limit whose validity is not demonstrated. The significance therefore hinges on whether the authors can justify these assumptions and provide numerical or analytic evidence.","major_comments":[{"comment":"The central existence threshold is explicitly conditional: 'A monotonic electrostatic potential profile, assumed in this work.' The critical angle is a threshold within the restricted class of monotonic profiles, not a proven property of physical sheaths. If non-monotonic or oscillatory potential structures exist for α < α_c, as is possible in kinetic sheath treatments, then α_c would not mark the disappearance of sheaths but rather the limit of the assumed ansatz. The abstract gives no evidence that sheaths in this regime are indeed monotonic, nor any reason to expect this. This is a load-bearing gap.","section":"Abstract"},{"comment":"The asymptotic limit λ_D/ρ_S → 0 is used to separate the Debye sheath and the magnetic presheath, while γ = ρ_e/λ_D is retained as a finite parameter. No equations or convergence checks are shown to demonstrate that this limit is well defined or that the two-scale expansion does not introduce uncontrolled errors. For real divertor conditions λ_D/ρ_S is finite and small, not zero; the abstract provides no quantitative estimate of the error incurred by this idealization when comparing to divertor angles.","section":"Abstract"},{"comment":"The statement that the critical angle 'is still typically smaller than the magnetic field angle at divertor targets' is an unquantified assertion. The abstract reports no values of α_c, γ, or the relevant divertor angles, and no comparison data are visible. This claim is central to the paper's practical relevance and must be supported with explicit numbers and/or a figure.","section":"Abstract"},{"comment":"The paper claims that GYRAZE 'simultaneously solves' for both the presheath and the Debye sheath, but no numerical method, discretization, convergence tests, or benchmark comparisons are provided. Without such details, the correctness of the iterative scheme and the reliability of the reported critical-angle behavior cannot be assessed. At minimum, the authors should state the governing equations and the iterative procedure, and show that the solution converges in λ_D/ρ_S and in numerical resolution.","section":"Abstract"}],"minor_comments":[{"comment":"The notation ρ_S is not defined; it should be stated as the ion sound Larmor radius (or equivalent) for clarity.","section":"Abstract"},{"comment":"The phrase 'monotonic electrostatic potential profile' should specify the coordinate with respect to which monotonicity is assumed (e.g., distance from the wall) and whether it includes both the presheath and Debye sheath regions.","section":"Abstract"},{"comment":"The reference [3] is cited for the existence of a critical angle; it would help to state what was established there and what is new in this work regarding the γ dependence.","section":"Abstract"},{"comment":"The abstract mentions 'energy-angle distribution of ions at the wall and the velocity distributions of electrons reflected by the wall' but does not explain how these are used to determine the critical angle; a sentence connecting these outputs to the threshold would improve readability.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The referee received only the abstract; the full text is not available in the submission. The report is based solely on the abstract. The central claim is clear but conditional on an explicit monotonicity assumption and an asymptotic limit. If the full manuscript contains derivations, convergence studies, and comparisons, the authors should be able to address the major comments. If not, the paper is not yet ready for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The abstract gives a clear picture: a new combined presheath-sheath solver (GYRAZE) that works in the grazing-incidence limit, retaining finite electron gyroradius relative to Debye length. The genuinely new result is that the critical angle for a monotonic potential increases with gamma = rho_e / lambda_D, and it remains below typical divertor field angles. That is a sharp, testable claim, and it is directly relevant to tokamak edge modeling.\n\nI think this is a solid extension, not a conceptual breakthrough. It builds on the authors' own iterative scheme, which is fine—the citation pattern is not a problem when the prior work is theirs and they are extending it honestly. The abstract is also refreshingly transparent: it states outright that a monotonic profile is assumed. That is exactly the load-bearing assumption, and they do not hide it.\n\nThat said, the stress-test worry is legitimate. If GYRAZE enforces monotonicity in its iteration, then the critical angle is a threshold within a restricted class of solutions, not necessarily a physical boundary. Non-monotonic sheaths could exist below alpha_c, and then the comparison to divertor angles would be misleading. The abstract does not show evidence that physical sheaths are monotonic in this regime. This is a real soft spot, but it is the kind of thing a good referee can probe. The asymptotic lambda_D/rho_S -> 0 is standard and not a concern in itself.\n\nI cannot judge soundness beyond the abstract because there is no full text here: no equations, no convergence checks, no comparisons. But what I see is internally coherent, and the authors are not fitting data or hiding free parameters. The monotonicity issue is the main thing I would want resolved—does the solver actually allow non-monotonic profiles, and if not, what is the justification for assuming they are absent?\n\nWho is this for? Plasma sheath and divertor physicists, especially anyone building reduced models for SOL/divertor codes. It deserves a serious referee. I would send it out, with a request to focus on the monotonic assumption and to compare against independent kinetic simulations at grazing incidence.\n\nBottom line: worth engaging with, worth citing once the full derivation is verified, and the monotonic-profile caveat should be taken seriously rather than treated as a deal-breaker.","headline":"A transparent, incremental but useful extension of the authors' own presheath-sheath solver; the new critical-angle dependence on gamma is provocative, but the monotonic-profile assumption needs scrutiny in the full paper.","tokens_in":1462,"tokens_out":1416,"would_cite":true,"duration_ms":16711,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.40.Kh"],"model":"deepseek-v4-flash","headline":"This paper establishes that monotonic plasma sheaths near a wall exist only above a critical magnetic-field angle that grows with the ratio of electron gyroradius to Debye length, and that this angle typically stays below fusion divertor fi","keywords":["magnetised plasma sheath","grazing incidence","Debye sheath","magnetic presheath","critical angle","monotonic potential","gyrokinetic","divertor"],"falsifier":"A kinetic simulation that does not impose monotonicity, run at fixed $\\gamma$ with $\\alpha$ below the predicted $\\alpha_{\\rm c}^*(\\gamma)$, should fail to find any steady monotonic potential; finding one would disprove the predicted boundary, while observing a non-monotonic steady state would support it.","tokens_in":732,"feed_emoji":"🧲","tokens_out":7882,"duration_ms":75123,"temperature":0.7,"pith_summary":"This paper shows that a monotonic electrostatic potential in a magnetised plasma sheath near an absorbing wall exists only when the magnetic field incidence angle exceeds a critical value, and that this critical angle grows with the ratio of electron gyroradius to Debye length, $\\gamma = \\rho_{\\rm e}/\\lambda_{\\rm D}$. The authors solve simultaneously for the quasineutral magnetic presheath and the thin non-neutral Debye sheath in the grazing-incidence limit, and their scheme yields the ion energy-angle distribution at the wall and the reflected-electron velocity distribution. The result matters for fusion experiments because magnetic field lines hit divertor targets at shallow angles; the predicted critical angle remains below typical divertor angles, but the margin narrows as $\\gamma$ grows.","feed_headline":"Plasma sheath turns non-monotonic below a critical angle","feed_subtitle":"New solver shows the cutoff rises with electron gyroradius relative to Debye length, yet stays below divertor angles.","key_machinery":"The central machinery is an iterative solution scheme that matches the quasineutral magnetic presheath (width ~$\\rho_{\\rm S}$) with the non-neutral Debye sheath (width ~$\\lambda_{\\rm D}$) in the double limit $\\alpha \\ll 1$ and $\\lambda_{\\rm D}/\\rho_{\\rm S}\\to 0$, while keeping $\\gamma = \\rho_{\\rm e}/\\lambda_{\\rm D}$ finite. The scheme exploits the scale separation to impose quasineutrality in the presheath and space-charge balance in the sheath, and it resolves the full energy-angle dependence of particle orbits. The output is the existence domain for monotonic solutions as a function of $\\gamma$ and wall potential, from which the critical angle is read as the minimum angle for steady monoto","core_discovery":"The central claim is that the existence of a monotonic potential profile in the steady-state sheath is bounded below by a critical angle $\\alpha_{\\rm c}^*(\\gamma)$. In the idealized limit $\\lambda_{\\rm D}/\\rho_{\\rm S}\\to 0$, the presheath and Debye sheath are solved together, retaining finite $\\gamma$ through a gyrokinetic description valid at grazing incidence. As $\\gamma$ increases, $\\alpha_{\\rm c}^*$ increases substantially, so a larger electron gyroradius (relative to Debye length) makes the monotonic solution harder to sustain. Even so, $\\alpha_{\\rm c}^*$ is typically smaller than the field-line angles at fusion divertor targets, so the monotonic assumption used in many edge-plasma mode","pith_inferences":["One testable extension is to run a particle-in-cell simulation without the monotonicity constraint at $\\alpha$ just below the predicted critical angle; if a steady non-monotonic potential appears, the boundary marks a real transition, not a mathematical artifact.","The $\\gamma$ dependence suggests that the safety margin for monotonic sheaths in a fusion device could be quantified in terms of local temperature and density; a divertor design could use $\\gamma$ as a dimensionless monitor.","If non-monotonic profiles exist just below $\\alpha_{\\rm c}^*$, they may trap electrons and alter the heat flux to the wall, making the critical angle relevant to power-load predictions."],"forward_implications":["In edge-plasma models that assume a monotonic sheath, the angle between the magnetic field and the wall must remain above $\\alpha_{\\rm c}^*(\\gamma)$; below this, a different treatment is needed.","Because $\\alpha_{\\rm c}^*$ grows with $\\gamma$, devices with hotter electrons or smaller Debye lengths have a more restrictive range of angles for which the standard monotonic-sheath picture holds.","The GYRAZE output provides ion and electron velocity distributions at the wall for a given wall potential, which can be used directly as boundary conditions for fluid or kinetic simulations of the scrape-off layer.","The same matching procedure can be applied to determine how the plasma-wall interaction changes when the monotonicity assumption is relaxed."],"supporting_citations":[],"fun_headline_variants":["Critical angle rises with gyroradius for monotonic plasma sheaths","Thicker electron orbits shrink range of monotonic wall sheaths","Grazing-field sheath monotonicity: a sharper cutoff with γ","New solver GYRAZE maps when plasma wall sheaths stay monotonic","Monotonic sheaths vanish below a γ-dependent critical angle"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The analysis assumes the electrostatic potential profile is monotonic; if a real grazing-incidence sheath can be non-monotonic, the critical-angle boundary does not apply.","fun_headline_variants_meta":{"raw":{"variants":["Critical angle rises with gyroradius for monotonic plasma sheaths","Thicker electron orbits shrink range of monotonic wall sheaths","Grazing-field sheath monotonicity: a sharper cutoff with γ","New solver GYRAZE maps when plasma wall sheaths stay monotonic","Monotonic sheaths vanish below a γ-dependent critical angle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000354,"raw_usage":{"total_tokens":1807,"prompt_tokens":832,"completion_tokens":975,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":896}},"tokens_in":576,"tokens_out":975,"duration_ms":8074,"temperature":1.0,"reasoning_tokens":896,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T21:12:31.583299+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A kinetic simulation that does not impose monotonicity, run at fixed $\\gamma$ with $\\alpha$ below the predicted $\\alpha_{\\rm c}^*(\\gamma)$, should fail to find any steady monotonic potential; finding one would disprove the predicted boundary, while observing a non-monotonic steady state would support it.","supporting_citations":[],"review_version":1}