  We have developed a new method for post-Newtonian, high-precision integrationof stellar systems containing a super-massive black hole (SMBH), splitting theforces on a particle between a dominant central force and perturbations. Weused this method to perform fully collisional N-body simulations ofinspiralling intermediate-mass black holes (IMBHs) in the centre of the MilkyWay. We considered stellar cusps of different power-law indices and analysedthe effects of IMBHs of different masses, all starting from circular orbits atan initial distance of 0.1 pc.  Our simulations show how IMBHs deplete the central cusp of stars, leavingbehind a flatter cusp with slope consistent with what has recently beenobserved. If an additional IMBH spirals into such a flat cusp, it can take 50Myr or longer to merge with the central SMBH, thus allowing for directobservation in the near future. The final merger of the two black holesinvolves gravitational wave radiation which may be observable with plannedgravitational wave detectors.  Furthermore, our simulations reveal detailed properties of the hypervelocitystars (HVSs) created, and how generations of HVSs can be used to trace IMBHs inthe Galactic Centre. We find that significant rotation of HVSs (which would beevidence for an IMBH) can only be expected among very fast stars (v > 1000km/s). Also, the probability of creating a hypervelocity binary star is foundto be very small.
We have developed a new method for post-Newtonian, high-precision integration of stellar systems containing a super-massive black hole (SMBH), splitting the forces on a particle between a dominant central force and perturbations. We used this method to perform fully collisional N-body simulations of inspiralling intermediate-mass black holes (IMBHs) in the centre of the Milky Way.