Doppler Broadening of Spectral Line Shapes in Relativistic Plasmas

Transcription

1 atoms Article Doppler Broadening of Spectral Line Shapes in Relativistic Plasmas Mohammed Tayeb Meftah 1,2, *, Hadda Gossa 1, Kamel Ahmed Touati 1,3, Keltoum Chenini 1,4 and Amel Naam 1,2 1 Laboratoire de Recherche de Physique des Plasmas et Surfaces (LRPPS), UKMO Ouargla 30000, Algerie; hadda.gossa@gmail.com (H.G.); ktouati@yahoo.com (K.A.T.); k1_chenini@yahoo.fr (K.C.); naamnaam10@gmail.com (A.N.) 2 Département de Physique, Faculté de Mathématiques et Sciences de la matière, Université Kasdi-Merbah, Ouargla 30000, Algerie 3 Lycée professionnel les Alpilles, Rue des Lauriers, Miramas, France 4 Département des Sciences et Technologies, Faculté des Sciences et Technologies, Université de Ghardaia, Ghardaia 47000, Algerie * Correspondence: mewalid@yahoo.com or meftah.tayeb@univ-ouargla.dz; Tel.: Received: 5 January 2018; Accepted: 27 March 2018; Published: 4 April 2018 Abstract: In this work, we report some relativistic effects on the spectral line broadening. In particular, we give a new Doppler broadening in extra hot plasmas that takes into account the possible high velocity of the emitters. This suggests the use of an appropriate distribution of the velocities for the emitters. Indeed, the Juttner-Maxwell distribution of the velocities is more adequate for relativistic velocities of the emitters when the latter are in plasma with an extra high temperature. We find an asymmetry in the Doppler line shapes unlike the case of the traditional Doppler effect. Keywords: plasmas; Maxwell; Juttner-Maxwell; relativistic; Doppler effect; asymmetry 1. Introduction The Doppler effect, discovered by physicist and mathematician Christian Doppler in the nineteenth century, is the modification of the frequency of a wave when the emitting source and the receiver are in relative motion. The frequency change also implies that of the period and the wavelength. This effect concerns both mechanical waves and electromagnetic waves. In plasmas, the neutral atoms, molecules or ions moving inside the plasma are similar to the moving antennae. Atoms or ions subjected to the Doppler effect, exhibit the well-known phenomenon: the Doppler broadening of the line profile. The investigation in recent decades of the derivation and illustration of the Doppler effect, especially the generalized relativistic Doppler effect, is still being actively pursued today [1 5]. In the following we will illustrate this method to formulate the classical Doppler effect first, and then the same method is formulated to get the relativistic Doppler effect on the Doppler broadening of the line profile. In our work, we will present a better derivation allowing quick and exact expressions of the classic and relativistic Doppler effect on the broadening of the spectral lines observed in the plasmas. In the formulation of the classical Doppler effect, we used the Maxwell velocity distribution for the emitters, while in the formulation of the relativistic Doppler effect we used the relativistic Juttner-Maxwell distribution. The latter is justified for the case of very high temperatures (in the range K) such as that encountered in fusion plasmas, in astrophysics, in cosmology (primordial Universe) and in unstable Z Pinch experiments [5]. Indeed, the Juttner-Maxwell distribution remains valid for all temperatures since it is more general than the Maxwell distribution. Atoms 2018, 6, 16; doi: /atoms

2 Atoms 2018, 6, 16 2 of 7 2. Doppler Broadening 2.1. Classical Doppler Broadening: Non Relativistic Case Often the emission (or absorption) of radiation by a particle (atom, ion, etc.) occurs during the movement. By the Doppler effect, the observed frequency in the observer (at rest) frame (see Figure 1) is different from the frequency emitted in the atom frame. The mean particle velocity at thermodynamic equilibrium is related to the temperature of the medium. Hence the broadening of the statistical Doppler effect is related to the distribution of the velocities of the emitter at the temperature T of the medium and the mass m of the emitter. Figure 1. The fixed frame where the emitter moves with a velocity V forming an angle θ with the observation direction Ox. One can assume a motionless observer, looking at an emitting atom moving with a velocity V in a direction forming an angle θ with the direction of observation (Ox) (see Figure 1), records a shifted angular frequency ω with respect to the angular eigenfrequency ω 0 of the emitter assumed to be stationary. This angular frequency is given by ω(v x ) = ω 0 (1 V c cos θ) = ω 0(1 V x c ) (1) where c is the velocity of the light in vacuum. The normalized intensity (normalized to one) of the line at the angular frequency ω is given by the average over the normalized Maxwell distribution (normalized to one). of the Dirac delta distribution as the following f Maxwell (V x ) = (m/(2πk B T)) 1/2 exp( mv 2 x /(2k B T)) (2) m I(ω) =< δ(ω ω(v x )) > Maxwell = 2πk B T exp( m 2k B T V2 x )δ(ω ω(v x ))dv x (3) Using the integral representation of the Dirac delta distribution (u is the integration variable whose unit is the second) δ(ω ω(v x )) = 1 + exp(iu(ω ω(v x )))du (4) 2π

3 Atoms 2018, 6, 16 3 of 7 We find the normalized intensity (normalized to one) as I(ω) = 1 m 2π 2πk B T exp(iuω)du = 1 2π exp( m 2k B T V2 x iuω(v x ))dv x (5) 2πmc 2 k B Tω 2 0 ) exp ( mc2 ( ω 1)2 2k B T where ω = ω ω 0, k B is the Boltzmann constant and m is the emitter mass. This is the formula of the intensity of the line in the non relativistic case. It is symmetric (Gaussian) around the central angular frequency ω 0. We note that the integrals in Formula (6) are strongly convergent because we deal with purely Gaussian integrals. The full width at the half maximum (FWHM) is given by the well known formula (in angular frequency unit) ω Doppler = ω 0 ( 8k BT ln(2) T mc 2 ) = ω 0 M where M is the mass of the emitter in atomic mass unit whereas T is the temperature in Kelvin Relativistic Doppler Broadening When an observer at a rest, recording the emitted radiation from a moving atom (or ion) with relativistic velocity V, they find that the angular frequency of this radiation is equal to [6]: ω(β) = ω 0 γ(1 + β cos θ) (8) (6) (7) where β = V/c, (9) γ = 1/ (1 β 2 ) (10) and ω 0 is the angular eigenfrequency and θ is the angle between the velocity of the emitter and the observation direction (Ox) (see Figure 1). By using the normalized Juttner-Maxwell distribution (normalized to one) [7] W J M (β)dβ = λ γ5 β 2 dβ exp ( λγ) (11) K 2 (λ) where λ = mc 2 /(k B T) (12) and K 2 (X) is the modified Bessel function of order two, we obtained the normalized relativistic intensity (normalized to one) of the line profile I(ω) =< δ(ω ω(β)) > Juttner Maxwell = W J M (β)dβ δ(ω ω(β)) (13) = 1 4π du W J M (β)dβ exp(iu(ω ω(β))) sin θdθdφ (14) Here, we have introduced the integral over the spherical angles that makes the emitter velocity with the fixed frame axis (see Figure 1). We have replaced the Dirac delta distribution by its integral representation by integrating over the variable u. Replacing the Juttner-Maxwell distribution W J M (β)

4 Atoms 2018, 6, 16 4 of 7 given by (11) and ω(β) given by (8) in Formula (14), we reach a more suitable expression of the relativistic intensity of the line profile I(ω) = + λ 4π K 2 (λ) or after integration on θ between zero and π; exp(iuω)du 1 0 γ5 β 2 exp( λγ) exp( iuω 0 γ)dβ π 0 exp(iuω 0γβ cos θ) sin θdθ (15) I(ω) = λ 1 2 K 2 (λ) λ 1 2 K 2 (λ) 0 γ4 β exp( λγ)dβ 0 γ4 β exp( λγ)dβ + + exp(iuω iuω 0 γ+iuω 0 γβ) iuω 0 du exp(iuω iuω 0 γ iuω 0 γβ) iuω 0 du (16) Finally, the integration over u, allows us to get the relativistic intensity of the line profile I( ω) = λ ( ) γdγ exp( λγ) S( ω γ + γ 2 K 2 (λ) 2 1) S( ω γ γ 2 1) 1 (17) where S(t) = +1 if t > 0 and S(t) = 1 if t < 0 and ω = ω ω 0 is the reduced angular frequency. We can manage the formula to be more suitable for the numerical treatment: I( ω) = λ 2 K 2 (λ) exp( λ) 1 (S( γdγ exp( λ(γ 1)) ω γ + γ 2 1) S( ω γ ) γ 2 1) (18) if we put λ(γ 1) = y, (λ is given by formula (12)) we find I( ω) = exp( λ) 2 K 2 (λ) 0 ( y λ + 1)dy exp( y) ) (S( ω y λ ( 1 + y λ )2 + 2 y λ ) S( ω y λ ( 1 y λ )2 + 2 y λ ) (19) We note that the integral in the last formula is convergent because we deal with the integral in distribution sense [8]. Unlike the classical Doppler effect, the relativistic one has a property: an asymmetric broadening as it is shown clearly in the following table corresponding to the temperatures in the range K. We remark also that in the relativistic case, as in the classical case, the central frequency is unchanged (see Figure 2). The maximum value of the intensity is at ω = 1 both for classical and relativistic case but the maximum value of the relativistic case is smaller than the maximum of the classical case (the asymmetry at ω = 1) is negative, see Figure 3. We mention that in the Table 1, we have denoted by ˆω L,R the value of the reduced angular frequency at the left and the right of the peak of the line (centred at ˆω = 1). ˆω L and ˆω R are chosen to be symmetrical with respect the centre of the line at ˆω = 1 and giving intensities very close to the half of the maximum of the intensity. Strictly speaking, we have considered ˆω R = 1 + g and ˆω R = 1 g with g = for Fe +25, for W +73, for Fm +99 (produced in nuclear reactions) and for Cn +111 (it is synthesized in laboratories for use in nuclear reactions). As we see in this table, the value of the intensity at the right I( ˆω R ) is greater than the intensity at the left I( ˆω L ). This remark shows clearly that the line profile has an asymmetry as defined by [9] (see the definition at the last line in the above table). Another feature in this study is that we have not specified the line profile, because we have used the reduced angular frequency ω = ω ω 0 : for each specific line (specific transition), we must multiply ˆω (the x-axis) by the corresponding angular eigenfrequency ω 0 to obtain the intensity I(ω). If we define the asymmetry as [10] Asym = I(ω, relativistic) I(ω, classical) (20)

5 Atoms 2018, 6, 16 5 of 7 which is the difference between the normalized relativistic intensity given by (18) and the normalized classical intensity given by (6), we obtain the following figure representing the asymmetry for Cn +111 at K. Figure 2. Relativistic intensity as defined by Formula (19) for Iron at T = K. Figure 3. Asymmetry as defined by Formula (20) for Cn +111 at T = K. Table 1. Asymmetry percentages for different hydrogen-like ions. T = 10 8 K, Fe +25 T = K, W +73 T = K, Fm +99 T = K, Cn +111 ˆω L ˆω R I( ˆω L ) I( ˆω R ) I( ˆω R) I( ˆω L) 100 I( ˆω R)+I( ˆω L)

6 Atoms 2018, 6, 16 6 of 7 In Figure 3, we see clearly that, in the left of ˆω = ω/ω 0 = , the intensity of the relativistic profile is lower than of the classical profile, whereas it is higher in the right of ˆω = ω/ω 0 = It can be seen clearly in this figure that the asymmetry is a function of ( ˆω = ω/ω 0 ) and that means that for any line, the asymmetry is as indicated in this figure. To obtain the asymmetry, for a specific line centred at w 0, we must multiply ˆω by w 0. The same remark holds for Figure 4 for the hydrogen-like Iron, but with a more pronounced asymmetry since the maximum of the asymmetry is equal to 2.40 for the iron (Fe +25 ) whereas for the Copernicium (Cn +111 ) is equal to Figure 4. Asymmetry as defined by Formula (20) for Fe +25 at T = K. 3. Conclusions In this work, we report some relativistic effects on the spectral line broadening. In particular, we obtained a new expression for the Doppler broadening that takes into account the possible high velocity of the emitters. This suggests the use of an appropriate distribution of the velocities for emitters. We find, an asymmetry in the Doppler broadening unlike the well known classical Gaussian Doppler broadening. Acknowledgments: We wish to acknowledge the support of LRPPS laboratory and its director Pr: Fethi Khelfaoui, by offering us the encouragement, and some technical materials for developing this work. Author Contributions: All authors M.T. Meftah, H. Gossa, K.A. Touati, K. Chenini and A. Naam were participated equivalently to this theoretical work. Conflicts of Interest: We declare no conflict of interest to signal. References 1. Huang, Y.-S. Formulation of relativistic Doppler-broadened absorption line profile. Europhys. Lett. 2012, 97, Huang, Y.-S. Formulation of the classical and the relativistic Doppler effect by a systematic method. Can. J. Phys. 2004, 82, Kichenassamy, S.; Krikorian, R.; Nikogosian, A. The relativistic Doppler broadening of the line absorption profile. J. Quant. Spectrosc. Radiat. Transf. 1982, 27, McKinley, J.M. Relativistic transformations of light power. Am. J. Phys. 1979, 47,

7 Atoms 2018, 6, 16 7 of 7 5. Haines, M.G.; LePell, P.D.; Coverdale, C.A.; Jones, B.; Deeney, C.; Apruzese, J.P. Ion Viscous Heating in a Magnetohydrodynamically Unstable Z Pinch at Over Kelvin. Phys. Rev. Lett. 2006, 96, Jackson, J.D. Special Theory of Relativity. In Classical Electrodynamics, 3rd ed.; John Wiley: New York, NY, USA, 1962; Chapter 11, pp Zenitani, S. Loading relativistic Maxwell distributions in particle simulations. Phys. Plamas 2015, 22, Stehlé, C.; Gilles, D.; Demura, A.V. Asymmetry of Stark profiles: The microfield point of view. Eur. Phys. J. D 2000, 12, Schwartz, L. Théorie des Distributions; Editions Hernmann: Paris, France, Huang, Y.-S.; Chiue, J.-H.; Huang, Y.-C.; Hsiung, T.-C. Relativistic formulation for the Doppler-broadened line profile. Phys. Rev. A 2010, 82, (R). c 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (