Comparison of Model ENAs Produced from Heliospheric Multi-fluid MHD with the First All-Sky ENA Maps

Using a multi-fluid MHD plus neutrals model of the heliosphere, we produce and compare model energetic neutral atom (ENA) maps to the first observations of the Interstellar Boundary Explorer (IBEX) mission, specifically comparing to data from the IBEX - Hi sensor at energies of 0.5 keV to 6 keV. We explore how model ENA maps are affected by varying several parameters including the Very Local Interstellar Medium (VLISM) magnetic field strength and orientation, the assumed plasma distribution of the inner heliosheath, and the plasma distribution of the outer heliosheath. If a suprathermal population exists in the outer heliosheath, then at energies < 2 keV the high density region of VLISM plasma and neutrals along the nose of the heliopause produces as many, or more, ENAs as the inner heliosheath. We provide a system of metrics for comparison with the IBEX data and discuss the extent to which our model agrees or disagrees with the first complete ENA sky maps. Based on the metric results, we suggest future improvements on our model.