Perpendicular Ion Heating by Reduced Magnetohydrodynamic Turbulence

Recent theoretical studies argue that the rate of stochastic ion heating in low-frequency Alfv\'en-wave turbulence is given by $Q_\perp = c_1 [(\delta u)^3 /\rho] \exp(-c_2/\epsilon)$, where $\delta u$ is the rms turbulent velocity at the scale of the ion gyroradius $\rho$, $\epsilon = \delta u/v_{\perp \rm i}$, $v_{\perp \rm i}$ is the perpendicular ion thermal speed, and $c_1$ and $c_2$ are dimensionless constants. We test this theoretical result by numerically simulating test particles interacting with strong reduced magnetohydrodynamic (RMHD) turbulence. The heating rates in our simulations are well fit by this formula. The best-fit values of $c_1$ are $\sim 1$. The best-fit values of $c_2$ decrease (i.e., stochastic heating becomes more effective) as the grid size and Reynolds number of the RMHD simulations increase. As an example, in a $1024^2 \times 256$ RMHD simulation with a dissipation wavenumber of order the inverse ion gyroradius, we find $c_2 = 0.21$. We show that stochastic heating is significantly stronger in strong RMHD turbulence than in a field of randomly phased Alfv\'en waves with the same power spectrum, because coherent structures in strong RMHD turbulence increase orbit stochasticity in the regions where ions are heated most strongly. We find that $c_1$ increases by a factor of $\sim 3$ while $c_2$ changes very little as the ion thermal speed increases from values $\ll v_{\rm A}$ to values $\sim v_{\rm A}$, where $v_{\rm A}$ is the Alfv\'en speed. We discuss the importance of these results for perpendicular ion heating in the solar wind.