Search

We present observations from X-ray to mid-infrared wavelengths of the most energetic non-quasar transient ever observed, AT2021lwx. Our data show a single optical brightening by a factor $>100$ to a luminosity of $7\times10^{45}$ erg s$^{-1}$, and a total radiated energy of $1.5\times10^{53}$ erg, both greater than any known optical transient. The decline is smooth and exponential and the ultra-violet - optical spectral energy distribution resembles a black body with temperature $1.2\times10^4$ K. Tentative X-ray detections indicate a secondary mode of emission, while a delayed mid-infrared flare points to the presence of dust surrounding the transient. The spectra are similar to recently discovered optical flares in known active galactic nuclei but lack some characteristic features. The lack of emission for the previous seven years is inconsistent with the short-term, stochastic variability observed in quasars, while the extreme luminosity and long timescale of the transient disfavour the disruption of a single solar-mass star. The luminosity could be generated by the disruption of a much more massive star, but the likelihood of such an event occurring is small. A plausible scenario is the accretion of a giant molecular cloud by a dormant black hole of $10^8 - 10^9$ solar masses. AT2021lwx thus represents an extreme extension of the known scenarios of black hole accretion., 11 pages, 5 figures, Accepted for publication in MNRAS

Publisher Copyright: © 2023 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society. The bright blazar OJ 287 routinely parades high brightness bremsstrahlung flares, which are explained as being a result of a secondary supermassive black hole (SMBH) impacting the accretion disc of a more massive primary SMBH in a binary system. The accretion disc is not rigid but rather bends in a calculable way due to the tidal influence of the secondary. Next, we refer to this phenomenon as a variable disc level. We begin by showing that these flares occur at times predicted by a simple analytical formula, based on general relativity inspired modified Kepler equation, which explains impact flares since 1888. The 2022 impact flare, namely flare number 26, is rather peculiar as it breaks the typical pattern of two impact flares per 12-yr cycle. This is the third bremsstrahlung flare of the current cycle that follows the already observed 2015 and 2019 impact flares from OJ 287. It turns out that the arrival epoch of flare number 26 is sensitive to the level of primary SMBH's accretion disc relative to its mean level in our model. We incorporate these tidally induced changes in the level of the accretion disc to infer that the thermal flare should have occurred during 2022 July-August, when it was not possible to observe it from the Earth. Thereafter, we explore possible observational evidence for certain pre-flare activity by employing spectral and polarimetric data from our campaigns in 2004/05 and 2021/22. We point out theoretical and observational implications of two observed mini-flares during 2022 January-February.

X-ray luminosity functions (XLFs) of Active Galactic Nuclei (AGN) trace the growth and evolution of supermassive black hole populations across cosmic time, however, current XLF models are poorly constrained at redshifts of z>6, with a lack of spectroscopic constraints at these high redshifts. In this work we \redit{place limits} on the bright-end of the XLF at z=5.7-6.4 using high-redshift AGN identified within the Extragalactic Serendipitous Swift Survey (ExSeSS) catalogue. Within ExSeSS we find one serendipitously X-ray detected z>6 AGN, ATLAS J025.6821-33.4627, with an X-ray luminosity of $L_\mathrm{X}=8.47^{+3.40}_{-3.13}\times10^{44}\mathrm{erg.s^{-1}}$ and $z=6.31\pm0.03$ making it the highest redshift, spectroscopically confirmed, serendipitously X-ray detected quasar known to date. We also calculate an upper limit on the space density at higher luminosities where no additional sources are found, enabling us to place constraints on the shape of the XLF. Our results are consistent with the rapid decline in the space densities of high-luminosity AGN toward high redshift as predicted by extrapolations of existing parametric models of the XLF. We also find that our X-ray based measurements are consistent with estimates of the bolometric quasar luminosity function based on UV measurements at $z\gtrsim6$, although they require a large X-ray to bolometric correction factor at these high luminosities., 10pages, 6 figures. Resubmitted to MNRAS, following referee comments

We propose a solution to the problem of accretion disc sizes in active galactic nuclei being larger when measured by reverberation mapping than predicted by theory. Considering blackbody reprocessing on a disc with thickness profile $H(r)$, our solution invokes a steep rim or rippled structures irradiated by the central lamp-post. We model the continuum lags and the faint and bright disc spectral energy distribution (SED) in the best-studied case NGC 5548 (black hole mass $M = 7\times10^{7} M_\odot$, disc inclination $i=45^\circ$). With the lamp-post off, the observed disc SED requires a low accretion rate ($\dot{M} \sim 0.0014 M_\odot$/yr) and high prograde black hole spin ($a \sim 0.93$). Reprocessing on the thin disc gives time lags increasing with wavelength but 3 times smaller than observed. Introducing a steep $H(r)$ rim, or multiple crests, near $r = 5$ light days, reprocessing on their steep centre-facing slopes increases temperatures from $\sim1500$ K to $\sim6000$ K and this increases optical lags to match the lag data. Most of the disc surface maintains the cooler $T\propto r^{-3/4}$ temperature profile that matches the SED. The bright lamp-post may be powered by magnetic links tapping the black hole spin. The steep rim occurs near the sublimation radius for dust in the disc, as in the "failed disc wind model" for broad-line clouds. Lens-Thirring torques aligning the disc and black hole spin may also raise a warp and associated waves. In both scenarios, the small density scale height implied by the inferred value of $H(r)$ suggests possible marginal gravitational instability in the disc., 17 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Aims. We intend to use the impact of microlensing on the Fe III lambda lambda 2039 2113 emission line blend along with a measure of its gravitational redshift to estimate the mass of the quasar's central supermassive black hole (SMBH). Methods. We fit the Fe III feature in multiple spectroscopic observations between 2008 and 2016 of the gravitationally lensed quasar Q 0957 +561 with relatively high signal-to-noise ratios (at the adequate wavelength). Based on the statistics of microlensing magnifications, we used a Bayesian method to derive the size of its emitting region. Results. The Fe III lambda lambda 2039 2113 spectral feature appears systematically redshifted in all epochs of observation by a value of Delta lambda similar to 17 angstrom on average. We find clear di fferences in the shape of the Fe III line blend between images A and B. Measuring the strength of those magnitude di fferences, we conclude that this blend may arise from a region of half-light radius of R-1/2 similar to 15 ltdays, which is in good agreement with the accretion disk dimensions for this system. We obtain a mass for the central SMBH of M-BH = 1.5(-0.5)(+0.5) x 10(9) M-circle dot, consistent within uncertainties with previous mass estimates based on the virial theorem. The relatively small uncertainties in the mass determination (, Tel Aviv University, University of Haifa through a DFG grant HA3555-14/1, Spanish Government AYA2016- 79104-C3-1-P AYA2016-79104-C3-3-P, Center for Forestry Research & Experimentation (CIEF), European Commission Prometeo/2020/085, Spanish Government AYA2014-53506-P, European Commission Spanish Government AYA2014-53506-P, Junta de Andalucia FQM-108, Centro de Astrofisica de Valparaiso (CAV)

We study the demographics of z similar to 6 broad-line quasars in the black hole (BH) mass-luminosity plane using a sample of more than 100 quasars at 5.7 < z < 6.5. These quasars have well-quantified selection functions and nearly one-third of them also have virial BH masses estimated from near-IR spectroscopy. We use forward modelling of parametrized intrinsic distributions of BH masses and Eddington ratios, and account for the sample flux limits and measurement uncertainties of the BH masses and luminosities. We find significant differences between the intrinsic and observed distributions of the quantities due to measurement uncertainties and sample flux limits. There is also marginal evidence that the virial BH masses are susceptible to a positive luminosity-dependent bias (BH mass is overestimated when luminosity is above the average), and that the mean Eddington ratio increases with BH mass. Our models provide reliable constraints on the z similar to 6 BH mass function at M-BH > 10(8.5) M-circle dot, with a median 1 sigma uncertainty of similar to 0.5 dex in abundance. The intrinsic Eddington ratio distribution of M-BH > 10(8.5) M-circle dot quasars can be approximated by a mass-dependent Schechter model, with a broad peak around log (L-bol/L-Edd)similar to -0.9. We also find that, at 4.5 less than or similar to z less than or similar to 6, the number densities of more massive BHs tend to decline more rapidly with increasing redshift, contrary to the trend at 2.5 less than or similar to z less than or similar to 4.5 reported previously.

Pulsar timing arrays (PTAs) are searching for gravitational waves from supermassive black hole binaries (SMBHBs). Here we show how future PTAs could use a detection of gravitational waves from individually resolved SMBHB sources to produce a purely gravitational wave-based measurement of the Hubble constant. This is achieved by measuring two separate distances to the same source from the gravitational wave signal in the timing residual: the luminosity distance $D_L$ through frequency evolution effects, and the parallax distance $D_\mathrm{par}$ through wavefront curvature (Fresnel) effects. We present a generalized timing residual model including these effects in an expanding universe. Of these two distances, $D_\mathrm{par}$ is challenging to measure due to the pulsar distance wrapping problem, a degeneracy in the Earth-pulsar distance and gravitational wave source parameters that requires highly precise, sub-parsec level, pulsar distance measurements to overcome. However, in this paper we demonstrate that combining the knowledge of two SMBHB sources in the timing residual largely removes the wrapping cycle degeneracy. Two sources simultaneously calibrate the PTA by identifying the distances to the pulsars, which is useful in its own right, and allow recovery of the source luminosity and parallax distances which results in a measurement of the Hubble constant. We find that, with optimistic PTAs in the era of the Square Kilometer Array, two SMBHB sources within a few hundred Mpc could be used to measure the Hubble constant with a relative uncertainty on the order of 10 per cent., 23 pages, 11 figures. Published by the Monthly Notices of the Royal Astronomical Society. arXiv admin note: text overlap with arXiv:2109.07603

The gamma-ray emitting galaxy SBS 0846+513 has been classified as a Narrow-Line Seyfert 1 from its spectroscopy, and on that basis it was thought likely to have a small central black hole hosted in a spiral galaxy. But very few of the gamma-ray Narrow-Line Seyfert 1s have high-resolution imaging of their hosts, so it is unknown how those expectations hold up for the gamma-emitting class. We have observed this galaxy in the J-band with the Large Binocular Telescope's LUCI1 camera and the ARGOS adaptive optics system. We estimate its black hole mass to lie between $7.70 \leq \log \frac{\text{M}}{\text{M}_\odot} \leq 8.19$, using the correlation with bulge luminosity, or $7.96 \leq \log \frac{\text{M}}{\text{M}_\odot} \leq 8.16$ using the correlation with S\'{e}rsic index, putting its mass at the high end of the Narrow Line Seyfert 1 range. These estimates are independent of the Broad Line Region viewing geometry and avoid underestimates due to looking down the jet axis. Its host shows evidence of a bulge + disc structure, both from two-dimensional modeling and isophote shape, in keeping with the expectations. Mergers and interactions appear to be common among the gamma-ray Narrow-Line Seyfert 1s, and we see some circumstantial evidence for companion galaxies or disturbed features in the host., Comment: 10 pages, 6 figures. Submitted to MNRAS

Context. The strong asymmetry in the optical [O III] λ5007 emission line is one of the best signatures of active galactic nuclei (AGN) driven warm (∼104 K) ionized gas outflows on host galaxy scales. While large spectroscopic surveys such as the sloan digital sky survey (SDSS) have characterized the kinematics of [O III] for large samples of AGN, estimating the associated energetics requires spatially resolving these outflows with, for example, integral field unit (IFU) studies. Aims. As part of the Close AGN Reference Survey, we obtained spatially resolved IFU spectroscopy for a representative sample of 39 luminous type 1 AGN at 0.01 80% of the flux associated with a point-like source. We measured < 100 pc offsets in the spatial location of the outflow from the AGN nucleus using the spectro-astrometry technique for these sources. For the other 13 AGN, the [O III] wing emission is resolved and possibly extended on several kiloparsec scales. Conclusions. We conclude that [O III] wing emission can be compact or extended in an unbiased luminous AGN sample, where both cases are likely to appear. Electron density in the compact [O III] wing regions (median ne ∼ 1900 cm−3) is nearly a magnitude higher than in the extended ones (median ne ∼ 500 cm−3). The presence of spatially extended and compact [O III] wing emission is unrelated to the AGN bolometric luminosity and to inclination effects, which means other features such as time delays, or mechanical feedback (radio jets) may shape the ionized gas outflow properties. © ESO 2022., The work of MS was supported in part by the University of Manitoba Faculty of Science Graduate Fellowship (Cangene Award), and by the University of Manitoba Graduate Enhancement of Tri-Council Stipends (GETS) program. BH is greatful to financial support by the DFG grant GE625/17-1 and DLR grant 50OR1911. BH and ISP are additionally supported by the DAAD travel grant 57509925. TU acknowledges funding by the Competitive Fund of the Leibniz Association through grants SAW-2013-AIP-4, SAW-2015-AIP-2 and SAW-2016-IPHT-2. JS is supported by the international Gemini Observatory, a program of NSF’s NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation, on behalf of the Gemini partnership of Argentina, Brazil, Canada, Chile, the Republic of Korea, and the United States of America. MS, CO and SB acknowledge partial support from the Natural Sciences and Engineering Research Council (NSERC) of Canada. MG acknowledges partial support by NASA Chandra GO8-19104X/GO9-20114X and HST GO-15890.020-A. TR is supported by the Science and Technology Facilities Council (STFC) through grant ST/R504725/1. TAD acknowledges support from the UK Science and Technology Facilities Council through grant ST/S00033X/1. BAT was supported by the Harvard Future Faculty Leaders Postdoctoral Fellowship. Based on observations collected at the European Southern Observatory under ESO programme(s) 083.B-0801, 094.B-0345(A) and 095.B-0015(A). GRT acknowledges support from NASA through grant numbers HST-GO-15411 and HST-GO-15440 from the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. MPT acknowledges financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709) and through grant PGC2018-098915-B-C21 (MCI/AEI/FEDER, UE). The Science, Technology and Facilities Council is acknowledged by JN for support through the Consolidated Grant Cosmology and Astrophysics at Portsmouth, ST/S000550/1. Based on observations collected at the Centro Astronómico Hispano-Alemán (CAHA) at Calar Alto, operated jointly by Junta de Andalucía and Consejo Superior de Investigaciones Científicas (IAA-CSIC). This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018), SciPy (Virtanen et al. 2020), NumPy (van der Walt et al. 2011); and the plotting packages Matplotlib (Hunter 2007), CMasher (van der Velden 2020).

We present a study of 9242 spectroscopically-confirmed quasars with multi-epoch ugriz photometry from the SDSS Southern Survey. By fitting a separable linear model to each quasar's spectral variations, we decompose their five-band spectral energy distributions into variable (disc) and non-variable (host galaxy) components. In modelling the disc spectra, we include attenuation by dust on the line of sight through the host galaxy to its nucleus. We consider five commonly used attenuation laws, and find that the best description is by dust similar to that of the Small Magellanic Cloud, inferring a lack of carbonaceous grains from the relatively weak 2175AA absorption feature. We go on to construct a composite spectrum for the quasar variations spanning 700 to 8000AA. By varying the assumed power-law $L_{\nu}\propto\nu^\alpha$ spectral slope, we find a best-fit value $\alpha=0.71\pm0.02$, excluding at high confidence the canonical $L_{\nu}\propto\nu^{1/3}$ prediction for a steady-state accretion disc with a $T\propto r^{-3/4}$ temperature profile. The bluer spectral index of the observed quasar variations instead supports the model of Mummery & Balbus in which a steeper temperature profile, $T\propto r^{-7/8}$, develops as a result of finite magnetically-induced stress at the innermost stable circular orbit extracting energy and angular momentum from the black hole spin., Comment: 18 pages, 13 figures, accepted for publication in MNRAS

Improving our understanding of the nuclear properties of high-Eddington-ratio (λEdd) active galactic nuclei (AGN) is necessary since at this regime the radiation pressure is expected to affect the structure and efficiency of the accretion disc-corona system. This may cause departures from the typical nuclear properties of low-λEdd AGN, which have been largely studied so far. We present here the X-ray spectral analysis of 14 radio-quiet, λEdd ≳ 1 AGN at 0.4 ≤ z ≤ 0.75, observed with XMM-Newton. Optical/UV data from simultaneous Optical Monitor observations have also been considered. These quasars were selected to have relatively high values of black hole mass (MBH ∼ 108 − 8.5 M⊙) and bolometric luminosity (Lbol ∼ 1046 erg s−1) in order to complement previous studies of high-λEdd AGN at lower MBH and Lbol. We studied the relation between λEdd and other key X-ray spectral parameters, such as the photon index (Γ) of the power-law continuum, the X-ray bolometric correction (kbol, X), and the optical/UV-to-X-ray spectral index (αox). Our analysis reveals that, despite the homogeneous optical and supermassive black hole accretion properties, the X-ray properties of these high-λEdd AGN are quite heterogeneous. We indeed measured values of Γ between 1.3 and 2.5, at odds with the expectations based on previously reported Γ − λEdd relations, for which Γ ≥ 2 would be a ubiquitous hallmark of AGN with λEdd ∼ 1. Interestingly, we found that ∼30% of the sources are X-ray weak, with an X-ray emission about a factor of ∼10 − 80 fainter than that of typical AGN at similar UV luminosities. The X-ray weakness seems to be intrinsic and not due to the presence of absorption along the line of sight to the nucleus. This result may indicate that high-λEdd AGN commonly undergo periods of intrinsic X-ray weakness. Furthermore, results from follow-up monitoring with Swift of one of these X-ray weak sources suggest that these periods can last for several years. © ESO 2022., This work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. Part of this work is based on archival data, software and online services provided by the Space Science Data Center – ASI. This work has been partially supported by the ASI-INAF program I/004/11/4. ML acknowledges financial support from the PhD programme in Astronomy, Astrophysics and Space Science supported by MIUR (Ministero dell’Istruzione, dell’Università e della Ricerca). EP, LZ, FT, SB, AL and GV acknowledge financial support under ASI/INAF contract 2017-14-H.0. EP, SB, MB and GV acknowledge support from PRIN MIUR project “Black Hole winds and the Baryon Life Cycle of Galaxies: the stone-guest at the galaxy evolution supper”, contract #2017PH3WAT. AdO acknowledges financial support from the State Agency for Research of the Spanish MCIU through the project PID2019-106027GB-C41 and the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). RM acknowledges the financial support of INAF (Istituto Nazionale di Astrofisica), Osservatorio Astronomico di Roma, ASI (Agenzia Spaziale Italiana) under contract to INAF: ASI 2014-049-R.0 dedicated to SSDC. CR acknowledges support from the Fondecyt Iniciacion grant 11190831.

This is an Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Context. Active galactic nuclei (AGN) are thought to be intimately connected with their host galaxies through feeding and feedback processes. A strong coupling is predicted and supported by cosmological simulations of galaxy formation, but the details of the physical mechanisms are still observationally unconstrained. Aims. Galaxies are complex systems of stars and a multiphase interstellar medium (ISM). A spatially resolved multiwavelength survey is required to map the interaction of AGN with their host galaxies on different spatial scales and different phases of the ISM. The goal of the Close AGN Reference Survey (CARS) is to obtain the necessary spatially resolved multiwavelength observations for an unbiased sample of local unobscured luminous AGN. Methods. We present the overall CARS survey design and the associated wide-field optical integral-field unit (IFU) spectroscopy for all 41 CARS targets at z, BH is grateful for the financial support from the DFG grant GE625/17-1, DLR grant 50OR1911 and DAAD grant 57509925. The work of MS was supported in part by the University of Manitoba Faculty of Science Graduate Fellowship (Cangene Award), and by the University of Manitoba Graduate Enhancement of Tri-Council Stipends (GETS) program. JS is supported by the international Gemini Observatory, a program of NSF’s NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation, on behalf of the Gemini partnership of Argentina, Brazil, Canada, Chile, the Republic of Korea, and the United States of America. The Science, Technology and Facilities Council is acknowledged by JN for support through the Consolidated Grant Cosmology and Astrophysics at Portsmouth, ST/S000550/1. ISP acknowledges funding from the DLR grant 50OR2006. SMC and RM acknowledge support from the Australian Research Council through grant DP190102714. VNB gratefully acknowledges assistance from a National Science Foundation (NSF) Research at Undergraduate Institutions (RUI) grant AST-1909297. We note that findings and conclusions do not necessarily represent views of the NSF. TAD acknowledges support from STFC grant ST/S00033X/1. TR is supported by the Science and Technology Facilities Council (STFC) through grant ST/R504725/1. MK acknowledges support by DFG grant KR 3338/4-1. MG acknowledges partial support by NASA Chandra GO8-19104X/GO9-20114X and HST GO-15890.020-A grants. MPT acknowledges financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709) and through grant PGC2018-098915-B-C21 (MCI/AEI/FEDER, UE). GRT acknowledges support from NASA through grant numbers HST-GO-15411-.002-A, HST-GO-15440.002-A, and HST-GO-16173.001-A from the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. Based on observations collected at the European Southern Observatory under ESO program(s) 083.B-0801(A), 094.B-0345(A), 095.B-0015(A), 099.B-0242(B), and 099.B-0249(A). Based on observations collected at the Centro Astronómico Hispano-Alemán (CAHA) at Calar Alto, operated jointly by Junta de Andalucía and Consejo Superior de Investigaciones Científicas (IAA-CSIC). The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation. The national facility capability for SkyMapper has been funded through ARC LIEF grant LE130100104 from the Australian Research Council, awarded to the University of Sydney, the Australian National University, Swinburne University of Technology, the University of Queensland, the University of Western Australia, the University of Melbourne, Curtin University of Technology, Monash University and the Australian Astronomical Observatory. SkyMapper is owned and operated by The Australian National University’s Research School of Astronomy and Astrophysics. The survey data were processed and provided by the SkyMapper Team at ANU. The SkyMapper node of the All-Sky Virtual Observatory (ASVO) is hosted at the National Computational Infrastructure (NCI). Development and support the SkyMapper node of the ASVO has been funded in part by Astronomy Australia Limited (AAL) and the Australian Government through the Commonwealth’s Education Investment Fund (EIF) and National Collaborative Research Infrastructure Strategy (NCRIS), particularly the National eResearch Collaboration Tools and Resources (NeCTAR) and the Australian National Data Service Projects (ANDS).

[Abridged] We classify AT 2020zso as a TDE based on the blackbody evolution inferred from UV/optical photometric observations, and spectral line content and evolution. We identify transient, double-peaked Bowen (N III), He I, He II and Halpha emission lines. We model medium resolution optical spectroscopy of the He II (after careful deblending of the N III contribution) and Halpha lines during the rise, peak and early decline of the light curve using relativistic, elliptical accretion disk models. We find that the spectral evolution before peak can be explained by optical depth effects consistent with an outflowing, optically thick Eddington envelope. Around peak the envelope reaches its maximum extent (approximately 10^15 or 3000-6000 gravitational radii for an inferred black hole mass of 5-10 10^5) and becomes optically thin. The Halpha and He II emission lines at and after peak can be reproduced with a highly inclined (i=85+-5 degrees), highly elliptical (e=0.97+-0.01) and relatively compact (Rin = several 100 Rg and Rout = several 1000 Rg ) accretion disk. Overall, the line profiles suggest a highly elliptical geometry for the new accretion flow, consistent with theoretical expectations of newly formed TDE disks. We quantitatively confirm, for the first time, the high inclination nature of a Bowen (and X-ray dim) TDE, consistent with the unification picture of TDEs where the inclination largely determines the observational appearance. Rapid line profile variations rule out the binary SMBH hypothesis as the origin of the eccentricity; these results thus provide a direct link between a TDE in an AGN and the eccentric accretion disk. We illustrate for the first time how optical spectroscopy can be used to constrain the black hole spin, through (the lack of) disk precession signatures (changes in inferred inclination) - and rule out high black hole spin values (a < 0.8)., 20 pages, 15 figures and 5 tables. Accepted for publication in A&A. The spectra will be made publicly available through WISErep

GRAVITY+ is the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) with wide-separation fringe tracking, new adaptive optics, and laser guide stars on all four 8~m Unit Telescopes (UTs), for ever fainter, all-sky, high contrast, milliarcsecond interferometry. Here we present the design and first results of the first phase of GRAVITY+, called GRAVITY Wide. GRAVITY Wide combines the dual-beam capabilities of the VLTI and the GRAVITY instrument to increase the maximum separation between the science target and the reference star from 2 arcseconds with the 8 m UTs up to several 10 arcseconds, limited only by the Earth's turbulent atmosphere. This increases the sky-coverage of GRAVITY by two orders of magnitude, opening up milliarcsecond resolution observations of faint objects, and in particular the extragalactic sky. The first observations in 2019 - 2022 include first infrared interferometry of two redshift $z\sim2$ quasars, interferometric imaging on the binary system HD 105913A, and repeated observations of multiple star systems in the Orion Trapezium Cluster. We find the coherence loss between the science object and fringe-tracking reference star well described by the turbulence of the Earth's atmosphere. We confirm that the larger apertures of the UTs result in higher visibilities for a given separation due to larger overlap of the projected pupils on sky and give predictions for visibility loss as a function of separation to be used for future planning., Comment: 15 pages, 12 figures, 5 tables. Accepted by A&A

In early type galaxies, externally accreted gas is thought to be the main source of gas replenishment at late times. We use MUSE integral field spectroscopy data to study the active S0 galaxy NGC 5077, known to have disturbed dynamics, indicative of a past external interaction. We confirm the presence of a stellar kinematically distinct core with a diameter of 2.8 kpc, counter-rotating with respect to the main stellar body of the galaxy. We find that the counter-rotating core consists of an old stellar population, not significantly different from the rest of the galaxy. The ionised gas is strongly warped and extends out to 6.5 kpc in the polar direction and in a filamentary structure. The gas dynamics is complex, with significant changes in the position angle as a function of radius. The ionised gas line ratios are consistent with LINER excitation by the AGN both in the nucleus and at kiloparsec scales. We discover a nuclear outflow with projected velocity V ~ 400 km/s, consistent with a hollow outflow cone intersecting the plan of the sky. The properties of the misaligned gas match predictions from numerical simulations of misaligned gas infall after a gas-rich merger. The warp and change in the gas orientation as a function of radius are consistent with gas relaxation due to stellar torques, that are stronger at small radii where the gas aligns faster than in the outer regions, driving gas to the nucleus. The stellar and gas dynamics indicate that NGC 5077 has had at least two external interactions, one that resulted in the formation of the counter-rotating core followed by late time external gas accretion. NGC 5077 illustrates the importance of external interactions in the replenishment of the galaxy gas reservoir and the nuclear gas content available for black hole fuelling., Comment: Accepted for publication in A&A

The past year has seen numerous publications underlining the importance of a space mission to the ice giants in the upcoming decade. Proposed mission plans involve a $\sim$10 year cruise time to the ice giants. This cruise time can be utilized to search for low-frequency gravitational waves (GWs) by observing the Doppler shift caused by them in the Earth-spacecraft radio link. We calculate the sensitivity of prospective ice giant missions to GWs. Then, adopting a steady-state black hole binary population, we derive a conservative estimate for the detection rate of extreme mass ratio inspirals (EMRIs), supermassive- (SMBH) and stellar mass binary black hole (sBBH) mergers. We link the SMBH population to the fraction of quasars $f_\rm{bin}$ resulting from galaxy mergers that pair SMBHs to a binary. For a total of ten 40-day observations during the cruise of a single spacecraft, $\mathcal{O}(f_\rm{bin})\sim0.5$ detections of SMBH mergers are likely, if Allan deviation of Cassini-era noise is improved by $\sim 10^2$ in the $10^{-5}-10^{-3}$ Hz range. For EMRIs the number of detections lies between $\mathcal{O}(0.1) - \mathcal{O}(100)$. Furthermore, ice giant missions combined with the Laser Interferometer Space Antenna (LISA) would improve the localisation by an order of magnitude compared to LISA by itself., Accepted for publication in MNRAS: Letters

Accreting supermassive black holes (SMBHs), also known as active galactic nuclei (AGN), are generally surrounded by large amounts of gas and dust. This surrounding material reprocesses the primary X-ray emission produced close to the SMBH and gives rise to several components in the broadband X-ray spectra of AGN, including a power-law possibly associated with Thomson-scattered radiation. In this work, we study the properties of this scattered component for a sample of 386 hard-X-ray-selected, nearby ($z \sim 0.03$) obscured AGN from the 70-month Swift/BAT catalog. We investigate how the fraction of Thomson-scattered radiation correlates with different physical properties of AGN, such as line-of-sight column density, X-ray luminosity, black hole mass, and Eddington ratio. We find a significant negative correlation between the scattering fraction and the column density. Based on a large number of spectral simulations, we exclude the possibility that this anti-correlation is due to degeneracies between the parameters. The negative correlation also persists when considering different ranges of luminosity, black hole mass, and Eddington ratio. We discuss how this correlation might be either due to the angle dependence of the Thomson cross-section or to more obscured sources having a higher covering factor of the torus. We also find a positive correlation between the scattering fraction and the ratio of [OIII] $\lambda$5007 to X-ray luminosity. This result is consistent with previous studies and suggests that the Thomson-scattered component is associated with the narrow-line region., Comment: Accepted for publication in MNRAS, 18 pages, 9 figures, 2 tables, 1 equation

We publicly release the spectroscopic and photometric redshift catalog of the sources detected with Chandra in the field of the $z$=6.3 quasar SDSS J1030+0525. This is currently the fifth deepest X-ray field, and reaches a 0.5-2 keV flux limit $f_{\rm 0.5-2}$=6$\times$10$^{-17}$ erg s$^{-1}$ cm$^{-2}$. By using two independent methods, we measure a photometric redshift for 243 objects, while 123 (51%) sources also have a spectroscopic redshift, 110 of which coming from an INAF-Large Binocular Telescope (LBT) Strategic Program. We use the spectroscopic redshifts to determine the quality of the photometric ones, and find it in agreement with that of other X-ray surveys which used a similar number of photometric data-points. In particular, we measure a sample normalized median absolute deviation $\sigma_{NMAD}$=1.48||$z_{phot}$-$z_{spec}$||/(1+$z_{spec}$)=0.065. We use these new spectroscopic and photometric redshifts to study the properties of the Chandra J1030 field. We observe several peaks in our spectroscopic redshift distribution between $z$=0.15 and $z$=1.5, and find that the sources in each peak are often distributed across the whole Chandra field of view. This evidence confirms that X-ray selected AGN can efficiently track large-scale structures over physical scales of several Mpc. Finally, we computed the Chandra J1030 $z>$3 number counts: while the spectroscopic completeness at high-redshift of our sample is limited, our results point towards a potential source excess at $z\geq$4, which we plan to either confirm or reject in the near future with dedicated spectroscopic campaigns., Comment: 21 pages, 15 figures. Accepted for publication in Astronomy and Astrophysics

Low-mass compact galaxies (ultracompact dwarfs, UCDs, and compact ellipticals, cEs) populate the stellar size-mass plane between globular clusters and early-type galaxies. Known to be formed either in situ with an intrinsically low mass or resulting from the stripping of a more massive galaxy, the presence of a supermassive or an intermediate-mass black hole (BH) could help discriminate between these possible scenarios. With this aim, we have performed a multiwavelength search of active BH activity, i.e. active galactic nuclei (AGN), in a sample of 937 low-mass compact galaxies (580 UCDs and 357 cEs). This constitutes the largest study of AGN activity in these types of galaxies. Based on their X-ray luminosity, radio luminosity, and morphology, and/or optical emission line diagnostic diagrams, we find a total of 11 cEs that host an AGN. We also study for the first time the location of both low-mass compact galaxies (UCDs and cEs) and dwarf galaxies hosting AGN on the BH-galaxy scaling relations, finding that low-mass compact galaxies tend to be overmassive in the BH mass-stellar mass plane but not as much in the BH mass-stellar velocity dispersion correlation. This, together with available BH mass measurements for some of the low-mass compact galaxies, supports a stripping origin for the majority of these objects that would contribute to the scatter seen at the low-mass end of the BH-galaxy scaling relations. However, the differences are too large to be explained solely by this scatter, and thus our results suggest that a flattening at such low masses is also plausible, happening at a velocity dispersion of ~20-40 km s-1., We thank Soo-Chang Rey and Suk Kim for providing their list of cEs from Kim et al. (2020) and Aaron Romanowsky for insightful discussions. AFM has received financial support through the Postdoctoral Junior Leader Fellowship Programme from ‘La Caixa’ Banking Foundation (LCF/BQ/LI18/11630007). MM acknowledges support from the Beatriu de Pinos fellowship (2017-BP-00114) and the Ramon y Cajal fellowship (RYC2019-027670-I). AFM and MM both acknowledge support from the RAVET project PID2019-107427GB-C32 from the Spanish Ministry of Science and Innovation.

The quasar main sequence appears to be an incredibly powerful tool to organize the diversity in large samples of type-1 quasars but the most important physical parameters governing it are still unclear. Here, we investigate the origin of the broadening and of a defining feature of Population B sources: a strong redward asymmetry of the Balmer emission lines. We focus on a prototypical source, Fairall 9. Spectropolarimetric data of the Fairall 9 broad H β and H α profiles allowed for a view of the geometric and dynamical complexity of the line emitting regions. Measurements (1) provided evidence of rotational motion; (2) were helpful to test the presence of polar and equatorial scatterers, and their association with non-virial motions. However, we suggest that the polarization properties appear to be more consistent with a warped disc geometry induced by Lense–Thirring precession. © 2021 The Author(s)., BWJ and JMW acknowledge financial support from the National Natural Science Foundation of China (11833008 and 11991054), from the National Key Research and Development Program of China (2016YFA0400701), from the Key Research Program of Frontier Sciences, Chinese Academy of Sciences (CAS; QYZDJ-SSW-SLH007), and from the CAS Key Research Program (KJZD-EW-M06). VLA and ES thank the grant of Russian Science Foundation project number 20-12-00030 ‘Investigation of geometry and kinematics of ionized gas in active galactic nuclei by polarimetry methods’, which supported the spectropolarimetric data analysis. ĐS and LČP acknowledge funding provided by the Astronomical Observatory Belgrade (the contract 451-03-68/2020-14/200002) through the grants by the Ministry of Education, Science, and Technological Development of the Republic of Serbia. This work was supported by the F.R.S.–FNRS under grant PDR T.0116.21. AdO acknowledges financial support from the Spanish grants MCI PID2019-106027GB-C41 and the State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award for the IAA (SEV-2017-0709). Under the same award, PM acknowledges the Hypatia of Alexandria visiting grant.

Narrow-line Seyfert 1 (NLS1) galaxies are a class of active galactic nuclei (AGN) that, in some cases, can harbor powerful relativistic jets. One of them, PKS 2004-447, shows gamma-ray emission, and underwent its first recorded multifrequency flare in 2019. However, past studies revealed that in radio this source can be classified as a compact steep-spectrum source (CSS), suggesting that, unlike other gamma-ray sources, the relativistic jets of PKS 2004-447 have a large inclination with respect to the line of sight. We present here a set of spectroscopic observations of this object, aimed at carefully measuring its black hole mass and Eddington ratio, determining the properties of its emission lines, and characterizing its long term variability. We find that the black hole mass is $(1.5\pm0.2)\times10^7$ M$_\odot$, and the Eddington ratio is 0.08. Both values are within the typical range of NLS1s. The spectra also suggest that the 2019 flare was caused mainly by the relativistic jet, while the accretion disk played a minor role during the event. In conclusion, we confirm that PKS 2004-447 is one of the rare examples of gamma-ray emitting CSS/NLS1s hybrid, and that these two classes of objects are likely connected in the framework of AGN evolution., Comment: 14 pages, 8 figures, 4 tables, resubmitted to A&A after minor comments from the referee

arXiv:2106.10678v1, We use X-ray active galactic nuclei (AGN) observed by the Chandra X-ray Observatory within the 9.3 deg2 Boötes field of the NDWFS to study whether there is a correlation between X-ray luminosity (LX) and star formation rate (SFR) of the host galaxy, at 0.5, GM acknowledges support by the Agencia Estatal de Investigación, Unidad de Excelencia María de Maeztu, ref. MDM-2017-0765. MB acknowledges FONDECYT regular grant 1170618 The project has received funding from Excellence Initiative of Aix-Marseille University - AMIDEX, a French ‘Investissements d’Avenir’ programme. KM has been supported by the National Science Centre (UMO2018/30/E/ST9/00082).

In order to gauge the role that active galactic nuclei play in the evolution of galaxies via the effect of kinetic feedback in nearby QSO 2’s (z ∼ 0.3), we observed eight such objects with bolometric luminosities Lbol∼1046ergs−1 using Gemini Multi-Object Spectrograph-integral field units. The emission lines were fitted with at least two Gaussian curves, the broadest of which we attributed to gas kinetically disturbed by an outflow. We found that the maximum extent of the outflow ranges from ∼1 to 8 kpc, being ∼0.5±0.3 times the extent of the [OIII] ionized gas region. Our ‘default’ assumptions for the gas density (obtained from the [SII] doublet) and outflow velocities resulted in peak mass outflow rates of M˙defout∼ 3–30 M⊙yr−1 and outflow power of E˙defout∼1041–1043ergs−1⁠. The corresponding kinetic coupling efficiencies are εdeff=E˙defout/Lbol∼7×10−4–0.5 per cent, with the average efficiency being only 0.06 per cent (0.01 per cent median), implying little feedback powers from ionized gas outflows in the host galaxies. We investigated the effects of varying assumptions and calculations on M˙out and E˙out regarding the ionized gas densities, velocities, masses, and inclinations of the outflow relative to the plane of the sky, resulting in average uncertainties of 1 dex. In particular, we found that better indicators of the [OIII] emitting gas density than the default [SII] line ratio, such as the [ArIV] λλ4711,40 line ratio, result in almost an order of magnitude decrease in the εf. Peerreview

In order to gauge the role that active galactic nuclei play in the evolution of galaxies via the effect of kinetic feedback in nearby QSO 2’s (z ∼ 0.3), we observed eight such objects with bolometric luminosities Lbol∼1046ergs−1 using Gemini Multi-Object Spectrograph-integral field units. The emission lines were fitted with at least two Gaussian curves, the broadest of which we attributed to gas kinetically disturbed by an outflow. We found that the maximum extent of the outflow ranges from ∼1 to 8 kpc, being ∼0.5±0.3 times the extent of the [OIII] ionized gas region. Our ‘default’ assumptions for the gas density (obtained from the [SII] doublet) and outflow velocities resulted in peak mass outflow rates of M˙defout∼ 3–30 M⊙yr−1 and outflow power of E˙defout∼1041–1043ergs−1⁠. The corresponding kinetic coupling efficiencies are εdeff=E˙defout/Lbol∼7×10−4–0.5 per cent, with the average efficiency being only 0.06 per cent (0.01 per cent median), implying little feedback powers from ionized gas outflows in the host galaxies. We investigated the effects of varying assumptions and calculations on M˙out and E˙out regarding the ionized gas densities, velocities, masses, and inclinations of the outflow relative to the plane of the sky, resulting in average uncertainties of 1 dex. In particular, we found that better indicators of the [OIII] emitting gas density than the default [SII] line ratio, such as the [ArIV] λλ4711,40 line ratio, result in almost an order of magnitude decrease in the εf.

The most standard scenario for the evolution of massive galaxies across cosmic time assumes a correspondence based on the interplay between active galactic nuclei (AGN) feedback, which injects large amounts of energy into the host environment, and galaxy mergers, with their ability to trigger massive star formation events and accretion onto supermassive black holes. Interacting systems hosting AGN are useful laboratories for obtaining key insights into both phenomena. In this context, we present an analysis of the optical spectral properties of IRAS 20210+1121 (I20210), a merging system at z = 0.056. According to X-ray data, this object comprises two interacting galaxies, each hosting an obscured AGN. The optical spectra confirm the presence of AGN features in both galaxies. In particular, we are able to provide a Seyfert classification for I20210 North. The spectrum of I20120 South shows broad blueshifted components associated with the most intense emission lines that indicate the presence of an ionized outflow, for which we derive a maximum velocity of ∼2000 km s-1, an extension of ∼2 kpc, and a mass rate of ∼0.6 M⊙ yr-1. We also report the existence of an ionized nebular component with v ∼ 1000 km s-1 at ∼6.5 kpc southwards of I20210 South, which can be interpreted as disrupted gas ejected from the host galaxy by the action of the outflow. I20120 therefore exhibits a double obscured AGN, with one of them showing evidence of ongoing events for AGN-powered outflows. Future spatially resolved spectroscopy will allow for an accurate mapping of the gas kinematics in this AGN pair and evaluate the impact of the outflow on both the interstellar medium and the galaxy environment. © ESO 2021., GV, EP, CV, MB, CF and FF acknowledge support from PRIN MIUR project "Black Hole winds and the Baryon Life Cycle of Galaxies: the stone-guest at the galaxy evolution supper", contract #2017PH3WAT. GV also acknowledges financial support from Premiale 2015 MITic (PI: B. Garilli). CRA acknowledges financial support from the Spanish Ministry of Science, Innovation and Universities (MCIU) under grant with reference RYC-2014-15779, from the European Union's Horizon 2020 research and innovation programme under Marie Skodowska-Curie grant agreement No 860744 (BiD4BESt), from the State Research Agency (AEI-MCINN) of the Spanish MCIU under grants "Feeding and feedback in active galaxies" with reference PID2019-106027GB-C42, "Feeding, feedback and obscuration in active galaxies" with reference AYA2016-76682-C3-2-P, and "Quantifying the impact of quasar feedback on galaxy evolution (QSOFEED)" with reference EUR2020-112266. CRA also acknowledges support from the Consejeria de Economia, Conocimiento y Empleo del Gobierno de Canarias and the European Regional Development Fund (ERDF) under grant with reference ProID2020010105 and from IAC project P/301404, financed by the Ministry of Science and Innovation, through the State Budget and by the Canary Islands Department of Economy, Knowledge and Employment, through the Regional Budget of the Autonomous Community. Based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundacion Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

We thank the anonymous referee for the helpful comments and the constructive remarks on this manuscript. C.F. acknowledges the financial support the Tel Aviv University and University of Haifa through a DFG grant HA3555-14/1. E.M. and J.A.M are supported by the Spanish MINECO with the grants AYA2016- 79104-C3-1-P and AYA2016-79104-C3-3-P. J.A.M. is also supported from the Generalitat Valenciana project of excellence Prometeo/2020/085. J.J.V. is supported by the project AYA2017-84897-P financed by the Spanish Ministerio de Economia y Competividad and by the Fondo Europeo de Desarrollo Regional (FEDER), and by project FQM-108 financed by Junta de Andalucia. V.M. acknowledges partial support from Centro de Astrofisica de Valparaiso. Research by D.C. has been partially supported by an Israeli Science Foundation grant 2398/19., Aims. We aim to study the structure and kinematics of the broad line region (BLR) of a sample of 27 gravitationally lensed quasars with up to five different epochs of observation. This sample is composed of similar to 100 spectra from the literature plus 22 unpublished spectra of 11 systems. Methods. We measure the magnitude differences in the broad emission line (BEL) wings and statistically model the distribution of microlensing magnifications to determine a maximum likelihood estimate for the sizes of the C IV, C III], and Mg II emitting regions. Results. The BELs in lensed quasars are expected to be magnified differently owing to the different sizes of the regions from which they originate. Focusing on the most common BELs in our spectra (C IV, C III], and Mg II), we find that the low-ionization line Mg II is only weakly affected by microlensing. In contrast, the high-ionization line C IV shows strong microlensing in some cases, indicating that its emission region is more compact. Thus, the BEL profiles are deformed differently depending on the geometry and kinematics of the corresponding emitting region. We detect microlensing in either the blue or the red wing (or in both wings with different amplitudes) of C IV in more than 50% of the systems and find outstanding asymmetries in the wings of QSO 0957+561, SDSS J1004+4112, SDSS J1206+4332, and SDSS J1339+1310. This observation indicates that the BLR is, in general, not spherically symmetric and supports the existence of two regions in the BLR, one insensitive to microlensing and another that only shows up when it is magnified by microlensing. Disregarding the existence of these two regions, our estimate for Mg II, R-1/2 = 67.3(-15.7)(+3.8) root M/M-circle dot light-days, is in good agreement with previous results from smaller samples, while we obtain smaller sizes for the C III] and CIV lines, R-1/2 = 31.0(-4.0)(+1.9) root M/M-circle dot light-days and R-1/2 = 15.5(-3.9)(+0.8) root M/M-circle dot light-days, respectively., Tel Aviv University, University of Haifa through a DFG HA3555-14/1, Spanish MINECO AYA2016- 79104-C3-1-P AYA2016-79104-C3-3-P, Generalitat Valenciana European Commission, General Electric Prometeo/2020/085, Spanish Ministerio de Economia y Competividad AYA2017-84897-P, European Commission, Junta de Andalucia FQM-108, Centro de Astrofisica de Valparaiso, Israel Science Foundation 2398/19

Context. Sources at the brightest end of the quasi-stellar object (QSO) luminosity function, during the peak epoch in the history of star formation and black hole accretion (z ∼ 2-4, often referred to as "Cosmic noon") are privileged sites to study the cycle of feeding & feedback processes in massive galaxies. Aims. We aim to perform the first systematic study of cold gas properties in the most luminous QSOs, by characterising their host-galaxies and environment. These targets exhibit indeed widespread evidence of outflows at nuclear and galactic scales. Methods. We analyse ALMA, NOEMA and JVLA observations of the far-infrared continuum, CO and [CII] emission lines in eight QSOs (bolometric luminosity LBol ≳ 3 × 1047 erg s-1) from the WISE-SDSS selected hyper-luminous (WISSH) QSOs sample at z ∼ 2.4-4.7. Results. We report a 100% emission line detection rate and a 80% detection rate in continuum emission, and we find CO emission to be consistent with the steepest CO ladders observed so far. Sub-millimetre data reveal presence of (one or more) bright companion galaxies around ∼80% of WISSH QSOs, at projected distances of ∼6-130 kpc. We observe a variety of sizes for the molecular gas reservoirs (∼1.7-10 kpc), mostly associated with rotating disks with disturbed kinematics. WISSH QSOs typically show lower CO luminosity and higher star formation efficiency than infrared matched, z ∼ 0-3 main-sequence galaxies, implying that, given the observed SFR ∼170-1100 M⊙ yr-1, molecular gas is converted into stars in ≲ 50 Myr. Most targets show extreme dynamical to black-hole mass ratios Mdyn/MBH ∼ 3-10, two orders of magnitude smaller than local relations. The molecular gas fraction in the host-galaxies of WISSH is lower by a factor of ∼10-100 than in star forming galaxies with similar M∗. Conclusions. Our analysis reveals that hyper-luminous QSOs at Cosmic noon undergo an intense growth phase of both the central super-massive black hole and of the host-galaxy. These systems pinpoint the high-density sites where giant galaxies assemble, where we show that mergers play a major role in the build-up of the final host-galaxy mass. We suggest that the observed low molecular gas fraction and short depletion timescale are due to AGN feedback, whose presence is indicated by fast AGN-driven ionised outflows in all our targets. © ESO 2020., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00417.S, ADS/JAO.ALMA#2019.1.01070.S, ADS/JAO.ALMA#2015.1.01602.S, ADS/JAO.ALMA#2016.1.00718.S, ADS/JAO. ALMA#2016.1.01515.S, ADS/JAO.ALMA#2018.1.01806.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work is based on observations carried out under project number S17BW and W17DT with the IRAM NOEMA Interferometer. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). We are grateful to the anonymous referee for useful feedback which helped us to improve the paper. We thank Dr. Y. Wang for providing us with SUBARU images of J0209-0005 and Dr. P. Santini for useful discussion about SFE in high-z galaxies. MBi, MBr, CF, FF, AM, EP acknowledge support from PRIN MIUR project "Black Hole winds and the Baryon Life Cycle of Galaxies: the stone-guest at the galaxy evolution supper", contract #2017PH3WAT. MBi, CF and FF acknowledge support from INAF under PRIN SKA/CTA FORECaST and PRIN MAINSTREAM 2018 "Black hole winds and the baryon cycle". EP, LZ and MBi acknowledge financial support under ASI-INAF contract 2017-14-H.0. GB acknowledges financial support under the INTEGRAL ASI-INAF agreement 2019-35-HH.0 MPT acknowledges financial support from the State Agency for Research of the Spanish MCIU through the Center of Excellence Severo Ochoa award to the Instituto de Astrofisica de Andalucia (SEV-2017-0709) and through grant PGC2018-098915-B-C21 (MCI/AEI/FEDER, UE). RM acknowledges ERC Advanced Grant 695671 "QUENCH" and support by the Science and Technology Facilities Council (STFC). GVe acknowledges support from CONICYT Basal-CATA AFB-170002 and FONDECYT Postdoctorado 3200802 grants. CC acknowledges support from the Royal Society. MP is supported by the Programa Atraccion de Talento de la Comunidad de Madrid via grant 2018-T2/TIC-11715. GVi acknowledges financial support from Premiale 2015 MITiC (PI B. Garilli).

Results from a few decades of reverberation mapping (RM) studies have revealed a correlation between the radius of the broad-line emitting region (BLR) and the continuum luminosity of active galactic nuclei. This "radius-luminosity" relation enables survey-scale black-hole mass estimates across cosmic time, using relatively inexpensive single-epoch spectroscopy, rather than intensive RM time monitoring. However, recent results from newer reverberation mapping campaigns challenge this widely used paradigm, reporting quasar BLR sizes that differ significantly from the previously established radius-luminosity relation. Using simulations of the radius--luminosity relation with the observational parameters of SDSS-RM, we find that this difference is not likely due to observational biases. Instead, it appears that previous RM samples were biased to a subset of quasar properties, and the broader parameter space occupied by the SDSS-RM quasar sample has a genuinely wider range of BLR sizes. We examine the correlation between the deviations from the radius-luminosity relation and several quasar parameters; the most significant correlations indicate that the deviations depend on UV/optical SED and the relative amount of ionizing radiation. Our results indicate that single-epoch black-hole mass estimates that do not account for the diversity of quasars in the radius-luminosity relation could be overestimated by an average of ~0.3 dex., Comment: 11 pages, 12 figures. Submitted to ApJ. Current version includes revisions after referee report

arXiv:2004.04167v3, The leading contenders for the seeds of the first quasars are direct-collapse black holes (DCBHs) formed during catastrophic baryon collapse in atomically cooled halos at z ~ 20. The discovery of the Lyα emitter CR7 at z = 6.6 was initially held to be the first detection of a DCBH, although this interpretation has since been challenged on the grounds of Spitzer IRAC and Very Large Telescope X-Shooter data. Here we determine if radio flux from a DCBH in CR7 could be detected and discriminated from competing sources of radio emission in the halo such as young supernovae and H ii regions. We find that a DCBH would emit a flux of 10–200 nJy at 1.0 GHz, far greater than the sub-nJy signal expected for young supernovae but on par with continuum emission from star-forming regions. However, radio emission from a DCBH in CR7 could be distinguished from free–free emission from H ii regions by its spectral evolution with frequency and could be detected by the Square Kilometre Array in the coming decade., D. J. W. was supported by STFC New Applicant Grant ST/P000509/1 and the Ida Pfeiffer Professorship at the Institute of Astrophysics at the University of Vienna. M. M. acknowledges support from the Beatriu de Pinos fellowship (2017-BP-00114). A. M. acknowledges support from the UK Science and Technology Facilities Council Consolidated Grant ST/R000972/1. T. H. was supported by JSPS KAKENHI Grant Number 17F17320. M. L. acknowledges funding from UAEU via UPAR grant No. 31S372.

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$., 24+2 pages, 18+1 figures, published in Astronomy & Astrophysics

Following the recent discovery of X-ray quasi-periodic eruptions (QPEs) coming from the nucleus of the galaxy GSN 069, here we report on the detection of QPEs in the active galaxy named RX J1301.9+2747. QPEs are rapid and recurrent increases of the X-ray count-rate by more than one order of magnitude with respect to a stable quiescent level. During a XMM-Newton observation lasting 48 ks that was performed on 30 and 31 May 2019, three strong QPEs lasting about half an hour each were detected in the light curves of RX J1301.9+2747. The first two QPEs are separated by a longer recurrence time (about 20 ks) compared to the second and third (about 13 ks). This pattern is consistent with the alternating long-short recurrence times of the GSN 069 QPEs, although the difference between the consecutive recurrence times is significantly smaller in GSN 069. Longer X-ray observations will better clarify the temporal pattern of the QPEs in RX J1301.9+2747 and will allow a detailed comparison with GSN 069 to be performed. The X-ray spectral properties of QPEs in the two sources are remarkably similar, with QPEs representing fast transitions from a relatively cold and likely disk-dominated state to a state that is characterized by a warmer emission similar to the so-called soft X-ray excess, a component that is almost ubiquitously seen in the X-ray spectra of unobscured, radiatively efficient active galaxies. Previous X-ray observations of RX J1301.9+2747 in 2000 and 2009 strongly suggest that QPEs have been present for at least the past 18.5 years. The detection of QPEs from a second galactic nucleus after GSN 069 rules out contamination by a Galactic source in both cases, such that QPEs ought to be considered a novel extragalactic phenomenon associated with accreting supermassive black holes., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

Context. The quasar class of extreme Population A (xA) (also known as super-Eddington accreting massive black holes, SEAMBHs) has been hailed as potential distance indicators for cosmology. Aims. The aim of this paper is to define tight criteria for their proper identification, starting from the main selection criterion RFeII > 1, and to identify potential intruders that do not meet the selection criteria, which nonetheless have been selected as xA because of the coarseness of automatic searches. The inclusion of the spurious xA sources may dramatically increase the dispersion in the Hubble diagram of quasars, which were obtained from virial luminosity estimates. Methods. We studied a sample of 32 low-z quasars that were originally selected from the seventh data release of the Sloan Digital Sky Survey as xA or SEAMBHs, which have been proved to be almost certainly misclassified sources. All of them show moderate to strong Fe II emission and the large majority show strong absorption features in their spectra which are typical of fairly evolved stellar populations. We performed a simultaneous fit of a host galaxy spectrum, active galactic nucleus (AGN) continuum, FeII template, and emission lines to spectra, using the fitting technique based on ULySS, the full spectrum fitting package. We derived the main accretion parameters (i.e., luminosity, black hole mass, and Eddington ratio) and investigate the relation between host galaxy properties and AGN. Results. For sources in our sample with spectral types that correspond to a relatively low Eddington ratio, we find an overall consistency between HβNC, [O III]λλ4959,5007 line shifts, and the mean stellar velocity obtained from the host galaxy fit (within|60| km s-1). Only one source in our sample qualifies as a xA source. Conclusions. The correct classification of spectra that were contaminated by heavy absorption requires careful determination of the host galaxy spectrum. Contamination and misclassification are not usual in the identification of the xAs, nor at low z or at high z. We find a high fraction of host galaxy spectrum; in half of the sample this is even higher than 40%. When absorption lines are prominent, and the fraction of the host galaxy is high, host galaxy spectrum mimics FeII, which may result in a mistaken identification of FeII spectral features. We have identified several stellar absorption lines that, along with the continuum shape, may lead to an overestimate of RFeII, and therefore to the misclassification of sources as xA sources. © ESO 2020., This research is part of the projects 176001 "Astrophysical spectroscopy of extragalactic objects" and 176003 "Gravitation and the large scale structure of the Universe", funded by Ministry of Education, Science and Technological Development of the Republic of Serbia. PM and MDO acknowledge funding from the INAF PRIN-SKA 2017 program 1.05.01.88.04. PM also acknowledges the Programa de Estancias de Investigacion (PREI) No. DGAP/DFA/2192/2018 of UNAM. A.d.O. acknowledges financial support from the Spanish Ministry of Economy and Competitiveness through grant AYA2016-76682-C3-1-P and from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709). ML. M. A. acknowledges financial support of National Science Centre, Poland, grant No. 2017/26/A/ST9/00756 (Maestro 9). DD acknowledges support from grants PAPIIT, UNAM 113719, Mexico D.F. 04510, Mexico. AN acknowledges support from grant CONACyT research fellow -Instituto de Astronomia, UNAM, Mexico D.F. 04510, Mexico. We thank to Pu Du for his help and constructive comments. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/.The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.

Aims. The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) was designed to conduct a blind search for AGN-driven outflows on X-ray-selected AGNs at redshift z ∼ 2 with high (∼2 kpc) spatial resolution, and to correlate them with the properties of their host galaxy and central black hole. The main aims of this paper are: (a) to derive reliable estimates for the masses of the black holes and accretion rates for the Type-1 AGNs in this survey; and (b) to characterise the properties of the AGN-driven winds in the broad line region (BLR). Methods. We analysed rest-frame optical and UV spectra of 21 Type-1 AGNs. We used Hα, Hβ, and MgII line profiles to estimate the masses of the black holes. We used the blueshift of the CIV line profile to trace the presence of winds in the BLR. Results. We find that the Hα and Hβ line widths are strongly correlated, as is the line continuum luminosity at 5100 Å with Hα line luminosity, resulting in a well-defined correlation between black hole masses estimated from Hα and Hβ. Using these lines, we estimate that the black hole masses for our objects are in the range Log (MBH/M⊙) = 8.4–10.8 and are accreting at λEdd = 0.04–1.3. Furthermore, we confirm the well-known finding that the CIV line width does not correlate with the Balmer lines and the peak of the line profile is blueshifted with respect to the [OIII]-based systemic redshift. These findings support the idea that the CIV line is tracing outflowing gas in the BLR for which we estimated velocities up to ∼4700 km s−1. We confirm the strong dependence of the BLR wind velocity on the UV-to-X-ray continuum slope, the bolometric luminosity, and Eddington ratio. We infer BLR mass outflow rates in the range 0.005–3 M⊙ yr−1, revealing a correlation with the bolometric luminosity consistent with that observed for ionised winds in the narrow line region (NLR), and X-ray winds detected in local AGNs, and kinetic power ∼10−7 − 10−4 × LBol. The coupling efficiencies predicted by AGN-feedback models are much higher than the values reported for the BLR winds in the SUPER sample; although it should be noted that only a fraction of the energy injected by the AGN into the surrounding medium is expected to become kinetic power in the outflow. Finally, we find an anti-correlation between the equivalent width of the [OIII] line and the CIV velocity shift, and a positive correlation between this latter parameter and [OIII] outflow velocity. These findings, for the first time in an unbiased sample of AGNs at z ∼ 2, support a scenario where BLR winds are connected to galaxy-scale detected outflows, and are therefore capable of affecting the gas in the NLR located at kiloparsec scale distances. We thank the anonymous referee for the useful comments that improved the paper. The scientific results reported in this article are based on observations collected at the European organisation for Astronomical Research in the Southern Hemisphere under ESO programs 196.A-0377 (SINFONI data) and LP175.A-0839 (zCOSMOS data), and SDSS. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the US Department of Energy Office of Science. The SDSS-III web site is http://www.sdss3.org/. SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University. This research project was supported by the DFG Cluster of Excellence “Origin and Structure of the Universe” (universe-cluster.de). G. V. acknowledges financial support from Premiale 2015 MITic (PI B. Garilli). G. C., M. Bi, C. F. and F. F. acknowledge support from PRIN MIUR project “Black Hole winds and the Baryon Life Cycle of Galaxies: the stone-guest at the galaxy evolution supper”, contract #2017PH3WAT. M. P. is supported by the Programa Atracción de Talento de la Comunidad de Madrid via grant 2018-T2/TIC-11715.. With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737) Peer review

Context. The existence of tight correlations between supermassive black holes (BHs) and their host galaxies' properties in the local Universe suggests a closely linked evolution. Investigating these relations up to the high redshifts (z ≳ 6) is crucial in order to understand the interplay between star formation and BH growth across the cosmic time and to set constraints on galaxy formation and evolution models. In this work, we focus on the relation between BH mass (MBH) and the dynamical mass (Mdyn) of the host galaxy. Aims. Previous works suggest an evolution of the MBH-Mdyn relation with redshift indicating that BH growth precedes the galaxy mass assembly during their co-evolution at z > 3. However, dynamical galaxy masses at high redshift are often estimated through the virial theorem, thus introducing significant uncertainties. Within the scope of this work, our aim is to study the MBH-Mdyn relation of a sample of 2 < z < 7 quasars by constraining their galaxy masses through a full kinematical modelling of the cold gas kinematics, thus avoiding all possible biases and effects introduced by the rough estimates usually adopted so far. Methods. For this purpose, we retrieved public observations of 72 quasar host galaxies observed in [CII]158 μm or CO transitions with the Atacama Large Millimeter Array (ALMA). We then selected those quasars whose line emission is spatially resolved, and performed a kinematic analysis on ALMA observations. We estimated the dynamical mass of the systems by modelling the gas kinematics with a rotating disc, taking into account geometrical and instrumental effects. Our dynamical mass estimates, combined with MBH obtained from literature and our own new CIVλ1550 observations allowed us to investigate the MBH/Mdyn in the early Universe. Results. Overall, we obtained a sample of ten quasars at z ∼ 2-7, in which line emission is detected with high S/N (≳ 5-10) and the gas kinematics are spatially resolved and dominated by ordered rotation. The estimated dynamical masses place six out of ten quasars above the local relation yielding to MBH/Mdyn ratios ∼10× higher than those estimated in low-z galaxies. On the other hand, we found that four quasars at z ∼ 4-6 have dynamical-to-BH-mass ratios consistent with what is observed in early-type galaxies in the local Universe., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. We use the same method to analyse the two simulations, post-processing MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies., MNRAS in press

Aims. The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) was designed to conduct a blind search for AGN-driven outflows on X-ray-selected AGNs at redshift z∼2 with high (∼2 kpc) spatial resolution, and to correlate them with the properties of their host galaxy and central black hole. The main aims of this paper are: (a) to derive reliable estimates for the masses of the black holes and accretion rates for the Type-1 AGNs in this survey; and (b) to characterise the properties of the AGN-driven winds in the broad line region (BLR). Methods. We analysed rest-frame optical and UV spectra of 21 Type-1 AGNs. We used Hα, Hβ, and MgII line profiles to estimate the masses of the black holes. We used the blueshift of the CIV line profile to trace the presence of winds in the BLR. Results. We find that the Hα and Hβ line widths are strongly correlated, as is the line continuum luminosity at 5100 Å with Hα line luminosity, resulting in a well-defined correlation between black hole masses estimated from Hα and Hβ. Using these lines, we estimate that the black hole masses for our objects are in the range Log (MBH/M·) = 8.4-10.8 and are accreting at λEdd = 0.04-1.3. Furthermore, we confirm the well-known finding that the CIV line width does not correlate with the Balmer lines and the peak of the line profile is blueshifted with respect to the [OIII]-based systemic redshift. These findings support the idea that the CIV line is tracing outflowing gas in the BLR for which we estimated velocities up to ∼4700 km s-1. We confirm the strong dependence of the BLR wind velocity on the UV-to-X-ray continuum slope, the bolometric luminosity, and Eddington ratio. We infer BLR mass outflow rates in the range 0.005-3 M· yr-1, revealing a correlation with the bolometric luminosity consistent with that observed for ionised winds in the narrow line region (NLR), and X-ray winds detected in local AGNs, and kinetic power ∼10-7-10-4 × LBol. The coupling efficiencies predicted by AGN-feedback models are much higher than the values reported for the BLR winds in the SUPER sample; although it should be noted that only a fraction of the energy injected by the AGN into the surrounding medium is expected to become kinetic power in the outflow. Finally, we find an anti-correlation between the equivalent width of the [OIII] line and the CIV velocity shift, and a positive correlation between this latter parameter and [OIII] outflow velocity. These findings, for the first time in an unbiased sample of AGNs at z∼2, support a scenario where BLR winds are connected to galaxy-scale detected outflows, and are therefore capable of affecting the gas in the NLR located at kiloparsec scale distances., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

We present a description of the Australian Dark Energy Survey (OzDES) and summarize the results from its 6 years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-m Anglo-Australian Telescope, OzDES has monitored 771 active galactic nuclei, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30 000 sources, some as faint as rAB = 24 mag, and 375 000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to z ∼ 1.2 and the masses of hundreds of black holes out to z ∼ 4. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs.

We report on the spectroscopic confirmation of a large-scale structure around the luminous z = 6.31 quasi-stellar object (QSO) SDSS J1030+0524, powered by a one billion solar mass black hole. The structure is populated by at least six members, namely, four Lyman-break galaxies (LBGs), and two Lyman alpha emitters (LAEs). The four LBGs were identified among a sample of 21 i-band dropouts with zAB < 25.5 selected up to projected separations of 5 physical Mpc (15 arcmin) from the QSO. Their redshifts were determined through multi-object spectroscopic observations at 8-10 m class telescopes lasting up to eight hours. The two LAEs were identified in a 6 h VLT/MUSE observation centered on the QSO. The redshifts of the six galaxies cover the range between 6.129-6.355. Assuming that the peculiar velocities are negligible, this range corresponds to radial separations of ±5 physical Mpc from the QSO, that is comparable to the projected scale of the observed LBG distribution on the sky. We conservatively estimate that this structure is significant at a level > 3.5σ and that the level of the galaxy overdensity is at least 1.5-2 within the large volume sampled (∼780 physical Mpc3). The spectral properties of the six member galaxies (Lyα strength and UV luminosity) are similar to those of field galaxies at similar redshifts. This is the first spectroscopic identification of a galaxy overdensity around a supermassive black hole in the first billion years of the Universe. Our finding lends support to the idea that the most distant and massive black holes form and grow within massive (>1012 M) dark matter halos in large-scale structures and that the absence of earlier detections of such systems is likely due to observational limitations., We acknowledge funding from the INAF main-stream (1.05.01.86.31) and from the agreement ASI-INAF n. 2017-14-H.O. This work is based on ESO program 0102.A-0263(A), on observations obtained at the W. M. Keck Observatory, and on data acquired using the Large Binocular Telescope (ID 2017/2018 #18).

We present a description of the Australian Dark Energy Survey (OzDES) and summarize the results from its 6 years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-m Anglo-Australian Telescope, OzDES has monitored 771 active galactic nuclei, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30 000 sources, some as faint as r = 24 mag, and 375 000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to z ∼1.2 and the masses of hundreds of black holes out to z ∼4. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs., Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Fi-nanciadora de Estudos e Projetos, Fundac¸ão Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemein-schaft, and the Collaborating Institutions in the DES. LG was funded by the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement no. 839090. The DES data management system is supported by the National Science Foundation under Grant Numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2015-71825, ESP2015-66861, FPA2015-68048, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciência e Tecnologia (INCT) e-Universe (CNPq grant 465376/2014-2). Parts of this research were supported by the Australian Research Council under grants DP160100930, FL180100168, and FT140101270.

Accretion of stars on massive black holes (MBHs) can feed MBHs and generate tidal disruption events (TDEs). We introduce a new physically motivated model to self-consistently treat TDEs in cosmological simulations, and apply it to the assembly of a galaxy with final mass $3\times 10^{10}\mathrm{\, M}_\odot$ at $z=6$. This galaxy exhibits a TDE rate of $\sim 10^{-5}\mathrm{\, yr}^{-1}$, consistent with local observations but already in place when the Universe was one billion year old. A fraction of the disrupted stars participate in the growth of MBHs, dominating it until the MBH reaches mass $\sim 5 \times 10^5 \mathrm{\, M}_\odot$, but their contribution then becomes negligible compared to gas. TDEs could be a viable mechanism to grow light MBH seeds, but fewer TDEs are expected when the MBH becomes sufficiently massive to reach the luminosity of, and be detected as, an active galactic nucleus. Galaxy mergers bring multiple MBHs in the galaxy, resulting in an enhancement of the global TDE rate in the galaxy by $\sim 1$ order of magnitude during $100\mathrm{\, Myr}$ around mergers. This enhancement is not on the central MBH, but caused by the presence of MBHs in the infalling galaxies. This is the first self consistent study of TDEs in a cosmological environment and highlights that accretion of stars and TDEs are a natural process occurring in a Milky~Way-mass galaxy at early cosmic times., Comment: submitted to MNRAS

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh–Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work, we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh–Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh–Mstar relation requires a mean radiative efficiency ε ≳ 0.15, in line with theoretical expectations for accretion on to spinning black holes. However, matching the ‘raw’ observed relation for inactive black holes requires ε ∼ 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve ε ∼ 0.12–0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized., FS acknowledges partial support from a Leverhulme Trust Research Fellowship. RC acknowledges financial support from CONICYT Doctorado Nacional N° 21161487 and CONICYT PIA ACT172033. DMA thanks the Science and Technology Facilities Council (STFC) for support from grant no. ST/L00075X/1. ID is supported by the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement no. 788679. MM acknowledges support from the Beatriu de Pinos fellowship (2017-BP-00114).

69 pages, 21 figures. Accepted for publication in New Astronomy Reviews © 2020 Elsevier B.V., The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.© 2020 Elsevier B.V., All authors acknowledge the hospitality of the Lorentz Center for international workshops (Leiden, The Netherlands), where the idea of this work was born, and acknowledge the support of the International Space Science Institute (ISSI Bern, Switzerland), where the collaboration was originated. We thank the reviewers for having provided valuable and constructive comments that improved the clarity of the manuscript, and J.E. Barnes, L. Blecha, M. Eracleous, B.D. Farris, H. Fu, X. Liu, R. Pfeifle, L.C. Popovic, C. Ricci, S. Rodriguez, S. Tang, G. B. Taylor for their kind permission to reuse figures from their publications. ADR, CV and SB acknowledge financial support from ASI under grant ASI-INAF I/037/12/0 , and from the agreement ASI-INAF n. 2017-14-H.O . TB acknowledges support by the National Aeronautics and Space Administration under Grant No. NNX15AK84G and 80NSSC19K0319 issued through the Astrophysics Theory Program and by the Research Corporation for Science Advancement through a Cottrell Scholar Award. SF and KÉG thank the Hungarian National Research, Development, and Innovation Office ( OTKA NN110333 ) for support. KÉG was supported by the János Bolyai Research Scholarship of the Hungarian Academy of Sciences and by the ÚNKP-19-4 New National Excellence Program of the Ministry for Innovation and Technology. BH acknowledges financial support by the DFG grant GE625/17-1 . EL is supported by a European Union COFUND/Durham Junior Research Fellowship (under EU grant agreement no. 609412). DL acknowledges support from the European Research Council (ERC) under grant 647208 (PI Jonker). AS is supported by the ERC CoG grant 818691 (B Massive). MGi is supported by the “Programa de Atracción de Talento” of the Comunidad de Madrid grant 2018-T1/TIC-11733 for the project: “Unveiling Black Hole Winds from Space”, and by the Spanish State Research Agency (AEI) grant MDM-2017-0737 Unidad de Excelencia “María de Maeztu” - Centro de Astrobiología (INTA-CSIC). PRC and LM acknowledge support from the Swiss National Science Foundation under the Grant 200020_178949 . MC acknowledges support from the National Science Foundation (NSF) NANOGrav Physics Frontier Center, award number 1430284 . NHR acknowledges support from the BMBF Verbundforschung under FKZ: 05A17PC1 and FKZ: 05A17PC2. ZH acknowledges support from NASA grants NNX17AL82G and 80NSSC19K0149 and NSF grant 1715661 . KI acknowledges support by the Spanish MINECO under grant AYA2016-76012-C3-1-P and MDM-2014-0369 of ICCUB (Unidad de Excelencia ’María de Maeztu’). MPT acknowledges financial support from the Spanish MCIU through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía ( SEV-2017-0709 ) and through the MINECO grants AYA2012-38491-C02-02 and AYA2015-63939-C2-1-P .

Supermassive primordial stars forming in atomically-cooled halos at z∼15−20 are currently thought to be the progenitors of the earliest quasars in the Universe. In this picture, the star evolves under accretion rates of 0.1−1 M⊙ yr−1 until the general relativistic instability triggers its collapse to a black hole at masses of ∼105 M⊙. However, the ability of the accretion flow to sustain such high rates depends crucially on the photospheric properties of the accreting star, because its ionising radiation could reduce or even halt accretion. Here we present new models of supermassive Population III protostars accreting at rates 0.001−10 M⊙ yr−1, computed with the GENEVA stellar evolution code including general relativistic corrections to the internal structure. We use the polytropic stability criterion to estimate the mass at which the collapse occurs, which has been shown to give a lower limit of the actual mass at collapse in recent hydrodynamic simulations. We find that at accretion rates higher than 0.001 M⊙ yr−1 the stars evolve as red, cool supergiants with surface temperatures below 104 K towards masses >105M⊙, and become blue and hot, with surface temperatures above 105 K,only for rates ≲0.001 M⊙ yr−1. Compared to previous studies, our results extend the range of masses and accretion rates at which the ionising feedback remains weak, reinforcing the case for direct collapse as the origin of the first quasars.

Recent results have suggested that active galactic nuclei (AGN) could provide enough photons to reionise the Universe. We assess the viability of this scenario using a semi-numerical framework for modeling reionisation, to which we add a quasar contribution by constructing a Quasar Halo Occupation Distribution (QHOD) based on Giallongo et al. observations. Assuming a constant QHOD, we find that an AGN-only model cannot simultaneously match observations of the optical depth $\tau_e$, neutral fraction, and ionising emissivity. Such a model predicts $\tau_e$ too low by $\sim 2\sigma$ relative to Planck constraints, and reionises the Universe at $z\lesssim 5$. Arbitrarily increasing the AGN emissivity to match these results yields a strong mismatch with the observed ionising emissivity at $z\sim 5$. If we instead assume a redshift-independent AGN luminosity function yielding an emissivity evolution like that assumed in Madau & Haardt model, then we can match $\tau_e$ albeit with late reionisation, however such evolution is inconsistent with observations at $z\sim 4-6$ and poorly motivated physically. These results arise because AGN are more biased towards massive halos than typical reionising galaxies, resulting in stronger clustering and later formation times. AGN-dominated models produce larger ionising bubbles that are reflected in $\sim\times 2$ more 21cm power on all scales. A model with equal parts galaxies and AGN contribution is still (barely) consistent with observations, but could be distinguished using next-generation 21cm experiments HERA and SKA-low. We conclude that, even with recent claims of more faint AGN than previously thought, AGN are highly unlikely to dominate the ionising photon budget for reionisation., Comment: 16 pages, 9 figures, matches the accepted version for publication in MNRAS, 2017

We present near-infrared 1.1-1.3 and 1.3-1.6 $\mu$m spectra of the Ly$\alpha$-luminous quasar PSO J006.1240+39.2219 at $z=6.617$ obtained with the NIRSPEC spectrograph at the Keck-II telescope. The spectra cover the CIV $\lambda$1549, CIII] $\lambda$1909 emission lines and part of the UV continuum of the quasar. From the NIRSPEC observations of PSO J006.1240+39.2219, we constrain the spectral slope of its UV continuum to be $\alpha_{\lambda}=-1.35\pm0.26$ and measure an absolute magnitude of $M_{1450}=-25.60$. Using the scaling relation between black hole mass, width of the CIV line and ultraviolet continuum luminosity, we derive a black hole mass of $(2.19\pm0.30)\times 10^8 M_{sun}$, which is consistent but somewhat smaller than the typical black hole masses of $z\gtrsim6$ quasars of similar luminosities. The inferred accretion rate of $L_{bol}/L_{edd}\gtrsim2$ indicates that PSO J006.1240+39.2219 is in the phase of the rapid growth of its supermassive black hole characterized by the high NV/CIV line ratio, NV/CIV$>1$, and lower level of ionization of its circumnuclear gas than in other high-redshift luminous quasars. The NV/CIV line ratio of PSO J006.1240+39.2219 implies relatively high abundance of nitrogen in its circumnuclear gas. This abundance might be produced by the post-starburst population of stars that provide the fuel for black hole accretion., Comment: 13 pages, 8 figures, submitted to ApJ

We present near-infrared 1.1-1.3 and 1.3-1.6 μm spectra of the Ly-luminous quasar PSO J006.1240+39.2219 at z = 6.617 ± 0.003 obtained with the NIRSPEC spectrograph at the Keck II telescope. The spectra cover the C iv λ1549, C iii] λ1909 emission lines and part of the UV continuum of the quasar. From the NIRSPEC observations of PSO J006.1240+39.2219, we constrain the spectral slope of its UV continuum to be λ =-1.35 ± 0.26 and measure an absolute magnitude of M 1450 =-25.60. Using the scaling relation between black hole mass, width of the C iv line and ultraviolet continuum luminosity, we derive a black hole mass of (2.19 ± 0.30) × 108 M o, which is consistent but somewhat smaller than the typical black hole masses of z 6 quasars of similar luminosities. The inferred accretion rate of L Bol/L Edd 2 indicates that PSO J006.1240+39.2219 is in the phase of the rapid growth of its supermassive black hole characterized by the high N v/C iv line ratio, N v/C iv > 1, and the lower level of ionization of its circumnuclear gas than in other high-redshift luminous quasars. The N v/C iv line ratio of PSO J006.1240+39.2219 implies a relatively high abundance of nitrogen in its circumnuclear gas. This abundance might be produced by the post-starburst population of stars that provide the fuel for black hole accretion.

One of the ideas to explain the existence of supermassive black holes (SMBH) that are in place by z~7 is that there was an earlier phase of very rapid accretion onto direct collapse black holes (DCBH) that started their lives with masses ~ 10^4-10^5 M_solar. Working in this scenario, we show that the mass function of SMBH after such a limited time period with growing formation rate paired with super-Eddington accretion can be described as a broken power-law with two characteristic features. There is a power-law at intermediate masses whose index is the dimensionless ratio {\alpha} = {\lambda}/{\gamma}, where {\lambda} is the growth rate of the number of DCBH during their formation era, and {\gamma} is the growth rate of DCBH masses by super-Eddington accretion during the DCBH growth era. A second feature is a break in the power law profile at high masses, above which the mass function declines rapidly. The location of the break is related to the dimensionless number \b{eta} = {\gamma} T, where T is the duration of the period of DCBH growth. If the SMBH continue to grow at later times at an Eddington-limited accretion rate, then the observed quasar luminosity function can be directly related to the tapered power-law function derived in this paper., Comment: 6 pages, 3 figures, accepted to ApJ

Context. The flat-spectrum radio quasar 4C +71.07 is a high-redshift (z = 2.172), γ-loud blazar whose optical emission is dominated by thermal radiation from the accretion disc. Aims. 4C +71.07 has been detected in outburst twice by the AGILE γ-ray satellite during the period from the end of October to mid-November 2015, when it reached a γ-ray flux of the order of F(E > 100 MeV)=(1.2 ± 0.3)×10 photons cm s and F(E > 100 MeV)=(3.1 ± 0.6)×10 photons cm s, respectively, allowing us to investigate the properties of the jet and the emission region. Methods. We investigated its spectral energy distribution by means of almost-simultaneous observations covering the cm, mm, near-infrared, optical, ultraviolet, X-ray, and γ-ray energy bands obtained by the GASP-WEBT Consortium and the Swift, AGILE, and Fermi satellites. Results. The spectral energy distribution of the second γ-ray flare (whose energy coverage is more dense) can be modelled by means of a one-zone leptonic model, yielding a total jet power of about 4 × 10 erg s. Conclusions. During the most prominent γ-ray flaring period our model is consistent with a dissipation region within the broad-line region. Moreover, this class of high-redshift, flat-spectrum radio quasars with high-mass black holes might be good targets for future γ-ray satellites such as e-ASTROGAM. © ESO 2019., AGILE is an ASI space mission developed with programmatic support by INAF and INFN. We acknowledge partial support through the ASI grant no. I/028/12/0. SV and PR acknowledge contract ASI-INAF I/004/11/0 and INAF/IASF Palermo where most of the work was carried out. SV acknowledges financial contribution from the agreement ASI-INAF no. 2017-14-H.0. Part of this work is based on archival data, software, or online services provided by the ASI SPACE SCIENCE DATA CENTER (ASI-SSDC). SV and PR thank Leonardo Barzaghi and Sara Baitieri for useful discussions. The Osservatorio di Torino team acknowledges the financial contribution from the agreement ASI-INAF No. 2017-14-H.0 and from the contract PRIN-SKA-CTA-INAF 2016. OMK acknowledges financial support from the Shota Rustaveli National Science Foundation under contract FR/217950/16 and grants NSFC11733001, NSFCU1531245. IA acknowledges support from a Ramón y Cajal grant of the Ministerio de Economía y Compet-itividad (MINECO) of Spain. The research at the IAA–CSIC was supported in part by the MINECO through grants AYA2016–80889–P, AYA2013–40825–P, and AYA2010–14844, and by the regional government of Andalucía through grant P09–FQM–4784. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain). Calar Alto Observatory is jointly operated by the MPIA and the IAA-CSIC. This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grant DN 08-1/2016. The St. Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. AZT-24 observations are made within an agreement among the Pulkovo, Rome, and Teramo observatories. GD and OV gratefully acknowledge the observing grant support from the Institute of Astronomy and Rozhen National Astronomical Observatory, Bulgaria Academy of Sciences, via bilateral joint research project “Observations of ICRF radio-sources visible in optical domain” (PI G. Damljanovic). This work is a part of Project No. 176011 (“Dynamics and kinematics of celestial bodies and systems”), No. 176004 (“Stellar physics”) and No. 176021 (“Visible and invisible matter in nearby galaxies: theory and observations”) supported by the Ministry of Education, Science, and Technological Development of the Republic of Serbia. The Maidanak Observatory team acknowledges support from Uzbekistan Academy of Sciences grants No. F2-FA-F027 and F.4-16.

Tables A.1–A.3 are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/630/A110, Context. When can an active galactic nucleus (AGN) be considered radio loud (RL)? Following the established view of the AGNs inner workings, an AGN is RL if associated with relativistic ejections emitting a radio synchrotron spectrum (i.e., it is a “jetted” AGN). In this paper we exploit the AGN main sequence that offers a powerful tool to contextualize radio properties. Aims. If large samples of optically-selected quasars are considered, AGNs are identified as RL if their Kellermann’s radio loudness ratio RK > 10. Our aims are to characterize the optical properties of different classes based on radio loudness within the main sequence and to test whether the condition RK > 10 is sufficient for the identification of RL AGNs, since the origin of relatively strong radio emission may not be necessarily due to relativistic ejection. Methods. A sample of 355 quasars was selected by cross-correlating the Very Large Array Faint Images of the Radio Sky at Twenty-Centimeters survey (FIRST) with the twelfth release of the Sloan Digital Sky Survey Quasar Catalog published in 2017. We classified the optical spectra according to their spectral types along the main sequence of quasars. For each spectral type, we distinguished compact and extended morphology (providing a FIRST-based atlas of radio maps in the latter case), and three classes of radio loudness: detected ( specific flux ratio in the g band and at 1.4 GHz, R′K < 10), intermediate (10 ≤ R′K < 70), and RL (R′K ≥ 70). Results. The analysis revealed systematic differences between radio-detected (i.e., radio-quiet), radio-intermediate, and RL classes in each spectral type along the main sequence. We show that spectral bins that contain the extreme Population A sources have radio power compatible with emission by mechanisms ultimately due to star formation processes. RL sources of Population B are characteristically jetted. Their broad Hβ profiles can be interpreted as due to a binary broad-line region. We suggest that RL Population B sources should be preferential targets for the search of black hole binaries, and present a sample of binary black hole AGN candidates. Conclusions. The validity of the Kellermann’s criterion may be dependent on the source location along the quasar main sequence. The consideration of the main sequence trends allowed us to distinguish between sources whose radio emission mechanisms is jetted from the ones where the mechanism is likely to be fundamentally different.© ESO 2019, PM and MDO acknowledge funding from the INAF PRIN-SKA 2017 program 1.05.01.88.04. PM also acknowledges the Programa de Estancias de Investigación (PREI) No. DGAP/DFA/2192/2018 of UNAM. AdO acknowledges financial support from the Spanish Ministry of Economy and Competitiveness through grant AYA2016-76682-C3-1-P and from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). The relevant research is part of the project 176001 “Astrophysical spectroscopy of extragalactic objects” supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia

We present calculations on the formation of massive black holes with 10^5 Msun at z > 6 that can be the seeds of supermassive black holes at z > 6. Under the assumption of compact star cluster formation in merging galaxies, star clusters in haloes of 10^8 ~ 10^9 Msun can undergo rapid core-collapse leading to the formation of very massive stars (VMSs) with ~1000 Msun which directly collapse into black holes with similar masses. Star clusters in halos of > 10^9 Msun experience type-II supernovae before the formation of VMSs due to long core-collapse time scales. We also model the subsequent growth of black holes via accretion of residual stars in clusters. 2-body relaxation efficiently re-fills the loss cones of stellar orbits at larger radii and resonant relaxation at small radii is the main driver for accretion of stars onto black holes. As a result, more than ninety percent of stars in the initial cluster are swallowed by the central black holes before z=6. Using dark matter merger trees we derive black hole mass functions at z=6-20. The mass function ranges from 10^3 to 10^5 Msun at z ~ 4*10^8 Msun at z ~ 20 successfully leads to the formation of >~ 10^5 Msun BHs by z >~ 10 which can be the potential seeds of supermassive black holes seen today., Comment: 11 pages, 10 figures, accepted for publication in MNRAS