1. Introduction
Electron Cyclotron Resonance Ion Sources (ECRIS) are some of the most widespread devices used to generate highly charged ion beams of variable intensity to accelerators across the world. They are based on the dual concepts of ECR heating and magnetic confinement, whereby plasma electrons gyrating about a longitudinal magnetic field
B gain energy through resonance with circularly polarised electromagnetic (EM) radiation, and a min-B configuration traps them long enough to sequentially ionise atoms to high-charge states. The resultant plasma is composed of multi-charged ions immersed in a cloud of electrons of density
–
and temperature
–
, which makes it ideal for not just ion beam generation, but also for research in other applied disciplines. For all its advantages, however, the system is quite difficult to study. The complex energy transfer process on account of the peculiar magnetostatic field profile and multi-modal nature of EM wave launched into the chamber [
1], complicated transport phenomena [
2,
3], and the presence of instabilities [
4] render the plasma non-homogeneous and anisotropic. Additionally, ECR plasmas are also known to support multiple electron populations, of which the three most important classes are the cold (
–
), warm (
–
), and hot electrons (
-beyond).
Fundamental and application-oriented research require information on intermediate energy electrons spanning the boundary between warm and hot (
–
) because their properties govern the sequential ionisation process that forms the backbone of ECRIS operation. The PANDORA project is one such application which aims to measure
-decay rates of radioisotopes modified by the plasma environment, confining them in an ECR magnetic trap [
5]. The idea is to then use these experimentally measured rates to verify the theory put forward by Takahashi and Yokoi [
6] and in case of satisfactory match, extrapolate the same to stellar plasmas. The theory requires detailed inputs on the atomic level population and charge state distribution (CSD) of the ions which in case of ECR plasmas implies evaluating them as a function of position since the system is spatially non-homogeneous. The first step in achieving this is by characterising the spatial distribution of the aforementioned intermediate energy electrons.
To this effect, we present here a comprehensive analysis of the space-resolved properties of warm plasma electrons confined in compact magnetic traps, studied using a robust electron kinetics model and benchmarked with suitable experiments. In this article, we intend to elaborate on both aspects, but focus more on phenomenological plasma emission models connected with the latter, in an effort to demonstrate the utility of energy dispersive X-ray spectroscopy [
7,
8] when it comes to investigating intermediate energy electrons.
The contents of the article are divided into different sections as follows: in
Section 2, the self-consistent electron kinetics simulations are briefly summarised and a quick overview is provided on determination of space-resolved electron energy distribution functions (EEDF). This effectively resulted in a theoretical 3D map of warm electrons whose properties were to be corroborated by experiments.
Section 3 introduces the idea behind these validation experiments, and elaborates on the different kinds of measurements and their individual applicability. The focus is then shifted onto two particular experiments, namely 2D space-resolved X-ray imaging in photon counting mode and volumetric fluorescence/low energy bremsstrahlung emission spectroscopy, and corresponding theoretical models are presented, respectively, in
Section 4 and
Section 5. By comparing the model-predicted results with those obtained experimentally, the two methods, respectively, furnished information about the spatial distribution of warm electrons, and estimates on the absolute number density of confined electrons and characteristics of escaping electrons. We conclude with
Section 6, summarising the methodology employed and preliminary results on ECR plasma electron characterisation. The correctness of the latter is discussed with regards to the physics, and future improvements to the model are outlined in keeping with the requirements for the PANDORA project.
2. Space-Resolved Electron Kinetics: Theoretical Modelling
Numerical simulations are practical methods for probing microscopic properties of ECR plasmas because they offer reliable solutions to the several coupled differential equations describing such systems, which may otherwise remain unsolved. They are fast, accurate, and flexible with regards to model complexity. An iterative procedure to solve the collisional Vlasov–Boltzmann Equation [
9] was developed to obtain a self consistent picture of stationary ECR plasmas. The simulations aimed to solve the plasma particles’ equation of motion in the presence of an EM field self-consistently updated with the trajectories of the particles. Such an approach was essential because of the very nature of wave-particle interaction—the EM field dictates the motion of charged particles in the ECR plasma and the electrons’ energy through resonant interactions, but the field profile itself is affected by the 3D dielectric tensor [
10] calculated according to the spatial distribution of the particles, and hence to the electron density. The simulation model was generated using COMSOL Multiphysics
® as an FEM solver and MATLAB
® as a particle mover, based on a schematic outlined in
Figure 1.
Details about the simulation scheme and associated algorithm can be found in [
11], and experimental validation in [
12,
13,
14]. The routine was first applied to warm electrons alone; the phase space trajectories of
N = 40,000 macroparticles were followed for a fixed simulation time
with time step
and after each step, the code saved position and energy of particles in a 3D matrix corresponding to the simulation domain sliced into cells of
, producing occupation and energy maps, respectively, [
11]. The microwave frequency was taken as
while the power was
, matching the operating conditions of experiments as described in
Section 3. At the end of the iteration routine (after achieving steady-state) the occupation maps were scaled to density maps by assuming a total number of particles spread out according to a plasmoid/halo structure [
9], where plasmoid density ∼
and halo density one hundredth of it. The result was concise data on electron number and energy as a function of position in the plasma, expressed as a set of 7 number and 7 energy density matrices of dimensions
, corresponding to the energy intervals most occupied by warm electrons
,
,
,
,
,
, and
. Some of the XY-projection maps (along the Z axis) are shown in
Figure 2.
The discrete data thus obtained was subjected to post-processing for determining a phenomenological, space-resolved EEDF. First, the occupation matrices
and energy density matrices
were used to calculate average electron energy (AVE) in each cell according to the expression
and the plasma was divided into finer
regions of interest (ROIs) by grouping together cells with similar AVE. This facilitated analysis and helped resolving the inherent anisotropy in the plasma.
Figure 3 shows an isometric view of ROIs 1, 3 and 4, corresponding, respectively, to
= 0–0.1,
= 0.2–0.3, and
= 0.3–0.4
.
The second step involved an exhaustive comparison of the performances of different EEDFs in the various ROIs based on statistical metrics like MSE and
, after which a two-component distribution composed of Maxwell and Druyvesteyn distribution as shown in Equation (
2) was deduced as the phenomenologically and globally correct function.
Here,
–
refers to the temperature of the Maxwell distribution function which described cold electrons accumulated in the
interval, while
–
is the same for the Druyvesteyn function representing the warm electrons in the remaining energy intervals
.
are, respectively, the normalisation coefficients of the two components, while the energy intervals chosen could not sufficiently isolate the two populations, the phenomenological model adopted did manage to resolve them. It should be noted that though the form of the EEDF was seemingly uniform throughout the plasma, the defining parameters of the components, temperature and normalisation coefficients, varied spatially, reflecting the plasma non-homogeneity. More details on the choice of the EEDF, rationale behind multiple components, and statistical analysis can be found in [
15,
16].
Figure 4 shows a crude verification of the deduced EEDF where the aggregated number and energy density for a few ROIs are plotted against the same calculated using Equation (
2)—the degree of match can be appreciated.
3. Energy Dispersive Soft X-ray Spectroscopy
ECR plasmas emit radiation across the EM spectrum, with each type of radiation reflecting the properties of some element of the plasma. As such, soft X-ray emission spectroscopy is a well-established tool for studying intermediate energy electrons [
7,
8,
17,
18,
19] because the energy of the emitted photons lies in the range 2–30
. The spectra include both discrete fluorescence line emissions arising from the ionisation of the atoms, and continuous bremsstrahlung from the deceleration of electrons in the Coulomb field of the ions. Additionally, the experimental setup can be modified to study both global properties of the electrons (through volumetric spectroscopy) and local structure (through 2D space-resolved images).
The present work focuses on the analysis of soft X-rays from an Ar ECR plasma measured during an experiments at ATOMKI, Debrecen in 2014 [
8], by developing a plasma emission model to explain the outputs from the silicon drift detector (SDD) and charge-coupled device (CCD) camera. The former was used for volumetric measurements alone, while the latter was intended for local structural studies by incorporating it into a more elaborate pinhole camera setup. The specific schematics for both setups are shown in
Figure 5. The plasma chamber was of the same dimensions as the simulation domain, i.e.,
.
The SDD was coupled to a collimator that filtered out all photons outside of a narrow near-axis zone, allowing determination of plasma properties in a highly localised region. The CCD-pinhole setup could be operated in two different ways. The first mode of operation was the spectrally integrated mode wherein the CCD chip was exposed to the plasma for tens of seconds and lost spectral resolution, but recorded the total photon energy impinging on it. This was useful for analysing the coarse shape and structure of the plasma as well as the local energy content.
The second mode of operation, and more relevant to present calculations, was the spectrally resolved or photon counting mode. Here, the CCD chip functioned as a fast camera, capturing a large number of
-duration frames. This allowed each pixel to retain information about individual photons, resulting in a full spatially
and spectrally resolved X-ray map showing the finer details about the plasma structure and local sources of fluorescence. More details about the aforementioned modes of operation and their importance in plasma studies can be found in [
8].
4. The 2D Space-Resolved X-ray Imaging
In order to benchmark the theoretical model of warm electrons detailed in
Section 2, the CCD-pinhole setup was operated in photon counting mode for an Ar plasma heated at operating frequency of
and with RF power of
. A total of 2000 frames were captured, each of duration
.
Experimental validation was attempted by comparing a theoretical fluorescence emission model based on the simulated
map and Equation (
2) with the Ar fluorescence map. The calculation was rather straightforward. The starting point was the volumetric reaction rate (in units
) defined as
where
is the electron density,
is the Ar ion density,
E is the collision energy in centre of mass (CM) frame of reference,
is the K-shell ionisation cross-section as a function of collision energy,
is the collision speed, and
I is the K-shell binding energy (
for Ar). The volumetric reaction rate could be converted into an emission map by calculating the measured photon counts from each cell (dimensionless units), given by the simple expression
where 300 is the total exposure time in
,
Y is the fluorescence yield describing the percentage of K-shell ionisation converted into K
photons,
is the geometrical efficiency of the detector setup, and
is the quantum efficiency of the CCD chip for detecting a
photon.
V is the volume of each cell which in this case was simply
. The electron density was nothing but
.
While in principle Equation (
4) fully described the spatially-resolved Ar K
emission map, there were some missing quantities and uncertainties. The ion density
(and its spatial distribution) was unknown, but as a first approximation, it was taken as
assuming the plasma to be locally neutral and the ions to be in
state. There was also uncertainty in the contribution of the warm electrons to the total emission, which could manifest as a mismatch between predicted and measured photon counts in the X-ray images. Furthermore, finally, due to the peculiar configuration of the imaging setup and the finite size of the plasma chamber,
was expected to be position-dependent as well, which entailed a more rigorous calculation of the spatially-resolved geometrical efficiency. These missing aspects will be addressed in the following subsections.
4.1. Assessing the Warm Electron Contribution—EEDF Integrated Cross Section
To understand better which electrons contribute most to the photon emission, the most straightforward way was calculation of the EEDF-averaged cross section defined as
where
is the EEDF. Assessing the expected contribution of warm electrons from among the intermediate energy electrons is an important task—if hotter electrons with higher temperatures constitute the chief source of radiation, using only warm electrons in the theoretical emission model would eventually lead to underestimation of photon counts, or a corresponding overestimation in
if theoretical and experimental maps were scaled to same order of magnitude. This analysis has already been detailed in [
15], so only a brief overview will be provided here.
To evaluate Equation (
5), the semi-empirical Lotz formula was used as the Ar ionisation cross Section [
20] while the EEDF was taken as a pure Maxwell distribution. The upper limit of the integral in Equation (
5) was truncated to
for sake of brevity.
Figure 6a shows an overlap plot of the Lotz cross section vs. Maxwell EEDFs with different temperatures, while
Figure 6b shows the EEDF-integrated cross section. It can be easily concluded that the overlap between warm electrons as described in
Section 2 and the cross section is weak as compared to slightly hotter electrons with
and thus the emission map of Equation (
4) will likely predict fewer photons than seen experimentally.
4.2. Evaluation of Local Geometrical Efficiency (LGE)
When dealing with isotropic radiation from a point-sized source, the geometrical efficiency of the setup is given by the simple expression
where
is the solid-angle subtended by the detector on the source and
is the full emission solid angle. This formalism cannot be adopted for the present case because the (1) the plasma is not one single point-sized source but rather a collection of point-sources of finite size, (2) the solid angle of emission is subtended by the
source on the detector and not the other way round, and (3) each cell (point source) behaves differently with regards to
as defined by the optics of the cell and the pinhole setup.
Figure 7 summarises the key ideas presented in these points, while imperceptible in the schematic, the emission space and detection probability of a photon varies spatially with the position of the source cell in the plasma.
To address the aforementioned issues, a ray-tracing Monte Carlo method was used to calculate the geometrical efficiency of each cell in a cuboidal simulation domain of the same format as the 3D density and energy density matrices, thus called the
local geometrical efficiency (LGE). The strategy was quite simple and based on a two-step approach. First, assuming each cell to be an independent source, the limits of the polar and azimuthal angles were calculated such that the photon when emitted within this defined emission space, would surely pass through the pinhole. Normally this would be an unnecessary step and one could straight away proceed with simulating a
N photons emitted isotropically from each cell and check how many passed through the pinhole. This is a computationally expensive procedure here since the pinhole blocks a majority of the photons emitted, and thus the resultant emission space is severely limited. Knowing the minimum and maximum of the polar/azimuthal angles, the emission space could be calculated as
and the
ideal LGE would then be given by Equation (
6).
Figure 8a shows a rough sketch for calculating the emission space and
Figure 8b shows a 3D sliced view of
for each plasma cell. The subtle variation across the simulation domain can be appreciated.
However, simply passing through the pinhole is not a guarantee of detection—the photon could very well be hindered by some absorbing element or simply not intercepted by the CCD chip. To account for this,
N photons
within the above calculated emission space were simulated, and their trajectories traced along their journey through the setup. The number of photons
which made it till the CCD were checked, and a correction
was factored into Equation (
6) to obtain the final LGE.
Figure 9a demonstrates the idea of ray-tracing while
Figure 9b shows the final LGE. The setup of the experiment was modelled based on data provided in [
18]—the
diameter pinhole was drilled into a stack of tungsten and lead plates of thickness 1 and
, respectively, and an aluminium window of thickness
was placed after the pinhole to block UV photons. Together with the mesh at the end of the plasma chamber, they acted as X-ray absorbers. The impact of the ray-tracing correction can be immediately confirmed by looking at the dark band of cells at both edges of the simulation domain—these are the cells from which emitted photons pass through the pinhole but never make it to the CCD. The lack of cylindrical symmetry is due to a combination of factors like assumed cuboidal shape of the simulation domain and limited resolution in terms of cell size. Furthermore, it is important to note the equalising effect of the absorbing elements on the LGE (since the final
looks more uniformly distributed as compared to
Figure 8) and this is an artefact of the low statistics in the Monte Carlo routine. To retain the original distribution, a minimum of 10,000 photons should be simulated but currently only 500 were simulated. This does not have a large impact on the absolute results, however.
4.3. Experimental Benchmarking
Using Equation (
4), Lotz formula for Ar K-shell ionisation cross section, a fluorescence yield
,
and the calculated LGE, a theoretical emission map was obtained. The integral was evaluated with a Trapezoidal integration method and truncating the upper limit to
(the EEDF, and thus the integral, practically vanished after this limit). The global electron density was crudely re-scaled to
, just so a match could be made with experimental data.
Figure 10a shows the sliced emission map and
Figure 10b shows the more relevant longitudinally integrated image which should be compared to the real map in
Figure 10c.
Comparing the theoretical and experimental maps, it can be seen that the rough shape and structure of the plasma have been reproduced, as well as the “hole” in the near-axis region also seen in other experiments. There are of course differences, which can be attributed to a number of reasons like incorrect assumptions in ion distribution (they may not be locally neutral with the electrons [
2]), negligence of photon scattering, and incomplete modelling of CCD readout. In addition, reconstruction of the number of counts came at the expense of re-scaling
from
to
which arises from evaluating the reaction rate using warm electrons alone. These issues, however, only seem to imply that our model is yet incomplete but does hold potential for studying spatial distribution of electrons in ECR plasmas. Corrections to the LGE evaluation scheme, improvement to ion dynamics simulation and a model of CCD action are still underway, while the problems associated with uncertainties in electron density/temperature will be addressed through volumetric spectroscopy.
5. Volumetric X-ray Spectroscopy
To resolve the uncertainties about which electrons constitute the bulk of the fluorescence emission, analysis of the volumetric spectrum was done. The experiment was performed on the same Ar plasma heated with
RF at
power. The setup was as shown in
Figure 5b and additional details about the measurement scheme/results can be found in [
19]. Since many details about the the present analysis have already been published in [
15,
19], only a broad overview will be given here.
The bare spectrum from the SDD/ADC was first calibrated with Fe lines, corrected for quantum efficiency (QE) and dead time, and then converted to emissivity density . The presence of K and K emissions from Ar, Cr and Fe was confirmed from the experimental emissivity density.
The emission model was built on the premise that a single component Maxwell distribution
represented the entire emission zone. Comparison of model-predicted results with experimental data would validate the assumption.
was composed of three different components—bremsstrahlung from confined plasma electrons, fluorescence from plasma ions (Ar lines), and fluorescence from extraction plate atoms generated by escaping electrons (Cr and Fe lines). The bremsstrahlung emissivity density
equation was taken directly from [
19], expressed as
Here, is the electron mass, is the fine structure constant, is the photon energy, is the atomic number of Ar, is the temperature parameter of the Maxwell EEDF, and is the product of the electron and ion number density.
The Ar K
and K
line emissivity densities were calculated according to the expression
where
is the fluorescence factor associated with the transition,
is the energy per channel,
is the K-shell binding energy of Ar, and
is the semi-empirical Lotz formula. The fluorescence factors are the same as fluorescence yield
Y in Equation (
4) and
are connected to each other through the expression
where
is the line intensity ratio.
Finally, the emission from the extraction plate was calculated as
where
X represents either Cr or Fe,
is the loss electron density,
is the plasma volume introduced to convert the
total extraction plate fluorescence emissivity density into a volume-averaged value,
is the target atom number density, and
A is the area of the extraction plate interacting with the escaping electrons calculated as
,
l being the separation between the extraction plate and detection cone vertex and
d the diameter of the extraction hole. This area was a a rather thin ring-like shape intercepted by the collimator on the extraction plate.
Figure 11 shows the schematic for calculating the area. The quantity
represents the electron stopping power in AISI 316 steel which constituted the extraction plate, and was inserted to account for the constant modification of the ionising power of lost electrons as they traversed the target material.
For
, the Lotz cross section was replaced with the Deutsch-Märk formalism [
21]. Analogous to
Section 4.1, the relative contribution of warm and hot electrons was also checked for Cr and Fe by plotting the overlap plot of Deutsch-Märk cross section with Maxwell EEDFs of different temperatures, as well as the EEDF-integrated values. This is also shown in
Figure 6 and here too electron population with
look to contribute more strongly. These results align well with calculations done in [
19] on bremsstrahlung from
Figure 12 alone, where electron temperature was estimated around
.
The fluorescence emissivity density from Equations (
9) and (
11) were modulated with a Gaussian profile to account for the line broadening and added to Equation (
8) to generate the full model comprising of all three components. This was fit to the experimental spectrum using a Trust Region Selective least squares fitting routine, and parameters of interest, namely
,
and
were estimated as
,
and
, respectively. The result of the fit is shown in
Figure 12.
The FWHM of the fluorescence peaks was estimated in the range 0.117–0.212
and the escaping electron current as calculated from
was around 2–5
which is the same order of magnitude as extracted ion current [
19]. Just as importantly, the combined charge particle density
implies
(if the ions are in
charge state) which is a fully valid result. This serves to prove that volumetric soft-X-ray spectroscopy is a powerful method to probe densities and temperatures of relevant electron populations and by coupling it to the 2D X-ray imaging analysis described in
Section 4, the plasma can be completely characterised. For the moment, what it tells us is that warm electrons with in energy range 0–30
intervals are not sufficient to reproduce the degree of ionisation present in the plasma, and thus hotter species with
should be simulated to match the experimentally obtained maps, both in terms of structure and absolute photon counts.