{"id":"8a328517-7969-4a55-841a-aab4f3279707","arxiv_id":"2504.16190","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"First GRILLIX turbulence simulations reproduce detached X-point radiator states in ASDEX Upgrade and reveal a turbulent radiator with intermittent recombination cores and strong radial flows.","lead":"Using a turbulence code, the authors simulated the detached 'X-point radiator' state in the ASDEX Upgrade tokamak, matching measured temperatures and radiation patterns. The simulations show the radiating layer is highly unstable and turbulent, with implications for controlling heat exhaust and edge instabilities in future fusion devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"High-XPR height claim rests on a confounded parameter change: final High-XPR differs from Low-XPR in c_imp and T_N,min while N_div returns to baseline, so neutral-gas height modulation is not cleanly established.","rationale":"I read the paper as a credible first demonstration of XPR turbulence: the GRILLIX model is described transparently, the new implicit resistivity scheme is verified by the Method of Manufactured Solutions, and the Low-XPR and High-XPR states are compared against OMP profiles, DTS divertor profiles, and bolometry. The strongest claim is not that the XPR height is predicted from first principles, but that a turbulence code can reproduce a detached XPR and expose its turbulent structure. The reader's CONDITIONAL verdict already captures the main limitation: the front position and detached state are set by imposed parameters (c_imp, N_div, T_N,min, and a temporary neutral seed) rather than by self-consistent impurity and neutral transport. My stress-test sharpens one specific aspect of that limitation: the High-XPR case used for the most favorable validation is not a controlled test of neutral-gas height modulation, because its final state differs from Low-XPR in two parameters (c_imp halved and T_N,min raised) while N_div returns to its original value. The inference in the abstract and Section 3.3 that neutral gas is critical for modulating XPR height is therefore less secure than the rest of the paper. This does not overturn the central demonstration, but it does mean the quantitative 12 cm height and the resulting transport/fluctuation analysis should be regarded as scenario reproduction with tuning, not as predictive confirmation. A simple four-branch restart scan would settle whether the height is due to N_div, T_N,min, or c_imp. Since the reader's verdict was already CONDITIONAL, my read does not change it; it only gives a sharper reason for the condition.","tokens_in":16444,"tokens_out":5995,"duration_ms":59095,"concrete_test":"Branch four simulations from the Low-XPR state at t=1.7 ms (Sec 3.3): (A) raise N_div to 2e19 and hold it there, with c_imp=10% and T_N,min=30; (B) only increase T_N,min to 50 eV; (C) only reduce c_imp to 5%; (D) reproduce the published High-XPR recipe. If (A) alone reaches and holds a ~12 cm front, the neutral-gas modulation claim survives; if (B) or (C) alone reproduces the high front, the height is controlled by the imposed minimum neutral diffusivity or impurity fraction rather than by divertor neutral density, and the High-XPR bolometry match is a tuned result rather than a neutral-gas prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's validation of a 12 cm XPR and its conclusion that neutral gas modulates XPR height rest on an uncontrolled comparison. In Section 3.3, the High-XPR case is created from the Low-XPR state at t=1.7 ms by transiently raising N_div from 1e19 to 2e19 m^-3 and T_N,min from 30 to 50 eV, then restoring N_div to 1e19 m^-3 while simultaneously lowering c_imp from 10% to 5%. The final High-XPR state therefore differs from Low-XPR in c_imp and T_N,min, not in the supposedly decisive divertor neutral density. Because Eq. (1) makes radiated power proportional to c_imp, and Section 3.2 shows c_imp strongly controls XPR persistence (reducing it from 10% to 2% makes the structure drift below the separatrix and vanish), the lower c_imp is a plausible alternative driver of the 12 cm front. Likewise, the higher T_N,min changes neutral penetration by imposing a minimum charge-exchange diffusivity. As a result, the attribution of the high front to neutral gas is not supported by the presented simulations, and the bolometry agreement of the tuned High-XPR case is weaker evidence for predictive capability than it appears.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents first-of-kind turbulence simulations of detached X-point radiator (XPR) conditions in ASDEX Upgrade using the GRILLIX code. The plasma model is extended with a fixed-impurity-fraction radiation term, an implicit parallel resistivity scheme, and a lower bound on the neutral diffusivity. Two XPR states are reported, with radiation fronts at 5 cm and 12 cm above the X-point, radiating about 80% of the input heating power; they are compared with an attached reference case and validated against outboard-midplane, divertor, and bolometry measurements. The paper further analyzes fluctuation statistics, radial transport coefficients, the radial electric field, and E×B flows, and argues that neutral gas is critical for detachment and for modulating the XPR height, and that the observed turbulence and flows may relate to ELM suppression.","tokens_in":16697,"tokens_out":5668,"duration_ms":54753,"significance":"If the conclusions hold, this is a valuable step: it is the first global turbulence simulation of a detached XPR, and it demonstrates the feasibility of coupling drift-fluid turbulence with neutral and impurity effects in a diverted geometry. The numerical extensions are verified with the Method of Manufactured Solutions, and the paper is generally transparent about the modeling assumptions. The reported intermittent recombination cores and the large (400%) fluctuation amplitudes are interesting and potentially relevant to the ELM-suppressed XPR regime. The main caveat is that the XPR height and radiation fraction are strongly influenced by tunable parameters (c_imp, N_div, T_N,min, and a temporary neutral source), so the validation is partly a consistency check rather than an independent predictive test.","major_comments":[{"comment":"The claim that neutral gas modulates the XPR height is not supported by the presented simulations, because the High-XPR case differs from the Low-XPR case in three parameters simultaneously. The text states that N_div is transiently raised and T_N,min is increased, and then 'the impurity concentration parameter was reduced from c_imp = 10%→5%.' After N_div is restored to 1e19 m^-3, the final High-XPR state differs from Low-XPR in both T_N,min (30→50 eV) and c_imp (10%→5%). Since Eq. (1) makes p_rad proportional to c_imp, and §3.2 shows that reducing c_imp strongly affects XPR persistence and position, the lower c_imp is a plausible alternative driver of the 12 cm front. To support the stated conclusion, a control simulation that varies only N_div (or a clear statement that the height is a combined effect of all three changes) is required.","section":"§3.3 (and abstract)"},{"comment":"The bolometric validation is partially circular. The DLX diode data is scaled 'such that the experimental peak aligns approximately with the peak value obtained in the High-XPR case,' and the High-XPR state itself is obtained by tuning N_div, T_N,min, and c_imp. Moreover, the radiation fraction of 0.8 in Table 1 is a direct consequence of the chosen c_imp and the imposed neutral sources; it is not an emergent prediction. The sentence 'The agreement of the High-XPR case is remarkable' therefore overstates the evidential value: the comparison demonstrates that the model can localize a strong radiator at the tuned height, but it does not independently validate the predicted height or the 80% radiation fraction.","section":"§4.3 and Table 1"}],"minor_comments":[{"comment":"The statement that 'whenever simulation data is shown ... it has been averaged over ... the last 50 snapshots, spanning 0.1 ms' should be accompanied by a sensitivity check or a statement on statistical convergence, since 0.1 ms is short compared to the 1 ms quasi-steady interval quoted in §3.2.","section":"§3.1"},{"comment":"The temporary neutral source is a significant ad hoc element. The paper says the front persists after removal, but no quantitative diagnostic is shown to demonstrate that the quasi-steady state is unaffected by the source's memory; a brief comparison of fields before and after removal would be useful.","section":"§3.2"},{"comment":"The phrase 'though note that some overestimation should be expected due to higher radiation compared to the experiment' is vague; it would be clearer to state that the simulated heating power and radiated power are higher than the experimental values (Table 1) by specific amounts.","section":"§4.3"},{"comment":"The histograms are presented with a count axis (log scale) but not as normalized probability densities; please clarify in the caption how the counts are obtained and whether the 400% figure corresponds to the sample maximum or a percentile.","section":"Figure 10"}],"recommendation":"major_revision","confidential_remarks":"The confounded parameter change in §3.3 is the main stumbling block. The authors should be encouraged to run a control simulation with only N_div varied, or to rewrite the abstract and conclusions to attribute the height change to the combined effect. The paper is otherwise within the journal's scope and of interest to the plasma-edge community."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First turbulence-resolving simulations of an X-point radiator in a diverted tokamak; that is the real news. The GRILLIX extensions (impurity radiation, implicit resistivity, neutrals diffusivity limiter) are described concretely, the Ohm's law change is MMS-verified, and the validation against OMP/divertor profiles and bolometry is honest and mostly convincing. The turbulent front morphology, fluctuation PDFs with density fluctuations above 400%, and the E×B vortex around the X-point are new and worth taking seriously. As a capability demonstration, the paper largely succeeds.\n\nThe soft spots are concentrated in the height-modulation claim. The High-XPR case starts from Low-XPR at 1.7 ms with N_div raised from 1e19 to 2e19 m^-3 and T_N,min from 30 to 50 eV; after 0.24 ms N_div returns to 1e19 m^-3 while c_imp is simultaneously lowered from 10% to 5%. The final high state therefore differs from the low state in c_imp and T_N,min, not in the divertor neutral density that the text calls decisive. The stress-test note frames lower c_imp as an alternative driver, but from Sec 3.2 lowering c_imp from 10% to 2% makes the front drift below the separatrix and vanish; here lowering to 5% appears to pin the front at 12 cm rather than lift it. Either way the comparison is confounded, and \"neutral gas modulates XPR height\" is not cleanly established by these runs. The bolometry agreement of the tuned High-XPR case is good, but it is partly circular since c_imp, N_div, and T_N,min were chosen to make the front sit where the experiment sees it.\n\nOther limitations are proportionate: fixed c_imp in coronal equilibrium is a strong simplification (no impurity transport, no non-coronal states); the temporary neutral seed needed to initiate detachment means the simulation does not predict onset; the 0.1 ms averaging window is short; no code or data are posted. None of this breaks the central demonstration that a global turbulence code can hold a detached XPR-like state and produce plausible profiles. The physics claims about intermittency and enhanced radial transport are supported by the simulations as simulations; the link to ELM suppression is explicitly speculative.\n\nWho is this for: edge/SOL turbulence and power-exhaust people. It deserves a serious referee. The referee should ask for a cleaner height-modulation experiment (vary N_div alone, or vary c_imp alone) and for more transparency on the scan outcomes, not for a rewrite of the whole paper.","headline":"First turbulence-resolving XPR simulation, worth a serious referee, but the neutral-gas height-modulation claim is confounded by simultaneous changes in c_imp and T_N,min.","tokens_in":17328,"tokens_out":2440,"would_cite":true,"duration_ms":23412,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Global turbulence simulations reproduce a detached X-point radiator in a tokamak and show that the radiating front is strongly turbulent, a result the paper links to ELM suppression.","keywords":["X-point radiator","divertor detachment","edge turbulence simulation","GRILLIX","impurity radiation","neutral gas physics","ASDEX Upgrade","ELM suppression"],"falsifier":"A turbulence simulation that evolves the nitrogen charge-state distribution self-consistently, or uses a kinetic treatment of neutrals, should be run for the same discharge; if the front height no longer sits near 12 cm or the radiated fraction moves away from about 80%, the fixed-fraction and fixed-neutral-density assumptions are doing the work attributed to the physics.","tokens_in":16165,"feed_emoji":"⚛️","tokens_out":10620,"duration_ms":90425,"temperature":0.7,"pith_summary":"Predicting fusion power exhaust requires simulations that resolve both the turbulence in the tokamak edge and the atomic physics of detachment. This paper claims that the GRILLIX edge turbulence code, extended with impurity-radiation cooling, can sustain a detached X-point radiator (XPR) and reproduce measurements from an ASDEX Upgrade discharge. Two simulated states place a dense nitrogen radiation front 5 and 12 cm above the X-point, radiating about 80% of the injected heating power, with the taller front matching the measured bolometry best. The paper's new content is what turbulence does to the front: density and temperature fluctuations exceed 400% of the background, the front becomes a set of intermittent cold recombining spots wrapped in fluctuating ionization and radiation mantles, and a stationary E×B vortex appears with radial flows an order of magnitude above normal outboard transport. Because tall X-point radiators are observed to suppress edge-localized modes, the paper connects these turbulence-driven flows to ELM suppression.","feed_headline":"First turbulence simulations reproduce a detached X-point radiator","feed_subtitle":"Two simulated fronts at 5 and 12 cm above the X-point radiate 80% of input power and match measurements.","key_machinery":"The load-bearing element is the impurity-radiation power density $p_{\\rm rad}^{\\rm imp} = -L_{\\rm imp}(T_e)\\,n^2\\,c_{\\rm imp}$, added to the electron temperature equation with a coronal-equilibrium nitrogen rate coefficient $L_{\\rm imp}(T_e)$ that peaks near 7–15 eV. Around it sit two numerical and physics additions: a fluid neutral model whose perpendicular diffusivity has a temperature floor $T_{N,\\min}$ meant to mimic kinetic non-local neutral spreading, and an implicit treatment of parallel resistivity that removes the $T_e^{-3/2}$ stiffness of Ohm's law at low temperature. The fixed impurity fraction determines whether a radiation front can persist; the divertor neutral density and the neutral diffusivity floor control its height. On top of this, the turbulence is generated by GRILLIX's full-$f$ drift-fluid model, which evolves the whole profile rather than only fluctuations, letting the extreme fluctuation levels develop self-consistently.","core_discovery":"The paper's central claim is that a global, full-f turbulence code can reproduce a detached X-point radiator state in a tokamak. With nitrogen modeled as a fixed fraction of electron density, $c_{\\rm imp}=10\\%$ (later $5\\%$), radiating in coronal equilibrium, GRILLIX forms dense radiation fronts at 5 cm and 12 cm above the X-point that radiate about 80% of the input heating power. The taller front matches the experimental divertor temperature and bolometry profiles of ASDEX Upgrade discharge #40333 especially well. The discovery is that the XPR is not a smooth radiating shell but a strongly turbulent structure: intermittent cold recombining cores sit inside fluctuating ionization and radiation mantles, density fluctuations reach more than 400% of the background (40% in the attached reference), and the averaged electrostatic potential becomes poloidally asymmetric, producing a stationary E×B vortex whose radial particle fluxes exceed outboard turbulent transport by up to an order of magnitude. The paper argues these effects may explain the ELM suppression seen when the XPR is tall.","pith_inferences":["If the fixed-fraction and coronal assumptions were replaced by self-consistent nitrogen transport, the front height and radiation fraction might shift; that calculation would separate physics from prescription and is a direct next step.","The negative and space-dependent effective diffusivities found near the XPR imply that mean-field transport codes used for reactor extrapolation may need qualitatively new closure models in detached regimes, not just higher diffusivity values.","The E×B vortex intensity is a testable predictor for ELM suppression: scanning simulated XPR height and correlating vortex-driven radial transport with the experimentally observed height threshold around 7 cm would check the causal link directly.","Since averaged atomic rates differ from rates evaluated on averaged profiles, synthetic diagnostics applied to the instantaneous turbulent fields could quantify how much the intermittency changes predicted ionization, recombination, and radiation balances."],"forward_implications":["Detached X-point radiators can now be studied in global turbulence simulations, so the radiating front, neutral gas, and turbulent transport no longer have to be packaged into heuristic transport coefficients.","Neutral gas acts as the control knob for XPR height: raising the divertor neutral density and the minimum neutral diffusivity moves the front from 5 cm to 12 cm above the X-point.","The High-XPR simulation reproduces the measured divertor temperature and peaked bolometry profile, supporting the experimental picture of a tall radiating front in discharge #40333.","The time-averaged structure of the turbulent XPR matches earlier transport simulations, but only the turbulence simulation shows the intermittency: cold recombining cores and fluctuating ionization and radiation mantles that a mean-field picture misses.","The XPR shifts the radial electric field well inward, breaks poloidal potential symmetry, and drives an E×B vortex with radial flows an order of magnitude above outboard transport, which the paper links to the observed ELM suppression."],"supporting_citations":[{"why":"Supplies the fluid neutral model with charge-exchange diffusion and the implicit pressure-diffusion scheme used to reach detachment.","marker":"[19]"},{"why":"Defines the global electromagnetic drift-fluid plasma model that GRILLIX evolves, including the free-streaming heat-flux closure.","marker":"[51]"},{"why":"Provides the coronal-equilibrium nitrogen radiation rate coefficient $L_{\\rm imp}(T_e)$ used in the impurity cooling term.","marker":"[61]"},{"why":"The SOLPS-ITER transport simulation of the ASDEX Upgrade XPR whose mantle-core structure and height-control picture the turbulence results are compared with.","marker":"[36]"},{"why":"The analytical power-balance model for XPR stability and height that motivates treating the front as a radiation-conduction balance.","marker":"[59]"},{"why":"Documents the experimental X-point radiator, its feedback control, and the ELM-suppressed regime that the paper's turbulence results are linked to.","marker":"[28]"},{"why":"Describes the fully detached L-H transition discharge #40333 and the experimental XPR height used as the simulation target.","marker":"[30]"},{"why":"An earlier detached-plasma turbulence simulation in TCV whose inward $E_r$ shift the paper confirms and extends.","marker":"[20]"},{"why":"The original GRILLIX paper providing the flux-coordinate-independent grid and numerical framework the simulations run on.","marker":"[45]"},{"why":"Established the role of neutral gas in validated GRILLIX edge turbulence simulations, supporting the conclusion that neutrals control detachment.","marker":"[52]"}],"fun_headline_variants":["Turbulent X-point radiator simulated in detached divertor","First full-f sim of XPR: turbulence and ELM link","Simulated XPR: 80% power radiated, front turbulent","X-point radiator turbulence: cold spots, strong flows","Detached XPR sim matches experiment, reveals turbulence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole result leans on treating nitrogen as a fixed slice of the electron density whose cooling is precomputed from temperature, and on setting the neutral gas density in the divertor and a minimum neutral diffusivity by hand; if real impurity transport, non-equilibrium charge states, or kinetic neutral spreading shift the front's height or brightness, the agreement with the experiment is not predictive.","fun_headline_variants_meta":{"raw":{"variants":["Turbulent X-point radiator simulated in detached divertor","First full-f sim of XPR: turbulence and ELM link","Simulated XPR: 80% power radiated, front turbulent","X-point radiator turbulence: cold spots, strong flows","Detached XPR sim matches experiment, reveals turbulence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00026,"raw_usage":{"total_tokens":1644,"prompt_tokens":1055,"completion_tokens":589,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":671,"completion_tokens_details":{"reasoning_tokens":505}},"tokens_in":671,"tokens_out":589,"duration_ms":6183,"temperature":1.0,"reasoning_tokens":505,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:10:27.005289+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A turbulence simulation that evolves the nitrogen charge-state distribution self-consistently, or uses a kinetic treatment of neutrals, should be run for the same discharge; if the front height no longer sits near 12 cm or the radiated fraction moves away from about 80%, the fixed-fraction and fixed-neutral-density assumptions are doing the work attributed to the physics.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the fluid neutral model with charge-exchange diffusion and the implicit pressure-diffusion scheme used to reach detachment."},{"cited_title":"Zholobenko, K","cited_arxiv_id":null,"evidence_quote":"Defines the global electromagnetic drift-fluid plasma model that GRILLIX evolves, including the free-streaming heat-flux closure."},{"cited_title":"https://open.adas.ac.uk/","cited_arxiv_id":null,"evidence_quote":"Provides the coronal-equilibrium nitrogen radiation rate coefficient $L_{\\rm imp}(T_e)$ used in the impurity cooling term."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The SOLPS-ITER transport simulation of the ASDEX Upgrade XPR whose mantle-core structure and height-control picture the turbulence results are compared with."},{"cited_title":"Stroth, M","cited_arxiv_id":null,"evidence_quote":"The analytical power-balance model for XPR stability and height that motivates treating the front as a radiation-conduction balance."},{"cited_title":"Bernert, F","cited_arxiv_id":null,"evidence_quote":"Documents the experimental X-point radiator, its feedback control, and the ELM-suppressed regime that the paper's turbulence results are linked to."},{"cited_title":"Bernert, T","cited_arxiv_id":null,"evidence_quote":"Describes the fully detached L-H transition discharge #40333 and the experimental XPR height used as the simulation target."},{"cited_title":"Mancini, P","cited_arxiv_id":null,"evidence_quote":"An earlier detached-plasma turbulence simulation in TCV whose inward $E_r$ shift the paper confirms and extends."},{"cited_title":"Stegmeir, A","cited_arxiv_id":null,"evidence_quote":"The original GRILLIX paper providing the flux-coordinate-independent grid and numerical framework the simulations run on."},{"cited_title":"Zholobenko, A","cited_arxiv_id":null,"evidence_quote":"Established the role of neutral gas in validated GRILLIX edge turbulence simulations, supporting the conclusion that neutrals control detachment."}],"review_version":1}