A&A 405, 397-403 (2003)
DOI: 10.1051/0004-6361:20030612
V. Milosavljevic1,2 - S. Djenize1,2,3
1 - Faculty of Physics, University of Belgrade, PO Box 368,
Belgrade, Serbia
2 -
Isaac Newton Institute of Chile, Yugoslavia Branch,
Belgrade, Serbia
3 -
Hungarian Academy of Sciences, Budapest, Hungary
Received 30 January 2003 / Accepted 15 April 2003
Abstract
On the basis of five accurately recorded neutral argon
(Ar I) line shapes (in the 4s-5p transition) we have
recovered the basic plasma parameters i.e. electron temperature (T) and electron density (N) using our new line deconvolution
procedure in the case of three different plasmas created in a
linear, low-pressure, pulsed arc discharge. The mentioned plasma
parameters have also been measured using independent experimental
diagnostic techniques. An excellent agreement has been found among
the two sets of obtained parameters. This recommends our
deconvolution procedure for plasma diagnostic purposes, especially
in astrophysics where direct measurements of the main plasma
parameters (T and N) are not possible. On the basis of the
observed asymmetry of the Stark broadened line profile we have
obtained not only its ion broadening parameter (A) which is
caused by the influence of the ion microfield on the line
broadening mechanism but also the influence of the ion-dynamical
effect (D) on the line shape. The separate electron ()
and ion (
)
contributions to the total Stark width,
which have not been measured so far, have also been obtained and
represent the first experimental data in this field. We have found
a stronger influence of the ion contribution to the Ar I line profiles than the existing theoretical approximation
provides. This is of importance for astrophysical plasma modeling
and diagnostics.
Key words: plasmas - line: profiles - atomic data
The presence of the neutral argon (Ar I) spectral lines has
been discovered in various cosmic light sources in the last few
years. Recently, Weaver et al. (2002) have referred to the presence of
Ar I lines in the spectra of long-period comets. The
Ar I absorption lines have been detected in the spectra of
the quasar Q0347-3819 (Levshakov et al. 2002) and PG 1259 + 593
(Richter et al. 2001). In the studies by Friedman et al. (2000), Jenkins et al. (2000),
Kruk et al. (2002) and Lehner et al. (2002) the absorptions in the Ar I lines were presented. Argon is detected in the spectrum of the
damped Ly system of IZw 18 (Levshakov et al. 2001). Mallouris et al. (2001)
refer to the presence of Ar I lines in the spectrum of the
Wolf-Rayet binary SK 108. Thus, the Ar I spectral line
shapes represent important sources of information about the
physical conditions in the place of birth of the radiation,
especially since the launch of the Hubble space telescope.
In our work we applied our line deconvolution procedure
(Milosavljevic & Poparic 2001) to five precisely recorded Ar I line
profiles. The basic plasma parameters, i.e. electron temperature ()
and electron density (
)
have been
obtained, using our line deconvolution procedure, in the case of
three different plasmas created in a linear, low-pressure, pulsed
arc discharge in helium-argon and hydrogen-argon mixtures. To
our knowledge, our results of the T and N values have been the
first published data obtained directly from the line profile,
using a deconvolution procedure. Plasma parameters have also been
measured (
and
)
using independent,
well-known, experimental diagnostical techniques. Excellent
agreement has been found within the two sets of the obtained
parameters (
and
;
and
and
). This recommends our deconvolution procedure for
plasma diagnostical purposes, especially in astrophysics where
direct measurements of the plasma parameters (T and N) are not
possible. The method is applicable for optically thin plasmas. For
the optically thick plasmas the situation is more complicated
because of very weak asymmetry of the spectral line profile caused
by self-absorption.
In plasmas with electron densities higher than 10 m-3, where the Stark effect begins to play an
important role in the Ar I spectral line broadening, the
Stark broadening characteristics can also be used for plasma
diagnostical purposes. However, the existing experimental
Ar I Stark width data
(Lesage & Fuhr 1999; NIST 2002; Konjevic et al. 2002, and references
therein) present only the sum of the electron
(
)
and ion (
)
contributions to the total
Stark width (
)
without any possibility for their
separation. The
and
Ar I values
presented here are the first data, in this field, separated from
the measured total Stark width by using the line deconvolution
procedure described by Milosavljevic & Poparic (2001) which have already been
applied for some He I, Ne I, Ar I and Kr I lines (Milosavljevic & Djenize 2003,2002b,c,a). A
significant number of experimental studies have been dedicated to
the Ar I total Stark FWHM (full-width at half intensity
maximum,
)
investigations of the 4s-5p transition.
However, only one study (Griem 1974) contains theoretical
,
and
values in the 4s-5p transition.
In this paper we are presenting the measured Stark broadening
parameters of the 415.86 nm, 416.42 nm, 419.83 nm, 420.07 nm
and 426.63 nm Ar I spectral lines (in the 4s-5p transition) at about 16 000 K electron temperature and at about
m
electron density. The used
T values are typical for many cosmic light sources. On the basis
of the observed Ar I line profile asymmetry, the
characteristics of the ion contribution to the total Stark FWHM
(
), expressed as a function of the ion contribution
parameter (A) and ion-dynamical effect (D) have also been
obtained. As an optically thin plasma source we have used a
linear, low-pressure, pulsed arc operated in three various
discharge conditions. Our
,
,
and A values have been compared to all available theoretical and
experimental Stark broadening parameters.
The total line Stark FWHM ()
is given as
![]() ![]() |
For the purpose of the deconvolution iteration process we need to
know the value of K (Eq. (6)) as a function of
for every group of parameters (
,
,
,
,
R, A).
is
defined in Eq. (2.3) of Milosavljevic & Poparic (2001). The used numerical
procedure for solution of Eq. (6) is described in
earlier publications (Milosavljevic & Poparic 2001; Milosavljevic 2001; Milosavljevic & Djenize 2002a).
Working gases | Exp. | C | U | H | ![]() |
P |
![]() |
![]() |
![]() |
![]() |
72% Ar + 28% He | a | 14 | 1.5 | 7.2 | 5 | 133 | ![]() |
![]() |
![]() |
![]() |
97% Ar + 3% H2 | b | 14 | 1.5 | 7.2 | 5 | 67 | ![]() |
![]() |
![]() |
![]() |
97% Ar + 3% H2 | c | 14 | 1.5 | 7.2 | 5 | 133 | ![]() |
![]() |
![]() |
![]() |
From Eqs. (1-6) it is possible to obtain
the plasma parameters (N and T) and the line broadening
characteristics (,
,
,
A and D). One can see that the ion contribution, expressed in terms of
the A and D parameters directly determine the ion width (
)
component in the total Stark width
(Eqs. (1) and (2)).
This sophisticated deconvolution method, which allows direct
determination of all six parameters by fitting the theoretical
K-profile (6) to experimental data, requires a
sufficient number of experimental points per line, and small
statistical errors. The upper limits of numerical conditionality
of this method are a minimum of twenty experimental points per
line (within the range
),
and maximal statistical indeterminacy in intensity is 5% at every
experimental point. Poor experimental measurements weaken the
conditionality of the system of equations, and lead to
non-applicability of this method. This has been concluded by
testing the sensitivity of the algorithm by generating random
statistical noise with Gaussian distribution in every point
convolved by theoretical profiles. The fitting procedure with the
K-convolution integral has also been tested using another set of
experimental data (see Milosavljevic & Djenize 2002b,c,a). The K convolution integral is used for the analysis of our new data for
many spectral lines of neutral rare gases. By comparing the
different spectral lines obtained under the same plasma
conditions, we tested the physical stability of the deconvolution
procedure. The obtained parameters, which are tied to plasma
conditions, such as T and N, are independent from the analyzed
lines. Our calculated values of temperature from each spectral
line and values obtained by Boltzmann and Saha equations are in
very good agreement, within
7%. The electron density
calculated from each spectral line shows even better agreement
with the values measured by interferometry, the agreement being
within
5%.
Taking into account the uncertainties of the line profile
measurements and above mentioned, we estimate errors 12% for
the
and
,
15% for the A parameter and
20%
for D.
The modified version of the linear low pressure pulsed arc
(Djenize et al. 2002,1991,1998; Milosavljevic et al. 2000,2001) has been used as a plasma
source. Pulsed discharge was performed in a quartz discharge tube.
The working gases were helium-argon (28% He + 72% Ar) and
hydrogen-argon (3% H + 97% Ar) mixtures. The used
tube geometry and corresponding discharge conditions are presented
in Table 1.
The spectroscopic observation of spectral lines has been made end-on along the axis of the discharge tube.
The line profiles were recorded by a step-by-step technique using
a photomultiplier (EMI 9789 QB and EMI 9659B) and a grating
spectrograph (Zeiss PGS-2, reciprocal linear dispersion 0.73 nm mm-1in the first order) system. The instrumental FWHM of 8 pm was
obtained by using narrow spectral lines emitted by the hollow
cathode discharge. The spectrograph exit slit (10 m) with the
calibrated photomultipliers was micrometrically traversed along
the spectral plane in small wavelength steps (7.3 pm). The
averaged photomultiplier signal (five shots in each position) was
digitized using an oscilloscope, interfaced to a computer.
Plasma reproducibility was monitored by the Ar I and
Ar II line radiation and, also, by the discharge current
using a Rogowski coil signal (it was found to be within 5%).
The used deconvolution procedure in its details is described in
Milosavljevic & Poparic (2001) and Milosavljevic (2001). It includes a new advanced
numerical procedure for deconvolution of theoretical asymmetric
convolution integral of a Gaussian and a plasma broadened spectral
line profile jA,R()
for spectral lines. This method
gives complete information on the plasma parameters from a single
recorded spectral line. The method determines all broadening
(
,
,
,
A and D) and plasma
parameters (N and T) self-consistently and directly from the
shape of spectral lines without any assumptions or prior
knowledge, making it useful in astrophysics. All one needs to know
is the instrumental width of the spectrometer. The measured
profiles are the results of convolution with the Lorentzian Stark
and Gaussian profiles caused by Doppler and instrumental
broadening (Griem 1974). Van der Waals and resonance broadenings
(Griem 1974) were estimated to be smaller by more than an order
of magnitude in comparison to Stark, Doppler and instrumental
broadenings. The deconvolution procedure was computed using the
least Chi-square function (Milosavljevic & Poparic 2001).
The absence of self-absorption was checked using the method presented in Djenize & Bukvic (2001).
The plasma parameters were determined independently using standard
diagnostics methods. Thus, the electron temperature was determined
from the ratios of the relative line intensities of seven
Ar I spectral lines (415.859 nm, 416.418 nm, 419.103 nm,
419.832 nm, 420.067 nm, 425.936 nm and 426.627 nm) to the five
Ar II spectral lines (335.093 nm, 420.197 nm,
426.653 nm, 487.986 nm, 488.903 nm) with an estimated error of 11%, assuming the existence of LTE (Griem 1974). The
necessary atomic data have been taken from NIST (2002). The
electron density decay was measured using a well-known single
wavelength He-Ne laser interferometer technique for the 632.8 nm
transition with an estimated error of
7%. The electron
densities (
)
and temperatures (
),
obtained at the moment when the line profiles were analyzed, are
presented in Table 1 together with the
and
values obtained using the deconvolution procedure.
![]() |
Figure 1:
Electron temperature (T) and density (N) decays.
Full lines represent measured data using independent experimental
techniques. Dashed lines represent plasma parameters obtained
using our line deconvolution procedure in various plasmas (see a,
b and c in Table 1). Error bars represent estimated
accuracies of the measurements (![]() ![]() ![]() |
The measured
and
decays are presented
in Fig. 1 together with the averaged
and
values obtained using the line profile deconvolution
procedure for five Ar I lines. One can conclude that the
agreement among
and
values is excellent
(within 3% on average in the three plasmas investigated). This
fact confirms the homogeneity of the investigated plasmas in the
linear part of our light source (see Fig. 1 in
Djenize et al. 1998). In
the case of the electron density, the situation is similar. The
agreement among the two sets of the electron density decays
(
and
)
is within the experimental
accuracy of
7% and uncertainties (
12%) of the results
obtained with the deconvolution procedure.
The plasma broadening parameters (
,
,
,
,
)
obtained using our deconvolution procedure of the recorded
line profiles at measured
and
values
are presented in Table 2 together with those of other
authors. Theoretical (Griem 1974) predictions (index G) of the
and A are also given. For the normalization of the
values to our electron density the well-known
numerical factor (Griem 1974) was used.
Multiplet | ![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Ref. |
![]() |
![]() |
![]() |
419.83 | 15.6 | 6.7 | 130 | 112 | 18 | 0.143 | 1.347 | Tw | 223 | 0.113 |
16.0 | 7.0 | 137 | 117 | 20 | 0.145 | 1.656 | Tw | 235 | 0.114 | ||
16.2 | 7.1 | 139 | 119 | 20 | 0.145 | 1.649 | Tw | 239 | 0.114 | ||
11.1 | 0.4 | 3.3 | KM | ||||||||
12.4 | 7.3 | 120 | B | ||||||||
14.0 | 1.0 | 23 | G | ||||||||
![]() |
415.86 | 15.6 | 6.7 | 132 | 114 | 18 | 0.143 | 1.332 | Tw | ||
16.0 | 7.0 | 139 | 119 | 20 | 0.146 | 1.638 | Tw | ||||
16.2 | 7.1 | 143 | 122 | 21 | 0.147 | 1.624 | Tw | ||||
11.9 | 6.2 | 121 | 0.127 | J | |||||||
11.9 | 6.2 | 112 | 0.125 | HW | |||||||
13.5 | 10.0 | 246 | A | ||||||||
12.7 | 9.4 | 173 | B | ||||||||
11.4 | 4.6 | 110 | G | ||||||||
14.0 | 1.0 | 18.3 | P | ||||||||
12.5 | 9.2 | 190 | MM | ||||||||
11.9 | 6.2 | 123 | 0.127 | JP | |||||||
13.4 | 4.5 | 107 | M | ||||||||
![]() |
416.42 | 15.6 | 6.7 | 126 | 109 | 17 | 0.140 | 1.357 | Tw | ||
16.0 | 7.0 | 134 | 115 | 19 | 0.142 | 1.664 | Tw | ||||
16.2 | 7.1 | 134 | 115 | 19 | 0.142 | 1.663 | Tw | ||||
11.9 | 6.2 | 97 | 0.196 | HW | |||||||
13.8 | 14.5 | 332 | Ch | ||||||||
11.1 | 0.4 | 3.7 | KM | ||||||||
12.7 | 9.4 | 171 | B | ||||||||
11.4 | 4.6 | 100 | G | ||||||||
14.0 | 1.0 | 18.1 | P | ||||||||
10.6 | 2.6 | 59.6 | MM | ||||||||
11.5 | 5.1 | 96 | M | ||||||||
![]() |
426.63 | 15.6 | 6.7 | 127 | 110 | 27 | 0.132 | 1.380 | Tw | ||
16.0 | 7.0 | 130 | 113 | 27 | 0.132 | 1.712 | Tw | ||||
16.2 | 7.1 | 136 | 118 | 28 | 0.134 | 1.689 | Tw | ||||
11.9 | 6.2 | 126 | 0.076 | HW | |||||||
11.4 | 4.6 | 110 | G | ||||||||
14.0 | 1.0 | 17.9 | P | ||||||||
12.6 | 9.9 | 240 | M | ||||||||
![]() |
420.07 | 15.6 | 6.7 | 138 | 119 | 19 | 0.140 | 1.325 | Tw | 153 | 0.085 |
16.0 | 7.0 | 143 | 123 | 20 | 0.142 | 1.636 | Tw | 162 | 0.085 | ||
16.2 | 7.1 | 147 | 126 | 21 | 0.142 | 1.624 | Tw | 165 | 0.085 | ||
13.8 | 14.5 | 324 | Ch | ||||||||
11.1 | 0.4 | 2.8 | KM | ||||||||
12.7 | 9.4 | 155 | B | ||||||||
11.4 | 4.6 | 100 | G | ||||||||
12.6 | 9.2 | 180 | GR | ||||||||
13.4 | 4.5 | 94 | M |
![]() |
Figure 2:
Ratios of the experimental total Stark FWHM (
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
In order to make the comparison among the measured (
)
and calculated (
)
total
(electron + ion) width values easier, the dependence of the ratio
/
on the electron
temperature is presented graphically in Figs. 2 and 3. The
(Griem 1974) values are
calculated using Eq. (226) from Griem (1974).
Our broadening parameters (
)
represent the
first measured values at electron temperatures higher than 14 000 K.
The comparison of our
,
and
values with uniquely theoretical data
(Griem 1974) was possible for only two transitions. Thus, in the
case of the 419.832 nm Ar I spectral line Griem's (1974)
theoretical
values are about twice as high as
the experimental
values (see
Fig. 2). Also, the
values are about
twice as high as our
values. We think that
the discrepancy between measured and theoretical values indicates
that the atomic data used for calculations (Griem 1974) were
unreliable. Our ion broadening parameters (
)
overvalue the theoretical (
) values by about 26% and
are multiplied 1.3-1.6 times with the ion-dynamical effect,
depending on the discharge conditions. In the case of the
420.068 nm Ar I line all experimental
values lie below Griem's (1974) values by about 20% on
average (see Fig. 3). Our
values
are about 20% smaller than the theoretical
data. Our ion broadening parameters (
)
overvalue for
about 65% the theoretical (
)
values and are multiplied
1.3-1.6 times with the ion-dynamical effect, depending on the
discharge conditions. Generally, the electron contribution to the
total Stark width found experimentally is about 85% (on average)
at about 16 000 K electron temperature.
![]() |
Figure 3:
Ratios of the experimental total Stark FWHM (
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
It turns out that
values obtained by
Jones et al. (1986,1987) and Hahn & Woltz (1990) for the 415.86 nm line
agree well with ours. In the case of the 416.42 nm line,
Hahn & Woltz (1990) have found
higher than ours (see
Table 2). Direct comparison of our
values with other experimental data for the 415.86 nm, 416.42 nm
and 426.63 nm lines is impossible, because of the various plasma
conditions in the experiments. But, taking into account the
approximative normalization factor (
)
we have found a
tolerable scatter for the mentioned transitions among
experimental values (within 15%-38%). The only exceptions
are very low
values from Kusz & Mazur (1996) which are about
twice as small as ours and the
values of other
authors.
We have found clear influence of the quasi-static ion and ion-dynamical effects on the investigated spectral line shapes. They play a much more important role than the semiclassical theory provides. The observed ion-dynamical effect, at our plasma conditions, multiplies the quasi-static ion contribution by up to a factor of 1.6. This is of importance for the use of these Ar I lines for astrophysical plasma modeling or for diagnostics.
Acknowledgements
This work is a part of the project "Determination of the atomic parameters on the basis of the spectral line profiles'' supported by the Ministry of Science, Technologies and Development of the Republic of Serbia. S. Djenize is grateful to the Foundation "Arany János Közalapitvány'' Budapest, Hungary.