Figure 9

Figure 9: Evolution of the binary semi-major axis (a) and hardening rate (b) in a set of high accuracy N-body simulations; the initial galaxy model was a low-central-density Plummer sphere [20]. Units are G = Mgal = 1, E = - 1/4, with E the total energy. (a) Dashed lines are simulations with binary mass M1 = M2 = 0.005 and solid lines are for M1 = M2 = 0.02, in units where the total galaxy mass is one. (b) Filled (open) circles are for M = M = 0.005 1 2 (0.02). Crosses indicate the hardening rate predicted by a simple model in which the supply of stars to the binary is limited by the rate at which they can be scattered into the binary’s influence sphere by gravitational encounters. The simulations with largest (M1, M2) exhibit the nearly N - 1 dependence expected in the “empty loss cone” regime that is characteristic of real galaxies.