Testing the relativistic Doppler boost hypothesis for supermassive black hole binary candidates

Context. Recently we discovered an unprecedented supermassive black hole binary (SMBHB) candidate in the nearby Seyfert galaxy SDSS J1430+2303, which is predicted to merge within three years. X-ray spectroscopy may bring unique kinematic evidence for the last inspiraling stage, when the binary is too close to allow each of them to hold an individual broad line region. Aims. We try to confirm the unique SMBHB merger event and understand the associated high-energy processes from a comprehensive X-ray view. Methods. We observed SDSS J1430+2303 with XMM-Newton, NuSTAR, Chandra, and Swift spanning the first ∼200 days since its discovery. Results. X-ray variability, up to a factor of 7, has been detected on a timescale of a few days. The broadband spectrum from 0.2–70 keV can be well fitted with a model consisting of a power law and a relativistic reflection covered by a warm absorber. The properties of the warm absorber changed dramatically, for example, with a decrease in the line-of-sight velocity from ∼0.2c to ∼0.02c, between the two XMM-Newton observations separated by only 19 days, which can be naturally understood in the context of the SMBHB; although, the clumpy wind scenario cannot be completely excluded. Broad Fe Kα emission has been robustly detected, though its velocity shift or profile change is not yet measurable. Further longer X-ray observations are highly encouraged to detect the expected orbital motion of the binary.

Graham et al. discovered a supermassive black hole binary (SMBHB) candidate and identified the detected 5.2 yr period of the optical variability as the orbital period of the binary. Hydrodynamical simulations predict multiple periodic components for the variability of SMBHBs, thus raising the possibility that the true period of the binary is different from 5.2 yr. We analyse the periodogram of PG1302 and find no compelling evidence for additional peaks. We also point out that, despite the 5.2 yr peak being significant if a single source is considered, further analysis is required to account for the fact that PG1302 was selected among a large sample of 247 000 quasars. We derive upper limits on any additional periodic modulations in the available data, by modelling the light curve as the sum of stochastic noise and the known 5.2 yr periodic component, and injecting additional sinusoidal signals. We find that, with the current data, we would be able to detect with high significance (false alarm probability <1 per cent) secondary periodic terms, with periods in the range predicted by the simulations, if the amplitude of the variability was at least ∼0.06 mag (compared to 0.14 mag for the main sinusoid). A three-year follow-up monitoring campaign with weekly observations can increase the sensitivity for detecting secondary peaks by ≈50 per cent, and would allow a more robust test of predictions from hydrodynamical simulations.

The elusive supermassive black hole binaries (SMBHBs) are thought to be the penultimate stage of galaxy mergers, preceding a final coalescence phase. SMBHBs are sources of continuous gravitational waves, possibly detectable by pulsar timing arrays; the identification of candidates could help in performing targeted gravitational wave searches. Due to SMBHBs’ origin in the innermost parts of active galactic nuclei (AGN), X-rays are a promising tool for unveiling their presence, by means of either double Fe Kα emission lines or periodicity in their light curve. Here we report on a new method for selecting SMBHBs by means of the presence of a periodic signal in their Swift Burst Alert Telescope (BAT) 105 month light curves. Our technique is based on Fisher’s exact g-test and takes into account the possible presence of colored noise. Among the 553 AGN selected for our investigation, only the Seyfert 1.5 galaxy Mrk 915 emerges as a candidate SMBHB; from subsequent analysis of its light curve we find a period P 0 = 35 ± 2 months, and the null hypothesis is rejected at the 3.7σ confidence level. We also present a detailed analysis of the BAT light curve of the only previously X-ray-selected binary candidate source in the literature, the Seyfert 2 galaxy MCG+11-11-032. We find P 0 = 26.3 ± 0.6 months, consistent with the one inferred from previously reported double Fe Kα emission lines.

We present an intensive multiwavelength monitoring campaign of the quasar PG 1302−102 with Swift and the Las Cumbres Observatory network telescopes. At z ∼ 0.3, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time the RM technique has been applied to test disk structures predicted in the SMBH binary model that has been suggested for this source. PG 1302−102 was observed at a ∼daily cadence for ∼9 months in 14 bands spanning from X-ray to UV and optical wavelengths, and it shows moderate to significant levels of variability correlated between wavelengths. We measure the interband time lags, which are consistent with a τ ∝ λ 4/3 relation as expected from standard disk reprocessing, albeit with large uncertainties. The disk size implied by the lag spectrum is consistent with the expected disk size for its black hole mass within uncertainties. While the source resembles other reverberation-mapped active galactic nuclei in many respects, and we do not find evidence supporting the prevalent hypothesis that it hosts an SMBH binary, we demonstrate the feasibility of studying SMBH binaries from this novel angle and suggest possibilities for the LSST Deep Drilling Fields.

Supermassive black hole binaries (SBHBs) are a natural outcome of galaxy mergers, and they are expected to be among the loudest gravitational-wave sources at low frequencies. The source SDSS J2320+0024 was recently proposed as a promising SBHB candidate due to a possible periodicity in its light curve and variability in the MgII emission line. In this work, we reanalysed the optical ( gr , and i bands) light curves of J2320+0024 within the framework of Bayesian model selection. When periodicity was searched for together with red noise, analysis of the g -band light curve reveals a peak in the posterior of the period at ∼290 days. The posterior profile is too broad to yield a preference for periodic models over models that include only red noise. Furthermore, the same peak is not present in the analysis of the r -band and i -band light curve. A periodic model without red noise identified a different (∼1100 days) periodicity, but this model is statistically significantly disfavoured relative to the other models tested. In summary, we find no significant evidence in favour of a true periodic signal over red-noise variability. Our analysis questions the robustness of the previously proposed periodicity and emphasises the importance of rigorous statistical treatment. While our findings challenge the binary interpretation for J2320+0024, they do not rule it out. A statistically robust joint analysis of the photometric light curves and evolving broad-line profiles would shed further light on the true nature of this object.

Here we present a short overview and main results of our investigations of several effects which can induce shifts in the broad Fe K$\alpha$ line emitted from relativistic accretion disks around single and binary supermassive black holes. We used numerical simulations based on ray-tracing method in the Kerr metric to study the role of classical Doppler shift, special relativistic transverse Doppler shift and Doppler beaming, general relativistic gravitational redshift, and perturbations of the disk emissivity in the formation of the observed Fe K$\alpha$ line profiles. Besides, we also investigated whether the observed line profiles from the binary systems of supermassive black holes could be affected by the Doppler shifts due to dynamics of such systems. The presented results demonstrate that all these effects could have a significant influence on the observed profiles of the broad Fe K$\alpha$ line emitted from relativistic accretion disks around single and binary supermassive black holes.

We investigate the time-varying electromagnetic emission of a low-mass-ratio supermassive black hole binary (SMBHB) embedded in a circumprimary disc, with a particular interest in variability of shocks driven by the binary. We perform a 2D, locally isothermal hydrodynamics simulation of an SMBHB with mass ratio q = 0.01 and separation a = 100 Rg, using a physically self-consistent steady disc model. We estimate the electromagnetic variability from the system by monitoring accretion on to the secondary and using an artificial viscosity scheme to capture shocks and monitor the energy dissipated. The SMBHB produces a wide, eccentric gap in the disc, previously only observed for larger mass ratios, which we attribute to our disc model being much thinner (H/R ≈ 0.01 near the secondary) than is typical of previous works. The eccentric gap drives periodic accretion on to the secondary SMBH on a time-scale matching the orbital period of the binary, $t_{\rm {bin}}\approx 0.1\,\,\rm {yr}$, implying that the variable accretion regime of the SMBHB parameter space extends to lower mass ratios than previously established. Shocks driven by the binary are periodic, with a period matching the orbital period, and the shocks are correlated with the accretion rate, with peaks in the shock luminosity lagging peaks in the accretion rate by 0.43 tbin. We propose that the correlation of these quantities represents a useful identifier of SMBHB candidates, via observations of correlated variability in X-ray and ultraviolet monitoring of active galactic nuclei, rather than single-waveband periodicity alone.

We review the results of multi‐scale, hydrodynamical simulations of major mergers between galaxies with or without central supermassive black holes (SMBHs) to investigate the orbital decay of SMBH pairs in galactic nuclei and the formation of massive SMBH seeds via direct gas collapse. Both SPH simulations and AMR simulations are carried out. The complex balance between heating and cooling is modeled via an effective EOS with varying adiabatic index γ apparopriate for the conditions of an intense nuclear starburst such as that expected during and after the merger. Prominent gas inflows due to tidal torques produce nuclear disks at the centers of merger remnants whose properties depend sensitively on the details of gas thermodynamics. In parsec‐scale resolution simulations starting with two SMBHs originally at the centers of the two galaxies, a SMBH binary forms very rapidly, less than a million years after the merger of the two galaxies, owing to the drag exerted by the surrounding gaseous nuclear disk. Binary formation is significantly suppressed if heating, by e.g. radiative feedback from the accreting SMBHs, renders cooling negligible.The nuclear disk rearranges its mass distribution in response to a second, internal gas inflow occurring while the binary sinks. The inflow is driven by spiral instabilities imprinted by the final collision between the two galactic cores. In simulations with 0.1 pc resolution, the gas inflow continues all the way down to the center and peaks at >104 M⊙/yr, producing a Jeans‐unstable supermassive central cloud (with mass a few times 108 M⊙) only 105 yr after the merger. If the collapse continues, the cloud could form a massive black hole seed (Mseed>105 M⊙) after prior formation of a supermassive star or quasi‐star. The massive SMBH seed can grow up to a billion solar masses in less than a billion years by accreting the surrounding nuclear gas. If the gas‐rich merger occurs at z>8, this is then a new, attractive way to explain the rapid emergence of the bright QSOs discovered by the Sloan Digital Sky survey at z>6, which does not require the assumption of primordial gas composition in order to suppress cooling below 104 K and star formationas in models starting from unstable, isolated protogalactic disks. If there is a pre‐existing pair of SMBHs their orbital decay stalls at parsec scales because, as a result of the formation of the supermassive cloud, the nuclear disk density decreases outside the cloud, yielding much weaker dynamical friction. We envision a new scenario in which direct formation of massive black hole seeds and SMBH binary formation are mutually exclusive; if a SMBH is already present in the nuclear disk it can stabilize it and weaken the secondary inflow via its energetic feedback, maintaining a high enough density where the SMBHs are located and assisting their sinking.

The couplings between supermassive black hole binaries (SMBHBs) and their environments within galactic nuclei have been well studied as part of the search for solutions to the final parsec problem. The scattering of stars by the binary or the interaction with a circumbinary disk may efficiently drive the system to sub-parsec separations, allowing the binary to enter a regime where the emission of gravitational waves can drive it to merger within a Hubble time. However, these interactions can also affect the orbital parameters of the binary. In particular, they may drive an increase in binary eccentricity which survives until the system’s gravitational-wave (GW) signal enters the pulsar-timing array (PTA) band. Therefore, if we can measure the eccentricity from observed signals, we can potentially deduce some of the properties of the binary environment. To this end, we build on previous techniques to present a general Bayesian pipeline with which we can detect and estimate the parameters of an eccentric SMBHB system with PTAs. Additionally, we generalize the PTA -statistic to eccentric systems, and show that both this statistic and the Bayesian pipeline are robust when studying circular or arbitrarily eccentric systems. We explore how eccentricity influences the detection prospects of single GW sources, as well as the detection penalty incurred by employing a circular waveform template to search for eccentric signals, and conclude by identifying important avenues for future study.

A technique is described that permits the robust decomposition of the bulge and disk components of a sample of Seyfert galaxies, as well as a (control) sample of non-active galaxies matched to the Seyferts in the distributions of redshift, luminosity and morphological classification. The structural parameters of the host galaxies in both samples are measured. No statistically significant differences at greater than the 95% level are found in these parameters according to a K-S test. ``Companion galaxies'' --- defined as any galaxy within a projected separation of 200 h^{-1} kpc from the center of the host --- are identified and their basic properties measured. A comparison between the active and control samples in the distributions of apparent R magnitude, absolute R magnitude (assuming the ``companions'' are at the distance of the host), projected separation from the host, position angle relative to the host, magnitude difference between the companion and host, and strength of the tidal parameter, show no statistically significant differences. Similarly, no statistically significant differences are found between the control and active sample host galaxies in terms of light asymmetries --- bars, rings, isophotal twisting, etc. The implications for a model in which interactions and mergers are responsible for inciting activity in galactic nuclei are discussed briefly.

A star that wanders too close to a massive black hole (BH) is shredded by the BH's tidal gravity. Stellar gas falls back to the BH, releasing a flare of energy. In anticipation of upcoming transient surveys, we predict the light curves and spectra of tidal flares as a function of time, highlighting the unique signatures of tidal flares in the optical and near-IR. Some of the gas initially bound to the BH is likely blown away when the fallback rate is super-Eddington at early times. This outflow produces an optical luminosity comparable to that of a supernova (Figure 1, left panel); such events have durations of ~ 10 days and may have been missed in supernova searches that exclude the nuclear regions of galaxies. When the fallback rate subsides below Eddington, the gas accretes onto the BH via a thin disk whose emission peaks in the UV to soft X-rays. Some of this emission is reprocessed by the unbound stellar debris, producing a spectrum of very broad emission lines, with no corresponding narrow forbidden lines (center panel). These lines are strongest for BHs with MBH ~ 105–106M⊙ and thus optical surveys are particularly sensitive to the lowest mass BHs in galactic nuclei. Calibrating our models to ROSAT and GALEX observations, we predict detection rates for Pan-STARRS, Palomar Transit Factory, and LSST (right panel) and highlight observational challenges in the optical. Pan-STARRS should detect at least several events per year — many more if current theoretical models of super-Eddington outflows are correct. These surveys will significantly improve our knowledge of stellar dynamics in galactic nuclei, the physics of super-Eddington accretion, the demography of intermediate mass BHs, and the role of tidal disruption in the growth of massive BHs.

Context. Galaxy interactions and mergers can lead to supermassive black hole (SMBH) binaries, which become active galactic nucleus (AGN) pairs when the SMBHs start accreting mass. If there is a third galaxy involved in the interaction, then a triple-AGN system can form. Aims. Our goal is to investigate the nature of the nuclear emission from the galaxies in the interacting pair NGC 7733–NGC 7734 using archival VLT/MUSE integral field spectrograph data and study its relation to the stellar mass distribution traced by near-infrared (NIR) observations from the South African Astronomical Observatory (SAAO). Methods. We conducted NIR observations using the SAAO and identified the morphological properties of bulges in each galaxy. We used MUSE data to obtain a set of ionized emission lines from each galaxy and studied the ionization mechanism. We also examined the relation of the galaxy pair with any nearby companions with far-ultraviolet observations using the UVIT. Results. The emission line analysis from the central regions of NGC 7733 and NGC 7734 shows Seyfert and low ionization nuclear emission-line regions type AGN activity. The galaxy pair NGC 7733−34 also shows evidence of a third component, which has Seyfert-like emission. Hence, the galaxy pair NGC 7733−34 forms a triple-AGN system. We also detected an extended narrow-line region associated with the nucleus of NGC 7733.

The observed relation between the X-ray radiation from AGNs, originating in the corona, and the optical/UV radiation from the disk is usually described by the anticorrelation between the UV to X-ray slope alpha_ox and the UV luminosity. Many factors can affect this relation, including: enhanced X-ray emission associated with the jets of radio-loud AGNs; X-ray absorption associated with the UV Broad Absorption Line (BAL) outflows; other X-ray absorption not associated with BALs; intrinsic X-ray weakness; UV and X-ray variability, and non-simultaneity of UV and X-ray observations. The separation of these effects provides information about the intrinsic alpha_ox-L_UV relation and its dispersion, constraining models of disk-corona coupling. We extract simultaneous data from the second XMM-Newton serendipitous source catalogue and the XMM-Newton Optical Monitor Serendipitous UV Source Survey Catalog, and derive the single-epoch alpha_ox indices. We use ensemble structure functions to analyse multi-epoch data. We confirm the anticorrelation of alpha_ox with L_UV, and do not find any evidence of a dependence of alpha_ox on z. The dispersion in our simultaneous data (0.12) is not significantly smaller than in previous non-simultaneous studies, suggesting that "artificial alpha_ox variability" introduced by non-simultaneity is not the main cause of dispersion. "Intrinsic alpha_ox variability", i.e., the true variability of the X-ray to optical ratio, is instead important, and accounts for ~30% of the total variance, or more. "Inter-source dispersion", due to intrinsic differences in the average alpha_ox values from source to source, is also important. The dispersion introduced by variability is mostly caused by the long timescale variations, which are expected to be driven by the optical variations.

We calculate the local UV flux density in the GALEX MIS FUV and NUV bands using redshifts provided by SDSS DR7. Luminosity functions are calculated for the overlapping MIS and SDSS sample, allowing flux densities to be measured and the local star formation rate (SFR) to be calculated using volumes much larger than previous FUV based estimates. We calculate flux densities for a dust corrected low redshift (0.013 < z < 0.1) sample of f[FUV] = 22.24 \pm 3.13 \times 10^25 h ergs s^-1 Hz^-1 Mpc^-3, f[NUV] = 38.54 \pm 5.30 \times 10^25 h ergs s^-1 Hz^-1 Mpc^-3. The star formation rate density found is 0.0312\pm0.0045 h M\odot yr^-1 Mpc^-3. This is larger than published rates recently found using the UV implied SFR, though the major discrepancy is the correction made for dust attenuation and once this is dealt with consistently the results agree well. These values are also consistent with recent Halpha derived SFRs. Once cosmic variance is taken into account most of the recent SFRs at low redshift (z < 0.3) found in the literature can be brought into agreement, however the lowest redshift values (z < 0.045) do appear to be significantly lower.

Images of the central and southern parts of the local group Sb spiral galaxy M31 were obtained by the Ultraviolet Imaging Telescope during the Astro-1 Spacelab mission. Stellar photometry is obtained for 611 stars in 59 van den Bergh associations in the near-UV A1 band (central wavelngth 2490 A), and for 130 of these stars in the far UV B1 band (central wavlength 1520 A). In addition, B-band magnitudes for 456 of the stars in approximately 30 associations, together with U- and R-band photometry for approximately 290 stars in approximately 15 associations are determined from CCD images. B-magnitudes from the catalog of Magnier et al. (1992) are given for 137 stars and R-magnitudes for 233 stars. Stellar positions are obtained from the ground based CCD images and from UIT images. Total fluxes in the UV bands are measured for 63 associations . Finding charts are presented in the NUV and FUV bands, as well as, one optical band. Tables give the images utilized in this investigation, the associations observed, with positions, total UV fluxes, extinctions estimated from color-magnitude diagrams, and stellar positions with photometry in up to five bands. UV color magnitude diagrams are presented with extinctions estimated by comparing the observed association color magnitude diagrams with colors and magnitudes of early-type stars computed from the IUE spectral atlas of Fanelli et al. (1992). Comparison with evolutionary models suggests a maximum initial mass approx. 85 solar mass. Ultraviolet extinctions within most M31 OB associations is found to follow an extinction curve similar to that of Hutchings et al. (1992). The variation in the fraction of near-UV stars also measured in the far-UV among the associations is probably caused by variations in age, extinction, and the degree of foreground contamination.