Galaxy Merger Remnants Research Articles

Massive black hole binaries as sources of low-frequency gravitational waves and X-shaped radio galaxies

ABSTRACT We present the study of multimessenger signatures of massive black hole (MBH) binaries residing in the centres of galaxy merger remnants. In particular, we first focus on the gravitational wave background (GWB) produced by an ensemble of MBH binary inspirals in the frequency range probed by the Pulsar Timing Array (PTA) experiments. The improved estimates of the characteristic strain were obtained with the inclusion of environmental effects on the MBH binary orbital decay within the galaxy merger remnants, added in post-processing to the semi-analytical model of galaxy formation and evolution SHARK. Secondly, we explore two, intriguing in terms of the MBH binary evolution studies, hypotheses aiming to explain the origins of X-shaped radio galaxies – a peculiar type of objects with double lobe structures, constituting approximately 6–10 per cent of known radio loud galaxies. The two considered scenarios involve either an abrupt change in the jet direction after an MBH merger (a spin-flip) or an unresolved close binary, where each of the two components produces a jet. We find that the estimated GWB amplitude at the reference frequency $f_0=1 \, {\rm yr}^{-1}$ is in the range of $A_{\rm { yr^{-1}}} = 1.20\times 10^{-15}{\!-\!}1.46\times 10^{-15}$, which is 50 per cent lower than the strain of the signal detected by the PTA experiments. We also show that the spin-flip scenario considered in gas-poor mergers reproduces the observed properties of X-shaped radio galaxies well in terms of flip angle, redshift, and luminosity distributions.

Galaxy mergers can rapidly shut down star formation

ABSTRACT Galaxy mergers trigger both star formation and accretion on to the central supermassive black hole. As a result of subsequent energetic feedback processes, it has long been proposed that star formation may be promptly extinguished in galaxy merger remnants. However, this prediction of widespread, rapid quenching in late stage mergers has been recently called into question with modern simulations and has never been tested observationally. Here we perform the first empirical assessment of the long-predicted end phase in the merger sequence. Based on a sample of ∼500 post-mergers identified from the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), we show that the frequency of post-merger galaxies that have rapidly shutdown their star formation following a previous starburst is 30–60 times higher than expected from a control sample of non-merging galaxies. No such excess is found in a sample of close galaxy pairs, demonstrating that mergers can indeed lead to a rapid halt to star formation, but that this process only manifests after coalescence.

The complex globular cluster system of the S0 galaxy NGC 4382 in the outskirts of the Virgo Cluster

ABSTRACT NGC 4382 is a merger-remnant galaxy that has been classified as morphological type E2, S0, and even Sa. In this work, we performed a photometric and spectroscopic analysis of the globular cluster (GC) system of this peculiar galaxy in order to provide additional information about its history. We used a combination of photometric data in different filters, and multiobject and long-slit spectroscopic data obtained using the Gemini/GMOS instrument. The photometric analysis of the GC system, using the Gaussian Mixture Model algorithm in the colour plane, reveals a complex colour distribution within Rgal < 5 arcmin (26.1 kpc), showing four different groups: the typical blue and red subpopulations, a group with intermediate colours, and the fourth group towards even redder colours. From the spectroscopic analysis of 47 GCs, confirmed members of NGC 4382 based on radial velocities, we verified 3 of the 4 photometric groups from the analysis of their stellar populations using the ULySS code. NGC 4382 presents the classic blue (10.4 ± 2.8 Gyr, [Fe/H] = −1.48 ± 0.18 dex), and red (12.1 ± 2.3 Gyr, [Fe/H] = −0.64 ± 0.26 dex) GCs formed earlier in the lifetime of the galaxy, and a third group of young GCs (2.2 ± 0.9 Gyr; [Fe/H] = −0.05 ± 0.28 dex). Finally, analysis of long-slit data of the galaxy reveals a luminosity-weighted mean age for the stellar population of ∼2.7 Gyr, and an increasing metallicity from [Fe/H] = −0.1 to +0.2 dex in Rgal < 10 arcsec (0.87 kpc). These values, and other morphological signatures in the galaxy, are in good agreement with the younger group of GCs, indicating a common origin as a result of a recent merger.

The combined and respective roles of imaging and stellar kinematics in identifying galaxy merger remnants

ABSTRACT One of the central challenges to establishing the role of mergers in galaxy evolution is the selection of pure and complete merger samples in observations. In particular, while large and reasonably pure interacting galaxy pair samples can be obtained with relative ease via spectroscopic criteria, automated selection of post-coalescence merger remnants is restricted to the physical characteristics of remnants alone. Furthermore, such selection has predominantly focused on imaging data – whereas kinematic data may offer a complimentary basis for identifying merger remnants. Therefore, we examine the theoretical utility of both the morphological and kinematic features of merger remnants in distinguishing galaxy merger remnants from other galaxies. Deep classification models are calibrated and evaluated using idealized synthetic images and line-of-sight stellar velocity maps of a heterogeneous population of galaxies and merger remnants from the TNG100 cosmological hydrodynamical simulation. We show that even idealized stellar kinematic data have limited utility compared to imaging and underperforms by $2.1 \pm 0.5 {{\ \rm per\ cent}}$ in completeness and $4.7 \pm 0.4 {{\ \rm per\ cent}}$ in purity for our fiducial model architecture. Combining imaging and stellar kinematics offers a small boost in completeness (by $1.8 \pm 0.4 {{\ \rm per\ cent}}$, compared to $92.7 \pm 0.2 {{\ \rm per\ cent}}$ from imaging alone) but no change in purity ($0.1\pm 0.3 {{\ \rm per\ cent}}$ improvement compared to $92.7 \pm 0.2 {{\ \rm per\ cent}}$, evaluated with equal numbers of merger remnant and non-remnant control galaxies). Classification accuracy of all models is particularly sensitive to physical companions at separations ≲ 40 kpc and to time-since-coalescence. Taken together, our results show that the stellar kinematic data have little to offer in compliment to imaging for merger remnant identification in a heterogeneous galaxy population.

Cold Molecular Gas in Merger Remnants. II. The Properties of Dense Molecular Gas

Abstract We present the 3 mm wavelength spectra of 28 local galaxy merger remnants obtained with the Large Millimeter Telescope. Sixteen molecular lines from 14 different molecular species and isotopologues were identified, and 21 out of 28 sources were detected in one or more molecular lines. On average, the line ratios of the dense gas tracers, such as HCN (1–0) and HCO+(1–0), to 13CO (1–0) are 3–4 times higher in ultra/luminous infrared galaxies (U/LIRGs) than in non-LIRGs in our sample. These high line ratios could be explained by the deficiency of 13CO and high dense gas fractions suggested by high HCN (1–0)/12CO (1–0) ratios. We calculate the IR-to-HCN (1–0) luminosity ratio as a proxy of the dense gas star formation efficiency. There is no correlation between the IR/HCN ratio and the IR luminosity, while the IR/HCN ratio varies from source to source ((1.1–6.5) × 103 L ☉/(K km s−1 pc2)). Compared with the control sample, we find that the average IR/HCN ratio of the merger remnants is higher by a factor of 2–3 than those of the early/mid-stage mergers and nonmerging LIRGs, and it is comparable to that of the late-stage mergers. The IR-to-12CO (1–0) ratios show a similar trend to the IR/HCN ratios. These results suggest that star formation efficiency is enhanced by the merging process and maintained at high levels even after the final coalescence. The dynamical interactions and mergers could change the star formation mode and continue to impact the star formation properties of the gas in the postmerger phase.

The Detectability of Kiloparsec-scale Dual Active Galactic Nuclei: The Impact of Galactic Structure and Black Hole Orbital Properties

Abstract Observational searches for dual active galactic nuclei (dAGNs) at kiloparsec separations are crucial for understanding the role of galaxy mergers in the evolution of galaxies. In addition, kiloparsec-scale dAGNs may serve as the parent population of merging massive black hole (MBH) binaries, an important source of gravitational waves. We use a semi-analytical model to describe the orbital evolution of unequal mass MBH pairs under the influence of stellar and gaseous dynamical friction in post-merger galaxies. We quantify how the detectability of approximately 40,000 kpc-scale dAGNs depends on the structure of their host galaxies and the orbital properties of the MBH pair. Our models indicate that kiloparsec-scale dAGNs are most likely to be detected in gas-rich post-merger galaxies with smaller stellar bulges and relatively massive, rapidly rotating gas disks. The detectability is also increased in systems with MBHs of comparable masses following low eccentricity prograde orbits. In contrast, dAGNs with retrograde, low eccentricity orbits are some of the least detectable systems among our models. The dAGNs in models in which the accreting MBHs are allowed to exhibit radiative feedback are characterized by a significantly lower overall detectability. The suppression in detectability is most pronounced in gas-rich merger remnant galaxies, where radiation feedback is more likely to arise. If so, then large, relatively gas-poor galaxies may be the best candidates for detecting dAGNs.

The Pairing Probability of Massive Black Holes in Merger Galaxies in the Presence of Radiative Feedback

Abstract Dynamical friction (DF) against stars and gas is thought to be an important mechanism for orbital evolution of massive black holes (MBHs) in merger remnant galaxies. Recent theoretical investigations, however, show that DF does not always lead to MBH inspiral. For MBHs evolving in gas-rich backgrounds, the ionizing radiation that emerges from the innermost parts of their accretion flow can affect the surrounding gas in such a way to cause the MBHs to accelerate and gain orbital energy. This effect was dubbed “negative DF.” We use a semianalytic model to study the impact of negative DF on pairs of MBHs in merger remnant galaxies evolving under the combined influence of stellar and gaseous DF. Our results show that for a wide range of merger galaxy and MBH properties negative DF reduces the MBH pairing probability by ∼46%. The suppression of MBH pairing is most severe in galaxies with one or more of these properties: (1) a gas fraction of f g ≥ 0.1; (2) a galactic gas disk rotating close to the circular velocity; (3) MBH pairs in prograde, low eccentricity orbits; and (4) MBH pairs with mass <108 M ⊙. The last point is of importance because MBH pairs in this mass range are direct progenitors of merging binaries targeted by the future space-based gravitational wave observatory the Laser Interferometer Space Antenna (LISA).

Deviations from the Infrared-radio Correlation in Massive, Ultracompact Starburst Galaxies

Abstract Feedback through energetic outflows has emerged as a key physical process responsible for transforming star-forming galaxies into the quiescent systems observed in the local universe. To explore this process, this paper focuses on a sample of massive and compact merger remnant galaxies hosting high-velocity gaseous outflows ( km s−1), found at intermediate redshift (z ∼ 0.6). From their mid-infrared emission and compact morphologies, these galaxies are estimated to have exceptionally large star formation rate (SFR) surface densities (ΣSFR ∼ 103 M ⊙ yr−1 kpc−2), approaching the Eddington limit for radiation pressure on dust grains. This suggests that star formation feedback may be driving the observed outflows. However, these SFR estimates suffer from significant uncertainties. We therefore sought an independent tracer of star formation to probe the compact starburst activity in these systems. In this paper, we present SFR estimates calculated using 1.5 GHz continuum Jansky Very Large Array observations for 19 of these galaxies. We also present updated infrared (IR) SFRs calculated from WISE survey data. We estimate SFRs from the IR to be larger than those from the radio for 16 out of 19 galaxies by a median factor of 2.5. We find that this deviation is maximized for the most compact galaxies hosting the youngest stellar populations, suggesting that compact starbursts deviate from the IR-radio correlation. We suggest that this deviation stems either from free–free absorption of synchrotron emission, a difference in the timescale over which each indicator traces star formation, or exceptionally hot IR-emitting dust in these ultra-dense galaxies.

Pairing of Massive Black Holes in Merger Galaxies Driven by Dynamical Friction

Abstract Motivated by observational searches for massive black hole (MBH) pairs at kiloparsec separations we develop a semianalytic model to describe their orbital evolution under the influence of stellar and gaseous dynamical friction (DF). The goal of this study is to determine how the properties of the merger remnant galaxy and the MBHs affect the likelihood and timescale for formation of a close MBH pair with separation of ≲1 pc. We compute approximately 40,000 configurations that cover a wide range of host galaxy properties and investigate their impact on the orbital evolution of unequal mass MBH pairs. We find that the percentage for MBH pairing within a Hubble time is larger than 80% in remnant galaxies with a gas fraction <20% and in galaxies hosting MBH pairs with total mass >106 M ⊙ and mass ratios ≥1/4. Among these, the remnant galaxies characterized by the fastest formation of close, gravitationally bound MBHs have one or more of the following properties: (1) a large stellar bulge, (2) comparable mass MBHs, and (3) a galactic gas disk rotating close to the circular speed. In such galaxies, the MBHs with the shortest inspiral times, which are likely progenitors of coalescing MBHs, are either on circular prograde orbits or on very eccentric retrograde orbits. Our model also indicates that remnant galaxies with opposite properties, which host slowly evolving MBH pairs, are the most likely hosts of dual active galactic nuclei at kiloparsec separations.

Black hole fueling and coalescence in galaxy mergers

AbstractWe study the gas accretion and orbital decay of super‐massive black hole (SMBH) pairs in the nuclear region of galaxy merger remnants, where a large amount of gas is expected to accumulate. For that purpose, we perform numerical simulations with both smooth particle hydro‐dynamics (SPH) (Gadget) and adaptive mesh refinement (AMR) (Ramses) techniques. We study the places of the highest SMBH gas accretion, finding that peaks of mass accretion rates are near pericenter passages. We also study the conditions for fast migration and coalescence of SMBHs in galactic nuclei, finding that, in the typical conditions of relatively wet mergers, the SMBHs will experience a fast shrinking, down to the scales where their final coalescence will be driven by the emission of gravitational waves, unless other processes strongly affect the distribution of the gas near the SMBH binary, such as active galactic nuclei (AGN) feedback.

The Hot Gas Exhaust of Starburst Engines in Mergers: Testing Models of Stellar Feedback and Star Formation Regulation

Abstract Using archival data from the Chandra X-ray telescope, we have measured the spatial extent of the hot interstellar gas in a sample of 49 nearby interacting galaxy pairs, mergers, and merger remnants. For systems with SFR > 1 M ☉ yr−1, the volume and mass of hot gas are strongly and linearly correlated with the star formation rate (SFR). This supports the idea that stellar/supernovae feedback dominates the production of hot gas in these galaxies. We compared the mass of X-ray-emitting hot gas M X(gas) with the molecular and atomic hydrogen interstellar gas masses in these galaxies ( and M H i , respectively), using published carbon monoxide and 21 cm H i measurements. Systems with higher SFRs have larger M X(gas)/( + M H i ) ratios on average, in agreement with recent numerical simulations of star formation and feedback in merging galaxies. The M X(gas)/( + M H i ) ratio also increases with dust temperature on average. The ratio M X(gas)/SFR is anticorrelated with the Infrared Astronomical Satellite 60–100 μm flux ratio and with the Spitzer 3.6–24 μm color. These trends may be due to variations in the spatial density of young stars, the stellar age, the ratio of young to old stars, the initial mass function, and/or the efficiency of stellar feedback. Galaxies with low SFR (<1 M ☉ yr−1) and high K band luminosities may have an excess of hot gas relative to the relation for higher SFR galaxies, while galaxies with low K band luminosities (and therefore low stellar masses) may have a deficiency in hot gas, but our sample is not large enough for strong statistical significance.

Recoiling supermassive black holes in analytical and numerical galaxy potential

ABSTRACTWe follow trajectories of recoiling supermassive black holes (SMBHs) in analytical and numerical models of galaxy merger remnants with masses of 1011 M⊙ and 1012 M⊙. We construct various merger remnant galaxies in order to investigate how the central SMBH mass and the mass ratio of progenitor galaxies influence escape velocities of recoiling SMBHs. Our results show that static analytical models of major merger remnant galaxies overestimate the SMBHs escape velocities. During major mergers violent relaxation leads to the decrease of galaxy mass and lower potential at large remnant radii. This process is not depicted in static analytical potential but clearly seen in our numerical models. Thus, the evolving numerical model is a more realistic description of dynamical processes in galaxies with merging SMBHs. We find that SMBH escape velocities in numerical major merger remnant galaxies can be up to 25 per cent lower compared to those in analytical models. Consequently, SMBHs in numerical models generally reach greater galactocentric distances and spend more time on bound orbits outside of the galactic nuclei. Thus, numerical models predict a greater number of spatially offset active galactic nuclei (AGNs).

Tidal disruption event rates in galaxy merger remnants

ABSTRACT The rate of tidal disruption events (TDEs) depends sensitively on the stellar properties of the central galactic regions. Simulations show that galaxy mergers cause gas inflows, triggering nuclear starbursts, increasing the central stellar density. Motivated by these numerical results, and by the observed overrepresentation of post-starburst galaxies among TDE hosts, we study the evolution of the TDE rate in high-resolution hydrodynamical simulations of a galaxy merger, in which we capture the evolution of the stellar density around the massive black holes (BHs). We apply analytical estimates of the loss-cone theory, using the stellar density profiles from simulations, to estimate the time evolution of the TDE rate. At the second pericentre, a nuclear starburst enhances the stellar density around the BH in the least massive galaxy, leading to an enhancement of the TDE rate around the secondary BH, although the magnitude and the duration of the increase depend on the stochasticity of star formation on very small scales. The central stellar density around the primary BH remains instead fairly constant, and so is its TDE rate. After the formation of the binary, the stellar density decreases, and so does the TDE rate.

Star Clusters in the Elliptical Galaxy NGC 4589 Hosting a Calcium-rich SN Ib (SN 2005cz)

Abstract NGC 4589, a bright E2 merger-remnant galaxy, hosts the peculiar fast and faint calcium-rich SN Ib SN 2005cz. The progenitor of Ca-rich SNe Ib has been controversial: it could be (1) a young, massive star with 6–12 M ⊙ in a binary system, or (2) an old, low-mass star in a binary system that was kicked out from the galaxy center. Moreover, previous distance estimates for this galaxy have shown a large spread, ranging from 20 to 60 Mpc. Thus, using archival Hubble Space Telescope/ACS F435W, F555W, and F814W images, we search for star clusters in NGC 4589 in order to help resolve these issues. We find a small population of young star clusters with 25 < V ≤ 27 (−7.1 < M V ≤ −5.1) mag and age <1 Gyr in the central region at R < 0.′5 (<3.8 kpc), thus supporting the massive-star progenitor scenario for SN 2005cz. In addition to young star clusters, we also find a large population of old globular clusters. In contrast to previous results in the literature, we find that the color distribution of the globular clusters is clearly bimodal. The turnover (Vega) magnitude in the V-band luminosity functions of the blue (metal-poor) globular clusters is determined to be V 0(max) = 24.40 ± 0.10 mag. We derive the total number of globular clusters, N GC = 640 ± 50, and the specific frequency, S N = 1.7 ± 0.2. Adopting a calibration for the metal-poor globular clusters, M V (max) = −7.66 ± 0.14 mag, we derive a distance to this galaxy: (m − M)0 = 32.06 ± 0.10(ran) ± 0.15(sys) (d = 25.8 ± 2.2 Mpc).

CALIFA reveals prolate rotation in massive early-type galaxies: A polar galaxy merger origin?

We present new evidence for eight early-type galaxies (ETGs) from the CALIFA Survey that show clear rotation around their major photometric axis ("prolate rotation"). These are LSBCF560-04, NGC 0647, NGC 0810, NGC 2484, NGC 4874, NGC 5216, NGC 6173 and NGC 6338. Including NGC 5485, a known case of an ETG with stellar prolate rotation, as well as UGC 10695, a further possible candidate for prolate rotation, we report ten CALIFA galaxies in total that show evidence for such a feature in their stellar kinematics. Prolate rotators correspond to ~9% of the volume-corrected sample of CALIFA ETGs, a fraction much higher than previously reported. We find that prolate rotation is more common among the most massive ETGs. We investigate the implications of these findings by studying N-body merger simulations, and show that a prolate ETG with rotation around its major axis could be the result of a major polar merger, with the amplitude of prolate rotation depending on the initial bulge-to-total stellar mass ratio of its progenitor galaxies. Additionally, we find that prolate ETGs resulting from this formation scenario show a correlation between their stellar line-of-sight velocity and higher order moment h_3, opposite to typical oblate ETGs, as well as a double peak of their stellar velocity dispersion along their minor axis. Finally, we investigate the origin of prolate rotation in polar galaxy merger remnants. Our findings suggest that prolate rotation in massive ETGs might be more common than previously expected, and can help towards a better understanding of their dynamical structure and formation origin.

Infalling clouds on to supermassive black hole binaries – II. Binary evolution and the final parsec problem

The formation of massive black hole binaries (MBHBs) is an unavoidable outcome of galaxy evolution via successive mergers. However, the mechanism that drives their orbital evolution from parsec separations down to the gravitational wave (GW) dominated regime is poorly understood, and their final fate is still unclear. If such binaries are embedded in gas-rich and turbulent environments, as observed in remnants of galaxy mergers, the interaction with gas clumps (such as molecular clouds) may efficiently drive their orbital evolution. Using numerical simulations, we test this hypothesis by studying the dynamical evolution of an equal-mass, circular MBHB accreting infalling molecular clouds. We investigate different orbital configurations, modelling a total of 13 systems to explore different possible impact parameters and relative inclinations of the cloud-binary encounter. We focus our study on the prompt, transient phase during the first few orbits when the dynamical evolution of the binary is fastest, finding that this evolution is dominated by the exchange of angular momentum through gas capture by the individual black holes and accretion. Building on these results, we construct a simple model for evolving a MBHB interacting with a sequence of clouds, which are randomly drawn from reasonable populations with different levels of anisotropy in their angular momenta distributions. We show that the binary efficiently evolves down to the GW emission regime within a few hundred million years, overcoming the 'final parsec' problem regardless of the stellar distribution.

A NEW CHANNEL FOR THE FORMATION OF KINEMATICALLY DECOUPLED CORES IN EARLY-TYPE GALAXIES

We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting from a major merger simulation of two disk galaxies. We show that although the two progenitor galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their final coalescence. This results in a central kinematic decoupling and the formation of a large-scale (~2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant (M*=1.3x10^11 Msun), while its outer parts keep the rotation direction of the initial orbital spin. The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical off-centered peaks, comparable to the observed "2-sigma galaxies". The KDC/CRC consists mainly of old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde mergers is in addition to previously known formation scenarios from retrograde mergers and can help towards explaining the substantial fraction of KDCs observed in early-type galaxies.

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.

MASSIVE BLACK HOLE PAIRS IN CLUMPY, SELF-GRAVITATING CIRCUMNUCLEAR DISKS: STOCHASTIC ORBITAL DECAY

We study the dynamics of massive black hole pairs in clumpy gaseous circumnuclear disks. We track the orbital decay of the light, secondary black hole $M_{\bullet2}$ orbiting around the more massive primary at the center of the disk, using $N$-body/smoothed particle hydrodynamic simulations. We find that the gravitational interaction of $M_{\bullet2}$ with massive clumps $M_{\rm cl}$ erratically perturbs the otherwise smooth orbital decay. In close encounters with massive clumps, gravitational slingshots can kick the secondary black hole out of the disk plane. The black hole moving on an inclined orbit then experiences the weaker dynamical friction of the stellar background, resulting in a longer orbital decay timescale. Interactions between clumps can also favor orbital decay when the black hole is captured by a massive clump which is segregating toward the center of the disk. The stochastic behavior of the black hole orbit emerges mainly when the ratio $M_{\bullet2}/M_{\rm cl}$ falls below unity, with decay timescales ranging from $\sim1$ to $\sim50$ Myr. This suggests that describing the cold clumpy phase of the inter-stellar medium in self-consistent simulations of galaxy mergers, albeit so far neglected, is important to predict the black hole dynamics in galaxy merger remnants.

DUAL SUPERMASSIVE BLACK HOLE CANDIDATES IN THE AGN AND GALAXY EVOLUTION SURVEY

Dual supermassive black holes (SMBHs) with kiloparsec scale separations in merger-remnant galaxies are informative tracers of galaxy evolution, but the avenue for identifying them in large numbers for such studies is not yet clear. One promising approach is to target spectroscopic signatures of systems where both SMBHs are fueled as dual active galactic nuclei (AGNs), or where one SMBH is fueled as an offset AGN. Dual AGNs may produce double-peaked narrow AGN emission lines, while offset AGNs may produce single-peaked narrow AGN emission lines with line-of-sight velocity offsets relative to the host galaxy. We search for such dual and offset systems among 173 Type 2 AGNs at z<0.37 in the AGN and Galaxy Evolution Survey (AGES), and we find two double-peaked AGNs and five offset AGN candidates. When we compare these results to a similar search of the DEEP2 Galaxy Redshift Survey and match the two samples in color, absolute magnitude, and minimum velocity offset, we find that the fraction of AGNs that are dual SMBH candidates increases from z=0.25 to z=0.7 by a factor of ~6 (from 2/70 to 16/91, or 2.9% to 18%). This may be associated with the rise in the galaxy merger fraction over the same cosmic time. As further evidence for a link with galaxy mergers, the AGES offset and dual AGN candidates are tentatively ~3 times more likely than the overall AGN population to reside in a host galaxy that has a companion galaxy (from 16/173 to 2/7, or 9% to 29%). Follow-up observations of the seven offset and dual AGN candidates in AGES will definitively distinguish velocity offsets produced by dual SMBHs from those produced by narrow-line region kinematics, and will help sharpen our observational approach to detecting dual SMBHs.