Free Access
Issue
A&A
Volume 508, Number 1, December II 2009
Page(s) 491 - 500
Section Atomic, molecular, and nuclear data
DOI https://doi.org/10.1051/0004-6361/200912046
Published online 08 October 2009

A&A 508, 491-500 (2009)

Line broadening in the Si I, Si II, Si III, and Si IV spectra in the helium plasma

S. Bukvic - S. Djenize - A. Sreckovic

Faculty of Physics, University of Belgrade, Studentski trg 14, 11000 Belgrade, Serbia

Received 12 March 2009 / Accepted 10 September 2009

Abstract
Context. The neutral and ionized silicon spectral line shapes have been investigated in the laboratory helium plasma at electron densities ranging between $3.7\times 10^{22}$ m-3 and $1.1\times 10^{23}$ m-3 and electron temperatures between 12 500 K and 19 000 K, both interesting for astrophysics.
Aims. The aim of this work is to present experimental Stark FWHM (full-width at half of the maximum line intensity, W) for number of spectral lines from neutral (Si I), singly (Si II), doubly (Si III), and triply (Si IV) ionized silicon spectra emitted by the pulsed helium discharge, which is optically thin at the wavelengths of the investigated ionic silicon lines. A specific method for estimating self-absorbtion is presented in detail. For investigated Si I spectral lines, applying the proposed method, an optical depth of less than 0.38 is found. Appropriate corrections of the Si I Stark widths were made. The Stark widths of different ionic species, presented in this paper, are measured for the first time in the essentially same laboratory plasma.
Methods. The silicon atoms were evaporated from the walls of the specially designed pyrex discharge tube in the pulsed helium discharge at a pressure of 665 Pa in a flowing regime. The Si I, Si II, Si III, and Si IV spectral line profiles were recorded using the McPherson model 209 spectrograph and the Andor ICCD camera as the detection system.
Results. The Stark FWHMs of 13 Si I, 15 Si II, 28 Si III, and 9 Si IV spectral lines were measured in the wavelength interval between 206 nm and 640 nm. Five Si I, four Si II, eleven Si III, and one Si IV W values from the above set not had measured or calculated. Our W values are compared with the existing theoretical and experimental data.
Conclusions. At the mentioned plasma parameters tolerable agreement was found (within the accuracy of the experiment and uncertainties of the theoretical approaches used) between measured and calculated Stark FWHM values. We recommend the Stark FWHMs of the intense 254.182 nm, 308.624 nm, and 309.342 nm Si III, and 314.956 nm and also 316.571 nm Si IV lines for the plasma diagnostic purposes.

Key words: plasmas - line: profiles - atomic data - atomic processes

1 Introduction

The silicon (Si) atoms and ions are present in many kinds of cosmic plasmas. The neutral (Si I), singly (Si II), doubly (Si III), and triply (Si IV) ionized silicon spectral lines are discovered in the emission and absorption spectra emitted by various cosmic light sources. We mention a few recent articles emphasizing the role of the prominent Si I-IV spectral lines in various fields of the astrophysical plasma diagnostics.

Statistical equilibrium of neutral and ionized silicon in the solar photosphere has been investigated by Shi et al. (2008). The Si I line was discovered in the emission spectra of a K giants (Sundqvist et al. 2008). The equivalent widths of Si II 400 nm lines have been used for the determining diversity of supernovae Ia (Arsenijevic et al. 2008). In the 500-700 nm wavelength interval, the SN Ia - defining Si II feature occurs (Bongard et al. 2008). Prominent Si II emission lines have been observed in the first helium nova V445 Puppis (Iijima & Nakanishi 2008). Strong Si IV ion absorption has been identified in seven GRB host galaxies (z=2-4) (Fox et al. 2008). The presence of Si III and Si IV is identified in the interstellar medium seen towards the nearby (170 pc) star HD 102065 (Nehme et al. 2008). Kawka et al. (2008) refer to the detection of Si IV in the spectra of the hot white dwarfs in the post-common-envelope binaries, Feige 24, EUVE J00720-317, and EUVE J2013+400, and to detection of Si III and Si IV in the spectra of BPM 6502. Some Si III lines have been observed in the spectrum of the rapidly rotating classical T Tauri star RY Tau (Gomez de Castro & Verdugo 2007) and the wind kinematics were derived from the Si III line profiles. The Si III and Si IV ions, as absorbers, have been discovered in the low-z(z<0.4) intergalactic medium (Danforth & Shull 2008). The emergence of a high-velocity broad-line outflow in a luminous quasar J105400.40+034801.2 is evident in ultraviolet Si IV absorption lines (Hamann et al. 2008). The Si III triplet at 455 nm are used by Catanzaro & Leone (2008) to determine the spectroscopic variability of the effective temperature ( $T_{\rm eff}$) in the hot pulsating star Beta Cephei. These lines, relatively insensitive to the variation in $T_{\rm eff}$ through a pulsation cycle, show small equivalent width variations (Catanzaro et al. 2008).

Thesee references clearly confirm the recent astrophysical interest for the neutral and ionized silicon spectral line characteristics. In many models of astrophysical plasmas (Lesage 1994; Dimitrijevic et al. 2007) the Stark width, beside the Doppler caused broadening, is essential. A summary of the papers dedicated to the Si I-IV Stark width calculations is available in NIST (2009), while the critical reviews of the measured values are given in Lesage (2009) and Konjevic et al. (2002), and references therein.

The aim of this work is to present experimental Stark widths for a number of spectral lines from neutral, singly, doubly, and triply ionized silicon spectra emitted by the pulsed helium discharge, which is optically thin at the investigated ionic silicon wavelengths. After applying procedure proposed in the Appendix, we have found a weak self-absorption in Si I spectra with optical depth <0.38. Results for the Si I lines are corrected in the amount of 2-22%. The specific design of the discharge tube provides the Si I, Si II, Si III, and Si IV spectra within the single pulse at the different stages of the discharge. Therefore, the Stark widths of different ionic species presented in this paper are for the first time measured in the essentially same laboratory plasma. The Stark FWHM (full-width at half of the maximum line intensity, W) values are obtained in the wavelength interval (206-640) nm at the electron densities between $3.7\times 10^{22}$ m-3 and $1.1\times 10^{23}$ m-3 and at the electron temperatures (T) ranging from 12 500 K to 19 000 K, both of which are interesting in astrophysics.

2 Experiment

A modified version of the linear low-pressure pulsed arc (Djenize et al. 1992, 2002) has already been used as a plasma source (see Fig. 1 in Djenize 2009). A pulsed discharge was produced in a pyrex discharge tube of 5 mm inner diameter with a plasma length of 12 cm. The silicon atoms (O, B and Na atoms, also) are evaporated from the axial part of the discharge tube. As a working gas helium (90% He + 7% N + 3%O) was used at 665 Pa flowing pressure. The helium metastables favor the ionization of the silicon atoms due to the Penning effect (Kruithof & Penning 1937) ensuring a high concentration of the Si II ions.

The discharge is created using the capacitor of 14 $\mu $F charged up to 44 J bank energy. Plasma reproducibility was monitored through the radiation originating from the chosen Si I, Si II, and Si III lines and discharge current using a Rogowski coil signal (Rogowski & Steinhaus 1912). It was found that scatter of the peak value of the discharge current is within ${\pm }4 \%$.

For the first time, we applied the McPherson model 209 spectrograph (1.33 m focal length) equipped with 2400 grooves/mm holographic grating, resulting in reciprocal linear dispersion of 0.28 nm/mm in the first order. As a detection system the Andor DH740-18F-03 iStar intensified CCD camera is employed. The system was carefully calibrated using a set of pen-light (Ne, Ar and Hg) sources produced by the LOT-Oriel. We found that instrumental profile of the spectrograph itself in the first order, conventionally measured (with 9789 QB EMI photomultiplier), corresponds to the Gaussian function of ${\it FWHM} = 3.0$ pm, while the overall profile (spectrograph + ICCD camera) can be approximated by the Voigt function with FWHM of 8.5 pm at 614.3 nm. Standard light guide and focusing optics, employed at some stages of measurement, have no influence on the shape and width of the instrumental profile.

The spectroscopic observations were made end-on along the axis of the discharge tube. The camera was triggered at the specified moment with exposure time adapted to our discharge conditions (0.5 $\mu $s). To reduce thermal noise the ICCD detector is always kept on -25 $^{\circ}$C. Some of the recorded Si I, Si II, Si III, and Si IV spectral line profiles are presented in Figs. 1-4. It should be mentioned that, at our experimental conditions, the Si II 505.598/504.103 lines intensity ratio is more than twice higher than the value which follows from NIST (see Fig. 2c), but agrees with one predicted by transition probabilities and statistical weights related to the mentioned transitions (see Eq. (2) in Djenize & Bukvic 2001, and data in NIST).

The helium plasma parameters were determined using standard diagnostic methods; the electron density (N) was estimated relying on the known Stark FWHM of the He II P$_{\alpha }$ (468.6 nm) spectral line. It should be mentioned that the intensity of the He II P$_{\alpha }$ drops down quickly, and after the 35th $\mu $s signal to noise ratio is not sufficient for accurate processing.

The electron temperature was obtained using the relative line intensity ratio method (Saha equation) (Griem 1964) between Si III (308.62, 309.34, and 309.68 nm) and Si IV (314.96 nm and 316.57 nm), and also between O II (397.33 nm) and O III (396.16 nm) spectral lines with an estimated error of ${\pm } 14\%$ assuming the local thermodynamic equilibrium (LTE) (Rompe & Stenbeeck 1967), and taking the lowering of the ionization energy into account (Drawin & Felenbok 1965). The necessary atomic data are taken from NIST (2009).

Temporal evolutions of the N and T are presented in Fig. 5. Recovery of the electron temperature and electron density at the 30th $\mu $s is due to second half period of the discharge current, caused by subcritically damped electrical circuit (see Fig. 5). Unexpectedly high value of T at the 30th $\mu $s could be attributed, to the some extent, to hot Penning electrons produced in collisions of two helium $\rm 2s~^{3}S_{1}$ metastables with 18.92 eV excitation energy (Silfvast 1971).

The Si I and Si II line profiles are investigated at lower, while the Si III and Si IV line profiles at higher electron temperatures. The Stark FWHMs are obtained at various electron densities ranging between $3.7\times 10^{22}$ m-3 and $1.1\times 10^{23}$ m-3. Their values, normalized at an $N=1.0\times 10^{23}$ m-3, are presented in Tables 2-5.

The spectral purity of the investigated silicon lines are checked using a list of the identified wavelengths given by NIST (2009).

\begin{figure}
\par\includegraphics[width=8cm,clip]{12046Fg1.eps}
\end{figure} Figure 1:

The recorded resonance Si I line profiles (see Table 2). The numbers in brackets denote relative line intensities tabulated by NIST (2009). The resonance B I line profiles are also presented.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8cm,clip]{12046F2.eps}
\end{figure} Figure 2:

Panels a-c) The recorded profiles of prominent Si II spectral lines present in many astrophysical spectra. The numbers in brackets denote relative line intensities tabulated by NIST.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8cm,clip]{12046F3.eps}
\end{figure} Figure 3:

Panels a) and b) The recorded profiles of prominent Si III spectral lines present in many astrophysical spectra. The numbers in brackets denote the relative line intensities tabulated by NIST.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8cm,clip]{12046F4.eps}
\end{figure} Figure 4:

The recorded profiles of the 314.956 nm and 316.571 nm Si IV spectral lines. The numbers in brackets denote relative line intensities tabulated by NIST.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8cm,clip]{12046Fg5.eps}
\end{figure} Figure 5:

Temporal evolutions of the electron density (N in 1023 m-3) and electron temperature (T in 104 K). The error bar represents ${\pm } 14\%$ estimated uncertainties. The gray line is the discharge current (in arbitrary units) captured by the Rogowski coil.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8cm,clip]{12046Fg6.eps}
\end{figure} Figure 6:

The 254.182 nm Si III line profile recorded at $N = 5.2 \times 10^{22}$ m-3 and T= 19 000 K. Circles represent measured values. The solid line represents the corresponding Voigt profile using 4.4 pm for the Gaussian parameter (Doppler + instrumental convolution) in the fitting procedure. The total FWHM ($W_{\rm t}$) (Lorentzian + Gaussian convolution) is 12.02 pm with the Stark component of $W_{\rm m}=3.4 $ pm (see Table 4). The baseline estimation is made using the procedure described by Bukvic & Spasojevic (2005) and Bukvic et al. (2008).

Open with DEXTER

3 Line profile deconvolution procedure

The measured Si II, Si III, and Si IV line profiles were of the Voigt type due to a convolution of the Lorentzian Stark, Gaussian profile caused by Doppler, and Voigt instrumental (spectrograph + ICCD) broadening. For the electron density and temperature in the presented experiment, the Lorentzian fraction was dominant. Van der Waals and resonance broadenings (Griem 1974) were estimated to be smaller by more than one order of magnitude than Stark, Doppler, and instrumental broadenings.

For estimation of spectral line parameters a deconvolution procedure (Davies & Vaughan 1963) is applied based on the least-squares algorithm. Estimation of the spectrum baseline is based on the procedure presented by Bukvic & Spasojevic (2005), and Bukvic et al. (2008). As an example a Voigt fit of 254.182 nm Si III spectral line is presented in Fig. 6. For this, a very-well defined line profile, error of the Lorentz component is $\sim$10%. In general, for ionic lines the error of the Lorentz component is between 10% and 20%, while it is slightly higher for atomic lines, up to 30%. Uncertainty of the Lorentz component for neutral silicon lines is additionally increased for a few percent due to corrections for self-absorption (see Tables 2-5).

\begin{figure}
\par\includegraphics[width=8cm,clip]{12046F7.eps}
\end{figure} Figure 7:

Panel a): the He II P$_{\alpha }$ line recorded in single pass (black) and in double pass (gray). Both profiles are taken at the 10th $\mu $s and same detection conditions. The light is captured only from the central part of the discharge, 1 mm in radius. Intensity fluctuations are caused by a huge number of very low-intensity spectral lines superimposed on the ICCD noise. Panel b): the ratio of the two P$_{\alpha }$ line profiles recorded in double (I2) and in single (I1) pass. In spite of noticeable scatter, specially on the line wings, it is easy to see that the main trend of the data is straight horizontal line. It means that the ratio of the two profiles is constant, i.e. self-absorption is not present.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8cm,clip]{12046F8.eps}
\vspace*{-3mm}
\end{figure} Figure 8:

The profile of the He II P$_{\alpha }$ line recorded axially (gray) at the 10th $\mu $s from the beginning of the discharge. The black profile originates from the light emitted only by the one contiguous layer. It is recorded side-on, at the same moment capturing emission of approximately the same volume that contributes in the axial view.

Open with DEXTER

4 Results and discussion

4.1 Electron density measurement

We mentioned above that the electron density (N) was estimated by relying on the known Stark FWHM of the He II P$_{\alpha }$(468.6 nm) spectral line. The inherent angle of view of the lens used to collect light was sufficient to cover whole cross section of the narrow part of the discharge tube. It means that the electron density evaluated in this way is an average value with respect to the actual density distribution along the tube radius. To estimate electron density on the axis of the tube, we reduced the angle of view to capture light only from the central part of the discharge. In this case the monitored area is no more than 1 mm in radius, along the whole discharge tube. A drawback of such a tiny angle of view is a significant decrease of signal/noise ratio. We chose 10th $\mu $s to record P$_{\alpha }$, the moment when the intensity of the line is still sufficient for reliable processing. According to our analysis, electron density measured on the axis of the tube is approximately 5% higher than the average value - the Lorentz FWHM of He II P$_{\alpha }$ at the axis is Wax=0.402 nm, while for the profile recorded integrally, the corresponding value is Wav=0.384 nm.

In the second stage of our analysis we discuss one possible issue of the self-absorption. To accomplish a check for self-absorption, we doubled the optical path length by placing a concave aluminum mirror behind the plasma (Hutchinson 2002). Two He II P$_{\alpha }$ profiles are subsequently recorded, with and without the auxiliary mirror. In Fig. 7 the corresponding P$_{\alpha }$profiles are presented together with their ratio after subtraction of the dark spectrum introduced by the ICCD. One can notice that there is no evidence of the self absorption and therefore optical depth for the $\lambda=468.6$ nm can be neglected at the 10th $\mu $s when P$_{\alpha }$ attains its maximum intensity[*].

Finally, we consider an influence of the contiguous layers of plasma, on both sides of the central narrow tube, on the procedure for electron density estimation. In Fig. 8 two P$_{\alpha }$profiles are shown, the first one in gray is recorded along the discharge tube axis with angle of view set to collect light integrally, and it includes contributions of the both contiguous layers. The second profile originates from the light emitted only by one contiguous layer. It is recorded side-on, capturing emission of approximately the same volume that contributes in the axial view. By subtracting contribution of both contiguous layers from the profile recorded axially, one can get corrected P$_{\alpha }$ profile, the basis for electron density estimation inside the narrow tube. The corrected profile is $\sim$6% wider, and the Lorentz FWHM is $W_{\rm corr}=0.436$ nm, which results in $\sim$7% higher electron density. Within this paper we have adopted corrected electron density values estimated from integral axial observations with uncertainties less than 15% due to errors related to the averaging over the tube radius and contribution of the contiguous plasma layers. Error of the best fit Lorentz FWHM of the P$_{\alpha }$ is low, about 1%, owing to a large number of data points in the line profile.

4.2 The Si I spectrum

Radial distribution of the Si I emitters should be close to the uniform, meaning the line profiles of Si I spectra are based on integral observations along the axis of the discharge tube. A certain amount of asymmetry in the Si I spectral line profiles is to be expected due to interaction of neutral silicon atoms with Si II, He II, and O II ions. In our experimental conditions, this effect is masked by overall noise present in the Si I spectra. Consequently, the Si I spectral lines are processed as symmetric lines. In Fig. 9 we present the best fit of the Si I 251.920 nm line to the Voigt profile, together with residuals. One can notice that data scatter completely mask the expected asymmetry of the line profile. Therefore, we were not able to isolate the ionic contribution to the line broadening. Because of that, values quoted in Table 2 represent the overall Stark FWHM ($W_{\rm m}$).

\begin{figure}
\par\includegraphics[width=7cm,clip]{12046F9.eps}
\end{figure} Figure 9:

The neutral silicon $\lambda =251.920$ nm spectral line, together with the best fit Voigt function. Expected asymmetry of the line is masked by the data scatter, so that the residuals are evenly distributed around the model function.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=7.45cm,clip]{12046F10.eps}
\end{figure} Figure 10:

Panels a-d). The measured and calculated Stark FWHM (in pm) versus electron temperature (T in K) at $1\times 10^{23}$ m-3 electron density for the most investigated Si II line profiles. The symbols mean: $\bullet $ this work; $\circ $(Gonzalez et al. 2002); $\square $ (Wollschläger et al. 1997); $\blacktriangle $ (Peréz et al. 1993); $\vartriangle $ (Lesage et al. 1983); + (Kusch & Schröder 1982); $\oplus $ (Chiang & Griem 1978); $\bigtriangledown $ (Lesage et al. 1977); $\blacksquare $ (Puric et al. 1974), and $\boxtimes $ (Konjevic et al. 1970). The SB, G and LDA denote values calculated by: Sahal-Bréchot (1968), Griem (1974), and Lanz et al. (1988) only for the electrons as perturbers. Error bars represent estimated accuracies in this work (see Table 3).

Open with DEXTER

Since the silicon is introduced in the working gas by sputtering of the glass one can expect generally low density of the silicon atoms. On the other hand, the investigated spectral lines of neutral silicon belong to the resonance transitions, so the issue of self-absorption has to be considered. Standard procedure for self-absorption check (Hutchinson 2002), applied in analysis of the He II P$_{\alpha }$, is inefficient here because the Si I lines are narrow, and corresponding profiles are recorded in no more than 10 points. Also, scatter of the data points here is more pronounced than for He I-II lines. Because of that it is hard to conclude whether two profiles, recorded with and without the auxiliary mirror, are proportional to the same factor in each point.

An independent procedure is employed to check for self-absorption. The essence is to calculate ratio of the measured intensities for two spectral lines, p and q, belonging to the same transition ( $R^{\rm m}_{p,q}=I_p/I_q$) and to compare it to the ratio $R^{\rm th}_{p,q}=A_pg_p/A_qg_q$ where A denotes the transition probability while g is the statistical weight. Ideally the both ratios should be the same, if the plasma is free of self-absorption. Because of the experimental errors associated to the line intensities (typically $\sim$6%) and uncertainties of the transition probabilities (mark C, up to the $\sim$25% in NIST 2009) the ratios will differ. Table 1 summarizes quotients of both ratios ( $Q_{p,q}=R^{\rm m}_{p,q}/R^{\rm th}_{p,q}$) for six Si I lines given in Fig. 1. One can notice that the worst disagreement is 33%, still within the margin imposed by experimental errors and uncertainties of the tabulated A values. Applying additional analysis presented in Appendix we have found that optical depth  $k_\lambda l$ for considered spectral lines is in the range $0.07<k_\lambda l<0.38$. Corresponding Stark widths are subsequently corrected to the limit of an optically thin layer (see Table 2).

Table 1:   Quotients Q for set of the six Si I spectral lines.

To check the contribution of the plasma layers beyond the narrow part of the discharge tube, we made the side-on recordings of the spectral regions with the most intense Si I lines. According to our findings, the intensity of the Si I lines originating from contiguous plasma layers is under the detectible limit, therefore their contributions to the profiles recorded axially can be neglected.

Table 2:   The Si I experimental Stark widths.

Our experimental results are presented in Table 2, together with existing results published by various authors. Three experiments are dedicated to investigating the resonance Si I lines W values (see Table 2), no one is made in helium plasma. Eight Si I Wvalues for high lying transitions (with $\lambda > 500$ nm), not researched in this work, are measured by Miller & Bengtson (1970). A review of calculated Si I W values is given in Griem (1974).

Table 3:   The measured Si II Stark widths.

Table 4:   The measured Si III Stark widths.

It is to be pointed out that the Si I spectrum has low intensity and is available only in a narrow window, between 22th $\mu $s and 24th $\mu $s after the beginning of the discharge, which corresponds to the noticeable minima in the electron temperature and electron density flow. We already mentioned that the minima of T and N appear at the moment when the discharge current changes its sign and has a low magnitude, see Fig. 5. Outside of this tiny window, the Si I line intensities are very low and the Si I spectrum is overpowered by Si III and Si IV spectra. Our new acquisition set-up provides accurate measurement of low-intensity spectral lines, and therefore Si I W values are more accurate than results by previous experimental technique.

4.3 The Si II spectrum

It should be mentioned that in the last four decades Si II was the fourth most experimentally investigated emitter, after hydrogen, helium, and argon. Therefore, the Si II Stark widths have been measured by many authors using different plasma sources (Konjevic et al. 1970; Puric et al. 1974; Lesage et al. 1977; Chiang & Griem 1978; Lesage et al. 1983; Kusch & Schröder 1982; Peréz et al. 1993; Wollschläger et al. 1997; Gonzalez et al. 2002; Lesage & Redon 2004). The critical reviews of the published experimental Stark widths are presented by Lesage (2009) and Konjevic et al. (2002, and references therein).

The calculated Si II W values (based on various theoretical approximations) are presented by Sahal-Bréchot (1968), Griem (1974) and Lanz et al. (1988). Our experimental Stark FWHMs $(W_{\rm m})$ are presented in Table 3, side by side with nine results published in the last forty years. The Stark FWHMs dependencies on the electron temperature for the most investigated Si II line profiles are presented in Fig. 10.

Here presented Si II W values agree with ones calculated by Lanz et al. (1988) for the most investigated transitions (see Fig. 10). It should be mentioned that at electron temperatures below 15 000 K values calculated by Lanz et al. (1988) and values presented by Griem (1974) differ for $25\%$ for the higher Si II multiplets. In the multiplet No. 1, the LDA W values are significantly below G values (see Fig. 10a), while our $W_{\rm m}$ values agree with the calculated LDA W values.

The Si II spectrum, as for Si I, is visible at the moment when Nand T have minimal values during the discharge. However, the time window for observation is almost twice wider, from 21th $\mu $s up to 26th $\mu $s.

4.4 The Si III spectrum

The Si III lines are the most intense spectral lines observed in our helium plasma source at this discharge condition. In particular, the well isolated 254.182 nm Si III line has been found to be the most intense during the whole plasma decay interval. It is also the most narrow among the investigated Si III lines (see Fig. 6 and Table 4). This line originates from the $\rm 3p^2~ ^1D_2$ Si III level with a 15.153 eV excitation energy, and represents a single allowed spontaneous de-excitation channel from this level (NIST 2009). Our measured Stark FWHMs $(W_{\rm m})$are presented in Table 4, together with other authors' results.

Some Si III W values were calculated by Dimitrijevic & Konjevic (1981) in the 3p-4s, 4s-4p, 3s-4p, 4p-4d, and 4p-5s transitions. Our Si III $W_{\rm m}$ values for ten lines show the best agreement (within 20%) with the calculated ones based on the modified semiempirical approximation.

4.5 The Si IV spectrum

The Si IV W values were experimentally investigated in two works and calculated by Dimitrijevic et al. (1991). Our measured Stark FWHMs $(W_{\rm m})$ are presented in Table 5, together with other authors' results.

Direct comparison of our results with calculated values presented by Dimitrijevic et al. (1991) is not possible due to the significant difference in electron temperatures. However, the values estimated by extrapolation indicate that, for the majority of lines, calculated Stark parameters are higher, about 27% in average, compared to our data. For two Si IV spectral lines (212.018 nm and 212.747 nm), though, agreement is encouraging - better than $\pm$5%. We would like to emphasize that existing experimental Si IV W values are also lower than the calculated by Dimitrijevic et al. (1991).

The Si III and Si IV spectra were recorded later, around the 30th $\mu $s from the beginning of the discharge. We have found that at this stage of the discharge Si III and Si IV spectral line profiles are well defined, while the continuum is inconsiderable (see Figs. 34, and 6).

5 Conclusions

The neutral and ionized silicon spectral line shapes were investigated in the laboratory helium plasma at electron densities ranging between $3.7\times 10^{22}$ m-3 and $1.1\times 10^{23}$ m-3 and electron temperatures between 12 500 K and 19 000 K. The spectral line profiles were recorded using a highly sensitive spectrograph (McPherson) and ICCD camera (Andor) detection system.

The Stark FWHMs of 13 Si I, 15 Si II, 28 Si III, and 9 Si IV spectral lines were measured in the wavelength interval between 206 nm and 640 nm. Five Si I, four Si II, eleven Si III values and one Si IV W value from the investigated set have not been measured or calculated so far. Our W values are compared with the existing theoretical and experimental data.

At the mentioned plasma parameters tolerable agreement was found (within the accuracy of the experiment and uncertainties of the theoretical approaches used) among our measured and calculated Stark FWHM values. We recommend the Stark FWHMs of the intense 254.182 nm, 308.624 nm and 309.342 nm Si III, and the 314.956 nm and 316.571 nm Si IV lines for purpose of plasma diagnostics.

Table 5:   The measured Si IV Stark widths.

Acknowledgements

This work is part of the ``Determination of atomic parameters on the basis of spectral line profiles'' (ON141010) project supported by the Ministry of Science and Technological Development of the Republic of Serbia.

Appendix A

It is interesting to estimate the maximum amount of self-absorption, which, if present, would still be consistent with our data. Several procedures for correcting spectral line profiles due to self-absorption are in use (Konjevic 1999; Escarguel et al. 2000). The fact that six Si I resonant spectral lines, considered in this paper, effectively have the same lower level gives us possibility for a specific analysis of the presented data.

It is well known that for homogenous plasma layer of length l intensity of the emitted radiation per unit wavelength (the radiance) is (Hilborn 1982):

\begin{displaymath}I_{\lambda}=F_{\lambda}[1-\exp(-k_{\lambda}l)]
\end{displaymath} (1)

where $F_{\lambda}$ denotes the radiance per unit length per unit wavelength, while $k_{\lambda}$ is the absorption coefficient at the wavelength $\lambda$. Product $k_\lambda l$ is the optical depth - the quantity which determines the amount of self-absorption. For very weak self-absorption, $k_{\lambda}l\ll
1$, after preserving only the first two terms in the Maclaurin expansion of the exponential function $\exp(k_{\lambda}l)\sim 1-
k_{\lambda}l$, the radiance of homogenous plasma layer simplifies to $I_{\lambda}=F_{\lambda}l\equiv I_{0\lambda}$. Suppose that self-absorption is slightly higher, so that we need first three terms in the expansion of the exponential function. Equation (1) becomes

\begin{displaymath}I_{\lambda}=I_{0\lambda}[1-k_{\lambda}l/2].
\end{displaymath} (2)

Coefficient $k_{\lambda}$ depends on oscillator strength f12and density of the lower level of the transition n, i.e. $k_{\lambda}=\frac{e^2\lambda^2}{4\epsilon_0mc^2}g_1f_{12}nz(\lambda)$, where e is the electron charge, $\epsilon_0$ the vacuum permittivity, c the speed of light, m the electron mass, g1the statistical weight of the lower level of the transition and $z(\lambda)$ the normalized profile of the spectral line, $\int_{-\infty}^{\infty}z(\lambda){\rm d}\lambda=1$. Inserting $k_{\lambda}$ into Eq. (2) and integrating over the line profile we have

\begin{displaymath}I=I_0\left[1-\frac{e^2\lambda^2l}{8\epsilon_0mc^2}g_1f_{12}n\int_{-\infty}^{\infty}z\:^2(\lambda){\rm d}\lambda\right].
\end{displaymath} (3)

In Eq. (3) I is the reduced spectral line intensity due to self-absorption, while I0 is the one that would be obtained in a case of negligible self-absorption. We introduce quantity $Z^{\:2}=\int_{-\infty}^{\infty}z\:^2(\lambda){\rm d}\lambda$. For common line profiles $Z^{\:2}=(C\cdot W)^{-1}$ where W is the halfwidth of the line, C is the dimensionless constant equal to $2\pi$ for Lorentz profile and $\sqrt{2\pi/\ln2}$ for Gauss. Usually $z(\lambda)$ is not known, but for a weak self absorbtion, one can take the experimental value for W sacrifiing a few percent in exactness. If we have two measured intensities, for spectral lines p and q with the same lower level, their ratio will be

\begin{displaymath}\frac{I_p}{I_q}=\frac{I_{0p}}{I_{0q}}\left[\frac
{1-\frac{e^2...
...a_q^2\:l}{8\epsilon_0mc^2}g_qf_{q}\;n\:Z_q^{\:2}}
\right]\cdot
\end{displaymath} (4)

Here we applied a simplified notation for oscillator strengths $f_p\equiv f^p_{12}$ and $f_q\equiv f^q_{12}$, in a similar way $g_p\equiv g^p_1$ and $g_q\equiv g^q_1$. For lines from the same multiplet, $\frac{I_{0p}}{I_{0q}}=
\frac{A_p\:g_p\:\lambda_q}{A_q\:g_q\:\lambda_p}$, where Ap, Aq are transition probabilities for the lines p and q. From Eq. (4) one can express density of the lower level of the transition, n, as

\begin{displaymath}n=\frac{Q_{p,q}-1}{\bigl(Q_{p,q}\frac{1}{W_q}g_qf_q\lambda_q^...
..._p}g_pf_p\lambda_p^2\bigr)}
\frac{8C\epsilon_0mc^2}{e^2l}\cdot
\end{displaymath} (5)

Quantity $Q_{p,q}=\frac{I_p}{I_q}\frac{A_q\:g_q\:\lambda_p}{A_p\:g_p\:\lambda_q}$is already defined in Sect. 4.2 and tabulated in Table 1. Taking into account that, for lines from the same multiplet $\lambda_1\sim\lambda_2\sim\lambda_3\cdots=\lambda$, Eq. (5) simplifies to

\begin{displaymath}n=\frac{Q_{p,q}-1}{\bigl(Q_{p,q}\frac{1}{w_q}g_qf_q-\frac{1}{w_p}g_pf_p\bigr)}
\frac{8CW_0\epsilon_0mc^2}{e^2l\;\lambda^2},
\end{displaymath} (6)

where wp, wq are dimensionless numbers expressing walfwidths of the lines p and q in picometers, i.e. $W_{i}=w_{i}\cdot
W_0$ where i = {pq} and $W_0=1\times10^{-12}$ m. There are N=6 independent values of quotient Q for the six Si I resonant lines[*] given in Table 1, which will result in 6 values for n. Ideally, all values should be the same. Due to uncertainties of measured intensities I and halfwidths of the line profiles W and uncertainties of tabulated gf data one will get a set of 6 different values for n. We can easily calculate the mean value $\overline{n}=\frac{\sum_i n_i}{N}$ and the error of the mean value $\Delta n=\sqrt{\frac{\sum_i(\overline{n}-n_i)^2}{N(N-2)}}$of the set. In practice, one should remove obvious outliers from the data set before calculating $\overline{n}$$\Delta n$. For the data presented in Table 1 and adopting[*] C=4, we found $\overline{n}=(0.95\pm0.83)\times 10^{17}\:$m-3. This result allows calculation of optical depth for each spectral line $\;\;k_{\lambda}l=\frac{e^2\lambda^2\:l}{4\epsilon_0mc^2}g_1f_{12}\overline{n}z(\lambda)$. Obviously, the amount of the self-absorption for a particular line is proportional to the product gf. It is common practice to apply adequate analytical or numerical procedure to correct spectral line profile if the optical depth $k_\lambda l<1$. For the Si I resonant spectral lines optical depth, at our discharge conditions, is in the range $0.07<k_\lambda l<0.38$. Therefore, by relying on Eq. (2) one can correct measured radiance $I_\lambda$ of the spectral line at the particular wavelength $\lambda$ to be close to the value $I_{0\lambda}$ which would be obtained in the case of the negligible self-absorption

\begin{displaymath}I_{o\lambda}=\frac{I_\lambda}{1-\frac{k_\lambda\:l}{2}}\cdot
\end{displaymath} (7)

Relying on Eq. (7) we evaluated corrections of the Stark widths given in Table 2.

References

  • Arsenijevic, V., Fabbro, S., Mourao, A. M., & Rica da Silva, A. J. 2008, A&A, 492, 535 [NASA ADS] [CrossRef] [EDP Sciences]
  • Bongard, S., Baron, E., Smadja, G., Branch, D., & Hauschildt, P. H. 2008, ApJ, 687, 456 [NASA ADS] [CrossRef]
  • Bukvic, S., & Spasojevic, D. 2005, Spectrochim. Acta B, Atomic Spectr., 60, 1308
  • Bukvic, S., Spasojevic, D., & Zigman, V. 2008, A&A, 477, 967 [NASA ADS] [CrossRef] [EDP Sciences]
  • Catanzaro, G., & Leone, F. 2008, MNRAS, 389, 1414 [NASA ADS] [CrossRef]
  • Catanzaro, G., Leone, F., Busa, I., & Romano, P. 2008, New Astron., 13, 113 [NASA ADS] [CrossRef]
  • Chiang, T., & Griem, H. R. 1978, Phys. Rev. A, 18, 1169 [NASA ADS] [CrossRef]
  • Danforth, C. W., & Shull, M. J. 2008, ApJ, 679, 194 [NASA ADS] [CrossRef]
  • Davies, J. T., & Vaughan, J. M. 1963, ApJ, 137, 1302 [NASA ADS] [CrossRef]
  • Dimitrijevic, M. S., & Konjevic, N. 1981, Spectral Line Shapes, ed. B. Wende (Berlin: W. de Gruyter)
  • Dimitrijevic, M. S., Sahal-Bréchot, S., & Bommier, V. 1991, A&ASS, 89, 591 [NASA ADS]
  • Dimitrijevic, M. S., Ryabchikova, T., Simic, Z., Popovic, L. C., & Dacic, M. 2007, A&A, 469, 681 [NASA ADS] [CrossRef] [EDP Sciences]
  • Djenize, S. 2009, Spectrochim. Acta B, Atomic Spectr., 64, 242
  • Djenize, S., & Bukvic, S. 2001, A&A, 365, 252 [NASA ADS] [CrossRef] [EDP Sciences]
  • Djenize, S., Sreckovic, A., Labat, J., Puric, J., & Platisa, M. 1992, J. Phys. B, 25, 785 [NASA ADS] [CrossRef]
  • Djenize, S., Dimitrijevic, M. S., Sreckovic, A., & Bukvic, S. 2002, A&A, 396, 331 [NASA ADS] [CrossRef] [EDP Sciences]
  • Drawin, H. W., & Felenbok, P. 1965, Data for Plasmas in Local Thermodynamic Equilibrium (Paris: Gauthier-Villars)
  • Escarguel, A., Ferhat, B., Lesage, A., & Richou, J. 2000, JQSRT, 64, 353 [NASA ADS]
  • Fox, A. J., Ledoux, C., Vreeswijk, P. M., Smette, A., & Jaunsen, A. O. 2008, A&A, 491, 189 [NASA ADS] [CrossRef] [EDP Sciences]
  • Gomez de Castro, A. I., & Verdugo, E. 2007, ApJ, 654, L91 [NASA ADS] [CrossRef]
  • Gonzalez, V. R., Aparicio, J. A., del-Val, A., & Mar, S. 2000, A&A, 363, 1177 [NASA ADS]
  • Gonzalez, V. R., Aparicio, J. A., del-Val, A., & Mar, S. 2002, J. Phys. B: At. Mol. Opt. Phys., 35, 3557 [NASA ADS] [CrossRef]
  • Griem, H. R. 1964, Plasma Spectroscopy (New York: McGraw Hill)
  • Griem, H. R. 1974, Spectral Line Broadening by Plasmas (New York: Acad. Press)
  • Hamann, F., Kaplan, K. F., Hidalgo, R. P., Prochaska, J. X., & Herbert-Fort, S. 2008, MNRAS, 391, L39 [NASA ADS] [CrossRef]
  • Hilborn, R. C. 1982, Am. J. Phys., 50, 982 [NASA ADS] [CrossRef]
  • Hutchinson, I. H. 2002, Principles of plasma diagnostics (Cambridge University Press), second edition
  • Iijima, T., & Nakanishi, H. 2008, A&A, 482, 865 [NASA ADS] [CrossRef] [EDP Sciences]
  • Kawka, A., Vennes, S., Dupuis, J., Chayer, P., & Lanz, T. 2008, ApJ, 675, 1518 [NASA ADS] [CrossRef]
  • Konjevic, N. 1999, Phys. Rep., 316, 339 [NASA ADS] [CrossRef]
  • Konjevic, N., Puric, J., Cirkovic, Lj., & Labat, J. 1970, J. Phys. B, 3, 999 [NASA ADS] [CrossRef]
  • Konjevic, N., Lesage, A., Fuhr, J. R., & Wiese, W. L. 2002, J. Phys. Chem. Ref. Data, 31, 819 [NASA ADS] [CrossRef]
  • Kruithof, A. A., & Penning, F. M. 1937, Physica, 4, 430 [NASA ADS] [CrossRef]
  • Kusch, H. J., & Schröder, K. 1982, A&A, 116, 255 [NASA ADS]
  • Lanz, T., Dimitrijevic, M. S., & Artru, M.-C. 1988, A&A, 192, 249 [NASA ADS]
  • Lesage, A. 1994, Astrophysical Applications of Powerful New Databases, ed. S. J. Adelman, & W. L. Wiese, Astron. Soc. Pacific Conf. Ser., 78
  • Lesage, A. 2009, New Astron. Rev., 52, 471 [NASA ADS] [CrossRef]
  • Lesage, A., & Redon, R. 2004, A&A, 418, 765 [NASA ADS] [CrossRef] [EDP Sciences]
  • Lesage, A., Sahal-Bréchot, S., & Miller, M. H. 1977, Phys. Rev. A, 16, 1617 [NASA ADS] [CrossRef]
  • Lesage, A., Rathore, B. A., Lakicevic, I. S., & Puric, J. 1983, Phys. Rev. A, 28, 2264 [NASA ADS] [CrossRef]
  • Meyer, J., & Beck, R. J. 1970, A&A, 8, 93 [NASA ADS]
  • Miller, M. H., & Bengtson, R. D. 1970, Phys. Rev. A, 1, 983 [NASA ADS] [CrossRef]
  • Milosavljevic, V., & Djenize, S. 2001, Eur. Phys. J. D, 15, 99 [NASA ADS] [CrossRef] [EDP Sciences]
  • Nehme, C., Gry, C., Boulanger, F., et al. 2008, A&A, 483, 471 [NASA ADS] [CrossRef] [EDP Sciences]
  • NIST 2009, Atomic Spectra Database Lines (wavelength order), http://www.physics.nist.gov
  • Peréz, C., de la Rosa, I., de Frutos, A. M., & Mar, S. 1993, Phys. Rev. E, 47, 756 [NASA ADS] [CrossRef]
  • Platisa, M., Dimitrijevic, M. S., Popovic, M., & Konjevic, N. 1977, J. Phys. B, 10, 2997 [NASA ADS] [CrossRef]
  • Puric, J., Djenize, S., Labat, J., & Cirkovic, Lj. 1974, Z. Phys., 267, 71 [NASA ADS] [CrossRef]
  • Rogowski, W., & Steinhaus, W. 1912, Archiv für Elektrotechnik, 1, 141 [CrossRef]
  • Rompe, R., & Stenbeeck, M. 1967, Ergebnisse der Plasmaphysik und der Gaselektronik (Berlin: Akademie Verlag)
  • Sahal-Bréchot, S. 1968, Z. Astrophys., 69, 74 [NASA ADS]
  • Shi, J. R., Gehren, T., Butler, K., Mashonkina, L. I., & Zhao, G. 2008, A&A, 486, 303 [NASA ADS] [CrossRef] [EDP Sciences]
  • Silfvast, W. T. 1971, Phys. Rev. Lett., 27, 1489 [NASA ADS] [CrossRef]
  • Sreckovic, A., Bukvic, S., & Djenize, S. 1998, Phys. Scr., 57, 225 [NASA ADS] [CrossRef]
  • Sundqvist, J. O., Ryde, N., Harper, G. M., Kruger, A., & Richter, M. J. 2008, A&A, 486, 985 [NASA ADS] [CrossRef] [EDP Sciences]
  • Wollschläger, F., Mitsching, J., Meiners, D., et al. 1997, JQSRT, 58, 135 [NASA ADS]

Footnotes

... intensity[*]
Uncertainty of the optical depth estimated by this method is proportional to $\sigma/\!\!\sqrt{N}$ where $\sigma$ is the standard deviation of the data scatter and N number of data points. Because of the nonsymmetric scatter we applied recently developed DLS method for data fitting (Bukvic et al. 2008) to estimate relevant quantities. It follows that the uncertainty of the optical depth is $\sim$0.009. In other words, the maximum optical depth consistent with our data is $\leq $0.009.
... lines[*]
In general for N spectral lines there are N!/(2(N-2)!) independent combinations; however, we have two groups of three spectral lines measured at two different temperatures. We suppose that the density of the effective lower level of the multiplet is the same for both temperatures.
... adopting[*]
Constant Cfor Voigt profile is less than $2\pi$ (Lorentz) and above $\sqrt{2\pi/\ln2}$ (Gauss), and it has to be evaluated separately for a specific Voigt profile.

All Tables

Table 1:   Quotients Q for set of the six Si I spectral lines.

Table 2:   The Si I experimental Stark widths.

Table 3:   The measured Si II Stark widths.

Table 4:   The measured Si III Stark widths.

Table 5:   The measured Si IV Stark widths.

All Figures

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12046Fg1.eps}
\end{figure} Figure 1:

The recorded resonance Si I line profiles (see Table 2). The numbers in brackets denote relative line intensities tabulated by NIST (2009). The resonance B I line profiles are also presented.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12046F2.eps}
\end{figure} Figure 2:

Panels a-c) The recorded profiles of prominent Si II spectral lines present in many astrophysical spectra. The numbers in brackets denote relative line intensities tabulated by NIST.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12046F3.eps}
\end{figure} Figure 3:

Panels a) and b) The recorded profiles of prominent Si III spectral lines present in many astrophysical spectra. The numbers in brackets denote the relative line intensities tabulated by NIST.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12046F4.eps}
\end{figure} Figure 4:

The recorded profiles of the 314.956 nm and 316.571 nm Si IV spectral lines. The numbers in brackets denote relative line intensities tabulated by NIST.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12046Fg5.eps}
\end{figure} Figure 5:

Temporal evolutions of the electron density (N in 1023 m-3) and electron temperature (T in 104 K). The error bar represents ${\pm } 14\%$ estimated uncertainties. The gray line is the discharge current (in arbitrary units) captured by the Rogowski coil.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12046Fg6.eps}
\end{figure} Figure 6:

The 254.182 nm Si III line profile recorded at $N = 5.2 \times 10^{22}$ m-3 and T= 19 000 K. Circles represent measured values. The solid line represents the corresponding Voigt profile using 4.4 pm for the Gaussian parameter (Doppler + instrumental convolution) in the fitting procedure. The total FWHM ($W_{\rm t}$) (Lorentzian + Gaussian convolution) is 12.02 pm with the Stark component of $W_{\rm m}=3.4 $ pm (see Table 4). The baseline estimation is made using the procedure described by Bukvic & Spasojevic (2005) and Bukvic et al. (2008).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12046F7.eps}
\end{figure} Figure 7:

Panel a): the He II P$_{\alpha }$ line recorded in single pass (black) and in double pass (gray). Both profiles are taken at the 10th $\mu $s and same detection conditions. The light is captured only from the central part of the discharge, 1 mm in radius. Intensity fluctuations are caused by a huge number of very low-intensity spectral lines superimposed on the ICCD noise. Panel b): the ratio of the two P$_{\alpha }$ line profiles recorded in double (I2) and in single (I1) pass. In spite of noticeable scatter, specially on the line wings, it is easy to see that the main trend of the data is straight horizontal line. It means that the ratio of the two profiles is constant, i.e. self-absorption is not present.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12046F8.eps}
\vspace*{-3mm}
\end{figure} Figure 8:

The profile of the He II P$_{\alpha }$ line recorded axially (gray) at the 10th $\mu $s from the beginning of the discharge. The black profile originates from the light emitted only by the one contiguous layer. It is recorded side-on, at the same moment capturing emission of approximately the same volume that contributes in the axial view.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7cm,clip]{12046F9.eps}
\end{figure} Figure 9:

The neutral silicon $\lambda =251.920$ nm spectral line, together with the best fit Voigt function. Expected asymmetry of the line is masked by the data scatter, so that the residuals are evenly distributed around the model function.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.45cm,clip]{12046F10.eps}
\end{figure} Figure 10:

Panels a-d). The measured and calculated Stark FWHM (in pm) versus electron temperature (T in K) at $1\times 10^{23}$ m-3 electron density for the most investigated Si II line profiles. The symbols mean: $\bullet $ this work; $\circ $(Gonzalez et al. 2002); $\square $ (Wollschläger et al. 1997); $\blacktriangle $ (Peréz et al. 1993); $\vartriangle $ (Lesage et al. 1983); + (Kusch & Schröder 1982); $\oplus $ (Chiang & Griem 1978); $\bigtriangledown $ (Lesage et al. 1977); $\blacksquare $ (Puric et al. 1974), and $\boxtimes $ (Konjevic et al. 1970). The SB, G and LDA denote values calculated by: Sahal-Bréchot (1968), Griem (1974), and Lanz et al. (1988) only for the electrons as perturbers. Error bars represent estimated accuracies in this work (see Table 3).

Open with DEXTER
In the text


Copyright ESO 2009

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.