Stark-Zeeman Line Shape Modeling for Magnetic White Dwarf and Tokamak Edge Plasmas : Common Challenges

The shape of atomic spectral lines in plasmas contains information on the plasma parameters, and can be used as a diagnostic tool. Under specific conditions, the plasma located at the edge of tokamaks has parameters similar to those in magnetic white dwarf stellar atmospheres, which suggests that the same line shape models can be used. A problem common to tokamak and magnetic white dwarfs concerns the modeling of Stark broadening of hydrogen lines in the presence of an external magnetic field and the related Zeeman effect. In this work, we focus on a selection of issues relevant to Stark broadening in magnetized hydrogen plasmas. Various line shape models are presented and discussed through applications to ideal cases.


Introduction
Tokamaks are devices of toroidal shape employed in controlled nuclear fusion research for the confinement of a burning hydrogen plasma, with ion and electron temperatures in the order of 10 8 K.The ITER project (www.iter.org)aims to demonstrate the principle of producing more energy from the fusion process than is used to initiate it; this large-scale tokamak is presently under construction in France, and it should be operational for first plasma discharge by the end of 2025 [1].In order to support the operation of the machine, an extensive set of spectroscopic measurements is planned [2]; for example, passive spectroscopy of the Balmer series will provide information on the isotopic proportion of the fueling gas (N T /N D , N H /N D ) in the divertor region [3], where a large amount of neutrals is expected.Under specific conditions (plasma in "detached" regime), the electron density in the divertor can attain values larger than 10 14 cm −3 , so that the Stark effect due to the plasma microfield becomes the dominant line-broadening mechanism on lines with a high principal quantum number.Recent observations of Dε (transition 7 → 2 of deuterium) in the divertor of ASDEX-Upgrade have been used for the determination of the electron density, using the Stark broadening of the line [4].Balmer lines with a low principal quantum number like Dα are also affected by the microfield and Stark broadening enters into competition with both Doppler broadening and Zeeman splitting [5].The need for accuracy in spectroscopic diagnostics of tokamak edge plasmas has prompted an interest in the development of line broadening models accounting for the simultaneous action of electric and magnetic fields on atomic energy levels, e.g., [6].Recently, a line shape database for the first Balmer lines Atoms 2017, 5, 36 2 of 10 accounting for the Stark and Zeeman effects has been devised for tokamak spectroscopy applications using computer simulations [7].Similar databases for Stark broadening exist in astrophysics for stellar atmosphere diagnostic applications, e.g., STARK-B (stark-b.obspm.fr)[8].Magnetic white dwarfs are of particular interest because they present Balmer spectra with Zeeman triplets [9,10] just as in tokamaks, and atmosphere models indicate zones where the electron density value is comparable to that in detached divertor plasmas [11].Within this context, modeling work by tokamak spectroscopists in conjunction with astrophysicists aimed at improving Stark line shape databases is presently ongoing.In this paper, we report on models and techniques used for the description of Stark-Zeeman line shapes.A special emphasis is put on ion dynamics and collision operator models.

Stark Broadening Formalism
The formalism used in Stark line shape modeling involves atomic physics and statistical physics.An atom immersed in a plasma and emitting a photon is considered; during photon emission, the energy levels are perturbed due to the presence of the charged particles surrounding the atom, and this perturbation results in a splitting and broadening of spectral lines.The Doppler effect, which stems from the atom's motion, is not considered, hereafter, for the sake of clarity (it can be accounted for through convolution).According to standard textbooks (e.g., [12]), a line shape I(ω) is given by the Fourier transform of the atomic dipole autocorrelation function, a quantity described in the framework of quantum mechanics Here, ρ is the restriction of the density operator to the atomic Hilbert space evaluated at initial time, d•ε is the dipole projected onto the polarization vector ε, the indices α, α' (resp.β, β') are used for matrix elements and denote states in the upper (resp.lower) level, the brackets { . . .} denote an average over the perturber trajectories, and U(t) is the evolution operator.It obeys the time-dependent Schrödinger equation Here H 0 is the Hamiltonian accounting for the atomic energy level structure and V(t) = −d.F(t) is the time-dependent Stark effect term resulting from the action of the microscopic electric field F(t) and written in the Schrödinger picture, here.When this term is neglected, the Schrödinger equation has the trivial solution U(t) = exp(−iH 0 t/h), which shows, using Equations ( 2) and (1), that I(ω) reduces to a set of delta functions (or Lorentzian functions if the natural broadening is retained).By contrast, the case where F(t) is significant is much trickier, because there is no general exact analytical solution.The time-dependent perturbation theory yields a formal expansion (Dyson series), which is not applicable in explicit calculations because of the non-commutation of the interaction term at different times (time-ordering problem).This concerns, in particular, the microfield due to ions.Several models, based on suitable approximations, have been developed in such a way to provide an analytical expression for the line shape (e.g., the impact and static approximations, the model microfield method, etc.).On the other hand, computer simulations involving the numerical integration of the Schrödinger Equation (3) [13,14] provide reference line profiles and can serve as a benchmark for testing models.In the following, we focus on the ion dynamics issue, which is important for hydrogen lines with a low principal quantum number.

Ion Broadening at the Impact Limit
Modeling work has been done over the last decade, in order to set up an analytical formula for the Lyman α and Lyman β line shapes for opacity investigations in tokamak divertor plasmas ( [15] and Refs.therein).Kinetic simulations of radiation transport using the EIRENE code (www.eirene.de)have indicated the possibility of an alteration of the ionization-recombination balance due to the Lyman opacity [16].The code employs a Monte Carlo method with random walks and related estimators; line shapes in the atomic frame of reference, i.e., without Doppler broadening, are involved in the sampling of photon frequencies.Because of the low principal quantum numbers (n = 2 and 3) involved in Lyman α and Lyman β, the ion dynamics is important, and a collision operator can be used for ions as a first approach for an estimate of the Stark broadening.The model of Griem, Kolb, and Shen (GKS, [17]), developed for electron broadening in hydrogen plasmas, can be adapted in a simple way for ions.Consider a single ion perturbing the emitter (binary interaction); the collision operator K reads where the integrals are done over the velocity (module) v and the impact parameter (module) b, N and f stand for the perturber's density and velocity distribution function, the brackets denote an angular average and S is the scattering matrix corresponding to a binary collision.The emitter motion is accounted for through a reduced mass model.The upper bound λ D is the Debye length, which accounts phenomenologically for the screening of the Coulomb field due to the electrons.In the GKS model, the integration domain is restricted in such a way that S can be expanded at the second order in the Stark perturbation V(t) = −d.F(t) ("weak collisions"); the integrals can be performed explicitly in the case of a Coulomb field and the resulting collision operator reads (notations are different from those in [17]) Here, v 0 = √ (2k B T/m red ) is the thermal velocity, accounting for reduced mass; b W = (2/3) 1/2 hn 2 /m e v 0 is the Weisskopf radius, which discriminates between weak (b > b W ) and strong (b < b W ) collisions (sometimes this is defined without the (2/3) 1/2 factor); r 2 is the square of the atomic electron position operator, in atomic units, restricted to the subspace relative to quantum number n, and can be formally replaced by its diagonal matrix element (9/4)n 2 (n 2 − l 2 − l − 1); and E 1 is the exponential integral function.Equation (5) represents the contribution of weak collisions only; these are usually dominant with respect to strong collisions, that the impact approximation is valid.An estimate of the strong collision contribution to K can be done by setting formally S ≡ 0 in Equation (4) (this function strongly oscillates with respect to b and v) and adjusting the integration domain in such a way that only impact parameters smaller than the Weisskopf radius are retained.An explicit calculation yields and the final expression for the collision operator is K = K weak + K strong .A Stark line shape described within the impact approximation can be written as a sum of Lorentzian functions with a half-width at half-maximum given by the matrix elements of the collision operator.In the framework of radiative transfer Monte Carlo simulations, this renders the sampling of photon frequencies straightforward.
Figure 1 shows a plot of Lyman α at conditions of tokamak edge plasmas, calculated without Doppler broadening and using the GKS collision operator.The result is in contrasted to the result of a computer simulation, where the Schrödinger equation is solved numerically for a given set of microfield histories.These histories are simulated using particles moving in a cube with periodic boundary conditions (see [7] for details on the code).As can be seen, the impact approximation agrees well with the simulation result.Two essential features of the impact model are the assumptions of binary collision and short collision time; these features are characterized by the ratios b W /r 0 and τ c /t i , which must be much smaller than unity.Here, r 0 ~N−1/3 is the mean interparticle distance, τ c ~r0 /v 0 is the collision time, and t i is the time of interest (which is estimated as the inverse line width).Here, the values of 3 × 10 12 cm −3 and 1 eV, assumed for density and temperature, respectively, yield values of 4% and 0.6% for these ratios.binary collision and short collision time; these features are characterized by the ratios bW/r0 and τc/ti, which must be much smaller than unity.Here, r0 ~ N −1/3 is the mean interparticle distance, τc ~ r0/v0 is the collision time, and ti is the time of interest (which is estimated as the inverse line width).Here, the values of 3 × 10 12 cm −3 and 1 eV, assumed for density and temperature, respectively, yield values of 4% and 0.6% for these ratios.broadening, in conditions such that the ratios bW/r0 and τc/ti are much smaller than unity (see text).As can be seen, the impact theory agrees well with the simulation.

Influence of Zeeman Effect on Line Broadening
In tokamaks, the magnetic field is strong enough that hydrogen lines exhibit a Zeeman triplet structure.This structure is a feature of energy degeneracy removal relative to the magnetic quantum number m; it can alter the Stark broadening if the Zeeman perturbation ħωZ = |μBB| (μB being the Bohr magneton) exceeds the characteristic Stark perturbation, estimated either by n 2 ea0F0 (static picture, F0 being the Holtsmark field) or ħNπbW 2 v0ln(λD/bW) (impact limit, the logarithm stems from the behavior of the exponential integral for small argument).An extension of the GKS impact collision operator that accounts for this degeneracy removal can be devised through the use of an alternative S operator in Equation ( 4) involving the Zeeman Hamiltonian.The method employs an adaptation of the GBKO model (Griem, Baranger, Kolb, Oertel [18]), initially developed for helium lines, for hydrogen lines with Zeeman effect; details can be found in [6,19].The result is a collision operator with a dependence on the m quantum number.The following relation holds where α stands for the quantum numbers n, l, m; the constants 4 involve the atomic electron position operator (here, by convention, x ± = (x ± iy)/√2); a(s) and b(s) are the GBKO functions writable in terms of the modified Bessel functions, bα st is the strong collision radius; and bm = v/ωZ is a distance relative to Zeeman degeneracy removal.The latter mitigates the characteristic strength of the Stark perturbation.If the ordering bWα << bm is satisfied, the a function evaluated at the upper bound can be replaced by a logarithm, and the length bm plays a role similar to the Debye length.Note that the Weisskopf radius here is a function of v; the definition given in the previous section corresponds to evaluation at the thermal velocity v = v0.In Equation ( 8), the dependence on velocity has not been written explicitly for As can be seen, the impact theory agrees well with the simulation.

Influence of Zeeman Effect on Line Broadening
In tokamaks, the magnetic field is strong enough that hydrogen lines exhibit a Zeeman triplet structure.This structure is a feature of energy degeneracy removal relative to the magnetic quantum number m; it can alter the Stark broadening if the Zeeman perturbation hω Z = |µ B B| (µ B being the Bohr magneton) exceeds the characteristic Stark perturbation, estimated either by n 2 ea 0 F 0 (static picture, F 0 being the Holtsmark field) or hNπb W 2 v 0 ln(λ D /b W ) (impact limit, the logarithm stems from the behavior of the exponential integral for small argument).An extension of the GKS impact collision operator that accounts for this degeneracy removal can be devised through the use of an alternative S operator in Equation ( 4) involving the Zeeman Hamiltonian.The method employs an adaptation of the GBKO model (Griem, Baranger, Kolb, Oertel [18]), initially developed for helium lines, for hydrogen lines with Zeeman effect; details can be found in [6,19].The result is a collision operator with a dependence on the m quantum number.The following relation holds where α stands for the quantum numbers n, l, m; the constants Note that the Weisskopf radius here is a function of v; the definition given in the previous section corresponds to evaluation at the thermal velocity v = v 0 .In Equation ( 8), the dependence on velocity has not been written explicitly for the sake of simplicity.The so-called "strong collision radius"-here, b α st -coincides with the Weisskopf radius at the limit B → 0, and is different at finite magnetic field.It can be determined by the solving of an integral equation; see [6,19] for details.The mitigation of Stark broadening due to the Zeeman effect is all the more important so that the Zeeman energy Atoms 2017, 5, 36 5 of 10 sublevels are separated.This effect is particularly important on the lateral components of Lyman α because they are not affected by ∆m = 0 couplings (i.e., K α // vanishes identically in Equation ( 8)).
Figure 2 shows an illustration of the Stark broadening mitigation on the lateral blue component.The same plasma conditions as in the previous section are considered, and a magnetic field of 2 T is assumed.The adaptation of the GBKO model yields a reduction of the line width by a factor of three.A slight shift towards the blue side is also visible; this corresponds to the imaginary part of the collision operator matrix element, and its amplitude is determined by the GBKO b function.
the sake of simplicity.The so-called "strong collision radius"-here, bα st -coincides with the Weisskopf radius at the limit B → 0, and is different at finite magnetic field.It can be determined by the solving of an integral equation; see [6,19] for details.The mitigation of Stark broadening due to the Zeeman effect is all the more important so that the Zeeman energy sublevels are separated.This effect is particularly important on the lateral components of Lyman α because they are not affected by ∆m = 0 couplings (i.e., Kα // vanishes identically in Equation ( 8)). Figure 2 shows an illustration of the Stark broadening mitigation on the lateral blue component.The same plasma conditions as in the previous section are considered, and a magnetic field of 2 T is assumed.The adaptation of the GBKO model yields a reduction of the line width by a factor of three.A slight shift towards the blue side is also visible; this corresponds to the imaginary part of the collision operator matrix element, and its amplitude is determined by the GBKO b function.

An Extension of the Impact theory to Non-Binary Interactions
The impact approximation for ions can be inaccurate under high density or low temperature conditions, namely, when the dimensionless parameters <α|K|α> × τc and bW/r0 are significant.Comparisons with ab initio simulations have indicated that the use of the impact approximation can result in a systematic overestimate of the Stark width, e.g., [6].Recently, a model has been developed in such a way to account for the deviations analytically [20].The model uses kinetic theory techniques and is inspired by the "unified theory", a formalism proposed in the late sixties as an extension of the impact approximation for electrons [21,22].The main result of the unified theory is a frequencydependent collision operator, which accounts for incomplete collisions (i.e., collisions not achieved during the time of interest) and provides a correct asymptotic line shape formula.If we apply the unified theory to ions, the collision operator reads where V stands for the Stark interaction corresponding to a binary collision with an ion, N is the number of ions, the brackets {…} denote statistical average, and Q(t) is the evolution operator of the atom under the influence of one collision in the interaction picture.It obeys the Schrödinger equation and admits the Dyson series as a formal expression.Several techniques exist for the evaluation of the integral, either based on the exact solution or on the perturbation theory, as done for the impact collision operator.The behavior of the collision operator at small and large frequency detuning is obtained by asymptotic analysis.This shows that K(ω) behaves as |∆ω| −1/2 in the wings, and reduces A magnetic field of 2 T is assumed and the observation is assumed perpendicular to B (θ = 90 • )."GBKO" refers to the Griem-Baranger-Kolb-Oertel 1962 model, adapted to hydrogen in the presence of Zeeman effect; "GKS" refers to the original Griem-Kolb-Shen 1959 model for hydrogen without degeneracy removal.As can be seen, the Zeeman degeneracy removal, retained in the GBKO model, strongly reduces the line width.A slight shift is also visible.

An Extension of the Impact theory to Non-Binary Interactions
The impact approximation for ions can be inaccurate under high density or low temperature conditions, namely, when the dimensionless parameters <α|K|α> × τ c and b W /r 0 are significant.Comparisons with ab initio simulations have indicated that the use of the impact approximation can result in a systematic overestimate of the Stark width, e.g., [6].Recently, a model has been developed in such a way to account for the deviations analytically [20].The model uses kinetic theory techniques and is inspired by the "unified theory", a formalism proposed in the late sixties as an extension of the impact approximation for electrons [21,22].The main result of the unified theory is a frequency-dependent collision operator, which accounts for incomplete collisions (i.e., collisions not achieved during the time of interest) and provides a correct asymptotic line shape formula.If we apply the unified theory to ions, the collision operator reads where V stands for the Stark interaction corresponding to a binary collision with an ion, N is the number of ions, the brackets { . . .} denote statistical average, and Q(t) is the evolution operator of the atom under the influence of one collision in the interaction picture.It obeys the Schrödinger equation and admits the Dyson series as a formal expression.Several techniques exist for the evaluation of the integral, either based on the exact solution or on the perturbation theory, as done for the impact collision operator.The behavior of the collision operator at small and large frequency detuning is obtained by asymptotic analysis.This shows that K(ω) behaves as |∆ω| −1/2 in the wings, and reduces to the impact collision operator at the center.For applications in tokamak edge and divertor conditions, it has been necessary to modify the original model Equation ( 9) in such a way as to account for N body interaction effects, referred to as "correlated collisions".Such correlations occur if, during one collision, the atom "feels" the presence of the other perturbers.Correlated collisions are important when the characteristic collision frequency <α|K|α> becomes of the same order as, or larger than, the inverse correlation time of the emitter-perturber interaction potential [23,24].The latter is of the order of the inverse plasma frequency ω pi −1 (due to Debye screening).In the model reported in [20], correlated collisions are retained through a resummation procedure applied to kinetic equations of BBGKY-type.The resulting collision operator obeys a nonlinear equation with a structure similar to Equation ( 9) In the exponential, K 0 = K(ω 0 ) is the collision operator evaluated at the central frequency of the line under consideration.Its presence denotes a non-Hermitian part in the atomic Liouvillian, which can be interpreted as a renormalization or "dressing" of the atomic energy levels, induced by the presence of other perturbers during a single collision and their correlation with the collision under consideration.This model presents similarities to the result of the resonance broadening theory used for plasma turbulence [25,26].In practice, a calculation of the collision operator from Equation ( 10) should be done by iterations.Tests have indicated that such a procedure yields a fast convergence, with typically no more than ten iterations.Figure 3 presents a plot of Lyman α broadened due to ions at N e = 10 15 cm −3 , T e = T i = 1 eV, which was obtained using impact approximation, the unified theory Equation ( 9), and its extension Equation ( 10) applied to ions.No magnetic field is retained here.A numerical result from a simulation is also shown in the figure .As can be seen, the impact approximation strongly overestimates the width.This stems from the inadequacy of this model as a regime where incomplete and correlated collisions are present.The extended unified theory gives a much better result, with an overestimate of the width no larger than 20%.We can also see that the standard unified theory (i.e., not accounting for correlated collisions) is not sufficient here, and yields a different (and incorrect) shape structure, with a dip at the center.Correlated collisions are important in this region, because they govern the average atomic evolution operator long times; hence, by virtue of the Fourier transform, at small frequency detuning.The typical range for the dip corresponds to frequencies smaller than the matrix elements of K 0 .
Atoms 2017, 5, 36 6 of 10 to the impact collision operator at the center.For applications in tokamak edge and divertor conditions, it has been necessary to modify the original model Equation ( 9) in such a way as to account for N body interaction effects, referred to as "correlated collisions".Such correlations occur if, during one collision, the atom "feels" the presence of the other perturbers.Correlated collisions are important when the characteristic collision frequency <α|K|α> becomes of the same order as, or larger than, the inverse correlation time of the emitter-perturber interaction potential [23,24].The latter is of the order of the inverse plasma frequency ωpi −1 (due to Debye screening).In the model reported in [20], correlated collisions are retained through a resummation procedure applied to kinetic equations of BBGKY-type.The resulting collision operator obeys a nonlinear equation with a structure similar to Equation ( 9) In the exponential, K0 = K(ω0) is the collision operator evaluated at the central of the line under consideration.Its presence denotes a non-Hermitian part in the atomic Liouvillian, which can be interpreted as a renormalization or "dressing" of the atomic energy levels, induced by the presence of other perturbers during a single collision and their correlation with the collision under consideration.This model presents similarities to the result of the resonance broadening theory used for plasma turbulence [25,26].In practice, a calculation of the collision operator from Equation ( 10) should be done by iterations.Tests have indicated that such a procedure yields a fast convergence, with typically no more than ten iterations.Figure 3 presents a plot of Lyman α broadened due to ions at Ne = 10 15 cm −3 , Te = Ti = 1 eV, which was obtained using impact approximation, the unified theory Equation ( 9), and its extension Equation ( 10) applied to ions.No magnetic field is retained here.A numerical result from a simulation is also shown in the figure .As can be seen, the impact approximation strongly overestimates the width.This stems from the inadequacy of this model as a regime where incomplete and correlated collisions are present.The extended unified theory gives a much better result, with an overestimate of the width no larger than 20%.We can also see that the standard unified theory (i.e., not accounting for correlated collisions) is not sufficient here, and yields a different (and incorrect) shape structure, with a dip at the center.Correlated collisions are important in this region, because they govern the average atomic evolution operator at long times; hence, by virtue of the Fourier transform, at small frequency detuning.The typical range for the dip corresponds to frequencies smaller than the matrix elements of K0.Profile of the Lyman α line at conditions such that neither the impact approximation nor the unified theory (UT) for ions are valid.The extension of the unified theory (UT), which accounts for correlated collisions, provides a good estimate of the line width.

Stark-Zeeman Line Shapes in Magnetic White Dwarfs
About 10% of white dwarfs are known to have a magnetic field strength of 10 5 to 5 × 10 8 G (10 to 5 × 10 4 T), as indicated by spectroscopic observations and models [9,10,27].An interpretation of the Figure 3. Profile of the Lyman α line at conditions such that neither the impact approximation nor the unified theory (UT) for ions are valid.The extension of the unified theory (UT), which accounts for correlated collisions, provides a good estimate of the line width.

Stark-Zeeman Line Shapes in Magnetic White Dwarfs
About 10% of white dwarfs are known to have a magnetic field strength of 10 5 to 5 × 10 8 G (10 to 5 × 10 4 T), as indicated by spectroscopic observations and models [9,10,27].An interpretation of the shape of absorption lines requires the Zeeman effect be accounted for in line broadening Atoms 2017, 5, 36 7 of 10 models, as done in tokamak edge plasma spectroscopy.Work is currently ongoing in order to improve already-existing databases by implementing the Zeeman effect in Stark line shape models.Figure 4 shows an illustration of the importance of the Zeeman effect in white dwarf spectra.A profile of Hα has been calculated with a numerical simulation assuming N e = 10 17 cm −3 , T e = T i = 0.5 eV, B = 500 T and observation perpendicular to B. These conditions are characteristic of the atmosphere of a white dwarf of DAH (magnetic DA) type, presenting hydrogen absorption lines.The spectrum exhibits a triplet structure characteristic to the Zeeman effect.The description of Stark broadening can be more challenging than in tokamak devices because particle correlations are important (here r 0 /λ D is of the order of unity) and the quasiparticle model used in the simulation may not be accurate enough.An extension of this work would consist of using full molecular dynamic simulations.Features specific to dense plasmas, such as quadrupole interactions, could also be important, and should be addressed in detail; see [28] for recent works.The presence of a strong magnetic field changes the plasma dynamics through the gyromotion effect and this can result in an alteration of electron collision operator formulas [29,30].Specific modifications to the atomic energy level structure (quadratic Zeeman effect, Lorentz electric field) could also be important [31].
Atoms 2017, 5, 36 7 of 10 shape of absorption lines requires the Zeeman effect be accounted for in line broadening models, as done in tokamak edge plasma spectroscopy.Work is currently ongoing in order to improve alreadyexisting databases by implementing the Zeeman effect in Stark line shape models.Figure 4 shows an illustration of the importance of the Zeeman effect in white dwarf spectra.A profile of Hα has been calculated with a numerical simulation assuming Ne = 10 17 cm −3 , Te = Ti = 0.5 eV, B = 500 T and observation perpendicular to B. These conditions are characteristic of the atmosphere of a white dwarf of DAH (magnetic DA) type, presenting hydrogen absorption lines.The spectrum exhibits a triplet structure characteristic to the Zeeman effect.The description of Stark broadening can be more challenging than in tokamak devices because particle correlations are important (here r0/λD is of the order of unity) and the quasiparticle model used in the simulation may not be accurate enough.An extension of this work would consist of using full molecular dynamic simulations.Features specific to dense plasmas, such as quadrupole interactions, could also be important, and should be addressed in detail; see [28] for recent works.The presence of a strong magnetic field changes the plasma dynamics through the gyromotion effect and this can result in an alteration of electron collision operator formulas [29,30].Specific modifications to the atomic energy level structure (quadratic Zeeman effect, Lorentz electric field) could also be important [31].

Non-Hydrogen Species
Tokamak plasmas contain a sufficient amount of impurities, due to plasma-wall interactions or vacuum chamber conditioning, to be observable on spectra (e.g., oxygen, nitrogen, carbon [32]).Specific atomic species of high Z can be found in some, but not all, tokamaks (e.g., tungsten [33]).Helium is also expected to be present in high amounts in ITER during the nuclear phase, where discharges with tritium will be done.As a rule, an analysis of the position and the intensity of impurity lines provides information on the plasma composition, the recycling and sputtering processes, which helps in the setting up of tokamak discharge scenarios.The design of fast routines employs line shape codes (e.g., PPP [34]), and can also use databases available in astrophysics such as STARK-B (stark-b.obspm.fr).An issue that is specific to non-hydrogen lines concerns the modeling of collisional electron broadening, accounting for degeneracy removal with quantum number l.This degeneracy removal yields a reduction of the Stark effect and can induce a shift on the line under consideration.An illustration of this effect is shown in Figure 5.The He 492 nm line (1s4d 1 D-1s2p 1 P) has been addressed, assuming Ne = 10 15 cm −3 and Te = 10,000 K, without ion broadening and without magnetic field for the sake of clarity.Two Lorentzian curves are displayed in the graph.One corresponds to the result of a hydrogen approximation (through the GKS impact model discussed above), and the other one corresponds to the model reported in the STARK-B database, which accounts for the degeneracy removal.As can be seen, the line width is sensitive to the model used.The shift visible in the STARK-B result is a feature of the degeneracy removal.It depends on the plasma parameters,

Non-Hydrogen Species
Tokamak plasmas contain a sufficient amount of impurities, due to plasma-wall interactions or vacuum chamber conditioning, to be observable on spectra (e.g., oxygen, nitrogen, carbon [32]).Specific atomic species of high Z can be found in some, but not all, tokamaks (e.g., tungsten [33]).Helium is also expected to be present in high amounts in ITER during the nuclear phase, where discharges with tritium will be done.As a rule, an analysis of the position and the intensity of impurity lines provides information on the plasma composition, the recycling and sputtering processes, which helps in the setting up of tokamak discharge scenarios.The design of fast routines employs line shape codes (e.g., PPP [34]), and can also use databases available in astrophysics such as STARK-B (stark-b.obspm.fr).An issue that is specific to non-hydrogen lines concerns the modeling of collisional electron broadening, accounting for degeneracy removal with quantum number l.This degeneracy removal yields a reduction of the Stark effect and can induce a shift on the line under consideration.An illustration of this effect is shown in Figure 5.The He 492 nm line (1s4d 1 D-1s2p 1 P) has been addressed, assuming N e = 10 15 cm −3 and T e = 10,000 K, without ion broadening and without magnetic field for the sake of clarity.Two Lorentzian curves are displayed in the graph.One corresponds to the result of a hydrogen approximation (through the GKS impact model discussed above), and the other one corresponds to the model reported in the STARK-B database, which accounts for the degeneracy removal.As can be seen, the line width is sensitive to the model used.The shift visible in the STARK-B result is a feature of the degeneracy removal.It depends on the plasma parameters, and can also be used for diagnostic purposes.An extension of already-existing databases to magnetized plasma conditions, where the Zeeman effect is visible in spectra, is under consideration.Work is presently ongoing for applications both in magnetic fusion and astrophysical plasmas.
Atoms 2017, 5, 36 8 of 10 and can also be used for diagnostic purposes.An extension of already-existing databases to magnetized plasma conditions, where the Zeeman effect is visible in spectra, is under consideration.Work is presently ongoing for applications both in magnetic fusion and astrophysical plasmas.Conditions relevant to high-density tokamak divertor plasmas are considered here, and no magnetic field is retained for the sake of clarity.As can be seen, the result significantly depends on the model used.

Conclusions
The shape of atomic spectral lines in plasmas contains information on the plasma parameters, and can be used as a diagnostic tool.In this work, we have examined a selection of problems involved in the modeling of Stark broadening for applications in magnetic fusion and in astrophysics.A specific issue concerns the description of ion dynamics on hydrogen lines, which is important if the upper principal quantum number is low as on Balmer α.The use of a collision operator for ions is convenient because the line shape can be described analytically.In this framework, we have shown that the early model by Griem et al. (GBKO) for atoms with nondegenerate levels can be adapted to hydrogen in the presence of the Zeeman effect.It coincides with numerical simulations at the socalled impact limit, i.e., when the microfield fluctuates at a time scale much smaller than the time of interest, and when strong collisions are rare events.We have also shown that deviations in the impact limit due to non-binary interactions can be retained through an adaptation of the so-called unified theory, a formalism based on kinetic equations and related statistical concepts (BBGKY hierarchy).In astrophysics, magnetic white dwarfs present similarities with tokamak edge plasmas because the hydrogen absorption lines in atmospheres can exhibit a Zeeman triplet structure.We have illustrated this point through calculation of an absorption line in conditions relevant to a DAH stellar atmosphere.In this framework, an extension of already-existing Stark databases that is devoted to accounting for the Zeeman effect is presently under preparation.We have also shown that the description of the broadening of non-hydrogen lines is a challenging issue, both in astrophysics and in magnetic fusion.Modeling work is also ongoing in this framework in order to improve the databases.Plot of the electron impact broadening of He 1s4d 1 D-1s2p 1 P calculated using the STARK-B database (solid line) and a hydrogen approximation for the energy level structure (dashed line).Conditions relevant to high-density tokamak divertor plasmas are considered here, and no magnetic field is retained for the sake of clarity.As can be seen, the result significantly depends on the model used.

Conclusions
The shape of atomic spectral lines in plasmas contains information on the plasma parameters, and can be used as a diagnostic tool.In this work, we have examined a selection of problems involved in the modeling of Stark broadening for applications in magnetic fusion and in astrophysics.A specific issue concerns the description of ion dynamics on hydrogen lines, which is important if the upper principal quantum number is low as on Balmer α.The use of a collision operator for ions is convenient because the line shape can be described analytically.In this framework, we have shown that the early model by Griem et al. (GBKO) for atoms with nondegenerate levels can be adapted to hydrogen in the presence of the Zeeman effect.It coincides with numerical simulations at the so-called impact limit, i.e., when the microfield fluctuates at a time scale much smaller than the time of interest, and when strong collisions are rare events.We have also shown that deviations in the impact limit due to non-binary interactions can be retained through an adaptation of the so-called unified theory, a formalism based on kinetic equations and related statistical concepts (BBGKY hierarchy).In astrophysics, magnetic white dwarfs present similarities with tokamak edge plasmas because the hydrogen absorption lines in atmospheres can exhibit a Zeeman triplet structure.We have illustrated this point through calculation of an absorption line in conditions relevant to a DAH stellar atmosphere.In this framework, an extension of already-existing Stark databases that is devoted to accounting for the Zeeman effect is presently under preparation.We have also shown that the description of the broadening of non-hydrogen lines is a challenging issue, both in astrophysics and in magnetic fusion.Modeling work is also ongoing in this framework in order to improve the databases.

Figure 1 .
Figure 1.Deuterium Lyman α line profile in the atomic frame of reference, i.e., without Doppler

Figure 1 .
Figure1.Deuterium Lyman α line profile in the atomic frame of reference, i.e., without Doppler broadening, in conditions such that the ratios b W /r 0 and τ c /t i are much smaller than unity (see text).As can be seen, the impact theory agrees well with the simulation.
the atomic electron position operator (here, by convention, x ± = (x ± iy)/ √ 2); a(s) and b(s) are the GBKO functions writable in terms of the modified Bessel functions, b α st is the strong collision radius; and b m = v/ω Z is a distance relative to Zeeman degeneracy removal.The latter mitigates the characteristic strength of the Stark perturbation.If the ordering b Wα << b m is satisfied, the a function evaluated at the upper bound can be replaced by a logarithm, and the length b m plays a role similar to the Debye length.

Figure 2 .
Figure 2.Lateral blue component of the Lyman α Zeeman triplet in the atomic frame of reference.A magnetic field of 2 T is assumed and the observation is assumed perpendicular to B (θ = 90°)."GBKO" refers to the Griem-Baranger-Kolb-Oertel 1962 model, adapted to hydrogen in the presence of Zeeman effect; "GKS" refers to the original Griem-Kolb-Shen 1959 model for hydrogen without degeneracy removal.As can be seen, the Zeeman degeneracy removal, retained in the GBKO model, strongly reduces the line width.A slight shift is also visible.

Figure 2 .
Figure 2. Lateral blue component of the Lyman α Zeeman triplet in the atomic frame of reference.A magnetic field of 2 T is assumed and the observation is assumed perpendicular to B (θ = 90 • )."GBKO" refers to the Griem-Baranger-Kolb-Oertel 1962 model, adapted to hydrogen in the presence of Zeeman effect; "GKS" refers to the original Griem-Kolb-Shen 1959 model for hydrogen without degeneracy removal.As can be seen, the Zeeman degeneracy removal, retained in the GBKO model, strongly reduces the line width.A slight shift is also visible.

Figure 3 .
Figure 3. Profile of the Lyman α line at conditions such that neither the impact approximation nor the unified theory (UT) for ions are valid.The extension of the unified theory (UT), which accounts for correlated collisions, provides a good estimate of the line width.

Figure 4 .
Figure 4. Balmer α line shape in conditions relevant to DAH white dwarf atmospheres.The spectrum exhibits a triplet structure due to the strong magnetic field.

Figure 4 .
Figure 4. Balmer α line shape in conditions relevant to DAH white dwarf atmospheres.The spectrum exhibits a triplet structure due to the strong magnetic field.

Figure 5 .
Figure 5. Plot of the electron impact broadening of He 1s4d 1 D-1s2p 1 P calculated using the STARK-B database (solid line) and a hydrogen approximation for the energy level structure (dashed line).Conditions relevant to high-density tokamak divertor plasmas are considered here, and no magnetic field is retained for the sake of clarity.As can be seen, the result significantly depends on the model used.

Figure 5 .
Figure 5. Plot of the electron impact broadening of He 1s4d 1 D-1s2p 1 P calculated using the STARK-B database (solid line) and a hydrogen approximation for the energy level structure (dashed line).Conditions relevant to high-density tokamak divertor plasmas are considered here, and no magnetic field is retained for the sake of clarity.As can be seen, the result significantly depends on the model used.