-
The effect of solar wind on the charged particles' diffusion coefficients
Authors:
J. F. Wang,
G. Qin
Abstract:
The transport of energetic charged particles through magnetized plasmas is ubiquitous in interplanetary space and astrophysics, and the important physical quantities are the along-field and cross-field spatial diffusion coefficients of energetic charged particles. In this paper, the influence of solar wind on particle transport is investigated. Using the focusing equation, we obtain along- and cro…
▽ More
The transport of energetic charged particles through magnetized plasmas is ubiquitous in interplanetary space and astrophysics, and the important physical quantities are the along-field and cross-field spatial diffusion coefficients of energetic charged particles. In this paper, the influence of solar wind on particle transport is investigated. Using the focusing equation, we obtain along- and cross-field diffusion coefficient accounting for the solar wind effect. For different conditions, the relative importance of solar wind effect to diffusion are investigated. It is shown that when energetic charged particles are close to the sun, for along-field diffusion the solar wind effect needs to be taken into account. These results are important for studying energetic charged particle transport processes in the vicinity of the sun.
△ Less
Submitted 13 October, 2023;
originally announced October 2023.
-
Relationship of transport coefficients with statistical quantities of charged particles
Authors:
J. F. Wang,
G. Qin
Abstract:
In the previous studies, from the Fokker-Planck equation the general spatial transport equation, which contains an infinite number of spatial derivative terms $T_n=κ_{nz}\partial^n{F}/ \partial{z^n}$ with $n=1, 2, 3, \cdots$, was derived. Due to the complexity of the general equation, some simplified equations with finite spatial derivative terms have been used in astrophysical researches, e.g., t…
▽ More
In the previous studies, from the Fokker-Planck equation the general spatial transport equation, which contains an infinite number of spatial derivative terms $T_n=κ_{nz}\partial^n{F}/ \partial{z^n}$ with $n=1, 2, 3, \cdots$, was derived. Due to the complexity of the general equation, some simplified equations with finite spatial derivative terms have been used in astrophysical researches, e.g., the diffusion equation, the hyperdiffusion one, subdiffusion transport one, etc. In this paper, the simplified equations with the highest order spatial derivative terms up to the first-, second-, third-, fourth-, and fifth-order are listed, and their transport coefficient formulas are derived, respectively. We find that most of the transport coefficients are determined by the corresponding statistical quantities. In addition, we find that the well-known statistical quantities, skewness $\mathcal{S}$ and kurtosis $\mathcal{K}$, are determined by some transport coefficients. The results can help one to use different transport coefficients determined by the statistical quantities, including many that are relatively new found in this paper, to study charged particle parallel transport processes.
△ Less
Submitted 17 July, 2023; v1 submitted 23 June, 2023;
originally announced June 2023.
-
Predicting Sunspot Numbers for Solar Cycles 25 and 26
Authors:
S. -S. Wu,
G. Qin
Abstract:
The prediction of solar activity is important for advanced technologies and space activities. The peak sunspot number (SSN), which can represent the solar activity, has declined continuously in the past four solar cycles (21$-$24), and the Sun would experience a Dalton-like minimum, or even the Maunder-like minimum, if the declining trend continues in the following several cycles, so that the pred…
▽ More
The prediction of solar activity is important for advanced technologies and space activities. The peak sunspot number (SSN), which can represent the solar activity, has declined continuously in the past four solar cycles (21$-$24), and the Sun would experience a Dalton-like minimum, or even the Maunder-like minimum, if the declining trend continues in the following several cycles, so that the predictions of solar activity for cycles 25 and 26 are crucial. In Qin & Wu, 2018, ApJ, we established an SSN prediction model denoted as two-parameter modified logistic prediction (TMLP) model, which can predict the variation of SSNs in a solar cycle if the start time of the cycle has been determined. In this work, we obtain a new model denoted as TMLP-extension (TMLP-E), which can predict the solar cycle nearly two cycles in advance, so that the predictions of cycles 25 and 26 are made. It is found that the predicted solar maximum, ascent time, and cycle length are 115.1, 4.84 yr, and 11.06 yr, respectively, for cycle 25, and 107.3, 4.80 yr, and 10.97 yr, respectively, for cycle 26. The solar activities of cycles 25 and 26 are predicted to be at the same level as that of cycle 24, but will not decrease further. We therefore suggest that the cycles 24$-$26 are at a minimum of Gleissberg cycle.
△ Less
Submitted 7 October, 2021; v1 submitted 11 February, 2021;
originally announced February 2021.
-
The efficiency of electron acceleration by ICME-driven shocks
Authors:
G. Qin,
F. -J. Kong,
S. -S. Wu
Abstract:
We present a study of the acceleration efficiency of suprathermal electrons at collisionless shock waves driven by interplanetary coronal mass ejections (ICMEs), with the data analysis from both the spacecraft observations and test-particle simulations. The observations are from the 3DP/EESA instrument onboard \emph{Wind} during the 74 shock events listed in Yang et al. 2019, ApJ, and the test-par…
▽ More
We present a study of the acceleration efficiency of suprathermal electrons at collisionless shock waves driven by interplanetary coronal mass ejections (ICMEs), with the data analysis from both the spacecraft observations and test-particle simulations. The observations are from the 3DP/EESA instrument onboard \emph{Wind} during the 74 shock events listed in Yang et al. 2019, ApJ, and the test-particle simulations are carried out through 315 cases with different shock parameters. A total of seven energy channels ranging from 0.428 to 4.161 keV are selected. In the simulations, using a backward-in-time method, we calculate the average downstream flux in the $90^\circ$ pitch angle. On the other hand, the average downstream and upstream fluxes in the $90^\circ$ pitch angle can also be directly obtained from the 74 observational shock events. In addition, the variation of the event number ratio with downstream to upstream flux ratio above a threshold value in terms of the shock angle (the angle between the shock normal and upstream magnetic field), upstream Alfv$\acute{\text e}$n Mach number, and shock compression ratio is statistically obtained. It is shown from both the observations and simulations that a large shock angle, upstream Alfv$\acute{\text e}$n Mach number, and shock compression ratio can enhance the shock acceleration efficiency. Our results suggest that shock drift acceleration is more efficient in the electron acceleration by ICME-driven shocks, which confirms the findings of Yang et al. 2018.
△ Less
Submitted 23 November, 2022; v1 submitted 20 December, 2020;
originally announced December 2020.
-
Study of momentum diffusion with the effect of adiabatic focusing
Authors:
J. F. Wang,
G. Qin
Abstract:
Momentum diffusion of the energetic charged particles is an important mechanism of the transport process in astrophysics, physics of the fusion devices, and laboratory plasmas. In addition to the uniform field momentum diffusion, we obtain the modifying term due to the focusing effect of the large-scale magnetic field. After evaluating the modifying term, we find that it is determined by the sign…
▽ More
Momentum diffusion of the energetic charged particles is an important mechanism of the transport process in astrophysics, physics of the fusion devices, and laboratory plasmas. In addition to the uniform field momentum diffusion, we obtain the modifying term due to the focusing effect of the large-scale magnetic field. After evaluating the modifying term, we find that it is determined by the sign of the focusing characteristic length and the Fokker-Planck coefficients $D_{μμ}$, $D_{μp}$, $D_{pμ}$, and $D_{pp}$. It is shown that we get a new second order acceleration mechanism in this work.
△ Less
Submitted 7 December, 2021; v1 submitted 30 October, 2020;
originally announced December 2020.
-
The invariance of the diffusion coefficient with the iterative operations of charged particles' transport equation
Authors:
J. F. Wang,
G. Qin
Abstract:
The Spatial Parallel Diffusion Coefficient (SPDC) is one of the important quantities describing energetic charged particle transport. There are three different definitions for the SPDC, i.e., the Displacement Variance definition $κ_{zz}^{DV}=\lim_{t\rightarrow t_{\infty}}dσ^2/(2dt)$, the Fick's Law definition $κ_{zz}^{FL}=J/X$ with $X=\partial{F}/\partial{z}$, and the TGK formula definition…
▽ More
The Spatial Parallel Diffusion Coefficient (SPDC) is one of the important quantities describing energetic charged particle transport. There are three different definitions for the SPDC, i.e., the Displacement Variance definition $κ_{zz}^{DV}=\lim_{t\rightarrow t_{\infty}}dσ^2/(2dt)$, the Fick's Law definition $κ_{zz}^{FL}=J/X$ with $X=\partial{F}/\partial{z}$, and the TGK formula definition $κ_{zz}^{TGK}=\int_0^{\infty}dt \langle v_z(t)v_z(0) \rangle$. For constant mean magnetic field, the three different definitions of the SPDC give the same result. However, for focusing field it is demonstrated that the results of the different definitions are not the same. In this paper, from the Fokker-Planck equation we find that different methods, e.g., the general Fourier expansion and perturbation theory, can give the different Equations of the Isotropic Distribution Function (EIDFs). But it is shown that one EIDF can be transformed into another by some Derivative Iterative Operations (DIOs). If one definition of the SPDC is invariant for the DIOs, it is clear that the definition is also an invariance for different EIDFs, therewith it is an invariant quantity for the different Derivation Methods of EIDF (DMEs). For the focusing field we suggest that the TGK definition $κ_{zz}^{TGK}$ is only the approximate formula, and the Fick's Law definition $κ_{zz}^{FL}$ is not invariant to some DIOs. However, at least for the special condition, in this paper we show that the definition $κ_{zz}^{DV}$ is the invariant quantity to the kinds of the DIOs. Therefore, for spatially varying field the displacement variance definition $κ_{zz}^{DV}$, rather than the Fick's law definition $κ_{zz}^{FL}$ and TGK formula definition $κ_{zz}^{TGK}$, is the most appropriate definition of the SPDCs.
△ Less
Submitted 2 June, 2020;
originally announced June 2020.
-
Energy resolution and linearity of XENON1T in the MeV energy range
Authors:
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
V. C. Antochi,
E. Angelino,
J. Angevaare,
F. Arneodo,
D. Barge,
L. Baudis,
B. Bauermeister,
L. Bellagamba,
M. L. Benabderrahmane,
T. Berger,
P. A. Breur,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon
, et al. (113 additional authors not shown)
Abstract:
Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV due to the saturation effects. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{ββ}\simeq$ 2.46 MeV. For the XEN…
▽ More
Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV due to the saturation effects. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{ββ}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $σ/μ$ is as low as (0.80$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{ββ}$. We also present a new signal correction method to rectify the saturation effects of the signal readout system, resulting in more accurate position reconstruction and indirectly improving the energy resolution. The very good result achieved in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.
△ Less
Submitted 9 September, 2020; v1 submitted 8 March, 2020;
originally announced March 2020.
-
Diffusion coefficient with displacement variance of energetic particles with adiabatic focusing
Authors:
J. F. Wang,
G. Qin
Abstract:
The equation $κ_{zz}=dσ^2/(2dt)$ (hereafter DCDV) is a well-known formula of energetic particles describing the relation of parallel diffusion coefficient $κ_{zz}$ with the parallel displacement variance $σ^2$. In this study, we find that DCDV is only applicable to two kinds of transport equations of isotropic distribution function, one is without cross terms, the other is without convection term.…
▽ More
The equation $κ_{zz}=dσ^2/(2dt)$ (hereafter DCDV) is a well-known formula of energetic particles describing the relation of parallel diffusion coefficient $κ_{zz}$ with the parallel displacement variance $σ^2$. In this study, we find that DCDV is only applicable to two kinds of transport equations of isotropic distribution function, one is without cross terms, the other is without convection term. Here, by employing the more general transport equation, i.e., the variable coefficient differential equation derived from the Fokker-Planck equation, a new equation of $κ_{zz}$ as a function of $σ^2$ is obtained. We find that DCDV is the special case of the new equation. In addition, another equation of $κ_{zz}$ as a function of $σ^2$ corresponding to the telegraph equation is also investigated preliminarily.
△ Less
Submitted 25 September, 2019;
originally announced September 2019.
-
A Model of Sunspot Number with Modified Logistic Function
Authors:
G. Qin,
S. -S. Wu
Abstract:
Solar cycles are studied with the Version 2 monthly smoothed international sunspot number, the variations of which are found to be well represented by the modified logistic differential equation with four parameters: maximum cumulative sunspot number or total sunspot number $x_m$, initial cumulative sunspot number $x_0$, maximum emergence rate $r_0$, and asymmetry $α$. A two-parameter function is…
▽ More
Solar cycles are studied with the Version 2 monthly smoothed international sunspot number, the variations of which are found to be well represented by the modified logistic differential equation with four parameters: maximum cumulative sunspot number or total sunspot number $x_m$, initial cumulative sunspot number $x_0$, maximum emergence rate $r_0$, and asymmetry $α$. A two-parameter function is obtained by taking $α$ and $r_0$ as fixed value. In addition, it is found that $x_m$ and $x_0$ can be well determined at the start of a cycle. Therefore, a prediction model of sunspot number is established based on the two-parameter function. The prediction for cycles $4-23$ shows that the solar maximum can be predicted with average relative error being 8.8\% and maximum relative error being 22\% in cycle 15 at the start of solar cycles if solar minima are already known. The quasi-online method for determining solar minimum moment shows that we can obtain the solar minimum 14 months after the start of a cycle. Besides, our model can predict the cycle length with the average relative error being 9.5\% and maximum relative error being 22\% in cycle 4. Furthermore, we predict the sunspot number variations of cycle 24 with the relative errors of the solar maximum and ascent time being 1.4\% and 12\%, respectively, and the predicted cycle length is 11.0 (95\% confidence interval is 8.3$-$12.9) years. The comparison to the observation of cycle 24 shows that our prediction model has good effectiveness.
△ Less
Submitted 24 October, 2018; v1 submitted 10 April, 2018;
originally announced April 2018.
-
Effects of shock and turbulence properties on electron acceleration
Authors:
G. Qin,
F. -J. Kong,
L. -H. Zhang
Abstract:
Using test particle simulations we study electron acceleration at collisionless shocks with a two-component model turbulent magnetic field with slab component including dissipation range. We investigate the importance of shock normal angle $θ_{Bn}$, magnetic turbulence level $\left(b/B_0\right)^2$, and shock thickness on the acceleration efficiency of electrons. It is shown that at perpendicular s…
▽ More
Using test particle simulations we study electron acceleration at collisionless shocks with a two-component model turbulent magnetic field with slab component including dissipation range. We investigate the importance of shock normal angle $θ_{Bn}$, magnetic turbulence level $\left(b/B_0\right)^2$, and shock thickness on the acceleration efficiency of electrons. It is shown that at perpendicular shocks the electron acceleration efficiency is enhanced with the decreasing of $\left(b/B_0\right)^2$, and at $\left(b/B_0\right)^2=0.01$ the acceleration becomes significant due to strong drift electric field with long time particles staying near the shock front for shock drift acceleration (SDA). In addition, at parallel shocks the electron acceleration efficiency is increasing with the increasing of $\left(b/B_0\right)^2$, and at $\left(b/B_0\right)^2=10.0$ the acceleration is very strong due to sufficient pitch-angle scattering for first-order Fermi acceleration, as well as due to large local component of magnetic field perpendicular to shock normal angle for SDA. On the other hand, the high perpendicular shock acceleration with $\left(b/B_0\right)^2=0.01$ is stronger than the high parallel shock acceleration with ($\left(b/B_0\right)^2=10.0$), the reason might be the assumption that SDA is more efficient than first-order Fermi acceleration. Furthermore, for oblique shocks, the acceleration efficiency is small no matter the turbulence level is low or high. Moreover, for the effect of shock thickness on electron acceleration at perpendicular shocks, we show that there exists the bend-over thickness, $L_{\text{diff,b}}$. The acceleration efficiency does not change evidently if the shock thickness is much smaller than $L_{\text{diff,b}}$. However, if the shock thickness is much larger than $L_{\text{diff,b}}$, the acceleration efficiency starts to drop abruptly.
△ Less
Submitted 30 April, 2018; v1 submitted 22 February, 2018;
originally announced February 2018.
-
Parallel and perpendicular diffusion coefficients of energetic charged particles with adiabatic focusing
Authors:
J. F. Wang,
G. Qin
Abstract:
It is very important to understand stochastic diffusion of energetic charged particles in non-uniform background magnetic field in plasmas of astrophysics and fusion devices. Using different methods considering along-field adiabatic focusing effect, various authors derived parallel diffusion coefficient $κ_\parallel$ and its correction $T$ to $κ_{\parallel 0}$, where $κ_{\parallel 0}$ is the paral…
▽ More
It is very important to understand stochastic diffusion of energetic charged particles in non-uniform background magnetic field in plasmas of astrophysics and fusion devices. Using different methods considering along-field adiabatic focusing effect, various authors derived parallel diffusion coefficient $κ_\parallel$ and its correction $T$ to $κ_{\parallel 0}$, where $κ_{\parallel 0}$ is the parallel diffusion coefficient without adiabatic focusing effect. In this paper, using the improved perturbation method developed by He \& Schlickeiser and iteration process, we obtain a new correction $T'$ to $κ_{\parallel 0}$. Furthermore, by employing the isotropic pitch-angle scattering model $D_{μμ}=D(1-μ^2)$, we find that $T'$ has the different sign as that of $T$. In this paper the spatial perpendicular diffusion coefficient $κ_\bot$ with the adiabatic focusing effect is also obtained.
△ Less
Submitted 5 October, 2018; v1 submitted 25 July, 2017;
originally announced July 2017.
-
Magnetic field line random walk in two-dimensional dynamical turbulence
Authors:
J. F. Wang,
G. Qin,
Q. M. Ma,
T. Song,
S. B. Yuan
Abstract:
The field line random walk (FLRW) of magnetic turbulence is one of the important topics in plasma physics and astrophysics. In this article by using the field line tracing method mean square displacements (MSD) of FLRW is calculated in all possible length scales for pure two-dimensional turbulence with damping dynamical model. We demonstrate that in order to describe FLRW with damping dynamical mo…
▽ More
The field line random walk (FLRW) of magnetic turbulence is one of the important topics in plasma physics and astrophysics. In this article by using the field line tracing method mean square displacements (MSD) of FLRW is calculated in all possible length scales for pure two-dimensional turbulence with damping dynamical model. We demonstrate that in order to describe FLRW with damping dynamical model a new dimensionless quantity $R$ is needed to be introduced. In different length scales dimensionless MSD shows different relationship with the dimensionless quantity $R$. Although temporal effect impacts MSD of FLRW and even changes regimes of FLRW, it dose not affect the relationship between the dimensionless MSD and dimensionless quantity $R$ in all possible length scales.
△ Less
Submitted 30 June, 2017; v1 submitted 17 April, 2017;
originally announced April 2017.
-
Numerical simulations of particle acceleration at interplanetary quasi-perpendicular shocks
Authors:
F. -J. Kong,
G. Qin,
L. -H. Zhang
Abstract:
Using test particle simulations we study particle acceleration at highly perpendicular ($θ_{Bn}\geq 75^\circ$) shocks under conditions of modeling magnetic turbulence. We adopt a backward-in-time method to solve the Newton-Lorentz equation using the observed shock parameters for quasi-perpendicular interplanetary shocks, and compare the simulation results with $ACE$/EPAM observations to obtain the…
▽ More
Using test particle simulations we study particle acceleration at highly perpendicular ($θ_{Bn}\geq 75^\circ$) shocks under conditions of modeling magnetic turbulence. We adopt a backward-in-time method to solve the Newton-Lorentz equation using the observed shock parameters for quasi-perpendicular interplanetary shocks, and compare the simulation results with $ACE$/EPAM observations to obtain the injection energy and timescale of particle acceleration. With our modeling and observations we find that a large upstream speed is responsible for efficient particle acceleration. Our results also show that the quasi-perpendicular shocks are capable of accelerating thermal particles to high energies of the order of MeV for both kappa and Maxwellian upstream distributions, which may originate from the fact that in our model the local background magnetic field has a component parallel to the shock normal.
△ Less
Submitted 5 June, 2017; v1 submitted 10 April, 2017;
originally announced April 2017.
-
Perpendicular diffusion coefficient of cosmic rays in the presence of weak adiabatic focusing
Authors:
J. F. Wang,
G. Qin,
Q. M. Ma,
T. Song,
S. B. Yuan
Abstract:
The influence of adiabatic focusing on particle diffusion is an important topic in astrophysics and plasma physics. In the past several authors have explored the influence of along-field adiabatic focusing on parallel diffusion of charged energetic particles. In this paper by using the Unified NonLinear Transport (UNLT) theory developed by Shalchi (SH2010) and the method of He and Schlickeiser (HS…
▽ More
The influence of adiabatic focusing on particle diffusion is an important topic in astrophysics and plasma physics. In the past several authors have explored the influence of along-field adiabatic focusing on parallel diffusion of charged energetic particles. In this paper by using the Unified NonLinear Transport (UNLT) theory developed by Shalchi (SH2010) and the method of He and Schlickeiser (HS2014) we derive a new nonlinear perpendicular diffusion coefficient for non-uniform background magnetic field. This formula demonstrates that particle perpendicular diffusion coefficient is modified by along-field adiabatic focusing. For isotropic pitch-angle scattering and weak adiabatic focusing limit the derived perpendicular diffusion coefficient is independent of the sign of adiabatic focusing characteristic length. For two-component model we simplify the perpendicular diffusion coefficient up to second order of the power series of adiabatic focusing characteristic quantity. We find that the first order modifying factor is equal to zero and the sign of the second one is determined by the energy of particles.
△ Less
Submitted 11 July, 2017; v1 submitted 7 March, 2017;
originally announced March 2017.
-
Study of time evolution of the bend-over energy in the energetic particle spectrum at a parallel shock
Authors:
F. -J. Kong,
G. Qin,
S. -S. Wu,
L. -H. Zhang,
H. -N. Wang,
T. Chen,
P. Sun
Abstract:
Shock acceleration is considered one of the most important mechanisms for the acceleration of astrophysical energetic particles. In this work, we calculate the trajectories of a large number of test charged particles accurately in a parallel shock with magnetic turbulence. We investigate the time evolution of the accelerated-particle energy spectrum in the downstream of the shock in order to under…
▽ More
Shock acceleration is considered one of the most important mechanisms for the acceleration of astrophysical energetic particles. In this work, we calculate the trajectories of a large number of test charged particles accurately in a parallel shock with magnetic turbulence. We investigate the time evolution of the accelerated-particle energy spectrum in the downstream of the shock in order to understand the acceleration mechanism of energetic particles. From simulation results we obtain power-law energy spectra with a bend-over energy, $E_0$, increasing with time. With the particle mean acceleration time and mean momentum change during each cycle of the shock crossing from diffusive shock acceleration model (following Drury), a time-dependent differential equation for the maximum energy, $E_{acc}$, of particles accelerated at the shock, can be approximately obtained. We assume the theoretical bend-over energy as $E_{acc}$. It is found that the bend-over energy from simulations agrees well with the theoretical bend-over energy using the non-linear diffusion theory, NLGCE-F, in contrast to that using the classic quasi-linear theory (QLT).
△ Less
Submitted 15 April, 2019; v1 submitted 14 February, 2017;
originally announced February 2017.
-
Numerical simulations of solar energetic particle event timescales associated with ICMES
Authors:
S. -Y. Qi,
G. Qin,
Y. Wang
Abstract:
Recently, S.W. Kahler studied the solar energetic particle (SEP) event timescales associated with coronal mass ejections (CMEs) from spacecraft data analysis. They obtained different timescales of SEP events, such as TO, the onset time from CME launch to SEP onset, TR, the rise time from onset to half the peak intensity (0.5Ip), and TD, the duration of the SEP intensity above 0.5Ip. In this work,…
▽ More
Recently, S.W. Kahler studied the solar energetic particle (SEP) event timescales associated with coronal mass ejections (CMEs) from spacecraft data analysis. They obtained different timescales of SEP events, such as TO, the onset time from CME launch to SEP onset, TR, the rise time from onset to half the peak intensity (0.5Ip), and TD, the duration of the SEP intensity above 0.5Ip. In this work, we solve SEPs transport equation considering ICME shocks as energetic particle sources. With our modeling assumptions, our simulations show similar results to Kahler's spacecraft data analysis that the weighted average of TD increases with both CME speed and width. Besides, from our simulation results, we suggest TD is directly dependent on CME speed, but not dependent on CME width, which were not achieved from the observation data analysis.
△ Less
Submitted 2 November, 2016; v1 submitted 28 July, 2015;
originally announced July 2015.
-
Simulations of a gradual solar energetic particle event observed by Helios 1, Helios 2, and IMP 8
Authors:
G. Qin,
Y. Wang
Abstract:
In this work, a gradual solar energetic particle (SEP) event observed by multispacecraft has been simulated. The time profiles of SEP fluxes accelerated by an interplanetary shock in the three-dimensional interplanetary space are obtained by solving numerically the Fokker-Planck focused transport equation. The interplanetary shock is modeled as a moving source of energetic particles. By fitting th…
▽ More
In this work, a gradual solar energetic particle (SEP) event observed by multispacecraft has been simulated. The time profiles of SEP fluxes accelerated by an interplanetary shock in the three-dimensional interplanetary space are obtained by solving numerically the Fokker-Planck focused transport equation. The interplanetary shock is modeled as a moving source of energetic particles. By fitting the 1979/03/01 SEP fluxes observed by Helios 1, Helios 2, and IMP 8 with our simulations, we obtain the best parameters for the shock acceleration efficiency model. And we also find that the particle perpendicular diffusion coefficient with the level of $ \sim 1\% - 3\% $ of parallel diffusion coefficient at 1 AU should be included. The reservoir phenomenon is reproduced in the simulations, and the longitudinal gradient of SEP fluxes in the decay phase, which is observed by three spacecraft at different locations, is more sensitive to the shock acceleration efficiency parameters than that is to the perpendicular diffusion coefficient.
△ Less
Submitted 6 July, 2015; v1 submitted 12 May, 2015;
originally announced May 2015.
-
Improvements of the shock arrival times at the Earth model STOA
Authors:
H. -L. Liu,
G. Qin
Abstract:
Prediction of the shocks' arrival times (SATs) at the Earth is very important for space weather forecast. There is a well-known SAT model, STOA, which is widely used in the space weather forecast. However, the shock transit time from STOA model usually has a relative large error compared to the real measurements. In addition, STOA tends to yield too much `yes' prediction, which causes a large numb…
▽ More
Prediction of the shocks' arrival times (SATs) at the Earth is very important for space weather forecast. There is a well-known SAT model, STOA, which is widely used in the space weather forecast. However, the shock transit time from STOA model usually has a relative large error compared to the real measurements. In addition, STOA tends to yield too much `yes' prediction, which causes a large number of false alarms. Therefore, in this work, we work on the modification of STOA model. First, we give a new method to calculate the shock transit time by modifying the way to use the solar wind speed in STOA model. Second, we develop new criteria for deciding whether the shock will arrive at the Earth with the help of the sunspot numbers and the angle distances of the flare events. It is shown that our work can improve the SATs prediction significantly, especially the prediction of flare events without shocks arriving at the Earth.
△ Less
Submitted 21 April, 2015; v1 submitted 3 February, 2015;
originally announced February 2015.
-
Simulations of the spatial and temporal invariance in the spectra of gradual solar energetic particle events
Authors:
Yang Wang,
Gang Qin
Abstract:
The spatial and temporal invariance in the spectra of energetic particles in the gradual solar events is reproduced in the simulations. Based on a numerical solution of the focused transport equation, we obtain the intensity time profiles of solar energetic particles (SEPs) accelerated by an interplanetary shock in the three-dimensional interplanetary space. The shock is treated as a moving source…
▽ More
The spatial and temporal invariance in the spectra of energetic particles in the gradual solar events is reproduced in the simulations. Based on a numerical solution of the focused transport equation, we obtain the intensity time profiles of solar energetic particles (SEPs) accelerated by an interplanetary shock in the three-dimensional interplanetary space. The shock is treated as a moving source of energetic particles with a distribution function. The time profiles of particle flux with different energies are calculated in the ecliptic at $1$ AU. According to our model, we find that shock acceleration strength, parallel diffusion and adiabatic cooling are the main factors in forming the spatial invariance in SEP spectra, and perpendicular diffusion is a secondary factor. In addition, the temporal invariance in SEP spectra is mainly due to the effect of adiabatic cooling. Furthermore, a spectra invariant region, which agrees with observations but is different than the one suggested by Reames and co-workers, is proposed based on our simulations.
△ Less
Submitted 4 May, 2015; v1 submitted 13 January, 2015;
originally announced January 2015.
-
A numerical simulation of solar energetic particle dropouts during impulsive events
Authors:
Y. Wang,
G. Qin,
M. Zhang,
S. Dalla
Abstract:
This paper investigates the conditions for producing rapid variations of solar ener- getic particle (SEP) intensity commonly known as dropouts. In particular, we use numerical model simulations based on solving the focused transport equation in the 3-dimensional Parker interplanetary magnetic field to put constraints on the proper- ties of particle transport coefficients both in the direction perp…
▽ More
This paper investigates the conditions for producing rapid variations of solar ener- getic particle (SEP) intensity commonly known as dropouts. In particular, we use numerical model simulations based on solving the focused transport equation in the 3-dimensional Parker interplanetary magnetic field to put constraints on the proper- ties of particle transport coefficients both in the direction perpendicular and parallel to the magnetic field. To simulate the effect of field lines that may alternatively connect and disconnect an observer from a compact SEP source on the solar surface, we place several smaller SEP sources that release particles impulsively and simultaneously at separate longitudes on the solar surface. We let magnetic flux tubes filled or devoid of energetic particles pass the observer. The perpendicular particle diffusion tends to smooth out the intensity variation. Our calculations of the temporal intensity profile of 0.5 and 5 MeV protons at the Earth show that the perpendicular diffusion must be small enough while the parallel mean free path should be long in order to reproduce the phenomenon of SEP dropouts. When the parallel mean free path is a fraction of 1 AU and the observer is located at 1 AU, the perpendicular to parallel diffusion ratio must be below $10^{-5}$ if we want to see the particle flux dropping by at least several times within three hours. When the observer is located at a larger solar radial distance, the perpendicular to parallel diffusion ratio for reproducing the dropouts should be even lower than that in the case of 1 AU distance. A shorter parallel mean free path or a larger radial distance from the source to observer will cause the particles to arrive later, making the effect of perpendicular diffusion more prominent and SEP dropouts disappear.
△ Less
Submitted 22 May, 2014; v1 submitted 25 March, 2014;
originally announced March 2014.
-
The Modification of the Nonlinear Guiding Center Theory
Authors:
G. Qin,
L. -H. Zhang
Abstract:
We modify the NonLinear Guiding Center (NLGC) theory (Matthaeus et al. 2003) for perpendicular diffusion by replacing the spectral amplitude of the two-component model magnetic turbulence with the 2D component one (following Shalchi 2006), and replacing the constant $a^2$, indicating the degree particles following magnetic field line, with a variable $a^{\prime 2}$ as a function of the magnetic tu…
▽ More
We modify the NonLinear Guiding Center (NLGC) theory (Matthaeus et al. 2003) for perpendicular diffusion by replacing the spectral amplitude of the two-component model magnetic turbulence with the 2D component one (following Shalchi 2006), and replacing the constant $a^2$, indicating the degree particles following magnetic field line, with a variable $a^{\prime 2}$ as a function of the magnetic turbulence. We combine the modified model with the NonLinear PArallel (NLPA) diffusion theory (Qin 2007) to solve perpendicular and parallel diffusion coefficients simultaneously. It is shown that the new model agrees better with simulations. Furthermore, we fit the numerical results of the new model with polynomials, so that parallel and perpendicular diffusion coefficients can be calculated directly without iteration of integrations, and many numerical calculations can be reduced.
△ Less
Submitted 20 March, 2014; v1 submitted 9 January, 2014;
originally announced January 2014.
-
Study of Different Solar Cycle Variations of Solar Energetic Particles and Cosmic Rays by Despiking ACE/SIS Heavy-Ion Fluxes
Authors:
G Qin,
L. -L. Zhao
Abstract:
Cosmic Rays (CRs) include Galactic Cosmic Rays (GCRs) and Anomalous Cosmic Rays (ACRs). The CR flux data of protons and heavy-ions observed with spacecraft are often seriously contaminated by Solar Energetic Particle (SEP) events. In this work, we separate SEPs from CRs of ACE/SIS spacecraft observations with an automatic despiking algorithm, so we are able to study the different variations of SEP…
▽ More
Cosmic Rays (CRs) include Galactic Cosmic Rays (GCRs) and Anomalous Cosmic Rays (ACRs). The CR flux data of protons and heavy-ions observed with spacecraft are often seriously contaminated by Solar Energetic Particle (SEP) events. In this work, we separate SEPs from CRs of ACE/SIS spacecraft observations with an automatic despiking algorithm, so we are able to study the different variations of SEPs and CRs over a solar cycle. In particular, we study the elemental ratio, first ionization potential dependence and elemental dependence, and information entropy of SEPs and CRs. So that we can gain new insights into energetic particles' different compositions, origins, and transport processes, etc.
△ Less
Submitted 8 December, 2013;
originally announced December 2013.
-
Estimation of the release time of solar energetic particles near the Sun
Authors:
Y. Wang,
G. Qin
Abstract:
This paper investigates the onset time of Solar Energetic Particle (SEP) events with numerical simulations, and analyses the accuracy of the Velocity Dispersion Analysis (VDA) method. Using a 3-dimensional focused transport model, we calculate the fluxes of protons observed in the ecliptic at 1 AU in the energy range between 10 MeV and 80 MeV. In particular, three models are used to describe diffe…
▽ More
This paper investigates the onset time of Solar Energetic Particle (SEP) events with numerical simulations, and analyses the accuracy of the Velocity Dispersion Analysis (VDA) method. Using a 3-dimensional focused transport model, we calculate the fluxes of protons observed in the ecliptic at 1 AU in the energy range between 10 MeV and 80 MeV. In particular, three models are used to describe different SEP sources produced by flare or coronal shock, and the effects of particle perpendicular diffusion in the interplanetary space are also studied. We have the following findings: When the observer is disconnected from the source, the effects of perpendicular diffusion in the interplanetary space and particles propagating in the solar atmosphere have significant influence on the VDA results. As a result, although the VDA method is valid with impulsive source duration, low background, and weak scattering in the interplanetary space or fast diffusion in the solar atmosphere, the method is not valid with gradual source duration, high background, or strong scattering.
△ Less
Submitted 18 November, 2014; v1 submitted 29 November, 2013;
originally announced November 2013.
-
Modulation of galactic cosmic rays during the unusual solar minimum between cycles 23 and 24
Authors:
L. -L. Zhao,
G. Qin,
M. Zhang,
B. Heber
Abstract:
During the recent solar minimum between cycles 23 and 24 (solar minimum $P_{23/24}$) the intensity of Galactic Cosmic Rays (GCRs) measured at the Earth was the highest ever recorded since space age. It is the purpose of this paper to resolve the most plausible mechanism for this unusually high intensity. A GCR transport model in three-dimensional heliosphere based on a simulation of Markov stochas…
▽ More
During the recent solar minimum between cycles 23 and 24 (solar minimum $P_{23/24}$) the intensity of Galactic Cosmic Rays (GCRs) measured at the Earth was the highest ever recorded since space age. It is the purpose of this paper to resolve the most plausible mechanism for this unusually high intensity. A GCR transport model in three-dimensional heliosphere based on a simulation of Markov stochastic process is used to find the relation of cosmic ray modulation to various transport parameters, including solar wind (SW) speed, distance of heliospheric boundary, magnitude of interplanetary magnetic field (IMF) at the Earth, tilt angle of heliospheric current sheet (HCS), values of parallel and perpendicular diffusion coefficients. We calculate GCR proton energy spectra at the Earth for the last three solar minima $P_{21/22}$, $P_{22/23}$, and $P_{23/24}$, with the transport parameters obtained from observations. Besides weak IMF magnitude and slow SW speed, we find that a possible low magnetic turbulence, which increases the parallel diffusion and reduces the perpendicular diffusion in the polar direction, might be an additional possible mechanism for the high GCR intensity in the solar minimum $P_{23/24}$.
△ Less
Submitted 25 March, 2014; v1 submitted 26 October, 2013;
originally announced October 2013.