Improving the Method of Measuring the Electron Density via the Asymmetry of Hydrogenic Spectral Lines in Plasmas by Allowing for Penetrating Ions

There was previously proposed and experimentally implemented a new diagnostic method for measuring the electron density Ne using the asymmetry of hydrogenic spectral lines in dense plasmas. Compared to the traditional method of deducing Ne from the experimental widths of spectral lines, the new method has the following advantages. First, the traditional method requires measuring widths of at least two spectral lines (to isolate the Stark broadening from competing broadening mechanisms), while for the new diagnostic method it is sufficient to obtain the experimental profile of just one spectral line. Second, the traditional method would be difficult to implement if the center of the spectral lines was optically thick, while the new diagnostic method could still be used even in this case. In the theory underlying this new diagnostic method, the contribution of plasma ions to the spectral line asymmetry was calculated only for configurations where the perturbing ions were outside the bound electron cloud of the radiating atom/ion (non-penetrating configurations). In the present paper, we take into account the contribution to the spectral line asymmetry from penetrating configurations, where the perturbing ion is inside the bound electron cloud of the radiating atom/ion. We show that in high-density plasmas, the allowance for penetrating ions can result in significant corrections to the electron density deduced from the spectral line asymmetry.


Introduction
In medium-density plasmas, profiles of hydrogenic spectral lines look symmetric, but in high-density plasmas, they become asymmetric.This asymmetry is caused primarily by the nonuniformity of the ion microfield, as noted by Sholin and his co-workers in papers [1][2][3]-for the latest advances in the theory of the asymmetry we refer to papers [4,5] and the references therein, of which we especially note papers [6,7].(There are also secondary sources of the asymmetry, as discussed in more detail below in the first paragraph of Section 2).Often, the blue maximum of the spectral line is higher than the red maximum, and the positions of the intensity maxima are asymmetrical with respect to the unperturbed line center.
A new diagnostic method for measuring the electron density using the asymmetry of hydrogenic spectral lines in dense plasmas was proposed and implemented in paper [8].In that paper, in particular, from the experimental asymmetry of the C VI Lyman-delta line emitted by a vacuum spark discharge, the electron density was deduced to be N e = 3 × 10 20 cm −3 .This value of N e was in good agreement with the electron density determined from the experimental widths of C VI Lyman-beta and Lyman-delta lines.
Later, this diagnostic method was also employed in the experiment presented in paper [9].In that laser-induced breakdown spectroscopy experiment, the electron density N e ~3 × 10 17 cm −3 was determined from the experimental asymmetry of the H I Balmer-beta (H-beta) line.
This new diagnostic method has the following advantages compared to the method of deducing N e from the experimental widths of spectral lines.First, the latter, traditional, method requires measuring widths of at least two spectral lines-because the widths are affected not only by the Stark broadening, but also by competing broadening mechanisms, such as, e.g., Doppler broadening.In distinction, to use the new diagnostic method, it is sufficient to obtain the experimental profile of just one spectral line-because the Doppler broadening does not cause the asymmetry.
Second, the traditional method based on experimental widths would be difficult to implement if the center of the spectral lines were optically thick.In distinction, the new diagnostic method can still be used even if the spectral line is optically thick in its central part.This is because the overwhelming contribution to the asymmetry originates from the wings of the spectral line, the wings usually being optically thin.More details can be found in Section 1.6 of [10]  1 .
In the theory underlying this new diagnostic method, the contribution of plasma ions to the spectral line asymmetry was calculated only for configurations where the perturbing ions were outside the "atomic sphere", i.e., outside the bound electron cloud of the radiating atom/ion (non-penetrating configurations).In the present paper, we take into the contribution to the spectral line asymmetry from penetrating configurations, i.e., from configurations where the perturbing ion is inside the bound electron cloud of the radiating atom/ion (hereafter, radiator).We show that, in high-density plasmas, the allowance for penetrating ions can result in significant corrections to the electron density deduced from the spectral line asymmetry.

Allowance for Penetrating Ions
Let us first present a brief overview of the underlying theory for non-penetrating configurations.The dipole interaction of the radiator with perturbing ions outside the bound electron cloud, being calculated in the first order of the perturbation theory, splits the spectral line into Stark components symmetrically with respect to the unperturbed frequency or wavelength-in terms of both positions and intensities of the Stark components.The quadrupole interactions of the radiator with perturbing ions outside the bound electron cloud, being calculated in its first nonvanishing order, causes the asymmetry of the Stark splitting-in terms of both positions and intensities of the Stark components.The latter is the primary source of the asymmetry: other sources of the asymmetry-such as, but not limited to, e.g., the dipole interaction in the second order (known as the quadratic Stark effect), the quadrupole interaction in the second order, the octupole interaction in the first order-add to the asymmetry only higher-order corrections in terms of the corresponding small parameter n 2 a 0 N e 1/3 /Z 4/3 , where n is the principal quantum number of the upper level involved in the radiative transition, a 0 is the Bohr radius, Z is the charge of plasma ions and the nuclear charge of the radiating ion.We note that Ref.
[40] from chp. 1 of [10] on the paper referred to here as [8] has typographic errors.The correct one is our Ref.[8] here. 2 The Boltzmann factor exp(−h ∆ω/T) also contributes to the asymmetry (here ∆ω is the detuning from the unperturbed frequency of the spectral line and T is the temperature).For quasistatic wings, h ∆ω/T scales with the electron density as (a 0 N e 1/3 ) 2 .Therefore, for plasmas of the electron densities N e << 6.7 − 10 24 cm −3 (the right side being the atomic unit of the electron density), the Boltzmann factor contribution to the asymmetry is much smaller than the quadrupole interaction contribution to the asymmetry that scales as a 0 N e 1/3 .Additionally, there is also the factor (1 + ω/ω 0 ) 4 caused by the scaling of the dipole radiation intensity (ω 0 being the unperturbed frequency).The asymmetry contributions of this factor and of the Boltzmann factor essentially cancel each other out (see, e.g., Section 5.11 of paper [5]).(Continued at the bottom of the next page).There is also so-called trivial contribution to the asymmetry caused by the conversion from the frequency scale to the wavelength scale.This consists of two factors (see, e.g., Section 5.8 of paper [5]): the transformation of the argument ∆ω (given by Equation (24) from [5]) and the transformation of the intensity (given by Equation (25) from [5]).These two factors essentially cancel each other out (as shown in [5]), so that the resulting trivial contribution to the asymmetry is much smaller than the quadrupole interaction contribution to the asymmetry.However, in paper [11], it was shown that the quadrupole interaction, despite causing the asymmetric splitting of the spectral line into Stark components, does not shift the center of gravity of the line profile.Therefore, in the new diagnostic method presented in paper [8], first the center of gravity of the experimental profile was determined, and then it was taken as the reference point.Then, with respect to this point, the integrated intensities of the blue (I B ) and red (I R ) wings of the experimental profile were found.After that, the experimental degree of asymmetry, defined as was determined and then compared with the corresponding theoretical value given below.
The theoretical intensities of the blue and red wings, resulting from dipole and quadrupole interactions of the radiator with perturbing ions outside the bound electron cloud, can be expressed as follows (see paper [8]): and where Z p is the charge of perturbing ions, Z r is the nuclear charge of the radiator, a o is the Bohr radius, and R o = [(4π/3)N p ] −1/3 is the mean interionic distance, N p = N e /Z p being the perturbing ion density.Here, I k are the unperturbed intensity and the quadrupole correction to the intensity, respectively, the subscript k being the label of Stark components of the spectral line; k > 0 and k < 0 correspond to the blue-shifted and red-shifted components, respectively (the values of k for several Lyman and Balmer lines were tabulated in paper [2]).The quantity <R 0 /R> is the scaled inverse distance between the perturbing ion and the radiator averaged over the distribution of such distances.
Finally, the theoretical degree of asymmetry was presented in paper [8] in the form: Then the electron density N e was determined in paper [8] by substituting the experimental degree of asymmetry into the left side of Equation ( 4).
In the present paper, we consider the contribution of penetrating ions to the spectral line asymmetry in order to refine this diagnostic method.For simplicity, we limit ourselves below to the practically important case Z p = Z r = Z.The energy shifts due to penetrating ions can be calculated by the perturbation theory on the basis of the spherical wave functions of the so-called "united atom" of the nuclear charge 2Z.
The perturbed energy shifts (counted from the unperturbed energies) for the orbital quantum number l > 0 are given by (see, e.g., Equations ( 6) and (7) from paper [12] or Equations (5.11) and (5.12) from book [13]): For the case of l = 0, the calculated energy shift is: We note that Equation ( 6) can also be obtained from Equation ( 5), first by setting m = 0, and then by canceling out l(l + 1) in the numerator and denominator, and by setting l = 0. (This was mentioned in book [12], but in Equation (5.11) from [12] corresponding to our Equation ( 6), there was a typographic error in the sign.) The frequency change of an individual Stark component is thus given by where For the specific case where either l = 0 or l = 0, Equation ( 8) reduces to Then, the quasi-static profile of each Stark component can be represented in the form: Here, u ≡ R 2 , and the probability of finding the perturbing ion a distance u away from the radiating atom is taken to be the binary distribution.For simplifying the integration, we use the expansion of the distribution in powers u/R 0 2 and keep the terms up to ~u2 : For the case of a hydrogenic radiator under the presence of a penetrating ion, the relative intensities of each line component can be best calculated analytically using the robust perturbation theory developed by Oks and Uzer [14].A more detailed explanation of this procedure is outlined in Appendix A. The relative intensities of each component can be written as where ∆ 0 Ik and ∆ 1 Ik are tabulated in Appendix B for each component of the spectral line Balmer-alpha, considered here as an example.These coefficients represent corrections to the intensity of the line, calculated from the perturbation theory briefly mentioned above.
The upper limit u max of the integration in Equation ( 10) should be the smallest of the following two "candidates".One candidate for u max is the root mean square size of the bound electron cloud, which depends on the sublevel in consideration: The other candidate for u max is defined by the limit of the applicability of the perturbation theory.Of course, this would ensure that formally calculated corrections to the energy and intensities of the spectral line would remain relatively small.The allowance for penetrating ions shifts the center of gravity of the spectral line, as shown in paper [15].(This is the only contribution to the shift of the center of gravity, since the dipole and quadrupole interactions of the radiator with perturbing ions outside the bound electron cloud do not shift the center of gravity, as shown in paper [11] and mentioned above).For the He II Balmer-alpha line, which we use as an example, the center of gravity shift due to penetrating ions was calculated analytically in paper [15] to be So, with the allowance for penetrating ions, the reference point for calculating the integrated intensities of the blue and red wings must be shifted by the amount given by Equation (14).
After carrying out the integration in Equation ( 10), the profile reduces to where Θ[. ..] is the Heaviside step function and u 0 is the root of the delta function, given by Thus, for the contributions of the penetrating ions to the integrated intensities of the blue and read parts of the line profile, we get and Here which is obtained by equating to zero the argument of the Heaviside step function.Additionally, what is meant in Equations ( 17) and ( 18) by k < 0 (or k > 0) is the inclusion of only those components which involve corrections to the energy that are positive (or negative), implying blue-shifted (or red-shifted) components of the spectral line.By combining the above result with the contribution of the quadrupole interaction (the interaction of the radiator with perturbing ions outside the bound electron cloud) to the integrated intensities of the blue and read parts of the profile, we obtain our final result for the degree of asymmetry where subscript act stands for actual-in distinction to ρ quad .The combination of Equations ( 4) and ( 20) connects the degree of asymmetry with the electron density N e and thus allows a more accurate determination of the electron density from the experimental asymmetry.We illustrate this below by the example of the He II Balmer-alpha line.
Table 1 presents the following quantities for the He II Balmer-alpha line at five different values of the actual electron density: the theoretical degree of asymmetry ρ act calculated with the allowance for penetrating ions, -the theoretical degree of asymmetry ρ quad calculated without the allowance for penetrating ions, -the electron density N e,quad that would be deduced from the experimental asymmetry degree while disregarding the contribution of the penetrating ions, -the relative error |N e,quad -N e,act |/N e,act in determining the electron density from the experimental asymmetry degree while disregarding the contribution of the penetrating ions.
Table 1.The relative error in determining the electron density N e from the experimental asymmetry degree while disregarding the contribution of the penetrating ions for the He II Balmer-alpha line.The physical quantities in Table 1 are explained in the text directly above Table 1.
N e,act /(10 It is seen that in high-density plasmas, the allowance for penetrating ions can indeed result in significant corrections to the electron density deduced from the spectral line asymmetry.

Conclusions
To improve the diagnostic method for measuring the electron density using the asymmetry of spectral lines in dense plasmas, we took into account the contribution to the spectral line asymmetry from penetrating configurations, i.e., from the configurations where the perturbing ion is inside the bound electron cloud of the radiating atom/ion.After performing the corresponding analytical calculations, we demonstrated that in high-density plasmas, the allowance for penetrating ions can result in significant corrections to the electron density deduced from the spectral line asymmetry.
It is worth clarifying why we took into account the shift of the line as a whole due to penetrating ions, but did not take into account other mechanisms shifting the line as a whole, such as, e.g., plasma polarization shift and the shift by plasma electrons.The experimental integrated intensities of the blue (I B ) and red (I R ) parts of the profile are calculated with respect to the experimental center of gravity of the profile.The latter shifts of the line as a whole do not contribute to the asymmetry, and thus should not affect the experimental values of I B and I R .The reason why we took into account the shift of the line as a whole by penetrating ions is that penetrating ions contribute simultaneously to both the asymmetry and the shift of the line as the whole.Since these two effects of penetrating ions are two sides of the same coin, both of them should be taken into account.
We mention in passing that the potential transition of electrons into the quasistatic regime is practically irrelevant to the asymmetry.Indeed, as Demura and Sholin wrote in paper [3], for the quadrupole interaction U ~Q/R 3 , which controls the asymmetry of hydrogenic spectral lines, electrons can become quasistatic at the frequency detuning from the line center ∆ω ~vTe 3/2 /Q 1/2 , where v Te is the mean thermal velocity of plasma electrons-according to Holstein [16] (see also Sobelman book [17]).However, such detuning significantly exceeds the mean separation between spectral lines-both for the Lyman and Balmer series, as noted by Demura and Sholin [3], thus making the potential transition of electron into the quasistatic regime irrelevant to the problem of asymmetry of hydrogenic spectral lines.A similar conclusion was drawn also in paper [18].Finally, we note that the electron densities N e ~(10 18 -10 19 ) cm −3 , which we used in the illustrative example of the He II Balmer-alpha line, are achievable in plasma spectroscopy.Examples include experiment [19], with a hydrogen plasma, and experiment [20], with a helium plasma.