Distinguishing Cause from Correlation in Tokamak Experiments to Trigger Edge Localised Plasma Instabilities

The generic question is considered: How can we determine the probability of an otherwise quasirandom event, having been triggered by an external influence? A specific problem is the quantification of the success of techniques to trigger, and hence control, edge-localised plasma instabilities (ELMs) in magnetically confined fusion (MCF) experiments. The development of such techniques is essential to ensure tolerable heat loads on components in large MCF fusion devices, and is necessary for their development into economically successful power plants. Bayesian probability theory is used to rigorously formulate the problem and to provide a formal solution. Accurate but pragmatic methods are developed to estimate triggering probabilities, and are illustrated with experimental data. These allow results from experiments to be quantitatively assessed, and rigorously quantified conclusions to be formed. Example applications include assessing whether triggering of ELMs is a statistical or deterministic process, and the establishment of thresholds to ensure that ELMs are reliably triggered.

💡 Research Summary

The paper addresses a fundamental problem in tokamak plasma physics: how to rigorously determine whether an observed edge‑localized mode (ELM) has been caused by an external “kick” (a rapid vertical magnetic perturbation) or has occurred naturally. This question is crucial because controlled triggering of ELMs is a leading candidate for reducing the heat load on plasma‑facing components in future fusion power plants such as ITER and DEMO.

The authors adopt a Bayesian probability framework. They define three events: A = “an ELM occurs”, B = “the kick successfully triggers an ELM”, and C = “the kick starts to influence the plasma at time τₘ”. Using the law of total probability they write

 P(A|C) = P(A|B,C) P(B|C) + P(A|¬B,C) P(¬B|C).

Here P(B|C) ≡ P(K|τₘ) is the unknown trigger probability, while P(¬B|C)=1−P(K|τₘ). The key insight is that the kick can only affect the plasma for a short interval Δτ (the duration of the magnetic perturbation). By integrating the probability density over this interval they obtain

 P_Δτ = P(K) + P(¬K)