Authors: Senbei Du (Boston University), Merav Opher (Boston University), Joe Giacalone (University of Arizona), Fan Guo (Los Alamos National Laboratory), John D. Richardson (Massachusetts Institute of Technology), Bertalan Zieger (Boston University)

The distribution of ions in the heliosheath is important for understanding remote observations of energetic neutral atoms (ENAs). The ion distribution has been estimated based on hybrid simulations of the heating and evolution of solar wind and interstellar pickup ions across the solar wind termination shock, but this only provides an estimate of the distribution near the shock. It remains unclear how the distribution evolves in the heliosheath. We use hybrid kinetic simulations to investigate the effects of turbulence on ion distributions in the heliosheath. The simulations are compared against Voyager observations, constraining the feasible amplitude and compressibility of turbulence. We find that the heating due to turbulent dissipation can lead to a significant increase of cold solar wind ion temperature. Both turbulent velocity fluctuations and the heating of solar wind ions increase the charge-exchange source for ENAs at low energies (around 100 eV), where current ENA models underpredict observations by more than an order of magnitude. However, the effects of turbulence are likely not strong enough to fully explain these discrepancies.