Pressure Broadening of Some He I Lines

Quantum statistical approach is adopted for calculating the spectral line shapes of neutral helium in dense plasmas. Stark broadening of isolated He I lines 5048 Å (21P − 41S), 3889 Å (23S − 33P), and 3188Å (23S − 43P) is presented. Based on thermodynamic Green’s function, the electronic contribution to the shift and width is considered. The participation of ions to the line broadening is treated in a quasistatic approximation, by taking both quadratic Stark effect and quadrupole interaction into account. The calculated shifts and widths are compared with existing data.


Introduction
Plasma spectroscopy deals with the characteristics of radiation emitted from a plasma. In dense plasmas, damping of the emitted radiation occurs by means of several mechanisms; the most effective one is pressure broadening (Stark broadening). The interaction between a radiating atom and surrounding perturbing particles leads to Stark broadening. High-speed electrons perturb the emitter by collisions, causing the interruption of the spontaneous emission and altering the emitter energy levels [1][2][3][4].
Line profile calculation is an interesting tool for both laboratory and astrophysical plasma diagnostics, for example, to determine the internal parameters, to understand the microscopic processes within the plasma, and to check the quality of the predicted experimental and theoretical parameters [1,2].
The emission spectra of helium and He-like ions with their simple atomic structures are interesting for plasma diagnostics such as in shock wave tube or pulsed arc plasmas [5,6] and in the astrophysical context, for example, stellar atmospheres of hot stars and white dwarfs [7][8][9]. Helium is used as a carrier gas in many laboratories and is weekly interacting with materials and less harmful for plasma facing components than hydrogen and its isotopes [5]. The Helike ions may exist even at extremely high temperatures and densities. Even ITER is started with helium discharge, also He can be observed in discharge of JET. Spectroscopic measurements of tokamak plasmas are not free from helium [10].
Thermodynamic Green's function approach is a powerful tool to describe the Stark broadening [27][28][29][30]. In the last two decades, a quantum statistical approach has been developed, taking into account the medium effects by using Green's function [31][32][33][34][35][36]. In principle, this approach is able to describe dynamical screening and strong collisions by electrons, as well as the dynamic ion microfield, in a systematic way. This quantum statistical approach has been successfully applied to calculate the spectral line shapes of hydrogen, helium, H-, and He-like ions in dense plasmas [37][38][39][40][41][42].
In this paper thermodynamic Green's function approach is considered to calculate the pressure broadening of some selected neutral helium lines. In Section 2 the definitions and spectral properties of Green's function approach to spectral 2 International Journal of Spectroscopy line shapes of nonideal plasmas are presented. A review of the basic formalism is introduced and extended to helium lines. Section 3 provides the calculated shifts and full widths at half maximum (FWHM) for nonoverlapping, isolated (nondegenerate) He I lines 5048Å (2 1 P −4 1 S), 3889Å (2 3 S−3 3 P), and 3188Å (2 3 S−4 3 P) in dense plasmas. Finally, conclusions are given in Section 4.

Theoretical Calculations
Microscopic formation of the spectral line shapes in dense plasmas arises from perturbation of the radiative atom by collective and many-body effects [1,4]. Thus, the interaction with the surrounding particles must be taken into account. For example, screening is considered as an important collective effect in plasmas. The influence of electrons and ions may be treated separately due to the difference in mass and mobility. Green's function methods provide a perturbative approach to correlation functions and quantum effects of many-body systems. Two-particle Green's function is used to calculate the line shape and electron broadening from the self-energy and the vertex function. From statistical properties of the system, the current-current correlations determine the absorption spectrum, which is utilized by Ross [27]. Recently, further improvements of this approach have been made by Röpke et al. [32], Hitzschke et al. [33], Günter et al. [34], and Omar et al. [41].
Optical properties of many-particle systems are specified by the dielectric function ε(q, ω) based on Green's function theory, which is the response of the medium to an external electromagnetic field. The longitudinal dielectric function is related to the polarization function Π(q, ω) where V (q) = e 2 /( 0 q 2 ) is the Fourier transformed Coulomb potential. In terms of the absorption coefficient α(ω) and the index of refraction n(ω), the transverse dielectric function ε t (q, ω) in the long wavelength limit q → 0 reads In the visible region where the wavelength λ is large compared with the atomic dimension a B , the transverse and the longitudinal part of dielectric function coincide. The absorption coefficient is proportional to the imaginary part of the dielectric function where A denotes (Imε(q, ω)) 2 . In thermal equilibrium, the absorption coefficient is related to the emission coefficient by Kirchoff 's law [43]. In optically thin plasma, the emission coefficient is proportional to the line emission. As mentioned above, the dielectric function is related to the polarization function Π(q, ω). Then the medium modifications of spectral line shapes can be addressed to the bound-bound twoparticle polarization function, which concerns to dipoledipole autocorrelation function [36]. The perturber-radiator interaction leads to pressure broadening, which contains electronic and ionic contributions. Describing the ionic contribution in the quasistatic approximation by averaging over the ionic microfield at radiating atom [31,36,44], we get Here, the ionic microfield distribution function P(β) is taken according to the Hooper field distribution, and β = E/E 0 is the normalized field strength [45]. So Dω i f = E i − E f is the unperturbed transition energy between the initial i and the final f states; i and f are the corresponding intermediate states where i|r| f is identified as a dipole matrix-element for the transition between i and f states. The line profile itself is determined by the vertex correction Γ V i f for the overlapping lines and by the self-energy corrections Σ of the initial and final states Both electronic and ionic contributions occur in the selfenergy Σ n (ω, β), which is assumed to be diagonal in the atomic state n; The electronic self-energy is obtained by performing a Born approximation with respect to the perturber-radiator interaction [36] Here, the sum over α runs from n − 2 to n + 2 discrete bound states for virtual transitions, and the level splitting due to the International Journal of Spectroscopy 3 ion microfield E α (β) has been neglected [35]. So n B (ω) = [exp (Dω/k B T) − 1] −1 is the Bose distribution function, and M nα (q) are the transition matrix-elements, given below. The inverse dielectric function ε −1 (q, ω) contains many particle effects which account for the dynamical screening of the interaction in the plasma: The random phase approximation (RPA) for the dielectric function is used: where E p = D 2 p 2 /2m e is the kinetic energy of electrons, and f e (E p ) is the Fermi distribution function of the electrons, approximated to the Boltzmann distribution function in the degenerate limit The full expression of the inverse dielectric function has to be used if the transition frequency ω nα becomes comparable to the electron plasma frequency ω pl = (n e e 2 / 0 m e ) 1/2 . However, in the high-frequency limit ω nα ω pl , the inverse dielectric function can be approximated by This binary collision approximation leads to a linear behavior of the electronic shift contribution with respect to the electron density, whereas a nonlinear dependence of the electronic shift with increasing electron density is expected if the full expression of the inverse dielectric function is used [33,35]. The vertex function for the coupling between the upper and the lower state is given by The transition matrix-elements M nα (q) describe the coupling between free charges and bound states. In lowest order, they are determined by the atomic eigenfunctions ψ n (p) of the radiating electron and depend on the momentum transfer Dq [34,46] assuming that the ion with charge Z is much heavier than the electron m i m e : Expanding the plane wave into spherical harmonics where j l (qr) is the Bessel function, a multipole expansion can be derived; for example, l = 0, 1, 2 gives the monopole, dipole, and quadrupole contribution of the radiator-electron interaction, respectively. The radial part of helium wave function is calculated based on Coulomb approximation method of Bates and Damgaard [47][48][49]. For more detail see [41]. By using the Born approximation, the electronic selfenergy is overestimated. To avoid this we apply a cutoff procedure and add the strong collision term in contrast to partial summation of the three-particle T-matrix, which is quite suitable for treating short-range interactions between particles [50], where the result might be slightly modified. According to Griem, the cutoff parameter for the qintegration is the inverse of the minimum limiting impact parameter (q max = 1/ρ min ) [22,36,49,51].
To determine the ionic self-energy, we approximate the time-dependent microfield fluctuation by its static value. In general, dynamic ionic microfield is important for overlapping lines and at low electron density in the line center [52,53]. Due to the slow movement of heavy ions, the ion microfield is assumed to be constant during the time of interest for the radiation process. The static ionic contribution to the ionic self-energy is treated by means of the microfield concept including both quadratic Stark effect and quadrupole effects. The first-order perturbation term vanishes for nonhydrogenic like atoms because of nondegeneracy with respect to the orbital quantum number l. According to second-order perturbation theory, the quadratic Stark effect is proportional to the square of the microfield [54] where E is the microfield strength; n, l, and m are the wellknown principal, orbital, and magnetic atomic quantum numbers, respectively. The quadrupole Stark effect is due to the inhomogeneity of the ionic microfield. We use the expression derived by Halenka [55]: Here, B ρ (β) is the mean field gradient at a given field strength, and the screening parameter ρ = r 0 /r D is taken as the ratio between the mean particles distance r 0 and the Debye radius r D .

Results and Discussions
The Stark broadening parameters for the transition (2 1 P − 4 1 S) 5048Å are given in Table 1 [66], measurements of Kelleher [56], and measurements of Diatta [67]. Our calculations agree better with the measured value by Kelleher [56] in contrast to the results given by Bassalo et al. [66]. The measured shift by Diatta [67] is obviously smaller than the calculated shifts given in Table 1.
The Stark width and shift of the line (2 3 S − 3 3 P) 3889Å are measured by Pérez et al. [6]. The measured values were in the plasma density range of (1 − 6) × 10 22 m −3 and temperature interval of (0.8−3)×10 4 K with a mean value of 2 × 10 4 K. The error bar in the case of n e was ±10%, and the uncertainty in the temperature evaluation was about 20%. Recently, the FWHM of this line is measured by Gao et al. [61] for a helium arc for density range (0.5 − 4) × 10 22 m −3 . Figures 1 and 2 include other available experimental [56][57][58][59][60]63] and theoretical results [22,62]. The MD simulation results of Gigosos et al. [62] have been performed for independent as well as interacting particles in nonquenching approximation. Our calculations are also presented, the width shows a good agreement especially with the MD simulations data of Gigosos et al. [62], where no Doppler broadening is included. Nonlinearity can be seen at very high electron density. The discrepancy between the measured and calculated line broadening may be related to self absorption [61]. However, our results for the shift are overestimated, on the other hand better agreement can be seen by comparing our results in non-quenching approximation with the results of Gigosos et al. [62], which give lower values of shift at high electron density.
The Stark parameters of the line (2 3 S − 4 3 P) 3188Å are measured by Peláez et al. [64]. They made a spectroscopic and interferometric analysis of a pulsed plasma. The electron temperature (from 1900 K to 2300 K) was obtained from the intensity ratio of He II lines. Electron density was determined by interferometry, ranging from 1.25 × 10 22 m −3 to 6.22 × 10 22 m −3 . From these experimental results an empirical calibration for the Stark parameters was obtained in a broad range of electron densities [64]. The experimental and various theoretical Stark parameters are reported by Peláez et al. [64], and our results are included in Figures 3 and  4 for Stark width and shift as a function of electron density, respectively, at the given plasma temperatures. Theoretical predictions of Bassalo et al. (BCW) [24] can also be seen in Figures 3 and 4. Further experimental results are presented, carried out by Kelleher [56], Mijatovic et al. [65], Berg et al. [58], and Soltwisch and Kusch [60]. The measured value of Berg et al. [58] was compared with Griem's theory [58] [22,24,56,58,60,64,65]. at electron density (1.5 ± 0.08) × 10 23 m −3 and temperature (2.9 ± 0.2) × 10 4 K.

Conclusions
The quantum statistical approach has been developed to calculate spectral line shapes in dense plasmas. By using thermodynamic Green's function, a systematic perturbative treatment of the polarization function has been performed [33,34,46]. In contrast to the molecular dynamics (MD) simulations, consistent quantum description is applied here to calculate the Stark parameters by using the formalism presented above. The calculated Stark shift and full width at half maximum (FWHM) of He I lines 5048Å (2 1 P − 4 1 S), 3889Å (2 3 S − 3 3 P), and 3188Å (2 3 S − 4 3 P) have been calculated in the density range (10 21 − 10 24 ) m −3 and for temperatures between (0.5 − 6) × 10 4 K. In dense plasmas the binary (few-particle) collision approximation is not appropriate to describe the collective interactions; therefore at high densities the dynamical screening effect will be more important seen (12), which reduces the linearity of Stark parameters with increasing electron density. This affects mostly the shift than the width, and it can be seen in Figures 2 and 4. The Hooper microfield distribution function is applicable for plasmas in the regime Γ ≤ 1, in which the correlations between plasma particles are considered as small perturbations. This distribution is not applicable in strongly coupling regime, where the correlation effect is important at high density, and therefore a new distribution function should be adopted. Generally, our calculated line widths show good agreement by comparing with other results, while the shift is slightly overestimated for the lines 3889Å and 3188Å. However, by performing nonquenching approximation better agreement can be seen. 6 International Journal of Spectroscopy Moreover, our quantum statistical approaches can be applied not only to investigate the line shapes of twoelectron atom but also to complex atoms. Furthermore, the deformation of spectral line profiles can be investigated in a strong magnetic field for stellar diagnostic.