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 $\kappa_{zz}^{DV}=\lim_{t\rightarrow t_{\infty}}d\sigma^2/(2dt)$, the Fick's Law definition $\kappa_{zz}^{FL}=J/X$ with $X=\partial{F}/\partial{z}$, and the TGK formula definition $\kappa_{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 $\kappa_{zz}^{TGK}$ is only the approximate formula, and the Fick's Law definition $\kappa_{zz}^{FL}$ is not invariant to some DIOs. However, at least for the special condition, in this paper we show that the definition $\kappa_{zz}^{DV}$ is the invariant quantity to the kinds of the DIOs. Therefore, for spatially varying field the displacement variance definition $\kappa_{zz}^{DV}$, rather than the Fick's law definition $\kappa_{zz}^{FL}$ and TGK formula definition $\kappa_{zz}^{TGK}$, is the most appropriate definition of the SPDCs.

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. 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_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.79$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We achieve this by using a signal correction method to rectify the saturation effects of the signal readout system. The very good energy resolution from keV to MeV energies demonstrated in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

The equation $\kappa_{zz}=d\sigma^2/(2dt)$ (hereafter DCDV) is a well-known formula of energetic particles describing the relation of parallel diffusion coefficient $\kappa_{zz}$ with the parallel displacement variance $\sigma^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 $\kappa_{zz}$ as a function of $\sigma^2$ is obtained. We find that DCDV is the special case of the new equation. In addition, another equation of $\kappa_{zz}$ as a function of $\sigma^2$ corresponding to the telegraph equation is also investigated preliminarily.

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).

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 $\alpha$. A two-parameter function is obtained by taking $\alpha$ 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.

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 $\kappa_\parallel$ and its correction $T$ to $\kappa_{\parallel 0}$, where $\kappa_{\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 $\kappa_{\parallel 0}$. Furthermore, by employing the isotropic pitch-angle scattering model $D_{\mu\mu}=D(1-\mu^2)$, we find that $T'$ has the different sign as that of $T$. In this paper the spatial perpendicular diffusion coefficient $\kappa_\bot$ with the adiabatic focusing effect is also obtained.

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 $\theta_{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.

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.

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.

Using test particle simulations we study particle acceleration at highly perpendicular ($\theta_{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.

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.

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.

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.

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.

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.

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.

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}$.

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.

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.