#### Abstract

Stark broadening of hydrogen lines is investigated in low-density magnetized plasmas, at typical conditions of magnetic fusion experiments. The role of time ordering is assessed numerically, by using a simulation code accounting for the evolution of the microscopic electric field generated by the charged particles moving at the vicinity of the atom. The Zeeman effect due to the magnetic field is also retained. Lyman lines with a low principal quantum number are first investigated, for an application to opacity calculations; next Balmer lines with successively low and high principal quantum numbers are considered for diagnostic purposes. It is shown that neglecting time ordering results in a dramatic underestimation of the Stark effect on the low- lines. Another conclusion is that time ordering becomes negligible only when ion dynamics effects vanish, as shown in the case of high- lines.

#### 1. Introduction

In magnetic fusion, detailed line shapes are of interest for accurate diagnostics or radiative transfer simulations. For plasma conditions and magnetic fields encountered in the divertor of present and future tokamaks, an accurate model for the line shape of the hydrogen isotopes should include Zeeman and Stark effects, and retain the dynamics of the ion-emitter interaction. Since we then have to solve a quantum time-dependent problem, understanding the role of time ordering becomes an important issue both from the fundamental and computational points of view (note, this problem is also investigated in other contexts, e.g., [1–3]). Time ordering has already been studied in the Stark broadening literature, but generally for the electron broadening [4–7]. Our aim here is to investigate the role of time ordering for the ion perturbation on hydrogen lines for plasmas with temperature in the eV range, and densities of about , conditions which are expected in the divertor of the future ITER tokamak. We recall in Section 2 the basic formalism used for line shape calculations in the presence of Stark and Zeeman effects, and briefly introduce the issue of time ordering. Line shapes in the atom’s frame of reference are considered, that is, in the Doppler free case. We present in Section 3 an ab initio simulation technique able to provide accurate line shapes including all the effects of time ordering. Calculations of hydrogen line shapes of Lyman and Balmer series are presented in Section 4, with and without the effect of time ordering, and compared to calculations performed in the static ion limit. The role of time ordering and the issue of retaining it in a line shape calculation are discussed in the conclusion for lines with low and high principal quantum number of the upper state.

#### 2. Formalism

According to classical textbooks or review articles (e.g., [8–10]), a line shape in the atom’s rest frame at a frequency , is given by the Fourier transform of the dipole autocorrelation function : Here, and denote the principal quantum number with respect to the initial and final states, respectively; is the dipole operator projected into the polarization plane; (resp., ) is the evolution operator projected into the subspace related to the initial (resp., final) states; is the density operator related to the initial states; the brackets stand for a statistical average over the perturbers’ states and the trace is performed over the atomic states. The evolution operator obeys the time-dependent Schrödinger equation where is the dipole projected into the subspace related to is the atomic Hamiltonian accounting for the Zeeman effect, that is, with being, respectively, the magnetic moment and the magnetic field; and is the time-dependent electric field created at the atom’s location by both electrons () and ions (). In our conditions of interest the electron contribution can be safely described by an impact collision operator . We will therefore formally replace by . We will also use the interaction representation, that is, consider the evolution operator which obeys the following equation where is the dipole operator in the interaction representation.

The Schrödinger equation in the interaction representation (3), associated with the initial condition , admits as a formal solution the so-called Dyson series expansion where . Following Baranger [11, 12], we write this solution a time-ordered exponential of a matrix Here, is the time ordering operator, defined by its action on a product of time-dependent operators as follows: Because of the presence of the -operator in (5), the evolution operator cannot in general be described by a pure exponential of a matrix since the -operator does not commute with itself at different times. The time ordering operator can only be dropped out in the limiting case where the electric field does not evolve significantly during the time of interest (quasistatic limit). In the following we will address the general case by using a simulation code.

#### 3. Assessment of Time Ordering by Computer Simulations

##### 3.1. The Numerical Simulation Method

The purpose of ab initio simulations is to numerically reproduce the motion of the charged particles in the plasma so as to obtain the time-dependent electric microfield . Essentially, a numerical simulation consists of (i) the calculation of a set of realizations for the electric field; (ii) the numerical integration of the Schrödinger equation for each realization; (iii) the average of the evolution operators on a set of realizations and the Fourier transform of the autocorrelation function. In the simulations performed for this work, we use a code [13] developed according to the method reported in [14]. We consider that the ions move along straight line trajectories with constant velocities, sampled among the particles according to an equilibrium Maxwell distribution function. The electrons (whose contribution to line broadening is mostly negligible for our conditions of interest) are not simulated here, but are described with an impact collision operator (see previous section). The treatment of the correlations between ions and electrons is retained by using Debye screened fields. The largeness of the ratio between the ion Larmor radius and the Debye length (of about 40 at T, i.e., typical conditions expected in the ITER divertor) ensures the assumption of straight line trajectories. A cubic cell with periodic boundary conditions is considered. For each history of the electric field, the code solves the time-dependent Schrödinger equation for the evolution operator according to the algorithm , with being the infinitesimal evolution operator between times and . The latter operator is not affected by time ordering if the time step *∆*t is sufficiently small, and hence can be evaluated by a matrix exponential
In the code, we use this property and calculate the exponential according to the scaling squaring method (e.g., [15]). This method is also used for the calculation of the expression of not accounting for time ordering, given by the exponential of .

##### 3.2. Line Shape Calculations

Time ordering should play a role on lines which are affected by ion dynamics, that is, with a low upper principal quantum number . To illustrate this point, we have successively calculated (whose shapes are involved in the calculations of divertor plasma opacity), and (which is commonly observed in tokamak experiments and routinely used for diagnostic purposes). Figures 1(a), 1(b), and 1(c), show, for each line profile, the result obtained with simulations retaining time ordering or not, at T, and in perpendicular observation. The profiles resulting from the quasistatic approximation are also plotted. In each case, the Zeeman effect is important and leads either to a Lorentz triplet structure () or to an additional broadening (). As can be seen, neglecting time ordering leads to a dramatic underestimation of the Stark effect. The Zeeman components of and are much narrower than those obtained both with the exact solution and the quasistatic approximation (note, by a factor of ~3-4 in the case of the lateral components of ), and the Zeeman-Stark pattern becomes visible on when time ordering is neglected. Another gain of structure is provided by the apparition of Stark components, in particular on where a splitting of the lateral Zeeman components due to Stark effect is clearly visible.

**(a)**

**(b)**

**(c)**

The underestimation of the Stark effect can be explained by noting that, in the solution of Schrödinger’s equation neglecting time ordering, the matrix exponential essentially involves the time average of the electric field . This quantity vanishes on a time scale smaller than the time of interest, a quantity usually defined as the inverse of the line width. Therefore, except for very large frequencies (which correspond to short times), the emitter-ion interaction potential is very small and the resulting line shape becomes very narrow. The line shape is even narrower than that obtained using the quasistatic approximation because the cancellation of the average electric field, equivalent to a cancellation of the evolution operator oscillations, leads to a slower decorrelation of the atomic dipole. The deviation should be smaller for lines with a higher upper principal quantum number , because ion dynamics is less important in this case. This is indeed illustrated in Figure 2, where a plot of (a) , (b) , and (c) are presented. For these lines, the ratio between the time of interest at half maximum and the collision time is 14%, 9%, and 0.2%, respectively. As shown, the role of time ordering becomes less and less important as this ratio decreases, that is, as the ions become static.

**(a)**

**(b)**

**(c)**

#### 4. Conclusion

We have addressed the role of time ordering on hydrogen Zeeman-Stark profiles in low-density plasmas, for typical conditions of tokamak divertors. With numerical simulations, we have shown that neglecting time ordering on lines with a low upper principal quantum number leads to strong deviations, with a systematic underestimate of the Stark width of the Zeeman components. This is interpreted in terms of the time averaged electric field, namely, the latter rapidly vanishes during the decorrelation of the atomic dipole so that the resulting effective Stark effect is reduced. Conversely, we have shown that the deviations are weak on lines with a higher upper quantum number, merely because they are much less affected by ion dynamics. This result, of interest for spectroscopy of magnetic fusion experiments, shows that: (i) the development of line shape models including ion dynamics for Monte Carlo investigations of radiative transfer (e.g., [16]) requires to account for time ordering; (ii) in a similar way, time ordering must be accounted for in the Stark broadening models used for passive spectroscopy diagnostics; (iii) obviously, the effect of time ordering disappears as one approaches the validity conditions of the static ion approximation, and this is clearly the case for high- lines in our divertor conditions. A possible extension of our work would consist of a parameterization of the role of time ordering on low- lines, for diagnostic as well as for radiative transfer calculation purposes.

#### Acknowledgments

This work is supported by an EFDA Fellowship contract, by the French Federation on Magnetic Fusion Research (project “Radiation Absorption Effects”), by the French Research National Agency (Project “PHOTONITER”, Contract ANR-07-BLAN-0187-01), and by the collaboration PIIM/CEA Cadarache (Contract LRC DSM 99-14).