`
`
`A
`applicabletto argcndischarges over at
`
`
`hollow-
`ll : Application;t10a
`cathodearc‘andslowvpressureglowe V
`dischafgesrr
`
`
`
`
`
`
`
`
`J Vléekt and V Pelikéni.
`T Department of Physics, Institute of Mechanical and Electrical Engineering.
`30614 Plzen. Nejedlého sady 14. Czechoslovakia
`3; Department of Automatic Control Systems. Skoda National Corporation, 31600
`Plzen, Czechoslovakia
`
`Received 12 July 1988, in final form 28 November 1988
`
`Abstract. The extensive collislonal—radiative model for an argon atom plasma is
`applied to a low-pressure, hollow—cathode arc discharge and to the positive column
`of a low-pressure glow discharge in order to clarify the mechanisms by which the
`excited levels in these discharges are populated, the results being compared with
`experimental investigations in the literature. Computations are carried out for
`various sets of input parameters. such as the electron kinetic temperature T5, the
`atom temperature T3. the ion temperature T, the electron number density ne, the
`ground state atom population n-_ the plasma column radius R and the escape
`factors Am and
`characterising the non-equilibrium plasmas under
`consideration. The predictions of our model. i.e. the populations in the excited
`levels as a function of the electron number density, the effective principal quantum
`number and the discharge current, are compared with the experimental results and
`in two cases also with the theoretical results of other authors. it is shown that all
`calculated dependences are fairly close to the corresponding experimental curves
`referring to both discharges. The results presented confirm the applicability of the
`so—called ‘analytical top model’ of van der Mullen er al and Walsh's formula for AV
`interpreted according to Mills and Hleftje.
`
`1. introduction
`
`In the preceding paper (Vlcek 1989). a collisional-
`radiative (CR) model applicable over a wider range of
`conditions than those described in the literature (Gian-
`naris and Incropera 1973. Katsonis 1976. Gomés 1983.
`van der Sijde er al 1984a. I-Iasegawa and Haraguchi
`1985) was established for an argon atom plasma.
`Atorn—atom inelastic collisions and diffusion losses
`of the metastable states.
`together with the electron—
`atom inelastic collisions and radiative processes usually
`included. are considered in this model.
`taking into
`account 65 effective levels.
`
`Analytical expressions used for the corresponding
`cross sections are in good agreement with currently
`available experimental and theoretical data for argon.
`In particular. the measurements of Chutjian and Cart-
`
`O022-3727x89/O50682 -*- 12 $02.50
`
`1989 IOP Publishing Ltd
`
`wright (1981) and the computations of Kimura er al
`(1985) have been used extensively.
`With the help of the previously established numeri-
`cal method (Vléek 1989). we can calculate the popu-
`lation coefficients determining the populations in all
`excited effective levels. This enables us to study the
`mechanisms by which these levels are populated under
`various conditions in a non-equilibrium argon plasma
`characterised (even in the case of a non-Maxwcllian
`electron energy distribution function) (EEDF) by the set
`of parameters TC. Tu. Ti.
`/23. 111. R. A,_,,,, and A,,,.
`The main aim of the present paper is to check
`the reliability of our extensive CR model under the
`conditions in the low—pressure hollow cathode arc.
`studied experimentally by van der Mullen et a[ (1978.
`1980) and by van der Sijde er (2! (1984a.b). and in the
`positive column of the low-pressure glow discharge.
`
`GlobalFoundries 1112
`
`GlobalFoundries 1112
`
`
`
`investigated experimentally by Kagan er al (1963b).
`A further motivation for this study is our interest in
`understanding the mechanisms by which the excited
`levels are populated in these two argon discharges of
`practical interest.
`For experimental verification of our CR model. a
`combination of these discharges seems to be suitable
`because the population mechanisms differ appreciably
`in them. Moreover, because the electron energy distri-
`bution function in the hollow-cathode are investigated
`is Maxwellian (Pots 1979) and the radiation trapping
`is negligibly small (van der Sijde er al 1984a,b), the
`calculated excited level populations obtained under
`these conditions cannot be affected by possible inac-
`curacies arising from the solution of the Boltzmann
`equation and from determination of the optical escape
`factors.
`
`2. Results and discussion
`
`Assuming that the quasi-stationary state model can be
`applied (see for example Cacciatore er al 1976, Biber-
`man er al 1982) and the fundamental mode of diffusion
`, of the metastables to the wall is dominant in the dis-
`charge tube (see for example Delcroix er :2! 1976, Fer-
`reira er al 1985). we obtain a set of coupled linear
`equations
`
`n=-
`2 any/inn :
`
`_ anzinl
`
`. 65. from which the unknown excited
`.
`.
`where m = 2, .
`level populations 11,, may be calculated. provided that
`the coefficients am and (Sm are known (Vlcek 1989)
`and the ground state atom population 111 has been
`determined experimentally.
`Owing to the possibility of investigating the effect
`of the upward ionisation flow of electrons from the
`ground state atom and their downward recombination
`flow from a continuum on the populating of the excited
`levels, the system (1) is solved. in spite of the fact that
`71; is not an independent parameter in our case. in the
`standard form
`
`12,; =21?’ —.— G_l,”r1;
`
`f0rIz=2.....65
`
`(2)
`
`where the population coefficients 11$” and G3,” represent
`the solutions of (1) with 11, = 0 or r11 = 1 and a,,, = 0.
`respectively inserted into their right hand sides.
`The expression (2) can be rewritten as
`
`77.»: = riff’-’n?. + rl'7’It§
`
`(3)
`
`where r2,5, and 21,‘? are the corresponding Saha population
`and Boltzmann population, respectively,
`r‘,,°‘ and rf,"
`are the so-called CR coefficients relating the actual popu-
`lations n,, to Iii and 223?. respectively.
`p
`In a general case. the numerical method developed
`allows us to calculate the population coefficients nff“
`and G‘,,” as functions of the following input parameters:
`Tc. Ta. Ti. 11:. 211. R. A,,,,, and Am with the possibility of
`neglecting the atom—atom inelastic collisions when their
`
`Collisional—radiative model for argon discharges, ll
`
`Table 1. Data characterising the excited effective levels
`including all actual 4s and 4p states.
`
`Statistical
`Excitation
`Level
`number
`Designation
`energy
`weight
`
`[7
`npqn[[KlJ
`51:1 (evl
`9n
`
`5
`1 1 .548
`4s[3.«-‘2],a_
`2
`3
`11 .624
`4s{3.«"2]1
`3
`1
`11.723
`Ash .1210
`4
`3
`11.828
`4s’{1/2],
`5
`3
`12.907
`4pl1 /'2]1
`6
`20
`13.116
`4p{3/21,2, [5/2}“
`7
`8
`13.295
`4p'[3./211,2
`8
`3
`13.328
`4p'[l/2]1
`9
`1
`13.273
`4p[‘l /210
`10
`
`11 1 4p’[1 /219 13.480
`
`
`
`influence on the population mechanism is-studied.
`Under
`a
`reasonable assumption that only the
`reabsorption of the resonance radiation may be impor-
`tant in the discharges investigated. we have used the
`analytical formulae for the escape factors (Mills and
`Hieftje 1984) in which A1,, are dependent only on T,,,
`n-_ and R, where rt = 3. 5,12,15,16.17. 20.21.26, 27
`and 33 (see table 1 of Vlcek 1989).
`In the special case. when the Boltzmann equation
`need not be solved due to the Maxwellian form of the
`
`EEDF in a plasma. computations become straightforward
`and relatively rapid. Furthermore. the ion temperature
`T, does not then appear among the input parameters
`(Vlcek and Pelikan 1985. Vlcek 1989).
`The basic data characterising the excited levels con-
`sidered in our CR model, including all individual '-is and
`4p states, are given in table 1 to help in comparing the
`calculated and experimental results.
`
`2.1. Hollow-cathode arc discharge
`
`Hollow-cathode arc discharges have been widely used.
`for example in the investigation of argon ion laser
`plasmas (Pots at (111978). in fusion-oriented and welding
`technology (Chall and Uhlenbusch 1982).
`in plasma
`centrifuges (Wijnakker er al 1979). or generally in
`plasmas in magnetic fields (Boeschoten er al 1979).
`Valuable results contributing to the elucidation of
`the mechanisms by which the excited levels in a low-
`pressure, hollow-cathode arc are populated have been
`obtained in recent years by van der Mullen er al (1978.
`1980, 1983) and van der Sijde er al (_l984a. b).
`In our case. the numerical results for the excited
`
`level populations n,, and for a CR coefficient rf,” are
`respectively compared with the corresponding values
`determined experimentally by van der Mullen ez al
`(1978. 1980) and van der Sijde er al
`(198421. b) who
`measured the absolute intensities of numerous Ar 1 lines
`as a function of the electron number density 21_. for
`known values of Tc. T3, In and R in the highly ionised.
`magnetically confined plasma of a low-pressure. hollow-
`cathode arc.
`
`It has previously been verified that the EEDF is
`
`633
`
`
`
`“*4
`
`l1l
`l1
`
`
`
`J Vléek and V Pelikan
`
`‘1
`
`’1
`
`ll
`
`Flgure 1. The excited level populations nn related to the corresponding Sana
`values nfi for various electron number densities ne at TE = 40600 K, Ta =
`11600 K, Fl = 1 cm and n, ranging from 2.35 x 10‘3crn“3 tor ne =
`5 x 10’? cm‘3 to 3.62 x 10” cm‘3 for ne = 10‘5cm‘3, Present work: curve A,
`n=6; curve 8, n= 11; curve C, n=7: curve D, n: 10; curve E, n=8;
`curve F. n = 9 (see table 1); curves G and H, the 5d and 6d groups,
`respectively, with jc = % (full curves) and jc = § (broken curves). Symbols 0,
`A. 0 show the experimental data of van der Mullen er al (1980) for the
`whole 4p. 5d and 6d groups, respectively.
`
`Maxwellian (Pots 1979) and that doubly ionised ions
`can be neglected (van der Mullen er al 1980, Pots 1979)
`under the conditions considered.
`
`In figure 1 we compare the ratios n,,/n;‘; measured
`by van der Mullen er al (1980) for the 4p, 5d and 6d
`groups with the corresponding Values calculated by us
`at the plasma parameters Ta, 218 and n1 presented in the
`above—mentioned paper and for Ta = 11600 K and R =
`
`1cm taken from similar experiments (van der Sijcle er
`al 198421) carried out with the help of the same set-up.
`According to van der Mullen et al (1978), the accuracy
`in the measurement of n,, is estimated to be 50%.
`As can be seen in figure 1, the uniform decrease of
`the calculated values of 22,,/n,§ as rz;3 in the so—called
`‘excitation saturation phase’ agrees well with the
`measured dependences, but
`the theoretical values.
`
`634
`
`
`
`Collisional—radiative model for argon discharges. ll
`
`coefficient as a function of the electron number density ne. Theoretical results: curve A, Pots (1979);
`Figure 2. The
`curve B, van der Mullen et al (1977); curve C, van der Sijde et al (1984a) (all with T9 = 34800 K); curves D,
`determined by us using the model of Katsonis (1976) at T5 = 35 000 K (full curve) and Te = 40000 K (broken curve);
`present work: curve A‘, n., = 3 x tO‘~2cm‘3, Ta = 3480 K, R = 1 cm. T9 = 34800 K (chain curve) and Te = 46400 K
`(dotted curve); B’, as for A’ but at n, = 9 x t0‘2cm“3; C’, n, = 3 X 10‘? cm‘3, Ta = 6960 K, R: 1 cm. Ta = 34800 K
`(full curve) and T9 = 48400 K (broken curve); D’. as for C’ but at n1 = 9 X 1O‘2cm‘3. A, experimental results (van der
`Sijde et al 1984b).
`
`those obtained for the individual 4p states
`except
`denoted by n = 6 and ll. are somewhat underestimated
`compared with the measurements. The values of He
`relating to the transition of the excited levels considered
`to the regime of the partial local thermodynamic equi-
`librium are in satisfactory agreement with those
`obtained by extrapolation of the measured data in the
`work of van der Mullen et al (1980).
`It is known (see for example van der Sijde et al
`1984b) that the values of the CR coefficients rff’, occur-
`ring in (3), may be important for spectroscopic diag-
`nostics of a plasma. e.g. for the determination of the
`electron temperature T: and the ground state atom
`population n(. if the nth level is in the complete satu-
`ration phase.
`In figure 2 we give our results for the coefficient
`r§,". where n = 9. as a function of ne for comparison
`with the corresponding experimental values obtained
`by van der Sijde et al (1984b) under the following
`conditions in the discharge: 34800 K $ Ti. S 46400 K,
`3480 K S Ta S 6960 K,
`10” cm‘-‘ s ne $10“ crn‘3,
`3 X 10” cm‘3 S rzl s 9 X 10” crn‘3. As van der Sijde et
`al (1984b) did not give the value for the plasma column
`radius in their work. we have used R = 1 cm again (van
`der Sijde er al 1984a) in our computations illustrating
`the dependence of the coefficient rfi,“ on the electron
`
`temperature T3 and also on the atom temperature Ta
`and the ground state atom population n-1. Assuming
`that the value of 711 changes with rig in a manner similar
`to the measurements of van der Mullen er al (1980), we
`have determined the most probable values of r‘,_" at nc =
`1013 cm’3 and 101"‘ cm'3. They‘ are denoted by vertical
`arrows in figure 2. Our calculations have proved that
`the plasma investigated is completely optically thin also
`for all radiative transitions excepting the resonance tran-
`sition from the 4s’[1_/2], state at
`211 = 9 X 10‘3cm’3
`where we have obtained the value of [\15 = 0.46 and
`0.83 for Ta = 3480 K and 6960 K. respectively.
`The numerical results obtained from several CR
`
`models under the assumption that the plasma is com-
`pletely optically thin are also shown in figure 2.
`When we used the extensive 65-level model of Kat-
`
`sonis (1976). in which all important excitations from the
`ground state of an atom to the levels lying above the
`45’ states are omitted, except the excitations to the
`effective levels with n = 15. 16 and 17 (Vlcek 1989). the
`model curves for r‘,," obtained at Te = 35000K and
`40000K are given. When the simplified models (van
`der Mullen et al 1977, Pots 1979) and the extensive 49—
`level model of van der Sijde er al (19843) have been
`applied at TC = 34800 K, the results for the 4p group as
`a whole are shown. Note that in the model of van der
`
`635
`
`
`
`J Vlcek and V Pelikan
`
`
`
`17,./n‘,",-1
`
`Figure 3. The population factor n_n/n;‘,3 — 1 as a function of the effective
`principal quantum number n;q,, obtained with the following input parameters: A,
`jg =% and A, jc =§tor Te = 58000 K, Ta = 11600 K, ne = 6.7 ><10‘3cm‘3,n.=
`10” cm”3 and R=1cm; V,/'6 =1; and V.jc =% for T9 = 34800 K, Ta:
`11600 K, ne = 6.7 x 1O‘3cm'3, n, = 10” cm'3 and F?=1cm.O, experimental
`results (van der Mullen et al 1980). The full line x = 6.0 and the broken line x =
`5.0 represent the slopes predicted on the basis of the analytical top model (van
`der Mullen er al 1983).
`
`Sijde er al (1984a) the semiempirical formulae of Vriens
`and Smeets (1980) proposed for neutral hydrogen and
`alkali excited states are employed, all 45 states are
`separated and the statistical weights of the f groups are
`increased artificially.
`Figure 3 shows our numerical results together with
`the corresponding experimental results (van der Mullen
`er a1 1980) for the value of rz,,/rzfi - 1 as a function
`of the effective principal quantum number
`n;"q,, =
`(ell/e,,)""3. where 8? and an are the ionisation energies
`for atomic hydrogen in the ground state and for the nth
`level of argon respectively.
`As can be seen. the measured values obtained at
`actual electron temperatures in the range from 34 800-
`
`636
`
`SSOOOK and at 213 = 6.7 X 1013 cm” are found to lie
`between the model curves determined by us at these
`two extreme values of Tc. Moreover. the slopes of the
`straight parts of our dependences agree well not only
`with the experimental result but also with that obtained
`on the basis of the so-called analytical top model (van
`der Mullen er al 1983). Using this model to describe the
`population mechanisms of high-lying excited levels in a
`real. collisionally dominated, ionising plasma with the
`Maxwellian EEDF, one can write the following simple
`power law for non-hydrogenic systems:
`
`“:2//n/Sr F1: bC*(n,fqzz,)"'~
`
`where bu is a constant and the value ofx is in the range
`
`
`
`Colllsiona|~radiative model for argon discharges. ll
`
`
`
`
`
`/1,,_/llff‘,1‘,‘,’n1(L'ln''1)
`
`
`
`
`
`(full curves),
`Figure 4. The calculated values of nl?’ (broken curves) and Gglln,
`togetherwith the measured values (0, O, A, A) of n,, (van der Mullen er al
`1978) as functions of ne at T5 = 34800 K, Ta = 11600 K, R1 = 10” cm‘3 and
`F? = 1 cm for selected excited levels (see table 1): curves A, n = 2; curves B,
`5d
`n = 3; curves C and A, n= 7; curves D and 0, 5p (jc =%); curves E and
`(/C = %); F and 0, 7s (jg = %); G. npqn =10 (jc = is). The dotted curve, Gf.,“n, for
`nm = .0 (jc = %) but with atom—atom inelastic collisions neglected.
`
`Table 2. The experimental discharge parameters (Kagan er al 1963b) for filling pressure of 5Torr together
`with the estimated values for the atom temperature Ta and the ground state atom population /1,.
`
`Glow
`discharge
`current
`I
`(mA)
`
`Electron
`number
`density
`ne
`(10“ cm‘3)
`
`Electron
`kinetic
`temperature
`7;
`(104 K)
`
`Second estimate
`First estimate
`—-———————-——--————-—-— —-——e——~—~—-——
`Atom
`Ground state
`Atom
`Ground state
`temperature
`population
`temperature
`pozpulatlon
`7;’
`ng‘
`Tif’
`Q
`(102 K)
`(10‘5 cm‘3)
`(102 K)
`(10‘5 cm'3)
`
`12.10
`4.00
`16.10
`3.00
`2.50
`2.70
`25
`10.70
`4.50
`15.10
`3.20
`2.20
`6.00
`50
`9.66
`5.00
`13.80
`3.50
`1.80
`12.00
`100
`200
`22.00
`1.60
`4.00
`12.10
`5.50
`8.78
`400
`46.00
`1.30
`5.00
`9.66
`6.50
`7.43
`
`
`637
`
`
`
`J Vléek and V Peiikén
`
`
`
`n,,(cm”’l
`
`m 8
`
`16
`
`2
`
`Figure 5. Populations in the 4s and 4p levels (see table 1) obtained with
`experimental discharge parameters (Kagan at al 1963b) p = 5Torr, R = 1.2 cm,
`I = 25 mA: curves A, theory; A, experimental results; I = 400 mA: curves B,
`theory; 0, experimental; full curve A. ne = 2.7 x 10“ cm‘3. Ta = 25000 K, Tg” =
`Ti“) = 300 K, nl,” = 1.61 X 10” cm“; broken curve A, T?’ = Tim = 400 K,
`ng2>=1.21 x 10” cm ‘3; chain curve B. ne = 4.6 x 1012 cm“3, Te =13ooo K,
`Tgl = Tf” = 500 K. nl," = 9.55 x 1015 cm‘3; dotted curve B, T9 = T5?’ =
`650 K. ngzl = 7.43 X 1O’5cm”3.
`
`from 5.5 to 6.5; the lowest possible value, x = 5.0,
`may be obtained if ionisation and recombination are
`neglected (van der Mullen et al 1983).
`In figure 4 we give the calculated values of nff” and
`G,‘,”tz1 and the measured values ofnn for selected excited
`levels as functions of the electron number density ne at
`the discharge parameters characterising the experiments
`of van der Mullen at al (1978).
`A comparison of the terms appearing in (2) helps us
`to understand the role played by the upward ionisation
`flow of electrons in the system of excited levels and
`their downward recombination flow in determining the
`excited level populations under the conditions con-
`sidered.
`
`As seen in figure 4, the indicated critical values of
`ne, for which both terms are balanced in (2), decrease
`rapidly with the decreasing ionisation energy of the
`excited levels investigated. When the electron number
`density is higher than these critical values. excited levels
`come into Saha equilibrium where rz,, is proportional to
`ng. As shown, the results presented are in better agree-
`ment with the measurements (van der Mullen er (211978)
`than in figure 1.
`The effect of the atom—atom inelastic collisions on
`the excited level populations has been found to be
`negligibly small for all sets of discharge parameters
`considered if neg/n,_ 2 1. At lower degrees of ionisation
`this effect is growing. e.g. at in = 10” cm“3 and with
`
`638
`
`
`
`Collisional—radiative model for argon discharges, H
`
`
`
`n,,inn"?J
`
`R l
`2
`
`l
`
`l
`1+
`
`l
`
`I
`
`l
`3
`
`i
`
`10
`
`Figure 6. Populations in the 4s and 4p levels obtained with experimental
`discharge parameters (Kagan eta! 1963b) p = 5Torr, F? = 1.2 cm, l= 50 mA:
`A, experimental results; full curve. theory with ne = 6.0 X 10“ cm“3, Te =
`_
`22000 K, Ta = 7} = 320 K and n, = 1.51 X 1O‘7cm‘3; broken curve, theory with
`ne = 5.4 x 10” cm’3, Ta = 22000 K, Ta = 7: = 320 K and 1.51 ><1O‘7cm‘3:
`dotted curve, theory with ne = 6.0 X 10” cm‘3, T9 = 20 700 K, Ta = T, = 320 K
`and n, = 1.51 X 10” cm‘3; I= 200 mA: 0, experimental results; chain curve,
`theory with ne = 2,2 x 10‘? cm“3, Te =18000 K, Ta = 7} = 400 K and n., =
`1.21 ><10‘7 cm‘?
`
`the other parameters determined by van der Mullen er
`al (1978) it causes positive or negative changes smaller
`than 3.4% and 30.9% in the population coefficients
`nlf” and GE,”, respectively. The deviations of G E,” are
`strongly dependent on the level number 22. For n S 21,
`i.e. when the ionisation energy 8,, 2 0.961 eV (Vlcek
`1989). all changes of G1,“ are smaller than 3.0%. In
`figure 4, we show the effect caused by neglecting the
`atom—atorn collisions
`in calculating the value of
`Gfgllnj for the excited effective level denoted by n =
`46. which includes all states with the core quantum
`number jc :% and 22”" = 10. For all the other excited
`levels given in this figure we have registered deviations
`smaller than 1.0% even at He = 10” cm“3.
`
`2.2. Positive column of the glow discharge
`
`In recent years numerous studies have been devoted to
`the modelling of a positive column of an argon glow
`discharge in which self-consistent calculations of EEDFS
`and electron excitation and ionisation rates have been
`included‘
`
`Effective analytical methods for solving this problem
`have been developed, in particular by Golubovskii er a1
`(1972, 1976), Smits and Prins (1979a. b) and Ferreira
`and Ricard (1983). Many other authors (see for example
`Morgan and Vriens 1980) based their approach on an
`accurate numerical solution of the Boltzmann equation.
`In order to account for the deviations of the actual EEDF
`
`689
`
`
`
`4p states on the fundamental plasma parameters 713, TC
`and the electric field strength E in an argon glow
`discharge with the discharge current I in the range
`25-400 mA and the filiing pressure p in the range
`O.18—l0 Torr. The discharge tube radius R was 1.2 cm.
`The experimental work of Kagan et al ( 1963b) has
`been Chosen for two reasons. First. a set of the discharge
`parameters measured by these authors is consistent with
`the input parameters in our CR model. Secondly,
`in
`addition to the excited level populations in the 45 states.
`those in higher-lying levels have also been determined.
`in contrast with other studies.
`
`In the same way as Vriens (1973) and Morgan and
`Vriens (1980). who dealt with the applicability of the
`two—electron group model. we use as input parameters
`those sets of quantities which correspond to the pressure
`p =STorr. As Kagan at al (l963b) did not give the
`values measured for T3 under the various conditions
`and as the estimate Ta = 293 K used by Vriens (1973)
`and Morgan and Vriens (1980) is unrealistic. the most
`probable values of Ta evaluated on the basis of com-
`parison with experimental results obtained under similar
`discharge conditions are utilised in our computations
`(Vlcek and Pelikén 1985).
`The chosen sets of experimental discharge par-
`ameters are listed in table 2 together with those esti-
`mations of Ta and :11 which define the intervals of their
`probable values for the corresponding measurements.
`Under the assumption that T3 —- Ti and taking the dis-
`charge tube radius to be R = 1.2 cm, we obtain complete
`sets of input parameters for which the numerical results
`for the excited level populations have been compared
`with measurements.
`
`In calculating the escape factors A,,,, the pressure-
`broadening term is dominant for all excited levels except
`those with n = 12 and 26 (Vléek 1989). All values of
`A1,, are in the range 4.94 X 10'4—l.78 X 10‘: under the
`conditions investigated.
`Figures 5 and 6 show a comparison between the
`calculated populations 21,. of the excited levels charac-
`terised in table 1 and the corresponding experimental
`data obtained at discharge currents of 25. 50. 200 and
`400 mA. In accordance with Kagan er al (l963a). an
`error of a factor of two was expected in the measure-
`ments of n,,.
`As can be seen from these figures. the calculated
`populations are in better agreement with the measure-
`ments for the 45 states than for the 4p levels. for which
`the behaviour of the functions n,,(n) agrees well with
`the corresponding experimental dependences. but the
`calculated values of nn are somewhat overestimated
`compared with the measurements. Better agreement
`between theory and experiment is obtained for two
`resonance 45 states than for two 43 metastables in con-
`trast to modelling of the low-pressure argon positive
`column carried out by Ferreira et al (1985).
`In figure 5 a different effect of the two choices of Z,
`and :2, given in table 2 on the populations n,, is shown
`for discharge currents of 25 and 400 mA. It may be
`explained by the fact that the tail of the EEDF is much
`
`J Vlcek and V Pelikan
`
` .l
`
`l
`
`
`
` $.:......_i,.l_..__J.__l__.1l
`
`.._.,l_.
`
`2 !e“~’I'
`
`Figure 7. The ratio of the EEDF f(s) to the corresponding
`Maxwelllan function fM(£) for the following input
`parameters: full curve A, ne = 2.7 X 10“ cm‘3.
`T9 = 25000 K, 7;” = T?) = 300 K, ng” = 1,51 ><10"7cm‘3;
`broken curve A, Tf’ = Tl” = 400 K, nf’ = 1.21 x
`10” cm'3; curve B, ne = 6.0 X 10“ cm‘3, T5 = 22000 K,
`Ta = 7: = 320 K,n,=1.51 x 10’7crn“3; curve O.
`ne = 2.2 ><10‘2cm'3, Te = 16000 K, Ta = Ti = 400 K,
`I11 =1.21 ><1O‘7cm‘3;curve D, ne = 4.6 ><1O‘2cm’3,
`Te =13000 K. Ta = T, = 500 K, n, = 9.66 ><10‘5cm-3.
`
`from the corresponding Maxwellian form in a simple
`way, Vriens (1973) and Morgan and Vriens (1980) pro-
`posed an approximate treatment utilising two groups of
`Maxwellian electrons. representing the body and the
`tail of the distribution, respectively.
`Recently, a realistic self—consistent discharge model
`for a low-pressure argon positive column based on a
`reliable set of electron cross sections (Ferreira and Lou-
`reiro 1983) has been worked out and tested experi-
`mentally by Ferreira et al (1985). The maintenance
`electric field strength together with the populations in
`the individual 4s states has been calculated from a
`set of coupled equations expressing the steady state
`discharge ionisation balance and rate balance for all 45
`states as a function of the gas pressure, the discharge
`current and the discharge tube radius.
`In our case the numerical results for the excited level
`populations are compared with the experimental values
`of Kagan et al (l963b). who investigated the dependence
`of the populations in the four 45 states and in the nine
`
`540
`
`
`
`Collisional-radiative model for argon discharges, ll
`
`6
`I
`
`§
`
`.3
`
`4°
`A
`I
`‘
`
`A2
`04
`Qa
`‘
`I
`-
`
`V
`
`v
`V
`V
`v
`
`V
`
`V
`
`w
`
`V
`'
`
`v
`
`'
`
`v
`'
`
`L.
`
`10°
`
`E
`
`M_
`
`__.
`
`V‘-:
`
`\
`E
`
`p...
`
`,
`F‘
`
`’)
`I.._
`
`_
`
`j
`--
`
`Q ano
`A
`l
`
`V
`
`V
`
`v
`r
`
`5
`o
`
`‘
`
`7
`
`,
`
`A
`
`3
`
`‘
`
`v
`
`v
`
`3
`0
`
`‘
`
`v
`
`V
`
`_..
`
`..__
`
`“l
`
`_
`
`_
`
`"'7
`_l
`
`*4... I
`2
`L
`
`l
`I’:
`
`l
`B
`
`_ rel
`‘:3
`
`___, ,,
`
`Figure 8. The calculated ratios n,,./n;°', as functions of the effective principal
`quantum number n;,q,, obtained with input parameters ne = 2.7 X 10“ cm‘3.
`Ta = 25000 K, Ta: '7}: 300 K, n-, = 1.61 X 10‘7cm‘3 and H=1cm: A,jc=‘;:
`A, jc = %; O, jc = l but without atom—atom collisions; O, jc = % but without atom-
`atom collisions; and with ne = 4.6 X 1012 cm‘3, Te = 13000 K; Ta = ‘E = 500 K,
`n, = 9.66 ><10‘5cm”9 and Fl=1cm:V,jc = '42; V, jc = %.
`
`more dependent on the value of the degree of ionisation
`nc/"n, in the case of I = 25 mA than for I = 400 mA. At
`25 mA. a decrease of the ground state atom population
`rzl connected with an increase in the atom temperature
`Ta leads to essentially enlarged values of the rate coef-
`ficients for the excitation by electrons from the ground
`state. which are not compensated for by the lowered
`value of n] in the term Gfflrzl determining the popu-
`lations n,_ of the levels considered under these con-
`ditions. Both effects are approximately balanced when
`I = 400 mA.
`-
`
`In the papers of Kagan er al (196321, b) there are no
`data concerning the accuracy in probe measurements of
`TC and He»
`In figure 6 the deviations of the populations ran
`
`caused by small changes in the values of Te and in are
`shown for a discharge current of 50 mA. for which the
`largest discrepancies between theory and experiment
`exist. As can be seen. a decrease of the electron number
`density by 10% leads to a drop in the populations rzn by
`l8~29%. whereas a decrease of
`the electron tem-
`perature only to the value of 20700 K, representing a
`reduction in the experimental value (Kagan er al 1963b)
`of 5.9%, results in the lowering of the populations
`investigated by 36-43%. This reduced value of T3 was
`obtained (Morgan and Vriens 1980) with the help of
`the standard Boltzmann calculations at the discharge
`parameters determined experimentally by Kagan er al
`(l963b) for I = 50 mA.
`In figure 7 we give the numerical results for the ratio
`
`641
`
`
`
`J Vlcek and V Pelikén
`
`to‘-3|‘ "
`r
`
`“
`
`A
`
`B
`
`fig
`
`'
`
`”
`
`‘
`
`i
`4
`
`*1
`
`
`
`
`
`“ii”.CIW71lam") 5 tr
`
` nu:
`
`(broken curves)
`ri‘,?’
`Figure 9. The calculated values of
`and Gifln, (full curves), together with the measured values
`(symbols with error bars) of n,, (Kagan er al 1963b), for
`selected excited levels as functions of the discharge
`current [with the first choices of Ta = T, and n., given in
`table 2: curves A and O, n = 2; curves B and A, n = 3;
`curves C and V, n= 7; curves D, Sp (1; = 3); curves E, 5d
`(1; = i); curves F, 7s (jc = 3); curves G, npqn =10(j, = i).
`
`f (s)/fM(e) to compare the deviations of the realistic
`EEDF f (.9) from the Corresponding Maxwellian function
`fM(e) at several sets of input parameters listed in table
`2 and used in calculations presented in figures 5 and 6.
`The calculated ratios 21,, _/n,5, as functions of the effec-
`tive principal quantum number n;‘q,, are shown in figure
`8 for plasmas with discharge currents of 25 and 400 mA
`and they illustrate a strong tendency to a decrease in
`the overpopulation in higher excited levels compared
`with the Saha values. In addition, the effect of the atom-
`atom inelastic collisions on the excited level populations
`is presented for the case of 1 = 25 mA.
`Figure 9 shows the calculated values of nE,°l and
`GSPI1, and the measured values of n,, for several excited
`levels as functions of the discharge current I with the
`first choices of atom temperature T,, and ground state
`atom population n1 given in table 2. In contrast to
`the situation in the low-pressure, hollow-cathode are
`discussed in §2.l (see figure 4), the upward ionisation
`flow of electrons from the ground state of the argon
`atom is a dominant populating mechanism for all but
`
`642
`
`the highest-lying excited levels over the whole range of
`the discharge conditions investigated. Compared with
`Gifinl , the terms iii?’ are much more dependent on the
`value of 1 corresponding to the set of plasma parameters
`Tc, Ta. 215 and n1 listed in table 2. When the discharge
`current I grows, the term G,(,”n1 is given by the balance
`between the effects caused by a drop in TE and n1.
`lowering the terms G§—,"n1, and those caused by a rise
`in ne. enlarging them, whereas both changes in Tc and
`ne lead to an increase in the value of the population
`coefficient n£,°:’.
`The measured values of Kagan eta! (1963b) available
`for the selected excited levels over the whole range of
`the discharge current are also shown in figure 9 for
`comparison.
`'
`
`3. Conclusions
`
`The extensive co1lisional—radiative model for an argon
`atom plasma has been applied to a low-pressure, hollow-
`cathode arc and to the positive column of a low-pressure
`glow discharge in order to clarify the mechanisms by
`which the excited levels in these discharges are popu-
`lated, these having been investigated experimentally by
`van der Mullen er al (1978, 1980). van der Sijde er al
`(19843, b) and Kagan et al (1963b), respectively.
`Our computations are carried out for various sets of
`the input parameters Te, Ta, Ti, 22,. n1. R. Am and
`Am which were measured directly in the experiments
`mentioned above or determined from the available
`experimental data.
`From our results the following main conclusions can
`be drawn.
`
`(i) Our model provides a qualitatively good picture
`of the processes determining the populating mechanism
`both in the hollow-cathode are considered. where the
`electron energy distribution function is Maxwellian and
`the radiation trapping is negligibly small, and in the
`glow discharge investigated, where the Boltzmann
`equation for the electron distribution function must be
`solved and the resonance radiation is strongly reab-
`sorbed.
`
`the
`In the low-pressure, hollow-cathode arc.
`(ii)
`predicted values of the populations in all the groups
`considered, i.e. 4p, 5p, 5d and 7s. agree well with the
`measurements of van der Mullen er al (1978). but all
`the groups considered. i.e. 4p, 5d and 6d, are somewhat
`underpopulated compared with the experimental values
`of van der Mullen et al (1980). In the positive column
`of a low-pressure glow discharge. reasonable agreement
`is found between the calculated and measured popu-
`lations in all individual 45 states, but compared with the
`experiment the values obtained for all 4p states are
`systematically overestimated, presumably owing to the
`lack of accuracy in the determination of the electron
`kinetic temperature used as an input parameter in our
`computations.
`(iii) The analy