We present a numerical framework for the variability of active galactic nuclei (AGN), which links the variability of AGN over a broad range of timescales and luminosities to the observed properties of the AGN population as a whole, and particularly the Eddington ratio distribution function (ERDF). We have implemented our framework on GPU architecture, relying on previously published time series generating algorithms. After extensive tests that characterise several intrinsic and numerical aspects of the simulations, we describe some applications used for current and future time domain surveys and for the study of extremely variable sources (e.g., "changing look" or flaring AGN). Specifically, we define a simulation setup which reproduces the AGN variability observed in the PTF/iPTF survey, and use it to forward model longer light curves of the kind that may be observed within the LSST main survey. Thanks to our effcient implementations, these simulations are able to cover for example over 1 Myr with a roughly weekly cadence. We envision that this framework will become highly valuable to prepare for, and best exploit, data from upcoming time domain surveys, such as for example LSST.

We present a multi wavelength analysis of 28 of the most luminous low-redshift narrow-line, ultra-hard X-ray selected active galactic nuclei (AGN) drawn from the 70 month Swift/BAT all-sky survey, with bolometric luminosities of log(L_bol/erg/s) > 45.25. The broad goal of our study is to determine whether these objects have any distinctive properties, potentially setting them aside from lower-luminosity obscured AGN in the local Universe. Our analysis relies on the first data release of the BAT AGN Spectroscopic Survey (BASS/DR1) and on dedicated observations with the VLT, Palomar, and Keck observatories. We find that the vast majority of our sources agree with commonly used AGN selection criteria which are based on emission line ratios and on mid-infrared colours. Our AGN are predominantly hosted in massive galaxies (9.8 < log(M_*/M_sun) < 11.7); based on visual inspection of archival optical images, they appear to be mostly ellipticals. Otherwise, they do not have distinctive properties. Their radio luminosities, determined from publicly available survey data, show a large spread of almost 4 orders of magnitude - much broader than what is found for lower X-ray luminosity obscured AGN in BASS. Moreover, our sample shows no preferred combination of black hole masses (M_BH) and/or Eddington ratio (lambda_Edd), covering 7.5 < log(M_BH/M_sun) < 10.3 and 0.01 < lambda_Edd < 1. Based on the distribution of our sources in the lambda_Edd-N_H plane, we conclude that our sample is consistent with a scenario where the amount of obscuring material along the line of sight is determined by radiation pressure exerted by the AGN on the dusty circumnuclear gas.

We investigate the abundance of Super-Massive Black Hole (SMBH) seeds in primordial galaxy halos. We explore the assumption that dark matter halos outgrowing a critical halo mass M_c have some probability p of having spawned a SMBH seed. Current observations of local, intermediate-mass galaxies constrain these parameters: For $M_c=10^{11}M_\odot$, all halos must be seeded, but when adopting smaller M_c masses the seeding can be much less efficient. The constraints also put lower limits on the number density of black holes in the local and high-redshift Universe. Reproducing z~6 quasar space densities depends on their typical halo mass, which can be constrained by counting nearby Lyman Break Galaxies and Lyman Alpha Emitters. For both observables, our simulations demonstrate that single-field predictions are too diverse to make definitive statements, in agreement with mixed claims in the literature. If quasars are not limited to the most massive host halos, they may represent a tiny fraction (~10^-5) of the SMBH population. Finally, we produce a wide range of predictions for gravitational events from SMBH mergers. We define a new diagnostic diagram for LISA to measure both SMBH space density and the typical delay between halo merger and black hole merger. While previous works have explored specific scenarios, our results hold independent of the seed mechanism, seed mass, obscuration, fueling methods and duty cycle.

Major galaxy mergers are thought to play an important part in fuelling the growth of supermassive black holes. However, observational support for this hypothesis is mixed, with some studies showing a correlation between merging galaxies and luminous quasars and others showing no such association. Recent observations have shown that a black hole is likely to become heavily obscured behind merger-driven gas and dust, even in the early stages of the merger, when the galaxies are well separated (5 to 40 kiloparsecs). Merger simulations further suggest that such obscuration and black-hole accretion peaks in the final merger stage, when the two galactic nuclei are closely separated (less than 3 kiloparsecs). Resolving this final stage requires a combination of high-spatial-resolution infrared imaging and high-sensitivity hard-X-ray observations to detect highly obscured sources. However, large numbers of obscured luminous accreting supermassive black holes have been recently detected nearby (distances below 250 megaparsecs) in X-ray observations. Here we report high-resolution infrared observations of hard-X-ray-selected black holes and the discovery of obscured nuclear mergers, the parent populations of supermassive-black-hole mergers. We find that obscured luminous black holes (bolometric luminosity higher than 2x10^44 ergs per second) show a significant (P<0.001) excess of late-stage nuclear mergers (17.6 per cent) compared to a sample of inactive galaxies with matching stellar masses and star formation rates (1.1 per cent), in agreement with theoretical predictions. Using hydrodynamic simulations, we confirm that the excess of nuclear mergers is indeed strongest for gas-rich major-merger hosts of obscured luminous black holes in this final stage.

The study of unobscured active galactic nuclei (AGN) and quasars depends on the reliable decomposition of the light from the AGN point source and the extended host galaxy light. The problem is typically approached using parametric fitting routines using separate models for the host galaxy and the point spread function (PSF). We present a new approach using a Generative Adversarial Network (GAN) trained on galaxy images. We test the method using Sloan Digital Sky Survey (SDSS) r-band images with artificial AGN point sources added which are then removed using the GAN and with parametric methods using GALFIT. When the AGN point source PS is more than twice as bright as the host galaxy, we find that our method, PSFGAN, can recover PS and host galaxy magnitudes with smaller systematic error and a lower average scatter ($49\%$). PSFGAN is more tolerant to poor knowledge of the PSF than parametric methods. Our tests show that PSFGAN is robust against a broadening in the PSF width of $\pm 50\%$ if it is trained on multiple PSF's. We demonstrate that while a matched training set does improve performance, we can still subtract point sources using a PSFGAN trained on non-astronomical images. While initial training is computationally expensive, evaluating PSFGAN on data is more than $40$ times faster than GALFIT fitting two components. Finally, PSFGAN it is more robust and easy to use than parametric methods as it requires no input parameters.

We present a framework to link and describe AGN variability on a wide range of timescales, from days to billions of years. In particular, we concentrate on the AGN variability features related to changes in black hole fuelling and accretion rate. In our framework, the variability features observed in different AGN at different timescales may be explained as realisations of the same underlying statistical properties. In this context, we propose a model to simulate the evolution of AGN light curves with time based on the probability density function (PDF) and power spectral density (PSD) of the Eddington ratio ($L/L_{\rm Edd}$) distribution. Motivated by general galaxy population properties, we propose that the PDF may be inspired by the $L/L_{\rm Edd}$ distribution function (ERDF), and that a single (or limited number of) ERDF+PSD set may explain all observed variability features. After outlining the framework and the model, we compile a set of variability measurements in terms of structure function (SF) and magnitude difference. We then combine the variability measurements on a SF plot ranging from days to Gyr. The proposed framework enables constraints on the underlying PSD and the ability to link AGN variability on different timescales, therefore providing new insights into AGN variability and black hole growth phenomena.

We present optical and near-IR Integral Field Unit (IFU) and ALMA band 6 observations of the nearby dual Active Galactic Nuclei (AGN) Mrk 463. At a distance of 210 Mpc, and a nuclear separation of $\sim$4 kpc, Mrk 463 is an excellent laboratory to study the gas dynamics, star formation processes and supermassive black hole (SMBH) accretion in a late-stage gas-rich major galaxy merger. The IFU observations reveal a complex morphology, including tidal tails, star-forming clumps, and emission line regions. The optical data, which map the full extent of the merger, show evidence for a biconical outflow and material outflowing at $>$600 km s$^{-1}$, both associated with the Mrk 463E nucleus, together with large scale gradients likely related to the ongoing galaxy merger. We further find an emission line region $\sim$11 kpc south of Mrk 463E that is consistent with being photoionized by an AGN. Compared to the current AGN luminosity, the energy budget of the cloud implies a luminosity drop in Mrk 463E by a factor 3-20 over the last 40,000 years. The ALMA observations of $^{12}$CO(2-1) and adjacent 1mm continuum reveal the presence of $\sim$10$^{9}$M$_\odot$ in molecular gas in the system. The molecular gas shows velocity gradients of $\sim$800 km/s and $\sim$400 km/s around the Mrk 463E and 463W nuclei, respectively. We conclude that in this system the infall of $\sim$100s $M_\odot$/yr of molecular gas is in rough balance with the removal of ionized gas by a biconical outflow being fueled by a relatively small, $<$0.01% of accretion onto each SMBH.

We present new Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the core of IC 2497, the galaxy associated with Hanny's Voorwerp. The combined fits of the Chandra (0.5-8 keV) and NuSTAR (3-24 keV) X-ray spectra, together with WISE mid-IR photometry, optical longslit spectroscopy and optical narrow-band imaging, suggest that the galaxy hosts a Compton-thick AGN ($N_{\rm H} \sim 2 \times 10^{24}$ cm$^{-2}$, current intrinsic luminosity $L_{\rm bol} \sim 2-5 \times 10^{44}$ erg s$^{-1}$) whose luminosity dropped by a factor of $\sim$50 within the last $\sim 100$ kyr. This corresponds to a change in Eddington ratio from $\rm \lambda_{Edd} \sim$ 0.35 to $\rm \lambda_{Edd} \sim$ 0.007. We argue that the AGN in IC 2497 should not be classified as a changing-look AGN, but rather we favour the interpretation where the AGN is undergoing a change in accretion state (from radiatively efficient to radiatively inefficient). In this scenario the observed drop in luminosity and Eddington ratio corresponds to the final stage of an AGN accretion phase. Our results are consistent with previous studies in the optical, X-ray and radio although the magnitude of the drop is lower than previously suggested. In addition, we discuss a possible analogy between X-ray binaries and an AGN.

We consider energy budgets and radiative history of 8 fading AGN, identified from mismatch between the ionizion of very extended (>10 kpc) ionized clouds and the luminosity of the nucleus viewed directly. All show significant fading on ~50,000-year timescales. We explore the use of minimum ionizing luminosity Q derived from photoionization balance in the brightest pixels in H-alpha at each projected radius. Tests using PG QSOs, and one target with detailed photoionization modeling, suggest that we can derive useful histories of individual AGN; the minimum ionizing luminosity is always an underestimate and subject to fine structure in the ionized material. These tests suggest that the underestimation from the upper envelope of Q values is roughly constant for a given object. These AGN show rapid drops and standstills; the common feature is a rapid drop in the last 20,000 years before our view of the nucleus. E-folding timescales are mostly thousands of years, with a few episodes as short as 400. In the limit of largely obscured AGN, we find additional evidence for fading, comparing lower limits from recombination balance and the maximum luminosities derived from from infrared fluxes. We compare these long-term light curves to simulations of AGN accretion; the strongest variations on these timespans are seen in models with strong and local feedback. Gemini integral-field optical spectroscopy shows a very limited role for outflows in these structures. While rings and loops of emission are common, their kinematic structure shows some to be in regular rotation. UGC 7342 exhibits local signatures of outflows <300 km/s, largely associated with very diffuse emission. Only in the Teacup AGN do we see outflow signatures of order 1000 km/s. Clouds around these fading AGN consist largely of tidal debris being externally illuminated but not displaced by AGN outflows. (Abridged)

In order to perform a more complete census of active galactic nuclei (AGN) in the local Universe, we investigate the use of the HeII~$\lambda 4685$ emission line diagnostic diagram by Shirazi & Brinchmann (2012) in addition to the standard methods based on other optical emission lines. The He II based diagnostics is more sensitive to AGN ionization in the presence of strong star formation than conventional line diagnostics. We survey a magnitude-limited sample of 63,915 galaxies from the Sloan Digital Sky Survey Data Release 7 at $0.02<z<0.05$ and use both the conventional BPT emission line diagnostic diagrams, as well as the HeII diagram to identify AGN. In this sample, 1,075 galaxies are selected as AGN using the BPT diagram, while an additional 234 galaxies are identified as AGN using the HeII diagnostic diagram, representing a 22% increase of AGN in the parent galaxy sample. We explore the host galaxy properties of these new HeII selected AGN candidates and find that they are most common in star-forming galaxies on the blue cloud and on the main sequence where ionization from star-formation is most likely to mask AGN emission in the BPT lines. We note in particular a high HeII AGN fraction in galaxies above the high-mass end of the main sequence where quenching is expected to occur. We use archival Chandra observations to confirm the AGN nature of candidates selected through HeII based diagnostic. Finally, we discuss how this technique can help inform galaxy/black hole co-evolution scenarios.

We present deep Chandra X-ray observations of the core of IC 2497, the galaxy associated with Hanny's Voorwerp and hosting a fading AGN. We find extended soft X-ray emission from hot gas around the low intrinsic luminosity (unobscured) AGN ($L_{\rm bol} \sim 10^{42}-10^{44}$ erg s$^{-1}$). The temperature structure in the hot gas suggests the presence of a bubble or cavity around the fading AGN ($\mbox{E$_{\rm bub}$} \sim 10^{54} - 10^{55}$ erg). A possible scenario is that this bubble is inflated by the fading AGN, which after changing accretion state is now in a kinetic mode. Other possibilities are that the bubble has been inflated by the past luminous quasar ($L_{\rm bol} \sim 10^{46}$ erg s$^{-1}$), or that the temperature gradient is an indication of a shock front from a superwind driven by the AGN. We discuss the possible scenarios and the implications for the AGN-host galaxy interaction, as well as an analogy between AGN and X-ray binaries lifecycles. We conclude that the AGN could inject mechanical energy into the host galaxy at the end of its lifecycle, and thus provide a source for mechanical feedback, in a similar way as observed for X-ray binaries.

We investigated AGN activity in low-mass galaxies, an important regime that can shed light onto BH formation and evolution, and their interaction with their host galaxies. We identified 336 AGN candidates from a parent sample of $\sim 48,000$ nearby low-mass galaxies ($M_{\rm \star} \leq 10^{9.5}M_\odot$, $z < 0.1$) in the SDSS. We selected the AGN using the classical BPT diagram, a similar optical emission line diagnostic based on the HeII$\lambda$4686 line, and mid-IR color cuts. Different criteria select host galaxies with different physical properties such as stellar mass and optical color, and only 3 out of 336 sources fulfill all three criteria. This could be in part due to selection biases. The resulting AGN fraction of $\sim 0.7 \%$ is at least one order of magnitude below the one estimated for more massive galaxies. At optical wavelengths, the HeII-based AGN selection appears to be more sensitive to AGN hosted in star-forming galaxies than the classical BPT diagram, at least in the low-mass regime. The archival X-ray and radio data available for some of the optically selected AGN candidates seem to confirm their AGN nature, but follow-up observations are needed to confirm the AGN nature of the rest of the sample, especially in the case of mid-IR selection. Our sample will be important for future follow-up studies aiming to understand the relation between BHs and host galaxies in the low-mass regime.

We present an observational constraint for the typical active galactic nucleus (AGN) phase lifetime. The argument is based on the time lag between an AGN central engine switching on and becoming visible in X-rays, and the time the AGN then requires to photoionize a large fraction of the host galaxy. Based on the typical light travel time across massive galaxies, and the observed fraction of X-ray selected AGN without AGN-photoionized narrow lines, we estimate that the AGN phase typically lasts $\sim10^{5}$ years. This lifetime is short compared to the total growth time of $10^{7}-10^{9}$ years estimated from e.g. the Soltan argument and implies that black holes grow via many such short bursts and that AGN therefore "flicker" on and off. We discuss some consequences of this flickering behavior for AGN feedback and the analogy of X-ray binaries and AGN lifecycles.