Pair Of Supermassive Black Holes Research Articles

Dynamics of supermassive black hole triples in the ROMULUS25 cosmological simulation

For a pair of supermassive black holes (SMBHs) in the remnant of a dual galaxy merger, well-known models exist to describe their dynamical evolution until the final coalescence accompanied by the emission of a low-frequency gravitational wave signal. In this article, we investigate the dynamical evolution of three SMBH triple systems recovered from the ROMULUS25 cosmological simulation to explore common dynamical evolution patterns and assess typical coalescence times. For this purpose, we construct initial conditions from the ROMULUS25 data and perform high-resolution gravitodynamical N-body simulations. We track the orbital evolution from the galactic inspiral to the formation of hard binaries at sub-parsec separation and use the observed hardening rates to project the time of coalescence. In all cases, the two heaviest black holes form an efficiently hardening binary that merges within fractions of the Hubble time. The lightest SMBH either gets ejected, forms a stable hierarchical triple system with the heavier binary, forms a hardening binary with the previously merged binary’s remnant, or remains on a wide galactic orbit. The coalescence times of the lighter black holes are thus significantly longer than for the heavier binary, as they experience lower dynamical friction and stellar hardening rates. We observe the formation of hierarchical triples when the density profile of the galactic nucleus is sufficiently steep.

A close quasar pair in a disk-disk galaxy merger at z = 2.17.

Galaxy mergers produce pairs of supermassive black holes (SMBHs), which may be witnessed as dual quasars if both SMBHs are rapidly accreting. The kiloparsec (kpc)-scale separation represents a physical regime sufficiently close for merger-induced effects to be important1 yet wide enough to be directly resolvable with the facilities currently available. Whereas many kpc-scale, dual active galactic nuclei-the low-luminosity counterparts of quasars-have been observed in low-redshift mergers2, no unambiguous dual quasar is known at cosmic noon (z ≈ 2), the peak of global star formation and quasar activity3,4. Here we report multiwavelength observations of Sloan Digital Sky Survey (SDSS) J0749 + 2255 as a kpc-scale, dual-quasar system hosted by a galaxy merger at cosmic noon (z = 2.17). We discover extended host galaxies associated with the much brighter compact quasar nuclei (separated by 0.46″ or 3.8 kpc) and low-surface-brightness tidal features as evidence for galactic interactions. Unlike its low-redshift and low-luminosity counterparts, SDSS J0749 + 2255 is hosted by massive compact disk-dominated galaxies. The apparent lack of stellar bulges and the fact that SDSS J0749 + 2255 already follows the local SMBH mass-host stellar mass relation, suggest that at least some SMBHs may have formed before their host stellar bulges. While still at kpc-scale separations where the host-galaxy gravitational potential dominates, the two SMBHs may evolve into a gravitationally bound binary system in around 0.22 Gyr.

Differential Interferometric Signatures of Close Binaries of Supermassive Black Holes in Active Galactic Nuclei. II. Merged Broad-line Regions

Pairs of supermassive black holes (SMBHs) at different stages are the natural results of galaxy mergers in the hierarchical framework of galaxy formation and evolution. However, identification of close binaries of SMBHs (CB-SMBHs) with sub-parsec separations in observations are still elusive. Recently, unprecedented spatial resolutions achieved by GRAVITY/GRAVITY+ on board the Very Large Telescope Interferometer through spectroastrometry (SA) provide new opportunities to resolve CB-SMBHs. Differential phase curves of CB-SMBHs with two independent broad-line regions (BLRs) are found to have distinguished characteristic structures from a single BLR. Once the CB-SMBH evolves to the stage where BLRs merge to form a circumbinary BLR, it will hopefully be resolved by the pulsar timing array in the near future as sources of nanohertz gravitational waves. In this work, we use a parameterized model for circumbinary BLRs to calculate line profiles and differential phase curves for SA observations. We show that both profiles and phase curves exhibit asymmetries caused by the Doppler boosting effect of accretion disks around individual black holes, depending on the orbital parameters of the binary and geometries of the BLR. We also generate mock SA data using the model and then recover orbital parameters by fitting the mock data. Degeneracies between parameters contribute greatly to uncertainties of parameters but can be eased through joint analysis of multiple-epoch SA observations and reverberation mappings.

UGC 4211: A Confirmed Dual Active Galactic Nucleus in the Local Universe at 230 pc Nuclear Separation

We present multiwavelength high-spatial resolution (∼0.″1, 70 pc) observations of UGC 4211 at z = 0.03474, a late-stage major galaxy merger at the closest nuclear separation yet found in near-IR imaging (0.″32, ∼230 pc projected separation). Using Hubble Space Telescope/Space Telescope Imaging Spectrograph, Very Large Telescope/MUSE+AO, Keck/OSIRIS+AO spectroscopy, and the Atacama Large Millimeter/submillimeter Array (ALMA) observations, we show that the spatial distribution, optical and near-infrared emission lines, and millimeter continuum emission are all consistent with both nuclei being powered by accreting supermassive black holes (SMBHs). Our data, combined with common black hole mass prescriptions, suggest that both SMBHs have similar masses, ∼ 8.1 (south) and ∼ 8.3 (north), respectively. The projected separation of 230 pc (∼6× the black hole sphere of influence) represents the closest-separation dual active galactic nuclei (AGN) studied to date with multiwavelength resolved spectroscopy and shows the potential of nuclear (<50 pc) continuum observations with ALMA to discover hidden growing SMBH pairs. While the exact occurrence rate of close-separation dual AGN is not yet known, it may be surprisingly high, given that UGC 4211 was found within a small, volume-limited sample of nearby hard X-ray detected AGN. Observations of dual SMBH binaries in the subkiloparsec regime at the final stages of dynamical friction provide important constraints for future gravitational wave observatories.

Unveiling the population of dual and lensed active galactic nuclei at sub-arcsec separations

All cosmological models of structure formation predict the existence of a widespread population of dual supermassive black holes in-spiralling inside their common host galaxy, eventually merging and giving rise to intense gravitational waves. These systems can be identified as dual AGNs at kilo parsec separations, but only very few have been confirmed at z>0.5. The appearance of multiple AGNs at small angular separations can also be due to gravitational lensing of single AGNs, which are themselves very important systems for many astrophysical topics. Here we present a novel technique, dubbed the Gaia Multipeak (GMP) method, to obtain large and reliable samples of dual/lensed AGN candidates with sub-arcsec separations by looking for AGNs showing multiple peaks in the light profiles observed by the Gaia satellite. All the GMP-selected sources with high resolution images (26 from the HST archive and 5 from dedicated adaptive-optics assisted imaging at the Large Binocular Telescope) show multiple components with sub-arcsec separation pointing toward a very high reliability of the method. By sampling separations down to ~2 kpc at z>1, this method allows us to probe the physical processes that drive the inspiralling of a pair of SMBHs inside a single galaxy.

First direct dynamical detection of a dual supermassive black hole system at sub-kiloparsec separation

We investigated whether the two recently discovered nuclei in NGC 7727 both host a super-massive black hole (SMBH). We used the high spatial resolution mode of the integral-field spectrograph MUSE on the VLT in adaptive optics mode to resolve the stellar kinematics within the sphere of influence of both putative black holes. We combined the kinematic data with an HST-based mass model and used Jeans models to measure their SMBH mass. We report the discovery of a dual SMBH system in NGC 7727. We detect a SMBH in the photometric center of the galaxy in Nucleus 1, with a mass of MSMBH = 1.54−0.15+0.18 × 108 M⊙. In the second nucleus, which is 500 pc offset from the main nucleus, we also find a clear signal for a SMBH with a mass of MBH = 6.33−1.40+3.32 × 106 M⊙. Both SMBHs are detected at high significance. The off-axis nature of Nucleus 2 makes modeling the system challenging; however, a number of robustness tests suggest that a black hole is required to explain the observed kinematics. The SMBH in the offset Nucleus 2 makes up 3.0% of its total mass, which means its SMBH is over-massive compared to the MBH − MBulge scaling relation. This confirms it as the surviving nuclear star cluster of a galaxy that has merged with NGC 7727. This discovery is the first dynamically confirmed dual SMBH system with a projected separation of less than a kiloparsec and the nearest dynamically confirmed dual SMBH at a distance of 27.4 Mpc. The second Nucleus is in an advanced state of inspiral, and it will eventually result in a 1:24 mass ratio SMBH merger. Optical emission lines suggest Nucleus 2 is a Seyfert galaxy, making it a low-luminosity Active Galactic Nuclei. There are likely many more quiescent SMBHs as well as dual SMBH pairs in the local Universe that have been missed by surveys that focus on bright accretion signatures.

A Catalog of 204 Offset and Dual Active Galactic Nuclei (AGNs): Increased AGN Activation in Major Mergers and Separations under 4 kpc

Abstract During galaxy mergers, gas and dust are driven toward the centers of merging galaxies, triggering enhanced star formation and supermassive black hole (SMBH) growth. Theory predicts that this heightened activity peaks at SMBH separations &lt;20 kpc; if sufficient material accretes onto one or both of the SMBHs for them to become observable as active galactic nuclei (AGNs) during this phase, they are known as offset and dual AGNs, respectively. To better study these systems, we have built the ACS-AGN Merger Catalog, a large catalog (N = 204) of uniformly selected offset and dual AGN observed by the Hubble Space Telescope at 0.2 &lt; z &lt; 2.5 with separations &lt;20 kpc. Using this catalog, we answer many questions regarding SMBH−galaxy coevolution during mergers. First, we confirm predictions that the AGN fraction peaks at SMBH pair separations &lt;10 kpc; specifically, we find that the fraction increases significantly at pair separations of &lt;4 kpc. Second, we find that AGNs in mergers are preferentially found in major mergers and that the fraction of AGNs found in mergers follows a logarithmic relation, decreasing as merger mass ratio increases. Third, we do not find that mergers (nor the major or minor merger subpopulations) trigger the most luminous AGNs. Finally, we find that nuclear column density, AGN luminosity, and host galaxy star formation rate have no dependence on SMBH pair separation or merger mass ratio in these systems, nor do the distributions of these values differ significantly from that of the overall AGN population.

Discovery of a Close-separation Binary Quasar at the Heart of a z ∼ 0.2 Merging Galaxy and Its Implications for Low-frequency Gravitational Waves

Abstract Supermassive black hole (SMBH) binaries with masses of ∼108–109 M ⊙ are expected to dominate the contribution to the as-yet undetected gravitational wave background (GWB) signal at the nanohertz frequencies accessible to pulsar timing arrays. We currently lack firm empirical constraints on the amplitude of the GWB due to the dearth of confirmed SMBH binaries in the required mass range. Using Hubble Space Telescope/Wide Field Camera 3 images, we have discovered a z ∼ 0.2 quasar hosted in a merger remnant with two closely separated (0.″13 or ∼430 pc) continuum cores at the heart of the galaxy SDSS J1010+1413. The two cores are spatially coincident with two powerful [O iii]-emitting point sources with quasar-like luminosities (L AGN ∼ 5 × 1046 erg s−1), suggesting the presence of a bound SMBH system, each with M BH &gt; 4 × 108 M ⊙. We place an upper limit on the merging timescale of the SMBH pair of 2.5 billion years, roughly the universe lookback time at z ∼ 0.2. There is likely a population of quasar binaries similar to SDSS J1010+1413 that contribute to a stochastic GWB that should be detected in the next several years. If the GWB is not detected this could indicate that SMBHs merge only over extremely long timescales, remaining as close separation binaries for many Hubble times, the so-called “final-parsec problem.”

Jet-powered Outflows in Supermassive Black Hole Binary Candidate SDSS J1048+0055

Abstract The search and study of close pairs of supermassive black holes (SMBHs) is important in the study of galaxy mergers that can possibly trigger active galactic nucleus (AGN) activity, and in the context of their evolution into the gravitational wave emitting regime. The quasar SDSS J1048+0055 was identified as an SMBH binary (SMBHB) candidate based on the observed double-peaked [O iii] λλ4959, 5007 emission lines and two distinct radio components separated by ∼20 pc (Zhou et al. 2004). To ascertain the binary nature of this source, we analyzed multifrequency, multiepoch very long baseline interferometry (VLBI) data to investigate its parsec-scale radio properties. The source shows double components with the western feature being brighter than the eastern one. This brighter component has a brightness temperature of ∼1010 K, a spectral index of α = −0.09 ± 0.09 (flat), and is indicative of mildly relativistic beaming. In contrast, the faint component has a lower brightness temperature of ∼108−9 K and a steep spectrum. These clues are consistent with a core–jet structure; moreover, the apparent separation speed between the two components is much higher than the expected orbital motion in a binary SMBH. Thus the present study excludes the association of the two VLBI components with the cores of an SMBHB; though, the SMBHB possibility (e.g., a pair of radio-loud and radio-quiet AGNs) is not fully ruled out. In the single AGN scenario, the double-peaked optical emission lines can originate from the jet interacting with the narrow-line region as indicated by a change in the jet direction at ∼140 pc.

Swift data hint at a binary supermassive black hole candidate at sub-parsec separation

Dual/binary Supermassive Black Hole (SMBH) systems are the inevitable consequence of the current Lambda Cold Dark Matter cosmological paradigm. In this context, we discuss here the properties of MCG+11-11-032, a local (z=0.0362) Seyfert 2 galaxy. This source was proposed as a dual AGN candidate on the basis of the presence of double-peaked [OIII] emission lines in its optical spectrum. MCG+11-11-032 is also an X-ray variable source and was observed several times by the Swift X-ray Telescope (XRT) on time scales from days to years. In this work, we analyze the SDSS-DR13 spectrum and find evidence for double-peaked profiles in all the strongest narrow emission lines. We also study the XRT light curve and unveil the presence of an alternating behavior of the intrinsic 0.3-10 keV flux, while the 123-month Swift BAT light curve supports the presence of almost regular peaks and dips almost every 25 months. In addition, the XRT spectrum suggests for the presence of two narrow emission lines with rest-frame energies of E~6.16 keV and E~6.56 keV. Although by considering only the optical emission lines, different physical mechanisms may be invoked to explain the kinematical properties, the X-ray results are most naturally explained by the presence of a binary SMBH in the center of this source. In particular, we evidence a remarkable agreement between the putative SMBH pair orbital velocity derived from the BAT light curve and the velocity offset derived by the rest-frame Delta_E between the two X-ray line peaks in the XRT spectrum (i.e. Delta_v~0.06c).

Dancing to changa: a self-consistent prediction for close SMBH pair formation time-scales following galaxy mergers

We present the first self-consistent prediction for the distribution of formation timescales for close Supermassive Black Hole (SMBH) pairs following galaxy mergers. Using ROMULUS25, the first large-scale cosmological simulation to accurately track the orbital evolution of SMBHs within their host galaxies down to sub-kpc scales, we predict an average formation rate density of close SMBH pairs of 0.013 cMpc^-3 Gyr^-1. We find that it is relatively rare for galaxy mergers to result in the formation of close SMBH pairs with sub-kpc separation and those that do form are often the result of Gyrs of orbital evolution following the galaxy merger. The likelihood and timescale to form a close SMBH pair depends strongly on the mass ratio of the merging galaxies, as well as the presence of dense stellar cores. Low stellar mass ratio mergers with galaxies that lack a dense stellar core are more likely to become tidally disrupted and deposit their SMBH at large radii without any stellar core to aid in their orbital decay, resulting in a population of long-lived 'wandering' SMBHs. Conversely, SMBHs in galaxies that remain embedded within a stellar core form close pairs in much shorter timescales on average. This timescale is a crucial, though often ignored or very simplified, ingredient to models predicting SMBH mergers rates and the connection between SMBH and star formation activity.

Supermassive black hole pairs in clumpy galaxies at high redshift: delayed binary formation and concurrent mass growth

Massive gas-rich galaxy discs at $z \sim 1-3$ host massive star-forming clumps with typical baryonic masses in the range $10^7-10^8$~M$_{\odot}$ which can affect the orbital decay and concurrent growth of supermassive black hole (BH) pairs. Using a set of high-resolution simulations of isolated clumpy galaxies hosting a pair of unequal-mass BHs, we study the interaction between massive clumps and a BH pair at kpc scales, during the early phase of the orbital decay. We find that both the interaction with massive clumps and the heating of the cold gas layer of the disc by BH feedback tend to delay significantly the orbital decay of the secondary, which in many cases is ejected and then hovers for a whole Gyr around a separation of 1--2 kpc. In the envelope, dynamical friction is weak and there is no contribution of disc torques: these lead to the fastest decay once the orbit of the secondary BH has circularised in the disc midplane. In runs with larger eccentricities the delay is stronger, although there are some exceptions. We also show that, even in discs with very sporadic transient clump formation, a strong spiral pattern affects the decay time-scale for BHs on eccentric orbits. We conclude that, contrary to previous belief, a gas-rich background is not necessarily conducive to a fast BH decay and binary formation, which prompts more extensive investigations aimed at calibrating event-rate forecasts for ongoing and future gravitational-wave searches, such as with Pulsar Timing Arrays and the future evolved Laser Interferometer Space Antenna.

BINARY BLACK HOLES, GAS SLOSHING, AND COLD FRONTS IN THE X-RAY HALO HOSTING 4C+37.11

ABSTRACT We analyzed deep Chandra ACIS-I exposures of the cluster-scale X-ray halo surrounding the radio source 4C+37.11. This remarkable system hosts the closest resolved pair of super-massive black holes and an exceptionally luminous elliptical galaxy, the likely product of a series of past mergers. We characterize the halo with r 500 ∼ 0.95 Mpc, M 500 = 2.5 ± 0.2 × 1014 M ⊙, kT = 4.6 ± 0.2 keV, and a gas mass of M g,500 = 2.2 ± 0.1 × 1013 M ⊙. The gas mass fraction within r 500 is f g = 0.09 ± 0.01. The entropy profile shows large non-gravitational heating in the central regions. We see several surface brightness jumps, associated with substantial temperature and density changes but approximate pressure equilibrium, implying that these are sloshing structures driven by a recent merger. A residual intensity image shows a core spiral structure closely matching that seen in the Perseus cluster, although at z = 0.055 the spiral pattern is less distinct. We infer that the most recent merger occurred 1–2 Gyr ago and that the event that brought the two observed super-massive black holes to the system core is even older. Under this interpretation, the black hole binary pair has, unusually, remained at a parsec-scale separation for more than 2 Gyr.

Orbital decay of supermassive black hole binaries in clumpy multiphase merger remnants

We simulate an equal-mass merger of two Milky Way-size galaxy discs with moderate gas fractions at parsec-scale resolution including a new model for radiative cooling and heating in a multi-phase medium, as well as star formation and feedback from supernovae. The two discs initially have a $2.6\times10^6\mathrm{~M_{\odot}}$ supermassive black hole (SMBH) embedded in their centers. As the merger completes and the two galactic cores merge, the SMBHs form a a pair with a separation of a few hundred pc that gradually decays. Due to the stochastic nature of the system immediately following the merger, the orbital plane of the binary is significantly perturbed. Furthermore, owing to the strong starburst the gas from the central region is completely evacuated, requiring $\sim10$~Myr for a nuclear disc to rebuild. Most importantly, the clumpy nature of the interstellar medium has a major impact on the the dynamical evolution of the SMBH pair, which undergo gravitational encounters with massive gas clouds and stochastic torquing by both clouds and spiral modes in the disk. These effects combine to greatly delay the decay of the two SMBHs to separations of a few parsecs by nearly two orders of magnitude, $\sim 10^8$ yr, compared to previous work. In mergers of more gas-rich, clumpier galaxies at high redshift stochastic torques will be even more pronounced and potentially lead to stronger modulation of the orbital decay. This suggests that SMBH pairs at separations of several tens of parsecs should be relatively common at any redshift.

Prospects for gravitational-wave detection and supermassive black hole astrophysics with pulsar timing arrays

Large-area sky surveys show that massive galaxies undergo at least one major merger in a Hubble time. Ongoing pulsar timing array (PTA) experiments are aimed at measuring the gravitational-wave (GW) emission from binary supermassive black holes (SMBHs) at the centres of galaxy merger remnants. In this paper, using the latest observational estimates for a range of galaxy properties and scaling relations, we predict the amplitude of the GW background generated by the binary SMBH population. We also predict the numbers of individual binary SMBH GW sources. We predict the characteristic strain amplitude of the GW background to lie in the range 5.1 × 10−16 < Ayr < 2.4 × 10−15 at a frequency of (1 yr)−1, with 95 per cent confidence. Higher values within this range, which correspond to the more commonly preferred choice of galaxy merger time-scale, will fall within the expected sensitivity ranges of existing PTA projects in the next few years. In contrast, we find that a PTA consisting of at least 100 pulsars observed with next-generation radio telescopes will be required to detect continuous-wave GWs from binary SMBHs. We further suggest that GW memory bursts from coalescing SMBH pairs are not viable sources for PTAs. Both the GW background and individual GW source counts are dominated by binaries formed in mergers between early-type galaxies of masses ≳5 × 1010 M⊙ at redshifts ≲1.5. Uncertainties in the galaxy merger time-scale and the SMBH mass–galaxy bulge mass relation dominate the uncertainty in our predictions.

Pair of supermassive black holes are orbiting each other at relativistic speeds

Pair of supermassive black holes are orbiting each other at relativistic speeds

Compact object mergers: observations of supermassive binary black holes and stellar tidal disruption events

AbstractThe capture and disruption of stars by supermassive black holes (SMBHs), and the formation and coalescence of binaries, are inevitable consequences of the presence of SMBHs at the cores of galaxies. Pairs of active galactic nuclei (AGN) and binary SMBHs are important stages in the evolution of galaxy mergers, and an intense search for these systems is currently ongoing. In the early and advanced stages of galaxy merging, observations of the triggering of accretion onto one or both BHs inform us about feedback processes and BH growth. Identification of the compact binary SMBHs at parsec and sub-parsec scales provides us with important constraints on the interaction processes that govern the shrinkage of the binary beyond the “final parsec”. Coalescing binary SMBHs are among the most powerful sources of gravitational waves (GWs) in the universe. Stellar tidal disruption events (TDEs) appear as luminous, transient, accretion flares when part of the stellar material is accreted by the SMBH. About 30 events have been identified by multi-wavelength observations by now, and they will be detected in the thousands in future ground-based or space-based transient surveys. The study of TDEs provides us with a variety of new astrophysical tools and applications, related to fundamental physics or astrophysics. Here, we provide a review of the current status of observations of SMBH pairs and binaries, and TDEs, and discuss astrophysical implications.

Gravitational-Wave Limits from Pulsar Timing Constrain Supermassive Black Hole Evolution

The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density (Ω(GW)) with 95% confidence to be Ω(GW)(H0/73 kilometers per second per megaparsec)(2) < 1.3 × 10(-9) (where H0 is the Hubble constant) at a frequency of 2.8 nanohertz, which is approximately a factor of 6 more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46 and 91%. For example, the standard galaxy formation model implemented in the Millennium Simulation Project is inconsistent with our limit with 50% probability.

STRUCTURE AND EVOLUTION OF CIRCUMBINARY DISKS AROUND SUPERMASSIVE BLACK HOLE BINARIES

It is generally believed that gaseous disks around supermassive black hole (SMBH) binaries in centers of galaxies can facilitate binary merger and give rise to observational signatures both in electromagnetic and gravitational wave domains. We explore general properties of circumbinary disks by reformulating standard equations for the viscous disk evolution in terms of the viscous angular momentum flux F_J. In steady state F_J is a linear function of the specific angular momentum, which is a generalization of (but is not equivalent to) the standard constant \dot M disk solution. If the torque produced by the central binary is effective at stopping gas inflow and opening a gap (or cavity) in the disk, then the inner part of the circumbinary disk can be approximated as a constant F_J disk. We compute properties of such disks in different physical regimes relevant for SMBH binaries and use these results to understand the gas-assisted evolution of SMBH pairs starting at separations 10^{-4}-10^{-2} pc. We find the following. (1) Pile-up of matter at the inner edge of the disk leads to continuous growth of the torque acting on the binary and can considerably accelerate its orbital evolution compared to the gravitational wave-driven decay. (2) Torque on the binary is determined non-locally and does not in general reflect the disk properties in the vicinity of the binary. (3) Binary evolution depends on the past history of the disk evolution. (4) Eddington limit can be important in circumbinary disks even if they accrete at sub-Eddington rates at late stages of binary evolution. (5) Circumbinary disk self-consistently evolved under the action of the binary torque emits more power and has spectrum different from the spectrum of constant \dot M disk - it is steeper (\nu F_\nu\propto \nu^{12/7}) and extends to shorter wavelength, facilitating its detection.

Double-peaked narrow-line signatures of dual supermassive black holes in galaxy merger simulations

We present a first attempt to model the narrow-line (NL) region of active galactic nuclei (AGN) in hydrodynamic simulations of galaxy mergers, using a novel physical prescription. This model is used to determine the origin of double-peaked NL AGN in merging galaxies and their connection to supermassive black hole (SMBH) pairs, motivated by recent observations of such objects. We find that double-peaked NLs induced by the relative motion of dual SMBHs are a generic but short-lived feature of gaseous major mergers. Double-peaked NL AGN should often be observed in late-stage galaxy mergers, during the kiloparsec-scale phase of SMBH inspiral or soon after the nuclear coalescence and subsequent SMBH merger. However, even within the kiloparsec-scale phase, only a minority of double-peaked NLs are directly induced by the relative motion of binary SMBHs; their lifetimes are typically a few Myr. The majority of double-peaked NLs result from gas kinematics near the SMBH, although prior to the SMBH merger up to ∼80 per cent of all double-peaked NL profiles may be influenced by SMBH motion via altered peak ratios or overall velocity offsets. The total lifetimes of double-peaked NL AGN depend strongly on viewing angle and on properties of the merging galaxies; gas-rich, nearly-equal-mass mergers have more NL AGN activity but may also be more obscured. Furthermore, in a typical merger, at least 10–40 per cent of the double-peaked NLs induced by SMBH motion have small projected separations, ∼0.1–1 kpc, making it difficult to clearly identify dual peaks of stellar surface brightness. Diffuse tidal features can indicate a late-stage merger, although they do not distinguish an SMBH pair from a merged SMBH. We demonstrate that double-peaked NL AGN spectra with large peak velocity splittings ( ≳ 500 km s−1) or with discernible overall velocity shifts are often associated with inspiraling SMBH pairs. Our results support the notion that selection of double-peaked NL AGN is a promising method for identifying dual SMBH candidates, but demonstrate the critical importance of high-resolution, multiwavelength follow-up observations, and the use of multiple lines of evidence, for confirming the dual nature of candidate SMBH pairs.