Effect of kappa distribution on the damping rate of the obliquely propagating magnetosonic mode
1.Introduction
Magnetosonic(MS)waves are low frequency and highly dispersive waves[1]which propagate perpendicularly to the ambient magnetic field[2].These waves are also sometimes called fast MS waves or equatorial noise when their frequency is above proton gyro frequency all the way up to the lower hybrid frequency[3].MS waves are of great importance because of their ubiquity in space plasmas.In the case of obliqueness,i.e.the presence of small parallel wavenumber,the waves exhibit wave–particle interaction(WPI),in which case the wave may either grow or become damped out depending on whether a free energy source is present or not.The WPI leads to mechanisms such as heating and particle acceleration in the solar wind[4],magneto sheath[5],plasma sheet[6],corona[7],magnetic equator[8],aurora[9],solar flares[10,11],inner magnetosphere[12]and Van Allen radiation belts[13].Based on the measurements of Van Allen Probe A,fast MS waves are believed to penetrate into the plasma sphere via the plasma pause and have been proven to be the best candidates for controlling the dynamics of electrons in the inner radiation belt and the unusual butter fly distribution of relativistic electrons in the geostationary orbit[3,13,14].
Depending on the nature of the plasma,WPI is affected by various effects such as the presence of dust[15–18],temperature anisotropy[19]and density inhomogeneity[20].In Maxwellian distributed plasmas,studies have shown that increasing the density inhomogeneity or temperature anisotropy withenhances the growth rate[20],whereas the presence of dust reduces it[15].
加大盾尾油脂注入量,持续观察盾尾油脂腔压力,确保盾尾密封良好。在左线盾构完成穿越施工并满足隧道安全距离要求后,进行右线盾构掘进穿越施工。
Numerous research studies have been conducted to study MS waves in non-Maxwellian plasma.The WPI of MS waves under the assumption ω ≪ ωi has been investigated in a loss cone distribution[21]as well as in an anti-loss cone one[22].For the loss-cone distribution,depending upon the choice of the loss-cone parameters,the wave growth may either be increased or decreased[21].However,in the case of the antiloss cone,the wave is always damped and is smaller as compared to the Maxwellian distribution[22].
Apart from the aforesaid distributions,the non-Maxwellian character,especially the presence of suprathermal particles,can be modeled effectively by kappa-type distribution functions.Several kappa-type distributions appear in the literature:non-relativistic kappa[23,24],relativistic kappa[25–28],kappa-Maxwellian[29]and many more.In the nonrelativistic limit,the values of kappa from 2 to 6 provide a good fit for many regions in the space plasma[23].However,at very high energy,the non-relativistic kappa distribution deviates from the observed particle distribution;a relativistic kappa distribution seems to fit the data very well[25–28].The kappa-Maxwellian distribution is advantageous both physically and mathematically.Usually,in space magnetoplasma,particles are accelerated along the magnetic field direction,leading to a power law distribution;however,in the perpendicular direction,there is equilibrium and is otropization [29,30].From the mathematical standpoint,the perpendicular integration in the generalized dielectric tensor can be handled with ease.Depending on the environment,each distribution has its own importance,and all the environments discussed in the first paragraph where the MS waves play an important role in WPI have been con firmed to fit with the kappa-type distribution.
In this paper,we discuss the damping rate of the obliquely propagating MS waves in a suprathermal plasma modeled by kappa-Maxwellian distribution.To the best of our knowledge,ours is the first report that addresses this issue.The layout of this paper is as follows.In section 2,we present the mathematical model and derive the dispersion relation along with the analytical expressions for the real frequency and the damping rate.Then,we discuss the exact numerical analysis in section 3.Finally,we summarize the results and conclude the paper in section 4.
To manipulate χ ,we choose the following kappa-Maxwellian distribution function[29]
and
2.Mathematical model
with and
![](https://qkk.lw71.com/uploads/image/2018/3/a07419b422c10ad0c33f805e9e51a745.png)
where we have assumed integer values of κ.With this choice ofZκ*(ζ0α),we extract the real part of the dielectric tensor from equation(1)as
![](https://qkk.lw71.com/uploads/image/2018/3/5ecae983f54499005d47222236719dcf.png)
We assume low frequency waves, ω ≪ Ωi ,in a uniform low beta plasma,in which the real wave vector k lies in the x-z plane.The dispersion relation in such a plasma can be written as[15,19,20,31,32]
![](https://qkk.lw71.com/uploads/image/2018/3/8d8eb8dfe099ed0e9842e5bfe5139849.png)
雾化吸入治疗过程中,确保患儿采取坐位或者半坐卧位,使膈肌下移,便于雾化吸入后肺部充分扩展,增加气体的交换量。
会计信息体统以计算机为基础,对当前的数据分析更加的快捷、便利。传统的管理模式无法适应当前信息系统的发展,阻碍了其加快的步伐。而云计算的优势便由此日益凸显。在提高管理信息时效性的同时,节省了大量的空间,缩小了内存,提高了信息的安全性能。不可避免的是,网络时代的发展使得众多资源共享。由此,不良商家为了顾及个人的利益,通过一些非法的手段以谋取个人利益。后果不堪设想。综上而言,要积极的定期对此类事件进行排查,可以在一定晨读上保障信息的安全状况,同时,完善健全监督机制,实现信息安全化的发展。
![](https://qkk.lw71.com/uploads/image/2018/3/3818a086b173d1bf753e8baef6b6be76.png)
whereThe perpendicular and parallel thermal velocities are related to their corresponding temperatures by
and
with
The spectral index κ controls the number of high energy particles in a system.The smaller the value of κ,the larger the number of high energy particles.Moreover,in the limit κ → ∞,the Maxwellian distribution can be retrieved.
管理会计是从传统会计中分离出来与财务会计并列的、着重为企业改善经营管理、提高经济效益服务的一个会计分支。管理会计主要服务于企业内部经营管理,通过对财务会计资料、统计资料及其他业务资料在内的各种数据进行加工,向企业内部管理者提供有用的经营信息以帮助他们正确确定经营目标,进行经营决策,对企业经营活动实施控制,使企业的经营活动处于最优状态。
Following the same procedure adopted by the authors in references[15,19,20,31,32],we get a simplified expression for∈as
![](https://qkk.lw71.com/uploads/image/2018/3/93b9890779d3f7602f2c1c65c63451f4.png)
whereis the gyro radius and ζ0α=
Generally, we assume ω to be complex,
satisfying the condition γ ≪ ωr .Consequently,the modified plasma dispersion function is given by[23,24]
![](https://qkk.lw71.com/uploads/image/2018/3/ca85900938aa2a11c21aa35e38580c57.png)
with the supposition that .As we see,
is a complex-valued function which,in turn,makes the dielectric tensor,equation(1),a complex-valued function that can be separated into real and imaginary parts: D( k,ω)=
Invoking Taylor’s series with the assumption
it can easily be shown that
evaluated at ω=ωr .
For MS waves,as we know,the parallel wavenumber is small compared to the perpendicular wavenumber which means ζ0α ≫ 1should be used in equation(4),in which case,the asymptotic form of the plasma dispersion function needs to be taken into account,i.e.,[23,24]
![](https://qkk.lw71.com/uploads/image/2018/3/0b88d68e42a507965df778ea00043259.png)
where
![](https://qkk.lw71.com/uploads/image/2018/3/e9f7ba0d269f014f0b8cfd407ef7d4f8.png)
where the sound velocity and Alfvén velocity are,respectively,given by and
IfD r= 0,then it is straightforward to obtain
![](https://qkk.lw71.com/uploads/image/2018/3/dde3eadee3d1ad1f0ec036244ac68cd7.png)
which means that the real frequency is independent of κ,in much the same way as the real frequency of electrostatic waves.This has to do with the structural similarities of the two dielectric tensors(equations(40)and(3)of references[24]and[33],respectively).
Substituting the above expression in we obtain
![](https://qkk.lw71.com/uploads/image/2018/3/1533d401b046637da3e97e2eab455316.png)
For the WPI,as is known,we take into account the imaginary part of the plasma dispersion function with which the resulting imaginary part of the dielectric tensor becomes
![](https://qkk.lw71.com/uploads/image/2018/3/1a48e3a781ac326ac6359851a2a63e48.png)
which clearly shows that the wave must always be damped;a negative sign cannot be generated from any values defined above.Also,it should be noted that the termgets pretty small for large
,ensuring a reduced damping rate.
Since an obliquely propagating mode and its WPI in the solar wind have been observed[34],hence we choose the following parameters appropriate for the solar wind at about 1 AU[35]:10-5G and β= 0.0057.We numerically solve equation(1),
3.Numerical results and discussion
keeping the plasma dispersion function unexpanded.For the exact numerical analysis,we consider low integer values of the spectral index κ,i.e., κ =2,3 and 4,which were experimentally observed in the aforementioned environment[36].The other reason for taking low integer values is that they give closed-form expressions for the plasma dispersion function[23]
(三)改善教学方法。对于高中生传统文化素养的培育,大多依赖于语文课课堂上教师对古诗词及文言文的讲授方式进行,受应试教育的影响,教师习惯在这一过程中落入翻译及要求背诵的教学模式中。
![](https://qkk.lw71.com/uploads/image/2018/3/12ad2f22386e531b5f4d6699f5fe3ba9.png)
血痂儿结在嘴上 Upon lips, blood scab is ruthlessly growing
![](https://qkk.lw71.com/uploads/image/2018/3/3331070a98836d073aee19f859778598.png)
![](https://qkk.lw71.com/uploads/image/2018/3/965a873a7ae2fd6619ea5693ccfb9014.png)
Figure 1.Panel(a)shows the exact numerical plots of equation(1)for the normalized real frequencyagainst normalized wavenumber
for different values of the spectral index κ.In panel(b),the frequency is in the large wavenumber regimes.The Maxwellian curve(solid line)is obtained using equation(6)to make Comparison.
These expressions make the numerical calculations much easier as compared to the cases of non-integer values of κ.We also note that there are no closed-form expressions for both the Maxwellian distributions and the non-integer values of κ.
As shown figure 1(a),the real frequency is almost independent of the parameter κ if we take to be small.Even for moderate values of
(see figure 1(b)),the real frequencies for different κ are indistinguishable.The reason is that for the MS mode,
which implies thatζ0α must be larger than 1.When we carefully observe the behavior of Zκ*(ζ0α)for large ζ0α (see equation(5)),we see that the expression in the parentheses in the real dispersion relation(see equation(4))for ζ0α ≫ 1is not a sensitive function of κ unless higher-order terms(greater than five)are invoked.
![](https://qkk.lw71.com/uploads/image/2018/3/60b7c7949e8b66719aa935f6c6f75152.png)
Figure 2.The normalized damping ratebased on equation(8)vs normalized wavenumber
for different values of the spectral index κ.
As is evident from equation(8),κ does not provide the free energy source;the wave is always damped,and the damping rate strongly depends on the spectral index κ.Figure 2 is the plot of equation(8),while figures 3(a)and 3(b)are based on the exact numerical analysis.As expected,for small values of the wave is heavily damped,when κ is low,although the difference is minuscule in magnitude.When
increases,the wave behaves otherwise,i.e.,it now becomes heavily damped for large values of κ,and the difference in the magnitude is relatively substantial.The underlying mechanism behind this is the participation of resonant particles.From the given data,we have
.As the values of
become small,ζ0α gets large and is shifted towards the tail of the distribution(see figure 4(c))where the number of resonant particles is large for a low spectral index.Moreover,the curves are quite close to each other,and thus the difference between the two damping rates is not appreciable.Table 1 depicts the numerical values of the damping rate for different values of the spectral index κ at the given small
and large ζ0e.Obviously,the values of
from 0.03 to 0.06 demand ζ0eto be in the range of 7.23–3.62.
课堂教学质量预警指标体系也应该分层、分类反映学校、专业、课程三个层面的基本情况。在设计课堂教学质量预警指标体系时,要充分考虑后期的数据分析,上一级的预警需求和目标决定了下一级预警要素的细分和关联,下一级的数据采集又影响着上一级数据状态的呈现。所以只有分层递阶考虑课堂教学预警指标设计,预警才能做到精确、及时、全面。
Whenincreases, ζ0e ,still being larger than 1,is shifted to the left where the resonant particles are more in number for high values of κ.Figure 4(b)shows the marked difference between the curves which explains Significant damping.From table 2,it is clear that for
or above,the values ofζ0α ranges from 1.44 to 1.14 wherein the large number of resonant particles for higher values of κ are seen from figure 4(b).
利用以上设置的相关参数,得到的内河水域中的流场分布见图2。图2中同时标示漂移物体在不同时刻可能位置的区间范围,当矩形区域与航道边界存在交点时,就存在停止继续漂移的可能性,假设物体的可能位置在相应时刻在矩形区域内为均匀分布,可得到物体最终漂移位置的概率分布情况。
4.Summary and conclusion
![](https://qkk.lw71.com/uploads/image/2018/3/f72f3d59912dd34c83ceb560a6fb79c7.png)
Figure 3.The exact numerical plot of equation(1)for the normalized damping ratevs normalized wavenumber
for different κ values.In(a),the damping rate is in the small wavenumber regimes which is too small to be seen in(b)which portrays both small as well as large wavenumber regimes.
Suprathermal particles have been regularly observed in laboratory as well as in space plasmas.That ambience affects the dynamics of the system has astonishing consequences.In the present work,we have considered a kappa-Maxwellian distribution with κ being an integer.We have shown that the real dispersion relation is not affected considerably by κ variation,however,the damping rate is.We have used the data from the solar wind at about 1AU and have found that for small values of the parallel wavenumber,the high energy resonant particles participate in the WPI,leading to a high damping rate.When parallel wavenumber increases,the resonant particles are mainly in the low energy part of the distribution function.But in the low energy part,there are more resonant particles for higher values of κ,contrary to what is observed in the high energy tail.
The results we have obtained may find fruitful applications in various suprathermal plasma environments,such as the solar corona,magnetosphere,aurora,plasma sheets and magnetosheath.Some of these environments may also contain suprathermal particles with non-integer values.However,following Mace[24]and Summer[37],we can extend our analysis by reconsidering the modified plasma dispersion function for non-integers.Furthermore,in addition to the presence of high energy particles in the aforesaid space environments,they may contain streaming particles or have temperature anisotropy or are inhomogeneous in density.By incorporating these effects,our calculations can be further extended,which will be the focus of our attention in our future research work.
![](https://qkk.lw71.com/uploads/image/2018/3/645be3088e23a4e3fa4db0de78662ba3.png)
Figure 4.In panel(a),parallel part of the normalized distribution functionagainst normalized velocity
for different κ values is shown.The middle panel(b)shows that when κ decreases,the number of suprathermal particles decreases when
is between 1 and 2.If
is greater than 3,then,as shown in panel(c),suprathermal particles are large in number for smallκ.
Table 1.Values of normalized damping rate for different
,ζ0e and κ.The normalized damping rates are expressed in
![](https://qkk.lw71.com/uploads/image/2018/3/8f1330dd19ad588eb0f7d954eb74d989.png)
k ∣∣k⊥ ζ0e k vAκ=γ ⊥, 2k vA κ=γ ⊥, 3k vA κ=γ ⊥, 4 0.03 7.23 −0.36 −0.094 −0.022 0.04 5.43 −1.46 −0.644 −0.246 0.05 4.43 −4.33 −2.758 −1.485 0.06 3.62 −10.00 −8.702 −5.982
Table 2.Values of normalized damping rate for different
ζ0e and κ.The normalized damping rates are expressed in10-3.
![](https://qkk.lw71.com/uploads/image/2018/3/3eec6ff6b5f929f2d1d8d0e38996522a.png)
k⊥ ζ0e k vA⊥k vA ⊥kvA⊥0.15 1.44 −0.57 −0.94 −1.11 0.16 1.33 −0.73 −1.19 −1.38 0.18 1.20 −1.11 −1.73 −1.98 0.19 1.14 −1.33 −2.02 −2.28
IAK is grateful to the Salam Chair in Physics for providing a vital research environment.
References
[1]Cooper F C,Nakariakov V M and Williams D R 2003 Astron.&Astrophys 409 325
[2]Cravens T E 1997 Physics of Solar System Plasmas(Cambridge,UK:Cambridge University Press)
[3]Xiao F L et al 2012 J.Geophys.Res.117 A06208
[4]Marsch E,Goertz C K and Richter K 1982 J.Geophys.Res.87 5030
[5]Pallocchia G 2013 J.Geophys.Res.118 331
[6]Tsurutani B T and Smith E J 1984 Geophys.Res.Lett.11 331
[7]Pekünlü E R,Çakırlı Ö and Özetken E 2001 Mon.Not.Roy.Astron.Soc.326 675
[8]Horne R B,Wheeler G V and Alleyne H S C K 2000 J.Geophys.Res.105 27597
[9]Stasiewicz K 2007 Plasma Phys.Control.Fusion 49 B621
[10]Miller J A 1991 Astrophys.J 376 342
[11]Miller J A 1997 Astrophys.J 491 939
[12]Ma C Y and Summers D 1998 Geophys.Res.Lett.25 4099
[13]Xiao F L et al 2015 Geophys.Res.Lett.42 7287
[14]Xiao F L et al 2015 Nat.Commun.6 8590
[15]Uzma C,Naim H and Murtaza G 2017 Phys.Plasmas 24 043701
[16]Salimullah M et al 2006 Phys.Scr.73 320
[17]Salimullah M,Zaman A and Amin M R 1996 Phys.Scr.54 530
[18]Salimullah M et al 2006 Phys.Plasmas 13 122102
[19]Naim H et al 2015 Phys.Plasmas 22 062117
[20]Bashir M F and Chen L J 2016 Phys.Plasmas 23 102107
[21]Sreekala G et al 2017 Adv.Space Res.59 2679
[22]Venugopal C et al 2013 Indian J.Phys.87 939
[23]Summers D and Thorne R M 1991 Phys.Fluids B 3 1835
[24]Mace R L and Hellberg M A 1995 Phys.Plasmas 2 2098
[25]Xiao F L et al 2008 Plasma Phys.Control.Fusion 50 062001
[26]Xiao F L,Chen L X and Li J F 2008 Plasma Phys.Control.Fusion 50 105002
[27]Xiao F L et al 2008 J.Geophys.Res.113 A05203
[28]Xiao F L 2006 Plasma Phys.Control.Fusion 48 203
[29]Hellberg M A and Mace R L 2002 Phys.Plasmas 9 1495
[30]Cattaert T,Hellberg M A and Mace R L 2007 Phys.Plasmas 14 082111
[31]Bashir M F and Murtaza G 2012 Braz.J.Phys.42 487
[32]Naim H,Bashir M F and Murtaza G 2014 Phys.Plasmas 21 102112
[33]Thorne R M and Summers D 1991 Phys.Fluids B Phys.Plasmas 3 2117
[34]Barne A 1979 Hydromagnetic waves and turbulence in the solar wind Solar System Plasma Physics Vol 1 ed E N Parker,C F Kennel and L J Lanzerotti(Amsterdam:North-Holland Publishing Co.)
[35]Schwenn R 2001 Solar wind:global properties Encyclopedia of Astronomy and Astrophysics ed P Murdin(Hampshire:Nature Publishing Group)
[36]Maksimovic M,Pierrard V and Riley P 1997 Geophys.Res.Lett.24 1151
[37]Summers D,Thorne R M and Matsumoto H 1996 Phys.Plasmas 3 2496