{"id":"847b41de-cc51-4403-a19f-b2945cda16ac","arxiv_id":"2603.16515","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"An added resistive-ballooning transport channel lets the Saarelma-Connor model reproduce EDA H-mode pedestal densities on C-Mod up to 3×10^20 m^-3 and lowers predicted SPARC pedestal density by ~20%.","lead":"This paper studies why the density pedestal behaves differently in two tokamak edge regimes (ELMy versus EDA H-modes) on Alcator C-Mod and extends a pedestal density model to reproduce the EDA regime by adding a resistive-ballooning transport channel. The work validates the model to higher densities and gives first pedestal-density projections for the SPARC tokamak.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"D_RBM may be an artifact of an overestimated KN1D neutral boundary: Section 4.2 leaves unresolved whether the high-density overprediction is a neutral-source error or missing transport, so the central RBM channel is not yet established.","rationale":"Reader identified the same key uncertainty, and I agree. The chain 'neutral boundary -> overprediction -> D_RBM' is exactly where the central claim is least secure, because the paper itself states the ambiguity. I did not find an internal inconsistency in the RBM scaling; the alternative forms (Eq. 13, 14) are plausible and the authors are appropriately tentative. The main gap is that no independent constraint on n_sep0 is available for these discharges, and the chosen constants come from earlier C-Mod studies, so the model has enough flexibility to fit the high-density residual. This does not invalidate the paper: the ELMy validation, EPED comparison, and SPARC sensitivity analysis remain valuable. It does mean the claim of an RBM-driven channel needs a decisive boundary-condition check before being accepted as predictive. Verdict remains CONDITIONAL.","tokens_in":42021,"tokens_out":8522,"duration_ms":90550,"concrete_test":"Rerun KN1D with a finer scan (≥30 points) over p_OMP0=0.03–0.6 mTorr, including a limiter-shadow density falloff and T_i≠T_e, and refit n_sep0 vs p_OMP0. Then recompute the Saarelma-Connor predictions (same transport settings as Fig. 10) with D_RBM=0. If the EDA points with n_ped>2×10^20 m^-3 no longer systematically overpredict experiment, the RBM channel is not required; if the overprediction persists, the D_RBM residual is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that D_RBM is required for EDA pedestal densities depends on the KN1D mapping n_sep0=38.5 p_OMP0 used in Section 4.2/Appendix A. The paper admits: 'Two possibilities exist – either the neutral source is overestimated or the plasma transport is underestimated.' The present analysis does not decide between them: the KN1D mapping is not an independent measurement. It is built from ten moving-average profiles, assumes T_e=T_i (explicitly flagged in Appendix A as an assumption to relax), uses a linear fit in Fig. A3 despite a hint of saturation at high p_OMP0, and applies a limiter-shadow plasma density that the authors note may overestimate particle content. If a more realistic neutral boundary lowers n_sep0 for EDA points, the standard model (D_RBM=0) may already match the high-density data; D_RBM would then be absorbing a boundary-condition bias. The fluctuation measurements do not fill this gap: the coherent QCM amplitude saturates and even weakens at the highest n_ped (Section 3.1), while D_RBM is largest there. So the 'good agreement' in Fig. 10 is a postdiction conditional on an unverified boundary condition, not an independent confirmation of RBM-driven transport.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes a set of Alcator C-Mod discharges spanning ELMy and EDA H-modes, using high-resolution Thomson scattering profiles and PCI fluctuation measurements. The authors show that the pedestal density is insensitive to neutral fueling in the EDA regime, while it is fueling-sensitive in the ELMy regime. They validate the Saarelma-Connor pedestal density prediction model on this dataset and find that the model overpredicts the EDA pedestal density at high density. To correct this, they add an ad hoc resistive-ballooning-mode (RBM) particle transport channel, D_RBM, with two alternative forms scaling with α_t and with 1/(k_RBM^2 q_cyl), and report improved agreement up to n_ped = 3×10^20 m^-3. They also perform EPED scans at different n_sep/n_ped ratios, compare to experiment, and make initial SPARC pedestal density predictions for an ELMy and an EDA/QCE-like scenario, finding that D_RBM lowers the predicted SPARC pedestal density by about 20%.","tokens_in":42423,"tokens_out":4265,"duration_ms":42993,"significance":"If the central claim holds, this paper would extend the validated range of a predictive pedestal density model into high-density, ELM-free regimes relevant to ITER and SPARC, and would identify RBM-driven particle transport as the mechanism setting the pedestal density in EDA/QCE H-modes. The manuscript has several concrete strengths: it uses a well-characterized C-Mod dataset with two independent Thomson-scattering fitting approaches; it tracks the QCM amplitude systematically across regimes; it validates the Saarelma-Connor model on a new device and regime; and it is transparent about the model's free parameters and limitations. The EPED scans at fixed n_sep/n_ped are a useful sensitivity study. However, the evidence for the new D_RBM channel is not yet independent of the neutral-boundary model used to infer the separatrix neutral density, and several transport coefficients are hand-tuned. As a result, the main claim is plausible but not fully established.","major_comments":[{"comment":"The conclusion that an RBM-driven transport channel is required in EDA H-modes is load-bearing and rests on the KN1D boundary mapping n_sep0[10^15 m^-3] = 38.5 p_OMP0[mTorr]. The paper itself states 'Two possibilities exist – either the neutral source is overestimated or the plasma transport is underestimated.' The present analysis does not resolve this ambiguity: the mapping is fit to moving-average profiles from the same dataset, assumes T_e=T_i (flagged in Appendix A), uses a linear fit in Fig. A3 despite visible saturation at high p_OMP0, and applies a limiter-shadow density that may overestimate particle content. If n_sep0 is overestimated for EDA points, the standard model may already match the high-density data. Please provide sensitivity tests (e.g., T_i/T_e variation, limiter density variation, alternative neutral model) or independent neutral-density constraints to show D_RBM i","section":"4.2 / Appendix A / Eqs. (12)-(14)"},{"comment":"The improved transport settings use C_KBM=0.01 and α_crit=3, whereas α_crit=2 was found suitable for other devices. The motivation for α_crit=3 from an average separatrix α_c=2.6 is suggestive but not a derivation; α_crit is a free parameter for the KBM onset inside the pedestal. Fig. 14 shows that lowering α_crit to 2 changes SPARC n_ped by roughly 15%, so the model is sensitive. Please report a systematic parameter scan (C_KBM, α_crit, (D/χ)_TG) with a quantitative goodness-of-fit metric and uncertainty estimates; otherwise the claimed validation is vulnerable to overfitting.","section":"4.1, Eq. (9) and Fig. 9"},{"comment":"The fluctuation data do not currently corroborate the D_RBM term where it matters most. The QCM amplitude B saturates and weakens for n_ped > 2.5×10^20 m^-3 (Fig. 6, right), yet D_RBM is largest at the highest n_ped. The background amplitude A continues to grow, but A is not directly linked to the radial particle transport coefficient D_RBM. Please either connect the fluctuation measurements quantitatively to the proposed transport channel (e.g., through a mixing-length estimate) or temper the claim that the RBM channel is independently supported by the PCI data.","section":"3.1 / Fig. 6 and Fig. 10"},{"comment":"The SPARC high-density predictions use very crude inputs—T_EDA = 0.5 T_PRD, n_EDA = 1.5 n_PRD, and ad hoc width adjustments—and the result is highly sensitive to the choice of n_sep0 and whether D_RBM is included (n_ped ranges from 5.1 to 8.8×10^20 m^-3 across the explored settings). The statement that the predictions are 'consistent with assumptions used in previous EPED modeling' should be qualified with these sensitivities; as it stands, the SPARC section is illustrative rather than predictive.","section":"6, Fig. 16"}],"minor_comments":[{"comment":"The linear fit n_sep0 = 38.5 p_OMP0 appears to be strongly influenced by the highest-pressure point, and the text notes a possible saturation. Show residuals and fit uncertainty, or use a saturating form and justify the linear choice.","section":"Fig. A3"},{"comment":"Define all symbols (k_RBM, q_cyl) and give units. Currently C_k_RBM is given only numerically, and the physical dimensions of the expression are not stated.","section":"Eq. (14)"},{"comment":"The transition at p_OMP0 ≈ 0.1 mTorr is central to the regime classification. Provide an uncertainty estimate for the pressure measurement and for how the transition value is determined.","section":"Section 2.2"},{"comment":"The sentence 'C_RBM = 0.039, taken empirically from the dataset in [33]' is ambiguous—was the coefficient calibrated on the same run day or on an independent dataset? Please clarify to avoid circularity concerns.","section":"Section 4.2"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern lands: the central RBM conclusion is not yet established because the KN1D boundary mapping is not independent of the validation dataset, and the paper explicitly leaves open the neutral-source vs. transport ambiguity. However, the authors are unusually transparent about this limitation, and the issue is addressable with additional sensitivity analysis and uncertainty quantification. This is not a reject-level error, but the main claim needs stronger support before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid experimental paper that does something genuinely new. It validates the Saarelma-Connor pedestal density model on C-Mod, extends it into the EDA regime, and adds a D_RBM transport term that makes the model match EDA densities up to 3e20 m^-3. The EPED scans at variable n_sep/n_ped are a useful addition, and the SPARC projections, while rough, are the first for an ELM-free scenario. The fluctuation analysis is careful, and the paper is honest about its own limitations.\n\nThe main soft spot is exactly where the stress-test lands: the D_RBM term is not independently established. The KN1D boundary mapping n_sep0 = 38.5 p_OMP0 is built from characteristic profiles from this same dataset, assumes Te=Ti, and is fit through a linear relation that shows a hint of saturation. The paper itself says 'Two possibilities exist – either the neutral source is overestimated or the plasma transport is underestimated.' That is not a minor caveat; it is the load-bearing uncertainty. If the neutral source is overestimated, the standard model without D_RBM might already match the EDA data, and D_RBM would compensate for a boundary-condition bias. The fluctuation data do not resolve this: the QCM amplitude saturates, and even weakens, at the highest n_ped, where D_RBM is largest. So the model agreement in Fig. 10 is best read as a postdiction that depends on an unverified boundary condition, not a confirmation of RBM-driven transport.\n\nThere are smaller issues. The transport coefficients (alpha_crit=3 vs 2 elsewhere, C_KBM=0.01) are hand-tuned to improve agreement. There are no error bars on the fluctuation amplitudes. The SPARC predictions rely on very crude profile estimates (T_e divided by two, n_e multiplied by 1.5) and change quite a bit with the neutral boundary condition. None of this kills the paper, but it means the central quantitative claim needs more work.\n\nWho gets value: pedestal modelers, C-Mod experimentalists, and anyone doing SPARC edge scenario work. If I were refereeing, I would ask for: (1) a sensitivity study that varies the KN1D mapping within plausible uncertainty, especially the Te=Ti assumption and the limiter-shadow density; (2) a cleaner test of D_RBM, perhaps by comparing discharges where the neutral boundary is independently measured; (3) error bars on B. Even with those, I would send this to peer review. The dataset is valuable, the model extension is a reasonable hypothesis, and the authors have been transparent about what is not known. This is the kind of paper that should be published with the uncertainty stated, not desk-rejected.","headline":"A careful C-Mod study that plausibly extends pedestal density prediction to EDA H-modes, but the new RBM transport channel rests on a boundary condition that could be absorbing a neutral-source error.","tokens_in":42946,"tokens_out":2574,"would_cite":true,"duration_ms":24602,"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":"The pedestal density of the ELM-free EDA H-mode is set by resistive-ballooning-mode particle transport, not by neutral fueling.","keywords":["pedestal density","EDA H-mode","ELMy H-mode","resistive ballooning mode","particle transport","neutral fueling","quasi-coherent mode","SPARC pedestal prediction"],"falsifier":"Measure the actual neutral ionization source in the pedestal (for example, with Ly-alpha emission or a calibrated neutral-density diagnostic) across the ELMy-to-EDA transition and compare it with the model's assumed source. If the separatrix neutral density is lower than the mapped value by enough to remove the overprediction, then D_RBM as added is an artifact; if the neutral source matches, the RBM transport channel is required. A gas-puff modulation experiment that varies neutral pressure while holding the separatrix density fixed could separate source effects from transport effects directl","tokens_in":41933,"feed_emoji":"🧲","tokens_out":7767,"duration_ms":70596,"temperature":0.7,"pith_summary":"Using a set of Alcator C-Mod discharges that cross the ELMy-to-EDA H-mode transition, the paper argues that the two regimes are governed by different density-control mechanisms. In the ELMy H-mode, the pedestal density climbs with neutral pressure and responds to fueling; in the EDA H-mode, it saturates and becomes insensitive to the neutral source, implying that a turbulent transport channel limits the density gradient. The paper identifies that channel as resistive ballooning mode (RBM) transport and shows that adding a diffusion term D_RBM, scaling with the collisionality parameter α_t and inversely with k_RBM^2 q_cyl, lets a recent pedestal-density prediction model reproduce EDA pedestals up to 3×10^20 m^-3. If correct, this would mean that in high-density, ELM-free regimes the pedestal density is transport-limited, not fueling-limited, a result with direct consequences for designing ELM-free operation in next-step devices such as SPARC.","feed_headline":"RBM transport, not fueling, sets the high-density pedestal","feed_subtitle":"RBM diffusion term matches high-density EDA pedestals and shaves 20% off SPARC density predictions.","key_machinery":"The load-bearing object is the RBM particle-diffusion coefficient D_RBM added to the pedestal transport sum D_ped = D_neo + D_KBM + D_TG + D_RBM. Two forms are tested: D_RBM = C_RBM α_t^a with a≈1.2, and D_RBM = C_k/(k_RBM^2 qhat_cyl), where α_t is a collisionality-like turbulence drive parameter, k_RBM is the characteristic resistive-ballooning wavenumber from the two-fluid model, and qhat_cyl is the cylindrical safety factor. This term provides the extra outward particle flux that clamps the density gradient at high density; the neutral boundary is supplied by a kinetic neutral simulation mapping the measured wall neutral pressure to the separatrix neutral density.","core_discovery":"The central claim is that the pedestal density in the high-density EDA H-mode is set by RBM-driven particle transport, not by neutral fueling. Experimentally, n_ped rises with neutral pressure in ELMy H-modes but saturates and becomes insensitive to the neutral source in EDAs, even as n_sep continues to grow. Fluctuation spectra show the quasi-coherent mode (QCM) strengthening across the transition, then saturating and weakening at the highest densities, while broadband fluctuation levels keep rising. The paper adds an RBM diffusivity D_RBM (either C_RBM α_t^a or C_k/(k_RBM^2 qhat_cyl)) to the pedestal transport sum and shows it reproduces EDA pedestals up to 3×10^20 m^-3. EPED scans at n_se","pith_inferences":["If the α_t/k_RBM^2 q scaling for D_RBM is universal, pedestal models for other high-density, ELM-free regimes (for example, the quasi-continuous exhaust regime on other tokamaks) may need the same term; this is testable with existing databases.","The neutral-boundary uncertainty could be reduced by comparing the kinetic neutral simulation's mapping to direct Ly-alpha measurements of the ionization source; such a test would either strengthen or remove the case for a separate RBM transport channel.","The prediction that n_ped approaches n_sep at high density in SPARC suggests that pedestal performance and divertor protection become coupled through the same RBM transport; if true, optimizing the separatrix density may be the shared lever for both.","The weakening of the QCM just before the density limit hints that the RBM channel may replace, rather than merely supplement, the kinetic ballooning mode; a fluid turbulence simulation resolving both instabilities could identify which one is active."],"forward_implications":["The density-pedestal model, extended with D_RBM, is validated for ELMy H-modes up to 2×10^20 m^-3 and for EDA H-modes up to 3×10^20 m^-3, extending the model's range to non-ELMing, high-density regimes.","In EDA H-modes, n_ped is essentially fixed by turbulent transport, so gas fueling cannot be used to raise pedestal density; density must be controlled through edge transport or plasma shape.","EPED scans show that rising n_sep/n_ped shifts the peeling-ballooning transition to lower n_ped, so a high separatrix-to-pedestal ratio (typical of EDA) makes pedestal pressure ballooning-limited at lower density.","For SPARC, including RBM transport lowers n_ped by about 20% and weakens the density gradient near the separatrix in the high-density EDA/QCE-like scenario, changing the expected edge profile for power handling and ELM avoidance.","The observed saturation of the QCM amplitude at high n_ped, while broadband fluctuations keep rising, points to additional turbulence beyond the QCM contributing to transport near the density limit."],"fun_headline_variants":["RBM transport, not neutrals, sets high-density pedestal","Pedestal density controlled by RBM transport, not fueling","RBM diffusion, not neutral sources, governs EDA pedestals","Transport beats fueling: RBM sets pedestal density","RBM transport shaves 20% off SPARC pedestal density"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The need for the extra RBM transport channel rests on the neutral boundary mapping from the kinetic neutral simulation, n_sep0[10^15 m^-3] = 38.5 p_OMP0[mTorr] with the assumption T_e=T_i; if the neutral source is overestimated, D_RBM may be compensating for a boundary-condition error rather than representing a real transport channel.","fun_headline_variants_meta":{"raw":{"variants":["RBM transport, not neutrals, sets high-density pedestal","Pedestal density controlled by RBM transport, not fueling","RBM diffusion, not neutral sources, governs EDA pedestals","Transport beats fueling: RBM sets pedestal density","RBM transport shaves 20% off SPARC pedestal density"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000284,"raw_usage":{"total_tokens":1652,"prompt_tokens":1028,"completion_tokens":624,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":772,"completion_tokens_details":{"reasoning_tokens":536}},"tokens_in":772,"tokens_out":624,"duration_ms":6621,"temperature":1.0,"reasoning_tokens":536,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T18:00:14.735939+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the actual neutral ionization source in the pedestal (for example, with Ly-alpha emission or a calibrated neutral-density diagnostic) across the ELMy-to-EDA transition and compare it with the model's assumed source. If the separatrix neutral density is lower than the mapped value by enough to remove the overprediction, then D_RBM as added is an artifact; if the neutral source matches, the RBM transport channel is required. A gas-puff modulation experiment that varies neutral pressure while holding the separatrix density fixed could separate source effects from transport effects directl","supporting_citations":[],"review_version":1}