Pith. sign in

REVIEW 4 major objections 4 minor 21 references

Langmuir probe and infrared thermography measurements of wide and narrow heat flux profiles in the ST40 tokamak

T0 review · 4 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read In ST40 H-mode plasmas, the sub-millimetre near-SOL heat flux is carried by an enhanced electron current to the grounded divertor target, not by steeper density or temperature gradients.

desk verdict Solid single-pulse evidence that ST40's narrow near-SOL heat flux tracks the grounded divertor current, but the mechanism needs more data before it's secure. read the letter →

arxiv 2607.21362 v1 pith:SH65TSPL submitted 2026-07-23 physics.plasm-ph

T. Brewer , X. Zhang , C. Marsden , A. Scarabosio , J. Bryant , P. F. Buxton , M. Fontana , T. Gray
show 6 more authors
N. Lopez L. Martinelli A. Rengle E. J. C. Tinacba E. Vekshina H. V. Willett
This is my paper · ORCID
classification physics.plasm-ph PACS 52.55.Fa52.70.-m52.40.Hf
keywords Langmuirprobesheatfluxdecaylengthscrape-offlayersphericaltokamakdivertorelectroncurrentsheathtransmissioncoefficientinfraredthermography
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper compares infrared thermography and flush-mounted Langmuir probe measurements of parallel heat flux in the upper outer divertor of the ST40 tokamak, for similar L-mode and ELM-free H-mode pulses. The IR data show a very narrow (<1 mm) heat flux decay length near the strike point in H-mode, while the density and temperature profiles remain comparatively broad. The probes reveal a ground current density Jgnd near the strike point that is carried mostly by electrons and whose exponential decay length matches the narrow heat flux decay length. The paper concludes that the narrow near-SOL heat flux feature is carried by enhanced electron current to the target, and that the sheath heat transmission coefficient must include the Jgnd term to reproduce the measured narrow profile.

What carries the argument

The flush-mounted Langmuir probes are swept to obtain I–V traces fitted to Eq. (3), yielding Isat, Te, and Vf; setting Vb = 0 gives the grounded-target current density Jgnd. The parallel heat flux is then computed with a sheath transmission coefficient γ_sheath (Eq. 6) that contains Jgnd/Jsat terms rather than the usual constant. The matching exponential decay length between Jgnd and the IR-derived near-SOL heat flux is the argument that electrons, not profile gradients, carry the narrow channel.

What would settle it

Instrument the divertor tiles with current shunts to directly measure the current to the grounded target and simultaneously record the IR heat flux profile; if the directly measured current decay length does not track λq,near, the electron-current carrier claim is wrong. Alternatively, electrically isolate or bias the target tiles to remove the ground-current path and check whether the narrow heat flux feature disappears.

Watch

Extended reading notes

Core claim

Using six flush-mounted Langmuir probes and an IR thermography system on the ST40 upper outer divertor, the authors find that in ELM-free H-mode the profile of electron current to the grounded target, Jgnd, decays with a length scale close to λq,near, the sub-millimetre decay length of the narrow near-SOL heat flux feature. Neither electron density nor temperature profiles show such a narrow scale. The inferred sheath transmission coefficient, which peaks near the strike point when the Jgnd-dependent term is retained, falls to the standard 7–8 when zero ground current is enforced. The authors therefore attribute the narrow heat flux feature to an enhanced electron current into the target, po

Load-bearing premise

The entire argument rests on the inferred grounded-current density Jgnd from just six small flush-mounted probes; if those probes do not represent the electron current to the whole grounded target, the matching decay length with the IR heat flux could be a coincidence.

Editorial extensions

If this is right

  • The standard assumption of a roughly constant sheath heat transmission coefficient (~7–8) can miss the narrow channel; γ_sheath should be treated as a profile-dependent quantity when target current is significant.
  • ST40 divertor tile heat loads in H-mode may be concentrated in a sub-millimetre channel even when density/temperature profiles are broad, with direct implications for power-handling design.
  • Predictive scalings based on multi-machine λq databases may underestimate power widths for high-field spherical tokamaks that exhibit this electron-current-driven narrow feature.
  • Measurements of Jgnd can serve as a practical proxy for locating the narrow heat flux channel on the target.
  • The source of the enhanced electron current (for example, a fish-scale tile-induced potential hill) should be identified because it determines whether the effect is intrinsic to the divertor geometry.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the mechanism is electron current driven by a positive potential hill created by ion gyro-orbit losses on fish-scale tiles, then changing tile shaping or toroidal field direction should alter or suppress the narrow feature; this is a testable design knob that goes beyond the paper's speculation.
  • The same Jgnd-to-γ_sheath coupling should apply to other divertor geometries and tokamaks, so re-analyzing past flush-probe data with Eq. (6) might reveal hidden narrow features in machines where they were not previously identified.
  • If a substantial fraction of SOL power is carried by electrons as current, the plasma-wall interaction and detachment physics near the strike point may differ from the usual ion-dominated picture; controlling target current could then become a detachment-control handle.
  • The three-parameter fit and the Ti=Te assumption may hide non-saturation effects; the planned four-parameter fit with a wider voltage range will test whether the inferred Jgnd is quantitative or an artifact.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. This paper presents Langmuir probe (LP) and infrared thermography (IR) measurements of divertor heat flux in ST40, comparing an L-mode pulse (#13548) and an ELM-free H-mode pulse (#13599). The H-mode near-SOL parallel heat flux has a narrow decay length (lambda_q,near < 1 mm). From six flush-mounted LPs, the grounded-target current density Jgnd is inferred and found to have an exponential decay length close to lambda_q,near. Including Jgnd in the sheath heat transmission coefficient, Eq. (6), is needed for LP-derived q_parallel to match the IR narrow profile. The authors conclude that the narrow near-SOL heat flux is carried by enhanced electron current to the divertor.

Significance. If correct, the result provides a rare experimental observation connecting the extremely narrow heat-flux feature in a diverted H-mode to an electron-current channel, with direct implications for reactor divertor heat loads and detachment. The paper is valuable for the ST40 diagnostic description and for the LP/IR comparison. Strengths include the use of two independent diagnostics, explicit statement of assumptions (Ti=Te, no secondary electron emission, three-parameter fit), open-source equilibrium code, and honest reporting of fit R^2 and data rejection criteria. However, because the central mechanistic claim rests on Jgnd derived from a simplified I-V model and on a single H-mode pulse, the conclusion is plausible but not yet established to the standard claimed.

major comments (4)
  1. [§2.2, Eqs. (3)–(6), Fig. 3] The grounded-current density Jgnd is not an independent measurement: it is obtained by evaluating the three-parameter fit of Eq. (3) at Vb=0, with Ti=Te and no secondary electron emission. The same Jgnd enters gamma_sheath in Eq. (6), which is then used to compute LP q_parallel compared with IR in Fig. 3. Therefore, the reported 'satisfactory agreement' cannot validate the Jgnd model. The manuscript itself notes (§2.2) that a four-parameter fit is needed but has not been applied. I request a sensitivity analysis (e.g., four-parameter fit, Ti/Te variation, SEE correction) and, at minimum, error bars that propagate fit-parameter uncertainties into Jgnd and lambda_Jgnd.
  2. [§3, Fig. 6, Fig. 4] The central comparison is based on one H-mode pulse (#13599) and only six LP points. The Jgnd profile in Fig. 6(b) has no error bars, the exponential fit range is not specified, and the profiles were shifted manually by 1.4 mm (#13599) and 2.3 mm (#13548) in psi_N. With six points and fit-range freedom, the statement that lambda_Jgnd is 'close to' lambda_q,near is not quantitatively supported. Please report fit uncertainties, confidence intervals on lambda_Jgnd and lambda_q,near, the chosen fit domain, and ideally additional H-mode pulses or multiple time windows.
  3. [§3, Fig. 6 and Summary] The mechanistic conclusion ('narrow near SOL heat flux is carried by an enhanced electron current') rests on a matching decay length. This is a correlation between a probe-inferred current and an IR-inferred heat flux. If probe sheath effects bias Jgnd near the strike point (e.g., non-saturation or overestimated Te), the matching lambda could be coincidental. A concrete test would be to compare Jgnd with direct measurement of target current (e.g., through tile grounding), or to vary the target bias/grounding, or to demonstrate that lambda_Jgnd and lambda_q,near scale together across several pulses and conditions.
  4. [§3, Fig. 5] The claim that the narrow heat flux is not dominated by density or temperature gradients relies on comparing lambda_q,near with lambda_ne and lambda_Te in Fig. 5. The Te profile is described as 'slightly broader in H-mode but due to the noisy temperature data'; with noisy data and no error bars, this comparison is weak. Please provide uncertainties or at least show the scatter in Fig. 5 so the reader can assess whether the Te gradient is indeed too broad to explain the narrow heat-flux feature.
minor comments (4)
  1. [Throughout] Language issues: 'Niether' should be 'Neither' in §3; 'dishcarges' in §1; 'he current density' in §2.2; 'possessed minor fluctuations' in §3 is awkward. Please copyedit.
  2. [Fig. 3 caption] 'a singular chord' is ambiguous. Clarify which chord was used and how the IR profile was constructed; also specify the spatial averaging in the LP comparison.
  3. [Fig. 6] The figure lacks error bars on all profiles, including Jsat and Vf. Since the paper explicitly discusses error propagation from the fit, the absence of error bars in the key figure is surprising and should be addressed at least in the caption or text.
  4. [§4, Ref. [23]] The suggested mechanism (positive potential hill in front of fish-scaled divertor plates) is speculative and relies on a reference that is 'submitted'. This is acceptable as a hypothesis, but it should be explicitly labeled as speculative rather than as a supporting result.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Jgnd and IR-derived λq,near are independent measurements, and the LP/IR agreement provides external grounding.

full rationale

The paper's central comparison is between λJgnd, obtained from Langmuir probe I-V fits by solving Eq. (3) at Vb=0, and λq,near, obtained from IR thermography via FAHF/IRRITANT double-exponential fits. These are independent measurement channels: Jgnd comes from the probe current balance with no IR input, and the IR λq,near comes from surface-temperature inversion with no probe input. Although the LP q∥ used for validation in Fig. 3 is computed with γsheath that includes Jgnd (Eqs. 5-6), so the LP and Jgnd channels are not fully independent of each other, the agreement of LP q∥ with the IR data provides external grounding for that calculation. No fitted parameter is defined in terms of the target λq,near: the exponential fit to Jgnd is made to probe data alone, and the closeness of λJgnd to λq,near is an empirical comparison, not a construction. Self-citations are not load-bearing in a circular sense: Refs. [1,2] document prior ST40 observations, and Ref. [23] is cited only as a speculative mechanism, not as proof of the electron-current claim. The paper explicitly discloses limitations—the three-parameter fit with a four-parameter fit not yet applied, unclear ion-saturation regions, only six probes, no error bars on the decay lengths, and manual profile shifts of 1.4–2.3 mm—but these are measurement and correctness risks, not definitional circularity. Therefore no specific circular step can be exhibited, and the derivation chain is self-contained in the relevant sense.

Assumptions & free parameters 5 free parameters · 7 assumptions · 1 invented entities

The central claim depends on several modeling assumptions: the Langmuir probe interpretation, the sheath transmission coefficient, exponential profile shapes, and the equilibrium mapping. The only invented entity is the potential-hill mechanism, which has no independent evidence in this paper.

free parameters (5)
  • Profile alignment shift = 2.3 mm (#13548), 1.4 mm (#13599) at OMP
    Profiles are manually shifted so that peak heat flux aligns with ψN=1 (Section 3); this can affect the inferred decay lengths.
  • Fit range for exponential profiles = Not stated exactly
    Exponential fits to ne, Te, Jgnd and q∥ over an unspecified radial range; with only 6 probes, the fit range determines λ values.
  • λJgnd = Close to λq,near (<1 mm in H-mode), exact value not given
    The central claim rests on this fitted decay length from Jgnd data.
  • λq,near = <1 mm in H-mode
    The IR-derived narrow decay length, fitted from the double exponential.
  • Ti/Te ratio = 1 (assumed)
    Used in Eq. (6) for γsheath; if not 1, the inferred heat flux changes.
assumptions (7)
  • domain assumption Langmuir probe current-voltage characteristic: I = Isat [1 - exp((Vb - Vf)/Te)]
    Used to fit Te and Isat from LP sweeps; standard model but assumes simple sheath behavior.
  • domain assumption n_i = n_e and T_i = T_e in the sheath
    Invoked in Section 2.2 to compute ne,se and γsheath; not measured.
  • domain assumption Plasma contains only deuterium ions and secondary electron emission is negligible
    Assumed in Eq. (6) for γsheath; may not hold at high heat flux.
  • domain assumption Profiles decay exponentially in the SOL (ne, Te, Jgnd, q∥)
    Equations (7)-(8) and fits in Figs. 5-7; the narrow feature is characterized by a single exponential decay length.
  • standard math Sheath heat transmission coefficient formula from Stangeby (Eq. 6)
    Adopted from Ref [19]; standard model for a grounded surface.
  • domain assumption Equilibrium reconstruction (GSFit) accurately maps target locations to ψN and gives field-line angles α
    Used to map probe data to OMP and compute A_eff; errors affect all profiles.
  • domain assumption FAHF inversion of IR data yields accurate surface heat flux
    Assumes 2-D heat conduction with known material properties; from Ref [13].
invented entities (1)
  • Positive potential hill in front of fish-scaled divertor plates
    purpose: To explain the enhanced electron current to the target and thus the narrow heat flux feature
    Proposed via Ref [23] (submitted), not measured in this paper; no direct evidence provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Langmuir probe and infrared thermography measurements of wide and narrow heat flux profiles in the ST40 tokamak." pith.science (2026). https://pith.science/paper/SH65TSPL

@misc{pith2026260721362,
  author       = {Pith},
  title        = {Pith review of: Langmuir probe and infrared thermography measurements of wide and narrow heat flux profiles in the ST40 tokamak},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SH65TSPL}},
  note         = {Machine review of arXiv:2607.21362}
}
abstract

Very narrow heat flux profiles have been observed on ST40 with decay lengths in the near scrape-off-layer (SOL) of $\lambda_{q,near}<$~1~mm in H-mode plasmas. Measurements from divertor Langmuir probes are compared with an infrared (IR) thermography system for upper-single null diverted plasmas in L-mode and ELM-free H-mode. SOL current to the grounded divertor targets is measured from the Langmuir probes with profiles exhibiting a similar exponential decay to $\lambda_{q,near}$ suggesting the near SOL heat flux is related to enhanced electron current to the divertor. Inclusion of electron current contributions to the heat flux transmission coefficient is crucial in capturing similar narrow profiles to those from the IR thermography system. Ongoing upgrades to ST40 will enable more investigations on SOL power decay lengths on inboard and outboard targets in double and single null configurations.

Figures

Figures reproduced from arXiv: 2607.21362 by the authors.

Figure 1
Figure 1. IR camera data (gradient coloring) from pulse #13599 at [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Example of fitting done to raw data for ST40 pulse [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Target parallel heat flux from Langmuir probe and IR [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: (a) Plasma current, (b) smoothed line-integrated density, [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]
Figure 7
Figure 7. Figure 7: Sheath transmission coefficient for H-mode pulse #13599 [PITH_FULL_IMAGE:figures/full_fig_p005_7.png]
Figure 6
Figure 6. Figure 6: Target profiles of Jsat, Jgnd and Vf measured by divertor Langmuir probes. A fit similar to Eqs. (7) and (8) is applied to Jgnd with the length scale shown corresponding to the profile at OMP. Profiles of Jsat, Jgnd and Vf are shown in [PITH_FULL_IMAGE:figures/full_fi…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

21 extracted references · 7 canonical work pages

  1. [2]

    C. Marsden, et al., Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry, Nuclear Materials and Energy 41 (2024) 101773.doi:https://doi.org/10.1016/j. nme.2024.101773

    arXiv 2024

  2. [3]

    Eich, et al., Scaling of the tokamak near the scrape-off layer h-mode power width and implica- tions for iter, Nuclear Fusion 53 (9) (2013) 093031

    T. Eich, et al., Scaling of the tokamak near the scrape-off layer h-mode power width and implica- tions for iter, Nuclear Fusion 53 (9) (2013) 093031. doi:10.1088/0029-5515/53/9/093031

    doi:10.1088/0029-5515/53/9/093031 2013

  3. [5]

    P. Stangeby, et al., Power deposition on the diii- d inner wall limiter, Journal of Nuclear Materi- als 463 (2015) 389–392, plasma-Surface Interactions 21.doi:https://doi.org/10.1016/j.jnucmat. 2014.09.051

    doi:10.1016/j.jnucmat 2015

  4. [6]

    Arnoux, et al., Scrape-off layer properties of iter- like limiter start-up plasmas in jet, Nuclear Fusion 53 (7) (2013) 073016.doi:10.1088/0029-5515/53/ 7/073016

    G. Arnoux, et al., Scrape-off layer properties of iter- like limiter start-up plasmas in jet, Nuclear Fusion 53 (7) (2013) 073016.doi:10.1088/0029-5515/53/ 7/073016

    work page doi:10.1088/0029-5515/53/ 2013

  5. [7]

    F. Nespoli, et al., Understanding and suppressing the near scrape-off layer heat flux feature in inboard- limited plasmas in tcv, Nuclear Fusion 57 (12) (2017) 126029.doi:10.1088/1741-4326/aa84e0

    work page doi:10.1088/1741-4326/aa84e0 2017

  6. [8]

    D. Carralero, et al., An experimental investiga- tion of the high density transition of the scrape- off layer transport in asdex upgrade, Nuclear Fusion 54 (12) (2014) 123005.doi:10.1088/0029-5515/54/ 12/123005

    work page doi:10.1088/0029-5515/54/ 2014

  7. [9]

    B. LaBombard, et al., Particle transport in the scrape-off layer and its relationship to discharge den- sity limit in alcator c-mod, Physics of Plasmas 8 (5) (2001) 2107–2117.doi:10.1063/1.1352596

    work page doi:10.1063/1.1352596 2001

  8. [10]

    J. Hecko, et al., Experimental evidence of very short power decay lengths in h-mode discharges in the com- pass tokamak, Plasma Physics and Controlled Fu- sion 66 (1) (2023) 015013.doi:10.1088/1361-6587/ ad08f0

    doi:10.1088/1361-6587/ 2023

Show all 21 references
  1. [11]

    Asunta, et al., Overview of st40 results and fu- ture: expanding the physics basis of high-field spher- ical tokamaks, Nuclear Fusion 66 (11) (2026) 116004

    O. Asunta, et al., Overview of st40 results and fu- ture: expanding the physics basis of high-field spher- ical tokamaks, Nuclear Fusion 66 (11) (2026) 116004. doi:10.1088/1741-4326/ae4a3c

    2026 doi

  2. [12]

    McNamara, et al., Overview of recent results from the st40 compact high-field spherical tokamak, Nu- clear Fusion 64 (11) (2024) 112020.doi:10.1088/ 1741-4326/ad6ba7

    S. McNamara, et al., Overview of recent results from the st40 compact high-field spherical tokamak, Nu- clear Fusion 64 (11) (2024) 112020.doi:10.1088/ 1741-4326/ad6ba7

    2024

  3. [13]

    Moscheni, et al., Measurement of divertor surface heat flux by infra-red thermographic inversion in st40, Review of Scientific Instruments 96 (1) (2025) 013503

    M. Moscheni, et al., Measurement of divertor surface heat flux by infra-red thermographic inversion in st40, Review of Scientific Instruments 96 (1) (2025) 013503. doi:10.1063/5.0244494

    2025 doi

  4. [14]

    A. Podolník, et al., Interpretation of flush-mounted probe current-voltage characteristics using four- parametric fits, Plasma Physics and Controlled Fusion 61 (10) (2019) 105011.doi:10.1088/ 1361-6587/ab3de8

    2019

  5. [15]

    A. Kuang, et al., Design and operation of a high- heat flux, flush-mounted ‘rail’ langmuir probe ar- ray on alcator c-mod, Nuclear Materials and Energy 12 (2017) 1231–1235, proceedings of the 22nd Inter- national Conference on Plasma Surface Interactions 2016, 22nd PSI.doi:htt...

    2017 doi

  6. [16]

    G. F. Matthews, et al., Investigation of the fluxes to a surface at grazing angles of incidence in the toka- mak boundary, Plasma Physics and Controlled Fu- sion 32 (14) (1990) 1301.doi:10.1088/0741-3335/ 32/14/004

    1990 doi

  7. [17]

    J. P. Gunn, et al., Flush-mounted probes in the divertor plates of tokamak de varennes, Review of Scientific Instruments 66 (1) (1995) 154–159.doi: 10.1063/1.1145249

    1995 doi

  8. [18]

    P. F. Buxton, Gsfit (2025). URLhttps://github.com/tokamak-energy/gsfit 6

    2025

  9. [19]

    P. C. Stangeby, The Plasma Boundary of Magnetic Fusion Devices, 1st Edition, Institute of Physics Pub- lishing, 2000

    2000

  10. [20]

    Brunner, B

    D. Brunner, B. LaBombard, Surface thermocouples for measurement of pulsed heat flux in the divertor of the alcator c-mod tokamak, Review of Scientific Instruments 83 (3) (2012) 033501.doi:10.1063/1. 3689770

    2012 doi

  11. [21]

    R. P. Canale, S. C. Chapra, Numerical Methods for Engineers, 7th Edition, McGraw-Hill Education, 2014

    2014

  12. [22]

    De Oliveira, et al., Langmuir probe electronics up- grade on the tokamak à configuration variable, Re- view of Scientific Instruments 90 (8) (2019) 083502

    H. De Oliveira, et al., Langmuir probe electronics up- grade on the tokamak à configuration variable, Re- view of Scientific Instruments 90 (8) (2019) 083502. doi:10.1063/1.5108876

    2019 doi

  13. [23]

    C. S. Chang, et al., Positive sheath generation in front of fishscaled divertor plates, submitted to the Journal of Nuclear Fusion (2026). 7

    2026

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.