Probing the Progenitor of High-z Short-duration GRB 201221D and its Possible Bulk Acceleration in Prompt Emission

The growing observed evidence shows that the long-and short-duration gamma-ray bursts(GRBs) originate from massive star core-collapse and the merger of compact stars,respectively.GRB 201221 D is a short-duration GRB lasting~0.1 s without extended emission at high redshift z=1.046.By analyzing data observed with the Swift/BAT and Fermi/GBM,we find that a cutoff power-law model can adequately fit the spectrum with a soft E=113keV,and isotropic energy E=1.36× 10erg.In order to reveal the possible physical origin of GRB 201221 D,we adopted multi-wavelength criteria(e.g.,Amati relation,ε-parameter,amplitude parameter,local event rate density,luminosity function,and properties of the host galaxy),and find that most of the observations of GRB 201221 D favor a compact star merger origin.Moreover,we find that α is larger than 2+βin the prompt emission phase which suggests that the emission region is possibly undergoing acceleration during the prompt emission phase with a Poynting-flux-dominated jet.

Early acceleration of electrons and protons at the nonrelativistic quasiparallel shocks with different obliquity angles

The early acceleration of protons and electrons in the nonrelativistic collisionless shocks with three obliquities are investigated through 1D particle-in-cell simulations. In the simulations, the charged particles possessing a velocity of 0.2c flow towards a reflecting boundary, and the shocks with a sonic Mach number of 13.4 and an Alfven Mach number of 16.5 in the downstream shock frame are generated.In these quasi-parallel shocks with the obliquity angles θ = 15°, 30° and 45°, some of the protons and the electrons can be injected into the acceleration processes, and their downstream spectra in the momentum space show a power law tail at a time of 1.89 × 10^5ω^-1pe, where ωpe is the electron plasma frequency.Moreover, the charged particles reflected at the shock excite magnetic waves upstream of the shock. The shock drift acceleration is more prominent with a larger obliquity angle for the shocks, but the accelerated particles diffuse parallel to the shock propagation direction more easily to participate in the diffusive shock acceleration. In the early acceleration stage, more energetic protons and electrons appear in the downstream of the shock for θ = 15° compared with the other two obliquities. Moreover, in the upstream region, the spectrum of the accelerated electrons is the hardest for θnB = 45° among the three obliquities, whereas the proton spectra for θnB = 15° and 45° are similar as a result of the competition of the effectiveness of the shock drift acceleration and the diffusive shock acceleration.

The Influence of Ion-Acoustic Turbulence on the Electron Acceleration in the Reconnecting Current Sheet

Through solving the single electron equation of motion and the Fokker Planck equation including the terms of electric field strength and ion-acoustic turbulence, we study the influence of the ion-acoustic wave on the electron acceleration in turbulent reconnecting current sheets. It is shown that the ion-acoustic turbulence which causes plasma heating rather than particle acceleration should be considered. With typical parameter values, the acceleration time scale is around the order of 10-6 s, the accelerated electrons may have approximately a power-law distribution in the energy range 20-100 keV and the spectral index is about 3-10, which is basically consistent with the observed hard X-ray spectra in solar flares.

Injection Spectra of Different Species of Cosmic Rays from AMS-02, ACE-CRIS and Voyager-1

Precise measurements of energy spectra of different cosmic ray(CR) species have been obtained in recent years, by particularly the AMS-02 experiment on the International Space Station. It has been shown that apparent differences exist in different groups of the primary CRs. However, it is not straightforward to conclude that the source spectra of different particle groups are different since they will experience different propagation processes(e.g., energy losses and fragmentations) either. In this work, we study the injection spectra of different nuclear species using the measurements from Voyager-1 outside the solar system, and ACR-CRIS and AMS-02 on the top of atmosphere, in a physical framework of CR transportation. Two types of injection spectra are assumed, the broken power-law(BPL) form and the non-parametric spline interpolation form. The non-parametric form fits the data better than the BPL form, implying that potential structures beyond the constrained spectral shape of BPL may exist. For different nuclei the injection spectra are overall similar in shape but do show some differences among each other. For the non-parametric spectral form, the helium injection spectrum is the softest at low energies and the hardest at high energies. For both spectral shapes, the low-energy injection spectrum of neon is the hardest among all these species, and the carbon and oxygen spectra have more prominent bumps in 1–10 GV in the R2 d N dRpresentation.Such differences suggest the existence of differences in the sources or acceleration processes of various nuclei of CRs.

Influence of coronal holes on CMEs in causing SEP events

The issue of the influence of coronal holes （CHs） on coronal mass ejections （CMEs） in causing solar energetic particle （SEP） events is revisited. It is a continuation and extension of our previous work, in which no evident effects of CHs on CMEs in generating SEPs were found by statistically investigating 56 CME events. This result is consistent with the conclusion obtained by Kahler in 2004. We extrapolate the coronal magnetic field, define CHs as the regions consisting of only open magnetic field lines and perform a similar analysis on this issue for 76 events in total by extending the study interval to the end of 2008. Three key parameters, CH proximity, CH area and CH relative position, are involved in the analysis. The new result confirms the previous conclusion that CHs did not show any evident effect on CMEs in causing SEP events.

Broken-up spectra of the loop-top hard X-ray source during a solar limb flare

Solar hard X-rays(HXRs) appear in the form of either footpoint sources or coronal sources. Each individual source provides its own critical information on acceleration of nonthermal electrons and plasma heating. Earlier studies found that the HXR emission in some events manifests a broken-up power-law spectrum, with the break energy around a few hundred keV based on spatially-integrated spectral analysis,and it does not distinguish the contributions from individual sources. In this paper, we report on the brokenup spectra of a coronal source studied using HXR data recorded by Reuven Ramaty High Energy Solar Spectroscopic Imager(RHESSI) during the SOL2017–09–10 T16:06(GOES class X8.2) flare. The flare occurred behind the western limb and its footpoint sources were mostly occulted by the disk. We could clearly identify such broken-up spectra pertaining solely to the coronal source during the flare peak time and after. Since a significant pileup effect on the RHESSI spectra is expected for this intense solar flare, we have selected the pileup correction factor, p = 2. In this case, we found the resulting RHESSI temperature(~30MK) to be similar to the GOES soft X-ray temperature and break energies of 45–60 keV. Above the break energy, the spectrum hardens with time from spectral index of 3.4 to 2.7, and the difference in spectral indices below and above the break energy increases from 1.5 to 5 with time. However, we note that when p = 2 is assumed, a single power-law fitting is also possible with the RHESSI temperature higher than the GOES temperature by ~10MK. Possible scenarios for the broken-up spectra of the loop-top HXR source are briefly discussed.

Dependence of anomalous resistivity on bulk drift velocity of electrons in the reconnecting current sheets in solar flares

Anomalous resistivity is critical for triggering fast magnetic reconnection in the nearly collisionless coronal plasma. Its nonlinear dependence on bulk drift velocity is usually assumed in MHD simulations. However, the mechanism for the production of anomalous resistivity and its evolution is still an open question. We numerically solved the one dimension Vlasov equation with the typical solar coronal parameters and realistic mass ratios to infer the relationship between anomalous resistivity and bulk drift velocity of electrons in the reconnecting current sheets as well as its non- linear characteristics. Our principal findings are summarized as follows： 1） the relationship between the anomalous resistivity and bulk drift velocity of electrons relative to ions may be described as ηmax=0.03724（vd/ve）^5.702Ωm for vd/ve in the range of 1.4-2.0 and ηmax=0.8746（vd/ve）^1.284Ωm for vd/ve in the range of 2.5-4.5;2）if drift velocity is just slightly larger than the threshold of ion-acoustic instability, the anomalous resistivity due to the wave-particle interactions is enhanced by about five orders as compared with classic resistivity due to Coulomb collisions, With the increase of drift velocity from 1.4ve to 4.5Ve, the anomalous resistivity continues to increase 100 times; 3） in the rise phase of unstable waves, the anomalous resistivity has the same order as the one estimated from quasi-linear theory; after saturation of unstable waves, the anomalous resistivity decreases at least about one order as com- pared with its peak value; 4） considering that the final velocity of electrons ejected out of the reconnecting current sheet （RCS） decreases with the distance from the neutral point in the neutral plane, the anomalous resistivity decreases with the distance from the neutral point, which is favorable for the Petschek-like reconnection to take place.

An explanation for 13 consecutive day activities of Mrk 421

It is surprising to find an instance of migration in the peak positions of synchrotron spectral energy distribution components during the activity epochs of Markarian 421(Mrk 421),accompanying an orphan flare at the X-ray and GeV-TeVγ-ray bands.A geometric interpretation and standard shock or stochastic acceleration models of blazar emission have difficulty reproducing these observed behaviors.The present paper introduces a linear acceleration by integrating the reconnection electric field into the particle transport model for the observed behaviors of Mrk 421.We note that strong evidence for evolution in characteristic of multi-wavelength spectral energy distribution including shifting the peak frequency,accompanying an orphan flare at the X-ray and GeV-TeVγ-ray bands provides an important electrostatic acceleration diagnostic in a blazar jet.Assuming suitable model parameters,we apply the results of the simulation to the 13-day flaring event in March 2010 of Mrk 421,concentrating on the evolution of multiwavelength spectral energy distribution characteristic by shifting the peak frequency.It is clear that the ratio of the electric field and magnetic field strength plays an important role in temporal evolution of the peak frequency of synchrotron spectral energy distribution component.We suggest it is reasonable that the electrostatic acceleration is responsible for the evolution of multi-wavelength spectral energy distribution characteristic by shifting the peak frequency.Based on the model results,we assert that the peak frequency of the synchrotron spectral energy distribution component may signify a temporary characteristic of blazars,rather than a permanent one.

Numerical simulation of superhalo electrons generated by magnetic reconnection in the solar wind source region

Superhalo electrons appear to be continuously present in the interplane- tary medium, even during very quiet times, with a power-law spectrum at energies above ～2 keV. Here we numerically investigate the generation of superhalo electrons by magnetic reconnection in the solar wind source region, using magnetohydrody- namics and test particle simulations for both single X-line reconnection and multiple X-line reconnection. We find that the direct current electric field, produced in the mag- netic reconnection region, can accelerate electrons from an initial thermal energy of T ～105 K up to hundreds of keV. After acceleration, some of the accelerated elec- trons, together with the nascent solar wind flow driven by the reconnection, propagate upwards along the newly-opened magnetic field lines into interplanetary space, while the rest move downwards into the lower atmosphere. Similar to the observed superhalo electrons at 1 AU, the flux of upward-traveling accelerated electrons versus energy dis- plays a power-law distribution at ～ 2-100 keV, f（E）～ E^-δ, with a 6 of ～1.5 - 2.4. For single （multiple） X-line reconnection, the spectrum becomes harder （softer） as the anomalous resistivity parameter a （uniform resistivity η） increases. These modeling results suggest that the acceleration in the solar wind source region may contribute to superhalo electrons.

The Role of Nuclei-Nuclei Interactions in the Production of Gamma-ray Lines in Solar Flares

Dramatic extensions of experimental possibilities (spacecraft RHESSI, CORONAS-F and others) in solar gamma-ray astronomy call for urgent, detailed theoretical consideration of a set of physical problems of solar activity and solar-terrestrial relationships that earlier may have only been outlined. Here we undertake a theoretical analysis of issues related to the production of gamma-radiation in the processes of interactions of energetic (accelerated) heavy and middle nuclei with the nuclei of the solar atmosphere (the so-called i-j interactions). We also make an estimate of the contribution of these interactions to the formation of nuclear and isotopic abundances of the solar atmosphere in the range of light and rare elements. The analysis is carried out for solar flares in the wide range of their intensities. We compare our theoretical estimates with RHESSI observations for the flare of 2002 July 23. It was shown that the 24Mg gamma-ray emission in this event was produced by the newly generated Mg nuclei. With a high probability, the gamma-ray line emission of 28Si nuclei from this flare was generated by the same processes.

Resonant Heating of Ions by Parallel Propagating Alfvén Waves in Solar Coronal Holes

Resonant heating of H, O+5, and Mg+9 by parallel propagating ion cyclotron Alfven waves in solar coronal holes at a heliocentric distance is studied using the heating rate derived from the quasilinear theory. It is shown that the particle-AlfVen-wave interaction is a significant microscopic process. The temperatures of the ions are rapidly increased up to the observed order in only microseconds, which implies that simply inserting the quasilinear heating rate into the fluid/MHD energy equation to calculate the radial dependence of ion temperatures may cause errors as the time scales do not match. Different species ions are heated by Alfven waves with a power law spectrum in approximately a mass order. To heat O+5 over Mg+9 as measured by the Ultraviolet Coronagraph Spectrometer (UVCS) in the solar coronal hole at a region ≥ 1.9.R, the energy density of Alfven waves with a frequency close to the O+5-cyclotron frequency must be at least double of that at the Mg+9-cyclotron frequency. With an appropriate wave-energy spectrum, the heating of H, O+5 and Mg+9 can be consistent with the UVCS measurements in solar coronal holes at a heliocentric distance.

Ensemble Numerical Simulations of Realistic SEP Events and the Inspiration for Space Weather Awareness

The solar energetic particle(SEP)event is a kind of hazardous space weather phenomena,so its quantitative forecast is of great importance from the aspect of space environmental situation awareness.We present here a set of SEP forecast tools,which consists of three components:(1)a simple polytropic solar wind model to estimate the background solar wind conditions at the inner boundary of 0.1 AU(about 20 R⊙);(2)an ice-cream-cone model to estimate the erupted coronal mass ejection(CME)parameters;and(3)the improved Particle Acceleration and Transport in the Heliosphere(i PATH)model to calculate particle fluxes and energy spectra.By utilizing the above models,we have simulated six realistic SEP events from 2010 August 14 to 2014 September 10,and compared the simulated results to the Geostationary Operational Environmental Satellites(GOES)spacecraft observations.The results show that the simulated fluxes of>10 Me V particles agree with the observations while the simulated fluxes of>100 Me V particles are higher than the observed data.One of the possible reasons is that we have adopted a simple method in the model to calculate the injection rate of energetic particles.Furthermore,we have conducted the ensemble numerical simulations over these events and investigated the effects of different background solar wind conditions at the inner boundary on SEP events.The results imply that the initial CME density plays an important role in determining the power spectrum,while the effect of varying background solar wind temperature is not significant.Naturally,we have examined the influence of CME initial density on the numerical prediction results for virtual SEP cases with different CME ejection speeds.The result shows that the effect of initial CME density variation is inversely associated with CME speed.

Possible Contribution of Mature γ-ray Pulsars to Cosmic-ray Positrons

We restudy the possible contribution of mature gamma-ray pulsars to cosmic ray positrons based on the new version of outer gap model. In this model, the inclination angle and average properties of the outer gap are taken into account, and more mature pulsars can have the outer gap and emit high energy photons. Half of the primary particles in the outer gaps will flow back toward the star surface and emit synchrotron photons, which can produce electron/positron pairs by the cascade of pair production. Some of these pairs will escape from the light cylinder and be accelerated to relativistic energies in the pulsar wind driven by low-frequency electromagnetic waves. Using a Monte Carlo method, we obtain a sample of mature gamma-ray pulsars and then calculate the production of the positrons from these pulsars. The observed excess of cosmic positrons can be well explained by this model.

Solar Impulsive Hard X-Ray Emission and Two-Stage Electron Acceleration

Heating and acceleration of electrons in solar impulsive hard X-ray （HXR） flares are studied according to the two-stage acceleration model developed by Zhang for solar ^3Herich events. It is shown that electrostatic H-cyclotron waves can be excited at a parallel phase velocity less than about the electron thermal velocity and thus can significantly heat the electrons （up to 40 MK） through landau resonance. The preheated electrons with velocities above a threshold are further accelerated to high energies in the flare-acceleration process. The flareproduced electron spectrum is obtained and shown to be thermal at low energies and power law at high energies. In the non-thermal energy range, the spectrum can be double power law if the spectral power index is energy dependent or related. The electron energy spectrum obtained by this study agrees quantitatively with the result derived from the Reuven Ramaty High Energy Solar Spectroscopic Imager （RHESSI） HXR observations in the flare of 2002 July 23. The total flux and energy flux of electrons accelerated in the solar flare also agree with the measurements.

Probing shock geometry via the charge to mass ratio dependence of heavy ion spectra from multiple spacecraft observations of the 2013 November 4 event

In large Solar Energetic Particle （SEP） events, ions can be accelerated at coronal mass ejection （CME）-driven shocks to very high energies. The spectra of heavy ions in many large SEP events show features such as roll-overs or spectral breaks. In some events when the spectra are plotted in terms of energy/nucleon, they can be shifted relative to each other to make the spectral breaks align. The amount of shift is charge to mass ratio （Q/A） dependent and varies from event to event. This can be understood if the spectra of heavy ions are organized by the diffusion coefficients （Cohen et al. 2005）. In the work of Li et al. （2009）, the Q/A dependence of the scaling is related to shock geometry when the CME-driven shock is close to the Sun. For events where multiple in-situ spacecraft observations exist, one may expect that different spacecraft are connected to different portions of the CME-driven shock that have different shock geometries, therefore yielding different Q/A dependence. In this work, we examine one SEP event which occurred on 2013 November 4. We study the Q/A dependence of the energy scaling for heavy ion spectra using helium, oxygen and iron ions. Observations from STEREO-A, STEREO-B and ACE are examined. We find that the scalings are different for different spacecraft. We suggest that this is because ACE, STEREO- A and STEREO-B are connected to different parts of the shock that have different shock geometries. Our analysis indicates that studying the Q/A scaling of in-situ particle spectra can serve as a powerful tool to remotely examine the shock geometry for large SEP events.

Progenitor model of cosmic ray knee

The primary energy spectrum of cosmic rays exhibits a knee at about 3 PeV where a change in the spectral index occurs. Despite many efforts, the origin of such a feature in the spectrum is not satisfactorily solved yet. Here it is proposed that the steepening of the spectrum beyond the knee may be a consequence of the mass distribution of the progenitor of the cosmic ray source. The proposed speculative model can account for all the major observed features of cosmic rays without invoking any fine tuning to match flux or spectra at any energy point. The prediction of the proposed model regarding the primary composition scenario beyond the knee is quite different from most of the prevailing models of the knee, and thereby can be discriminated from precise experimental measurement of the primary composition.

Proton activity of the Sun in current solar cycle 24

We present a study of seven large solar proton events in the current solar cycle 24（from 2009 January up to the current date）. They were recorded by the GOES spacecraft with the highest proton fluxes being over 200 pfu for energies 〉10 Me V. In situ particle measurements show that：（1） The profiles of the proton fluxes are highly dependent on the locations of their solar sources, namely flares or coronal mass ejections（CMEs）, which confirms the ＂heliolongitude rules＂ associated with solar energetic particle fluxes;（2） The solar particle release（SPR） times fall in the decay phase of the flare emission, and are in accordance with the times when the CMEs travel to an average height of 7.9 solar radii; and（3） The time differences between the SPR and the flare peak are also dependent on the locations of the solar active regions. The results tend to support the scenario of proton acceleration by the CME-driven shock,even though there exists a possibility of particle acceleration at the flare site, with subsequent perpendicular diffusion of accelerated particles in the interplanetary magnetic field. We derive the integral time-of-maximum spectra of solar protons in two forms： a single power-law distribution and a power law roll-over with an exponential tail. It is found that the unique ground level enhancement that occurred in the event on 2012 May 17 displays the hardest spectrum and the largest roll-over energy which may explain why this event could extend to relativistic energies.

Synchrotron X-ray diagnostics of cutoff shape of nonthermal electron spectrum at young supernova remnants

Synchrotron X-rays can be a useful tool to investigate electron accelera- tion at young supemova remnants （SNRs）. At present, since the magnetic field con- figuration around the shocks of SNRs is uncertain, it is not clear whether electron acceleration is limited by SNR age, synchrotron cooling, or even escape from the ac- celeration region. We study whether the acceleration mechanism can be constrained by the cutoff shape of the electron spectrum around the maximum energy. We derive analytical formulae of the cutoff shape in each case where the maximum electron en- ergy is determined by SNR age, synchrotron cooling and escape from the shock. They are related to the energy dependence of the electron diffusion coefficient. Next, we discuss whether information on the cutoff shape can be provided by observations in the near future which will simply give the photon indices and the flux ratios in the soft and hard X-ray bands. We find that if the power-law index of the electron spectrum is independently determined by other observations, then we can constrain the cutoff shape by comparing theoretical predictions of the photon indices and/or the flux ratios with observed data which will be measured by NuSTAR and/or ASTRO-H. Such study is helpful in understanding the acceleration mechanism. In particular, it will supply another independent constraint on the magnetic field strength around the shocks of SNRs.

Spectral cutoffs of Fermi-LAT GRBs 080916C and 090926A

It is expected that there should be a spectral cutoff at the high energy end of emission from a prompt gamma-ray burst （GRB）, due to, e.g. γγ absorption and/or a high energy cutoff in the electron distribution. We analyze the spectral data of Fermi- LAT detected GRBs 080916C and 090926A, aiming at locating the spectral cutoff. By assuming that the prompt GRB spectrum at the high energy end is a power law with an exponential cutoff, our analysis finds that the cutoff energy Ecutoff depends on the photon index/3 and the cutoff occurs at very high energy, Ecutoff = 161＋533 GeV in GRB 080916C and Ecutoff ≥ 100 GeV （forβ ≈-2.3） in GRB 090926A. Such high energy photons, if they exist, may disfavor the synchrotron origin and need alternative generation mechanisms.

Analysis of the CME-driven shock from the SEP event that occurred on 2006 December 14

In a solar flare or coronal mass ejection （CME）, observations of the subse- quent interplanetary shock provide us with strong evidence of particle acceleration to energies of multiple MeV, even up to GeV. Diffusive shock acceleration is an efficient mechanism for particle acceleration. For investigating the shock structure, the energy injection and energy spectrum ofa CME-driven shock, we perform a dynamical Monte Carlo simulation of the CME-driven shock that occurred on 2006 December 14 using an anisotropic scattering law. The simulated results of the shock＇s fine structure, par- ticle injection, and energy spectrum are presented. We find that our simulation results give a good fit to the observations from multiple spacecraft.