arXiv:1901.04041v3 [physics.plasm-ph] 1 Feb 2019 utain sn h atL-E ytmon system Li-BES COMPASS fast the layer using scrape-off fluctuations and edge of Characterization .Bencze A. .Szuty´anyiM. 3 2 1 Abstract. 2019 January E-mail: o caeo ae lmnsadacmayn oe (densi holes accompanying and filaments the layer of scrape-off validity the for demonstr i verify automat different also we for -beam We diagnostic, at test non-perturbative ratio of measurement profiles. probability density physics the elect sequential electron generalized atomic of plasma different interpretation for the characterizing and the how in affect show system can We Li-BES adva COMPASS edge fluctuations. the investigate the to the is in study profiles of this of density purpose fast The u ultra normally plasmas. is r of it of and reconstruction terms scenarios in plasma for of power kind diagnostic any full in operation its reached has COMPASS on nttt fPam hsc fteCS rge zc Republ Czech Prague, CAS, the of Physics Plasma of Institute Hungary Budapest, Physics, Szech´enyi for Gy˝or, Istv´an Hungary University, Center Research Wigner [email protected] eetyteLtimBa msinSetocp (Li-BES) Spectroscopy Emission Lithium-Beam the Recently 1 .Berta M. , 2 n h OPS Team COMPASS the and 2 .Buz´as A. , 1 .Hacek P. , ceetdtcin Using detection. event ic tcatcGarcia-model stochastic tgsadlimitations and ntages e sasadr tool standard a as sed 3 eino COMPASS of region ydeficits). ty .Krbec J. , t h sblt of usability the ate uieautomatic outine ic ailpositions radial o density ron nteractions system 3 , Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 2

1. Introduction 2. Li-Beam Emission Spectroscopy

Tokamak plasma region defined as a narrow, few cm Neutral alkali beams are routinely used in a number of band around the last-closed-flux-surface (LCFS) is con- fusion-related devices as non perturbative diagnostics sidered to be of crucial importance for understanding for time evolving density profile measurements [10, 11, different transport phenomena determining both par- 12, 13]. ticle and energy confinement times and the strength Experimental data analyzed in the present article of the plasma-wall interaction [1, 2, 3, 4, 5]. This is have been acquired by the Li-BES diagnostic system the region where the important and interesting physics (see Figure 1) which consists of two main parts: the of the L-H transition, the generation of turbulence in- high energy (60 keV) diagnostic Li-beam [14] and the duced mesoscopic shear flows takes place and where Li-light detection system [15]. The Li-gun ion source the scrape-off layer (SOL) filaments, known as ’blobs’ emits Li+ ions which are further accelerated in two are born. This region is well known to be rich in stages (extraction and acceleration), in our case up to different stable and unstable modes like drift waves 60 keV. The ion beam is focused in the ion optics (DW), neoclassical tearing modes (NTMs), edge local- before entering the neutralization chamber where it ized modes (ELMs) etc. The unstable modes through gets neutralized by charge exchange processes induced nonlinear mode couplings can and will lead to coher- by collisions with sodium atoms. The neutral beam ent structure formation. The kinematics and dynam- reaches the plasma without any significant energy ics of these structures (which usually take the form of loss. During the plasma-beam interaction, through 3D filaments) is determined by the actual plasma pa- collisions with the plasma particles, Li-atoms are rameters both global and local, the underlying shear excited mainly to 2p state by the plasma particles. flows and the interaction of filaments with the mate- The excited 2p state decays with a τlife ≈ 27 ns rial components such as and limiter structures. lifetime with the emission of a photon of characteristic Theoretical descriptions of the coherent structure gen- wavelength (λ = 670.8 nm). These photons can be eration, structure propagation and detection statistics observed using an appropriate 2 nm wide interference have reached different levels of sophistication and de- filter and various detector systems (CCD camera, scriptive power. Some of the theories are competitive photomultiplier, photodiode or Avalanche photodiode [6, 7, 8, 9], therefore any experimental insight is of sig- (APD)). nificant value. For the present analysis we use the data acquired The present work is intended to be a comprehen- by a 18 channel APD-detector array [16] with sive contribution to the growing body of experimental observation volumes aligned along the major radius of data concerning the characterization of coherent struc- the at the outer midplane. The distance of tures: filaments or blobs in the scrape-off layer, holes neighboring observation volumes is 10 mm. The actual or density deficits and other fluctuating phenomena in radial resolution of the measurement strongly depends the plasma edge region. At the same time our intention on the nonlocal nature of the light emission due to the is to place the recently built COMPASS Lithium Beam finite lifetime of the exited Li(2p) states. The upper Emission Spectroscopy (Li-BES) system, on the map limit for this non-locality can be easily calculated as: of plasma turbulence studies, demonstrating its capa- 2E τ 2 bilities and uncover some of the interpretation issues ∆r = beam life , (1) that may arise during fluctuation data analysis. s mLi The rest of this paper is organized as follows: in where Ebeam is the beam energy (most of our section 2, the experimental setup and the principles of measurements were done using 60 keV beam energy), Li-BES diagnostic are discussed. subsection 2.1, qual- τlife is the lifetime of the excited Li(2p) atomic state itatively describes the raw data identifying the main and mLi is the mass of individual Li atoms. According noise sources. The effect of the plasma-beam interac- to the Eq. (1), in our measurement ∆r ≈ 35 mm. tion on the measured signal is described in details in This value sets a theoretical maximum for the smearing subsection 2.2. Subsection 2.3 is devoted to the fluctu- effect of the beam. The effective value of this non- ation amplitude profiles and the variation of the auto- locality is lower due to the collision induced de- correlation functions with the radial position. Section excitation and ionization depending on the actual 3 summarizes the analytic description of the Garcia- plasma density profiles [17]. In [18] a basic sensitivity model (subsection 3.1) and presents the results of the study has been performed concerning the response of synthetic diagnostic (subsection 3.2). Before analyzing the light emission profile to the local perturbations of the actual experimental data, two event detection tech- the density profile, concluding that the sensitivity is niques are compared (subsection 3.3). After presenting strongly suppressed at the maximum of the emission the experimental results in subsection 3.4, a summary profile. Later in this paper we are going to analyze the is given. Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 3

0 200 8

150 6 100 4

0 [kA] [m p e I

50 2 n 0 0 0 0.2 0 0.1 Li-BES

Hα 00 [a.u.] Intensity 0.0 1000 1100 1200 1300 t [ms] 0 Figure 2: Time traces of the main parameters of the discharge # 16072. In addition the Li-BES time signal 0 is shown for the brightest APD channel.

0 discharges in order to provide good statistics for in- depth plasma fluctuation studies. Medium density, L- 0 mode plasmas were produced with a stable 150 ms flat- 0 0 0 top. The main plasma parameters of the discharges are shown in Figure 2. For our practical purpose we separate the outer region of the toroidal plasma in two parts: the edge Figure 1: Experimental setup showing the Li-beam plasma and the SOL. In this article we define the path and the observation directions (APD line of edge plasma as the region of closed magnetic surfaces sights) relative to the tokamak vessel and the last just few cm inside the last closed flux surface (LCFS), closed surface (LCFS). meaning two observation channels (R12 = 709.7 mm, R13 = 719.4 mm) in the Li-BES system of COMPASS, while the SOL plasma is the region of open field lines effect of atomic physics on the statistical properties of and in this article we consider three Li-BES channels fluctuating Li-BES signals. (R15 = 738.7 mm, R16 = 748.3 mm, R17 = 758.0 mm) The detected light intensity is - in the first for our SOL analysis. approximation - proportional to the local electron It is quite obvious that the quality of the Li- density, therefore observing the time average intensity BES data is highly correlated to the reliability of and its fluctuations along the beam, it is possible the outcome results, therefore the assessment of data to reconstruct the density profile and the quasi-2D quality is very important and often neglected issue in correlation functions of the electron density fluctuation such data analysis. [19, 20]. Using a pair of deflection plates the beam In experimental physics the analyzed time traces can be either chopped or poloidally deflected (swept). are usually thought to consist of a series of instances Chopping the beam makes possible to correct for the taken from some statistical ensemble. The assumption background, while the poloidal deflection allows quasi- is that the statistics is determined by the – usually 2D measurements [21]. The switching frequency of nonlinear – dynamics of the physical system. In real the sweeping/deflection system presently works up to world various noise sources can corrupt and hide the 250 kHz. The poloidal resolution is limited by the interesting statistics of the system under investigation. ≈ 2 cm beam diameter. In the experimental campaign This is very much the case with the fluctuation data of which data are analyzed in this article, the quasi-2D detected by the Li-BES system. In our case the measurement was not available. interesting phenomena is the fluctuating turbulent plasma flow driven by local micro instabilities, while 2.1. Data characterization the main sources of measurement uncertainty can be identified as: Data analyzed in this work have been obtained in the experimental campaign CC18.18 dedicated to (i) Photon statistical noise acquiring large sets of fluctuation data in identical (ii) Amplifier noise Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 4

10 19 12 8 tanh 3. 10 saturation ) slope angle at value: n -3 s 6 8 inflection:

6 4

4 perturbated light 2 2. light profile Light Intensity (a.u.)

2 Density profile ( m 1. perturbation 0 density profile 0 0 2 4 6 8 10 r (cm)

Figure 4: Artificial density profile and calculated light profile with added perturbation and response. See detailed explanation in the text.

density is indirect. A natural question arises: what is the effect of plasma-beam interactions on statistical Figure 3: Power spectra of the fluctuating light signal quantities such as fluctuation amplitude profiles, from a SOL APD channel. correlation functions, conditionally averaged wave forms etc. To study the question a new reconstruction code has been developed at COMPASS [22] based on (iii) Pick-up oscillatory noise Schweinzer’s density profile reconstruction algorithm (iv) Fluctuations in the plasma background [23]. The method gives an approximate solution of the Figure 3 shows spectra of a typical Li-BES signal. The collisional-radiative model including electron-impact beam has been switched on and off with 20 ms period excitation, ionization and charge exchange processes: in order to asses the effect of plasma background and 5 d of the electronic noise. After the plasma discharge N (z)= [n (z)a (T (z)) + b ] N (z), (2) dz i e ij e ij j ends (after shot period, ASP) the beam observation j=1 X continues. When the beam is off in the ASP, the pure th where Ni is the occupation density of the i excitation electrical noise is detected (see the cyan line in Figure level, z is the coordinate along the lithium beam, where 3, this is also the sum of items (i)+(ii)+(iii) in the list z = 0 defines the position of the first measurement above), some sporadic narrow peaks are seen in the channel (called channel # 18). The coefficients spectrum. The same peaks are seen when the beam is aij (i 6= j) describe the transition rate from level on in ASP (see the red line in Figure 3), showing that i −→ j due to collisions with plasma electrons the presence of the beam does not introduce additional and ions. Attenuation of the ith atomic state due noise to the spectra. We should mention the resonance to charge exchange and ionization is included in around 650 kHz which gets amplified and broadened coefficient aii. The coefficients bij are the Einstein when the high photon flux hits the detector during the coefficients of spontaneous emission. ne and Te are the discharge (see the green and blue lines in Figure 3). electron density and temperature, respectively. This These high frequency features are irrelevant from the linear differential equation system describes the light point of view of edge turbulence studies, since we are emission of a Li-beam passing through a given plasma interested in the 1 − 100 kHz spectral range. It can be density and temperature profile. also seen, that the background fluctuation levels are at Figure 4 shows an artificially constructed density least one order of magnitude lower then the beam light profile (black curve) which consists of a tanh function fluctuations in the relevant spectral range, assuring smoothly merged to a linear function, matching the that the detected fluctuations are well localized at the observed shapes of real density profiles [24]. The beam position. shape of the profile is controlled by two parameters: (i) the gradient at the inflection point (α is the angle 2.2. Effects of plasma-beam interaction associated with the tagent line), (ii) the density value As the BES system detects the light emitted by plasma- where the tanh function saturates (ns). The density beam interactions, the information about plasma profiles constructed in such way have been perturbed Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 5

0

0 (a) 00 (d) 00 00 Li-light response[a.u] Li-light

0 (b) (e) 0

Shift [cm] Shift 0

(c) (f)

Broadening [cm] Broadening 00 0 0 0 δ0

Figure 5: The amplitude response (a), the offset of the response (b) and the widening of the response (c) as a function of perturbation amplitude at different radial locations with respect to the boundary of simulation domain. The amplitude response (d), the offset of the response (e) and the widening of the response (f) as a function of profile gradient at different saturation densities (see Figure 4). with localized Gaussian density pulses. Using the against the perturbation amplitude. The subplot (a) collisional-radiative model (2) the light emission can reproduces the findings in [18]: the amplitude response be calculated. The calculations shows some distinct is linear for small perturbations and a small deviation features of the light emission response (see Figure 4): is observed for larger δn. It is also clear that for a • amplification/attenuation of the perturbation given perturbation strength the response decreases as amplitude we travel uphill the profile. According to Figure 5.b • radial shift of maximum perturbed emission the maximum of the perturbation is shifted towards the center by a considerable amount of ≈ 0.5 cm • radial broadening of the perturbation for low perturbation amplitudes. Less radial shift is • decrease in the emission behind the light profile observed for larger amplitudes and for larger r values. maximum Depending on the amplitude and the initial position Recently it has been demonstrated that the of the perturbation, a well localized perturbation gets relative emission response to density perturbation broadened (smeared) due to the finite lifetime of the decreases with the radial position of the perturbation exited Li(2p) atomic states. Figure 5.c reveals that [18]. The present work expands on this idea as larger r values favour localization. we examine the characteristics of the light response Figure 5.d-f presents the situation when a fixed to various density perturbations of different density amplitude density perturbation is applied to the profiles. density profile at the radial position where the Let us first examine Figure 5.a-c where the tanh profile reaches its inflection point (maximum properties of light emission response are plotted Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 6 gradient). The response has been calculated for two ns saturation densities and different profile gradients L-mode 0.3 α (see Figure 4). As the density profile becomes H-mode 15 steeper the sensitivity to a given perturbation increases 0.2 (a) (b) independently of the saturation density (see Figure 10 5.d), while the radial shift of the perturbation does 0.1 5 not depend neither on α nor on ns as shown in Figure 5.e. For low densities the perturbation broadening is Ampl. Fluct. Rel. 0.0 quite large and it is basically independent of the profile 10 15 Ampl. [mV] Abs. Fluct. 10 15 gradient, on the other hand for higher density and Channel Channel shallow profiles the smearing is considerably smaller due to the collisional deexcitation processes.

2.3. Fluctuation amplitudes and correlations Before we turn to the detailed analysis of scrape-off layer structures (density filaments/blobs, holes), we shall describe the global characteristics of measured Li- light fluctuations (as a proxy for the electron density fluctuations). The relative and absolute fluctuation levels are plotted in Figure 6.a-b respectively, and are calculated as:

δIabs = σa − σb, (3)

δIabs δIrel = , (4) µa − µb Figure 6: Relative (a) and absolute (b) fluctuation where σa and µa are respectively the standard amplitude profile of Li-light signal in L (red) and inter- deviation and the time average of the signal Sa(t) = ELM H-mode (blue). Radial variation of the auto- Sbeam(t) + Sbackground(t). This is the superposition correlation function shape (c). The dashed black line of plasma background emission and the Li-beam depicts the separatrix position. emission. σb and µb are the standard deviation and the time average of plasma background signal Sb(t) = Sbackground(t) respectively. Applying a 3. SOL plasma fluctuations properly designed FIR bandpass filter, Sa(t) and Sb(t) contain only the relevant frequency components in the In this section we give a detailed description of range of f ∈ [1, 100] kHz. the observed statistics of scrape-off layer fluctuations A clear difference between L-mode and H- measured by Li-BES. Prior to the discussion of the mode fluctuation strength is visible in both absolute blob detection and statistical analysis, we present and relative fluctuation amplitude. The absolute the stochastic modeling of Garcia and a computer fluctuation level in L-mode is as much as five times simulation based on this model, simulating both the higher than in H-mode, while the relative fluctuation density and Li-light emission fluctuations. amplitude shows a factor of two difference. This behaviour is well known and generally observed in 3.1. The Garcia model all toroidal plasma devices (see e.g [25] and the references therein). The lower plot in Figure 6 shows In [27] Garcia presented a stochastic model describing the radial variation of auto-correlation function shapes. the fluctuations of scrape-off layer plasma as perceived In the scrape-off layer fluctuations with 10% − 40% by single point measurements. The observed behavior relative level are seen with auto-correlation functions is dominated by large-amplitude bursts, commonly smoothly decaying with a decay length τdecay ∈ referred to as ”blobs”. The model represents the [20, 50µs] decreasing towards the separatrix. In the measurement as a random sequence of bursts: last closed flux surface region (plasma edge) the Φ(t)= Akψ(t − tk), (5) fluctuations are different as can be seen in Figure 6.c. k The observed wave-like fluctuations with ≈ 70 kHz X where A is the burst amplitude, t is it’s arrival time frequency are similar to the earlier Li-BES observations k k and ψ(t) is a fixed waveform. It also states that at Wendelstein 7-AS [26]. bursts arrive according to a Poisson-process, which Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 7 means that the waiting time between them follows an exponential distribution with rate 1 : τw r = 3 cm (b) τ r = 4.5 cm 1 − τ Pτ (τ)= e w . (6) τw Calculating the cumulants of the signals probability density function (or PDF: PΦ(Φ)) yields formulas for its skewness and flatness, that show a parabolic relation, which is characteristic of systems that are Light intensity (a.u.) Light intensity (a.u.) dominated by intermittent fluctuations: perturbed light unperturbed light 2 4 (a) I2I4 hA ihA i 2 r0 F =3+ 2 3 2 S , (7) I3 hA i th where In is the integral of the n power of the wave form: Figure 7: Time traces of simulated Li-emission signals ∞ n at two different radial positions (a). Effect of atomic In = dt[ψ(t)] (8) −∞ physics on a single blob around r0 (b). Z The possible origin of such quadratic relation has been discussed in [28]. For experimental consideration the on blob amplitude, scattered by a normal distribution burst waveform is taken as a ”step” rise followed by an − t to simulate the findings in [29]. exponential decay: ψ(t)=Θ(t)e τd and exponentially Blob wave fronts are propagated radially outward distributed burst amplitudes. The connection between while their amplitudes are exponentially decaying. For the skewness and flatness thus simplifies to: each measurement point the passing through time 3 (blob duration) is calculated for every blob, than an F =3+ S2. (9) 2 exponential rise and fall is added to the signal, creating τd double exponential waveforms. In case of short waiting times and slow decay (γ = τw is large) PΦ(Φ) approaches a normal distribution while in The poloidal shape of individual blobs is described the opposite limit it approaches Gamma distribution: by a Gaussian function with FWHM = 3 cm. Then the perturbed signal is added to a radial density profile − Φ 1 1 A lim PΦ(Φ) = lim e h i . (10) reconstructed from Li-BES measurements [22]. The γ→0 γ→0 Γ(γ) Φ amplitudes are scaled to match the relative fluctuation In [29] a slightly refined version of the model was amplitudes seen in the experiment. The average blob presented according to experimental results from detection frequency is set so that larger blobs arrive at Langmuir-probe measurements on the TCV tokamak. around 1 kHz. The shape parameters and the decay of The shape of the bursts that dominate the signal the blobs is ascertained using conditionally averaged have a double exponential waveform with amplitude waveforms (conditional averaging will be explained in and waiting time both fitting exponential distributions. the next subsection) from COMPASS experiments as The PDF of the signal fits a Gamma distribution which well. The radial velocity of the blobs is set to be around is in agreement with previous results [27]. m 300 s as in e.g. [30]. The average blob amplitude is arbitrary as it can be scaled while maintaining 3.2. Simulation and synthetic diagnostic statistics. After simulating a fluctuating density signal, the atomic physics of plasma-beam interaction is taken We created a MATLAB simulation of intermittent into account as it has been described in subsection 2.2. SOL fluctuations based on the Garcia model described This way we get a synthetic signal similar to what in subsection 3.1. This simulation connected to the we can obtain from real Li-BES measurements (see atomic physics model described previously serves as Figure 7). At first glance two important features can a synthetic diagnostic helping the interpretation of be noted: (i) negative spikes appear at r > r , where the Li-BES measurements. It basically models a 0 r is the radial coordinate, where all blobs are created single point measurement in a two-dimensional (radial- 0 (in this simulation r = 4 cm). The negative spikes poloidal) plane. The basic parameters of the emerging 0 or holes are artifacts caused by the beam loss due to blobs are randomly generated: the arrival times at the increased ionization in the higher density zone. (ii) first radial observation point and the blob amplitudes the positive spikes at r > r , but closer to r are the are drawn from an exponential probability distribution. 0 0 result of the smearing of the density perturbation by Individual blobs are initialized at fixed radial position the finite lifetime of the excited Li(2p) states. In Figure r (r = 0 is the outer boundary of the simulation 0 7, the subplot (a) depicts two time traces at different domain) and moved with velocities linearly dependent Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 8

(a) (b) PDF PDF PDF

Amplitude (a.u.) Amplitude (a.u.)

(c) Sim. density Sim. light Garcia model

Amplitude (a.u.) PDF

Figure 8: The PDF of a simulated signal, where both blobs and artifact holes are present (r > r ). The 0 Amplitude (a.u.) red and magenta lines are fitted Gamma distribution according to the predictions in [27]. Figure 9: (a) PDF of simulated density fluctuations (red line: fitted Gamma function), (b) PDF of radial positions: the blue curve displays only positive light fluctuations with Gaussian noise (red line: fitted normal distribution, green line: fitted Gamma fluctuations at r < r0, while the red curve shows both positive and negative fluctuations relative to the distribution), (c) comparison of the density fluctuation average light profile. It is very important to realize PDF with the Li-light fluctuation PDF (red line: fitted that the appearance of holes as artifacts can corrupt Garcia model prediction). our interpretation of the experimental results, therefore extra care is required during BES data analysis. Their amplitude distributions are showed for the simulated contribution to the PDF (Probability Density Function noise-free density and light fluctuations. The Gamma approximated by the amplitude histograms) of the distribution predicted by the Garcia model accurately signal also fits a Gamma distribution predicted by [27] fits the simulation density as well as Li-light data. (as seen on Figure 8). We applied a broadband Gaussian random noise to the simulated ”density” signal with standard deviation 3.3. Event detection: two approaches proportional to the signal RMS. This accounts for the The most frequently used method in the literature different noises captured by the observation system, for SOL blob detection is based on a predetermined, such as the electronic noise and the photon statistical otherwise arbitrary threshold. This threshold is noise in the limit of high photon fluxes. As a usually given in terms of the standard deviation σ reference, in Figure 9.a the PDF of the ”pure” density [29]. The detection algorithm scans the signal and signal is shown. A Gamma distribution fits perfectly finds the peaks above some prescribed σ value. The the simulation. The resulting PDF of the noisy detected peaks are considered to be part of individual synthetic measurement (Li-light fluctuations) is shown filamentary events, therefore we can calculate the in Figure 9.b, and exhibit a PDF with Gaussian center detection frequency, the amplitude distribution and the (red curve) and Gamma-like tail (green curve). We waiting time distribution. This approach is used in this can conclude that applying the effect of collisional- article for conditional averaging (see section 3.3.2). radiative processes to the density does not change Another recently introduced approach to detect the amplitude distribution significantly, the change in Edge Localized Mode (ELMs) events in different diag- PDF shape (observed in the experiments) is introduced nostic signals is the generalized Sequential Probability by the detection noise. See Figure 9.c, where scaled Ratio Test (gSPRT) [31]. In the present work we ap- Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 9 plied, for the first time, the gSPRT technique to blob detection. It is a unique characteristic of the method that the event duration can be determined accurately according to the intrinsic statistics of the experimental data.

3.3.1. The gSPRT method. The generalized method described in [31] is based on the classical Sequential Probability Ratio Test [32], but adapt the local physics conditions and get the required parameters from the measured diagnostic signal. The SPRT algorithm Figure 10: The calculated gSPRT signal (blue) is based on the assumption that the signal consists together with the original signal (green) showing the of two separate random processes and the algorithm starting and ending points as calculated by the gSPRT can decide for each sample what is the logarithmic algorithm. probability ratio of the two situations where the sample belongs to one or the other processes. This parameter, usually denoted by λn (see its definition in [31]), depends on the probability distributions (PDF) of the two processes. Two decision levels are defined for λn, discriminating between two processes: β A = ln , (11) 1 − α 1 − β B = ln , (12) α where we set α - the probability of ’false alarm’ to 0.2 and β the probability of ’missed alarm’ to 0.05. We have to note that the gSPRT method is not sensitive to small variations of these values. When λn is below level A the sample belongs to the first process (here inter-blob fluctuations), while it is above level B the Figure 11: The observed linear correlation between the sample belongs to the second process (here the blob blob amplitude and the blob duration using the gSPRT event). method. For defining λn directly from the experimental data we use the best fit two-parameter amplitude distribution (ADF) as a proxy for the PDF. We have in the case of the 2nd and 3rd blob in the Figure found that the ADF of the inter-blob periods can be 10. We label the detected blobs according to their well described by the Gaussian distribution G(µ0, σ0). intensity (size). E.g. a ’size 2’ blob has its maximal On the other hand, the best fit for the empirical ADF amplitude in the corresponding gSPRT signal in the for conditionally averaged blob event was the Inverse interval [B +1 · ∆l,B +2 · ∆l], where ∆l = (M − B)/n Gauss (also known as Wald) distribution W (µ, Λ). We and M is the overall maxima for the gSPRT signal and have found a strong correlation (> 95%) between the n is the number of blob size labels we use (in the present maxima of the calculated gSPRT signal (λn) for each analysis we use n = 5, therefor we have 5 different blob blob and the maxima of the original signal. As it is amplitudes in the analysis). There has been observed shown in Fig. 10 the gSPRT can determine the start a strong correlation between the duration time of blobs and the end of the ’blob event’. as measured by the gSPRT method and the blob When the gSPRT signal (λn) exceeds the decision intensity as shown in Figure 11. We show the outlier level B, the blob event starts and when the signal point at the ’size 1’ blobs as a warning sign which has falls below the level the event ends. Figure 10 also to be taken into account when low amplitude blobs are illustrates the problem of merged events. When two or considered. As the blobs (filaments) are considered more events overlap, the gSPRT signal could not fall to be ’coherent structures’, the blob amplitudes, below the decision level B and the algorithm could not with good approximation are constant as they travel distinguish between the events. This problem can be radially across the detection channel, therefore the handled by finding the local minima of the smoothed measured blob duration is determined by the blob gSPRT signal within a given merged blob, and put velocity. The analysis presented in Figure 11 shows the end and the start points there as it can be seen that with increasing blob amplitude the blob duration Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 10

1 5 x σ a. b. 0.8

2.5 x σ 0.6

1 x σ 0.4

[ [

Intensity (a.u.) 0.2 50 150 250 350 450 10 30 50 70 Sim time [microsecond] Sim + at. ph. 0 Figure 12: (a) Blob detection using the threshold 0 10 20 30 40 50 60 70 method and time window selection for conditional t ( s) average. (b) Conditional average blob, averaged over 3 Figure 13: Conditionally averaged waveforms calcu- 420 events (3 · 10 blob/sec) with amplitude > 2.5 × σ. lated from simulated density signal (blue line) and af- ter adding the atomic physics effects (red line). increases therefore the velocity decreases. Nevertheless we should note that the above statement does not necessarily imply a relation between the blob velocity BES signal and the spatial (B⊥) blob size. PDF fit Simulation 3.3.2. The threshold method and conditional averaging. 10 -2 As it has been mentioned above, the threshold-based

event detection followed by conditional averaging (CA) PDF is the conventional approach in the field of scrape- off layer fluctuation studies. In this section we 10 -4 briefly summarize our algorithm before we dive into data analysis. The algorithm starts with setting a reasonable threshold level, usually given in terms of standard deviation (e.g. 2.5 × σ, see Figure 12), then 0 5 10 the signal is being scanned starting from the beginning Light Intensity (a.u.) until a peak is found above the threshold. Around the Figure 14: Comparison of the measured PDFs in the peak the signal is cut out within a ∆t time window (in SOL with the simulated synthetic signal and the result the case seen on Figure 12, ∆t = 70 µs) and saved for of the curve fitting procedure. further analysis. This procedure is repeated for all peaks above the threshold. As a result a collection of ∆t long time series After filtering the raw Li-BES data to the is created. Conditionally averaged blob is obtained relevant 1 − 100 kHz frequency range, we calculated after averaging over this ensemble of time series. the amplitude histogram (PDF) of the Li-light In Figure 13 we show the result of the conditional fluctuations. It can be seen in Figure14 that averaging procedure applied to simulated density time the PDF is clearly positively skewed suggesting the traces. The comparison of the ”density blob” and presence of large amplitude events (filaments, blobs) the ”light emission blob” is plotted in terms of the in our experimental data. It is also clear that the conditional average waveforms. The smearing effect of experimental PDF does not fit the expected theoretical the plasma-beam interaction can be observed in the Gamma distribution described in [27]. Nevertheless widening of the waveform. The rising exponent of the it was possible to develop a fitting function as a fitted exponential is increased by 72% while the falling combination of two two-parameter distributions: a by 37%. Gamma function and a Gaussian function which gave excellent fit to the data. It was also possible to 3.4. Experimental observation of SOL fluctuations deduce, directly from the experimental data, the In this section we present the results of basic input parameters for synthetic signal simulations, statistical analysis of intermittent SOL fluctuations namely the waiting time distribution and the average (blobs and holes) as measured by the Li-BES system blob duration. Running the simulation with such at COMPASS. parameters, the resulting PDF is in very good Characterization of edge and scrape-off layer fluctuations using the fast Li-BES system on COMPASS 11

the density drops locally in the SOL, it can cause a simultaneous density increase behind the emission (a) (b) profile. As the analyzed discharges are made to be as much identical as possible, we are able to collect large amount of events (few thousand). In Figure 16 we summarize the essential results about the statistical Intensity [a.u] Intensity − properties of these fluctuation events. In Figure 16.a- Cross-correlation b the conditional average waveforms are plotted for − − − blobs in the SOL and holes around the separatrix. τ μ τ μ The average blob duration is ≤ 20 µs while the hole duration time is ≤ 10 µs. We have to note though Figure 15: (a) Conditional average waveforms that these values contain contributions from different measured at different radial positions, reference being size blobs, and as it was shown on Figure 11 the channel # 15. (b) Measured cross-correlation functions higher amplitude blobs have longer duration times, with reference channel # 15. which is the reason for the relatively large errorbars. Fitting double exponential to the conditional average agreement with the actual measurement (see the red waveforms, the rising and falling times can be obtained. points in Figure 14), suggesting that our simulation As a general statement we can say that the blob/hole approach is reliable. We have to note that the waveform rises approximately two times faster than experimentally observed intermittency is lower than it decays. The Figure 16.c-d shows the amplitude the Garcia-limit, yet the final conclusion of the model, distribution for large events. According to Garcia that the signal PDF fits a Gamma distribution, still model the amplitude distribution can be described holds. by an exponential function, which indeed is the case As we have already showed in subsection 3.2, in for our Li-BES measurements. The above mentioned the simulation hole artifacts can appear behind the stochastic modeling predicts exponential distribution maximum of the Li-light emission profile. This will for the waiting times as well. Both the blobs and the present a challenge to experimental data analysis due holes detected in our experimental data have waiting to the fact that the real holes (density deficit) are times following exponential distribution as seen in expected around the separatrix, where the blobs are Figure 16.e-f. From the waiting time distribution the blob/hole detection rate can be inferred, this gives created. 3 −1 for blobs, λblob ≈ 3 · 10 s and for holes, λhole ≈ Indeed, hole artifacts can be found in the 3 −1 experimental data as well, as it can be seen in 13 · 10 s . Figure 15, where we show the conditional average waveforms for different radial channels (Figure 15.a). 4. Summary Here the reference channel # 15 is located in the scrape-off layer while channels # 10 and # 11 are We have demonstrated that the COMPASS Li-BES behind the light profile maximum (channel #12). system can be used as a non-perturbative means for There is no time delay observed between conditionally plasma fluctuation study, although the interpretation averaged waveform maxima, showing that these events of the data is affected by the plasma-beam interaction. appear at the same time in the signal. The same Our sensitivity study shows that for given density conclusion can be drawn from the cross-correlation profile shape with density perturbation applied at functions shown in Figure 15.b. It is worth mentioning a given radial position, the amplitude response is that reference channel is also # 15 and time shifts proportional to the perturbation amplitude. On the relative to zero time lag can be observed for the SOL other hand if we keep the perturbation amplitude channels (# 14, #16) indicating radial propagation of constant, the light response is more sensitive to outer the blobs, but there is no shift for channel # 11, which perturbations and for density profiles with larger is consistent with the presence of hole artifacts. gradients. We have also found that the localization In order to meet the challenge presented in the (maximum shift and broadening) of the response previous paragraph, the statistical analysis of the blob is better for outer channel perturbations and it and hole events should start with a procedure sorting practically does not depend on saturation density. out the artifacts. This has been done based on We constructed a synthetic diagnostic based on the following process: if the algorithm finds a hole two-dimensional numerical simulation of moving ran- in a given channel and simultaneously finds a blob dom structures with prescribed spatial and temporal somewhere else it qualifies the hole as a fake. Fake shapes. From the simulated data, using a collisional- blob appearance is also possible, if for some reason radiative model we have generated synthetic Li-BES REFERENCES 12

10 0 5 ch16 (a)( ¡ (c) ch16 #blob: 6404 10 3 #blob: 6404 4 f: 3.2 kHz f: 3.2 kHz Exp. fit Exp. fit 3 2 10 -2 10 2 PDF PDF

1 Intensity (a.u.) 10 1

0 CA waveform Rising Exp. 6.80 10 -4 -1 Falling Exp. -8.67 20 40 60 4 6 8 10 0.5 1 1.5 2 t ( s) Relative Blob Amplitude Waiting time (s) 10 -3

4 1 0 ch14 10 (d)10 (e)¢£¤ ch14 #hole: 26733 #hole: 26733 f: 13.4 kHz f: 13.4 kHz Exp. fit 0 Exp. fit

-1 10 2 PDF PDF -2 10 -2 Intensity (a.u.)

-3 0 CA waveform 10 Rising Exp. 4.79 -4 Falling Exp. -2.31 10 20 30 40 -5 -4 -3 -2 2 4 6 t ( s) Relative Hole Amplitude Waiting time (s) 10 -4

Figure 16: Fundamental statistical properties of blobs and holes measured by COMPASS Li-BES. signals. The analysis revealed the appearance of ar- Acknowledgement tificial events (holes and blobs) due to the plasma- beam interactions. These ”fake” events should not This work was partly supported by Czech Science be included in the analysis. We have shown that the Foundation Project No. GA16-25074S. PDFs, therefore the statistics, are not affected by the This work was also partialy supported by project atomic physics effects, but the photon noise and elec- No. CZ.02.1.01/0.0/0.0/16 019/0000768 and MEYS trical noises change the statistics. project LM2015045. We have shown that beside the generally used threshold method, the sequential probability ratio test References method can be successfully applied for blob detection, having the advantage of more accurate determination [1] Benjamin A. Carreras, “Plasma edge cross-field of event duration directly from the signal statistics. transport: experiment and theory”, Journal of We have also demonstrated that the Li-BES data Nuclear Materials, vol. 337-339, pp. 315 – 321, of intermittent fluctuations in the scrape-off layer 2005, PSI-16. can be described by the stochastic Garcia model for [2] V. Naulin, “Turbulent transport and the plasma convoluted PDFs assuming finite intermittency. It is edge”, Journal of Nuclear Materials, vol. 363-365, also interesting to note that we have found a change pp. 24 – 31, 2007, Plasma-Surface Interactions-17. in the correlation-functions of Li-light fluctuations [3] S. I. Krasheninnikov, D. A. D’Ippolito, and as we approach the separatrix from the scrape-off J. R. Myra, “Recent theoretical progress in layer. The origin of the wave-like correlations at understanding coherent structures in edge and sol the last closed flux surface is not yet clear. Quasi- turbulence”, Journal of Plasma Physics, vol. 74, two dimensional Li-BES measurements are planned no. 5, pp. 679717, 2008. for further investigation of turbulence phenomena at COMPASS. [4] O. E. Garcia, “Blob transport in the plasma edge: a review”, Plasma and Fusion Research, vol. 4, pp. 019–019, 2009. REFERENCES 13

[5] D. A. D’Ippolito, J. R. Myra, and S. J. Zweben, experiments”, Review of Scientific Instruments, “Convective transport by intermittent blob- vol. 89, no. 1, pp. 013503, 2018. filaments: Comparison of theory and experiment”, [14] G. Anda and et al., “Lithium beam diagnostic Physics of Plasmas, vol. 18, no. 6, pp. 060501, system on the COMPASS tokamak”, Fusion 2011. Engineering and Design, vol. 108, pp. 1–6, 2016. [6] B. A. Carreras, V. E. Lynch, and G. M. Zaslavsky, [15] M. Berta, G. Anda, A. Bencze, D. Dunai, “Anomalous diffusion and exit time distribution P. Hacek, M. Hron, A. Kovacsik, J. Krbec, of particle tracers in plasma turbulence model”, R. Panek, D. Refy, G. Veres, V. Weinzettl, Physics of Plasmas, vol. 8, no. 12, pp. 5096–5103, and S. Zoletnik, “Li-bes detection system for 2001. plasma turbulence measurements on the compass [7] R. S´anchez, B. Ph. van Milligen, D. E. Newman, tokamak”, Fusion Engineering and Design, vol. and B. A. Carreras, “Quiet-time statistics of 96-97, pp. 795 – 798, 2015, Proceedings of the electrostatic turbulent fluxes from the jet tokamak 28th Symposium On Fusion Technology (SOFT- and the w7-as and tj-ii ”, Phys. Rev. 28). Lett., vol. 90, pp. 185005, May 2003. [16] D. Dunai and et al., “Avalanche photodiode based [8] D Carralero, I Calvo, M Shoji, B A Carreras, detector for beam emission spectroscopy”, Rev. K Ida, S Ohdachi, S Sakakibara, H Yamada, and Sci. Instrum., vol. 81, pp. 103503, 2010. C Hidalgo, “Influence of on the self-similarity [17] O. Asztalos and et al., “Fluctuation response anal- properties of lhd edge fluctuations”, Plasma ysis of plasma turbulence with BES diagnostics”, Physics and Controlled Fusion, vol. 53, no. 9, pp. EPS Conference Proceedings, 2017. 095010, 2011. [18] M. Willensdorfer and et al., “Characterization of [9] P. Devynck, P. Ghendrih, and Y. Sarazin, “The the Li-BES at ASDEX Upgrade”, Plasma Phys. origin of the long time correlations of the density Control. Fusion, vol. 56, pp. 025008, 2014. fluctuations in the scrape-off layer of the tore [19] R. J. Fonck and et al., “Plasma fluctuation supra tokamak”, Physics of Plasmas, vol. 12, no. measurements in using beam-plasma 5, pp. 050702, 2005. interactions”, Rev. Sci. Instrum., vol. 61(11), pp. [10] S. Fiedler, R. Brandenburg, J. Baldzuhn, K. Mc- 3487, 1990. Cormick, F. Aumayr, J. Schweinzer, and H.P. [20] S. Zoletnik and et al., “Determination of Winter, “Edge plasma diagnostics on w7-as and electron density fluctuation via beam emission asdex-upgrade using fast li beams”, Journal of spectroscopy”, Plasma Phys. Control. Fusion, vol. Nuclear Materials, vol. 266-269, pp. 1279 – 1284, 40, pp. 1399–1416, 1998. 1999. [21] S Zoletnik, L Bardoczi, A Krmer-Flecken, Y Xu, [11] S. Zoletnik, G. Anda, M. Aradi, O. Asztalos, I Shesterikov, S Soldatov, G Anda, D Dunai, and S. Bat´o, A. Bencze, M. Berta, G. Demeter, G Petravich and, “Methods for the detection D. Dunai, P. Hacek, S. Heged˝ us, G. H. of zonal flows using one-point and two-point Hu, T. Krizsanczi, M. Lampert, D. Nagy, turbulence measurements”, Plasma Physics and J. N´emeth, M. Otte, G. Petravich, G. I. Pokol, Controlled Fusion, vol. 54, no. 6, pp. 065007, may D. R´efy, B. T´al, and M. V´ecsei, “Advanced 2012. neutral alkali beam diagnostics for applications in fusion research (invited)”, Review of Scientific [22] J. Krbec, P. Hek, M. Berta, J. Seidl, M. Hron, Instruments, vol. 89, no. 10, pp. 10D107, 2018. and R. Pnek, “Fast density reconstruction of li- bes signal on the compass tokamak”, Review of [12] D. I. R´efy, M. Brix, R. Gomes, B. T´al, Scientific Instruments, vol. 89, no. 11, pp. 113504, S. Zoletnik, D. Dunai, G. Kocsis, S. K´alvin, and 2018. T. Szabolics, “Sub-millisecond electron density profile measurement at the jet tokamak with the [23] J. Schweinzer and et al., “Reconstruction of fast lithium beam emission spectroscopy system”, plasma edge density profiles from Li I (2s-2p) Review of Scientific Instruments, vol. 89, no. 4, emission profiles”, Plasma Phys. Control. Fusion, pp. 043509, 2018. vol. 34, pp. 1173, 1992. [13] G. Anda, D. Dunai, M. Lampert, T. Krizsanczi, [24] E. Stefanikova, M. Peterka, P. Bohm, P. Bilkova, J. N´emeth, S. Bat´o, Y. U. Nam, G. H. Hu, and M. Aftanas, M. Sos, J. Urban, M. Hron, and S. Zoletnik, “Development of a high current R. Panek, “Fitting of the 60 kev neutral lithium beam injector for beam density and temperature profiles on the compass emission spectroscopy measurements on fusion tokamak”, Review of Scientific Instruments, vol. 87, no. 11, pp. 11E536, 2016. REFERENCES 14

[25] S J Zweben, J A Boedo, O Grulke, C Hidalgo, B LaBombard, R J Maqueda, P Scarin, and J L Terry, “Edge turbulence measurements in toroidal fusion devices”, Plasma Physics and Controlled Fusion, vol. 49, no. 7, pp. S1, 2007. [26] S. Zoletnik, M. Anton, M. Endler, S. Fiedler, M. Hirsch, K. McCormick, and J. Schweinzer, “Density fluctuation phenomena in the scrape- off layer and edge plasma of the wendelstein 7-as stellarator”, Physics of Plasmas, vol. 6, no. 11, pp. 4239–4247, 1999. [27] O. E. Garcia, “Stochastic Modeling of Inter- mittent Scrape-Off Layer Plasma Fluctuations”, Physical Review Letters, vol. 108, no. 26, pp. 265001, June 2012. [28] D. Guszejnov and et al., “On the effect of in- termittency of turbulence on the parabolic rela- tion between skewness and kurtosis in magnetized plasmas”, Physics of Plasmas, vol. 20, pp. 112305, 2013. [29] A. Theodorsen, O. E. Garcia, J. Horacek, R. Kube, and R. A. Pitts, “Scrape-off layer turbulence in TCV: evidence in support of stochastic modelling”, Plasma Phys. Control. Fusion, vol. 58, pp. 044006, 2016. [30] M. Lampert, S. Zoletnik, J. G. Bak, and Y. U. Nam, “2d scrape-off layer turbulence measurement using deuterium beam emission spectroscopy on ”, Physics of Plasmas, vol. 25, no. 4, pp. 042507, 2018. [31] M. Berta and et al., “Automatic ELM detection using gSPRT on the COMPASS tokamak”, Fusion Engineering and Design, vol. 123, pp. 950–954, 2017. [32] A. Wald, “Sequential Test of Statistical Hypotheses”, Ann. Math. Statist., vol. 16, pp. 117–186, 1945.