Binary supermassive black holes (SMBH) are expected to form naturally during galaxy mergers. After the dynamical friction phase, when the two SMBHs become gravitationally bound to each other, and a brief stage of initial rapid hardening, the orbit gradually continues to shrink due to three-body interactions with stars that enter the loss cone of the binary. Using the stellar-dynamical Monte Carlo code Raga, we explore the co-evolution of the binary SMBH and the nuclear star cluster in this slow stage, for various combinations of parameters (geometry of the star cluster, primary/secondary SMBH mass, initial eccentricity, inclusion of stellar captures/tidal disruptions). We compare the rates of stellar captures in galactic nuclei containing a binary to those of galaxies with a single SMBH. At early times, the rates are substantially higher in the case of a binary SMBH, but subsequently they drop to lower levels. Only in triaxial systems both the binary hardening rates and the capture rates remain sufficiently high during the entire evolution. We find that the hardening rate is not influenced by star captures, nor does it depend on eccentricity; however, it is higher when the difference between the black hole masses is greater. We confirm that the eccentricity of the binary tends to grow, which may significantly shorten the coalescence time due to earlier onset of gravitational-wave emission. We also explore the properties of the orbits entering the loss cone, and demonstrate that it remains partially full throughout the evolution in the triaxial case, but significantly depleted in the axisymmetric case. Finally, we study the distribution of ejected hypervelocity stars and the corresponding mass deficits in the central parts of the galaxies hosting a binary, and argue that the missing mass is difficult to quantify observationally.

We explore the rates of tidal disruption events (TDEs) of stars by supermassive black holes (SBHs) in galactic nuclei formed in mergers followed by a formation and coalescence of a binary SBH. Such systems initially have a deficit of stars on low-angular-momentum orbits caused by the slingshot process during the binary SBH stage, which tends to reduce the flux of stars into the SBH compared to the steady-state value. On the other hand, a newly formed galactic nucleus has a non-spherical shape which enhances the mixing of stars in angular momentum and thus the TDE rate. In galaxies with relatively low SBH masses (<10^7 Msun), relaxation times are short enough to wash out the anisotropy in initial conditions, and for more massive SBH the enhancement of flux due to non-sphericity turns out to be more important than its suppression due to initial anisotropy. Therefore, the present-day TDE rates generally exceed conventional steady-state estimates based on a spherical isotropic approximation. We thus conjecture that the lower observationally inferred TDE rates compared to theoretical predictions cannot be attributed to the depletion of low-angular-momentum stars by SBH binaries.

We compute the rates of capture of stars by supermassive black holes, using time-dependent Fokker--Planck equation with initial conditions that have a deficit of stars on low-angular-momentum orbits. One class of initial conditions has a gap in phase space created by a binary black hole, and the other has a globally tangentially-anisotropic velocity distribution. We find that for galactic nuclei that are younger than ~0.1 relaxation times, the flux of stars into the black hole is suppressed with respect to the steady-state value. This effect may substantially reduce the number of observable tidal disruption flares in galaxies with black hole masses M>10^7 Msun.

We analize the quasi-spherical accretion in the presence of axisymmetric vortex turbulence. It is shown that in this case the turbulence changes mainly the effective gravity potential but not the effective pressure.