O4 -Symmetry-Based Non-Perturbative Analytical Calculations of the Effect of the Helical Trajectories of Electrons in Strongly Magnetized Plasmas on the Width of Hydrogen/Deuterium Spectral Lines

Abstract

The effects of the helical trajectories of the perturbing electrons in magnetized plasmas on the dynamical Stark width of hydrogen or deuterium spectral lines have been studied analytically in our previous two papers—specifically in the situation where the magnetic field B is so strong that the dynamical Stark width of these lines reduces to the so-called adiabatic Stark width because the so-called nonadiabatic Stark width is completely suppressed. This situation corresponds, for example, to DA and DBA white dwarfs. We obtained those analytical results by using the formalism of the so-called conventional (or standard) theory of the impact Stark broadening: namely, by performing calculations in the second order of the Dyson perturbation expansion. The primary outcome was that the dynamical Stark broadening was found to not depend on the magnetic field B (for sufficiently strong B). In the present paper, we use the O₄ symmetry of hydrogen atoms for performing the corresponding non-perturbative analytical calculations equivalent to accounting for all orders of the Dyson perturbation expansion. The results, obtained by using the O₄ symmetry of hydrogen atoms, differ from our previous ones not only quantitatively, but—most importantly—qualitatively. Namely, the dynamical Stark broadening does depend on the magnetic field B, even for strong B. These results should be important for revising the interpretation of the hydrogen Balmer lines observed in DA and DBA white dwarfs. We also address confusion in the literature on this subject.

1. Introduction

Among the analytical theories of the dynamical Stark broadening of hydrogen (or deuterium) spectral lines in plasmas, the most practically useful are the semiclassical theories that started from the pioneering works by Baranger [1,2]. These were then further developed in paper [3] into what was later called the conventional or standard theory, where it was assumed that the ion microfield was quasistatic (from the point of view of radiating atoms), while the electron microfield was dynamical and treated in the impact approximation. The primary feature of the impact approximation is treating the electron microfield as a sequence of completed binary collisions of the electrons with the radiating atom. The allowance for incomplete collisions was later made in papers [4,5,6], calling this version the unified theory.

In the intervening dozens of years, many further advances have been made, especially in treating the dynamical part of the ion microfield. We refer the readers to books [7,8,9,10,11,12,13,14,15,16], listed in chronological order, and the references therein.

The effects of the helical trajectories of the perturbing electrons in magnetized plasmas on the dynamical Stark width of hydrogen or deuterium spectral lines has been studied analytically for the first time in papers [17,18]—specifically in the situation where the magnetic field B is so strong that the dynamical Stark width of these lines reduces to the so-called adiabatic Stark width because the so-called nonadiabatic Stark width is completely suppressed. This situation corresponds, for example, to DA and DBA white dwarfs. The analytical results were obtained in papers [17,18] by using the formalism of the so-called conventional or standard theory of the impact Stark broadening [3,19]: namely, by performing calculations in the second order of the Dyson perturbation expansion. The primary outcome was that the dynamical Stark broadening was found not to depend on the magnetic field B (for sufficiently strong B).

In the present paper, we use the O₄ symmetry of hydrogen atoms for performing the corresponding non-perturbative analytical calculations equivalent to accounting for all orders of the Dyson perturbation expansion. The results, obtained by using the O₄ symmetry of hydrogen atoms, differ from those obtained in papers [17,18] not only quantitatively, but—most importantly—qualitatively: namely, the dynamical Stark broadening does depend on the magnetic field B, even for strong B.

2. Summary of Starting Formulas from Papers [17,18]

1. The radius vector R(t) of a perturbing electron traveling on a helical trajectory in a magnetized plasma is

R(t) = v_ztB/B + ρ[1 + (r_Bp/ρ) cos(ω_Bt + φ)] + ρ × B [r_Bp/(ρB)] sin(ω_Bt + φ).(1)

In Equation (1), ρ is the impact parameter vector and ρ × B is its cross-product (vector product) with the magnetic field B; v_z is the component of the electron velocity parallel to B (the z-axis is chosen parallel to B). Other notations are presented as follows:

r_Bp = v_p/ω_B, ω_B = eB/(m_ec).(2)

In Equation (2), ω_B is the Larmor frequency and v_p is the component of the electron velocity perpendicular to B.

2. The electric field created by the perturbing electron at the location of the radiating hydrogen or deuterium atom:

E(t) = eR(t)/[R(t)]³.(3)

The z-projection of this electric field is

E_z(t) = e(v_zt)/[ρ² + v_z²t² + v_p²/ω_B² + 2(ρv_p/ω_B) cos(ω_Bt + φ)]^3/2,(4)

3. The nonadiabatic contribution to the dynamical Stark width, i.e., the contribution caused by the component of the field E(t) perpendicular to B, is completely suppressed if

B > B_cr = 2T_e/(3|X_αβ|eλ_c) = 9.2 × 10²T_e(eV)/|X_αβ| Tesla,(5)

where λ_c = ħ/(m_ec) = 2.426 × 10⁻¹⁰ cm is the Compton wavelength of electrons. In Equation (5), T_e is the electron temperature, and

X_αβ = |n_a(n₁ − n₂)_α − n_b(n₁ − n₂)_β|.(6)

In Equation (6), n is the principal quantum number, and n₁ and n₂ are the parabolic quantum numbers of the upper (aα) and lower (bβ) Stark states involved in the radiative transition.

4. The adiabatic part Φ_ad of the electron broadening operator Φ_ab (and the corresponding Stark width) is controlled by E_z(t), given by Equation (4).

5. For the helical trajectories of perturbing electrons, the starting formula for the adiabatic Stark width Γ_αβ = − Re[_αβ(Φ_ad)_βα] for the line component, corresponding to the radiative transition between the upper Stark sublevel α and the lower Stark sublevel β, is

(7)${\Gamma }_{\mathsf{\alpha }\mathsf{\beta },\mathrm{h}\mathrm{e}\mathrm{l}}={\mathrm{N}}_{\mathrm{e}}{\displaystyle \underset{-\infty }{\overset{\infty }{\int }}\mathrm{d}{\mathrm{v}}_{\mathrm{z}}{\mathrm{f}}_1\left({\mathrm{v}}_{\mathrm{z}}\right)}{\displaystyle \underset{0}{\overset{\infty }{\int }}\mathrm{d}{\mathrm{v}}_{\mathrm{p}}{\mathrm{f}}_2\left({\mathrm{v}}_{\mathrm{p}}\right)}{\left({{\mathrm{v}}_{\mathrm{z}}}^2+{{\mathrm{v}}_{\mathrm{p}}}^2\right)}^{1/2}\mathsf{\sigma }\left({\mathrm{v}}_{\mathrm{z}},{\mathrm{v}}_{\mathrm{p}}\right).$

In Equation (7), f₂(v_p) is the two-dimensional Maxwell distribution and f₁(v_z) is its one-dimensional counterpart. The operator σ(v_z, v_p) discussed below has the physical meaning of the cross-section of the “optical” collisions causing the dynamical Stark broadening of the spectral line.

3. Non-Perturbative Analytical Calculations

Let us start in a general way. We consider an atomic system described by the following Hamiltonian:

H(t) = H₀ + V(t), V(t) = −dE(t).(8)

In Equation (8), H₀ is the unperturbed Hamiltonian, d is the electric dipole moment operator, and dE(t) is its scalar product (also known as the dot-product) with a time-dependent electric field E(t).

The central idea is to split the interaction V(t) into two parts

V(t) = V₁(t) + V₂(t)(9)

and to rewrite the Hamiltonian from Equation (8) as follows:

H(t) = H₁(t) + V₁(t), H₁(t) = H₀ + V₁(t),(10)

where the extended “unperturbed” Hamiltonian H₁(t) can be diagonalized analytically at any instant of time (at least, in a subspace of a fixed principal quantum number n). In particular, for hydrogenic atoms or ions, the analytical diagonalization of the extended “unperturbed” Hamiltonian H₁(t) is possible due to the higher than geometrical symmetry (O₄ symmetry) of the system described by H₁(t).

We remind that one of the consequences of the O₄ symmetry of hydrogenic atoms or ions is that the separation of variables is possible in several different types of coordinates: in the spherical coordinates, in the parabolic coordinates, and in the elliptical coordinates. For the present paper, the most convenient is the utilization of the parabolic coordinates, in which the projection d_z of the dipole moment operator on the magnetic field can be diagonalized within any subspace of the fixed principal quantum number n.

We choose the z-axis along the magnetic field B and use the following specific breakdown of the interaction V(t) into two parts: the adiabatic part

V₁(t) = −d_zE_z(t),(11)

and the nonadiabatic part

V₂(t) = −d_perpE_perp(t),(12)

where d_perp and E_perp are the components of the electric dipole moment and of the electric field, respectively, perpendicular to the magnetic field B. Under the condition (5), the nonadiabatic part of the interaction is completely suppressed, so that only the adiabatic part V₁(t) remains.

We calculate the adiabatic contribution to the dynamical Stark broadening of hydrogenic spectral lines exactly, nonperturbatively (rather than in the second order of the Dyson perturbation expansion) [20,21,22]. This is achieved along the lines of the so-called old adiabatic theory—see, e.g., papers [23,24] and Appendix A of the present paper. In the formalism of the old adiabatic theory, the cross-section σ of the optical collisions has the following form:

(13)$\mathsf{\sigma }=2\mathsf{\pi }{\displaystyle \underset{0}{\overset{\infty }{\int }}\mathrm{d}\varrho \varrho }⟨1-\cos \left[{\displaystyle \underset{-\infty }{\overset{\infty }{\int }}{\mathrm{dtd}}_{\mathrm{z}}{\mathrm{E}}_{\mathrm{z}}(\mathrm{t})/\mathrm{ħ}}\right]⟩.$

In Equation (13), the symbol <…> stands for the average over angular or phase variables; d_z is the matrix element of the z-projection of the electric dipole moment:

d_z/ħ = 3X_αβħ/(2m_ee).(14)

The integral over time in Equation (13) vanishes for the odd part of E_z(t). Thus, it suffices to perform the integration only for the even part E_z,even(t) of E_z(t):

E_z(t)_even = [E_z(t) + E_z(−t)]/2.(15)

The latter was calculated in papers [17,18] as follows.

The integral over the impact parameter ρ was broken into two parts

σ = σ₁ + σ₂,(16)

σ₁ corresponding to the integral from ρ₀ to infinity and σ₂ corresponding to the integral from zero to ρ₀, where

ρ₀ = v_p/ω_B.(17)

For calculating σ₁, E_z,even(t) was expanded in terms of the small parameter ρ₀/ρ = v_p/(ω_Bρ). After keeping the first non-vanishing term of the expansion, the following was obtained:

E_z(t)_even = (sin φ) (3eρv_pv_z/ω_B) t[sin(ω_Bt)]/(ρ² + v_z²t²)^5/2.(18)

For calculating σ₂, E_z,even(t) was expanded in terms of the small parameter ρ/ρ₀ = ω_Bρ/v_p. After keeping the first non-vanishing term of the expansion, the following was obtained:

E_z(t)_even = (sin φ) (3eρv_pv_z/ω_B) t[sin(ω_Bt)]/(v_p²/ω_B² + v_z²t²)^5/2.(19)

Using expressions (18) and (19), obtained in papers [17,18], we proceed now to calculating

(20)${\mathsf{\sigma }}_1=2\mathsf{\pi }{\displaystyle \underset{{\varrho }_0}{\overset{\infty }{\int }}\mathrm{d}\varrho \varrho }⟨1-\cos \left[{\displaystyle \underset{-\infty }{\overset{\infty }{\int }}{\mathrm{dtd}}_{\mathrm{z}}{\mathrm{E}}_{\mathrm{z}}{(\mathrm{t})}_{\mathrm{even}}/\mathrm{ħ}}\right]⟩$

and

(21)${\mathsf{\sigma }}_2=2\mathsf{\pi }{\displaystyle \underset{0}{\overset{{\varrho }_0}{\int }}\mathrm{d}\varrho \varrho }⟨1-\cos \left[{\displaystyle \underset{-\infty }{\overset{\infty }{\int }}{\mathrm{dtd}}_{\mathrm{z}}{\mathrm{E}}_{\mathrm{z}}{(\mathrm{t})}_{\mathrm{even}}/\mathrm{ħ}}\right]⟩.$

For σ₁, after calculating the integral over time, the expression within the symbol <…> in Equation (20) becomes

1 − cos[(sinφ)3X_αβħv_pω_BK₁(ω_Bρ/|v_z|)/(m_e|v_z|³)],(22)

where K₁(s) is the modified Bessel function of the 2nd kind. Then, the averaging of the expression (22) over φ yields

<…> = 1 − J₀[3X_αβħv_pω_BK₁(ω_Bρ/|v_z|)/(m_e|v_z|³)],(23)

where J₀(u) is the Bessel function. For the range of ρ under consideration, the Bessel function K₁(ω_Bρ/|v_z|) is exponentially small, so that the argument of the Bessel function J₀[…] is much smaller than unity. Therefore, Equation (23) can be simplified to the following:

<…> ≈ [3X_αβħv_pω_BK₁(ω_Bρ/|v_z|)/(2m_e|v_z|³)]².(24)

Then, the integration over impact parameters in Equation (20) yields the following:

σ₁ ≈ (π/4)^3/2[3X_αβħv_p/(m_ev_z²)]² MeijerG[{{},{3/2}}, {{0, 0, 2},{}}, v_p²/v_z²],(25)

where MeijerG[…] is the Meijer G-function.

Equation (24) for σ₁ coincides with σ₁ from papers [17,18]. It does not depend on the magnetic field B (provided that B satisfies the condition (5)). However, below it is shown that our nonperturbative result for σ₂ significantly differs—both quantitatively and qualitatively—from the perturbatively obtained result for σ₂ from papers [17,18]. In particular, our nonperturbative result depends on the magnetic field B, while the perturbative result from papers [17,18] did not depend on B.

For σ₂, after calculating the integral over time, the expression within the symbol <…> in Equation (21) becomes

1 − cos[(sinφ)3X_αβħv_pω_B²K₁(v_p/|v_z|) ρ/(m_e|v_z|³)].(26)

Then, the averaging of the expression (26) over φ yields

<…> = 1 − J₀[3X_αβħω_B²K₁(v_p/|v_z|) ρ/(m_e|v_z|³)].(27)

The subsequent integration over impact parameters in Equation (27) leads to the following result for σ₂:

σ₂ = (πv_p²/ω_B²)[1 − J₁(w)/w], w = 3X_αβħv_pω_BK₁(v_p/|v_z|)/(m_e|v_z|³).(28)

Obviously, σ₂ depends on ω_B and thus on the magnetic field B.

For presenting the results in plots, we calculate the effective average values <v_p> and <|v_z|>, as well as their ratio, as follows (compare to Equation (7)):

(29)$<{\mathrm{v}}_{\mathrm{p}}>={\displaystyle \underset{-\infty }{\overset{\infty }{\int }}\mathrm{d}{\mathrm{v}}_{\mathrm{z}}{\mathrm{f}}_1\left({\mathrm{v}}_{\mathrm{z}}\right)}{\displaystyle \underset{0}{\overset{\infty }{\int }}\mathrm{d}{\mathrm{v}}_{\mathrm{p}}{\mathrm{f}}_2\left({\mathrm{v}}_{\mathrm{p}}\right)}{\left({{\mathrm{v}}_{\mathrm{z}}}^2+{{\mathrm{v}}_{\mathrm{p}}}^2\right)}^{1/2}{\mathrm{v}}_{\mathrm{p}},$

(30)$<{\mathrm{v}}_{\mathrm{p}}>={\displaystyle \underset{-\infty }{\overset{\infty }{\int }}\mathrm{d}{\mathrm{v}}_{\mathrm{z}}{\mathrm{f}}_1\left({\mathrm{v}}_{\mathrm{z}}\right)}{\displaystyle \underset{0}{\overset{\infty }{\int }}\mathrm{d}{\mathrm{v}}_{\mathrm{p}}{\mathrm{f}}_2\left({\mathrm{v}}_{\mathrm{p}}\right)}{\left({{\mathrm{v}}_{\mathrm{z}}}^2+{{\mathrm{v}}_{\mathrm{p}}}^2\right)}^{1/2}|{\mathrm{v}}_{\mathrm{z}}|.$

As a result, we find

<v_p> = 1.60 <v_T>, <|v_z|> = 0.665 <v_T>, <v_p>/<|v_z|> = 2.41,(31)

where v_T = (2T_e/m_e)^1/2 is the mean thermal velocity of plasma electrons. Then, we use these effective average values in Equation (25) for σ₁ and in Equation (28) for σ₂ in the subsequent plots where all quantities will be in atomic units.

Figure 1 shows the dependence of the ratio Γ/(N_ev_T), where Γ is the adiabatic Stark width calculated by combining Equations (7), (25), and (28), on the Larmor frequency ω_B and on the mean thermal velocity of electrons v_T for X_αβ = 6, X_αβ being defined in Equation (9). We note that X_αβ = 6 corresponds to either n_a = 3, (n₁ − n₂)_α = 2, n_b = 1 (which is a component of the Ly-beta line), or to n_a = 4, (n₁ − n₂)_α = 2, n_b = 2, (n₁ − n₂)_β = 1, (which is a component of the Balmer-beta line), or to n_a = 4, (n₁ − n₂)_α = 3, n_b = 3, (n₁ − n₂)_β = 2, (which is a component of the Paschen-alpha line), or to n_a = 6, (n₁ − n₂)_α = 1, n_b = 1, (which is a component of Lyman-epsilon line), or to n_a = 6, (n₁ − n₂)_α = 2, n_b = 3, (n₁ − n₂)_β = 2, (which is a component of the Paschen-gamma line), or to n_a = 4, (n₁ − n₂)_α = 2, n_b = 3, (n₁ − n₂)_β = 2, (which is a component of the Paschen-alpha line), or to n_a = 6, (n₁ − n₂)_α = 3, n_b = 4, (n₁ − n₂)_β = 3, (which is a component of the Brackett-beta line). Thus, the value of X_αβ = 6 represents Stark components of seven spectral lines, which is why we chose this value for Figure 1.

Figure 2 displays the dependence of the ratio Γ/(N_ev_T) on the Larmor frequency ω_B for the mean thermal velocity of electrons v_T = 0.2 a.u., corresponding to 2.7 eV, for X_αβ = 6 (solid line). The corresponding perturbative result from papers [17,18] is shown by the dashed line.

From Figure 2, it is seen that the perturbative calculation from papers [17,18] overestimated the Stark width. It also illustrates that the perturbative calculation from papers [17,18] did not produce any dependence of the Stark width on the magnetic field, while the nonperturbative (more accurate calculation) demonstrates a significant dependence of the Stark width on the magnetic field.

The physical reason for the decrease in the dynamical Stark width as the magnetic field becomes stronger is the following. The dynamical Stark width is Γ~N_ev_Tσ, where σ is the cross-section of the optical collisions responsible for Stark broadening of hydrogen or deuterium lines [9]. In its turn, σ~(v_T/Ω)², where Ω is the characteristic frequency of the variation of the electric field of the perturbing ions at the location of the radiating atom. Thus, the dynamical Stark width Γ is inversely proportional to Ω². In non or weakly magnetized plasmas, Ω~ω_W~T_e/(n²ħ), where ω_W is the Weisskopf frequency [24]. In strongly magnetized plasmas, as the Larmor frequency ω_B starts exceeding ω_W, then the characteristic frequency of the variation of the electric field of the perturbing ions at the location of the radiating atom becomes Ω~ω_B and the dynamical Stark width decreases with the growth of the magnetic field.

As for the dependence of the dynamical Stark width Γ on the mean thermal velocity of plasma electrons v_T, physically it is the competition of the following two factors. On the one hand, at v_T = 0, the dynamical Stark width vanishes because it is dynamical, while there no dynamics at v_T = 0. On the other hand, at relatively large v_T, the dynamical Stark width Γ decreases as v_T increases. This is because at relatively large v_T, so that the Weisskopf frequency exceeds the Larmor frequency, the dynamical Stark width Γ, being proportional to v_Tσ(v_T), becomes inversely proportional to v_T. Consequently, somewhere between v_T = 0 and relatively large v_T, Γ reaches a maximum, as it is seen in Figure 1 at a fixed ω_B.

4. Discussion of the Results

We showed how far the analytical semiclassical theory of the dynamical Stark broadening of hydrogen or deuterium lines in plasmas has progressed. While the pioneering works [1,2,3] presented the most simple theory, lots of subsequent improvements were made in the intervening dozens of years, as reflected, e.g., in works [4,5,6,7,8,9,10,11,12,13,14,15,16]. However, it was not until the year 2017 that the effect of helical rather than rectilinear trajectories of perturbing electrons in strongly magnetized plasmas on the dynamical Stark width of hydrogen or deuterium spectral lines was studied analytically for the first time, in paper [17].

Specifically, in papers [17,18], the situation was studied where the magnetic field B is so strong that the dynamical Stark width of these lines reduces to the so-called adiabatic Stark width because the so-called nonadiabatic Stark width is completely suppressed. The analytical calculations in papers [17,18] were perturbative: they were performed in the second order of the Dyson perturbation expansion.

In distinction, in the present paper we used the O₄ symmetry of hydrogen or deuterium atoms for performing non-perturbative analytical calculations equivalent to accounting for all orders of the Dyson perturbation expansion in the spirit of the generalized theory of the dynamical Stark broadening [20,21,22].

The results of the present paper differ from those obtained in papers [17,18], not only quantitatively, but—most importantly—qualitatively. Namely, the dynamical Stark broadening does depend on the magnetic field B even for strong B, while in papers [17,18] the dynamical Stark broadening did not depend on B. In addition to this, the non-perturbative calculations yield the dynamical Stark width by about an order of magnitude smaller than the perturbative calculations from papers [17,18], as seen, e.g., in Figure 2.

The results of the present paper should be important for revising the interpretation of the hydrogen Balmer lines observed from DA and DBA white dwarfs. According to observations, the magnetic field in plasmas of white dwarfs can range from 10³ Tesla to 10⁵ Tesla (see, e.g., papers [25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42], listed in chronological order and the references therein), thus easily exceeding the critical value from Equation (5). Indeed, for the typical electron temperature T_e~1 eV in the white dwarfs plasma emitting hydrogen lines, Equation (5) yields B_cr~10³/|X_αβ| Tesla < 10³ Tesla.

Finally, we note that in the literature there is a confusion concerning the effect that the helical trajectories of the perturbing electrons have on the dynamical Stark width of hydrogen spectral lines. We clarify this confusion in Appendix B of the present paper.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflicts of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures

Figure 1. Dependence of the ratio Γ/(NevT), where Γ is the adiabatic Stark width calculated by combining Equations (7), (25), and (28), on the Larmor frequency ωB and on the mean thermal velocity of electrons vT for Xαβ = 6, Xαβ being defined in Equation (9). All quantities are in atomic units.

Figure 2. Dependence of the ratio Γ/(NevT) on the Larmor frequency ωB for the mean thermal velocity of electrons vT = 0.2 a.u., corresponding to 2.7 eV, for Xαβ = 6 (solid line). The corresponding perturbative result from papers [1,2] is shown by the dashed line. All quantities are in atomic units.

Appendix A. Some Details on the Old Adiabatic Theory of the Stark Broadening of Spectral Lines in Plasmas

The adiabatic line broadening theory proceeds from the assumption that the electric field E(t) of the perturbing charged particles causes only the phase modulation of the atomic oscillator (i.e., a change in the phase of the wave function of the radiating atom)—see, e.g., review [34]. Then, the correlation function has the following form:(A1)$$\mathrm{K}(\mathsf{\tau })=<\mathrm{e}\mathrm{x}\mathrm{p}\left[\mathrm{i}\underset{0}{\overset{\mathsf{\tau }}{\int }}\mathsf{\kappa }(\mathrm{t})\mathrm{d}\mathrm{t}\right]>.$$

In Equation (A1), κ(t) = CE(t)/e,(A2) where C is the Stark constant of the line component in the parabolic quantization.

Next, it is assumed that the phase change κ(t) is the sum of the phase changes κ_j(t) due to the individual independent perturbing particles: j = 1, 2, …, N. Then, the correlation function factorizes into the product of the corresponding contributions from the individual particles:(A3)$$\mathrm{K}(\mathsf{\tau })=<\exp \left[\mathrm{i}{\displaystyle \sum _{\mathrm{j}=1}^{\mathrm{N}}\underset{0}{\overset{\mathrm{\tau }}{\int }}{\mathsf{\kappa }}_{\mathrm{j}}(\mathrm{t})\mathrm{d}\mathrm{t}}\right]{>}_{\mathrm{N}}=<\underset{\mathrm{j}=1}{\overset{\mathrm{N}}{\mathsf{\Pi }}}\mathrm{e}\mathrm{x}\mathrm{p}\left[\mathrm{i}\underset{0}{\overset{\mathsf{\tau }}{\int }}{\mathsf{\kappa }}_{\mathrm{j}}(\mathrm{t})\mathrm{d}\mathrm{t}\right]{>}_{\mathrm{N}}=<\mathrm{e}\mathrm{x}\mathrm{p}\left[\mathrm{i}\underset{0}{\overset{\mathsf{\tau }}{\int }}\mathsf{\kappa }(\mathrm{t})\mathrm{d}\mathrm{t}\right]{{>}_1}^{\mathrm{N}}.$$

In Equation (A3), the symbol <…>_N stands for the average over the phase volume of all particles, while the symbol <…>₁ stands for the average over the phase volume of one particle. The latter average incorporates the integration with respect to the coordinates dr₀/V and with respect to the velocity distribution f(v)dv of the perturbing particles (V being the spatial volume):(A4)$$\mathrm{K}(\mathsf{\tau })={\left\{\left[1-(1/\mathrm{V}){\displaystyle \int \mathbf{d}{\mathbf{r}}_0}{\displaystyle \int \mathrm{f}(\mathrm{v})\mathrm{d}\mathrm{v}}\right]\left[1-\mathrm{e}\mathrm{x}\mathrm{p}\left(\mathrm{i}\underset{0}{\overset{\mathsf{\tau }}{\int }}\mathsf{\kappa }(\mathrm{t})\mathrm{d}\mathrm{t}\right)\right]\right\}}^{\mathrm{N}=\mathrm{N}\mathrm{p}\mathrm{V}},$$ where N_p is the density of the perturbing particles.

In the limit where N → ∞, V → ∞, but N_p = N/V = const, one obtains K(τ) = exp{−N_pV(τ)],(A5) where V(τ) is the so-called “collision volume:(A6)$$\mathrm{R}\mathrm{e}[\mathrm{V}(\mathsf{\tau })]={\displaystyle \int \mathbf{d}\mathbf{r}\mathrm{f}(\mathrm{v})}\underset{0}{\overset{\infty }{\int }}2\mathsf{\pi }\varrho \mathrm{d}\varrho \underset{-\infty }{\overset{\infty }{\int }}\mathrm{d}\mathrm{z}<1-\cos [\underset{-\infty }{\overset{\infty }{\int }}\mathrm{d}\mathrm{t}\mathrm{C}{\mathrm{E}}_{\mathrm{z}}(\mathrm{t})]>.$$

In Equation (A6), the transformation to the cylindrical coordinates dr₀ = 2πρdρdz has been made.

For the relatively large time τ >> Ω, where Ω is the characteristic frequency of the variation of the electric field of the perturbing particle at the location of the radiating atom, one obtains Re[V(τ)] ≈ (dV/dτ)|_τ=∞ = vστ,(A7) where σ is given by Equation (13).

Appendix B. Clarification of the Confusion in the Literature on This Subject

The analytical results from our paper [17] were obtained in frames of the so-called conventional (or standard) theory of the impact Stark broadening [3,19]. Despite this, in paper [43] it was incorrectly stated that the results in paper [1] were obtained in the frames of the so-called generalized theory of Stark broadening.

Most importantly, in paper [43] it was incorrectly stated that presumably in our paper [17] it was predicted analytically that the allowance for the helical trajectories of the perturbing electrons (HTPE) leads to a dramatic width increase in the lines Balmer-beta, Balmer-delta, and Balmer-epsilon at high densities, while simulations in paper [43] yielded a decrease in the corresponding widths. For supporting this statement, in paper [43] there were specific examples of the Balmer-beta, Balmer-delta, and Balmer-epsilon lines at the electron temperature T_e = 1 eV and the electron density N_e = 2 × 10¹⁷ cm⁻³.

However, in reality, according to Equations (18) and (47), and Figure 1 from our paper [17], whether the allowance for HTPE increases or decreases the width of the Stark components of any hydrogen line depends on the value of the following dimensionless parameter D = 5.57 × 10⁻¹¹|X_αβ|[N_e(cm⁻³)]^1/2/T_e(eV),(A8) where X_αβ = n_aq_α − n_bq_β, q_α = (n₁ − n₂)_α, q_α = (n₁ − n₂)_α.(A9)

In Equation (A9), n is the principal quantum number, n₁ and n₂ are the parabolic quantum numbers (while q is often called the electric quantum number); we noted in paper [1] that X_αβ is the standard label of the Stark component of hydrogen or deuterium spectral lines corresponding to the radiative transition between the upper (α) and lower (β) Stark sublevels. According to Equation (47) and Figure 1 from paper [17], the allowance for HTPE increases the width of a Stark component if D > 0.44, but decreases its width if D < 0.44.

For the plasma parameters chosen in paper [43], Equation (A8) simplifies to D = 0.025 |X_αβ|.(A10)

The critical value of D = 0.44 corresponds in this case to the critical value of |X_αβ| = 18, so that according to paper [1], the allowance for HTPE increases the width of Stark components having |X_αβ| > 18, but decreases the width of Stark components having |X_αβ| < 18.

For the Balmer-beta line, all intense Stark components have a |X_αβ| of no more than 10. So, the actual prediction from paper [17] for the Balmer-beta line at the plasma parameters chosen in paper [43] is the decrease in the Stark width (which can be also seen in Figures 2 and 3 from our paper [17]), rather than the increase in the Stark width (incorrectly stated in paper [43] with respect to the predictions from paper [17].

For the Balmer-delta line, the most intense Stark component has |X_αβ| = 6. So, based on the results from paper [17], for the plasma parameters chosen in paper [43], one should not expect the increase in the Stark width—contrary to the incorrect statement from paper [43] with respect to the predictions from our paper [17].

For the Balmer-epsilon line, the most intense Stark component has |X_αβ| = 14. So, again, based on the results from paper [17], for the plasma parameters chosen in paper [28], one should not expect the increase in the Stark width—contrary to the incorrect statement from paper [43] with respect to the predictions from our paper [17].

The statements from our paper [17] concerning the increase in the Stark width of the Balmer-delta and higher lines due to HTPE related to the electron densities N_e > 10¹⁸ cm⁻³, i.e., electron densities much higher than the only one value of N_e chosen in simulations from paper [43]. This situation is yet another demonstration of the superiority of analytical results over simulations: the analytical results are valid for a broad range of the electron density, while the simulations from paper [43] were performed for only one value of the electron density N_e = 2 × 10¹⁷ cm⁻³.

References

1. Baranger, M. Simplified Quantum-Mechanical Theory of Pressure Broadening. Phys. Rev.; 1958; 111, pp. 481-493. [DOI: https://dx.doi.org/10.1103/PhysRev.111.481]

2. Baranger, M. Problem of Overlapping Lines in the Theory of Pressure Broadening. Phys. Rev.; 1958; 111, pp. 494-504. [DOI: https://dx.doi.org/10.1103/PhysRev.111.494]

3. Kolb, A.; Griem, H.R. Theory of Line Broadening in Multiplet Spectra. Phys. Rev.; 1958; 111, pp. 514-521. [DOI: https://dx.doi.org/10.1103/PhysRev.111.514]

4. Smith, E.; Cooper, J.; Vidal, C. Unified Classical-Path Treatment of Stark Broadening in Plasmas. Phys. Rev.; 1969; 185, pp. 140-151. [DOI: https://dx.doi.org/10.1103/PhysRev.185.140]

5. Vidal, C.; Cooper, J.; Smith, E. Hydrogen Stark Broadening Calculations with the Unified Classical Path Theory. J. Quant. Spectrosc. Radiat. Transf.; 1970; 10, pp. 1011-1063. [DOI: https://dx.doi.org/10.1016/0022-4073(70)90121-4]

6. Vidal, C.; Cooper, J.; Smith, E. Unified Theory Calculations of Stark Broadened Hydrogen Lines Including Lower State Interactions. J. Quant. Spectrosc. Rad. Transf.; 1971; 11, pp. 263-281. [DOI: https://dx.doi.org/10.1016/0022-4073(71)90013-6]

7. Griem, H.R. Plasma Spectroscopy; McGraw-Hill: New York, NY, USA, 1964.

8. Griem, H.R. Spectral Line Broadening by Plasmas; Academic: New York, NY, USA, 1974.

9. Oks, E. Plasma Spectroscopy. The Influence of Microwave and Laser Fields; Springer: New York, NY, USA, 1995.

10. Griem, H.R. Principles of Plasma Spectroscopy; Cambridge University Press: Cambridge, UK, 1997.

11. Fujimoto, T. Plasma Spectroscopy; Clarendon: Oxford, UK, 2004.

12. Oks, E. Stark Broadening of Hydrogen and Hydrogenlike Spectral Lines in Plasmas: The Physical Insight; Alpha Science International: Oxford, UK, 2006.

13. Fujimoto, T.; Iwamae, A. Plasma Polarization Spectroscopy; Springer: Berlin/Heidelberg, Germany, 2008.

14. Ochkin, V.N. Spectroscopy of Low Temperature Plasma; Wiley: Weinheim, Germany, 2009.

15. Oks, E. Diagnostics of Laboratory and Astrophysical Plasmas Using Spectral Lines of One-, Two-, and Three-Electron Systems; World Scientific: Singapore, 2017.

16. Oks, E. Advances in X-ray Spectroscopy of Laser Plasmas; IOP Publishing: Bristol, UK, 2020.

17. Oks, E. Effect of Helical Trajectories of Electrons in Strongly Magnetized Plasmas on the Width of Hydrogen/Deuterium Spectral Lines: Analytical Results and Applications to White Dwarfs. Intern. Rev. At. Mol. Phys.; 2017; 8, pp. 61-72.

18. Oks, E. Review of Recent Advances in the Analytical Theory of Stark Broadening of Hydrogenic Spectral Lines in Plasmas: Applications to Laboratory Discharges and Astrophysical Objects. Atoms; 2018; 6, 50.Erratum in Atoms 2019, 7, 68 [DOI: https://dx.doi.org/10.3390/atoms6030050]

19. Kepple, P.; Griem, H.R. Improved Stark Profile Calculations for the Hydrogen Lines Hα, Hβ, Hγ and Hδ. Phys. Rev.; 1968; 173, 317. [DOI: https://dx.doi.org/10.1103/PhysRev.173.317]

20. Ispolatov, Y.; Oks, E. A Convergent Theory of Stark Broadening of Hydrogen Lines in Dense Plasmas. J. Quant. Spectr. Rad. Transf.; 1994; 51, pp. 129-138. [DOI: https://dx.doi.org/10.1016/0022-4073(94)90073-6]

21. Oks, E.; Derevianko, A.; Ispolatov, Y. A Generalized Theory of Stark Broadening of Hydrogen-Like Spectral Lines in Dense Plasmas. J. Quant. Spectr. Rad. Transf.; 1995; 54, pp. 307-316. [DOI: https://dx.doi.org/10.1016/0022-4073(95)00065-S]

22. Oks, E. Stark Widths of Hydrogen Spectral Lines in Plasmas: A Highly-Advanced Non-Simulative Semiclassical Theory and Tables. AIP Conf. Proc.; 2006; 874, pp. 19-34.

23. Anderson, P.W. A Method of Synthesis of the Statistical and Impact Theories of Pressure Broadening. Phys. Rev.; 1952; 86, 809. [DOI: https://dx.doi.org/10.1103/PhysRev.86.809]

24. Lisitsa, V.S. Stark Broadening of Hydrogen Lines in Plasmas. Sov. Phys. Uspekhi; 1977; 122, pp. 603-630. [DOI: https://dx.doi.org/10.1070/PU1977v020n07ABEH005446]

25. Angel, J.R.P. Magnetic White Dwarfs. Ann. Rev. Astron. Astrophys.; 1978; 16, pp. 487-519. [DOI: https://dx.doi.org/10.1146/annurev.aa.16.090178.002415]

26. Angel, J.R.P.; Liebert, J.; Stockman, H.S. The Optical Spectrum of Hydrogen at 160–350 Million Gauss in the White Dwarf GRW +70 deg 8247. Astrophys. J.; 1985; 292, pp. 260-266. [DOI: https://dx.doi.org/10.1086/163156]

27. Wunner, G.; Roesner, W.; Herold, H.; Ruder, H. Stationary Hydrogen Lines in White Dwarf Magnetic Fields and the Spectrum of the Magnetic Degenerate GRW +70 8247. Astron. Astrophys.; 1985; 149, pp. 102-108.

28. Henry, R.J.W.; Oconnell, R.F. Hydrogen Spectrum in Magnetic White Dwarfs: Hα, Hβ, and Hγ Transitions. Publ. Astron. Soc. Pac.; 1985; 97, pp. 333-339. [DOI: https://dx.doi.org/10.1086/131540]

29. Bergeron, P.; Wesemael, F.; Fontaine, G.; Liebert, J. On the Surface Composition of Cool, Hydrogen-Line White Dwarfs: Discovery of Helium in the Atmospheres of Cool DA Stars and Evidence for Convective Mixing. Astrophys. J.; 1990; 351, pp. L21-L24. [DOI: https://dx.doi.org/10.1086/185670]

30. MacDonald, J.; Vennes, S. How Much Hydrogen is There in a White Dwarf?. Astrophys. J.; 1991; 371, pp. 719-738. [DOI: https://dx.doi.org/10.1086/169937]

31. Bergeron, P.; Saffer, R.A.; Liebert, J. A Spectroscopic Determination of the Mass Distribution of DA White Dwarfs. Astrophys. J.; 1992; 394, pp. 228-247. [DOI: https://dx.doi.org/10.1086/171575]

32. Fassbinder, P.; Schweizer, W. Stationary Hydrogen Lines in Magnetic and Electric Fields of White Dwarf Stars. Astron. Astrophys.; 1996; 314, pp. 700-706.

33. Reimers, D.; Jordan, S.; Koester, D.; Bade, N. Discovery of Four White Dwarfs with Strong Magnetic Fields by the Hamburg/ESO Survey. Astron. Astrophys.; 1996; 311, 572.

34. Kowalski, P.M.; Saumon, D. Found: The Missing Blue Opacity in Atmosphere Models of Cool Hydrogen White Dwarfs. Astrophys. J.; 2006; 651, pp. L137-L140. [DOI: https://dx.doi.org/10.1086/509723]

35. Kawka, A.; Vennes, S.; Schmidt, G.D.; Wickramasinghe, D.T.; Koch, R. Spectropolarimetric Survey of Hydrogen-Rich White Dwarfs Stars. Astrophys. J.; 2007; 654, pp. 499-520. [DOI: https://dx.doi.org/10.1086/509072]

36. Tremblay, P.-E.; Bergeron, P. Spectroscopic Analysis of DA White Dwarfs: Stark Broadening of Hydrogen Lines Including Nonideal Effects. Astrophys. J.; 2009; 696, pp. 1755-1770. [DOI: https://dx.doi.org/10.1088/0004-637X/696/2/1755]

37. Külebi, B.; Jordan, S.; Euchner, F.; Gänsicke, B.T.; Hirsch, H. Analysis of Hydrogen-Rich Magnetic White Dwarfs Detected in the Sloan Digital Sky Survey. Astron. Astrophys.; 2009; 506, pp. 1341-1350. [DOI: https://dx.doi.org/10.1051/0004-6361/200912570]

38. Gianninas, A.; Bergeron, P.; Ruiz, M.T. A Spectroscopic Survey and Analysis of Bright, Hydrogen-Rich White Dwarfs. Astrophys. J.; 2011; 743, 138. [DOI: https://dx.doi.org/10.1088/0004-637X/743/2/138]

39. Franzon, B.; Schramm, S. Effects of Magnetic Fields in White Dwarfs. J. Phys. Conf. Ser.; 2017; 861, 012015. [DOI: https://dx.doi.org/10.1088/1742-6596/861/1/012015]

40. Rolland, B.; Bergeron, P.; Fontaine, G. On the Spectral Evolution of Helium-atmosphere White Dwarfs Showing Traces of Hydrogen. Astrophys. J.; 2018; 857, 56. [DOI: https://dx.doi.org/10.3847/1538-4357/aab713]

41. Bohlin, R.C.; Hubeny, I.; Rauch, T. New Grids of Pure-hydrogen White Dwarf NLTE Model Atmospheres and the HST/STIS Flux Calibration. Astronom. J.; 2020; 160, 21. [DOI: https://dx.doi.org/10.3847/1538-3881/ab94b4]

42. Hardy, F.; Dufour, P.; Jordan, S. Spectrophotometric Analysis of Magnetic White Dwarf—I. Hydrogen-rich Compositions. Mon. Not. R. Astron. Soc.; 2023; 520, pp. 6111-6134. [DOI: https://dx.doi.org/10.1093/mnras/stad196]

43. Alexiou, S. Effects of Spiralling Trajectories on White Dwarf Spectra: High Rydberg States. Atoms; 2023; 11, 141. [DOI: https://dx.doi.org/10.3390/atoms11110141]