Globular clusters host a variety of lower-luminosity ($L_X<10^{35}$ erg s$^{-1}$) X-ray sources, including accreting neutron stars and black holes, millisecond pulsars, cataclysmic variables, and chromospherically active binaries. In this paper, we provide a comprehensive catalog of more than 1100 X-ray sources in 38 Galactic globular clusters observed by the Chandra X-ray Observatory's ACIS detector. The targets are selected to complement the MAVERIC survey's deep radio continuum maps of Galactic globular clusters. We perform photometry and spectral analysis for each source, determine a best-fit model, and assess the possibility of it being a foreground/background source based on its spectral properties and location in the cluster. We also provide basic assessments of variability. We discuss the distribution of X-ray binaries in globular clusters, their X-ray luminosity function, and carefully analyze systems with $L_X > 10^{33}$ erg s$^{-1}$. Among these moderately bright systems, we discover a new source in NGC 6539 that may be a candidate accreting stellar-mass black hole or a transitional millisecond pulsar. We show that quiescent neutron star LMXBs in globular clusters may spend ~2% of their lifetimes as transitional millisecond pulsars in their active ($L_X>10^{33}$ erg s$^{-1}$) state. Finally, we identify a substantial under-abundance of bright ($L_X>10^{33}$ erg s$^{-1}$) intermediate polars in globular clusters compared to the Galactic field, in contrast with the literature of the past two decades.

Supernova spectral time series contain a wealth of information about the progenitor and explosion process of these energetic events. The modeling of these data requires the exploration of very high dimensional posterior probabilities with expensive radiative transfer codes. Even modest parametrizations of supernovae contain more than ten parameters and a detailed exploration demands at least several million function evaluations. Physically realistic models require at least tens of CPU minutes per evaluation putting a detailed reconstruction of the explosion out of reach of traditional methodology. The advent of widely available libraries for the training of neural networks combined with their ability to approximate almost arbitrary functions with high precision allows for a new approach to this problem. Instead of evaluating the radiative transfer model itself, one can build a neural network proxy trained on the simulations but evaluating orders of magnitude faster. Such a framework is called an emulator or surrogate model. In this work, we present an emulator for the TARDIS supernova radiative transfer code applied to Type Ia supernova spectra. We show that we can train an emulator for this problem given a modest training set of a hundred thousand spectra (easily calculable on modern supercomputers). The results show an accuracy on the percent level (that are dominated by the Monte Carlo nature of TARDIS and not the emulator) with a speedup of several orders of magnitude. This method has a much broader set of applications and is not limited to the presented problem.

Resolved kinematical information, such as from molecular gas in star forming regions, is obtained from spectral line observations. However, these observations often contain multiple line-of-sight components, making estimates harder to obtain and interpret. We present a fully automatic method that determines the number of components along the line of sight, or the spectral multiplicity, through Bayesian model selection. The underlying open-source framework, based on nested sampling and conventional spectral line modeling, is tested using the large area ammonia maps of NGC 1333 in Perseus molecular cloud obtained by the Green Bank Ammonia Survey (GAS). Compared to classic approaches, the presented method constrains velocities and velocity dispersions in a larger area. In addition, we find that the velocity dispersion distribution among multiple components did not change substantially from that of a single fit component analysis of the GAS data. These results showcase the power and relative ease of the fitting and model selection method, which makes it a unique tool to extract maximum information from complex spectral data.

We look to provide a detailed description of the SPectroscopic IDentification of ERosita Sources (SPIDERS) survey, an SDSS-IV programme aimed at obtaining spectroscopic classification and redshift measurements for complete samples of sufficiently bright X-ray sources. We describe the SPIDERS X-ray Point Source Spectroscopic Catalogue, considering its store of 11,092 observed spectra drawn from a parent sample of 14,759 ROSAT and XMM sources over an area of 5,129 deg$^2$ covered in SDSS-IV by the eBOSS survey. This programme represents the largest systematic spectroscopic observation of an X-ray selected sample. A total of 10,970 (98.9\%) of the observed objects are classified and 10,849 (97.8\%) have secure redshifts. The majority of the spectra (10,070 objects) are active galactic nuclei (AGN), 522 are cluster galaxies, and 294 are stars. The observed AGN redshift distribution is in good agreement with simulations based on empirical models for AGN activation and duty cycle. Forming composite spectra of type 1 AGN as a function of the mass and accretion rate of their black holes reveals systematic differences in the H-beta emission line profiles. This study paves the way for systematic spectroscopic observations of sources that are potentially to be discovered in the upcoming eROSITA survey over a large section of the sky.

We present a statistical analysis of the optical properties of an X-ray selected Type 1 AGN sample, using high signal to-noise ratio ($S/N>20$) spectra of the counterparts of the ROSAT/2RXS sources in the footprint of the SDSS-IV/SPIDERS (Spectroscopic IDentification of eROSITA Sources) programme. The final sample contains 2100 sources. It significantly extends the redshift and luminosity ranges ($\rm z \sim 0.01 - 0.80$ and $\rm L_{0.1-2.4 \,keV}$ ~ $\rm 2.0 \times 10^{41}$ - $\rm 1.0 \times 10^{46} \, erg s^{-1}$ ) used so far in this kind of analysis. By means of a Principal Component Analysis, we derive Eigenvector (EV) 1 and 2 in an eleven dimensional optical and X-ray parameter space, which are consistent with previous results. The validity of the correlations of the Eddington ratio $\rm L/L_{Edd}$ with EV1 and the black hole mass with EV2 are strongly confirmed These results imply that $\rm L/L_{Edd}$ and black hole mass are related to the diversity of the optical properties of Type 1 AGN. Investigating the relation of the width and asymmetry of $\rm H\beta$ and the relative strength of the iron emission $\rm r_{FeII}$ , we show that our analysis supports the presence of a distinct kinematic region: the Very Broad Line Region. Furthermore, comparing sources with a red-asymmetric broad $\rm H\beta$ emission line to sources for which it is blue-asymmetric, we find an intriguing difference in the correlation of the FeII and the continuum emission strengths. We show that this contrasting behaviour is consistent with a flattened, stratified model of the Broad Line Region, in which the FeII emitting region is shielded from the central source.

We present results from a data challenge posed to the radial velocity (RV) community: namely, to quantify the Bayesian "evidence" for n={0,1,2,3} planets in a set of synthetically generated RV datasets containing a range of planet signals. Participating teams were provided the same likelihood function and set of priors to use in their analysis. They applied a variety of methods to estimate Z, the marginal likelihood for each n-planet model, including cross-validation, the Laplace approximation, importance sampling, and nested sampling. We found the dispersion in Z across different methods grew with increasing n-planet models: ~3 for 0-planets, ~10 for 1-planet, ~100-1000 for 2-planets, and >10,000 for 3-planets. Most internal estimates of uncertainty in Z for individual methods significantly underestimated the observed dispersion across all methods. Methods that adopted a Monte Carlo approach by comparing estimates from multiple runs yielded plausible uncertainties. Finally, two classes of numerical algorithms (those based on importance and nested samplers) arrived at similar conclusions regarding the ratio of Zs for n and (n+1)-planet models. One analytic method (the Laplace approximation) demonstrated comparable performance. We express both optimism and caution: we demonstrate that it is practical to perform rigorous Bayesian model comparison for <=3-planet models, yet robust planet discoveries require researchers to better understand the uncertainty in Z and its connections to model selection.

We present a unification model for a clumpy obscurer in active galactic nuclei (AGN) and investigate the properties of the resulting X-ray spectrum. Our model is constructed to reproduce the column density distribution of the AGN population and cloud eclipse events in terms of their angular sizes and frequency. We developed and release a generalised Monte Carlo X-ray radiative transfer code, XARS, to compute X-ray spectra of obscurer models. The geometry results in strong Compton scattering, causing soft photons to escape also along Compton-thick sight lines. This makes our model spectra very similar to the Brightman & Nandra TORUS model. However, only if we introduce an additional Compton-thick reflector near the corona, we achieve good fits to NuSTAR spectra. This additional component in our model can be interpreted as part of the dust-free broad-line region, an inner wall or rim, or a warped disk. It cannot be attributed to a simple disk because the reflector must simultaneously block the line of sight to the corona and reflect its radiation. We release our model as an Xspec table model and present corresponding CLUMPY infrared spectra, paving the way for self-consistent multi-wavelength analyses.

In the context of the upcoming SRG/eROSITA survey, we present an N-body simulation-based mock catalogue for X-ray selected AGN samples. The model reproduces the observed hard X-ray AGN luminosity function (XLF) and the soft X-ray logN-logS from redshift 0 to 6. The XLF is reproduced to within $\pm5\%$ and the logN-logS to within $\pm20\%$. We develop a joint X-ray -- optical extinction and classification model. We adopt a set of empirical spectral energy distributions to predict observed magnitudes in the UV, optical and NIR. With the latest eROSITA all sky survey sensitivity model, we create a high-fidelity full-sky mock catalogue of X-ray AGN. It predicts their distributions in right ascension, declination, redshift and fluxes. Using empirical medium resolution optical spectral templates and an exposure time calculator, we find that $1.1\times10^6$ ($4\times10^5$) fiber-hours are needed to follow-up spectroscopically from the ground the detected X-ray AGN with an optical magnitude $21<r<22.8$ ($22.8<r<25$) with a 4-m (8-m) class multi-object spectroscopic facility. We find that future clustering studies will measure the AGN bias to the percent level at redshift $z<1.2$ and should discriminate possible scenarios of galaxy-AGN co-evolution. We predict the accuracy to which the baryon acoustic oscillation standard ruler will be measured using X-ray AGN: better than 3\% for AGN between redshift 0.5 to 3 and better than 1\% using the Ly$\alpha$ forest of X-ray QSOs discovered between redshift 2 and 3. eROSITA will provide an outstanding set of targets for future galaxy evolution and cosmological studies.

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.

As matter accretes onto the central supermassive black holes in active galactic nuclei (AGN), X-rays are emitted. We present a population synthesis model that accounts for the summed X-ray emission from growing black holes; modulo the efficiency of converting mass to X-rays, this is effectively a record of the accreted mass. We need this population synthesis model to reproduce observed constraints from X-ray surveys: the X-ray number counts, the observed fraction of Compton-thick AGN [$\log$ (N$_{\rm H}$/cm$^{\rm -2}$) $>$ 24] and the spectrum of the cosmic X-ray background (CXB), after accounting for selection biases. Over the past decade, X-ray surveys by {\it XMM-Newton}, {\it Chandra}, \textit{NuSTAR} and \textit{Swift}-BAT have provided greatly improved observational constraints. We find that no existing X-ray luminosity function (XLF) consistently reproduces all these observations. We take the uncertainty in AGN spectra into account, and use a neural network to compute an XLF that fits all observed constraints, including observed Compton-thick number counts and fractions. This new population synthesis model suggests that, intrinsically, 50$\pm$9\% (56$\pm$7\%) of all AGN within z $\simeq$ 0.1 (1.0) are Compton-thick.

The short gamma-ray burst (GRB) 170817A was the first GRB associated with a gravitational-wave event. Due to the exceptionally low luminosity of the prompt $\gamma$-ray and the afterglow emission, the origin of both radiation components is highly debated. The most discussed models for the burst and the afterglow include a regular GRB jet seen off-axis and the emission from the cocoon encompassing a "choked" jet. Here, we report low radio-frequency observations at 610 and 1390~MHz obtained with the Giant Metrewave Radio Telescope (GMRT). Our observations span a range of $\sim7$ to $\sim152$ days after the burst. The afterglow started to emerge at these low frequencies about 60~days after the burst. The $1390$~MHz light curve barely evolved between 60 and 150 days, but its evolution is also marginally consistent with a $F_\nu\propto t^{0.8}$ rise seen in higher frequencies. We model the radio data and archival X-ray, optical and high-frequency radio data with models of top-hat and Gaussian structured GRB jets. We performed a Markov Chain Monte Carlo analysis of the structured-jet parameter space. Though highly degenerate, useful bounds on the posterior probability distributions can be obtained. Our bounds of the viewing angle are consistent with that inferred from the gravitational wave signal. We estimate the energy budget in prompt emission to be an order of magnitude lower than that in the afterglow blast-wave.

Methods. We analyzed broadband XMM-Newton and NuSTAR observations of NGC 300 ULX1, performing phase-averaged and phase-resolved spectroscopy. We compared two physically motivated models for the source spectrum: Non-thermal accretion column emission modeled by a power law with high-energy exponential roll-off (AC model) vs multicolor thermal emission from an optically thick accretion envelope plus a hard power-law tail (MCAE model). We combine the findings of our Bayesian analysis with qualitative physical considerations to evaluate the suitability of each model. Results. The low-energy part (<2 keV) of the source spectrum is dominated by non-pulsating, multicolor thermal emission. The (pulsating) high energy continuum is more ambiguous. If modelled with the AC model a residual structure is detected, that can be modeled using a broad Gaussian absorption line centered at ~12 keV. However, the same residuals can be successfully modeled using the MCAE model, without the need for the absorption-like feature. Model comparison, using the Bayesian approach strongly indicates that the MCAE model -- without the absorption line -- is the preferred model. Conclusions. All models considered strongly indicate the presence of an accretion disk truncated at a large distance from the central object, as has been recently suggested for a large fraction of both pulsating and non-pulsating ULXs. The hard, pulsed emission is not described by a smooth spectral continuum. If modelled by a broad Gaussian absorption line, the fit residuals can be interpreted as a cyclotron scattering feature (CRSF) compatible with a ~10^12 G magnetic field. However, the MCAE model describes the spectral and temporal characteristics of the source emission, without the need for an absorption feature and yields physically meaningful parameter values. Therefore strong doubts are cast on the presence of a CRSF in NGC 300 ULX1.

The data torrent unleashed by current and upcoming astronomical surveys demands scalable analysis methods. Many machine learning approaches scale well, but separating the instrument measurement from the physical effects of interest, dealing with variable errors, and deriving parameter uncertainties is often an after-thought. Classic forward-folding analyses with Markov Chain Monte Carlo or Nested Sampling enable parameter estimation and model comparison, even for complex and slow-to-evaluate physical models. However, these approaches require independent runs for each data set, implying an unfeasible number of model evaluations in the Big Data regime. Here I present a new algorithm, collaborative nested sampling, for deriving parameter probability distributions for each observation. Importantly, the number of physical model evaluations scales sub-linearly with the number of data sets, and no assumptions about homogeneous errors, Gaussianity, the form of the model or heterogeneity/completeness of the observations need to be made. Collaborative nested sampling has immediate application in speeding up analyses of large surveys, integral-field-unit observations, and Monte Carlo simulations.

We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS survey and $\sim$200 ksec in the remainder of the field. In this paper, we outline the survey's scientific goals, describe our observing strategy, and detail our data reduction and point source detection algorithms. Our analysis has resulted in a total of 868 band-merged point sources detected with a false-positive Poisson probability of $<1\times10^{-4}$. In addition, we present the results of an X-ray spectral analysis and provide best-fitting neutral hydrogen column densities, $N_{\rm H}$, as well as a sample of 51 Compton-thick active galactic nucleus candidates. Using this sample, we find the intrinsic Compton-thick fraction to be 30-35% over a wide range in redshift ($z=0.1-3$), suggesting the obscured fraction does not evolve very strongly with epoch. However, if we assume that the Compton-thick fraction is dependent on luminosity, as is seen for Compton-thin sources, then our results are consistent with a rise in the obscured fraction out to $z\sim3$. Finally, an examination of the host morphologies of our Compton-thick candidates shows a high fraction of morphological disturbances, in agreement with our previous results. All data products described in this paper are made available via a public website.

We report on the discovery of a dramatic X-ray spectral variability event observed in a $z\sim 1$ broad line type-1 QSO. The XMM-Newton spectrum from the year 2000 is characterized by an unobscured power-law spectrum with photon index of $\Gamma\sim 2$, a column density of $N_{\mathrm{H}}\sim 5\times 10^{20}\,\mathrm{cm^{-2}}$, and no prominent reflection component. Five years later, Chandra captured the source in a heavily-obscured, reflection-dominated state. The observed X-ray spectral variability could be caused by a Compton-thick cloud with $N_{\mathrm{H}}\sim 2\times 10^{24}\,\mathrm{cm^{-2}}$ eclipsing the direct emission of the hot corona, implying an extreme $N_{\mathrm{H}}$ variation never before observed in a type-1 QSO. An alternative scenario is a corona that switched off in between the observations. In addition, both explanations require a significant change of the X-ray luminosity prior to the obscuration or fading of the corona and/or a change of the relative geometry of the source/reflector system. Dramatic X-ray spectral variability of this kind could be quite common in type-1 QSOs, considering the relatively few datasets in which such an event could have been identified. Our analysis implies that there may be a population of type-1 QSOs which are Compton-thick in the X-rays when observed at any given time.

We report the detection of relativistically broadened iron K alpha emission in the X-ray spectra of AGN detected in the 4Ms CDF-S. Using the Bayesian X-ray analysis (BXA) package, we fit 199 hard band (2-7 keV) selected sources in the redshift range z=0.5--4 with three models: (i) an absorbed power-law, (ii) the first model plus a narrow reflection component, and (iii) the second model with an additional relativistic broadened reflection. The Bayesian evidence for the full sample of sources selects the model with the additional broad component as being 10^5 times more probable to describe the data better than the second model. For the two brightest sources in our sample, CID 190 (z=0.734) and CID 104 (z=0.543), BXA reveals the relativistic signatures in the individual spectra. We estimate the fraction of sources containing a broad component to be 54^{+35}_{-37}% (107/199 sources). Considering that the low signal-to-noise ratio of some spectra prevents the detection of the broad iron K alpha line, we infer an intrinsic fraction with broad emission of around two thirds. The detection of relativistic signatures in the X-ray spectra of these sources suggests that they are powered by a radiatively efficient accretion disk. Preliminary evidence is found that the spin of the black hole is high, with a maximally spinning Kerr BH model (a=1) providing a significantly better fit than a Schwarzschild model (a=0). Our analysis demonstrate the potential of X-ray spectroscopy to measure this key parameter in typical SMBH systems at the peak of BH growth.

Context: Redshifts are fundamental for our understanding of extragalactic X-ray sources. Ambiguous counterpart associations, expensive optical spectroscopy and/or multimission multiwavelength coverage to resolve degeneracies make estimation often difficult in practice. Aims: We attempt to constrain redshifts of obscured Active Galactic Nuclei (AGN) using only low-resolution X-ray spectra. Methods: Our XZ method fits AGN X-ray spectra with a moderately complex spectral model incorporating a corona, torus obscurer and warm mirror. Using the Bayesian X-ray Astronomy (BXA) package, we constrain redshift, column density, photon index and luminosity simultaneously. The redshift information primarily comes from absorption edges in Compton-thin AGN, and from the Fe K$\alpha$ fluorescent line in heavily obscured AGN. A new generic background fitting method allows us to extract more information from limited numbers of source counts. Results: We derive redshift constraints for 74/321 hard-band detected sources in the Chandra deep field South. Comparing with spectroscopic redshifts, we find an outlier fraction of 8%, indicating that our model assumptions are valid. For three Chandra deep fields, we release our XZ redshift estimates. Conclusions: The independent XZ estimate is easy to apply and effective for a large fraction of obscured AGN in todays deep surveys without the need for any additional data. Comparing to different redshift estimation methods, XZ can resolve degeneracies in photometric redshifts, help to detect potential association problems and confirm uncertain single-line spectroscopic redshifts. With high spectral resolution and large collecting area, this technique will be highly effective for Athena/WFI observations.

The basic unified model of active galactic nuclei (AGN) invokes an anisotropic obscuring structure, usually referred to as a torus, to explain AGN obscuration as an angle-dependent effect. We present a new grid of X-ray spectral templates based on radiative transfer calculations in neutral gas in an approximately toroidal geometry, appropriate for CCD-resolution X-ray spectra (FWHM > 130 eV). Fitting the templates to broadband X-ray spectra of AGN provides constraints on two important geometrical parameters of the gas distribution around the supermassive black hole: the average column density and the covering factor. Compared to the currently available spectral templates, our model is more flexible, and capable of providing constraints on the main torus parameters in a wider range of AGN. We demonstrate the application of this model using hard X-ray spectra from NuSTAR (3-79 keV) for four AGN covering a variety of classifications: 3C 390.3, NGC 2110, IC 5063, and NGC 7582. This small set of examples was chosen to illustrate the range of possible torus configurations, from disk-like to sphere-like geometries with column densities below, as well as above, the Compton-thick threshold. This diversity of torus properties challenges the simple assumption of a standard geometrically and optically thick toroidal structure commonly invoked in the basic form of the unified model of AGN. Finding broad consistency between the our constraints and those from infrared modeling, we discuss how the approach from the X-ray band complements similar measurements of AGN structures at other wavelengths.

The SUperluminous Supernova Host galaxIES (SUSHIES) survey aims to provide strong new constraints on the progenitors of superluminous supernovae (SLSNe) by understanding the relationship to their host galaxies. We present the photometric properties of 53 H-poor and 16 H-rich SLSN host galaxies out to $z\sim4$. We model their spectral energy distributions to derive physical properties, which we compare with other galaxy populations. At low redshift, H-poor SLSNe are preferentially found in very blue, low-mass galaxies with high average specific star-formation rates. As redshift increases, the host population follows the general evolution of star-forming galaxies towards more luminous galaxies. After accounting for secular evolution, we find evidence for differential evolution in galaxy mass, but not in the $B$-band and the far UV luminosity ($3\sigma$ confidence). Most remarkable is the scarcity of hosts with stellar masses above $10^{10}~M_\odot$ for both classes of SLSNe. In the case of H-poor SLSNe, we attribute this to a stifled production efficiency above $\sim0.4$ solar metallicity. However, we argue that, in addition to low metallicity, a short-lived stellar population is also required to regulate the SLSN production. H-rich SLSNe are found in a very diverse population of star-forming galaxies. Still, the scarcity of massive hosts suggests a stifled production efficiency above $\sim0.8$ solar metallicity. The large dispersion of the H-rich SLSNe host properties is in stark contrast to those of gamma-ray burst, regular core-collapse SN, and H-poor SLSNe host galaxies. We propose that multiple progenitor channels give rise to this sub-class.

Binary neutron-star mergers (BNSMs) are among the most readily detectable gravitational-wave (GW) sources with LIGO. They are also thought to produce short $\gamma$-ray bursts (SGRBs), and kilonovae that are powered by r-process nuclei. Detecting these phenomena simultaneously would provide an unprecedented view of the physics during and after the merger of two compact objects. Such a Rosetta Stone event was detected by LIGO/Virgo on 17 August 2017 at a distance of $\sim 44$ Mpc. We monitored the position of the BNSM with ALMA at 338.5 GHz and GMRT at 1.4 GHz, from 1.4 to 44 days after the merger. Our observations rule out any afterglow more luminous than $3\times 10^{26}~{\rm erg\,s}^{-1}\,{\rm Hz}^{-1}$ in these bands, probing $>$2--4 dex fainter than previous SGRB limits. We match these limits, in conjunction with public data announcing the appearance of X-ray and radio emission in the weeks after the GW event, to templates of off-axis afterglows. Our broadband modeling suggests that GW170817 was accompanied by a SGRB and that the GRB jet, powered by $E_{\rm AG,\,iso}\sim10^{50}$~erg, had a half-opening angle of $\sim20^\circ$, and was misaligned by $\sim41^\circ$ from our line of sight. The data are also consistent with a more collimated jet: $E_{\rm AG,\,iso}\sim10^{51}$~erg, $\theta_{1/2,\,\rm jet}\sim5^\circ$, $\theta_{\rm obs}\sim17^\circ$. This is the most conclusive detection of an off-axis GRB afterglow and the first associated with a BNSM-GW event to date. Assuming a uniform top-hat jet, we use the viewing angle estimates to infer the initial bulk Lorentz factor and true energy release of the burst.

We release the AllWISE counterparts and Gaia matches to 106,573 and 17,665 X-ray sources detected in the ROSAT 2RXS and XMMSL2 surveys with |b|>15. These are the brightest X-ray sources in the sky, but their position uncertainties and the sparse multi-wavelength coverage until now rendered the identification of their counterparts a demanding task with uncertain results. New all-sky multi-wavelength surveys of sufficient depth, like AllWISE and Gaia, and a new Bayesian statistics based algorithm, NWAY, allow us, for the first time, to provide reliable counterpart associations. NWAY extends previous distance and sky density based association methods and, using one or more priors (e.g., colors, magnitudes), weights the probability that sources from two or more catalogues are simultaneously associated on the basis of their observable characteristics. Here, counterparts have been determined using a WISE color-magnitude prior. A reference sample of 4524 XMM/Chandra and Swift X-ray sources demonstrates a reliability of ~ 94.7% (2RXS) and 97.4% (XMMSL2). Combining our results with Chandra-COSMOS data, we propose a new separation between stars and AGN in the X-ray/WISE flux-magnitude plane, valid over six orders of magnitude. We also release the NWAY code and its user manual. NWAY was extensively tested with XMM-COSMOS data. Using two different sets of priors, we find an agreement of 96% and 99% with published Likelihood Ratio methods. Our results were achieved faster and without any follow-up visual inspection. With the advent of deep and wide area surveys in X-rays (e.g. SRG/eROSITA, Athena/WFI) and radio (ASKAP/EMU, LOFAR, APERTIF, etc.) NWAY will provide a powerful and reliable counterpart identification tool.

Multiple channels have been proposed to produce Type Ia supernovae, with many scenarios suggesting that the exploding white dwarf accretes from a binary companion pre-explosion. In almost all cases, theory suggests that this companion will survive. However, no such companion has been unambiguously identified in ancient supernova remnants -- possibly falsifying the accretion scenario. Existing surveys, however, have only looked for stars as faint as $\approx 0.1 L_\odot$ and thus would have missed a surviving white dwarf companion. In this work, we present very deep DECAM imaging $(u, g, r, z)$ of the Type Ia supernova remnant SN 1006 specifically to search for a potential surviving white dwarf companion. We find no object within the inner third of the SN 1006 remnant that is consistent with a relatively young cooling white dwarf. We find that if there is a companion white dwarf, it must have formed long ago and cooled undisturbed for $> 10^8$ yr to be consistent with the redder objects in our sample. We conclude that our findings are consistent with the complete destruction of the secondary (such as in a merger) or an unexpectedly cool and thus very dim surviving companion white dwarf.

We present X-ray bolometric correction factors, $\kappa_{Bol}$ ($\equiv L_{Bol}/L_X$), for Compton-thick (CT) active galactic nuclei (AGN) with the aim of testing AGN torus models, probing orientation effects, and estimating the bolometric output of the most obscured AGN. We adopt bolometric luminosities, $L_{Bol}$, from literature infrared (IR) torus modeling and compile published intrinsic 2--10 keV X-ray luminosities, $L_{X}$, from X-ray torus modeling of NuSTAR data. Our sample consists of 10 local CT AGN where both of these estimates are available. We test for systematic differences in $\kappa_{Bol}$ values produced when using two widely used IR torus models and two widely used X-ray torus models, finding consistency within the uncertainties. We find that the mean $\kappa_{Bol}$ of our sample in the range $L_{Bol}\approx10^{42}-10^{45}$ erg/s is log$_{10}\kappa_{Bol}=1.44\pm0.12$ with an intrinsic scatter of $\sim0.2$ dex, and that our derived $\kappa_{Bol}$ values are consistent with previously established relationships between $\kappa_{Bol}$ and $L_{Bol}$ and $\kappa_{Bol}$ and Eddington ratio. We investigate if $\kappa_{Bol}$ is dependent on $N_H$ by comparing our results on CT AGN to published results on less-obscured AGN, finding no significant dependence. Since many of our sample are megamaser AGN, known to be viewed edge-on, and furthermore under the assumptions of AGN unification whereby unobscured AGN are viewed face-on, our result implies that the X-ray emitting corona is not strongly anisotropic. Finally, we present $\kappa_{Bol}$ values for CT AGN identified in X-ray surveys as a function of their observed $L_X$, where an estimate of their intrinsic $L_{X}$ is not available, and redshift, useful for estimating the bolometric output of the most obscured AGN across cosmic time.

The wide-area XMM-XXL X-ray survey is used to explore the fraction of obscured AGN at high accretion luminosities, $L_X (\rm 2-10 \, keV) > 10^{44} \, erg \,s ^{-1}$, and out to redshift $z\approx1.5$. The sample covers an area of about $\rm14\,deg^2$ and provides constraints on the space density of powerful AGN over a wide range of neutral hydrogen column densities extending beyond the Compton-thick limit, $\rm N_H\approx10^{24}\,cm^{-2}$. The fraction of obscured Compton-thin ($\rm N_H=10^{22}-10^{24}\,cm^{-2}$) AGN is estimated to be $\approx0.35$ for luminosities $L_X(\rm 2-10\,keV)>10^{44}\,erg\,s^{-1}$ independent of redshift. For less luminous sources the fraction of obscured Compton-thin AGN increases from $0.45\pm0.10$ at $z=0.25$ to $0.75\pm0.05$ at $z=1.25$. Studies that select AGN in the infrared via template fits to the observed Spectral Energy Distribution of extragalactic sources estimate space densities at high accretion luminosities consistent with the XMM-XXL constraints. There is no evidence for a large population of AGN (e.g. heavily obscured) identified in the infrared and missed at X-ray wavelengths. We further explore the mid-infrared colours of XMM-XXL AGN as a function of accretion luminosity, column density and redshift. The fraction of XMM-XXL sources that lie within the mid-infrared colour wedges defined in the literature to select AGN is primarily a function of redshift. This fraction increases from about 20-30% at z=0.25 to about 50-70% at $z=1.5$.

SPIDERS (SPectroscopic IDentification of eROSITA Sources) is an SDSS-IV survey running in parallel to the eBOSS cosmology project. SPIDERS will obtain optical spectroscopy for large numbers of X-ray-selected AGN and galaxy cluster members detected in wide area eROSITA, XMM-Newton and ROSAT surveys. We describe the methods used to choose spectroscopic targets for two sub-programmes of SPIDERS: X-ray selected AGN candidates detected in the ROSAT All Sky and the XMM-Newton Slew surveys. We have exploited a Bayesian cross-matching algorithm, guided by priors based on mid-IR colour-magnitude information from the WISE survey, to select the most probable optical counterpart to each X-ray detection. We empirically demonstrate the high fidelity of our counterpart selection method using a reference sample of bright well-localised X-ray sources collated from XMM-Newton, Chandra and Swift-XRT serendipitous catalogues, and also by examining blank-sky locations. We describe the down-selection steps which resulted in the final set of SPIDERS-AGN targets put forward for spectroscopy within the eBOSS/TDSS/SPIDERS survey, and present catalogues of these targets. We also present catalogues of ~12000 ROSAT and ~1500 XMM-Newton Slew survey sources which have existing optical spectroscopy from SDSS-DR12, including the results of our visual inspections. On completion of the SPIDERS program, we expect to have collected homogeneous spectroscopic redshift information over a footprint of ~7500 deg$^2$ for >85 percent of the ROSAT and XMM-Newton Slew survey sources having optical counterparts in the magnitude range 17<r<22.5, producing a large and highly complete sample of bright X-ray-selected AGN suitable for statistical studies of AGN evolution and clustering.

An important constraint for galaxy evolution models is how much gas resides in galaxies, in particular at the peak of star formation z=1-3. We attempt a novel approach by letting long-duration Gamma Ray Bursts (LGRBs) x-ray their host galaxies and deliver column densities to us. This requires a good understanding of the obscurer and biases introduced by incomplete follow-up observations. We analyse the X-ray afterglow of all 844 Swift LGRBs to date for their column density $N_H$. To derive the population properties we propagate all uncertainties in a consistent Bayesian methodology. The $N_H$ distribution covers the $10^{20-23}\mathrm{cm}^{-2}$ range and shows no evolutionary effect. Higher obscurations, e.g. Compton-thick columns, could have been detected but are not observed. The $N_H$ distribution is consistent with sources randomly populating a ellipsoidal gas cloud of major axis $N_H^\text{major}=10^{23}\mathrm{cm}^{-2}$ with 0.22 dex intrinsic scatter between objects. The unbiased SHOALS survey of afterglows and hosts allows us to constrain the relation between Spitzer-derived stellar masses and X-ray derived column densities $N_H$. We find a well-constrained powerlaw relation of $N_H=10^{21.7}\mathrm{cm}^{-2}\times\left(M_{\star}/10^{9.5}M_{\odot}\right)^{1/3}$, with 0.5 dex intrinsic scatter between objects. The Milky Way and the Magellanic clouds also follow this relation. From the geometry of the obscurer, its stellar mass dependence and comparison with local galaxies we conclude that LGRBs are primarily obscured by galaxy-scale gas. Ray tracing of simulated Illustris galaxies reveals a relation of the same normalisation, but a steeper stellar-mass dependence and mild redshift evolution. Our new approach provides valuable insight into the gas residing in high-redshift galaxies.

The "torus" obscurer of Active Galactic Nuclei (AGN) is poorly understood in terms of its density, substructure and physical mechanisms. Large X-ray surveys provide model boundary constraints, for both Compton-thin and Compton-thick levels of obscuration, as obscured fractions are mean covering factors $f_{\text{cov}}$. However, a major remaining uncertainty is host galaxy obscuration. In Paper I we discovered a relation of $N_H \propto M_{\star}^{1/3}$ for the obscuration of galaxy-scale gas. Here we apply this observational relation to the AGN population, and find that galaxy-scale gas is responsible for a luminosity-independent fraction of Compton-thin AGN, but does not produce Compton-thick columns. With the host galaxy obscuration understood, we present a model of the remaining, nuclear obscurer which is consistent with a range of observations. Our radiation-lifted torus model consists of a Compton-thick component ($f_{\text{cov}}\sim35\%$) and a Compton-thin component ($f_{\text{cov}}\sim40\%$), which depends on both black hole mass and luminosity. This provides a useful summary of observational constraints for torus modellers who attempt to reproduce this behaviour. It can also be employed as a sub-grid recipe in cosmological simulations which do not resolve the torus. We also investigate host-galaxy X-ray obscuration inside cosmological, hydro-dynamic simulations (EAGLE, Illustris). The obscuration from ray-traced galaxy gas can agree with observations, but is highly sensitive to the chosen feedback assumptions.

We present the Data Release 12 Quasar catalog (DR12Q) from the Baryon Oscillation Spectroscopic Survey (BOSS) of the SDSS-III. This catalog includes all SDSS-III/BOSS objects that were spectroscopically targeted as quasar candidates during the full survey and that are confirmed as quasars via visual inspection of the spectra, have luminosities Mi[z=2]<-20.5 (in a LCDM cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega _L=0.7), and either display at least one emission line with a full width at half maximum (FWHM)larger than 500 km/s or, if not, have interesting/complex absorption features. The catalog also includes previously known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 297,301 quasars detected over 9,376 square degrees with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with z>2.15 is about an order of magnitude greater than the number of z>2.15 quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 29,580 broad absorption line quasars and lists their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry together with some information on the optical morphology and the selection criteria. When available, the catalog also provides information on the optical variability of quasars using SDSS and PTF multi-epoch photometry. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra, covering the wavelength region 3,600-10,500A at a spectral resolution in the range 1,300<R<2,500, can be retrieved from the SDSS Catalog Archive Server.

In this paper we describe and publicly release a catalogue consisting of 8445 point-like X-ray sources detected in the XMM-XXL north survey. For the 2512 AGN which have reliable spectroscopy, we present the X-ray spectral fitting which has been computed with a Bayesian approach. We have also applied an X-ray spectral stacking method to different sub-samples. We confirm the well-known Iwasawa-Taniguchi effect in our luminosity-redshift sub-samples, and argue that such an effect is due to a decrease in the covering factor of a distant obscuring torus with increasing X-ray luminosity. By comparing the distribution of the reflection fraction, the ratio of the normalization of the reflected component to the direct radiation, we find that the low-luminosity, low-redshift sub-sample had systematically higher reflection fraction values than the high-redshift, high-luminosity one. On the other hand, no significant difference is found between samples having similar luminosity but different redshift, suggesting that the structure of the torus does not evolve strongly with redshift. Contrary to previous works, we do not find evidence for an increasing photon index at high Eddington ratio. This may be an indication that the structure of the accretion disc changes as the Eddington ratio approaches unity. Comparing our X-ray spectral analysis results with the optical spectral classification, we find that $\sim20$ per cent of optical type-1 AGN show an X-ray absorbing column density higher than $10^{21.5}\,\rm{cm^{-2}}$, and about 50 per cent of type-2 AGN have an X-ray absorbing column density less than $10^{21.5}\,\rm{cm^{-2}}$. We suggest that the excess X-ray absorption shown in the high-luminosity optical type-1 AGN can be due to small-scale dust-free gas within (or close to) the broad line region, while in the low-luminosity ones it can be due to a clumpy torus with a large covering factor.

The active galactic nuclei X-ray luminosity function traces actively accreting supermassive black holes and is essential for the study of the properties of the active galactic nuclei (AGN) population, black hole evolution, and galaxy-black hole coevolution. Up to now, the AGN luminosity function has been estimated several times in soft (0.5-2 keV) and hard X-rays (2-10 keV). AGN selection in these energy ranges often suffers from identification and redshift incompleteness and, at the same time, photoelectric absorption can obscure a significant amount of the X-ray radiation. We estimate the evolution of the luminosity function in the 5-10 keV band, where we effectively avoid the absorbed part of the spectrum, rendering absorption corrections unnecessary up to NH=10^23 cm^-2. Our dataset is a compilation of six wide, and deep fields: MAXI, HBSS, XMM-COSMOS, Lockman Hole, XMM-CDFS, AEGIS-XD, Chandra-COSMOS, and Chandra-CDFS. This extensive sample of ~1110 AGN (0.01<z<4.0, 41<log L_x<46) is 98% redshift complete with 68% spectroscopic redshifts. We use Bayesian analysis to select the best parametric model from simple pure luminosity and pure density evolution to more complicated luminosity and density evolution and luminosity-dependent density evolution. We estimate the model parameters that describe best our dataset separately for each survey and for the combined sample. We show that, according to Bayesian model selection, the preferred model for our dataset is the luminosity-dependent density evolution (LDDE). Our estimation of the AGN luminosity function does not require any assumption on the AGN absorption and is in good agreement with previous works in the 2-10 keV energy band based on X-ray hardness ratios to model the absorption in AGN up to redshift three. Our sample does not show evidence of a rapid decline of the AGN luminosity function up to redshift four. [abridged]

This paper presents a survey of X-ray selected active galactic nuclei (AGN) with optical spectroscopic follow-up in a $\sim 18\, \rm{deg^2}$ area of the equatorial XMM-XXL north field. A sample of 8445 point-like X-ray sources detected by XMM-Newton above a limiting flux of $F_{\rm 0.5-10\, keV} > 10^{-15} \rm\,erg\, cm^{-2}\, s^{-1}$ was matched to optical (SDSS) and infrared (WISE) counterparts. We followed up 3042 sources brighter than $r=22.5$ mag with the SDSS BOSS spectrograph. The spectra yielded a reliable redshift measurement for 2578 AGN in the redshift range $z=0.02-5.0$, with $0.5-2\rm\, keV$ luminosities ranging from $10^{39}-10^{46}\rm\,erg\,s^{-1}$. This is currently the largest published spectroscopic sample of X-ray selected AGN in a contiguous area. The BOSS spectra of AGN candidates show a bimodal distribution of optical line widths allowing a separation between broad- and narrow-emission line AGN. The former dominate our sample (70 per cent) due to the relatively bright X-ray flux limit and the optical BOSS magnitude limit. We classify the narrow emission line objects (22 per cent of full sample) using standard BPT diagnostics: the majority have line ratios indicating the dominant source of ionization is the AGN. A small number (8 per cent of full sample) exhibit the typical narrow line ratios of star-forming galaxies, or only have absorption lines in their spectra. We term the latter two classes "elusive'' AGN. We also compare X-ray, optical and infrared color AGN selections in this field. X-ray observations reveal, the largest number of AGN. The overlap between the selections, which is a strong function of the imaging depth in a given band, is also remarkably small. We show using spectral stacking that a large fraction of the X-ray AGN would not be selectable via optical or IR colours due to host galaxy contamination.

We combine deep X-ray survey data from the Chandra observatory and the wide-area/shallow XMM-XXL field to estimate the AGN X-ray luminosity function in the redshift range z=3-5. The sample consists of nearly 340 sources with either photometric (212) or spectroscopic (128) redshift in the above range. The combination of deep and shallow survey fields provides a luminosity baseline of three orders of magnitude, Lx(2-10keV)~1e43-1e46erg/s at z>3. We follow a Bayesian approach to determine the binned AGN space density and explore their evolution in a model-independent way. Our methodology accounts for Poisson errors in the determination of X-ray fluxes and uncertainties in photometric redshift estimates. We demonstrate that the latter is essential for unbiased measurement of space densities. We find that the AGN X-ray luminosity function evolves strongly between the redshift intervals z=3-4 and z=4-5. There is also suggestive evidence that the amplitude of this evolution is luminosity dependent. The space density of AGN with Lx<1e45erg/s drops by a factor of 5 between the redshift intervals above, while the evolution of brighter AGN appears to be milder. Comparison of our X-ray luminosity function with that of UV/optical selected QSOs at similar redshifts shows broad agreement at bright luminosities, Lx>1e45erg/s. The faint-end slope of UV/optical luminosity functions however, is steeper than for X-ray selected AGN. This implies that the type-I AGN fraction increases with decreasing luminosity at z>3, opposite to trends established at lower redshift. We also assess the significance of AGN in keeping the hydrogen ionised at high redshift. Our X-ray luminosity function yields ionising photon rate densities that are insufficient to keep the Universe ionised at redshift z>4. A source of uncertainty in this calculation is the escape fraction of UV photons for X-ray selected AGN.

We aim to constrain the evolution of AGN as a function of obscuration using an X-ray selected sample of $\sim2000$ AGN from a multi-tiered survey including the CDFS, AEGIS-XD, COSMOS and XMM-XXL fields. The spectra of individual X-ray sources are analysed using a Bayesian methodology with a physically realistic model to infer the posterior distribution of the hydrogen column density and intrinsic X-ray luminosity. We develop a novel non-parametric method which allows us to robustly infer the distribution of the AGN population in X-ray luminosity, redshift and obscuring column density, relying only on minimal smoothness assumptions. Our analysis properly incorporates uncertainties from low count spectra, photometric redshift measurements, association incompleteness and the limited sample size. We find that obscured AGN with $N_{H}>{\rm 10^{22}\, cm^{-2}}$ account for ${77}^{+4}_{-5}\%$ of the number density and luminosity density of the accretion SMBH population with $L_{{\rm X}}>10^{43}\text{ erg/s}$, averaged over cosmic time. Compton-thick AGN account for approximately half the number and luminosity density of the obscured population, and ${38}^{+8}_{-7}\%$ of the total. We also find evidence that the evolution is obscuration-dependent, with the strongest evolution around $N_{H}\thickapprox10^{23}\text{ cm}^{-2}$. We highlight this by measuring the obscured fraction in Compton-thin AGN, which increases towards $z\sim3$, where it is $25\%$ higher than the local value. In contrast the fraction of Compton-thick AGN is consistent with being constant at $\approx35\%$, independent of redshift and accretion luminosity. We discuss our findings in the context of existing models and conclude that the observed evolution is to first order a side-effect of anti-hierarchical growth.

We present photometric redshifts and associated probability distributions for all detected sources in the Extended Chandra Deep Field South (ECDFS). The work makes use of the most up-to-date data from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) and the Taiwan ECDFS Near-Infrared Survey (TENIS) in addition to other data. We also revisit multi-wavelength counterparts for published X-ray sources from the 4Ms-CDFS and 250ks-ECDFS surveys, finding reliable counterparts for 1207 out of 1259 sources ($\sim 96\%$). Data used for photometric redshifts include intermediate-band photometry deblended using the TFIT method, which is used for the first time in this work. Photometric redshifts for X-ray source counterparts are based on a new library of AGN/galaxy hybrid templates appropriate for the faint X-ray population in the CDFS. Photometric redshift accuracy for normal galaxies is 0.010 and for X-ray sources is 0.014, and outlier fractions are $4\%$ and $5.4\%$ respectively. The results within the CANDELS coverage area are even better as demonstrated both by spectroscopic comparison and by galaxy-pair statistics. Intermediate-band photometry, even if shallow, is valuable when combined with deep broad-band photometry. For best accuracy, templates must include emission lines.

AGN are known to have complex X-ray spectra that depend on both the properties of the accreting SMBH (e.g. mass, accretion rate) and the distribution of obscuring material in its vicinity ("torus"). Often however, simple and even unphysical models are adopted to represent the X-ray spectra of AGN. In the case of blank field surveys in particular, this should have an impact on e.g. the determination of the AGN luminosity function, the inferred accretion history of the Universe and also on our understanding of the relation between AGN and their host galaxies. We develop a Bayesian framework for model comparison and parameter estimation of X-ray spectra. We take into account uncertainties associated with X-ray data and photometric redshifts. We also demonstrate how Bayesian model comparison can be used to select among ten different physically motivated X-ray spectral models the one that provides a better representation of the observations. Despite the use of low-count spectra, our methodology is able to draw strong inferences on the geometry of the torus. For a sample of 350 AGN in the 4 Ms Chandra Deep Field South field, our analysis identifies four components needed to represent the diversity of the observed X-ray spectra: (abridged). Simpler models are ruled out with decisive evidence in favour of a geometrically extended structure with significant Compton scattering. Regarding the geometry of the obscurer, there is strong evidence against both a completely closed or entirely open toroidal geometry, in favour of an intermediate case. The additional Compton reflection required by data over that predicted by toroidal geometry models, may be a sign of a density gradient in the torus or reflection off the accretion disk. Finally, we release a catalogue with estimated parameters such as the accretion luminosity in the 2-10 keV band and the column density, $N_{H}$, of the obscurer.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Observations in the last decade have provided strong evidence that the growth of supermassive black holes at the centres of galaxies is among the most influential processes in galaxy evolution. Open questions that relate to our current understanding of black hole growth and its relation to the build-up of galaxies at redshifts z=1-4, when most black holes and stars we see in present-day galaxies were put in place, include: what is the nature of AGN feedback and whether it plays a significant role in the evolution of galaxies? what is the dominant population of accreting AGN at that critical epoch? is it dominated by obscured objects as required by many current observations and models? The Athena+ mission concept will provide the technological leap required for a breakthrough in our understanding of AGN and galaxy evolution at the heyday of the Universe. The high throughput of Athena+ will allow the systematic study of the incidence, nature and energetics of AGN feedback processes to z~4 via the identification and measurement of blue-shifted X-ray absorption lines. The excellent survey and spectral capabilities of the Athena+ Wide Field Imager will complete the census of black hole growth by yielding samples of up to 100 times larger than is currently possible of the most heavily obscured, including Compton thick, AGN to redshifts z~4. The demographics of this population relative to their hosts is fundamental for understanding how major black hole growth events relate to the build-up of galaxies.