Traditionally, the solar magnetic field has been considered to have a negligible effect in the outer regions of the heliosphere. Recent works have shown that the solar magnetic field may play a crucial role in collimating the plasma in the heliosheath. Interstellar Boundary Explorer (IBEX) observations of the heliotail indicated a latitudinal structure varying with energy in the energetic neutral atom (ENA) fluxes. At energies ~1 keV, the ENA fluxes show an enhancement at low latitudes and a deficit of ENAs near the poles. At energies >2.7 keV, ENA fluxes had a deficit within low latitudes, and lobes of higher ENA flux near the poles. This ENA structure was initially interpreted to be a result of the latitudinal profile of the solar wind during solar minimum. We extend the work of Kornbleuth et al. (2018) by using solar minimum-like conditions and the recently developed SHIELD model. The SHIELD model couples the magnetohydrodynamic (MHD) plasma solution with a kinetic description of neutral hydrogen. We show that while the latitudinal profile of the solar wind during solar minimum contributes to the lobes in ENA maps, the collimation by the solar magnetic field is important in creating and shaping the two high latitude lobes of enhanced ENA flux observed by IBEX. This is the first work to explore the effect of the changing solar magnetic field strength on ENA maps. Our findings suggest that IBEX is providing the first observational evidence of the collimation of the heliosheath plasma by the solar magnetic field.

Neutral hydrogen has been shown to greatly impact the plasma flow in the heliopshere and the location of the heliospheric boundaries. We present the results of the Solar-wind with Hydrogen Ion Exchange and Large-scale Dynamics (SHIELD) model, a new, self-consistent, kinetic-MHD model of the outer heliosphere within the Space Weather Modeling Framework. The charge-exchange mean free path is on order of the size of the heliosphere; therefore, the neutral atoms cannot be described as a fluid. The SHIELD model couples the MHD solution for a single plasma fluid to the kinetic solution from for neutral hydrogen atoms streaming through the system. The kinetic code is based on the Adaptive Mesh Particle Simulator (AMPS), a Monte Carlo method for solving the Boltzmann equation. The SHIELD model accurately predicts the increased filtration of interstellar neutrals into the heliosphere. In order to verify the correct implementation within the model, we compare the results of the SHIELD model to other, well-established kinetic-MHD models. The SHIELD model matches the neutral hydrogen solution of these studies as well as the shift in all heliospheric boundaries closer to the Sun in comparison the the multi-fluid treatment of the neutral hydrogen atoms. Overall the SHIELD model shows excellent agreement to these models and is a significant improvement to the fluid treatment of interstellar hydrogen.

A recent study by Opher et al. (2015) suggested the heliosphere has a "croissant" shape, where the heliosheath plasma is confined by the toroidal solar magnetic field. The "croissant" heliosphere is in contrast to the classically accepted view of a comet-like tail. We investigate the effect of the "croissant" heliosphere model on energetic neutral atom (ENA) maps. Regardless of the existence of a split tail, the confinement of the heliosheath plasma should appear in ENA maps. ENA maps from the Interstellar Boundary Explorer (IBEX) have shown two high latitude lobes with excess ENA flux at higher energies in the tail of the heliosphere. These lobes could be a signature of the confinement of the heliosheath plasma, while some have argued they are caused by the fast/slow solar wind profile. Here we present ENA maps of the "croissant" heliosphere, focusing on understanding the effect of the heliosheath plasma collimation by the solar magnetic field while using a uniform solar wind. We incorporate pick-up ions (PUIs) into our model based on Malama et al. (2006) and Zank et al. (2010). We use the neutral solution from our MHD model to determine the angular variation of the PUIs, and include the extinction of PUIs in the heliosheath. In the presence of a uniform solar wind, we find that the collimation in the "croissant" heliosphere does manifest itself into two high latitude lobes of increased ENA flux in the downwind direction.