EDN: TSFYCI УДК 533.951; 550.385
Definition of Lifetime due to Wave-Particle Interactions for the Pitch Angle of 90 Degrees
Sergei V. Smolin*
Siberian Federal University Krasnoyarsk, Russian Federation
Received 30.12.2022, received in revised form 15.01.2023, accepted 03.03.2023
Abstract. One of very important, but still it is not enough the investigated problems in the theoretical physics of the Earth's magnetosphere is the definition of lifetime of the charged particles due to wave-particle interactions. Therefore for the pitch angle of 90 degrees as mathematical model is offered the ordinary differential equation (ODE) for the analytical description of a perpendicular differential flux of the charged particles in the Earth's magnetosphere which depends on time t and several parameters. Using the analytical solution of the ODE, the new simple formula for definition of lifetime due to wave-particle interactions for the pitch angle of 90 degrees for different geophysical conditions is received. For calculation under this formula it is used correlated observation of enhanced electromagnetic ion cyclotron waves and dynamic evolution of ring current energetic (5-30 keV) proton flux collected by Cluster satellite near the location L = 4.5 during March 26-27, 2003, a nonstorm period (Dst >-10 nT. In addition are found the perpendicular coefficients of the particle loss function, the particle source function and the pitch angle diffusion. For the first time the modeling dependences of lifetime due to wave-particle interactions for the pitch angle of 90 degrees from the local time LT and the geomagnetic activity Kp-index are received. Mathematical statement of a problem in the general view (the system of two ODEs of the first order with two boundary conditions) is offered, when parameters depend from time, which allows to define the lifetime due to wave-particle interactions for the pitch angle of 90 degrees numerically.
Keywords: Earth's magnetosphere, pitch angle diffusion, data of the Cluster satellite, lifetime due to wave-particle interactions.
Citation: S.V. Smolin, Definition of Lifetime due to Wave-Particle Interactions for the Pitch Angle of 90 Degrees, J. Sib. Fed. Univ. Math. Phys., 2023, 16(3), 340-347. EDN: TSFYCI.
Introduction
Electromagnetic ion cyclotron (EMIC) waves play an important role in precipitation losses of ring current ions [1] via wave-particle interactions during geomagnetic storms in the Earth's magnetosphere. Pitch angle diffusion [2] induced by EMIC waves is considered as an efficient mechanism responsible for the ring current decay during active geomagnetic periods. Variation of radiation belt dynamics is considered as the contribution from wave-particle interactions [3,4] and from drift resonance associated with enhanced ultra low frequency waves [5-7]. In more detail for the review of the literature look [1-7].
Therefore a lifetime due to wave-particle interactions has very much great value at carrying out of calculations and mathematical modeling. But the works, devoted only to definition of the lifetime due to wave-particle interactions, practically it is not enough. In this connection the purpose of the work is to propose a new variant (a new formula) of definition of the lifetime
* smolinsv@inbox.ru © Siberian Federal University. All rights reserved
due to wave-particle interactions for the 90° pitch angle for different geophysical conditions and corresponding mathematical statement of a problem in the general view, which allows to define Twp± numerically.
1. The mathematical model
The following non-stationary differential equation is used [8, 9] for phase space density, which describes the pitch angle diffusion and losses due to wave-particle interactions in the Earth's magnetosphere for a range of pitch angles from 0° up to 180°:
df 1 d ( df sin2 a cos a dL \ f 2
Ism aDaa —--—--—f ) - T--+ S± sin af, (1)
dt sin a da \ da 2L dt J Twp
where f is the phase space density (or distribution function); t is the time; a is the local pitch angle; Daa is the coefficient of pitch angle diffusion; L is the Mcllwain parameter; dL/dt is the radial rate; Twp is the average lifetime due to wave-particle interactions; S± is the perpendicular coefficient of the particle source function.
Equation (1) describes the pitch angle diffusion in the velocity space with losses due to wave-particle interactions. The loss function is conditioned by the fall of charged particles in the so-called "loss cone" as a result of wave-particle interactions. The particle source function can be related, for example, to charged particles that move from the tail of the magnetosphere toward the Earth when affected by magnetospheric convection.
The mathematical model of pitch angle diffusion of charged particles in the Earth's magnetosphere represented by (1) takes into account three physical mechanisms [4]. First, the wave-particle interactions are considered with Twp. Second, the physical mechanism of injection and particle drift is taken into account through the radial drift velocity dL/dt in equation (1). Third, since the electric field potential depends on the geomagnetic activity index Kp [8,9], we take into account the splitting effect of drift shells of the electric field on the pitch angle distribution of charged particles.
In the further the following relationship will be used between a differential flux of particles j and phase space density f (or function of distribution) j = 2mEf, where m is the mass of the charged particle (in the given work of a proton), E is the energy of the charged particle.
Using (1) [10], for the 90° pitch angle as mathematical model is offered the ordinary differential equation (ODE) for the analytical description of a perpendicular differential flux of the charged particles in the Earth's magnetosphere which depends on time and several parameters:
+ (- .Ldt + ) j. = 0, (2)
dj± i__1_ dL + Y^ - 7±Q
dt V 2L dt y±Q (Y±Q + 2) TwP±
where j± is the perpendicular differential flux of the charged particles, y±0 is the well-known (when j (a) = sin7(a) a, a is the local pitch angle) parameter of the charged particle pitch angle distribution (or the pitch angle distribution anisotropy index) taken for the 90° pitch angle at t = 0, j± is the average parameter of the charged particle pitch angle distribution on a time interval of calculation t = tend (it is supposed, that j± (t) « const), Twp± is the lifetime due to wave-particle interactions for the 90° pitch angle.
For calculation on experimental data of charged particle pitch angle distributions it are found j±Q (t = 0), y±q (t = 0) and j±end (t = tend), J±end (t = tend) and then
Y±0 + Y^end
Y± «-2-• (3)
Thus the found average value y± is approximately fair for rather small time intervals and for nonstorm time conditions.
The equation (2) can be solved analytically under the certain geophysical conditions and on a time interval approximately no more than three hours (when a geomagnetic activity index Kp = const) or on a greater time interval, when Kp « const.
We shall take the following approached equality, using the bounce-averaged radial drift velocity of charged particles in the Earth's magnetosphere (measured in 1/s) [8-10],
dL dL\ ^62 t4 , cos 6
— «( — ) = — u—L4 cos 6 =--, (4)
dt \dt / 6o 24 • 60 • 60 • 92' w
where 6 is the azimuthal angle (local time LT = 0 h at midnight) or geomagnetic eastern longitude in the magnetic equator plane, U is the Earth's angular velocity, 6o = 92 kV, and the dependence of 62 (measured in kV) on geomagnetic activity, i.e. on the Kp-index, is determined from the formula [11]
62 «_0045_3. (5)
(1 - 0.16Kp + 0.01Kp2)3
Believing all parameters in the equation (2) in view of (3)-(5) are constants (or approximately constants), we find the analytical solution of the equation (2) exact (or approached):
. f f n02L3 cos6 . Y, - Y,o A A (c,,
j±(t) = ^ exHH^94880^ + 7±o (7,0 + 2) TwpJ • V . (6)
Here, j,o is the perpendicular differential flux of the charged particles at t = 0, and j, (t) is the perpendicular differential flux of the charged particles at the moment of time t.
From the equation (6) we find the new simple formula for definition of the lifetime due to wave-particle interactions for the 90° pitch angle (measured in s) for different geophysical conditions
T = (Y, - 7±o) •t A j,o nfaL3 cos 6 \- (7)
lwp^ = Y,o (y,o + 2) \ j, (t) 7948800 V . U)
Always the condition Twp, > 0 should be satisfied.
The definition of the lifetime due to wave-particle interactions for the 90° pitch angle (7) also allows to find such very important for modeling quantitative characteristics as the perpendicular coefficient of the particle loss function [8, 9]
L, = T—, (8)
T wp,
the perpendicular coefficient of the particle source function [8, 9]
S = (Y,o + 3)
= 7-ro^r-' (9)
(Y,o + 2) Twp,
and the perpendicular coefficient of pitch angle diffusion is offered to be determined [8, 9] as follows
D, =-(-+2) T . (10)
Y,o(Y,o + 2) Twp,
If all parameters in the equation (2) in view of (3)-(5) are not constants, then the equation (2) in the general view should be written down so
dj, + (n62 (t) L (t)cos 6 (t) + Y, (t) - Y,o \. =0 (11)
dt V 7948800 Y,o (y,o + 2) Twp, ]j, . ()
j,
The equation (11) allows to enter into calculation a trajectory of the satellite, for example, in a parametrical form L (t), $ (t), change of potential of a magnetospheric convection electric field $2 (Kp (t)) and change of a anisotropy index of the charged particle pitch angle distribution Y± (t). But in this equation Twp± = const, therefore the definition of the lifetime due to wave-particle interactions for the 90° pitch angle can be offered as follows.
The problem connected with ODE of the first order (11) can be considered as a boundary value problem with undetermined parameter Twp±, if to set two boundary conditions for a perpendicular differential flux: j±Q and j±end. By introducing the parameter Twp± as dependent variable, the problem can be written as a boundary value problem in standard form. To find the solution j± (t) and the value of Twp±, just add the equation
This system of ODEs of the first order (11) and (12) solves a problem with two boundary conditions for definition j± (t) and Twp±.
Thus, at use of the offered system of the equations (11), (12) with two boundary conditions j±Q and j±end, when parameters L (t), $ (t), $2 (Kp (t)), y± (t) depend from time, the lifetime due to wave-particle interactions for the 90° pitch angle Twp± is defined numerically.
2. Results of calculations
Further it will be used correlated observation of enhanced electromagnetic ion cyclotron (EMIC) waves and dynamic evolution of ring current proton flux collected by Cluster satellite near the location L = 4.5 during March 26-27, 2003, a nonstorm period (Dst > -10 nT) [1]. It is shown, that energetic (5-30 keV) proton fluxes are found to drop rapidly (e.g., a half hour) at lower pitch angles, corresponding to intensified EMIC wave activities.
At the moment of time 00:00 RT = 23:56 UT on 26 March, 2003 (RT is the "running", current time of modeling) Kp-index of geomagnetic activity was equaled 2.66 or 3- and further within almost 3 hours remained constant. As Kp-index is the planetary three-hour index. Time of modeling for comparison with experimental data has been taken 30 minutes, i.e. 00:30 RT = 00:26 UT on 27 March, 2003.
The initial condition in the moment of time 00:00 RT = 23:56 UT on 26 March, 2003 we shall take from work [1], where the pitch angle distribution measured on the Cluster satellite for energy E = 17.1 keV, L = 4.17, LT = 22.58 h (nonstorm conditions) is presented. In this work the initial pitch angle distribution (Fig. 1) is designated squares [12]. In the same figure experimental data (circles) are shown for energy E = 17.1 keV also in the moment of time 00:30 RT = 00:26 UT on 27 March, 2003 (again nonstorm conditions).
For the approached analytical description of experimental data it is used the pitch angle distribution
For example, the good consent turns out, when the initial perpendicular differential flux of protons is equal j±Q « 560360 (cm2 s sr keV)-1, and the parameter of pitch angle distribution Y = Y±Q = 0.5157.
For calculation we shall take following data on the interval t = (0-1800) s, tend = 1800 s (Fig. 1): j±Q(t = 0 s) = 5^6036 • 105 (cm2 s sr keV)-1, j^end(tend = 1800 s) = 21165 • 105 (cm2 s sr keV)-1, Kp = 2^66 or 3-, L = 4.17, $ = LT = 22.58 h, Y±o(t = 0 s) = 0.5157, y± = 0.580, E = 17.1 keV. Then using (5) and (7) the lifetime due to wave-particle interactions for the 90° pitch angle is received Twp± = 92.5356 s. If necessary it is possible to present dependence (t) (6) on the interval t = (0-1800) s.
(12)
j (a) = sin7 a, y = Y± = const•
(13)
L = 4.17; LT = 22.58 h; E = 17.1 keV; k = 2.8 10-6 10 -,-,-,-,-,—
= 0 _____—-□—e- -B-
n □ n -
□ □
-er-
o "O-ÇT O
cv t = 30 min ,o
Jo
0 0.5 1 1.5 2 2.5 3
Pitch angle, rad
Fig. 1. The pitch angle distributions of protons measured of the Cluster satellite for energy E = 17.1 keV. Squares and circles specify a differential flux for 00:00 RT = 23:56 UT on 26 March, 2003 and for 00:30 RT = 00:26 UT on 27 March, 2003 (nonstorm conditions) respectively. Blue and red lines specify a modeling differential flux for 00:00 RT and 00:30 RT (L = 4.17, LT = 22.58 h, k = 2.8 • 1Q-6 [12]) respectively
Accuracy of definition of the lifetime due to wave-particle interactions for the 90° pitch angle Twp± is connected with accuracy of definition of the data entering into the formula (7).
In addition in this case we receive following concrete results: the perpendicular coefficient of the particle loss function (8) L± = 0.0108 1/s, the perpendicular coefficient of the particle source function (9) S± = 0.0151 1/s and the perpendicular coefficient of pitch angle diffusion (10) Dl = 0.0083 1/s.
Using the previous data and assuming, that they remain constants or approximately constants, under the formula (7) we can find a modeling dependence Twp^ on azimuthal angle (local time LT = 0 h at midnight) or geomagnetic eastern longitude in the magnetic equator plane $ = LT (Fig. 2).
As it is well visible (Fig. 2), the presented modeling dependence corresponds to the law cos $ (7).
Again using the previous data and assuming, that they remain constants or approximately constants, under the formula (7) we can find a modeling dependence Twp± from the geomagnetic activity Kp-index (Fig. 3).
As a result (Fig. 3) it is received nonlinear dependence the lifetime due to wave-particle interactions for the 90° pitch angle from the geomagnetic activity Kp-index with the maximal value Twp± at Kp = 8.
Dependences Twp^ from the magnetic local time LT and the geomagnetic activity Kp-index are received for the first time.
Thus for the offered variant (the formula (7)) it is enough to have only the experimental data on two charged particle pitch angle distributions to find the lifetime due to wave-particle interactions for the 90° pitch angle, and then it is possible to find the perpendicular coefficient of the particle loss function (8), the perpendicular coefficient of the particle source function (9)
Protons, E = 17.1 keV, L = 4.17, Kp = 3-
0 3 6 3 12 15 18 21 24 LT, h
Fig. 2. The modeling dependence the lifetime due to wave-particle interactions for the 90° pitch angle Twp, on 6 or the local time LT
Protons, E = 17.1 keV, L = 4.17, LT = 22.58 h
-1-1-1-1-1-1-1-1-
0123456789
Kp-indei;
Fig. 3. The modeling dependence the lifetime due to wave-particle interactions for the 90° pitch angle Twp, from the geomagnetic activity Kp-index
and the perpendicular coefficient of pitch angle diffusion (10) for the subsequent mathematical modeling.
Further we use other mathematical statement of a problem, namely system of ODEs of the first order (11) and (12) with two boundary conditions. For an example of such calculation we shall take the previous data on a time interval (0-1800) s and we shall find Twp,. As a result as well as should be numerically we receive (for example using the Mathematica software program) the same result Twp, = 92.5356 s.
Thus in the further by means of system of two equations (11) and (12), when parameters depend from time on experimental data, it will be possible to receive more exact dependences Twp, from the local time LT and the geomagnetic activity Kp-index.
Conclusion
1. For the 90° pitch angle as the mathematical model the ordinary differential equation for the analytical description of a perpendicular differential flux of the charged particles in the Earth's magnetosphere which depends on time and several parameters is offered.
2. The new simple formula for definition of the lifetime due to wave-particle interactions for the 90° pitch angle Twp, for different geophysical conditions is found.
3. The definition of the lifetime due to wave-particle interactions for the 90° pitch angle also allows to find such very important quantitative characteristics as the perpendicular coefficient of the particle loss function, the perpendicular coefficient of the particle source function and the perpendicular coefficient of pitch angle diffusion for the subsequent mathematical modeling.
4. For concrete calculations data of energetic (5-30 keV) proton pitch angle distributions collected by Cluster satellite near the location L = 4.5 during March 26-27, 2003, a nonstorm period (Dst > -10 nT) are used.
5. For the first time the modeling dependences Twp, from the local time LT and the geomagnetic activity Kp-index are received.
6. For the offered variant (the new simple formula (7)) it is enough to have only the experimental data on two charged particle pitch angle distributions to find the lifetime due to wave-particle interactions for the 90° pitch angle analytically.
7. Mathematical statement of a problem in the general view (the system of two ODEs of the first order (11) and (12) with two boundary conditions) is offered, when parameters depend from time, which allows to define Twp, numerically.
8. By means of system of two ODEs, when parameters depend from time on experimental data, it will be possible to receive more exact dependences Twp, from the local time LT and the geomagnetic activity Kp-index.
References
[1] F.Xiao, C.Yang, Q.Zhou, Z.He, Y.He, X.Zhou, L.Tang, Nonstorm time scattering of ring current protons by electromagnetic ion cyclotron waves, J. Geophys. Res., 117(2012), no. 8, A08204. DOI: 10.1029/2012JA017922
[2] Y.S.Miyoshi, V.K.Jordanova, A.Morioka, M.F.Thomsen, G.D.Reeves, D.S.Evans, J.C.Green, Observations and modeling of energetic electron dynamics during the October 2001 storm, J. Geophys. Res., 111(2006), no. A11, 502-513. DOI: 10.1029/2005JA011351
[3] D.Summers, R.M.Thorne, Relativistic electron pitch angle scattering by electromagnetic ion cyclotron waves during geomagnetic storms, J. Geophys. Res., 108(2003), no. A4, 1143-1154. DOI: 10.1029/2002JA009489
[4] D.G.Sibeck, R.W.McEntire, A.T.Y.Lui, R.E.Lopez, S.M.Krimigis, Magnetic field drift shell splitting: cause of unusual dayside particle pitch angle distributions during storms and sub-storms, J. Geophys. Res., 92(1987), no. A12, 13485-13497. DOI: 10.1029/JA092iA12p13485
[5] H.Zhao, X.Li, Modeling energetic electron penetration into the slot region and inner radiation belt, J. Geophys. Res., 118(2013), 6936-6945. DOI: 10.1002/2013JA019240
[6] H.Zhao, X.Li, D.N.Baker, S.G.Claudepierre, J.F.Fennell, J.B.Blake, et al., Ring current electron dynamics during geomagnetic storms based on the Van Allen Probes measurements, J. Geophys. Res., 121(2016), 3333-3346. DOI: 10.1002/2016JA022358
[7] H.Zhao, R.H.W.Friedel, Y.Chen, G.D.Reeves, D.N.Baker, X.Li, et al., An empirical model of radiation belt electron pitch angle distributions based on Van Allen probes measurements, J. Geophys. Res., 123(2018), 3493-3511. DOI: 10.1029/2018JA025277
[8] S.V.Smolin, Modeling of pitch angle distribution on the dayside of the Earth's magnetosphere, Journal of Siberian Federal University. Mathematics & Physics, 5(2012), no. 2, 269-275 (in Russian).
[9] S.V.Smolin, Modeling the pitch angle distribution on the nightside of the Earth's magnetosphere, Geomagnetism and Aeronomy, 55(2015), no. 2, 166-173.
DOI: 10.1134/S0016793215020152
[10] S.V.Smolin, The analytical description of the Earth's ring current proton flux for the pitch angle of 90 degrees, Space, Time and Fundamental Interactions, 27(2019), no. 2, 70-74. DOI: 10.17238/issn2226-8812.2019.2.70-74
[11] A.Nishida, Geomagnetic diagnosis of the magnetosphere, Springer-Verlag, New York, 1978.
[12] S.V.Smolin, Ring current proton dynamics driven by wave-particle interactions during a nonstorm period, Journal of Siberian Federal University. Mathematics & Physics, 14(2021), no. 1, 98-104. DOI: 10.17516/1997-1397-2021-14-1-98-104
Определение времени жизни вследствие взаимодействий волна-частица для питч-угла 90 градусов
Сергей В. Смолин
Сибирский федеральный университет Красноярск, Российская Федерация
Аннотация. Одной из очень важных, но еще недостаточно исследованных задач в теоретической физике магнитосферы Земли является определение времени жизни заряженных частиц вследствие взаимодействий волна-частица. Поэтому для питч-угла 90 градусов как математическая модель предлагается обыкновенное дифференциальное уравнение (ОДУ) для аналитического описания перпендикулярного дифференциального потока заряженных частиц в магнитосфере Земли, которое зависит от времени и нескольких параметров. Используя аналитическое решение ОДУ, получена новая простая формула для определения времени жизни вследствие взаимодействий волна-частица для питч-угла 90 градусов для разных геофизических условий. Для расчета по этой формуле используется коррелированное наблюдение усиленных электромагнитных ионно-циклотронных волн и динамической эволюции потока энергичных (5-30 кэВ) протонов кольцевого тока, собранное спутником Cluster около положения L = 4.5 в течение 26-27 марта 2003 г. в магнитоспокойный период (Dst > -10 nT). В дополнение найдены перпендикулярные коэффициенты функции потерь частиц, функции источника частиц и питч-угловой диффузии. Впервые получены модельные зависимости времени жизни вследствие взаимодействий волна-частица для питч-угла 90 градусов от местного времени LT и Кр-индекса геомагнитной активности. Предложена математическая постановка задачи в общем виде (система двух ОДУ первого порядка с двумя граничными условиями), когда параметры зависят от времени, которая позволяет определять время жизни вследствие взаимодействий волна-частица для питч-угла 90 градусов численно.
Ключевые слова: магнитосфера Земли, питч-угловая диффузия, данные спутника Cluster, время жизни вследствие взаимодействий волна-частица.