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.

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.

We investigate the use of deep convolutional neural networks (deep CNNs) for automatic visual detection of galaxy mergers. Moreover, we investigate the use of transfer learning in conjunction with CNNs, by retraining networks first trained on pictures of everyday objects. We test the hypothesis that transfer learning is useful for improving classification performance for small training sets. This would make transfer learning useful for finding rare objects in astronomical imaging datasets. We find that these deep learning methods perform significantly better than current state-of-the-art merger detection methods based on nonparametric systems like CAS and GM$_{20}$. Our method is end-to-end and robust to image noise and distortions; it can be applied directly without image preprocessing. We also find that transfer learning can act as a regulariser in some cases, leading to better overall classification accuracy ($p = 0.02$). Transfer learning on our full training set leads to a lowered error rate from 0.038 $\pm$ 1 down to 0.032 $\pm$ 1, a relative improvement of 15%. Finally, we perform a basic sanity-check by creating a merger sample with our method, and comparing with an already existing, manually created merger catalogue in terms of colour-mass distribution and stellar mass function.

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 use a phenomenological model which connects the galaxy and AGN populations to investigate the process of AGN triggering through major galaxy mergers at z~0. The model uses stellar mass functions as input and allows the prediction of AGN luminosity functions based on assumed Eddington ratio distribution functions (ERDFs). We show that the number of AGN hosted by merger galaxies relative to the total number of AGN increases as a function of AGN luminosity. This is due to more massive galaxies being more likely to undergo a merger and does not require the assumption that mergers lead to higher Eddington ratios than secular processes. Our qualitative analysis also shows that to match the observations, the probability of a merger galaxy hosting an AGN and accreting at a given Eddington value has to be increased by a factor ~10 relative to the general AGN population. An additional significant increase of the fraction of high Eddington ratio AGN among merger host galaxies leads to inconsistency with the observed X-ray luminosity function. Physically our results imply that, compared to the general galaxy population, the AGN fraction among merger galaxies is ~10 times higher. On average, merger triggering does however not lead to significantly higher Eddington ratios.

The American Astronomical Society's WorldWide Telescope (WWT) project enables terabytes of astronomical images, data, and stories to be viewed and shared among researchers, exhibited in science museums, projected into full-dome immersive planetariums and virtual reality headsets, and taught in classrooms from middle school to college levels. We review the WWT ecosystem, how WWT has been used in the astronomical community, and comment on future directions.

The large majority of the accreting supermassive black holes in the Universe are obscured by large columns of gas and dust. The location and evolution of this obscuring material have been the subject of intense research in the past decades, and are still highly debated. A decrease in the covering factor of the circumnuclear material with increasing accretion rates has been found by studies carried out across the electromagnetic spectrum. The origin of this trend has been suggested to be driven either by the increase in the inner radius of the obscuring material with incident luminosity due to the sublimation of dust; by the gravitational potential of the black hole; by radiative feedback; or by the interplay between outflows and inflows. However, the lack of a large, unbiased and complete sample of accreting black holes, with reliable information on gas column density, luminosity and mass, has left the main physical mechanism regulating obscuration unclear. Using a systematic multi-wavelength survey of hard X-ray-selected black holes, here we show that radiation pressure on dusty gas is indeed the main physical mechanism regulating the distribution of the circumnuclear material. Our results imply that the bulk of the obscuring dust and gas in these objects is located within the sphere of influence of the black hole (i.e., a few to tens of parsecs), and that it can be swept away even at low radiative output rates. The main physical driver of the differences between obscured and unobscured accreting black holes is therefore their mass-normalized accretion rate.

We use stellar mass functions to study the properties and the significance of quenching through major galaxy mergers. In addition to SDSS DR7 and Galaxy Zoo 1 data, we use samples of visually selected major galaxy mergers and post merger galaxies. We determine the stellar mass functions of the stages that we would expect major merger quenched galaxies to pass through on their way from the blue cloud to the red sequence: 1: major merger, 2: post merger, 3: blue early type, 4: green early type and 5: red early type. Based on the similar mass function shapes we conclude that major mergers are likely to form an evolutionary sequence from star formation to quiescence via quenching. Relative to all blue galaxies, the major merger fraction increases as a function of stellar mass. Major merger quenching is inconsistent with the mass and environment quenching model. At z~0 major merger quenched galaxies are unlikely to constitute the majority of galaxies that transition the green valley. Furthermore, between z~0-0.5 major merger quenched galaxies account for 1-5% of all quenched galaxies at a given stellar mass. Major galaxy mergers are therefore not a significant quenching pathway, neither at z~0 nor within the last 5 Gyr. The majority of red galaxies must have been quenched through an alternative quenching mechanism which causes a slow blue to red evolution.

We use a phenomenological model to show that black hole growth in the local Universe (z < 0.1) can be described by two separate, mass independent Eddington ratio distribution functions (ERDFs). We assume that black holes can be divided into two independent groups: those with radiatively efficient accretion, primarily hosted by optically blue and green galaxies, and those with radiatively inefficient accretion, which are mainly found in red galaxies. With observed galaxy stellar mass functions as input, we show that the observed AGN luminosity functions can be reproduced by using mass independent, broken power law shaped ERDFs. We use the observed hard X-ray and 1.4 GHz radio luminosity functions to constrain the ERDF for radiatively efficient and inefficient AGN, respectively. We also test alternative ERDF shapes and mass dependent models. Our results are consistent with a mass independent AGN fraction and AGN hosts being randomly drawn from the galaxy population. We argue that the ERDF is not shaped by galaxy-scale effects, but by how efficiently material can be transported from the inner few parsecs to the accretion disc. Our results are incompatible with the simplest form of mass quenching where massive galaxies host higher accretion rate AGN. Furthermore, if reaching a certain Eddington ratio is a sufficient condition for maintenance mode, it can occur in all red galaxies, not just the most massive ones.

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 a comprehensive method for determining stellar mass functions, and apply it to samples in the local Universe. We combine the classical 1/Vmax approach with STY, a parametric maximum likelihood method and SWML, a non-parametric maximum likelihood technique. In the parametric approach, we are assuming that the stellar mass function can be modelled by either a single or a double Schechter function and we use a likelihood ratio test to determine which model provides a better fit to the data. We discuss how the stellar mass completeness as a function of z biases the three estimators and how it can affect, especially the low mass end of the stellar mass function. We apply our method to SDSS DR7 data in the redshift range from 0.02 to 0.06. We find that the entire galaxy sample is best described by a double Schechter function with the following parameters: $\log (M^{*}/M_\odot) = 10.79 \pm 0.01$, $\log (\Phi^{*}_1/\mathrm{h^3\ Mpc^{-3}}) = -3.31 \pm 0.20$, $\alpha_1 = -1.69 \pm 0.10$, $\log (\Phi^{*}_2/\mathrm{h^3\ Mpc^{-3}}) = -2.01 \pm 0.28$ and $\alpha_2 = -0.79 \pm 0.04$. We also use morphological classifications from Galaxy Zoo and halo mass, overdensity, central/satellite, colour and sSFR measurements to split the galaxy sample into over 130 subsamples. We determine and present the stellar mass functions and the best fit Schechter function parameters for each of these subsamples.

We investigate early black hole (BH) growth through the methodical search for $z\gtrsim5$ AGN in the $Chandra$ Deep Field South. We base our search on the $Chandra$ 4-Ms data with flux limits of $9.1\times\ 10^{-18}$ (soft, 0.5 - 2 keV) and $5.5\times\ 10^{-17}\ \mathrm{erg}\ \mathrm{s}^{-1}\ \mathrm{cm}^{-2}$ (hard, 2 - 8 keV). At $z\sim5$ this corresponds to luminosities as low as $\sim10^{42}$ ($\sim10^{43}$) $\mathrm{erg}\ \mathrm{s}^{-1}$ in the soft (hard) band and should allow us to detect Compton-thin AGN with $M_\mathrm{BH}>10^7 M_{\odot}$ and Eddington ratios > 0.1. Our field ($0.03~\mathrm{deg}^2$) contains over 600 $z\sim5$ Lyman Break Galaxies. Based on lower redshift relations we would expect $\sim20$ of them to host AGN. After combining the $Chandra$ data with GOODS/ACS, CANDELS/WFC3 and $Spitzer$/IRAC data, the sample consists of 58 high-redshift candidates. We run a photometric redshift code, stack the GOODS/ACS data, apply colour criteria and the Lyman Break Technique and use the X-ray Hardness Ratio. We combine our tests and using additional data find that all sources are most likely at low redshift. We also find five X-ray sources without a counterpart in the optical or infrared which might be spurious detections. We conclude that our field does not contain any convincing $z\gtrsim5$ AGN. Explanations for this result include a low BH occupation fraction, a low AGN fraction, short, super-Eddington growth modes, BH growth through BH-BH mergers or in optically faint galaxies. By searching for $z\gtrsim5$ AGN we are setting the foundation for constraining early BH growth and seed formation scenarios.