Context. Solar Energetic Particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce Ground Level Enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation only parallel to the magnetic field lines of interplanetary space, i.e. it is spatially 1D. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the ~100 MeV-1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons through the Interplanetary Magnetic Field (IMF) using a 3D approach, to assess the effect of the Heliospheric Current Sheet (HCS) and drifts associated to the gradient and curvature of the Parker spiral. The latter are influenced by the IMF polarity. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used to compare with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings are strongly influenced by 3D effects, with significant differences between periods of A+ and A- polarities. The decay time constant of 1 AU intensity profiles varies depending on the polarity and position of the observer, and it is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, in many cases resulting in a roll-over in the fluence spectrum.

The Parker Solar Probe was launched on 2018 August 12 and completed its second orbit on 2019 June 19 with perihelion of 35.7 solar radii. During this time, the Energetic particle Instrument-Hi (EPI-Hi, one of the two energetic particle instruments comprising the Integrated Science Investigation of the Sun, ISOIS) measured seven proton intensity increases associated with stream interaction regions (SIRs), two of which appear to be occurring in the same region corotating with the Sun. The events are relatively weak, with observed proton spectra extending to only a few MeV and lasting for a few days. The proton spectra are best characterized by power laws with indices ranging from -4.3 to -6.5, generally softer than events associated with SIRs observed at 1 au and beyond. Helium spectra were also obtained with similar indices, allowing He/H abundance ratios to be calculated for each event. We find values of 0.016-0.031, which are consistent with ratios obtained previously for corotating interaction region events with fast solar wind < 600 km s-1. Using the observed solar wind data combined with solar wind simulations, we study the solar wind structures associated with these events and identify additional spacecraft near 1 au appropriately positioned to observe the same structures after some corotation. Examination of the energetic particle observations from these spacecraft yields two events that may correspond to the energetic particle increases seen by EPI-Hi earlier.

Observations at 1 au have confirmed that enhancements in measured energetic particle fluxes are statistically associated with "rough" magnetic fields, i.e., fields having atypically large spatial derivatives or increments, as measured by the Partial Variance of Increments (PVI) method. One way to interpret this observation is as an association of the energetic particles with trapping or channeling within magnetic flux tubes, possibly near their boundaries. However, it remains unclear whether this association is a transport or local effect; i.e., the particles might have been energized at a distant location, perhaps by shocks or reconnection, or they might experience local energization or re-acceleration. The Parker Solar Probe (PSP), even in its first two orbits, offers a unique opportunity to study this statistical correlation closer to the corona. As a first step, we analyze the separate correlation properties of the energetic particles measured by the \isois instruments during the first solar encounter. The distribution of time intervals between a specific type of event, i.e., the waiting time, can indicate the nature of the underlying process. We find that the \isois observations show a power-law distribution of waiting times, indicating a correlated (non-Poisson) distribution. Analysis of low-energy \isois data suggests that the results are consistent with the 1 au studies, although we find hints of some unexpected behavior. A more complete understanding of these statistical distributions will provide valuable insights into the origin and propagation of solar energetic particles, a picture that should become clear with future PSP orbits.

We draw a comparison between a solar energetic particle event associated with the release of a slow coronal mass ejection close to the sun, and the energetic particle population produced in high current density field-aligned current structures associated with auroral phenomena in planetary magnetospheres. We suggest that this process is common in CME development and lift-off in the corona, and may account for the electron populations that generate Type III radio bursts, as well as for the prompt energetic ion and electron populations typically observed in interplanetary space.

A solar energetic particle event was detected by the Integrated Science Investigation of the Sun (ISOIS) instrument suite on Parker Solar Probe (PSP) on 2019 April 4 when the spacecraft was inside of 0.17 au and less than 1 day before its second perihelion, providing an opportunity to study solar particle acceleration and transport unprecedentedly close to the source. The event was very small, with peak 1 MeV proton intensities of ~0.3 particles (cm^2 sr s MeV)^-1, and was undetectable above background levels at energies above 10 MeV or in particle detectors at 1 au. It was strongly anisotropic, with intensities flowing outward from the Sun up to 30 times greater than those flowing inward persisting throughout the event. Temporal association between particle increases and small brightness surges in the extreme-ultraviolet observed by the Solar TErrestrial RElations Observatory, which were also accompanied by type III radio emission seen by the Electromagnetic Fields Investigation on PSP, indicates that the source of this event was an active region nearly 80 degrees east of the nominal PSP magnetic footpoint. This suggests that the field lines expanded over a wide longitudinal range between the active region in the photosphere and the corona.

A series of solar energetic particle (SEP) events were observed at Parker Solar Probe (PSP) by the Integrated Science Investigation of the Sun (\ISOIS) during the period from April 18, 2019 through April 24, 2019. The PSP spacecraft was located near 0.48 au from the Sun on Parker spiral field lines that projected out to 1 au within $\sim 25^\circ$ of near Earth spacecraft. These SEP events, though small compared to historically large SEP events, were amongst the largest observed thus far in the PSP mission and provide critical information about the space environment inside 1 au during SEP events. During this period the Sun released multiple coronal mass ejections (CMEs). One of these CMEs observed was initiated on April 20, 2019 at 01:25 UTC, and the interplanetary CME (ICME) propagated out and passed over the PSP spacecraft. Observations by the Electromagnetic Fields Investigation (FIELDS) show that the magnetic field structure was mostly radial throughout the passage of the compression region and the plasma that followed, indicating that PSP did not directly observe a flux rope internal to the ICME, consistent with the location of PSP on the ICME flank. Analysis using relativistic electrons observed near Earth by the Electron, Proton and Alpha Monitor (EPAM) on the Advanced Composition Explorer (ACE) demonstrates the presence of electron seed populations (40--300 keV) during the events observed. The energy spectrum of the \ISOIS~ observed proton seed population below 1 MeV is close to the limit of possible stationary state plasma distributions out of equilibrium. \ISOIS~ observations reveal the \revise{enhancement} of seed populations during the passage of the ICME, which \revise{likely indicates a key part} of the pre-acceleration process that occurs close to the Sun.

The Integrated Science Investigations of the Sun (IS$\odot$IS) instrument suite includes two Energetic Particle instruments: EPI-Hi, designed to measure ions from ~1-200 MeV/nuc, and EPI-Lo, designed to measure ions from ~20 keV/nuc to ~15 MeV/nuc. We present an analysis of eight energetic proton events observed across the energy range of both instruments during PSP's first two orbits in order to examine their combined energy spectra. Background corrections are applied to help resolve spectral breaks between the two instruments and are shown to be effective. In doing so we demonstrate that, even in the early stages of calibration, IS$\odot$IS is capable of producing reliable spectral observations across broad energy ranges. In addition to making groundbreaking measurements very near the Sun, IS$\odot$IS also characterizes energetic particle populations over a range of heliocentric distances inside 1 au. During the first two orbits, IS$\odot$IS observed energetic particle events from a single corotating interaction region (CIR) at three different distances from the Sun. The events are separated by two Carrington rotations and just 0.11 au in distance, however the relationship shown between proton intensities and proximity of the spacecraft to the source region shows evidence of the importance of transport effects on observations of energetic particles from CIRs. Future IS$\odot$IS observations of similar events over larger distances will help disentangle the effects of CIR-related acceleration and transport. We apply similar spectral analyses to the remaining five events, including four that are likely related to stream interaction regions (SIRs) and one solar energetic particle (SEP) event.

Nonthermal loop-top sources in solar flares are the most prominent observational signature that suggests energy release and particle acceleration in the solar corona. Although several scenarios for particle acceleration have been proposed, the origin of the loop-top sources remains unclear. Here we present a model that combines a large-scale magnetohydrodynamic simulation of a two-ribbon flare with a particle acceleration and transport model for investigating electron acceleration by a fast-mode termination shock at the looptop. Our model provides spatially resolved electron distribution that evolves in response to the dynamic flare geometry. We find a concave-downward magnetic structure located below the flare termination shock, induced by the fast reconnection downflows. It acts as a magnetic trap to confine the electrons at the looptop for an extended period of time. The electrons are energized significantly as they cross the shock front, and eventually build up a power-law energy spectrum extending to hundreds of keV. We suggest that this particle acceleration and transport scenario driven by a flare termination shock is a viable interpretation for the observed nonthermal loop-top sources.

This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

We present numerical modelling of particle acceleration at coronal shocks propagating through a streamer-like magnetic field by solving the Parker transport equation with spatial diffusion both along and across the magnetic field. We show that the location on the shock where the high-energy particle intensity is the largest, depends on the energy of the particles and on time. The acceleration of particles to more than 100 MeV mainly occurs in the shock-streamer interaction region, due to perpendicular shock geometry and the trapping effect of closed magnetic fields. A comparison of the particle spectra to that in a radial magnetic field shows that the intensity at 100 MeV (200 MeV) is enhanced by more than one order (two orders) of magnitude. This indicates that the streamer-like magnetic field can be an important factor in producing large solar energetic particle events. We also show that the energy spectrum integrated over the simulation domain consists of two different power laws. Further analysis suggests that it may be a mixture of two distinct populations accelerated in the streamer and open field regions, where the acceleration rate differs substantially. Our calculations also show that the particle spectra are affected considerably by a number of parameters, such as the streamer tilt angle, particle spatial diffusion coefficient, and shock compression ratio. While the low-energy spectra agree well with standard diffusive shock acceleration theory, the break energy ranges from $\sim$1 MeV to $\sim$90 MeV and the high-energy spectra can extend to $\sim$1 GeV with a slope of $\sim$2-3.

The effects of turbulence in the very local interstellar medium (VLISM) have been proposed by Giacalone & Jokipii (2015) to be important in determining the structure of the Interstellar Boundary Explorer (IBEX) ribbon via particle trapping by magnetic mirroring. We further explore this effect by simulating the motion of charged particles in a turbulent magnetic field superposed on a large-scale mean field, which we have considered to be either spatially-uniform or a mean field derived from a 3D MHD simulation. We find that the ribbon is not double-peaked, in contrast to Giacalone & Jokipii (2015). However, the magnetic mirror force still plays an important role in trapping particles. Furthermore, the ribbon$'$s thickness is considerably larger if the large-scale mean field is draped around the heliosphere. Voyager 1 observations in the VLISM show a turbulent field component that is stronger than previously thought, which we test in our simulation. We find that the inclusion of turbulent fluctuations at scales ${\gtrsim}$100 au and power consistent with Voyager 1 observations produces a ribbon whose large-scale structure is inconsistent with IBEX observations. However, restricting the fluctuations to ${\sim}$10 au or smaller produces a smoother ribbon structure similar to IBEX observations. Different turbulence realizations produce different small-scale features ${\lesssim}10{\deg}$ in the ribbon, but its large-scale structure is robust if the maximum fluctuation size is ${\sim}$10 au. This suggests that the magnetic field structure at scales ${\lesssim}$10 au is determined by the heliosphere$'$s interaction with the VLISM and cannot entirely be represented by homogeneous interstellar turbulence.

Little is known about the origin of the high-energy and sustained emission from solar Long-Duration Gamma-Ray Flares (LDGRFs), identified with the Compton Gamma Ray Observatory (CGRO), the Solar Maximum Mission (SMM), and now Fermi. Though Fermi/Large Area Space Telescope (LAT) has identified dozens of flares with LDGRF signature, the nature of this phenomenon has been a challenge to explain both due to the extreme energies and long durations. The highest-energy emission has generally been attributed to pion production from the interaction of >300 MeV protons with the ambient matter. The extended duration suggests that particle acceleration occurs over large volumes extending high in the corona, either from stochastic acceleration within large coronal loops or from back precipitation from coronal mass ejection driven shocks. It is possible to test these models by making direct comparison between the properties of the accelerated ion population producing the gamma-ray emission derived from the Fermi/LAT observations, and the characteristics of solar energetic particles (SEPs) measured by the Payload for Matter-Antimatter Exploration and Light Nuclei Astrophysics (PAMELA) spacecraft in the energy range corresponding to the pion-related emission detected with Fermi. For fourteen of these events we compare the two populations -- SEPs in space and the interacting particles at the Sun -- and discuss the implications in terms of potential sources. Our analysis shows that the two proton numbers are poorly correlated, with their ratio spanning more than five orders of magnitude, suggesting that the back precipitation of shock-acceleration particles is unlikely the source of the LDGRF emission.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

We present an interpretation of the time variability of the $X$-ray flux recently reported from a multi-epoch campaign of $15$ years observations of the supernova remnant Cassiopeia A by {\it Chandra}. We show for the first time quantitatively that the $[4.2-6]$ keV non-thermal flux increase up to $50\%$ traces the growth of the magnetic field due to vortical amplification mechanism at a reflection inward shock colliding with inner overdensities. The fast synchrotron cooling as compared with shock-acceleration time scale qualitatively supports the flux decrease.

Recent observations have shown that coronal shocks driven by coronal mass ejections can develop and accelerate particles within several solar radii in large solar energetic particle (SEP) events. Motivated by this, we present an SEP acceleration study including the process that a fast shock propagates through a streamer-like magnetic field with both closed and open field lines in the low corona region. The acceleration of protons is modeled by numerically solving the Parker transport equation with spatial diffusion both along and across the magnetic field. We show that particles can be sufficiently accelerated to up to several hundred MeV within 2-3 solar radii. When the shock propagates through a streamer-like magnetic field, particles are more efficiently accelerated compared to the case with a simple radial magnetic field, mainly due to perpendicular shock geometry and the natural trapping effect of closed magnetic fields. Our results suggest that the coronal magnetic field configuration is an important factor for producing large SEP events. We further show that the coronal magnetic field configuration strongly influences the distribution of energetic particles, leading to different locations of source regions along the shock front where most of high energy particles are concentrated. This work may have strong implications to SEP observations. The upcoming Parker Solar Probe will provide in situ observations for the distribution of energetic particles in the coronal shock region, and test the results of the study.

Late on July 23, 2012, the STEREO-A spacecraft encountered a fast forward shock driven by a coronal mass ejection launched from the Sun earlier that same day. The estimated travel time of the disturbance ($\sim 20$ hrs), together with the massive magnetic field strengths measured within the ejecta ($> 100$nT), made it one of the most extreme events observed during the space era. In this study, we examine the properties of the shock wave. Because of an instrument malfunction, plasma measurements during the interval surrounding the CME were limited, and our approach has been modified to capitalize on the available measurements and suitable proxies, where possible. We were able to infer the following properties. First, the shock normal was pointing predominantly in the radial direction (${\bf n} = 0.97 {\bf e}_r -0.09 {\bf e}_t -0.23 {\bf e}_n$). Second, the angle between ${\bf n}$ and the upstream magnetic field, $\theta_{Bn}$, was estimated to be $\approx 34^{\circ}$, making the shock "quasi-parallel," and supporting the idea of an earlier "preconditioning" ICME. Third, the shock speed was estimated to be $\approx 3300$ km s$^{-1}$. Fourth, the sonic Mach number, $M_s$, for this shock was $\sim 28$. We support these results with an idealized numerical simulation of the ICME. Finally, we estimated the change in ram pressure upstream of the shock to be $\sim 5$ times larger than the pressure from the energetic particles, suggesting that this was not a standard "steady-state" cosmic-ray modified shock (CRMS). Instead it might represent an early, transient phase in the evolution of the CRMS.

Observations show that there is a proton spectral "break" with E$_{break}$ at 1-10MeV in some large CME-driven shocks. Theoretical model usually attribute this phenomenon to a diffusive shock acceleration. However, the underlying physics of the shock acceleration still remains uncertain. Although previous numerical models can hardly predict this "break" due to either high computational expense or shortcomings of current models, the present paper focuses on simulating this energy spectrum in converged two shocks by Monte Carlo numerical method. Considering the Dec 13 2006 CME-driven shock interaction with an Earth bow shock, we examine whether the energy spectral "break" could occur on an interaction between two shocks. As result, we indeed obtain the maximum proton energy up to 10MeV, which is the premise to investigate the existence of the energy spectral "break". Unexpectedly, we further find a proton spectral "break" appears distinctly at the energy $\sim$5MeV.

We investigate the spatial distribution of charged particles accelerated by non-relativistic oblique fast collisionless shocks using three-dimensional test-particle simulations. We find that the density of low-energy particles exhibit a localised enhancement at the shock, resembling the "spike" measured at interplanetary shocks. In contrast to previous results based on numerical solutions to the focused transport equation, we find a shock spike for any magnetic obliquity, from quasi-perpendicular to parallel. We compare the pitch-angle distribution with respect to the local magnetic field and the momentum distribution far downstream and very near the shock within the spike; our findings are compatible with predictions from the scatter-free shock drift acceleration (SDA) limit in these regions. The enhancement of low-energy particles measured by Voyager 1 at solar termination shock is comparable with our profiles. Our simulations allow for predictions of supra-thermal protons at interplanetary shocks within ten solar radii to be tested by Solar Probe Mission. They also have implications for the interpretation of ions accelerated at supernova remnant shocks.

We perform a numerical-simulation study of the acceleration of electrons at shocks that propagate through a prespecified, kinematically defined turbulent magnetic field. The turbulence consists of broadband magnetic fluctuations that are embedded in the plasma and cover a range of wavelengths, the smallest of which is larger than the gyroadii of electrons that are initially injected into the system. We find that when the variance of the turbulent component of the upstream magnetic field is sufficiently large -- $\sigma^2 \sim$ 10 $B_0^2$, where $B_0$ is the strength of the background magnetic field -- electrons can be efficiently accelerated at a collisionless shock regardless of the orientation of the mean upstream magnetic field relative to the shock-normal direction. Since the local angle between the incident magnetic-field vector and the shock-normal vector can be quite large, electrons can be accelerated through shock-drift acceleration at the shock front. In the upstream region, electrons are mirrored back to the shock front leading to multiple shock encounters. Eventually the accelerated electrons are energetic enough that their gyroradii are of the same order as the wavelength of waves that are included in our description of the turbulent magnetic field. Our results are consistent with recent \textit{in situ} observations at Saturn's bow shock. The study may help understand the acceleration of electrons at shocks in space and astrophysical systems.

Using numerical simulations of charged-particles propagating in the heliospheric magnetic field, we study small-scale gradients, or "dropouts", in the intensity of solar energetic particles seen at 1 AU. We use two turbulence models, the foot-point random motion model (Jokipii & Parker 1969; Giacalone et al. 2006) and two-component model (Matthaeus et al. 1990), to generate fluctuating magnetic fields similar to spacecraft observations at 1 AU. The turbulence models include a Kolmogorov-like magnetic field power spectrum containing a broad range of spatial scales from those that lead to large-scale field-line random walk to small scales leading to resonant pitch-angle scattering of energetic particles. We release energetic protons (20 keV - 10 MeV) from a spatially compact and instantaneous source. The trajectories of energetic charged particles in turbulent magnetic fields are numerically integrated. Spacecraft observations are mimicked by collecting particles in small windows when they pass the windows at a distance of 1 AU. We show that small-scale gradients in the intensity of energetic particles and velocity dispersions observed by spacecraft can be reproduced using the foot-point random motion model. However, no dropouts is seen in simulations using the two-component magnetic turbulence model. We also show that particle scattering in the solar wind magnetic field needs to be infrequent for intensity dropouts to form.

We present three-dimensional hybrid simulations of collisionless shocks that propagate parallel to the background magnetic field to study the acceleration of protons that forms a high-energy tail on the distribution. We focus on the initial acceleration of thermal protons and compare it with results from one-dimensional simulations. We find that for both one- and three-dimensional simulations, particles that end up in the high-energy tail of the distribution later in the simulation gained their initial energy right at the shock. This confirms previous results but is the first to demonstrate this using fully three-dimensional fields. The result is not consistent with the "thermal leakage" model. We also show that the gyrocenters of protons in the three-dimensional simulation can drift away from the magnetic field lines on which they started due to the removal of ignorable coordinates that exist in one- and two-dimensional simulations. Our study clarifies the injection problem for diffusive shock acceleration.

We present a study of the acceleration of electrons at a perpendicular shock that propagates through a turbulent magnetic field. The energization process of electrons is investigated by utilizing a combination of hybrid (kinetic ions and fluid electron) simulations and test-particle electron simulations. In this method, the motions of the test-particle electrons are numerically integrated in the time-dependent electric and magnetic fields generated by two-dimensional hybrid simulations. We show that large-scale magnetic fluctuations effect electrons in a number of ways and lead to efficient and rapid energization at the shock front. Since the electrons mainly follow along magnetic lines of force, the large-scale braiding of field lines in space allows the fast-moving electrons to interact with the shock front and get accelerated multiple times. Ripples in the shock front occurring at various scales will also contribute to the acceleration by mirroring the electrons. Our calculation shows that this process favors electron acceleration at perpendicular shocks. The acceleration efficiency is critically dependent on the turbulence amplitude and coherence length. We also discuss the implication of this study for solar energetic particles (SEPs) by comparing the acceleration of electrons with that of protons. Their correlation indicates that perpendicular shocks play an important role in SEP events.

Using test-particle simulations, we investigate the temporal dependence of the two-point velocity correlation function for charged particles scattering in a time-independent spatially fluctuating magnetic field derived from a three-dimensional isotropic turbulence power spectrum. Such a correlation function allowed us to compute the spatial coefficients of diffusion both parallel and perpendicular to the average magnetic field. Our simulations confirm the dependence of the perpendicular diffusion coefficient on turbulence energy density and particle energy predicted previously by a model for early-time charged particle transport. Using the computed diffusion coefficients, we exploit the particle velocity autocorrelation to investigate the time-scale over which the particles "decorrelate" from the solution to the unperturbed equation of motion. Decorrelation time-scales are evaluated for parallel and perpendicular motions, including the drift of the particles from the local magnetic field line. The regimes of strong and weak magnetic turbulence are compared for various values of the ratio of the particle gyroradius to the correlation length of the magnetic turbulence. Our simulation parameters can be applied to energetic particles in the interplanetary space, cosmic rays at the supernova shocks, and cosmic-rays transport in the intergalactic medium.

We examine the effects of streaming cosmic rays upstream of a strong, parallel collisionless shock. We include explicitly the inertia of the cosmic rays in our analysis, which was neglected in previous work. For parameters relevant to the acceleration of cosmic rays at a supernova blast wave, we find no MHD fluid instability that would lead to the amplification of the magnetic field above that given by the compression at the shock. We show how to recover, from our own analysis, the cosmic-ray-driven MHD fluid instability found by previous authors. We conclude that including the inertia and dynamics of the cosmic rays keeps the system stable. More over, the cosmic ray current leads to an additional Hall-like term in the magnetic evolution equation. The implications of this paper for acceleration of galactic cosmic rays at supernova remnants are briefly discussed.

We investigate the acceleration of charged particles (both electrons and protons) at collisionless shocks predicted to exist in the vicinity of solar flares. The existence of standing termination shocks has been examined by flare models and numerical simulations e.g., Shibata,Forbes. We study electron energization by numerically integrating the equations of motion of a large number of test-particle electrons in the time-dependent two-dimensional electric and magnetic fields generated from hybrid simulations (kinetic ions and fluid electron) using parameters typical of the solar flare plasma environment. The shock is produced by injecting plasma flow toward a rigid piston. Large-scale magnetic fluctuations -- known to exist in plasmas and known to have important effects on the nonthermal electron acceleration at shocks -- are also included in our simulations. For the parameters characteristic of the flaring region, our calculations suggest that the termination shock formed in the reconnection outflow region (above post-flare loops) could accelerate electrons to a kinetic energy of a few MeV within 100 ion cyclotron periods, which is of the order of a millisecond. Given a sufficient turbulence amplitude level ($\delta B^2/B_0^2 \sim 0.3$), about 10% of thermal test-particle electrons are accelerated to more than 15 keV. We find that protons are also accelerated, but not to as high energy in the available time and the energy spectra are considerably steeper than that of the electrons for the parameters used in our simulations. Our results are qualitatively consistent with the observed hard X-ray emissions in solar flares.

We have performed extensive two-dimensional magnetohydrodynamic simulations to study the amplification of magnetic fields when a supernova blast wave propagates into a turbulent interstellar plasma. The blast wave is driven by injecting high pressure in the simulation domain. The interstellar magnetic field can be amplified by two different processes, occurring in different regions. One is facilitated by the fluid vorticity generated by the ``rippled" shock front interacting with the background turbulence. The resulting turbulent flow keeps amplifying the magnetic field, consistent with earlier work \citep{Giacalone2007}. The other process is facilitated by the growth of the Rayleigh-Taylor instability at the contact discontinuity between the ejecta and the shocked medium. This can efficiently amplify the magnetic field and tends to produce the highest magnetic field. We investigate the dependence of the amplification on numerical parameters such as grid-cell size and on various physical parameters. We show the magnetic field has a characteristic radial profile that the downstream magnetic field gets progressively stronger away from the shock. This is because the downstream magnetic field needs a finite time to reach the efficient amplification, and will get further amplified in the Rayleigh-Taylor region. In our simulation we do not observe a systematic strong magnetic field within a small distance to the shock. This indicates that if the magnetic-field amplification in supernova remnants indeed occurs near the shock front, other processes such as three-dimensional instabilities, plasma kinetics and/or cosmic ray effect may need to be considered to explain the strong magnetic field in supernova remnants.

We study the physics of electron acceleration at collisionless shocks that move through a plasma containing large-scale magnetic fluctuations. We numerically integrate the trajectories of a large number of electrons, which are treated as test particles moving in the time dependent electric and magnetic fields determined from 2-D hybrid simulations (kinetic ions, fluid electron). The large-scale magnetic fluctuations effect the electrons in a number of ways and lead to efficient and rapid energization at the shock front. Since the electrons mainly follow along magnetic lines of force, the large-scale braiding of field lines in space allows the fast-moving electrons to cross the shock front several times, leading to efficient acceleration. Ripples in the shock front occuring at various scales will also contribute to the acceleration by mirroring the electrons. Our calculation shows that this process favors electron acceleration at perpendicular shocks. The current study is also helpful in understanding the injection problem for electron acceleration by collisionless shocks. It is also shown that the spatial distribution of energetic electrons is similar to in-situ observations (e.g., Bale et al. 1999; Simnett et al. 2005). The process may be important to our understanding of energetic electrons in planetary bow shocks and interplanetary shocks, and explaining herringbone structures seen in some type II solar radio bursts.

We present an alternate interpretation of recent STEREO/STE observations that were originally attributed to energetic neutral atoms (ENA) from the heliosheath. The signal attributed to the diffuse ENA source instead shows the characteristics of a point source. We point out that the peak intensity seen by STEREO/STE is centered at the ecliptic longitude of the bright X-ray source Sco X-1. The observed energy spectrum and intensity are also consistent with the X-rays from Sco X-1. The problem of energy dissipation at the solar wind termination shock remains unsolved while current understanding of the interaction between the solar wind and interstellar wind awaits future observations.