An outstanding problem in heliospheric physics is understanding the acceleration of solar energetic particles (SEP) in coronal mass ejections (CMEs) and flares. A fundamental question is whether the acceleration occurs in interplanetary space, or near the Sun. Recent work has shown that CME-driven shocks may produce SEPs while still below 5 solar radii. In this work we explore SEP acceleration during the onset of CMEs and shocks even lower in the corona, using realistic suprathermal spectra, for a selection of events. We have calculated quiet-time, pre-event suprathermal particle spectra from 1 AU observations, and scaled them back to the low corona to serve as seed spectra. For each event, AIA observations and the CASHeW framework were used to model the compressive/shock wave kinematics and its interaction with the corona. The proton acceleration was then modeled using an analytic diffusive shock acceleration model as the shock waves propagate between 1.05 and 1.3 solar radii. We demonstrate the capability of low coronal shock-related EUV waves to accelerate protons to multi-MeV energies in a matter of minutes, in the very early stages of the associated solar eruptions. We find that strong proton energization occurs for high values of the density jump, Alfven Mach number, and shock speed. In future work the results of this early-stage shock acceleration will be used to model the continued acceleration higher in the corona.

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.

The Interstellar Boundary Explorer (IBEX) has now operated in space for 7 years and returned nearly continuous observations that have led to scientific discoveries and reshaped our entire understanding of the outer heliosphere and its interaction with the local interstellar medium. Here we extend prior work, adding the 2014-2015 data for the first time, and examine, validate, initially analyze, and provide a complete 7-year set of Energetic Neutral Atom (ENA) observations from ~0.1 to 6 keV. The data, maps, and documentation provided here represent the 10th major release of IBEX data and include improvements to various prior corrections to provide the citable reference for the current version of IBEX data. We are now able to study time variations in the outer heliosphere and interstellar interaction over more than half a solar cycle. We find that the Ribbon has evolved differently than the globally distributed flux (GDF), with a leveling off and partial recovery of ENAs from the GDF, owing to solar wind output flattening and recovery. The Ribbon has now also lost its latitudinal ordering, which reflects the breakdown of solar minimum solar wind conditions and exhibits a greater time delay than for the surrounding GDF. Together, the IBEX observations strongly support a secondary ENA source for the Ribbon, and we suggest that this be adopted as the nominal explanation of the Ribbon going forward.

In large SEP events, ions can be accelerated at CME-driven shocks to very high energies. 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 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 dependences 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.

We fit the $\sim$0.1-500 MeV/nucleon H-Fe spectra in 46 large SEP events surveyed by Desai et al. (2016) with the double power-law Band function to obtain a normalization constant, low- and high-energy parameters $\gamma_a$ and $\gamma_b$; and break energy $E_B$. We also calculate the low-energy power-law spectral slope $\gamma_1$. We find that: 1) $\gamma_a$, $\gamma_1$, and $\gamma_b$ are species-independent within a given SEP event, and the spectra steepen with increasing energy; 2) $E_B$'s are well ordered by Q/M ratio, and decrease systematically with decreasing Q/M, scaling as (Q/M)$^\alpha$ with $\alpha$ varying between $\sim$0.2-3; 3) $\alpha$ is well correlated with Fe/O at $\sim$0.16-0.23 MeV/nucleon and CME speed; 4) In most events: $\alpha<$1.4, the spectra steepen significantly at higher energy with $\gamma_b$-$\gamma_a >$3; and 5) Seven out of 9 extreme SEP events (associated with faster CMEs and GLEs) are Fe-rich, have $\alpha >$1.4, have flatter spectra at low and high energies with $\gamma_b$-$\gamma_a <$3. The species-independence of $\gamma_a$, $\gamma_1$, and $\gamma_b$ and the systematic Q/M dependence of $E_B$ within an event, as well as the range of values for $\alpha$ suggest that the formation of double power-laws in SEP events occurs primarily due to diffusive acceleration at near-Sun CME shocks and not due to scattering in the interplanetary turbulence. In most events, the Q/M-dependence of $E_B$ is consistent with the equal diffusion coefficient condition while the event-to-event variations in $\alpha$ are probably driven by differences in the near-shock wave intensity spectra, which are flatter than the Kolmogorov turbulence spectrum but still weaker compared to that inferred for the extreme events.

The acceleration of protons and electrons to high (sometimes GeV/nucleon) energies by solar phenomena is a key component of space weather. These solar energetic particle (SEP) events can damage spacecraft and communications, as well as present radiation hazards to humans. In-depth particle acceleration simulations have been performed for idealized magnetic fields for diffusive acceleration and particle propagation, and at the same time the quality of MHD simulations of coronal mass ejections (CMEs) has improved significantly. However, to date these two pieces of the same puzzle have remained largely decoupled. Such structures may contain not just a shock but also sizable sheath and pileup compression regions behind it, and may vary considerably with longitude and latitude based on the underlying coronal conditions. In this work, we have coupled results from a detailed global three-dimensional MHD time-dependent CME simulation to a global proton acceleration and transport model, in order to study time-dependent effects of SEP acceleration between 1.8 and 8 solar radii in the 2005 May 13 CME. We find that the source population is accelerated to at least 100 MeV, with distributions enhanced up to six orders of magnitude. Acceleration efficiency varies strongly along field lines probing different regions of the dynamically evolving CME, whose dynamics is influenced by the large-scale coronal magnetic field structure. We observe strong acceleration in sheath regions immediately behind the shock.