For the first time, we explore the tightly coupled interior-magnetosphere system of Mercury by employing a three-dimensional ten-moment multifluid model. This novel fluid model incorporates the non-ideal effects including the Hall effect, inertia, and tensorial pressures that are critical for collisionless magnetic reconnection; therefore, it is particularly well suited for investigating $collisionless$ magnetic reconnection in Mercury's magnetotail and at the planet's magnetopause. The model is able to reproduce the observed magnetic field vectors, field-aligned currents, and cross-tail current sheet asymmetry (beyond the MHD approach) and the simulation results are in good agreement with spacecraft observations. We also study the magnetospheric response of Mercury to a hypothetical extreme event with an enhanced solar wind dynamic pressure, which demonstrates the significance of induction effects resulting from the electromagnetically-coupled interior. More interestingly, plasmoids (or flux ropes) are formed in Mercury's magnetotail during the event, indicating the highly dynamic nature of Mercury's magnetosphere.

We consider ion heating by turbulent Alfven waves (AWs) and kinetic Alfven waves (KAWs) with perpendicular wavelengths comparable to the ion gyroradius and frequencies smaller than the ion cyclotron frequency. When the turbulence amplitude exceeds a certain threshold, an ion's orbit becomes chaotic. The ion then interacts stochastically with the time-varying electrostatic potential, and the ion's energy undergoes a random walk. Using phenomenological arguments, we derive an analytic expression for the rates at which different ion species are heated, which we test by simulating test particles interacting with a spectrum of randomly phased AWs and KAWs. We find that the stochastic heating rate depends sensitively on the quantity epsilon = dv/vperp, where vperp is the component of the ion velocity perpendicular to the background magnetic field B0, and dv (dB) is the rms amplitude of the velocity (magnetic-field) fluctuations at the gyroradius scale. In the case of thermal protons, when epsilon << eps1, where eps1 is a constant, a proton's magnetic moment is nearly conserved and stochastic heating is extremely weak. However, when epsilon > eps1, the proton heating rate exceeds the cascade power that would be present in strong balanced KAW turbulence with the same value of dv, and magnetic-moment conservation is violated. For the random-phase waves in our test-particle simulations, eps1 is approximately 0.2. For protons in low-beta plasmas, epsilon is approximately dB/B0 divided by the square root of beta, and epsilon can exceed eps1 even when dB/B0 << eps1. At comparable temperatures, alpha particles and minor ions have larger values of epsilon than protons and are heated more efficiently as a result. We discuss the implications of our results for ion heating in coronal holes and the solar wind.

A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.