Processing math: 19%
Advanced Search+
Qianglin HU, Wen HU. A classical relativistic hydrodynamical model for strong EM wave-spin plasma interaction[J]. Plasma Science and Technology, 2022, 24(3): 035001. DOI: 10.1088/2058-6272/ac3985
Citation: Qianglin HU, Wen HU. A classical relativistic hydrodynamical model for strong EM wave-spin plasma interaction[J]. Plasma Science and Technology, 2022, 24(3): 035001. DOI: 10.1088/2058-6272/ac3985

A classical relativistic hydrodynamical model for strong EM wave-spin plasma interaction

More Information
  • Author Bio:

    Qianglin HU: E-mail: 9920070002@jgsu.edu.cn

  • Received Date: July 21, 2021
  • Revised Date: November 13, 2021
  • Accepted Date: November 13, 2021
  • Available Online: January 29, 2024
  • Published Date: March 08, 2022
  • Based on the covariant Lagrangian function and Euler–Lagrange equation, a set of classical fluid equations for strong EM wave-spin plasma interaction is derived. Analysis shows that the relativistic effects may affect the interaction processes by three factors: the relativistic factor, the time component of four-spin, and the velocity-field coupling. This set of equations can be used to discuss the collective spin effects of relativistic electrons in classical regime, such as astrophysics, high-energy laser-plasma systems and so on. As an example, the spin induced ponderomotive force in the interaction of strong EM wave and magnetized plasma is investigated. Results show that the time component of four-spin, which approaches to zero in nonrelativistic situations, can increase the spin-ponderomotive force obviously in relativistic situation.

  • In the last two decades, the spin effect in plasma has attracted a lot of attentions because of its potential applications [13]. Spin effect is an important property of quantum plasma, and researches have shown that it may exist even on a classic scale, which is longer than the thermal de Broglie wavelength [4]. Researches have also shown that under certain conditions, the spin contributions can be more important than the usual quantum plasma corrections [5]. In the nonlinear situation, the spin force appears to be more importance and cannot be ignored even in moderate-density high-temperature plasma, due to the fact that the electrons with different spin eigenstates can be separated by the spin induced ponderomotive force [6].

    Some models have been developed to study the spin plasma. In nonrelativistic regime, a fluid theory is developed from the Pauli Hamiltonian and Madelung decomposition of electron wave function [5, 7]. Based on this fluid theory, some new phenomena, such as new propagation modes [8], magnetosonic solitons [9], instabilities [10], and so on, are revealed. In relativistic regime [1116], the EM wave-spin plasma interactions have been considered on three levels: classical physics, relativistic quantum mechanics, and quantum field theory. On classical physics level, a classical relativistic model is developed directly from the Bargmann–Michel–Telegdi equation. However, in this model the spin does not affect the electron dynamics [17]. On relativistic quantum mechanics level, using the Foldy–Wouthuysen ransformation of Hamiltonian for positive energy states of electron, a phase-space scalar kinetic model based on the gauge-invariant Stratonovich–Wigner function [18] and a fluid model based on the Madelung decomposition of electron wave function with two spinor [19] are developed for weakly relativistic spin plasma. On quantum field theory level, a set of fully covariant hydrodynamic equations are obtained based on Dirac theory of electrons [20]. However, the covariant model is very complex and is rather difficult to implement in practical situations, especially for numerical applications [19].

    In addition to these, there are some other works on spin plasma [2123]. In the present work, we developed a classical hydrodynamical model for relativistic spin plasma based on the covariant Lagrangian function and Euler–Lagrange equation. The equations of this model are formally consistent with the equations of the model proposed in [20]. However, this model is mathematically much simpler and is easy to implement in practical situations. Admittedly, as a classical model, some quantum plasma corrections, such as quantum statistic effects and quantum diffraction effects, are not included in this model. Therefore, this model can be applied to discuss the spin effects of strong relativistic plasma in classical regime, such as in astrophysics, high-energy laser-plasma systems and so on. Analysis shows that the relativistic effects may affect the physical processes by three factors: the relativistic factor, the time component of four-spin, and the velocity-field coupling.

    In the end, this model is used to examine the spin contribution to the ponderomotive force of strong EM wave interacting with magnetized plasma, which has been discussed by some authors [6, 24, 25] in nonrelativistic regime. Results show that the time component of four-spin, which approaches to zero in nonrelativistic situations, can increase the spin-ponderomotive force obviously in relativistic situation.

    The Lagrangian function of electron in electromagnetic (EM) field can be written as

    =0+F+S, (1)

    where 0=-mc2/˜γ is the Lagrangian function of free electron, m, c, ˜γ are the rest mass of electron, the light speed in vacuum and the relativistic factor of electron, respectively, F=-e˜γcuαAα is the ordinary electron-EM field interaction Lagrangian [26], e is the electric quantity of an electron, uα is the 4-velocity of electron and Aα is the 4-vector potential, S=σμ(˙σμ-αFμνσν) is the spin-EM field interaction Lagrangian [27], α=e/mc, and the Landé factor g=2 is presumed, σμ is the 4-vector spin, and Fμν is Maxwell electromagnetic field tensor. The ions are viewed as a uniformly positive back-ground.

    To obtain the motion equation from equation (1), we use Lagrange equations

    ddτ˙q-q=0, (2)

    where τ is the proper time, and q is the generalized coordinates. From equations (1) and (2), we obtain

    mduαdτ=ecFαβuβ+xα(ασμFμνσν), (3)

    and

    dσαdτ=αFανσν. (4)

    Equations (3) and (4) are the motion equations of single electron. Considering ddτ=˜γddt=˜γ[t+(u·)], equations (3) and (4) can be rewritten as

    muββuα=ecFαβuβ+xα(ασμFμνσν), (5)

    and

    uββσα=αFανσν, (6)

    where β(/ct,/xi) is the 4-gradient operator.

    Assuming that the electron distribution function is f(x,p,σ), the fluid density n in the rest frame of fluid can be written as

    n(x)=fd (7)

    The statistical average of any physical quantity Rμ is defined as

    Rμ = 1n Rμ fd4 pd4 σ. (8)

    Therefore, the fluid four-velocity is Uμ= uμ and the four-spin density is Sμ= σμ . The microscopic four dimensional velocity and spin density in the rest frame are also defined as vμ=uμ-Uμ and Σμ=σμ-Sμ, respectively. It is easy to obtain that vμ =0= Σμ .

    Taking the ensemble average of equations (5) and (6), we obtain

    UββUα= emc FαβUβ+ em2 c xα SμFμνSν- vββvα + em2 c xα ΣμFμνΣν -Uβ βvα , (9)

    and

    UββSα= emc FανSν- vββΣα -Uβ βΣα . (10)

    Moreover, if the creation of electron pairs can be neglected, the conservation of electron number gives

    μ nUμ =0. (11)

    Equations (9)–(11), plus the macroscopic Maxwell equations

    μFμν=enUν, (12)

    form a complete and covariant description of relativistic spin plasma.

    It is useful to put equations (9)–(11) into the general three-dimensional vector forms. To do this, we use the usual procedure [20] and rewrite the four-velocity of fluid as

    Uμ=γf[c,U] , (13)

    where U is the spatial component of the fluid velocity, γf is the average relativistic factor. Considering the constraint uμuμ=c2 and the definition of vμ, it is easy to obtain

    γf=γ(1- vμvμ ) 1/2 , (14)

    where γ= 1-U2 /c2 -1/2 . The macroscopic four-spin density is

    Sμ= S0 , S. (15)

    Considering the covariant constraint UμSμ=0, the time component can be obtained as S0 = 1/cU·S.

    Therefore, equation (9) can be written in the vectorial form

    d γfmc2 dt =eU·E+ eγfm t Ψs+ Pcγf , (16)

    and

    d γfmUdt =e E+ Uc ×B+ eγfmc Ψs - 1γfn ·Π+ Fcγf , (17)

    where

    Ψs=γfS· B- Uc ×E- 2γf2 γf+1 S0 Uc ·B, (18)
    Πij=mn vi vj , (19)
    Pc= em t ΣμFμνΣν -mc β vβ v0 - βvβ v0 -mcUβ βv0 , (20)
    Fci = emc xi ΣμFμνΣν -m0 v0 vi +m βvβ vi -mUβ βvi . (21)

    Equations (16) and (17) are the time and special components of equation (9), respectively. Ψs is the spin potential (the derivation is shown in appendix), and in the nonrelativistic limit, it is equivalent to the corresponding item in [5]. Π is the pressure tensor. Pc and Fc are the nonlinear corrections of the power and force caused by the microscopic four dimensional velocity and spin density in the rest frame of fluid.

    Equation (10) can be written as

    dS0 dt = eγfmc S·E+ PSγf , (22)

    and

    dSdt = eγfmc S0 E+S×B- 1γfmn ·K+ ΞSγf , (23)

    where

    Kij=mn vi Σj , (24)
    PS=- β vβ Σ0 - βvβ Σ0 -Uβ βΣ0 , (25)
    ΞSi =-0 v0 Σi + βvβ Σi -Uβ βΣi . (26)

    Equations (22) and (23) are the time and special components of equation (10), respectively. K is the coupling tensor of thermal-spin. PS and ΞS are the nonlinear corrections to the spin evolution caused by the microscopic four dimensional velocity and spin density in the rest frame of fluid.

    The continuity equation (11) becomes

    t nγf+· nγfU=0. (27)

    Equations (16), (17), (22), (23), and (27) form a full set of macroscopic hydrodynamical equations. Comparing with the nonrelativistic limit results, it is easy to know that the relativistic effects mainly affect the interaction processes by three factors: firstly, the relativistic factor γf; secondly, the time component of four-spin S0 , as shown in equations (18) and (23); and lastly, the velocity-field coupling, as shown in equation (18). This set of equations can be used to discuss the spin effects of relativistic plasma, which exist widely in astrophysics, high-energy laser-plasma systems and so on.

    As an example, we will illustrate how the time component of four-spin S0 affects the EM wave-spin plasma interaction process in the following part. In order to do this, we apply the formalism developed in the present work to discuss the spin induced ponderomotive force in magnetized plasma, which has been studied by some authors [6, 24, 25] in nonrelativistic regime. Take the assumption that the back-ground magnetic field is B0 =B0 z^ . The electrons have two eigenstates, i.e. spin-up and -down relative to B0 . Though the spin states of electrons can be disturbed by the EM waves, the differentiation of spin-up and -down electrons is still well defined provided the physics associated with spin flips can be neglected [4, 6]. Therefore, the electrons are treated as two kinds of populations in the present work, and equations (17) and (23) are reduced to

    m d γfα Uαdt =e E+ Uαc ×B + emc Sα· B- Uαc ×E- 2γfαγfα+1 Sα0 Uαc ·B , (28)
    dSαdt = eγfαmc Sα0 E+Sα×B, (29)

    in the cold plasma, where the subscript α= or , denotes spin-up or spin-down. In the nonrelativistic limit, i.e. Uαc, equations (28) and (29) reduced to the equations which have been used in [6].

    In the present work, the derivation of the spin-ponderomotive force is also based on the perturbative analysis [6]. In the following part, we assume a slowly varying plane circularly polarized EM wave E= 1/2 E~ exp i kz-ωt +c.c., where c.c. denotes complex conjugate, and define the variables E± =Ex±iEy, U± =Ux±iUy, S± =Sx±iSy. Separating the EM wave magnetic field and the back-ground magnetic field by letting B=B0 +B, and writing Sα into Sα=Sα0 +Sα, the 1ω -frequency component of equation (29) can be written as

    t -iωSα= eγfαmc Sα0 E+Sα0 ×B+Sα×B0 , (30)

    where S0 = ħ/2 z^ =-S0 . Considering that S0 = 1/cU·S, we have

    Sα0 = 1c U·Sα= 1c U± * Sα± . (31)

    The lowest order of Uα can be obtained from equation (28) as

    Uα± = ieγfαm ω±Ω E± . (32)

    Equation (32) shows that in the lowest order of approximation, the velocity of electron is not influenced by the spin, and from now on, the subscript α in Uα and γfα will be dropped. Thus, the lowest order of Sα can be obtained from equation (30) as

    Sα± = eSα0 γfmc ω±Ω 1- a2 / γf2 1±Ω/ω2 B± , (33)

    where Ω= eB0 / γfmc is the cyclotron frequency, and a= eE/ mωc is the normalized potential of EM wave. Using γf2 =1+ a2 / 1±Ω/ω2 , equation (33) reduced to

    Sα± = γfeSα0 mc ω±Ω B± . (34)

    Therefore, substituting equation (34) into the correction term of equation (30) gives

    Sα± = γfeSα0 mc ω±Ω B± ± i ω±Ω B± t . (35)

    The spin induced ponderomotive force can be written as [6]

    Ƒαz= emc Ψs= emc S·B = emc Sα± z B± * + S± * z Bα± . (36)

    Substitution of equations (34) and (35) into equation (36) gives

    Ƒαz= γfe2 Sα0 mc2 ω±Ω z - k ω±Ω t B2 . (37)

    It is easy to verify that in the nonrelativistic limit, equation (37) is the same with the result in [6]. Equation (37) shows that in the relativistic situation, the spin-ponderomotive force increased with the relativistic factor, i.e. increased with the EM wave intensity. Letting ξ=z-vgt, equation (37) can be transformed to the commoving frame

    Ƒαz= γfe2 Sα0 mc2 ω±Ω 1+ kvg ω±Ω ξ B2 . (38)

    The ratio of the spin-ponderomotive force in relativistic and nonrelativistic situation is

    R= ƑαzFαz γf. (39)

    where Fαz is the spin-ponderomotive force for nonrelativistic situation, i.e. Equation (13) in [6]. Generally, the ponderomotive force is decreased with the increase of the effective mass of the particle. Equation (39) implies that in relativistic situation, in spite of the electron effective mass is increased, the spin-ponderomotive force is γf times larger than the spin-ponderomotive force in nonrelativistic situation due to the influence of the time component of four-spin. Figure 1 shows the variation of R versus a for left-hand and right-hand circularly polarized (LCP, RCP) waves, where the EM wave length λ=1.06 μm and the back-ground magnetic field B0 =1000 T are presumed. Figure 1 obviously shows that R increased with the increase of a. Researches [6] have shown that in the EM wave-spin plasma interaction, the spin-ponderomotive force can separates the spin-up and spin-down electrons. The reason is that the spin-ponderomotive forces acting on the two species of electrons are in opposite direction [6, 24]. Equation (39) means that in the relativistic situation, the time component of four-spin can increase the degree of the separation of these two kinds of electrons.

    Figure  1.  The variation of R versus a. EM wave length λ=1.06 μm, B0 =1000 T.

    Spin effect is an important property of quantum plasma and several models have been proposed on it. In this work, we developed a classical hydrodynamical model for relativistic spin plasma based on the covariant Lagrangian function of electrons. This model can be used to discuss the spin effects of strong relativistic plasma. Analysis shows that the relativistic effects may affect the interaction processes through three aspects: firstly, the relativistic factor γf; secondly, the time component of four-spin S0 ; and lastly, the velocity-field coupling. However, it needs further research to know what phenomenon can be triggered by these factors.

    The present model is different from the previous models [5, 7, 1720]. (1) The present model can be used for strong relativistic cases. The existing models developed from quantum mechanics are mainly applied to non-relativistic [5, 7] or weakly relativistic [18, 19] cases. (2) The present model is mathematically simple and is easy to implement in practical situations. Though the model developed from quantum field theory [20] is a covariant model and includes as many quantum effects as possible, it is very complex and is rather difficult to implement in practical situations, especially for numerical applications [19]. (3) The present model naturally includes the main spin effects. The model proposed in [17] is also a relativistic model, but some crucial spin effects, such as the spin potential, thermal-spin coupling, and so on, are ignored. However, the equations of the present model are formally consistent with the model proposed in [20]. Admittedly, as a classical model, some quantum effects, such as quantum statistic effects and quantum diffraction effects, are not included in the present model. Therefore, it is applicable for the conditions that the plasma temperature is much higher than the Fermi temperature and the force induced by Bohm potential is much smaller than the spin force. Actually, these conditions are satisfied in numerous cases such as astrophysics, high-energy laser-plasma systems and so on.

    The model has been applied to examine the spin contribution to the ponderomotive force of strong EM wave interacting with magnetized plasma in this work. Results show that the time component of the four-spin can increase the spin-ponderomotive force obviously in relativistic situation.

    Considering that SμFμνSν fulfills the requirements of covariance, we have

    Ψs=SμFμνSν=S' μF' μνS' ν=B'·S', (A1)

    where the quantities with prime denote the quantities in the electron rest frame and S' 0 =0. In the last equality, the commutation relation S' i, S' j=ϵijkS' k is used. Using the general Lorentz transformation, one can obtain [26]

    B' =γf B- Uc ×E- γf2 γf+1 Uc Uc ·B, (A2)
    S' =S- γfγf+1 S0 Uc . (A3)

    Therefore,

    Ψs=γfS· B- Uc ×E- 2γf2 γf+1 S0 Uc ·B, (A4)

    where the term including U/c3 is neglected.

    This work was supported by National Natural Science Foundation of China (No. 12065011), and Science and Technology Research Project of Jiangxi Provincial Department of Education (No. GJJ170642).

  • [1]
    Wu Y T et al 2020 Phys. Rev. Appl. 13 044064 doi: 10.1103/PhysRevApplied.13.044064
    [2]
    An X Y et al 2019 Phys. Plasmas 26 123106 doi: 10.1063/1.5118710
    [3]
    Sahu B, Choudhury S and Sinha A 2015 Phys. Plasmas 22 022304 doi: 10.1063/1.4907658
    [4]
    Brodin G, Marklund M and Manfredi G 2008 Phys. Rev. Lett. 100 175001 doi: 10.1103/PhysRevLett.100.175001
    [5]
    Marklund M and Brodin G 2007 Phys. Rev. Lett. 98 025001 doi: 10.1103/PhysRevLett.98.025001
    [6]
    Brodin G, Misra A P and Marklund M 2010 Phys. Rev. Lett. 105 105004 doi: 10.1103/PhysRevLett.105.105004
    [7]
    Brodin G and Marklund M 2007 New J. Phys. 9 277 doi: 10.1088/1367-2630/9/8/277
    [8]
    Misra A P 2007 Phys. Plasmas 14 064501 doi: 10.1063/1.2737765
    [9]
    Marklund M, Eliasson B and Shukla P K 2007 Phys. Rev. E 76 067401 doi: 10.1103/PhysRevE.76.067401
    [10]
    Misra A P and Shukla P K 2008 Phys. Plasmas 15 052105 doi: 10.1063/1.2913265
    [11]
    Hafez M G et al 2020 AIP Adv. 10 065234 doi: 10.1063/5.0011086
    [12]
    Hafez M G 2020 Astrophys. Space Sci. 365 78 doi: 10.1007/s10509-020-03791-9
    [13]
    Hafez M G 2019 IEEE Trans. Plasma Sci. 47 5314 doi: 10.1109/TPS.2019.2949254
    [14]
    Alam M S et al 2018 Phys. Plasmas 25 072904 doi: 10.1063/1.5037788
    [15]
    Hafez M G et al 2017 Plasma Sci. Technol. 19 015002 doi: 10.1088/1009-0630/19/1/015002
    [16]
    Hafez M G, Talukder M R and Hossain Ali M 2016 Pramana 87 70 doi: 10.1007/s12043-016-1275-x
    [17]
    Melrose D B and Mushtaq A 2011 Phys. Rev. E 83 056404 doi: 10.1103/PhysRevE.83.056404
    [18]
    Asenjo F A et al 2012 New J. Phys. 14 073042 doi: 10.1088/1367-2630/14/7/073042
    [19]
    Morandi O et al 2014 Phys. Rev. E 90 013103 doi: 10.1103/PhysRevE.90.013103
    [20]
    Asenjo F A et al 2011 Phys. Plasmas 18 012107 doi: 10.1063/1.3533448
    [21]
    Hakim R and Heyvaerts J 1978 Phys. Rev. A 18 1250 doi: 10.1103/PhysRevA.18.1250
    [22]
    Sivak H D 1985 Ann. Phys. 159 351 doi: 10.1016/0003-4916(85)90117-4
    [23]
    Kowalenko V, Frankel N E and Hines K C 1985 Phys. Rep. 126 109 doi: 10.1016/0370-1573(85)90097-3
    [24]
    Hu Q L et al 2020 Phys. Lett. A 384 126261 doi: 10.1016/j.physleta.2020.126261
    [25]
    Stefan M et al 2011 Phys. Rev. E 83 036410 doi: 10.1103/PhysRevE.83.036410
    [26]
    Jackson J D 1999 Classical Electrodynamics3rd ednNew YorkWiley 555
    [27]
    Fontana A Q 2011 Lat. Am. J. Phys. Educ. 5 406
  • Related Articles

    [1]Runhui WU (邬润辉), Song CHAI (柴忪), Jiaqi LIU (刘佳琪), Shiyuan CONG (从拾源), Gang MENG (孟刚). Numerical simulation and analysis of lithium plasma during low-pressure DC arc discharge[J]. Plasma Science and Technology, 2019, 21(4): 44002-044002. DOI: 10.1088/2058-6272/aafbc7
    [2]Jun DENG (邓俊), Liming HE (何立明), Xingjian LIU (刘兴建), Yi CHEN (陈一). Numerical simulation of plasma-assisted combustion of methane-air mixtures in combustion chamber[J]. Plasma Science and Technology, 2018, 20(12): 125502. DOI: 10.1088/2058-6272/aacdef
    [3]Guobao FENG (封国宝), Wanzhao CUI (崔万照), Lu LIU (刘璐). Dynamic characteristics of charging effects on the dielectric constant due to E-beam irradiation: a numerical simulation[J]. Plasma Science and Technology, 2018, 20(3): 35001-035001. DOI: 10.1088/2058-6272/aa9d0d
    [4]Gui LI (李桂), Muyang QIAN (钱沐杨), Sanqiu LIU (刘三秋), Huaying CHEN (陈华英), Chunsheng REN (任春生), Dezhen WANG (王德真). A numerical simulation study on active species production in dense methane-air plasma discharge[J]. Plasma Science and Technology, 2018, 20(1): 14004-014004. DOI: 10.1088/2058-6272/aa8f3c
    [5]R. KHOSHKHOO, A. JAHANGIRIAN. Numerical Simulation of Stall Flow Control Using a DBD Plasma Actuator in Pulse Mode[J]. Plasma Science and Technology, 2016, 18(9): 933-942. DOI: 10.1088/1009-0630/18/9/10
    [6]DUANMU Gang(端木刚), ZHAO Changming(赵长明), LIANG Chao(梁超), XU Yuemin(徐跃民). Numerical Simulation of Dual-Channel Communication of Column Plasma Antenna Excited by a Surface Wave[J]. Plasma Science and Technology, 2014, 16(11): 1059-1062. DOI: 10.1088/1009-0630/16/11/11
    [7]DUAN Yaoyong (段耀勇), GUO Yonghui (郭永辉), QIU Aici (邱爱慈). Shock Wave and Particle Velocities of Typical Metals on Shock Adiabats[J]. Plasma Science and Technology, 2013, 15(8): 727-731. DOI: 10.1088/1009-0630/15/8/02
    [8]ZHUANG Juan (庄娟), SUN Jizhong (孙继忠), SANG Chaofeng (桑超峰), WANG Dezhen (王德真). Numerical Simulation of VHF E®ects on Densities of Important Species for Silicon Film Deposition at Atmospheric Pressure[J]. Plasma Science and Technology, 2012, 14(12): 1106-1109. DOI: 10.1088/1009-0630/14/12/13
    [9]DENG Yongfeng(邓永锋), TAN Chang(谭畅), HAN Xianwei(韩先伟), TAN Yonghua(谭永华). Numerical Simulation of the Self-Heating Effect Induced by Electron Beam Plasma in Atmosphere[J]. Plasma Science and Technology, 2012, 14(2): 89-93. DOI: 10.1088/1009-0630/14/2/01
    [10]WU Junhui, WANG Xiaohua, MA Zhiying, RONG Mingzhe, YAN Jing. Numerical Simulation of Gas Flow during Arcing Process for 252kV Puffer Circuit Breakers[J]. Plasma Science and Technology, 2011, 13(6): 730-734.

Catalog

    Figures(1)

    Article views (125) PDF downloads (97) Cited by()

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return