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.

The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) suite of instruments operated throughout the over two years of the Rosetta mission operations in the vicinity of comet 67P/Churyumov-Gerasimenko. It measured gas densities and composition throughout the comet's atmosphere, or coma. Here we present two-years' worth of measurements of the relative densities of the four major volatile species in the coma of the comet, H2O. CO2, CO and O2, by one of the ROSINA sub-systems called the Double Focusing Mass Spectrometer (DFMS). The absolute total gas densities were provided by the Comet Pressure Sensor (COPS), another ROSINA sub-system. DFMS is a very high mass resolution and high sensitivity mass spectrometer able to resolve at a tiny fraction of an atomic mass unit. We have analyzed the combined DFMS and COPS measurements using an inversion scheme based on spherical harmonics that solves for the distribution of potential surface activity of each species as the comet rotates, changing solar illumination, over short intervals and as the comet changes distance from the sun and orientation of its spin axis over long time intervals. We also use the surface boundary conditions derived from the inversion scheme to simulate the whole coma with our fully kinetic Direct Simulation Monte Carlo model and calculate the production rates of the four major species throughout the mission. We compare the derived production rates with revised remote sensing observations by the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) as well as with published observations from the Microwave Instrument for the Rosetta Orbiter (MIRO). Finally we use the variation of the surface production of the major species to calculate the total mass loss over the mission and, for different estimates of the dust/gas ratio, calculate the variation of surface loss over the nucleus.

The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its two-year-long mission from 2014 to 2016, Rosetta remained in close proximity to 67P/Churyumov-Gerasimenko, studying its coma and nucleus in situ. Based on lobe-specific measurements of HDO and H2O performed with the ROSINA DFMS mass spectrometer on board Rosetta, the Deuterium-to-Hydrogen ratios in water from the two lobes could be compared. No appreciable difference was observed, suggesting that both lobes formed in the same region and are homogeneous in their Deuterium-to-Hydrogen ratios.

We study roles of the thermosphere and exosphere on the Martian ionospheric structure and ion escape rates in the process of the solar wind-Mars interaction. We employ a four-species multifluid MHD (MF-MHD) model to simulate the Martian ionosphere and magnetosphere. The $cold$ thermosphere background is taken from the Mars Global Ionosphere Thermosphere Model (M-GITM) and the $hot$ oxygen exosphere is adopted from the Mars exosphere Monte Carlo model - Adaptive Mesh Particle Simulator (AMPS). A total of four cases with the combination of 1D (globally averaged) and 3D thermospheres and exospheres are studied. The ion escape rates calculated by adopting 1D and 3D atmospheres are similar; however, the latter are required to adequately reproduce MAVEN ionospheric observations. In addition, our simulations show that the 3D hot oxygen corona plays an important role in preventing planetary molecular ions (O$_2^+$ and CO$_2^+$) escaping from Mars, mainly resulting from the mass loading of the high-altitude exospheric O$^+$ ions. The $cold$ thermospheric oxygen atom, however, is demonstrated to be the primary neutral source for O$^+$ ion escape during the relatively weak solar cycle 24.

In this Letter, we make use of sophisticated 3D numerical simulations to assess the extent of atmospheric ion and photochemical losses from Mars over time. We demonstrate that the atmospheric ion escape rates were significantly higher (by more than two orders of magnitude) in the past at $\sim 4$ Ga compared to the present-day value owing to the stronger solar wind and higher ultraviolet fluxes from the young Sun. We found that the photochemical loss of atomic hot oxygen dominates over the total ion loss at the current epoch whilst the atmospheric ion loss is likely much more important at ancient times. We briefly discuss the ensuing implications of high atmospheric ion escape rates in the context of ancient Mars, and exoplanets with similar atmospheric compositions around young solar-type stars and M-dwarfs.

Magnetohydrodynamics simulations have been carried out in studying the solar wind and cometary plasma interactions for decades. Various plasma boundaries have been simulated and compared well with observations for comet 1P/Halley. The Rosetta mission, which studies comet 67P/Churyumov-Gerasimenko, challenges our understanding of the solar wind and comet interactions. The Rosetta Plasma Consortium observed regions of very weak magnetic field outside the predicted diamagnetic cavity. In this paper, we simulate the inner coma with the Hall magnetohydrodynamics equations and show that the Hall effect is important in the inner coma environment. The magnetic field topology becomes complex and magnetic reconnection occurs on the dayside when the Hall effect is taken into account. The magnetic reconnection on the dayside can generate weak magnetic filed regions outside the global diamagnetic cavity, which may explain the Rosetta Plasma Consortium observations. We conclude that the substantial change in the inner coma environment is due to the fact that the ion inertial length (or gyro radius) is not much smaller than the size of the diamagnetic cavity.