The Bar--Halo Interaction--I. From Fundamental Dynamics to Revised   N-body Requirements
  Only through resonances can non-axisymmetric features such as spiral arms and bars exert torques over large scales and change the overall structure of a near-equilibrium galaxy. We describe the resonant interaction mechanism in detail and derive explicit criteria for the particle number required to simulate these dynamical processes accurately using N-body simulations and illustrate them with numerical experiments. To do this, we perform direct numerical solution of perturbation theory and make detailed comparisons with N-body simulations. The criteria include: sufficient particle coverage in phase space near the resonance and enough particles to minimize gravitational potential fluctuations that will change the dynamics of the resonant encounter. Some of our more surprising findings are as follows. First, the Inner-Lindblad-like resonance (ILR), responsible for coupling the bar to the central halo cusp, requires almost 10^9 equal mass particles within the virial radius for a Milky-Way-like bar in an NFW profile. Second, orbits that linger near the resonance receive more angular momentum than orbits that move through the resonance quickly. Small-scale fluctuations present in state-of-the-art particle-particle simulations can knock orbits out of resonance, preventing them from lingering and, thereby, decrease the torque. The required particle numbers are sufficiently high for scenarios of interest that apparent convergence in particle number is misleading: the convergence is in the noise-dominated regime. State-of-the-art simulations are not adequate to follow all aspects of secular evolution driven by the bar-halo interaction. We present a procedure to test the requirements for individual N-body codes for the actual problem of interest. [abridged] 
Martin D. Weinberg and Neal Katz