Evolution of stellar orbits around merging massive black hole binary

We study the long-term orbital evolution of stars around a merging massive or supermassive black hole binary (BHB), taking into account the general relativistic effect induced by the black hole (BH) spin. When the BH spin is significant compared to and misaligned with the binary orbital angular momentum, the orbital axis ($\hat{\boldsymbol {l}}$) of the circumbinary star can undergo significant evolution during the binary orbital decay driven by gravitational radiation. Including the spin effect of the primary (more massive) BH, we find that starting from nearly coplanar orbital orientations, the orbital axes $\hat{\boldsymbol {l}}$ of circumbinary stars preferentially evolve towards the spin direction after the merger of the BHB, regardless of the initial BH spin orientation. Such alignment phenomenon, i.e. small final misalignment angle between $\hat{\boldsymbol {l}}$ and the spin axis of the remnant BH $\hat{\boldsymbol {S}}$, can be understood analytically using the principle of adiabatic invariance. For the BHBs with extremely mass ratio (m₂/m₁ ≲ 0.01), $\hat{\boldsymbol {l}}$ may experience more complicated evolution as adiabatic invariance breaks down, but the trend of alignment still works reasonably well when the initial binary spin-orbit angle is relatively small. Our result suggests that the correlation between the orientations of stellar orbits and the spin axis of the central BH could provide a potential signature of the merger history of the massive BH.