*Article* **Polarization of Lyman-***α* **Line Due to the Anisotropy of Electron Collisions in a Plasma**

**Motoshi Goto 1,2,\* and Nilam Ramaiya <sup>3</sup>**


**\*** Correspondence: goto.motoshi@nifs.ac.jp

**Abstract:** We have developed an atomic model for calculating the polarization state of the Lyman-*α* line in plasma caused by anisotropic electron collision excitations. The model assumes a nonequilibrium state of the electron temperature between the directions parallel (*T*) and perpendicular (*T*⊥) to the magnetic field. A simplified assumption on the formation of an excited state population in the model is justified by detailed analysis of population flows regarding the upper state of the Lyman-*α* transition with the help of collisional-radiative model calculations. Calculation results give the polarization degree of several percent under typical conditions in the edge region of a magnetically confined fusion plasma. It is also found that the relaxation of polarization due to collisional averaging among the magnetic sublevels is effective in the electron density region considered. An analysis of the experimental data measured in the Large Helical Device gives *<sup>T</sup>*⊥/*T* <sup>=</sup> 7.6 at the expected Lyman-*α* emission location outside the confined region. The result is derived with the absolute polarization degree of 0.033, and *<sup>T</sup>*<sup>⊥</sup> <sup>=</sup> 32 eV and *<sup>n</sup>*<sup>e</sup> <sup>=</sup> 9.6 <sup>×</sup> <sup>10</sup><sup>18</sup> <sup>m</sup>−<sup>3</sup> measured by the Thomson scattering diagnostic system.

**Keywords:** plasma spectroscopy; polarization; Lyman-alpha; nuclear fusion

#### **1. Introduction**

In a magnetically confined fusion plasma, the velocity distribution function (VDF) of electrons and ions is thought to be more or less anisotropic. For example, energetic ions are unidirectionally introduced by the neutral beam, and the cyclotron motions of ions and electrons are selectively accelerated by the electron cyclotron resonance heating (ECRH) and the ion cyclotron resonance heating (ICRH), respectively. Because of the unavoidable magnetic field ripple, confinement characteristics are different between the particles having a large pitch angle and a small pitch angle with respect to the magnetic field. The former and the latter are called trapped particles and passing particles, respectively, and are sometimes treated separately when the particle transport is considered.

An example of the research relating to this topic is the influence of plasma pressure anisotropy on the MHD (Magnetohydrodynamic) equilibria, which has been intensively investigated [1]. The ITER experiment would also be influenced by the problem of anisotropy [2]. In the Large Helical Device (LHD), the so-called density clamping observed with strong ECRH is thought to be attributed to the difference in the confinement characteristics of the trapped and passing electrons [3]. The anisotropy also plays a role in the plasma edge region. The radial electric field formation in the edge stochastic region is thought to be related to the anisotropy in the electron VDF (EVDF) [4].

Although the anisotropy is regarded as an important subject for characterizing the plasma confinement as seen above, no reliable measurement method for the anisotropic VDF of electrons and ions has been established to date. Under such circumstances, polarization spectroscopy has been proposed as a technique to address the problem of anisotropy

**Citation:** Goto, M.; Ramaiya, N. Polarization of Lyman-*α* Line Due to the Anisotropy of Electron Collisions in a Plasma. *Symmetry* **2021**, *13*, 297. https://doi.org/10.3390/ sym13020297

Academic Editor: Eugene Oks Received: 22 January 2021 Accepted: 5 February 2021 Published: 9 February 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

in the EVDF [5]. This novel diagnostic method consists of two issues, i.e., the measurement of polarization in line emission from the plasma and the construction of an atomic model for analyzing the observation data. A critical problem in the measurement is a difficulty in detecting the polarization of line emission with accuracy, the degree of which is estimated in the order of one percent or smaller. Although some measurements have been attempted, no reliable results have been obtained [6,7].

Recently, the CLASP (Chromospheric Lyman-Alpha SpectroPolarimetry) project led by the NAOJ (National Astronomical Observatory of Japan) has successfully measured the polarization state of the hydrogen Lyman-*α* line in the solar atmosphere [8]. Although the polarization formation mechanism in the solar atmosphere is photoexcitation by the anisotropic radiation field, which is different from the mechanism in the fusion plasma, the observation technique itself can be transferred to the measurement for a fusion plasma. Actually, the same technique as CLASP has been attempted in LHD, and the polarization of the Lyman-*α* line has been successfully detected [9].

As for the atomic model, a sophisticated framework has been developed by Fujimoto [5], and some actual applications have been made [6,10]. We here report details of the implementation of Fujimoto's framework for the Lyman-*α* line with goal of utilizing it for the analysis of the polarization data taken in the LHD experiment. Furthermore, an analysis has been made of the LHD experimental data and derivation of an anisotropic EVDF is attempted.

#### **2. Theoretical Model for the Line Emission Polarization**

#### *2.1. Polarization Formation of the Lyman-α Line*

The polarization of an emission line originates in a population imbalance between the magnetic sublevels in the upper state of the transition. The Lyman-*α* line consists of two fine structure lines, i.e., 1 2S1/2 – 2 2P1/2 and 1 2S1/2 – 2 2P3/2. Figure 1 shows all transitions between magnetic sublevels composing the Lyman-*α* line. The Δ*mJ* = 0 and Δ*mJ* = ±1 transitions emit light linearly polarized in the quantization axis direction (*π* light) and circularly polarized on the plane perpendicular to the quantization axis (*σ* light), respectively, where *mJ* is the magnetic quantum number. The numbers next to the lines indicating the transitions are the relative values of the Einstein A coefficients. It is confirmed that if all the magnetic sublevels have the same population, the intensities of *π*, *σ*+, and *σ*<sup>−</sup> lights are identical, i.e., there is no polarization, and otherwise the line is polarized.

**Figure 1.** Line components included in the Lyman-*α* line. The solid and dashed lines represent the *π*- and *σ*-light, respectively. The numbers next to the lines indicate relative values of the Einstein A coefficient.

We assume axisymmetry with respect to the quantization axis which is taken in the magnetic field direction later. In this case, the population distribution has a "mirror symmetry", i.e., the populations of *mJ* and −*mJ* sublevels are identical in each of the 2 2P1/2 and 2 2P3/2 states. Because of this restriction, the line corresponding to the 1 2S1/2 – 2 2P1/2 transition is never polarized because *π*, *σ*+, and *σ*<sup>−</sup> light intensities are always identical. As for the 1 2S1/2 – 2 2P3/2 transition, populations of the *mJ* = |1/2| and *mJ* = |3/2| substates in the 2 2P3/2 state can be different, which could give rise to the line polarization because only the *mJ* = |1/2| substates are responsible for the *π* light. We first focus on deriving the polarization state of the 1 2S1/2 – 2 2P3/2 line and then incorporate the influence of the unpolarized 1 2S1/2 – 2 2P1/2 line into the result as explained below to enable a direct comparison of the model results with the observation results.

When the measurement is made from the direction perpendicular to the quantization axis with a linear polarizer, the polarization degree *P* of an emission line is generally defined as

*<sup>P</sup>* <sup>=</sup> *<sup>I</sup><sup>π</sup>* <sup>−</sup> *<sup>I</sup><sup>σ</sup> I<sup>π</sup>* + *I<sup>σ</sup>* , (1)

where *Iπ* and *Iσ* represent the line intensities observed when the polarizer is directed in the direction parallel and perpendicular to the quantization axis, respectively. Because the Lyman-*α* line includes the two fine structure components, *P* can be explicitly written as

$$P = \frac{\left(I\_{\pi}(3/2) + I\_{\pi}(1/2)\right) - \left(I\_{\sigma}(3/2) + I\_{\sigma}(1/2)\right)}{\left(I\_{\pi}(3/2) + I\_{\pi}(1/2)\right) + \left(I\_{\sigma}(3/2) + I\_{\sigma}(1/2)\right)'} \tag{2}$$

where *Iπ*(1/2) and *Iσ*(1/2) are the intensities of the *π*- and *σ*-components of the 1 2S1/2 – 2 2P1/2 line, respectively, and *Iπ*(3/2) and *Iσ*(3/2) are the same but of the 1 2S1/2 – 2 2P3/2 line, respectively. Because the 1 2S1/2 – 2 2P1/2 line is unpolarized under the present condition, the relation

$$I\_{\pi}(1/2) - I\_{\sigma}(1/2) = 0\tag{3}$$

should always hold. On the other hand, we assume that the population ratio of the 2 2P1/2 state to the 2 2P3/2 state follows the ratio of their statistical weights, i.e., the former is half the latter. In that case, the same is true for the line intensities of the 1 2S1/2 – 2 2P1/2 and 1 2S1/2 – 2 2P3/2 transitions, i.e.,

$$I\_{\pi}(1/2) + I\_{\sigma}(1/2) = \frac{1}{2} [I\_{\pi}(3/2) + I\_{\sigma}(3/2)].\tag{4}$$

By using the relations of Equations (3) and (4), we can rewrite Equation (2) as

$$\begin{split}P &= \frac{\left(I\_{\pi}(3/2) - I\_{\sigma}(3/2)\right) + \left(I\_{\pi}(1/2) - I\_{\sigma}(1/2)\right)}{\left(I\_{\pi}(3/2) + I\_{\sigma}(3/2)\right) + \left(I\_{\pi}(1/2) + I\_{\sigma}(1/2)\right)}\\ &= \frac{I\_{\pi}(3/2) - I\_{\sigma}(3/2)}{\frac{3}{2}\left[I\_{\pi}(3/2) + I\_{\sigma}(3/2)\right]}\\ &= \frac{2}{3}P(3/2),\end{split} \tag{5}$$

where *P*(3/2) is the polarization degree of the 1 2S1/2 – 2 2P3/2 line.

#### *2.2. Polarization Due to Anisotropic Electron Collisions*

The condition of an excited state which has a population imbalance among the magnetic sublevels is well represented by the density matrix [11]. Under an axisymmetric system, the spherical coordinate representation of the density matrix *ρ* of the state 2 2P3/2 can be expanded [5] as

$$
\rho(p) = \rho\_0^0(p) T\_0^{(0)}(p) + \rho\_0^2(p) T\_0^{(2)}(p),\tag{6}
$$

where *<sup>p</sup>* stands for the state 2 2P3/2 and *<sup>T</sup>*(*k*) *<sup>q</sup>* (*p*) is the so-called irreducible tensor operator. The coefficients *ρ*<sup>0</sup> <sup>0</sup>(*p*) and *<sup>ρ</sup>*<sup>2</sup> <sup>0</sup>(*p*) respectively correspond to the population and the alignment, the latter of which expresses the inhomogeneity over the magnetic sublevels in the state *p*. We hereafter use *a*(*p*) instead of *ρ*<sup>2</sup> <sup>0</sup>(*p*) for simplicity. The conventional population

is given as *<sup>n</sup>*(*p*) = <sup>√</sup>2*<sup>J</sup>* <sup>+</sup> <sup>1</sup>*ρ*<sup>0</sup> <sup>0</sup>(*p*), where *J* is the total angular momentum quantum number of the state *p*.

The population imbalance among the magnetic sublevels could be created by anisotropic electron collisions. We consider a simple atomic model for a quantitative calculation of *n*(*p*) and *a*(*p*) where the population inflow and outflow concerning the state *p* are balanced by the electron impact excitation of the ground state atoms and the spontaneous radiative decay. Such a condition can be expressed [5] as

$$C^{0,0}(1,p)n\_\text{e}n(1) = \sum\_s A(p,s)n(p),\tag{7}$$

where *C*0,0(1, *p*) is the rate coefficient of the electron impact excitation from the ground state denoted as "1" to the state *p*, *A*(*p*,*s*) is the Einstein A coefficient of the transition from *p* to a state *s* placed energetically lower than *p*, and *n*e is the electron density. The validity of this model under the conditions assumed here will be examined later.

The equilibrium condition of *a*(*p*) is similarly expressed [5] as

$$\mathcal{L}^{0,2}(1,p)n\_{\mathfrak{e}}n(1) = \left[\sum\_{s} A(p,s) + \mathcal{C}^{2,2}(p,p)n\_{\mathfrak{e}}\right]a(p),\tag{8}$$

where *C*0,2(1, *p*) is the alignment creation rate coefficient accompanying the excitation from the ground state to the state *p* and *C*2,2(*p*, *p*) is the alignment destruction rate coefficient in the state *p*. The population *n*(*p*) and the alignment *a*(*p*) are then expressed as

$$n(p) = \frac{\mathbb{C}^{0,0}(1, p)n\_{\text{c}}}{\sum\_{\text{s}} A(p, s)} n(1),\tag{9}$$

$$a(p) = \frac{\mathbb{C}^{0,2}(1,p)n\_{\text{e}}}{\sum\_{s} A(p,s) + \mathbb{C}^{2,2}(p,p)n\_{\text{e}}} n(1). \tag{10}$$

As discussed in Section 2.1, the measurement gives the polarization degree *P*. We here introduce another quantity, "longitudinal alignment", *A*L, which is defined slightly differently from *P* [5] as

$$\begin{split} A\_{\rm L} &= \frac{I\_{\rm \pi} - I\_{\rm \sigma}}{I\_{\rm \pi} + 2I\_{\rm \sigma}} \\ &= \frac{2P}{3 - P} . \end{split} \tag{11}$$

The longitudinal alignment for the transition from the state *p* to *s*, i.e., *A*L(*p*,*s*), is directly related to *a*(*p*)/*n*(*p*), the normalized alignment, [5] as

$$A\_{\mathcal{L}}(p,s) = (-1)^{l\_p + l\_s} \sqrt{\frac{3}{2}} (2l\_p + 1) \left\{ \begin{array}{ccc} l\_p & l\_p & 2 \\ 1 & 1 & l\_s \end{array} \right\} \frac{a(p)}{n(p)},\tag{12}$$

where {· · · } is the 6-*j* symbol, and *Jp* and *Js* are the total angular momentum quantum numbers of the states *p* and *s*, respectively.

Calculations of the coefficients *C*0,0(1, *p*) and *C*0,2(1, *p*) can be carried out under a certain EVDF. We assume that the EVDF is axisymmetric with respect to the quantization axis and is expressed by two temperatures, *T* and *T*⊥, which are in the directions parallel and perpendicular, respectively, to the quantization axis (or the magnetic field). Such EVDFs are explicitly given [5] as

$$f(v, \theta) = \left(\frac{m}{2\pi k}\right)^{3/2} \left(\frac{1}{T\_{\perp}^{2}T\_{\parallel}}\right)^{1/2} \exp\left[-\frac{mv^{2}}{2k}\left(\frac{\sin^{2}\theta}{T\_{\perp}} + \frac{\cos^{2}\theta}{T\_{\parallel}}\right)\right],\tag{13}$$

where *v* is the absolute velocity, *θ* is the pitch angle of the velocity with respect to the quantization axis, and *m* and *k* are the electron mass and the Boltzmann constant, respectively. The function *f*(*v*, *θ*) is here normalized as

$$2\pi \int \int f(v,\theta)v^2 \sin\theta \text{d}v \text{d}\theta = 1.\tag{14}$$

Figure <sup>2</sup> shows examples of *<sup>f</sup>*(*v*, *<sup>θ</sup>*) in the case of *<sup>T</sup>*<sup>⊥</sup> = 30 eV and *<sup>T</sup>* = 10 eV (a) and of *<sup>T</sup>*<sup>⊥</sup> = 30 eV and *<sup>T</sup>* = 100 eV (b). It is noted that when we focus our interest on a group of electrons having any constant absolute velocity, the number of electrons is larger in the higher temperature direction than in the lower temperature direction.

**Figure 2.** Examples of the electron velocity distribution function *<sup>f</sup>*(*v*, *<sup>θ</sup>*) for (**a**) *<sup>T</sup>*<sup>⊥</sup> = 30 eV and *<sup>T</sup>* <sup>=</sup> 10 eV, and (**b**) *<sup>T</sup>*<sup>⊥</sup> <sup>=</sup> 30 eV and *<sup>T</sup>* <sup>=</sup> 100 eV cases.

The rate coefficients *C*0,0(1, *p*) and *C*0,2(1, *p*) are calculated [5] as

$$\mathcal{L}^{0,0}(1,p) = \int Q\_0^{0,0}(1,p) 4\pi f\_0(v) v^3 dv \tag{15}$$

and

$$\mathbb{C}^{0,2}(1,p) = \int \mathbb{Q}\_0^{0,2}(1,p)[4\pi f\_2(v)/5]v^3 \mathrm{d}v,\tag{16}$$

respectively, where *Q*0,0 <sup>0</sup> (1, *<sup>p</sup>*) and *<sup>Q</sup>*0,2 <sup>0</sup> (1, *p*) are the excitation and alignment creation cross sections, respectively, for the excitation from the ground state to the state *p*, and *f*0(*v*) and *f*2(*v*) are the coefficients of the expansion of *f*(*v*, *θ*) by the Legendre polynomials *PK*(cos *θ*) as

$$f(\upsilon, \theta) = \sum\_{K} f\_K(\upsilon) P\_K(\cos \theta). \tag{17}$$

The coefficient *fK*(*v*) is explicitly given as

$$f\_K(v) = \frac{2K+1}{2} \int f(v, \theta) P\_K(\cos \theta) \sin \theta \,\mathrm{d}\theta. \tag{18}$$

The alignment creation cross section *Q*0,2 <sup>0</sup> (1, *<sup>p</sup>*) is derived from *<sup>Q</sup>*0,0 <sup>0</sup> (1, *p*) as [5].

$$Q\_0^{0,2}(1,p) = (-1)^{l\_p + l\_s} \sqrt{\frac{2}{3}} (2l\_p + 1)^{-1} \left\{ \begin{array}{ccc} l\_p & l\_p & 2 \\ 1 & 1 & l\_s \end{array} \right\}^{-1} A\_\mathcal{L}(p, 1) Q\_0^{0,0}(1, p), \tag{19}$$

with *A*<sup>L</sup> for the case when the excitation takes place with mono-energetic beam collisions. We adopt the cross section data by Bray [12] for *Q*0,0 <sup>0</sup> and by James [13] for *P* which is translated into *A*<sup>L</sup> by Equation (11). It is noted that the data found in Refences [12,13] include both the 1 2S1/2 – 2 2P1/2 and 1 2S1/2 – 2 2P3/2 transitions. We here assume that 2/3 of the total cross section is for the 1 2S1/2 – 2 2P3/2 transition, and the total polarization degree multiplied by 3/2 is for the 1 2S1/2 – 2 2P3/2 transition (cf. Equation (5)). Figure 3 shows these elemental quantities relating to the 1 2S1/2 – 2 2P3/2 transition. The opposite polarity between *A*<sup>L</sup> and *Q*0,2 <sup>0</sup> is due to the 6-*j* symbol which is negative in the present case.

**Figure 3.** *<sup>A</sup>*<sup>L</sup> values under an assumption of a mono-energetic beam collision experiment (**a**) and *<sup>Q</sup>*0,0 0 and *Q*0,2 <sup>0</sup> (**b**) for the 1 2S1/2 – 2 2P3/2 transition. The actual *<sup>Q</sup>*0,2 <sup>0</sup> labeled with (−) takes negative values.

The alignment destruction process is understood as the relaxation of the population imbalance among the magnetic sublevels. It is known that this process due to electron collisions has some correlation with the Stark broadening of the emission line from that state and its rate coefficient can be approximated by the half width of the Stark broadening [14]. Here, the Stark broadening data for the Lyman-*α* line by Stehlé [15] are adopted for evaluating *C*2,2(*p*, *p*)*n*e. The results are plotted in Figure 4. It is confirmed that the alignment destruction rate increases almost linearly with *n*e.

**Figure 4.** Alignment destruction rate *C*2,2(*p*, *p*)*n*<sup>e</sup> evaluated from the Stark broadening width [15].

#### **3. Results and Discussion**

We have calculated *A*<sup>L</sup> for the Lyman-*α* line for typical plasma conditions in LHD as an example. The quantization axis is taken in the direction of the magnetic field. In our previous experimental study, we found that the linearly polarized light intensity takes on a maximum (minimum) value in the direction perpendicular (parallel) to the magnetic field [9]. Therefore, it is natural to regard the magnetic field direction as the symmetry axis of the system. The calculation is made with *T*<sup>⊥</sup> fixed at 32 eV, where the experimental data analyzed later are borne in mind, while *T* is scanned in a range around the fixed *T*⊥. We adopt several *n*<sup>e</sup> values which cover a typical *n*<sup>e</sup> range for the edge region of a magnetic fusion plasma. The polarization degree *P* for the Lyman-*α* line is derived from the *A*<sup>L</sup> values with Equation (11), and the results are plotted in Figure 5.

**Figure 5.** Example of the calculation results for *P* with several *n*<sup>e</sup> values. *T*<sup>⊥</sup> is fixed at 32 eV and *T* is scanned. The open circle represents the combination of *<sup>n</sup>*<sup>e</sup> <sup>=</sup> 9.6 <sup>×</sup> 1019 <sup>m</sup>−<sup>3</sup> and *<sup>P</sup>* <sup>=</sup> <sup>−</sup>0.033 corresponding to the experimental value in Ref. [9], from which *<sup>T</sup>* <sup>=</sup> 4.2 eV is derived.

It is confirmed that the line is unpolarized when *<sup>T</sup>* = *<sup>T</sup>*⊥, and the absolute polarization degree decreases with increasing *n*e, which is caused by collisional averaging over the magnetic sublevels. The positive *<sup>P</sup>* values for *<sup>T</sup>* > *<sup>T</sup>*<sup>⊥</sup> indicate that the *<sup>π</sup>*-light intensity is larger than that of *<sup>σ</sup>*-light, and negative *<sup>P</sup>* values for *<sup>T</sup>* < *<sup>T</sup>*<sup>⊥</sup> mean the opposite condition. These results are inferred from the tendency of the elemental *P* data in Figure 3 as follows.

In the present *T*<sup>e</sup> range at around 30 eV, the collision energy lower than 100 eV is dominant, and such collisions give rise to a positive *A*<sup>L</sup> or *P* as seen in Figure 3a. A positive *A*<sup>L</sup> or *P* means a higher intensity of linearly polarized light in the direction parallel with the electron beam axis than in the perpendicular direction (cf. Equation (1)). The condition of *<sup>T</sup>*<sup>⊥</sup> > *<sup>T</sup>* means more electrons in the perpendicular direction than in the direction parallel to the magnetic field as seen in Figure 2, which causes a higher intensity in the perpendicular direction than in the direction parallel to the magnetic field. This tendency finally results in a negative *<sup>A</sup>*<sup>L</sup> or *<sup>P</sup>*. The negative *<sup>P</sup>* values in the case of *<sup>T</sup>* < *<sup>T</sup>*<sup>⊥</sup> can be understood similarly.

We have recently reported a value of *P* = −0.033 as an example in the actual measurement in LHD where the quantization axis is taken in the magnetic field direction [9]. In LHD, it is known that the radial location of neutral hydrogen emissions is almost fixed irrespective of the plasma condition [16–18], and local *T*e and *n*e at the emission location are obtained from the radial *T*e and *n*e profiles measured by the Thomson scattering diagnostic system. In the experiment where *P* = −0.033 is obtained, the *T*<sup>e</sup> and *n*<sup>e</sup> at the Lyman-*α* emission location are found to be 32 eV and 9.6 × 1018 <sup>m</sup><sup>−</sup>3, respectively. Because the Thomson scattering diagnostic system for LHD measures light scattered by electrons moving predominantly in the direction perpendicular to the magnetic field, this *T*<sup>e</sup> value can be regarded as *T*<sup>⊥</sup> in our model calculation. We now have *T*<sup>⊥</sup> and *n*<sup>e</sup> at the emission location

as well as *P*. Figure <sup>5</sup> indicates that *T* can be determined when *T*<sup>⊥</sup> and *n*<sup>e</sup> are known so as to give the measured *<sup>P</sup>*. In the present case, *<sup>T</sup>* = 4.2 eV is derived.

The atomic model developed in Section 2 assumes a simple mechanism for the formation of the excited level population. Here, the validity of the present model is examined with a collisional-radiative model (CR-model) for atomic hydrogen [19] which treats energy levels resolved only by the principal quantum number *p*.

The CR-model solves coupled rate equations for all of the excited levels considered in the model under the quasi-steady-state condition [19] for determining the population distribution over the excited levels. Because the present plasma is in the ionizing state [20], each excited level population is expressed as

$$n(p) = \mathcal{R}\_1(p) n\_\text{\textquotedblleft}(1),\tag{20}$$

where *R*1(*p*) is called the population coefficient of the level *p* and is a function of *n*<sup>e</sup> and *T*e, and *n*(1) stands for the ground state density. The CR-model derives *R*1(*p*) for all the excited levels considered.

By using the results of the CR-model, we have evaluated breakdowns of the population flows from and to the *p* = 2 level. Figure 6a shows the *n*<sup>e</sup> dependence of component fractions of the inflow to the *p* = 2 level from other levels at *T*<sup>e</sup> = 10 eV. It is found that more than 90% of the inflow is dominated by the electron impact excitation, and the cascades from higher levels account for the remaining part of the inflow in the *n*e range of our interest, i.e., from 1018 m−<sup>3</sup> to 10<sup>19</sup> m−3. We have confirmed that the results are hardly changed at *T*<sup>e</sup> = 30 eV. This result justifies the assumption regarding the populating process of the *p* = 2 level in the present model in Section 2.

**Figure 6.** Fractions of breakdowns for the (**a**) population inflow to the level *p* = 2 and (**b**) population outflow from the level *p* = 2 as a function of *n*<sup>e</sup> at *T*<sup>e</sup> = 10 eV. The numbers represent the principal quantum number of the levels and the arrows indicate the transition direction. The hatched and open areas indicate the radiative and collisional transitions, respectively. The label "ioniz". means the ionization.

Similar results concerning the outflow from the *p* = 2 are shown in Figure 6b. It is found that the radiative decay to the ground state predominates over other processes in the *n*e range of our interest, which supports the assumption for the depopulating process from the *p* = 2 level in the model. Because collisional transition rates between different *l*-levels and *j*-levels are generally small as compared to those between *p*-levels under the conditions assumed here [21,22], the present model is regarded as a good approximation for the plasma considered.

In this paper, we have developed an atomic model for analyzing polarization states of the Lyman-*α* line where observations of the magnetically confined fusion plasma are borne in mind. The simplified model adopted is confirmed to be adequate through analyses for the populating mechanism of the *p* = 2 level which is the upper state of the Lyman-*α* line emission. An analysis of experimental data has been attempted with the present model, and the anisotropy is derived in terms of the difference in *T*e between the magnetic field direction and the direction perpendicular to it. It is finally noted that the model developed in this paper is dedicated to the Lyman-*α* line, and therefore the anisotropy can be diagnosed in the limited narrow region where the Lyman-*α* line emission mainly takes place. However, the methodology can be easily transferred to other emission lines of other atoms and ions, so that an appropriate emission line can be chosen depending on the plasma region where the interest is focused.

**Author Contributions:** Conceptualization, methodology, and writing, M.G.; formal analysis, N.R.; investigation, M.G. and N.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partly supported by JSPS KAKENHI Grant Number 18K03588 and by the National Institute for Fusion Science grant administrative budgets (ULHH028).

**Acknowledgments:** We thank Takashi Fujimoto for his helpful suggestions concerning the collisional process of polarization relaxation.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

#### **references**

