Signatures of Massive Black Hole Merger Host Galaxies from Cosmological Simulations. II. Unique Stellar Kinematics in Integral Field Unit Spectroscopy

Abstract Secure methods for identifying the host galaxies of individual massive black hole (MBH) binaries and mergers detected by gravitational-wave experiments such as the Laser Interferometer Space Antenna and pulsar timing arrays are currently lacking, but will be critical to a variety of science goals. Recently, in Bardati et al., we used the Romulus25 cosmological simulation to show that MBH merger host galaxies have unique morphologies in imaging, due to their stronger bulges. Here, we use the same sample of simulated MBH merger host galaxies to investigate their stellar kinematics, as probed by optical integral field unit (IFU) spectroscopy. We perform stellar population synthesis and dust radiative transfer to generate synthetic 3D optical spectral data cubes of each simulated galaxy, and produce mock stellar kinematic maps. Based on a linear discriminant analysis of a combination of kinematic parameters derived from these maps, we show that this approach can identify MBH binary and merger host galaxies with accuracies that increase with chirp mass and mass ratio. For mergers with high chirp masses (≳10 8.2 M ⊙ ) and high mass ratios (≳0.5), the accuracies reach ≳85%, and their host galaxies are uniquely characterized by slower rotation and stronger stellar kinematic misalignments. These kinematic properties are commonly associated with massive early-type galaxies that have experienced major mergers, and naturally act as signposts for MBH binaries and mergers with high chirp masses and mass ratios. These results suggest that IFU spectroscopy should also play a role in telescope follow-up of future MBH binaries and mergers detected in gravitational waves.