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.

The interstellar neutral (ISN) gas enters the heliosphere and is detected at a few au from the Sun, as demonstrated by Ulysses and the nterstellar Boundary Explorer (IBEX) missions. Ulysses observed ISN gas from different vantage points in a polar orbit from 1994 to 2007, while IBEX has been observing in an Earth orbit in a fixed direction relative to the Sun from 2009. McComas et al. 2018 reported about an IMAP-Lo detector on board the Interstellar Mapping and Acceleration Probe (IMAP), with an ability to track the ISN flux in the sky. We present observation geometries for ISN gas for a detector with the capability to adjust the boresight direction along the Earth orbit over a year within a multichoice ISN observation scheme. We study science opportunities from the observations as a function of time during a year and the phase of solar activity. We identify observation geometries and determine the observation seasons separately for various ISN species and populations. We find that using an adjustable viewing direction allows for ISN gas observations in the upwind hemisphere, where the signal is not distorted by gravitational focusing, in addition to the viewing of ISN species throughout the entire year. Moreover, we demonstrate that with appropriately adjusted observation geometries, primary and secondary populations can be fully separated. Additionally, we show that atoms of ISN gas on indirect trajectories are accessible for detection, and we present their impact on the study of the ionization rates for ISN species.

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.

The Parker Solar Probe (PSP) spacecraft has flown into the most dense and previously unexplored region of our solar system's zodiacal cloud. While PSP does not have a dedicated dust detector, multiple instruments onboard are sensitive to the effects of meteoroid bombardment. Here, we discuss measurements taken during PSP's first two orbits and compare them to models of the zodiacal cloud's dust distribution. Comparing the radial impact rate trends and the timing and location of a dust impact to an energetic particle detector, we find the impactor population to be consistent with dust grains on hyperbolic orbits escaping the solar system. Assuming PSP's impact environment is dominated by hyperbolic impactors, the total quantity of dust ejected from our solar system is estimated to be 1-14 tons/s. We expect PSP will encounter an increasingly more intense impactor environment as its perihelion distance and semi-major axis are decreased.

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.

Energetic neutral atoms (ENA) are an important tool for investigating the structure of the heliosphere. Recently, it was observed that fluxes of ENAs (with energy $\le$ 55 keV) coming from the upwind and downwind regions of the heliosphere are similar in strength. This led the authors of these observations to hypothesize that the heliosphere is bubble-like rather than comet-like, meaning that it has no extended tail. We investigate the directional distribution of the ENA flux for a wide energy range (3--88 keV) including the observations from IBEX (Interstellar Boundary Explorer), INCA (Ion and Neutral Camera, on board Cassini), and HSTOF (High energy Suprathermal Time Of Flight sensor, on board SOHO, Solar and Heliospheric Observatory). An essential element is the model of pickup ion acceleration at the termination shock (TS) proposed by Zank. We use state of the art models of the global heliosphere, interstellar neutral gas density, and pickup ion distributions. The results, based on the "comet-like" model of the heliosphere, are close in flux magnitude to ENA observations by IBEX, HSTOF and partly by INCA (except for the 5.2-13.5 keV energy channel). We find that the ENA flux from the tail dominates at high energy (in agreement with HSTOF, but not INCA). At low energy, our comet-like model produces the similar strengths of the ENA fluxes from the upwind and downwind directions, which, therefore, removes this as a compelling argument for a bubble-like heliosphere.

In mid-2012, a GMIR observed by Voyager 2 crossed through the heliosheath and collided with the heliopause, generating a pressure pulse that propagated into the very local interstellar medium. The effects of the transmitted wave were seen by Voyager 1 just 93 days after its own heliopause crossing. The passage of the transient was accompanied by long-lasting decreases in galactic cosmic ray intensities that occurred from ~2012.55 to ~2013.35 and ~2012.91 to ~2013.70 at Voyager 2 and Voyager 1, respectively. Omnidirectional (>20 MeV) proton-dominated measurements from each spacecraft's Cosmic Ray Subsystem reveal a remarkable similarity between these causally-related events, with a correlation coefficient of 91.2% and a time-lag of 130 days. Knowing the locations of the two spacecraft, we use the observed time-delay to calculate the GMIR's average speed through the heliosheath (inside the heliopause) as a function of temperature in the very local interstellar medium. This, combined with particle, field, and plasma observations enables us to infer previously unmeasured properties of the heliosheath, including a range of sound speeds and total effective pressures. For a nominal temperature of ~20,000 K just outside the heliopause, we find a sound speed of 314 (+/-) 32 km/s and total effective pressure of 267 (+/-) 55 fPa inside the heliopause. We compare these results with the Interstellar Boundary Explorer's data-driven models of heliosheath pressures derived from energetic neutral atom fluxes (the globally distributed flux) and present them as additional evidence that the heliosheath's dynamics are driven by suprathermal energetic processes.

Interstellar neutral gas atoms penetrate the heliopause and reach 1~au, where they are detected by IBEX. The flow of neutral interstellar helium through the perturbed interstellar plasma in the outer heliosheath (OHS) results in creation of the secondary population of interstellar He atoms, the so-called Warm Breeze, due to charge exchange with perturbed ions. The secondary population brings the imprint of the OHS conditions to the IBEX-Lo instrument. Based on a global simulation of the heliosphere with measurement-based parameters and detailed kinetic simulation of the filtration of He in the OHS, we find the number density of interstellar He$^+$ population at $(8.98\pm 0.12)\times 10^{-3}$~cm$^{-3}$. With this, we obtain the absolute density of interstellar H$^+$ $5.4\times 10^{-2}$~cm$^{-3}$ and electrons $6.3\times 10^{-2}$~cm$^{-3}$, and ionization degrees of H 0.26 and He 0.37. The results agree with estimates of the Very Local Interstellar Matter parameters obtained from fitting the observed spectra of diffuse interstellar EUV and soft X-Ray background.

We describe the initial version of the Solar Particle Event (SPE) Threat Assessment Tool or STAT. STAT relies on elements of Corona-Heliosphere (CORHEL) and the Earth-Moon-Mars Radiation Environment Module (EMMREM), and allows users to investigate coronal mass ejection (CME) driven SPEs using coupled magnetohydrodynamic (MHD) and focused transport solutions. At the present time STAT focuses on modeling solar energetic particle (SEP) acceleration in and transport from the low corona, where the highest energy SEP events are generated. We illustrate STAT's capabilities with a model of the July 14, 2000 "Bastille Day" event, including innovative diagnostics for understanding the three-dimensional distribution of particle fluxes and their relation to the structure of the underlying CME driver. A preliminary comparison with NOAA GOES measurements is shown.

Starlight that becomes linearly polarized by magnetically aligned dust grains provides a viable diagnostic of the interstellar magnetic field (ISMF). A survey is underway to map the local ISMF using data collected at eight observatories in both hemispheres. Two approaches are used to obtain the magnetic structure: statistically evaluating magnetic field directions traced by multiple polarization position angles, and least-squares fits that provide the dipole component of the magnetic field. We find that the magnetic field in the circumheliospheric interstellar medium (CHM), which drives winds of interstellar gas and dust through the heliosphere, drapes over the heliopause and influences polarization measurements. We discover a polarization band that can be described with a great circle that traverses the heliosphere nose and ecliptic poles. A gap in the band appears in a region coinciding both with the highest heliosheath pressure, found by IBEX, and the center of the Loop I superbubble. The least-squares analysis finds a magnetic dipole component of the polarization band with the axis oriented toward the ecliptic poles. The filament of dust around the heliosphere and the warm helium breeze flowing through the heliosphere trace the same magnetic field directions. Regions along the polarization band near the heliosphere nose have magnetic field orientations within 15 degrees of sightlines. Regions in the IBEX ribbon have field directions within 40 degrees of the plane of the sky. Several spatially coherent magnetic filaments are within 15 pc. Most of the low frequency radio emissions detected by the two Voyager spacecraft follow the polarization band. The geometry of the polarization band is compared to the Local Interstellar Cloud, the Cetus Ripple, the BICEP2 low opacity region, Ice Cube IC59 galactic cosmic ray data, and Cassini results.

Direct-sampling observations of interstellar neutral (ISN) He by Interstellar Boundary Explorer (IBEX) provide valuable insight into the physical state of and processes operating in the interstellar medium ahead of the heliosphere. The ISN He atom signals are observed at the four lowest ESA steps of the IBEX-Lo sensor. The observed signal is a mixture of the primary and secondary components of ISN He and H. Previously, only data from one of the ESA steps have been used. Here, we extended the analysis to data collected in the three lowest ESA steps with the strongest ISN He signal, for the observation seasons 2009-2015. The instrument sensitivity is modeled as a linear function of the atom impact speed onto the sensor's conversion surface separately for each ESA step of the instrument. We found that the sensitivity increases from lower to higher ESA steps, but within each of the ESA steps it is a decreasing function of the atom impact speed. This result may be influenced by the hydrogen contribution, which was not included in the adopted model, but seems to exist in the signal. We conclude that the currently accepted temperature of ISN He and velocity of the Sun through the interstellar medium do not need a revision, and we sketch a plan of further data analysis aiming at investigating ISN H and a better understanding of the population of ISN He originating in the outer heliosheath.

We present a comprehensive study of energetic neutral atoms (ENAs) of 10 eV to 2.5 keV from the downwind hemisphere of the heliosphere. These ENAs are believed to originate mostly from pickup protons and solar wind protons in the inner heliosheath. This study includes all low-energy observations made with the Interstellar Boundary Explorer over the first 8 years. Since the protons around 0.1 keV dominate the plasma pressure in the inner heliosheath in downwind direction, these ENA observations offer the unique opportunity to constrain the plasma properties and dimensions of the heliosheath where no in-situ observations are available. We first derive energy spectra of ENA intensities averaged over time for 49 macropixels covering the entire downwind hemisphere. The results confirm previous studies regarding integral intensities and the roll-over around 0.1 keV energy. With the expanded dataset we now find that ENA intensities at 0.2 and 0.1 keV seem to anti-correlate with solar activity. We then derive the product of total plasma pressure and emission thickness of protons in the heliosheath to estimate lower limits on the thickness of the inner heliosheath. The temporally averaged ENA intensities support a rather spherical shape of the termination shock and a heliosheath thickness between 150 and 210 au for most regions of the downwind hemisphere. Around the nominal downwind direction of 76{\deg} ecliptic longitude, the heliosheath is at least 280 au thick. There, the neutral hydrogen density seems to be depleted compared to upwind directions by roughly a factor of 2.

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.

We have recently studied the development of an eruptive filament-driven, large-scale off-limb coronal bright front (OCBF) in the low solar corona (Kozarev et al. 2015), using remote observations from Solar Dynamics Observatory's Advanced Imaging Assembly EUV telescopes. In that study, we obtained high-temporal resolution estimates of the OCBF parameters regulating the efficiency of charged particle acceleration within the theoretical framework of diffusive shock acceleration (DSA). These parameters include the time-dependent front size, speed, and strength, as well as the upstream coronal magnetic field orientations with respect to the front's surface normal direction. Here we present an analytical particle acceleration model, specifically developed to incorporate the coronal shock/compressive front properties described above, derived from remote observations. We verify the model's performance through a grid of idealized case runs using input parameters typical for large-scale coronal shocks, and demonstrate that the results approach the expected DSA steady-state behavior. We then apply the model to the event of May 11, 2011 using the OCBF time-dependent parameters derived in Kozarev et al. (2015). We find that the compressive front likely produced energetic particles as low as 1.3 solar radii in the corona. Comparing the modeled and observed fluences near Earth, we also find that the bulk of the acceleration during this event must have occurred above 1.5 solar radii. With this study we have taken a first step in using direct observations of shocks and compressions in the innermost corona to predict the onsets and intensities of SEP events.

Melott et al. [2016] suggest that individual solar proton events (SPEs) are detectable as nitrate ion spikes in ice cores. They use the high fluence, high energy (hard spectrum) SPE of 23 February 1956 to calculate an enhancement of HNO3 from the surface to 46 km that is equivalent to a ~120 ng cm-2 nitrate ion spike observed in the GISP2H ice core. The Melott et al. [2016] approach is fundamentally flawed, since it considers only the absolute column burden of SPE-produced nitrate and not the pre-existing nitrate in the stratosphere. Modeling studies supported by extensive observations [Duderstadt et al., 2014, 2016, and this comment] show background HNO3 in the lower and middle stratosphere equivalent to 2000 to 3000 ng cm-2 nitrate. These high levels of background nitrate must also be included when estimating SPE enhancements to the deposition of nitrate ions that might eventually be preserved in an ice core. The 1956 SPE results in less than a 5% increase in the column burden of atmospheric HNO3, not large enough to explain the nitrate spike seen in the GISP2H ice core. Even extreme SPE enhancements cannot explain nitrate peaks (typically hundreds of percent increases) observed in the ice record [Duderstadt et al., 2016]. Realistic mechanisms linking nitrate ions in ice cores to SPEs have not been established. It is time to move the search for indicators of SPEs away from nitrate ions: Nitrate ions cannot be used as proxies for individual SPEs in the ice core record.

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.

Maps of Energetic Neutral Atom (ENA) fluxes obtained from Interstellar Boundary Explorer (IBEX) observations revealed a bright structure extending over the sky, subsequently dubbed the IBEX ribbon. The ribbon had not been expected from the existing models and theories prior to IBEX, and a number of mechanisms have since been proposed to explain the observations. In these mechanisms, the observed ENAs emerge from source plasmas located at different distances from the Sun. Since each part of the sky is observed by IBEX twice during the year from opposite sides of the Sun, the apparent position of the ribbon as observed in the sky is shifted due to parallax. To determine the ribbon parallax, we found the precise location of the maximum signal of the ribbon observed in each orbital arc. The obtained apparent positions were subsequently corrected for the Compton-Getting effect, gravitational deflection, and radiation pressure. Finally, we selected a part of the ribbon where its position is similar between the IBEX energy passbands. We compared the apparent positions obtained from the viewing locations on the opposite sides of the Sun, and found that they are shifted by a parallax angle of $0.41^\circ\pm0.15^\circ$, which corresponds to a distance of $140^{+84}_{-38}$ AU. This finding supports models of the ribbon with the source located just outside the heliopause.

With the velocity vector and temperature of the pristine interstellar neutral (ISN) He recently obtained with high precision from a coordinated analysis summarized by McComas et al.2015b, we analyzed the IBEX observations of neutral He left out from this analysis. These observations were collected during the ISN observation seasons 2010---2014 and cover the region in the Earth's orbit where the Warm Breeze persists. We used the same simulation model and a very similar parameter fitting method to that used for the analysis of ISN He. We approximated the parent population of the Warm Breeze in front of the heliosphere with a homogeneous Maxwell-Boltzmann distribution function and found a temperature of $\sim 9\,500$ K, an inflow speed of 11.3 km s$^{-1}$, and an inflow longitude and latitude in the J2000 ecliptic coordinates $251.6^\circ$, $12.0^\circ$. The abundance of the Warm Breeze relative to the interstellar neutral He is 5.7\% and the Mach number is 1.97. The newly found inflow direction of the Warm Breeze, the inflow directions of ISN H and ISN He, and the direction to the center of IBEX Ribbon are almost perfectly co-planar, and this plane coincides within relatively narrow statistical uncertainties with the plane fitted only to the inflow directions of ISN He, ISN H, and the Warm Breeze. This co-planarity lends support to the hypothesis that the Warm Breeze is the secondary population of ISN He and that the center of the Ribbon coincides with the direction of the local interstellar magnetic field. The common plane for the direction of inflow of ISN gas, ISN H, the Warm Breeze, and the local interstellar magnetic field %includes the Sun and is given by the normal direction: ecliptic longitude $349.7^\circ \pm 0.6^\circ$ and latitude $35.7^\circ \pm 0.6^\circ$ in the J2000 coordinates, with the correlation coefficient of 0.85.

Nitrate ion spikes in polar ice cores are contentiously used to estimate the intensity, frequency, and probability of historical solar proton events, quantities that are needed to prepare for potentially society-crippling space weather events. We use the Whole Atmosphere Community Climate Model to calculate how large an event would have to be to produce enough odd nitrogen throughout the atmosphere to be discernible as nitrate peaks at the Earth's surface. These hypothetically large events are compared with probability of occurrence estimates derived from measured events, sunspot records, and cosmogenic radionuclides archives. We conclude that the fluence and spectrum of solar proton events necessary to produce odd nitrogen enhancements equivalent to the spikes of nitrate ions in Greenland ice cores are unlikely to have occurred throughout the Holocene, confirming that nitrate ions in ice cores are not suitable proxies for historical individual solar proton events.

We analyzed observations of interstellar neutral helium (ISN~He) obtained from the Interstellar Boundary Explorer (IBEX) satellite during its first six years of operation. We used a refined version of the ISN~He simulation model, presented in the companion paper by Sokol_et al. 2015, and a sophisticated data correlation and uncertainty system and parameter fitting method, described in the companion paper by Swaczyna et al 2015. We analyzed the entire data set together and the yearly subsets, and found the temperature and velocity vector of ISN~He in front of the heliosphere. As seen in the previous studies, the allowable parameters are highly correlated and form a four-dimensional tube in the parameter space. The inflow longitudes obtained from the yearly data subsets show a spread of ~6 degree, with the other parameters varying accordingly along the parameter tube, and the minimum chi-square value is larger than expected. We found, however, that the Mach number of the ISN~He flow shows very little scatter and is thus very tightly constrained. It is in excellent agreement with the original analysis of ISN~He observations from IBEX and recent reanalyses of observations from Ulysses. We identify a possible inaccuracy in the Warm Breeze parameters as the likely cause of the scatter in the ISN~He parameters obtained from the yearly subsets, and we suppose that another component may exist in the signal, or a process that is not accounted for in the current physical model of ISN~He in front of the heliosphere. From our analysis, the inflow velocity vector, temperature, and Mach number of the flow are equal to lambda_ISNHe = 255.8 +/- 0.5 degree, beta_ISNHe = 5.16 +/- 0.10 degree, T_ISNHe = 7440 +/- 260 K, v_ISNHe = 25.8 +/- 0.4$ km/s, and M_ISNHe = 5.079 +/- 0.028, with uncertainties strongly correlated along the parameter tube.

This paper is one of three companion papers presenting the results of our in-depth analysis of the interstellar neutral helium (ISN He) observations carried out using the IBEX-Lo during the first six Interstellar Boundary Explorer (IBEX) observation seasons. We derive corrections for losses due to the limited throughput of the interface buffer and determine the IBEX spin-axis pointing. We develop an uncertainty system for the data, taking into account the resulting correlations between the data points. This system includes uncertainties due to Poisson statistics, background, spin-axis determination, systematic deviation of the boresight from the prescribed position, correction for the interface buffer losses, and the expected Warm Breeze (WB) signal. Subsequently, we analyze the data from 2009 to examine the role of various components of the uncertainty system. We show that the ISN He flow parameters are in good agreement with the values obtained by the original analysis. We identify the WB as the principal contributor to the global $\chi^2$ values in previous analyses. Other uncertainties have a much milder role and their contributions are comparable to each other. The application of this uncertainty system reduced the minimum $\chi^2$ value 4-fold. The obtained $\chi^2$ value, still exceeding the expected value, suggests that either the uncertainty system may still be incomplete or the adopted physical model lacks a potentially important element, which is likely an imperfect determination of the WB parameters. The derived corrections and uncertainty system are used in the accompanying paper by Bzowski et al. in an analysis of the data from six seasons.

The interstellar magnetic field (ISMF) near the heliosphere is a basic part of the solar neighborhood that can only be studied using polarized starlight. Results of an ongoing survey of polarized starlight are analyzed with the goal of linking the interstellar magnetic field that shapes the heliosphere to the nearby field in interstellar space. New results for the direction of the nearby ISMF, based on a merit function that utilizes polarization position angles, identify several magnetic components. The dominant interstellar field, B_pol, is aligned with the direction L,B= 36.2,49.0 (+/-16.0) degrees and is within 8 degrees of the IBEX Ribbon ISMF direction. Stars tracing B_pol have the same mean distance as stars that do not trace B_pol, but show weaker polarizations consistent with lower column densities of polarizing grains. The variations in the polarization position angle directions indicate a low level of magnetic turbulence. B_pol is found after excluding polarizations that trace a separate magnetic structure that apparently is due to interstellar dust deflected around the heliosphere. Local interstellar cloud velocities relative to the LSR increase with the angles between the LSR velocities and ISMF, indicating that the kinematics of local interstellar material is ordered by the ISMF. Polarization and color excess data are consistent with an extension of Loop I to the solar vicinity. Polarizations are consistent with previous findings of more efficient grain alignment in low column density sightlines. Optical polarization and color excess data indicate the presence of nearby interstellar dust in the BICEP2 field. Color excess E(B-V) indicates an optical extinction of A_V about 0.59 mag in the BICEP2 field, while the polarization data indicate that A_V is larger than 0.09 mag. The IBEX Ribbon ISMF extends to the boundaries of the BICEP2 region.

In this paper, we perform numerical modeling of the interstellar hydrogen fluxes measured by IBEX-Lo during orbit 23 (spring 2009) using a state-of-the-art kinetic model of the interstellar neutral hydrogen distribution in the heliosphere. This model takes into account the temporal and heliolatitudinal variations of the solar parameters as well as non-Maxwellian kinetic properties of the hydrogen distribution due to charge exchange in the heliospheric interface. We found that there is a qualitative difference between the IBEX-Lo data and the modeling results obtained with the three-dimensional, time-dependent model. Namely, the model predicts a larger count rate in energy bin~2 (20-41 eV) than in energy bin~1 (11-21 eV), while the data shows the opposite case. We perform study of the model parameter effects on the IBEX-Lo fluxes and the ratio of fluxes in two energy channels. We shown that the most important parameter, which has a major influence on the ratio of the fluxes in the two energy bins, is the solar radiation pressure. The parameter fitting procedure shows that the best agreement between the model result and the data occurs in the case when the ratio of the solar radiation pressure to the solar gravitation, $\mu_0$, is 1.26$^{+0.06}_{-0.076}$, and the total ionization rate of hydrogen at 1 AU is $\beta_{E,0}=3.7^{+0.39}_{-0.35}\times 10^{-7}$~s$^{-1}$. We have found that the value of $\mu_0$ is much larger than $\mu_0=0.89$, which is the value derived from the integrated solar Lyman-alpha flux data for the period of time studied. We discuss possible reasons for the differences.

A recently discovered filament of polarized starlight that traces a coherent magnetic field is shown to have several properties that are consistent with an origin in the outer heliosheath of the heliosphere: (1) The magnetic field that provides the best fit to the polarization position angles is directed within 6.7+-11 degrees of the observed upwind direction of the flow of interstellar neutral helium gas through the heliosphere. (2) The magnetic field is ordered; the component of the variation of the polarization position angles that can be attributed to magnetic turbulence is small. (3) The axis of the elongated filament can be approximated by a line that defines an angle of 80+/-14 degrees with the plane that is formed by the interstellar magnetic field vector and the vector of the inflowing neutral gas (the "BV" plane). We propose that this polarization feature arises from aligned interstellar dust grains in the outer heliosheath where the interstellar plasma and magnetic field are deflected around the heliosphere. The proposed outer heliosheath location of the polarizing grains requires confirmation by modeling grain-propagation through three-dimensional MHD heliosphere models that simultaneously calculate torques on asymmetric dust grains interacting with the heliosphere.

The journey of the Sun through space carries the solar system through a dynamic interstellar environment that is presently characterized by Mach 1 motion between the heliosphere and the surrounding interstellar medium (ISM). The interaction between the heliosphere and ISM is an evolving process due to the variable solar wind and to interstellar turbulence. Frisch et al. presented a meta-analysis of the historical data on the interstellar wind flowing through the heliosphere and concluded that temporal changes in the ecliptic longitude of the wind were statistically indicated by the data available in the refereed literature at the time of that writing. Lallement and Bertaux disagree with this result, and suggested, for instance, that a key instrumental response function of IBEX-Lo was incorrect and that the STEREO pickup ion data are unsuitable for diagnosing the flow of interstellar neutrals through the heliosphere. Here we show that temporal variations in the interstellar wind through the heliosphere are consistent with our knowledge of local ISM. The statistical analysis of the historical helium wind data is revisited, and a recent correction of a typographical error in the literature is incorporated into the new fits. With this correction, and including no newer IBEX results, these combined data still indicate that a change in the longitude of the interstellar neutral wind over the past forty years is statistically likely, but that a constant flow longitude is now also statistically possible. It is shown that the IBEX instrumental response function is known, and that the STEREO pickup ion data have been correctly utilized in this analysis.

The local interstellar magnetic field affects both the heliosphere and the surrounding cluster of interstellar clouds (CLIC). Measurements of linearly polarized starlight provide the only test of the magnetic field threading the CLIC. Polarization measurements of the CLIC magnetic field show multiple local magnetic structures, one of which is aligned with the magnetic field traced by the center of the "ribbon" of energetic neutral atoms discovered by the Interstellar Boundary Explorer (IBEX). Comparisons between the bulk motion of the CLIC through the local standard of rest, the magnetic field direction, the geometric center of Loop I, and the polarized dust bridge extending from the heliosphere toward the North Polar Spur direction all suggest that the CLIC is part of the rim region of the Loop I superbubble.

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.

We derive kinematic properties for two recent solar coronal transient waves observed off the western solar limb with the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) mission. The two waves occurred over $\sim10$-min intervals on consecutive days - June 12 and 13, 2010. For the first time, off-limb waves are imaged with a high 12-sec cadence, making possible detailed analysis of these transients in the low corona between $\sim1.1$-2.0 solar radii ($R_{s}$). We use observations in the 193 and 211 {\AA} AIA channels to constrain the kinematics of both waves. We obtain initial velocities for the two fronts of $\sim1287$ and $\sim736$ km s$^{-1}$, and accelerations of $-1170$ and $-800$ m s$^{-2}$, respectively. Additionally, differential emission measure analysis shows the June 13 wave is consistent with a weak shock. EUV wave positions are correlated with positions from simultaneous type II radio burst observations. We find good temporal and height association between the two, suggesting that the waves may be the EUV signatures of coronal shocks. Furthermore, the events are associated with significant increases in proton fluxes at 1 AU, possibly related to how waves propagate through the coronal magnetic field. Characterizing these coronal transients will be key to connecting their properties with energetic particle production close to the Sun.

Plasma wave observations from Voyager 1 have recently shown large increases in plasma density, to about 0.1 cm^{-3}, consistent with the density of the local interstellar medium. However, corresponding magnetic field observations continue to show the spiral magnetic field direction observed throughout the inner heliosheath. These apparently contradictory observations may be reconciled if Voyager 1 is inside an interstellar flux transfer event - similar to flux transfer events routinely seen at the Earth's magnetopause. If this were the case, Voyager 1 remains inside the heliopause and based on the Voyager 1 observations we can determine the polarity of the interstellar magnetic field for the first time.

The properties of interstellar clouds near the Sun are ordered by the Loop I superbubble and by the interstellar radiation field. Comparisons of the kinematics and magnetic field of the interstellar gas flowing past the Sun, including the Local Interstellar Cloud (LIC), indicate a geometric relation between Loop I as defined by radio synchrotron emission, and the interstellar magnetic field that polarizes nearby starlight. Depletion of Fe and Mg onto dust grains in the LIC shows a surprising relation to the far ultraviolet interstellar radiation field that is best explained by a scenario for the LIC to be extended, possibly filamentary, porous material drifting through space with the Loop I superbubble. The interstellar velocity and magnetic field measured by the Interstellar Boundary Explorer (IBEX) help anchor our understanding of the physical properties of the nearby interstellar medium.

We report on a numerical investigation of two coronal mass ejections (CMEs) which interact as they propagate in the inner heliosphere. We focus on the effect of the orientation of the CMEs relative to each other by performing four different simulations with the axis of the second CME rotated by 90 degrees from one simulation to the next. Each magneto-hydrodynamic (MHD) simulation is performed in three dimensions (3-D) with the Space Weather Modeling Framework (SWMF) in an idealized setting reminiscent of solar minimum conditions. We extract synthetic satellite measurements during and after the interaction and compare the different cases. We also analyze the kinematics of the two CMEs, including the evolution of their widths and aspect ratios. We find that the first CME contracts radially as a result of the interaction in all cases, but the amount of subsequent radial expansion depends on the relative orientation of the two CMEs. Reconnection between the two ejecta and between the ejecta and the interplanetary magnetic field (IMF) determines the type of structure resulting from the interaction. When a CME with a high inclination with respect to the ecliptic overtakes one with a low inclination, it is possible to create a compound event with a smooth rotation in the magnetic field vector over more than 180 degrees. Due to reconnection, the second CME only appears as an extended "tail", and the event may be mistaken for a glancing encounter with an isolated CME. This configuration differs significantly from the one usually studied of a multiple-magnetic cloud event, which we found to be associated with the interaction of two CMEs with the same orientation.

A key indicator of the galactic environment of the Sun is provided by the magnetic field in the interstellar medium (ISM), which influences the shape of the heliosphere. The direction of the nearby interstellar magnetic field (ISMF) is determined from starlight polarized in the ISM. The local ISMF direction is found from the ISMF direction that provides the best fit to the polarization position angles of nearby stars, using weighted fits to the data. New polarization observations are included in the analysis. The best-fitting ISMF is close to the magnetic field direction traced by the center of the Ribbon of energetic neutral atoms, discovered by the Interstellar Boundary Explorer spacecraft. Both the magnetic field and kinematics of the local ISM are consistent with a scenario where the local ISM is a fragment of the Loop I superbubble. An ordered component of the local ISMF is found in a region where PlanetPol data show that polarization increases with distance. It extends to within 8 parsecs of the Sun and implies a weak curvature in the nearby ISMF. Variations from the ordered component indicate turbulence of +/-23 deg. The local ISMF is generally uniform in direction over spatial scales of 8-200 parsecs so that it appears similar to interarm magnetic fields. The best-fitting ISMF direction also agrees with the position of tail-in spatial asymmetries in GeV-TeV galactic cosmic rays. The peculiar geometrical relation between the CMB dipole moment, the heliosphere nose, and local ISMF is supported by these new results. Radiative torques are not likely to play a role in grain alignment for these polarizations.

Post-launch boresight of the IBEX-Lo instrument onboard the Interstellar Boundary Explorer (IBEX) is determined based on IBEX-Lo Star Sensor observations. Accurate information on the boresight of the neutral gas camera is essential for precise determination of interstellar gas flow parameters. Utilizing spin-phase information from the spacecraft attitude control system (ACS), positions of stars observed by the Star Sensor during two years of IBEX measurements were analyzed and compared with positions obtained from a star catalog. No statistically significant differences were observed beyond those expected from the pre-launch uncertainty in the Star Sensor mounting. Based on the star observations and their positions in the spacecraft reference system, pointing of the IBEX satellite spin axis was determined and compared with the pointing obtained from the ACS. Again, no statistically significant deviations were observed. We conclude that no systematic correction for boresight geometry is needed in the analysis of IBEX-Lo observations to determine neutral interstellar gas flow properties. A stack-up of uncertainties in attitude knowledge shows that the instantaneous IBEX-Lo pointing is determined to within $\sim 0.1\degr$ in both spin angle and elevation using either the Star Sensor or the ACS. Further, the Star Sensor can be used to independently determine the spacecraft spin axis. Thus, Star Sensor data can be used reliably to correct the spin phase when the Star Tracker (used by the ACS) is disabled by bright objects in its field-of-view. The Star Sensor can also determine the spin axis during most orbits and thus provides redundancy for the Star Tracker.

Hydrogen gas is the dominant component of the local interstellar medium. However, due to ionization and interaction with the heliosphere, direct sampling of neutral hydrogen in the inner heliosphere is more difficult than sampling the local interstellar neutral helium, which penetrates deep into the heliosphere. In this paper we report on the first detailed analysis of the direct sampling of neutral hydrogen from the local interstellar medium. We confirm that the arrival direction of hydrogen is offset from that of the local Helium component. We further report the discovery of a variation of the penetrating Hydrogen over the first two years of IBEX observations. Observations are consistent with hydrogen experiencing an effective ratio of outward solar radiation pressure to inward gravitational force greater than unity ({\mu}>1); the temporal change observed in the local interstellar hydrogen flux can be explained with solar variability.

Neutral Interstellar Helium (NISHe) is almost unaffected at the heliospheric interface with the interstellar medium and freely enters the solar system. It provides some of the best information on the characteristics of the interstellar gas in the Local Interstellar Cloud. The Interstellar Boundary Explorer (IBEX) is the second mission to directly detect NISHe. We present a comparison between recent IBEX NISHe observations and simulations carried out using a well-tested quantitative simulation code. Simulation and observation results compare well for times when measured fluxes are dominated by NISHe (and contributions from other species are small). Differences between simulations and observations indicate a previously undetected secondary population of neutral helium, likely produced by interaction of interstellar helium with plasma in the outer heliosheath. Interstellar neutral parameters are statistically different from previous in situ results obtained mostly from the GAS/Ulysses experiment, but they do agree with the local interstellar flow vector obtained from studies of interstellar absorption: the newly-established flow direction is ecliptic longitude 79.2 deg, latitude -5.1 deg, the velocity is \sim 22.8 km/s, and the temperature is 6200 K. These new results imply a markedly lower absolute velocity of the gas and thus significantly lower dynamic pressure on the boundaries of the heliosphere and different orientation of the Hydrogen Deflection Plane compared to prior results from Ulysses. A different orientation of this plane also suggests a new geometry of the interstellar magnetic field and the lower dynamic pressure calls for a compensation by other components of the pressure balance, most likely a higher density of interstellar plasma and strength of interstellar magnetic field.

The Interstellar Boundary Explorer (IBEX) spacecraft is providing the first all-sky maps of the energetic neutral atoms (ENAs) produced by charge-exchange between interstellar neutral \HI\ atoms and heliospheric solar wind and pickup ions in the heliosphere boundary regions. The 'edge' of the interstellar cloud presently surrounding the heliosphere extends less than 0.1 pc in the upwind direction, terminating at an unknown distance, indicating that the outer boundary conditions of the heliosphere could change during the lifetime of the IBEX satellite. Using reasonable values for future outer heliosphere boundary conditions, ENA fluxes are predicted for one possible source of ENAs coming from outside of the heliopause. The ENA production simulations use three-dimensional MHD plasma models of the heliosphere that include a kinetic description of neutrals and a Lorentzian distribution for ions. Based on this ENA production model, it is then shown that the sensitivities of the IBEX 1.1 keV skymaps are sufficient to detect the variations in ENA fluxes that are expected to accompany the solar transition into the next upwind cloud. Approximately 20% of the IBEX 1.1 keV pixels appear capable of detecting the predicted model differences at the $ 3 \sigma$ level, with these pixels concentrated in the Ribbon region. Regardless of the detailed ENA production model, the success of the modeled \BdotR\ directions in reproducing the Ribbon locus, together with our results, indicate that the Ribbon phenomenon traces the variations in the heliosphere distortion caused by the relative pressures of the interstellar magnetic and gaseous components.

Variations in the spatial configuration of the interstellar magnetic field (ISMF) near the Sun can be constrained by comparing the ISMF direction at the heliosphere found from the Interstellar Boundary Explorer spacecraft (IBEX) observations of a 'Ribbon' of energetic neutral atoms (ENAs), with the ISMF direction derived from optical polarization data for stars within ~40 pc. Using interstellar polarization observations towards ~30 nearby stars within 90 deg of the heliosphere nose, we find that the best fits to the polarization position angles are obtained for a magnetic pole directed towards ecliptic coordinates of lambda, beta 263 deg, 37 deg (or galactic coordinates of L,B 38 deg, 23deg), with uncertainties of +/- 35 deg, based on the broad minimum of the best fits and the range of data quality. This magnetic pole is 33 deg from the magnetic pole that is defined by the center of the arc of the ENA Ribbon. The IBEX ENA ribbon is seen in sightlines that are perpendicular to the ISMF as it drapes over the heliosphere. The similarity of the polarization and Ribbon directions for the local ISMF suggest that the local field is coherent over scale sizes of tens of parsecs. The ISMF vector direction is nearly perpendicular to the flow of local ISM through the local standard of rest, supporting a possible local ISM origin related to an evolved expanding magnetized shell. The local ISMF direction is found to have a curious geometry with respect to the cosmic microwave background dipole moment.

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the May 12, 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance ($\approx 1R_\odot$) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

The brightest and most surprising feature in the first all-sky maps of Energetic Neutral Atoms (ENA) emissions (0.2-6 keV) produced by the Interstellar Boundary Explorer (IBEX) is an almost circular ribbon of a ~140{\deg} opening angle, centered at (l,b) = (33{\deg}, 55{\deg}), covering the part of the celestial sphere with the lowest column densities of the Local Interstellar Cloud (LIC). We propose a novel interpretation of the IBEX results based on the idea of ENA produced by charge-exchange between the neutral H atoms at the nearby edge of the LIC and the hot protons of the Local Bubble (LB). These ENAs can reach the Sun's vicinity because of very low column density of the intervening LIC material. We show that a plane-parallel or slightly curved interface layer of contact between the LIC H atoms (n_H = 0.2 cm^-3, T = 6000-7000 K) and the LB protons (n_p = 0.005 cm^-3, T ~ 10^6 K), together with indirect contribution coming from multiply-scattered ENAs from the LB, may be able to explain both the shape of the ribbon and the observed intensities provided that the edge is < (500-2000) AU away, the LIC proton density is (correspondingly) < (0.04-0.01) cm^-3, and the LB contains ~1% of non-thermal protons over the IBEX energy range. If this model is correct, then IBEX, for the first time, has imaged in ENAs a celestial object from beyond the confines of the heliosphere and can directly diagnose the plasma conditions in the LB.

We present new kinetic-gasdynamic model of the solar wind interaction with the local interstellar medium. The model incorporates several processes suggested by McComas et al. (2009) for the origin of the heliospheric ENA ribbon -- the most prominent feature seen in the all sky maps of heliospheric ENAs discovered by the Interstellar Boundary Explorer (IBEX). The ribbon is a region of enhanced fluxes of ENAs crossing almost the entire sky. Soon after the ribbon's discovery it was realized (McComas et al., 2009) that the enhancement of the fluxes could be in the directions where the radial component of the interstellar magnetic field around the heliopause is close to zero (Schwadron et al., 2009). Our model includes secondary charge exchange of the interstellar H atoms with the interstellar pickup protons outside the heliopause and is a further advancement of the kinetic-gasdynamic model by Malama et al. (2006) where pickup protons were treated as a separate kinetic component. Izmodenov et al. (2009) have shown in the frame of Malama's model that the interstellar pickup protons outside the heliopause maybe a significant source of ENAs at energies above 1 keV. The difference between the current work and that of Izmodenov et al. (2009) is in the assumption of no-scattering for newly created pickup protons outside the heliopause. In this limit the model produces a feature qualitatively similar to the ribbon observed by IBEX.

The Solar Wind Around Pluto (SWAP) instrument on New Horizons will measure the interaction between the solar wind and ions created by atmospheric loss from Pluto. These measurements provide a characterization of the total loss rate and allow us to examine the complex plasma interactions at Pluto for the first time. Constrained to fit within minimal resources, SWAP is optimized to make plasma-ion measurements at all rotation angles as the New Horizons spacecraft scans to image Pluto and Charon during the flyby. In order to meet these unique requirements, we combined a cylindrically symmetric retarding potential analyzer (RPA) with small deflectors, a top-hat analyzer, and a redundant/coincidence detection scheme. This configuration allows for highly sensitive measurements and a controllable energy passband at all scan angles of the spacecraft.

Motivated by recent experimental proposals to search for extragalactic cosmic rays (including anti-matter from distant galaxies), we study particle propagation through the intergalactic medium (IGM). We first use estimates of the magnetic field strength between galaxies to constrain the mean free path for diffusion of particles through the IGM. We then develop a simple analytic model to describe the diffusion of cosmic rays. Given the current age of galaxies, our results indicate that, in reasonable models, a completely negligible number of particles can enter our Galaxy from distances greater than $\sim 100$ Mpc for relatively low energies ($E$ $< 10^6$ GeV/n). We also find that particle destruction in galaxies along the diffusion path produces an exponential suppression of the possible flux of extragalactic cosmic rays. Finally, we use gamma ray constraints to argue that the distance to any hypothetical domains of anti-matter must be roughly comparable to the horizon scale.