Next Article in Journal
An Investigation of Laser Produced Lead-Tin Alloy Plasmas between 10 and 18 nm
Next Article in Special Issue
On the Stark Effect of the O I 777-nm Triplet in Plasma and Laser Fields
Previous Article in Journal
Calcium and Potassium Spectra in the EUV
Previous Article in Special Issue
Theoretical Stark Broadening Parameters for UV–Blue Spectral Lines of Neutral Vanadium in the Solar and Metal-Poor Star HD 84937 Spectra
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogen Line Shapes in Plasmas with Large Magnetic Fields

Laboratoire PIIM, UMR 7345 Aix-Marseille Université/CNRS, Case 232, F-13397 Marseille CEDEX 20, France
Submission received: 19 September 2020 / Revised: 6 October 2020 / Accepted: 12 October 2020 / Published: 15 October 2020

Abstract

:
We report on hydrogen line shape calculations in the presence of an external magnetic field, at conditions such that the quadratic Zeeman effect is important. The latter is described through a term proportional to B2 in the Hamiltonian, accounting for atomic diamagnetism. It provides a shift and an asymmetry on Lorentz triplets, and it leads to the occurrence of forbidden components. Motivated by investigations performed at the fifth edition of the Spectral Line Shape in Plasmas (SLSP5) code comparison workshop, we perform new calculations of hydrogen Lyman line profiles. Field values representative of magnetized white dwarf atmosphere conditions are taken. The calculations are done using a computer simulation technique, designed for Stark broadening modeling. A discussion of the results is done in the framework of plasma diagnostics.

1. Introduction

A significant number of white dwarfs contain a high magnetic field, with values up to 10 kT [1,2]. For such high values, the atoms emitting and absorbing radiation are subject to the quadratic Zeeman effect, which yields a perturbation of the energy levels proportional to B2 in addition to the linear Zeeman perturbation. A feature of the quadratic Zeeman effect is that the Lorentz triplets observed in spectra are strongly shifted and are no longer symmetric; furthermore, the spectral lines exhibit additional components corresponding to forbidden transitions. The quadratic Zeeman effect has already been investigated in the literature, in particular for hydrogen; e.g., see [3] for a comprehensive review. Recent works have been devoted to the modeling of the combined action of the quadratic Zeeman effect and the Stark effect associated with the plasma microfield, for applications to the diagnostic of magnetized white dwarf atmospheres [4,5]. In 2019, this topic was selected as a case problem at the Spectral Line Shape in Plasmas code comparison workshop (SLSP [6], fifth edition). We here report on new line shape calculations performed in the conditions close to those considered at the workshop. The first hydrogen Lyman lines are considered.

2. Line Broadening Modeling

We first give a summary of the formalism involved in line shape calculations. Consider an atom that emits or absorbs a photon (for the sake of clarity in discussions, it will always be referred to as the “emitter”) in the presence of charged particles located at its vicinity. The electric field they generate modifies the energy level structure of the emitter (Stark effect), and the resulting spectral lines are broader. The spectrum of a Stark-broadened line is proportional to the Fourier transform of the dipole autocorrelation function [7]:
I ( ω ) = 1 π Re 0 d t { T r [ ρ d ( 0 ) · d ( t ) ] } e i ω t .
Here, d(t) denotes the atomic dipole operator expressed in the Heisenberg picture, ρ is the density operator, Tr(...) denotes a trace over the atomic states and the brackets {...} correspond to an average over the perturbers. The dipole at time t is related to its initial value through d ( t ) =   U + ( t ) d ( 0 ) U ( t ) , where U ( t )     U ( t , 0 ) is the evolution operator from time 0 to t. The latter obeys the time-dependent Schrödinger equation:
i d U d t ( t ) = [ H 0 + V ( t ) ] U ( t ) .
Here, H0 is the Hamiltonian corresponding to the unperturbed atom and V ( t ) = d · E ( t ) is the Stark effect term; E(t) is the electric field generated at the emitter’s location by the perturbers. In the framework of the so-called standard model, the ions are assumed motionless and the electrons are described within a collisional picture in terms of a relaxation operator Ke. In practice, these two assumptions yield the substitution d · E ( t ) d · E i ( 0 ) i K e in the Stark effect term; the Schrödinger equation has an exponential solution and the line shape is obtained from a matrix inversion (e.g., [7] for details). The standard model provides a good estimate of the line broadening but it can be inaccurate in regimes such that the ions move significantly during the characteristic decorrelation time for the dipole (usually referred to as the “time-of-interest”; it is of the order of the inverse line width). More elaborate techniques using dedicated models and codes can be used (e.g., see reports of previous SLSP workshops [8,9]). In the following, we report on computer simulations performed using a code developed previously for tokamak plasma spectroscopy applications [10,11]. The Zeeman effect is retained by adding the term μ · B (with μ being the magnetic moment) to H0. An illustration is provided in Figure 1. The Lyman α line shape has been calculated in high-density tokamak divertor plasma conditions [12]. A magnetic field of 8 T has been assumed and the observation direction has been set perpendicular to B. The spectrum in the atomic frame of reference is shown, i.e., Doppler broadening due to the thermal motion of emitters is not retained. For the sake of clarity, only the contribution due to ions has been included here in the Stark broadening calculation. A density of 1015 cm−3 and a temperature of 1 eV have been assumed for the ions. The triplet structure visible in the spectrum is a feature of the Zeeman effect. The separation between components is proportional to B. The central component (π) is polarized parallel to the magnetic field and the lateral components (σ) are polarized perpendicular to the magnetic field. In the next section, we report on spectra calculations in magnetic white dwarf conditions. A focus is put on the case of large magnetic fields, such that the Zeeman splitting becomes a function nonlinear in B.

3. Investigation of Quadratic Zeeman Effect

A significant amount of white dwarf spectra exhibit lines with an asymmetric triplet structure, which is a feature of the presence of strong magnetic fields [1,2]. Simple estimates based on the separation between the components indicate that the magnetic field can attain several kiloteslas. The asymmetry stems from the contribution of the quadratic Zeeman effect to the energy level structure. This effect is retained by adding the term e2(B × r)2/8me (atomic diamagnetism) to the Hamiltonian H0 [3]. An estimate based on perturbation theory indicates that the Zeeman π and σ components are not equidistant and are shifted. The energy level correction due to the quadratic Zeeman effect is given by the following diagonal matrix element
n l m | e 2 8 m e ( B × r ) 2 | n l m = e 2 B 2 a 0 2 8 m e n 2 [ 5 n 2 + 1 3 l ( l + 1 ) ] ( l 2 + l 1 + m 2 ) ( 2 l 1 ) ( 2 l + 3 ) .
Applying this relation to the π and σ components, one obtains, for Lyman lines,
Δ E n π = e 2 B 2 a 0 2 8 m e [ n 4 n 2 2 ] ,
Δ E n σ ± = ± e B 2 m e + e 2 B 2 a 0 2 4 m e [ n 4 n 2 1 ] ,
where Δ E n π and Δ E n σ ± refer to the position of the π and σ± components, respectively, with respect to the energy level En of the unperturbed atom. These relations serve as an indicator of global shift (directly through Equation (4)) and asymmetry (by noting that the energy differences Δ E n σ + Δ E n π and Δ E n π Δ E n σ are different). The diagonal matrix element (3) does not provide a complete description of the quadratic Zeeman effect because the operator (B × r)2 involves couplings between states with different quantum numbers n and l. A consequence is the occurrence of forbidden components, violating the |Δl| = 1 selection rule. Figure 2 gives an illustration of this property. The π component of the hydrogen Lyman γ line shape has been calculated assuming a magnetic field of 5 kT. Couplings between different n have been neglected but couplings between different l have been retained. The resulting spectrum exhibits two components, which correspond to the 4p0-1s transition (allowed) and the 4f0-1s transition (forbidden). A homogeneous Lorentz broadening has been assumed for each component. The numerical value for the width has been set up in order to get clean components; in particular, it does not rely on any physical mechanism here. Furthermore, for the sake of clarity, the overlap with neighboring lines has not been retained in the calculation.
The quadratic Zeeman effect is currently the subject of investigations for the design of accurate line shape models. A part of the fifth edition of the SLSP workshop was dedicated to this issue. A cross-check between line shape codes simultaneously accounting for Stark broadening and the Zeeman effect has been performed. Two test cases have been defined in the following way: (i) H Lyman α line with Ne = 1019 cm−3, T = 5 eV, B = 0, 0.5, 2, 5 kT; (ii) H Lyman β with the same parameter set. A pure hydrogen plasma is considered in each case, and no couplings between different n are considered. We show in Figure 3 an example of the result. The π component of Lyman β, obtained assuming Stark broadening due to ions only, is plotted. The two curves visible in the figure correspond to one calculation done retaining the linear Zeeman effect only and one calculation done retaining both linear and quadratic Zeeman effects. The overlap with neighboring lines has not been retained, again for the sake of clarity. As can be seen in the figure, the spectrum is split in many components when the quadratic Zeeman is retained. These components are a consequence of the simultaneous action of the magnetic field and the plasma’s microscopic electric field on the atomic energy levels. Presently, no analytical formula is applicable to the modeling of such a spectrum in a general case. New calculations are presently underway in order to have a better understanding of this structure. They will involve both codes and models based on specific approximations (such as the standard model discussed above).

4. Conclusions

New hydrogen line shape calculations have been performed with the assumption that a high magnetic field is present. A focus has been put on the quadratic Zeeman effect. This effect is accounted for through an additional term, proportional to B2, in the atomic Hamiltonian. It provides a shift and an asymmetry on Lorentz triplets, and it leads to the occurrence of forbidden components. We have illustrated these properties through calculations of Lyman lines in ideal cases. In the framework of plasma spectroscopy modeling, we have shown that the simultaneous action of Stark and Zeeman effects results in spectra presenting a complex structure, with many components. New calculations are presently underway in order to have a better understanding of this issue. Applications to the diagnostic of magnetic white dwarf atmospheres will be performed. An additional issue that will be addressed is the modification of the particle trajectories due to the magnetic field, and how this alters spectral line shapes. New investigations will be done along the lines of previous works [13,14,15].

Funding

This research received no external funding.

Acknowledgments

This work was motivated by the 5th edition of the Spectral Line Shape in Plasmas (SLSP5) code comparison workshop. I thank all participants for helpful discussions during the workshop.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. 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, 1341–1350. [Google Scholar] [CrossRef]
  2. Kepler, S.O.; Pelisoli, I.; Jordan, S.; Kleinman, S.J.; Koester, D.; Külebi, B.; Peçanha, V.; Castanheira, B.G.; Nitta, A.; Costa, J.E.S.; et al. Magnetic white dwarf stars in the Sloan Digital Sky Survey. Mon. Not. R. Astron. Soc. 2013, 429, 2934–2944. [Google Scholar] [CrossRef]
  3. Garstang, R.H. Atoms in high magnetic fields. Rep. Prog. Phys. 1977, 40, 105–154. [Google Scholar] [CrossRef]
  4. Kieu, N.; Rosato, J.; Stamm, R.; Kovačević-Dojcinović, J.; Dimitrijević, M.S.; Popović, L.C.; Simić, Z. A New Analysis of Stark and Zeeman Effects on Hydrogen Lines in Magnetized DA White Dwarfs. Atoms 2017, 5, 44. [Google Scholar] [CrossRef] [Green Version]
  5. Rosato, J.; Hannachi, I.; Stamm, R. Hydrogen Stark broadening calculations in white dwarf atmosphere conditions. Contrib. Astron. Obs. Skalnaté Pleso 2020, 50, 105–110. [Google Scholar] [CrossRef]
  6. Spectral Line Shapes in Plasmas Workshops. Available online: http://plasma-gate.weizmann.ac.il/slsp/ (accessed on 19 September 2020).
  7. Griem, H.R. Spectral Line Broadening by Plasmas; Academic Press: London, UK, 1974. [Google Scholar]
  8. Stambulchik, E.; Calisti, A.; Chung, H.-K.; González, M.Á. Special Issue on Spectral Line Shapes in Plasmas. Atoms 2014, 2, 378–381. [Google Scholar] [CrossRef] [Green Version]
  9. Stambulchik, E.; Calisti, A.; Chung, H.-K.; González, M.Á. Special Issue on Spectral Line Shapes in Plasmas II. Atoms 2019, 7, 20. [Google Scholar] [CrossRef] [Green Version]
  10. Rosato, J.; Marandet, Y.; Capes, H.; Ferri, S.; Mossé, C.; Godbert-Mouret, L.; Koubiti, M.; Stamm, R. Stark broadening of hydrogen lines in low-density magnetized plasmas. Phys. Rev. E 2009, 79, 046408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Rosato, J.; Marandet, Y.; Stamm, R. Quantifying the statistical noise in computer simulations of Stark broadening. J. Quant. Spectrosc. Radiat. Transfer 2020, 249, 107002. [Google Scholar] [CrossRef]
  12. Welch, B.L.; Griem, H.R.; Terry, J.; Kurz, C.; LaBombard, B.; Lipschultz, B.; Marmar, E.; McCracken, G. Density measurements in the edge, divertor and X-point regions of Alcator C-Mod from Balmer series emission. Phys. Plasmas 1995, 2, 4246–4251. [Google Scholar] [CrossRef]
  13. Rosato, J.; Ferri, S.; Stamm, R. Influence of Helical Trajectories of Perturbers on Stark Line Shapes in Magnetized Plasmas. Atoms 2018, 6, 12. [Google Scholar] [CrossRef] [Green Version]
  14. Alexiou, S. Line Shapes in a Magnetic Field: Trajectory Modifications I: Electrons. Atoms 2019, 7, 52. [Google Scholar] [CrossRef] [Green Version]
  15. Alexiou, S. Line Shapes in a Magnetic Field: Trajectory Modifications II: Full Collision-Time Statistics. Atoms 2019, 7, 94. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Plot of hydrogen Lyman α line shape calculated in dense tokamak divertor plasma conditions. The spectrum in the atomic frame of reference is shown here, i.e., Doppler broadening due to thermal motion is not retained. See text for details on the calculation. The separation between the components is proportional to the magnetic field.
Figure 1. Plot of hydrogen Lyman α line shape calculated in dense tokamak divertor plasma conditions. The spectrum in the atomic frame of reference is shown here, i.e., Doppler broadening due to thermal motion is not retained. See text for details on the calculation. The separation between the components is proportional to the magnetic field.
Atoms 08 00074 g001
Figure 2. A feature of the quadratic Zeeman effect is the occurrence of forbidden components, namely such that the selection rule |Δl| = 1 is not satisfied. Here, this effect is illustrated on the π component of the Lyman γ line shape. The two components correspond to the 4p0-1s transition (allowed) and the 4f0-1s transition (forbidden). Other physical mechanisms of line broadening have not been retained and overlapping with neighboring lines has not been included for the sake of clarity. See text for details on the calculation.
Figure 2. A feature of the quadratic Zeeman effect is the occurrence of forbidden components, namely such that the selection rule |Δl| = 1 is not satisfied. Here, this effect is illustrated on the π component of the Lyman γ line shape. The two components correspond to the 4p0-1s transition (allowed) and the 4f0-1s transition (forbidden). Other physical mechanisms of line broadening have not been retained and overlapping with neighboring lines has not been included for the sake of clarity. See text for details on the calculation.
Atoms 08 00074 g002
Figure 3. Plot of the π component of Lyman β, obtained from one calculation neglecting the quadratic Zeeman effect (dashed line) and one calculation retaining the quadratic Zeeman effect (solid line). In both cases, the Stark broadening due to ions only is retained and the overlap with neighboring lines has not been included for the sake of clarity. The spectrum is split in many components when the quadratic Zeeman is retained.
Figure 3. Plot of the π component of Lyman β, obtained from one calculation neglecting the quadratic Zeeman effect (dashed line) and one calculation retaining the quadratic Zeeman effect (solid line). In both cases, the Stark broadening due to ions only is retained and the overlap with neighboring lines has not been included for the sake of clarity. The spectrum is split in many components when the quadratic Zeeman is retained.
Atoms 08 00074 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rosato, J. Hydrogen Line Shapes in Plasmas with Large Magnetic Fields. Atoms 2020, 8, 74. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms8040074

AMA Style

Rosato J. Hydrogen Line Shapes in Plasmas with Large Magnetic Fields. Atoms. 2020; 8(4):74. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms8040074

Chicago/Turabian Style

Rosato, Joël. 2020. "Hydrogen Line Shapes in Plasmas with Large Magnetic Fields" Atoms 8, no. 4: 74. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms8040074

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop