Solar wind proton precipitation to the surface of Mercury

We examine the effects of the interplanetary magnetic field (IMF) orientation and solar wind dynamic pressure on the solar wind proton precipitation to the surface of Mercury. We use the Amitis model, a three-dimensional GPU-based hybrid model of plasma (particle ions and fluid electrons), and explain a method we found necessary to accurately calculate plasma precipitation to the surface of Mercury through the highly dynamic Hermean magnetosphere. We use our model to explain ground-based telescope observations of Mercury's neutral sodium exosphere, and compare our simulation results with MESSENGER observations. For the typical solar wind dynamic pressure near the orbit of Mercury (i.e., ~7-8 nPa) our model shows a high solar wind proton flux precipitates through the magnetospheric cusps to the high latitudes on both hemispheres on the dayside with a higher precipitation rate to the southern hemisphere compared to the north, which is associated with the northward displacement of Mercury's intrinsic magnetic dipole. We show that this two peak pattern, which is also a common feature observed for neutral sodium exosphere, is controlled by the radial component (B_x) of the IMF and not the B_z component. Our model also suggests that the southward IMF and its associated magnetic reconnection do not play a major role in controlling plasma precipitation to the surface of Mercury through the magnetospheric cusps, in agreement with MESSENGER observations that show that, unlike the Earth, there is almost no dependence between the IMF angle and magnetic reconnection rate at Mercury. For the typical solar wind dynamic pressure, our model suggests that the solar wind proton precipitation through the cusps is longitudinally centered near noon with ~11^o latitudinal extent in the north and ~21^o latitudinal extent in the south, which is consistent with MESSENGER observations. We found an anti-correlation in the incidence area on the surface and the incidence particle rate between the northern and southern cusp precipitation such that the total area and the total rate through both of the cusps remain constant and independent of the IMF orientation. We also show that the solar wind proton incidence rate to the entire surface of Mercury is higher when the IMF has a northward component and nearly half of the incidence flux impacts the low latitudes on the nightside. During extreme solar events (e.g., Coronal Mass Ejections) a large area on the dayside surface of Mercury is exposed to the solar wind plasma, especially in the southern hemisphere. Our model suggests that over 70 nPa solar wind dynamic pressure is required for the entire surface of Mercury to be exposed to the solar wind plasma.