1. Introduction
Hydrogen spectral lines are of permanent interest both from the fundamental point of view and for applications in plasma diagnostics. Since the middle of the 20th century, a significant number of works and monographs (e.g., [
1,
2,
3,
4,
5,
6,
7,
8]) have been accumulated on the spectra of hydrogen plasma. However, the effect of a magnetic field on spectral lines remains a problem in plasma diagnostics. Difficulties connected with the analysis of the Stark–Zeeman spectra relate to a hydrogen in external crossed
electric and
magnetic fields.
Calculation of spectral line shape in plasmas is complicated by two serious problems. The first is the complex structure of the dipole matrix elements, as well as the rapid growth of the array of radiative transitions. The second is connected with the problem of the influence of an ion thermal motion on the intensity profile. This paper shows how one can get around these difficulties. The complex expressions for the transition probabilities can be simplified by using the semiclassical approximation for coordinate matrix elements obtained by S.A. Gulyaev [
9,
10] and specific properties of the Wigner d-functions [
11]. The thermal motion of ions can be taken into account using the FFM. It turns out that spectral line shape in a plasma with moving ions is a functional of the static profile [
12].
The first attempt to provide a solution to this problem was made within the framework of the classical mechanics [
13]. The perturbed motion of the electron was reduced to an independent precession of two vectors representing the combination of the angular momentum and the averaged coordinate around two different axes. The quantum approach to this problem has become possible thanks to the fundamental work of V.A. Fock [
14], who showed that an electron in the Coulomb field has enhanced O(4) (instead of common O(3)) symmetry. This fact leads to the existence of two additional constant of motion
, which are connected with the orbital momentum
by the simple relation
A rigorous quantum consideration of this problem is given in [
15]. In this work, the authors used the O(4) symmetry properties of the Coulomb field to obtain the spectra of a hydrogen atom in the external electric and magnetic fields. The Hamiltonian of this system is equal to
where
,
and
are the momentum, the coordinate and the angular momentum operators of the electron, correspondingly, Z is the charge of nuclei. This formula and every other in this paper is written in the atomic units. The perturbed part
can be rewritten in another way.
where
where 1 relates to + and 2 to −.
is the specific constant of motion in the Coulomb field—the Runge–Lenz vector.
We can do this, because in the Coulomb field there is connection between the Runge–Lenz vector and the coordinate (this is valid only within the manifold of a fixed principal quantum number
n):
The energy shift is equal to
where
and
are projections of (
3) on the vectors (
4).
Obviously, the vectors (
3) are constants of motion in the Coulomb field. Moreover, they have the angular momentum properties. Coming to the parabolic basis, there is the connection between the parabolic quantum numbers
,
and projections on the single direction [
16]
where
are projections of (
3) on the z direction (quantization axis) and m is the magnetic quantum number.
Using the angular momentum properties of vectors (
3), we can express the wave functions in the representation of
in terms of the parabolic states.
where
is the Wigner d-function.
In (
8),
are the angles between vectors
and
. We choose the reference frame in which the direction of the magnetic field coincides with the z axis.
where
is the angle between electric
and magnetic fields
.
The coordinate matrix elements in basis (
8) has the following form
where
(the intensity of radiation in the dipole approximation is proportional to the squared absolute value of the coordinate matrix element).
Accurate quantum expressions for the matrix elements
in (
11) were obtained by Gordon [
17,
18] and have a very complicated structure. They contain the hyper-geometric series which makes calculations for Rydberg atoms very cumbersome. The detailed analysis of computational complexities and ways to get around them are presented in [
19]. The array of radiative transitions grows proportionally to
. To carry out calculations for highly excited levels, the author of [
9,
10] obtained the approximation of
and developed the method of the distribution of atomic transitions into special groups. These results and the usage of the specific d-functions properties allow one to simplify the complicated expression (
11). The results for the
(
) and
(
) series are obtained in [
20,
21]. The basics of the approach to calculating these matrix elements and the results (see Formulas (
A11)–(
A16)) are presented in the
Appendix A.
In the present paper, we take into account thermal velocity of ions. This problem is closely related to stochastic processes in the Coulomb-like interacting medium. Statistical aspects of the collective motion in the Coulomb field are considered in [
22]. The frequency fluctuation model (FFM) consists in the dependence of the spectral line profile on the jumping ion frequency
.
where
is the density and
is the thermal velocity of ions.
Using the FFM, the intensity
can be analytically expressed as the functional of the normalized static profile
[
12].
A similar calculation, using the FFM, but without taking into account the influence of the magnetic field on the shape of spectral lines, is given in [
23]. Moreover, this work contains the detailed analysis of the influence of the thermal ion velocity on the plasma line shapes.
The spectral lines profiles in the presence of a magnetic field were calculated by Novikov et al. [
24]. They used the accurate analytical expressions for the dipole matrix elements. However, this work contains profiles only of
and
lines. In addition, the authors did not consider the thermal motion of ions. Application of the frequency-fluctuation model to Stark–Zeeman line shapes was demonstrated by Ferri et al. [
25]. Again, this work considers only the transitions between levels with a low value of the principal quantum number n. Computer modeling, taking into account the largest possible number of effects on the shape of spectral lines, is presented in the works of Rosato et al. [
26,
27].
In the present paper, we derive the analytical expressions for spectral lines profiles in plasma. These formulas are convenient for simply performing calculations with them. In fact, we provide the algorithm of spectral lines shapes calculations with given parameters of plasma. The expressions for transitions with and are presented.
2. Description of the Method
Firstly, we consider the general expression for the plasma spectral line static profile.
where
is the full set of all quantum numbers related to the initial and final states,
is the polarization,
is the distribution function of the absolute value of electric field,
defines the angle between the magnetic field and an ion microfield, and
is the dipole matrix element (expressions (
A11)–(
A13) for
and (
A14)–(
A16) for
are presented in
Appendix A).
Coming up to our notation, one obtains
Here, the values with bar relate to the final state.
In the summation, there are only two terms. For example, if one calculates the radiation intensity with direction of observation parallel to magnetic field, corresponds to X and Y directions. In other words, in this example, one has to put and .
In (
15),
is the Dirac delta-function and
is the energy shift that corresponds to the set of quantum numbers
. In our case, we have
Here, we use the absolute values of the vectors (
4). The fact that the system has the circular symmetry is also used, which means that there is no dependence of the energy shift and the matrix elements on the azimuthal angle.
As
, one can use the Holtsmark distribution
where
is the charge and
is the density of ions. The detailed analysis of the Holtsmark distribution is presented in [
5,
12].
To take into account the thermal velocity of ions in plasmas, we use Formulas (
13) and (
14). One can substitute expression (
15) into (
14). After that, it is possible to integrate over
and get rid of the delta functions.
Now, it is necessary to take into account the Doppler and the electron broadening mechanisms. To do this, one has to calculate the convolution of
from (
13) (with
from (
19)) with the Voigt profile
Expression (
21) relates to the Doppler broadening [
5]. Here, D is the Doppler parameter
where
, c is the speed of light in vacuum, and M and
are the mass and temperature of atoms, respectively.
That Lorentz distribution corresponds to the electron broadening. Generally, the calculation of the parameter
is a complicated process. We use the simplified approach to the electron broadening (see [
5]).
where
and
are the density and thermal velocity of electrons, correspondingly, and
is the Debye radius in plasmas.
where
is the coordinate operator,
is the transition intensity, and a and b denote different states referring to levels n,
. For example,
More accurate calculations for the electron broadening are discussed in [
28]. Moreover, modification of this broadening theory is presented in the Appendix of [
24].
Finally, we can obtain the formula for the intensity of radiation as the function of the energy shift (frequency). It is the convolution of expressions (
13) and (
20).
To sum up, we obtain the algorithm for plasma line shapes calculations. Expressions for the dipole matrix elements are presented in the
Appendix A. We take into account Stark–Zeeman, Doppler, and electron mechanisms of broadening. Moreover, using the FFM, the effect of thermal velocity is considered. For every step of the calculation, we use analytical expressions. The approximations for the transitions intensities are derived for lines with
and
. The main advantage this method is the opportunity of performing calculations for any large principle numbers.
3. Results and Discussion
Specific calculations for a deuterium plasma are demonstrated in this section. To estimate the accuracy of the method, we compare the results of other groups for the and spectral lines with our algorithm.
Comparison of semiclassical (universal approach) calculations and the results of Ferri et al. [
25] are presented in
Figure 1. One can observe the satisfactory correspondence between two approaches. Even for low levels
and
, the Rydberg approximation shows a high degree of accuracy. However, there is a slight discrepancy, mainly due to the imprecision in the calculations of the Stark shift (here,
). This calculation is extremely representative because in both approaches the thermal velocity of ions is taken into account in the same way. In fact, the main difference here is the choice of expressions for the transition probabilities.
One can observe a tolerable coincidence between the Rydberg approach to spectral lines in plasmas and the computer modeling in [
27]. Calculations for
line are presented in
Figure 2. The semiclassical approach has a distinct feature: the narrowing in the center. It is connected with the widening on the sides. The reason for these little “wings” is Zeeman components. The inaccuracy of the Rydberg approximation leads to a slight decrease in the Stark shift. However, this result should be considered quite satisfactory because distance between atomic levels is equal to principle quantum number of a lower state. Generally speaking, the condition for the applicability of our approach is
, but it turns out that it works well even for the Balmer series.
Obviously, with growth of n, the Rydberg approximation will be practically indistinguishable from the exact result. For large values of the principle quantum number, distance between atomic levels
, thus it will not influence the Stark shift, as well as the electron broadening parameter
(except specific cases of the interference of contributions to the electron width, the details can be found in [
4]). For highly excited levels, the electron impact on the broadening can blur spectral lines. Thus, the small inexactness connected with the Rydberg approximation can all the more be neglected.
To show how the shape of the spectral line changes,
Figure 3 is presented. Using the expression (
27), spectral profiles are obtained for the Paschen series (the transition 4-3). It is possible to trace how the Zeeman components disappear with increasing of an ion’s temperature due to the Doppler effect and an ion thermal motion. In the first graph, two Zeeman “wings” are clearly visible. We want to underline that, in this calculation, the temperatures and concentrations of ions and electrons are different. In the fourth graph, the Zeeman components are practically blurred.
In this method, we do not consider two effects connected with external magnetic field. The first one is that the environment can be spatially inhomogeneous. Therefore, this method is applicable only for comparison with local measurements, so that the temperature and density gradients cannot be large. The second effect is connected with the influence of an electron’s helical path related to the motion in an external magnetic field. Because of this, the argument of the logarithm in Formula (
24) can be greatly changed. The influence of a change in the trajectories of an electron on the shape of spectral lines is described in detail in [
29]. In the present paper, we focus on the problem of the great array of radiation transition. Our purpose is to show that neglecting of the big part of transitions does not strongly affect the shape of the spectral line. This can be seen by looking at
Figure 1 and
Figure 2. Comparison of our method with calculations of other groups clearly demonstrates that a significant part of radiative transitions does not contribute to the formation of the spectral profile. Because of that, there is a significant opportunity to simplify calculations of the spectral line shape in plasmas. This method is mainly applicable for the case when ion’s Stark broadening dominates over the electron one, in both static and impact limits for ions.
In this paper, a whole complex of problems in the theory of broadening of spectral lines in plasma is considered. However, individually, many of these topics have already been discussed previously. Semiclassical approximation for the transition probability (the intensity) is presented in [
30]. Mathematical aspects of the application of the O(4) symmetry to the theory of spectral line broadening are discussed in the work [
31]. However, this approach is not convenient for describing the transition array in the case of a magnetized plasma. The applicability of the FFM to plasma in the presence of a magnetic field is discussed in [
32].
The presented algorithm can be a convenient tool for diagnosing helium plasma in the ITER divertor [
33,
34]. An example of calculation for helium is shown in
Figure 4 (transition 5-4).
4. Conclusions
A semiclassical approach to the spectral lines in plasmas is presented. We demonstrate how one can use the analytical expressions for calculating
and
line shapes. Thanks to the presented method, it is possible to calculate the intensity profiles for transitions with large principle quantum numbers. Basically, this approach is suitable for the spectral line shapes of hydrogen-like ions. For calculations within the visible range, it is necessary to consider transitions with larger n. Moreover, it is shown that this approach can be used to achieve satisfactory agreement with the transitions related to the Balmer series (
Figure 1 and
Figure 2).
Using the FFM, one can relatively simply take into account the effect of the thermal motion of ions in the plasma on the shape of spectral lines. The FFM profile depends on the functions
(
14). The integrand in
depends on static profile
(
15). The function
is expressed in terms of delta-functions. It allows one do a simple integration and obtain Formula (
19) for
.
In [
12], the authors considered a Hydrogen atom in external crossed electric
and magnetic
fields. They performed calculations in the specific basis (
8). These states are closely related to the parabolic quantization on two different axes. Calculation of the transition probabilities in this basis leads one to Formula (
11). Using the symmetry properties of the Coulomb field, the Wigner d-functions recurrence relations, and the Rydberg asymptotic formulas for the coordinate matrix elements [
9,
10], the simple semiclassical approximations for the dipole matrix elements in representation of states (
8) are obtained in [
20,
21] and applied to the spectral line shape calculations. The approximations for the transitions probabilities in this basis are presented in
Appendix A (see (
A11)–(
A16)).
To demonstrate the power of the Rydberg atom approach, Balmer series line shapes are calculated and compared with the results of other groups. The shape of
line is presented in
Figure 1. One can observe a good correspondence with accurate consideration from [
25]. The comparison of Rydberg approach and the computer simulation [
27] is presented in
Figure 2. Even for the 4-2 transition, the analytical approximation shows a good correspondence with quantum calculations.
The specific calculation for the
series is demonstrated in
Figure 3. Using the presented algorithm, the shapes of spectral lines were calculated for various plasma parameters. This figure clearly illustrates how different mechanisms affect the shape of the spectral line.