On the Temperature and Plasma Distribution of an Inductively Driven Xe-I 2 -Discharge

: Inductively Coupled Plasma (ICP) discharges are part of intense research. Predicting different plasma parameters, like the distribution and temperature of the present species, is of great interest for many applications. Iodine-or halide-containing plasmas in particular have an important function, for example, in the development of mercury-free UV radiation sources. Therefore, a 2D simulation model of a xenon-and iodine-containing ICP was created by using the Finite Element Method (FEM) software COMSOL Multiphysics ® . The included species and the used reactions are presented in this paper. To verify the simulation in relation to the plasma distribution, the results were compared with measurements from literature. The temperature of the lamp vessel was measured in relation to the temperature distribution and also compared with the results of the simulation. It could be shown that the simulation reproduces the plasma distribution with a maximal deviation of ≈ 6.5% to the measured values and that the temperature distribution in the examined area can be predicted with deviations of up to ≈ 24% for long vessel dimensions and ≈ 3% for shorter dimensions. However, despite the deviating absolute values, the general plasma behaviour is reproduced by the simulation. The simulation thus offers a fast and cost-effective method to estimate an effective geometrical range of iodine-containing ICPs.


Introduction
In addition to use in plasma thrusters, a main application for iodine-or halidecontaining plasmas are mercury-free radiation sources [1,2].Besides the generation of radiation by using barrier discharges, ICPs for radiation generation were also investigated [3,4].Because of certain properties, such as the vapour pressure, halide compounds in particular are suitable as a substitute for mercury [5].Radiation sources containing halides are, therefore, still part of recent research [6].An electrodeless ICP operation of the lamps provide access to use the full bandwidth of halide compounds as lamp fillings, without the need to consider electrode chemistry as life-limiting process [7].
However, the ignition and operation of inductive mercury-free radiation sources can be difficult [8].Furthermore, the plasma and temperature distribution have important effects on the discharge.To achieve efficiencies similar to conventional mercury-containing systems, the lamp pressures must be in a comparable range [9].The lamp pressure is directly dependent on the temperature of the coldest point of the lamp vessel.As the temperature distribution is related to the plasma distribution, this correlation will be considered in the following.
A direct comparison of a pure Xe discharge with a discharge containing a Xe-I 2 -mixture shows a difference in the plasma distribution even from the optical appearance.Therefore, a simulation model was developed to investigate this behaviour in more detail.For the simulation, the Finite Element Method (FEM) software COMSOL ® was used [10].This software includes a module for ICP simulations.The simulation results were compared with the work of Barnes et al. [11] in relation to the plasma distribution.Furthermore, the results of the temperature distribution for the glass vessel were also compared with measurements on a lab system.The simulation model presented in the following shows a useful approximation of the plasma behaviour for future research.

Materials and Methods
The plasma model is represented by highly nonlinear coupled systems of differential equations.In order to approximate the plasma behaviour, it is possible to consider the system as a fluid model.This approach will be described more in detail below, as the used software applies the fluid model for simulation [12].In the fluid model, the plasma is described as a fluid by averaged macroscopic quantities such as density, velocity and energy of the included species.The solution for the macroscopic quantities results from the continuity, momentum, and energy conservation equations.The result are the moment equations of the Boltzmann transport equation [13].The Boltzmann transport equation assigns a distribution function f b = f b (r, t, ν) of the respective particle density to each species.This depends on the location r, the time t, as well as the velocity ν and is also dependent to the source or impact term S( f b ).Thus, the transport equation is given by Here, m corresponds to the particle mass and F to the external force acting on the particles, for example, through electromagnetic fields.If the equation is now integrated over the velocity space, the continuity equation already mentioned results in The number of changes due to collisions (e.g., ionization, recombination, etc.) as a function of time and place is given here by S and → Γ describes the particle flux.This can be given with help of the drift-diffusion approximation: To determine the energy conservation equation, the transport equation is multiplied by 1  2 mν 2 and integrated over the velocity space.For the use of the energy, this leads to Analogous to the particle flux, the electron energy flux → Γ can be calculated: where all quantities are related to energy.Consequently, n is the electron energy density, µ is the electron energy mobility and D is the corresponding diffusion coefficient of the electron energy.In the software, the quantities are related to the mean electron energy ¯ as follows [12]: n = n e ¯ ( 8) Here, T e describes the electron temperature and µ e the mobility of the electrons.As part of the electron transport calculation, Equation ( 2) is calculated by the software as follows: with the vector of fluid velocity u, and the electron flux: The source term of the electrons R e is obtained for all reactions j from the sum of their rate coefficients k j , their mass fraction x j , as well as the neutral particle density n n and the electron density n e : For this purpose, the rate coefficients can be calculated from the energy-dependent impact cross sections σ( ).COMSOL Multyphysics ® uses the approach: The energy-dependent impact cross sections σ( ) are available for many materials and reactions in various databases.The correspondence of the energy conservation Equation (4) for the electron energy density is given in the software as Here, Q and Q gen correspond to the energy input from an external and a general heat source.S n describes the energy gain or loss due to inelastic collisions.S n is formed by the sum of the energy loss ∆ j due to collisions occurring over all reactions j: It is obvious that the rate coefficients k j and thus the energy-dependent impact cross sections σ( ) play a decisive role in the simulation of the system.Even in a pure Xe-plasma, several hundreds to thousand reactions already take place, some simplifications have to be made.To limit the complexity of the model and still achieve a good approximation of the plasma behaviour, the modeling was initially limited to the most essential reactions.The reactions used are listed in Table 1.The reactions are linked in Table 1 with their data sources [14][15][16][17][18][19][20][21][22].It has to be noted that for a more compact nomenclature of the excited states of Xe the notation of Sommerer [23] was used.
With the used software it was not possible to include the evaporation of the halide in this simulation.For this reason, it was assumed that the material already is completely in the gas phase at the time of ignition.For this purpose, the I 2 -vapor pressure curve [24] as well as the lamp volume and the I 2 -filling quantity were used to determine the temperature at complete evaporation.The underlying geometries and filling parameters here are listed in Table 2. Subsequently, the pressure difference ∆p to the real initial pressure was determined with the temperature and the particle number of the starting gas as well as the I 2 -filling quantity according to the ideal gas law.This was used as initial pressure to simulate the discharge.An overview of the pressure curves obtained in this way can be seen in Figure 1.By means of this procedure, it is now possible to simulate the parameters of the discharge numerically without including the evaporation.The disadvantage of this approach is that the simulation can only be expected to provide reliable approximations of the plasma behaviour from the point at which an unsaturated operation is reached.Nevertheless, the simulation model provides a good estimation of the plasma parameters, as described below.The green marked area corresponds to the area of an unsaturated discharge.

Results and Discussion
The following section is divided into two subsections.In the first subsection, the results are presented and discussed in relation to the plasma distribution.The second subsection describes the effect of the plasma distribution on the temperature distribution on the glass vessel as well as gas temperature distribution in the plasma.

Plasma Distribution
By directly comparing an I 2 -containing Xe discharge to a pure Xe discharge, several differences are noticeable.As shown in Figure 2, the electron density behaves clearly different.For example, the electron density in the Xe discharge at 400 W input power reaches a mean value of n e = 6.8 • 10 14 cm −3 .The I 2 -containing discharge has a lower mean electron density with only n e = 4.8 • 10 10 cm −3 .This significant deviation at the same input power, geometry, and pressure can be justified by reaction 10 in Table 1.This reaction is due to the high electronegativity of iodine and the resulting possibility of dissociative attachment of electrons to form anions. Through this process, the electron density is significantly reduced compared to the pure Xe plasma.In the shown case, the mean electron density is reduced by approximately 4 orders of magnitude.A shift in the electron density distribution in contrast to the pure Xe discharge is also clear in Figure 2.This can be explained by the previously mentioned process too.The different behaviour becomes particularly clear by comparing the electron density distribution with the I − -and the I-distribution.Figure 3 shows the mentioned distributions.The I − distribution shows the increase in particle density on the edge regions.If the radial I 2 distribution is considered, the highest concentration is also reached on the edge regions.This behaviour in I 2 -containing inductive discharges is known and based on the recombination of dissociated iodine atoms at the wall [11].To verify the simulation, a model was created with the geometry and parameters as described by Barnes et al. in [11].A comparison of the radial I 2 -distribution between the simulation carried out and the measurements of is shown in Figure 4.The irregularly distributed simulated points in Figure 4 are based on the meshing of the model.The figure shows that the simulation approximates the measurements well.To determine the deviation, the measured values were interpolated and the maximum deviation was determined.The maximum deviation results in the curve slope with ≈6.5% to the measured values.A comparison with the simulation results of the geometry and parameters used in this paper is also shown.The results indicate that the general plasma behaviour with an increasing I 2 -density at the edges of the lamp vessel can also be assumed here.The described and simulated behaviour has direct effects on the plasma.This is particularly evident in the temperature distribution, which is explained more in detail in the next section.

Temperature Distribution
In relation to the temperature distribution in the lamp systems, the dimension of the lamp vessels play a decisive role.With a given frequency, power, pressure, and coil, only the geometry remains as a parameter to influence the temperature distribution.In lamps containing halides, the filling components condense at the coldest point of the system.To generate medium to high pressures, the aim is to create a temperature distribution that is as homogeneous as possible.However, this is only partly achievable due to the behaviour of the plasma.It has been observed that in the case of halide-containing discharges, the plasma tends to form a sphere which can already be seen in Figure 5. Here, the temperature behaviour is already visually observable by the plasma distribution.This behaviour means that the coldest point is always at the ends of the lamp due to the plasma behaviour.As the hottest point is therefore in the middle of the lamp, the coldest point is also determined by the lamp length.Therefore, temperature measurements were carried out at different lamp lengths.A thermographic camera was used to measure the lamp temperature (A325, FLIR Systems,Wilsonville, OR, USA).This method allows to monitor the temperature on the entire surface of the lamp vessel and to identify the hottest and the coldest point on the surface.The values used were measured after thermal stabilization of the lamp.These measurements can be seen in Figure 6.Note that with this method only the temperatures on the outer glass vessel can be recorded.
The measurement shows that the behaviour has a significant influence on the temperature distribution.From the hottest point in the center of the lamp, the temperature drops significantly towards the ends of the lamp.Despite the different lengths, the lamp bodies have a similar temperature distribution.For illustration, the quotient T q of the measured maximum temperature T max and the minimum temperature T min is compared.T q results as In order to achieve a homogeneous temperature distribution on the outer glass vessel, a geometry should be chosen, so that it has the smallest possible value of T q .Since the temperature decreases from the center along the length of the lamp, the shortest lamp length results in the smallest temperature ratio from the hottest to the coldest point.In order to verify the simulation model, the simulated temperature quotient was compared to the measured one.The results are shown in Figure 7.The comparison shows a large deviation of the simulation from the measurement at longer lamp lengths with a maximal deviation of ≈24%.For smaller lamp lengths, deviations of ≈3% are achieved.From the simulation and the measurement, it can be seen that the temperature quotient T q approaches a homogeneous distribution for shorter lamp lengths.As shown in Figure 8, the lamp length has only a small influence on the temperature distribution in the plasma.
The simulation shows that even with longer lamp geometries, the plasma is still concentrated in the center of the lamp, this behaviour considers to be consistent with the results in Figure 5.This means that it is not possible to select a lamp geometry of arbitrary length if a homogeneous temperature distribution over the entire glass vessel of the plasma is desired.

Conclusions
In this work, a 2D model for the simulation of Xe-I 2 plasmas was presented.The simulation results were compared with other publications in relation to the plasma distribution.The temperature distribution was also simulated and compared with own experimental results.It could be shown that the simulation can reproduce the measurements from the literature in terms of the plasma distribution with a maximum deviation of ≈6.5%.
The results suggest that also other plasma parameters such as the plasma impedance, which is influenced by the contained species, could be determined by the simulation with small deviations as well.Therefore, the simulation could be a useful tool for designing appropriate power supplies.
The simulation of the temperature distribution deviates significantly from the measured values for longer geometries.For small lamp lengths, the simulation results match the measurements with a deviation of ≈3%.A precise consideration of the validity limits beyond the areas analyzed needs further investigation.Nevertheless, the results show that the simulation reproduces the general temperature behaviour and could be used to estimate reasonable lamp geometries.

Figure 1 .
Figure 1.Pressure curve over the temperature of a discharge with 100 Pa of Xe as well as 1 mg of I 2 .The green marked area corresponds to the area of an unsaturated discharge.

Xe+I 2 :Figure 2 .
Figure 2. Comparison of the electron density distribution at 400 W input power with 3 MHz and an initial pressure of 190 Pa.Left: Xe discharge.Right: I 2 -Xe-discharge.

Figure 3 .
Figure 3.Comparison of different particle density distributions of the I 2 -Xe plasma at 400 W input power.Left: Electron density.Center: Atomic iodine density.Right: I − -density.

Figure 4 .
Figure 4.I 2 -particle density distribution along the lamp radius in the middle of the vessel.Comparison of the measurements of Barnes et al.[11], indicated with n B m ; a simulation using the same parameters and geometry, indicated with n B s ; and the simulation with the geometry and parameters presented in this paper, indicated as n s .

Figure 5 .
Figure 5.Comparison of the plasma distribution of pure gas-to halide-filled lamp systems at 400 W input power.Left: Xe-filled lamp system.Right: Xe-I 2 -filled lamp system.

Figure 6 .
Figure 6.Measurement of the temperature distributions for different lamp lengths.The zero was set at the hottest point.

Figure 7 .
Figure 7.Comparison of the simulation and measurements of the temperature quotient T q .

Figure 8 .
Figure 8.Comparison of plasma temperature distribution with different lamp lengths.Here, a coil with 15 windings was simulated to reproduce the real system for the measurement.

Table 1 .
Overview of the reactions and data sources used for the simulation.For this reaction the reaction rate of k j = 2.73 • 10 8 s −1 was used.

Table 2 .
Lamp and coil geometry for the I 2 simulation.