Search

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)

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).

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.

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).

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.

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 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

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

Aims. We present a multifrequency, dense radio monitoring program of the blazar OJ 287 using the 100-m E elsberg radio telescope. The program aims to test different binary supermassive black hole (SMBH) scenarios and studying the physical conditions in the central region of this bright blazar. Here, we analyze the evolution in total flux density, linear and circular polarization as a means to study the OJ 287 jet structure and its magnetic field geometry. Methods.We used a recently developed, high-precision data analysis methodology to recover all four Stokes parameters.We measured the total flux density of OJ 287 at nine bands from 2.64 GHz to 43 GHz, the linear polarization parameters at four bands between 2.64 GHz and 10.45 GHz, and the circular polarization at two bands, 4.85 GHz and 8.35 GHz. The mean cadence of our measurements is ten days. Results. Between December 2015 and January 2017 (MJD 57370-57785), OJ 287 showed flaring activity and complex linear and circular polarization behavior. The radio electric vector position angle (EVPA) showed a large clockwise (CW) rotation by 340 with a mean rate of-1.04°/day. Based on concurrent very long baseline interferometric (VLBI) polarization data at 15 GHz and 43 GHz, the rotation seems to originate within the jet core at 43 GHz (projected angular size 0.15 mas or 0.67 pc at the redshift of the source). Moreover, optical polarization data show a similar monotonic CW rotation with a rate of about-1.1°/day which is superposed with shorter and faster rotations that exhibit rates of about 7.8°/day, mainly in the CW sense. Conclusions. The flux density and polarization variability of the single dish, VLBI and optical data is consistent with a polarized emission component propagating on a helical trajectory within a bent jet. We constrained the helix arc length to 0.26 pc and radius to 0.04 pc as well as the jet bending arc length projected on the plane of the sky to ≤1.9-7.6 pc. A similar bending has been observed also in high angular resolution VLBI images of the OJ 287 jet at its innermost regions. The helical trajectory covers only a part of the jet width, possibly its spine. In addition, our results indicate the presence of a stable polarized emission component. Its EVPA (-10°) is oriented perpendicular to the large scale jet, suggesting dominance of the poloidal magnetic field component. Finally, the EVPA rotation begins simultaneously with an optical flare and hence the two might be physically connected. That optical flare has been suggested to be linked to the interaction of a secondary SMBH with the inner accretion disk or originating in the jet of the primary.© 2018 ESO., This program is supported by Fermi Guest Investigator grants NNX08AW56G, NNX09AU10G, NNX12AO93G, and NNX15AU81G.

We report the first results of a survey on 74 narrow-line Seyfert 1 galaxies (NLS1s) carried out in 2015 with the Karl G. Jansky Very Large Array (JVLA) at 5 GHz in A-configuration. So far, this is the largest survey aimed to image the radio continuum of NLS1s. We produced radio maps in order to compare the general properties of three different samples of objects: radio-quiet NLS1s (RQNLS1s), steep-spectrum radio-loud NLS1s (S-NLS1s), and flat-spectrum radio-loud NLS1s (F-NLS1s). We find that the three classes correspond to different radio morphologies, with F-NLS1s being more compact, and RQNLS1s often showing diffuse emission on kpc scales. We also find that F-NLS1s might be low-luminosity and possibly young blazars, and that S-NLS1s are part of the parent population of F-NLS1s. Dedicated studies to RQNLS1s are needed in order to fully understand their role in the unification pictures., 36 pages, 45 figures, 10 tables, accepted for publication by Astronomy & Astrophysics

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}, 21 pages, 12 figures, 9 tables. Published in The Astrophysical Journal. Please send any comments/questions to S. R. Taylor (srtaylor@caltech.edu)

The population of young stars near the Supermassive black hole (SMBH) in the Galactic Center (GC) has presented an unexpected challenge to theories of star formation. Kinematics measurements of these stars have revealed a stellar disk structure (with an apparent 20% disk membership) that has provided important clues to the origin of these mysterious young stars. However many of the apparent disk properties are difficult to explain, including the low disk membership fraction and the high eccentricities, given the youth of this population. Thus far, all efforts to derive the properties of this disk have made the simplifying assumption that stars at the GC are single stars. Nevertheless, stellar binaries are prevalent in our Galaxy, and recent investigations suggested that they may also be abundant in the Galactic Center. Here we show that binaries in the disk can largely alter the apparent orbital properties of the disk. The motion of binary members around each other adds a velocity component, which can be comparable to the magnitude of the velocity around the SMBH in the GC. Thus neglecting the contribution of binaries can significantly vary the inferred stars' orbital properties. While the disk orientation is unaffected the apparent disk's 2D width is increased to about 11.2deg, similar to the observed width. For a population of stars orbiting the SMBH with zero eccentricity, unaccounted for binaries will create a wide apparent eccentricity distribution with an average of 0.23.This is consistent with the observed average eccentricity of the stars' in the disk. We suggest that this high eccentricity value, which poses a theoretical challenge, may be an artifact of binary stars. Finally our results suggest that the actual disk membership might be significantly higher than the one inferred by observations that ignore the contribution of binaries, alleviating another theoretical challenge., Comment: Accepted to ApJ Lett, 5 figures, 8 pages

Context. Quasars radiating at extreme Eddington ratios (hereafter xA quasars) are likely a prime mover of galactic evolution and have been hailed as potential distance indicators. Their properties are still scarcely known. Aims. We aim to test the effectiveness of the selection criteria defined on the >4D Eigenvector 1> (4DE1) for identifying xA sources. We provide a quantitative description of their rest-frame UV spectra (1300-2200 Å) in the redshift range 2 ≤ z ≤ 2.9, with a focus on major emission features. Methods. Nineteen extreme quasar candidates were identified using 4DE1 selection criteria applied to SDSS spectra: Alλ1860Siλ1892 ≥ 0.5 and Cλ1909/Siλ1892 ≤ 1. The emission line spectra was studied using multicomponent fits of deep spectroscopic observations (S/N ≥ 40-50; spectral resolution ≈ 250 km s) obtained with the OSIRIS at Gran Telescopio Canarias (GTC). Results. GTC spectra confirm that almost all of these quasars are xA sources with very similar properties. We provide spectrophotometric and line profile measurements for the Siλ1397+O, Cλ1549+Heλ1640, and the 1900 Å blend. This last feature is found to be predominantly composed of Alλ1860, Siλ1892 and Fe emission features, with weak Cλ1909. The spectra can be characterized as very low ionization (ionization parameter, logU ∼ -3), a condition that explains the significant Fe emission observed in the spectra. xA quasars show extreme properties in terms of Cλ1549 equivalent width and blueshift amplitudes. Cλ1549 shows low equivalent width, with a median value of 15 Å (≤ 30 Å for the most sources), and high or extreme blueshift amplitudes (-5000 ≤ c(1/2) ≤ -1000 km s). Weak-lined quasars appear as extreme xA quasars and not as an independent class. The Cλ1549 high amplitude blueshifts coexists in all cases save one with symmetric and narrower Alλ1860 and Siλ1892 profiles. Estimates of the Eddington ratio using the Alλ1860 FWHM as a virial broadening estimator are consistent with the ones of a previous xA sample. Conclusions. xA quasars show distinguishing properties that make them easily identifiable in large surveys and potential >standard candles> for cosmological applications. It is now feasible to assemble large samples of xA quasars from the latest data releases of the SDSS. We provide evidence that Alλ1860 could be associated with a low-ionization virialized subsystem, supporting previous suggestions that Al is a reliable virial broadening estimator.© ESO 2018., AdO, MLMA and JWS acknowledge financial support from the Spanish Ministry for Economy and Competitiveness through grants AYA2013-42227-P and AYA2016-76682-C3-1-P. MLMA, PM and MDO acknowledge funding from the INAF PRIN-SKA 2017 program 1.05.01.88.04. JP acknowledge financial support from the Spanish Ministry for Economy and Competitiveness through grants AYA2013-40609-P and AYA2016-76682-C3-3-P.

We present an ALMA high-resolution observation of the 840 um continuum and [CII] line emission in the WISE-SDSS selected hyper-luminous (WISSH) QSO J1015+0020 at z~4.4. Our analysis reveals an exceptional overdensity of [CII]-emitting companions with a very small (, Accepted

Open Access.--Attribution 4.0 International (CC BY 4.0), We present preliminary results of a spectroscopic analysis for a sample of type 1 highly accreting quasars (L/L-Edd similar to 1.0) at high redshift, z similar to 2-3. The quasars were observed with the OSIRIS spectrograph on the GTC 10.4 m telescope located at the Observatorio del Roque de los Muchachos in La Palma. The highly accreting quasars were identified using the 4D Eigenvector 1 formalism, which is able to organize type 1 quasars over a broad range of redshift and luminosity. The kinematic and physical properties of the broad line region have been derived by fitting the profiles of strong UV emission lines such as AlIII lambda 1860, SiIII]lambda 1892 and CIII]lambda 1909. The majority of our sources show strong blueshifts in the high-ionization lines and high Eddington ratios which are related with the productions of outflows. The importance of highly accreting quasars goes beyond a detailed understanding of their physics: their extreme Eddington ratio makes them candidates standard candles for cosmological studies.© Copyright © 2018 Martínez-Aldama, Del Olmo, Marziani, Sulentic, Negrete, Dultzin, Perea and D'Onofrio., MM-A acknowledge the postdoctoral grant from the CONACyT. MM-A and AD acknowledge financial support from Spanish Ministry for Economy and Competitiveness through grants AYA2013-42227-P and AYA2016-76682-C3-3-1-P.

For a sample of star forming galaxies in the redshift interval 0.15$, Comment: 31 pages, 23 figures

We introduce massive black holes (BHs) in the Feedback In Realistic Environments project and perform high-resolution cosmological hydrodynamic simulations of quasar-mass halos ($M_{\rm halo}(z=2) \approx 10^{12.5}\,\rm{M}_{\odot}$) down to $z=1$. These simulations model stellar feedback by supernovae, stellar winds, and radiation, and BH growth using a gravitational torque-based prescription tied to resolved properties of galactic nuclei. We do not include BH feedback. We show that early BH growth occurs through short ($\lesssim 1\,$Myr) accretion episodes that can reach or even exceed the Eddington rate. In this regime, BH growth is limited by bursty stellar feedback continuously evacuating gas from galactic nuclei, and BHs remain under-massive relative to the local $M_{\rm BH}$-$M_{\rm bulge}$ relation. BH growth is more efficient at later times, when the nuclear stellar potential retains a significant gas reservoir, star formation becomes less bursty, and galaxies settle into a more ordered state, with BHs rapidly converging onto the scaling relation when the host reaches $M_{\rm bulge} \sim 10^{10}\,\rm{M}_{\odot}$. Our results are not sensitive to the details of the accretion model so long as BH growth is tied to the gas content within $\sim 100\,$pc of the BH. Our simulations imply that bursty stellar feedback has strong implications for BH and AGN demographics, especially in the early Universe and for low-mass galaxies., 5 pages, 3 figures, submitted to MNRAS

The relativistic jets created by some active galactic nuclei are important agents of AGN feedback. In spite of this, our understanding of what produces these jets is still incomplete. X-ray observations, which can probe the processes operating in the central regions in immediate vicinity of the supermassive black hole, the presumed jet launching point, are potentially particularly valuable in illuminating the jet formation process. Here, we present the hard X-ray NuSTAR observations of the radio-loud quasar 4C 74.26 in a joint analysis with quasi-simultaneous, soft X-ray Swift observations. Our spectral analysis reveals a high-energy cut-off of 183$_{-35}^{+51}$ keV and confirms the presence of ionized reflection in the source. From the average spectrum we detect that the accretion disk is mildly recessed with an inner radius of $R_\mathrm{in}=4-180\,R_\mathrm{g}$. However, no significant evolution of the inner radius is seen during the three months covered by our NuSTAR campaign. This lack of variation could mean that the jet formation in this radio-loud quasar differs from what is observed in broad-line radio galaxies., Accepted for publication in ApJ

We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispersion of $215.86 \pm 6.34$ km/s. From archival radio observations we see that this galaxy hosts a compact flat-spectrum radio source, and we conclude that it is the active galactic nucleus (AGN) responsible for generating the radio lobes. The radio luminosity distribution of the jets, and the broad-line classification of the host AGN, indicate this GRG is orientated well out of the plane of the sky, making its physical size one of the largest known for any GRG. Analysis of the infrared data suggests that the host is a lenticular type galaxy with a large stellar mass ($\log~\mathrm{M}/\mathrm{M}_\odot = 11.56 \pm 0.12$), and a moderate star formation rate ($1.2 \pm 0.3~\mathrm{M}_\odot/\mathrm{year}$). Spatially smoothing the SDSS images shows the system around UGC 9555 to be significantly disturbed, with a prominent extension to the south-east. Overall, the evidence suggests this host galaxy has undergone one or more recent moderate merger events and is also experiencing tidal interactions with surrounding galaxies, which have caused the star formation and provided the supply of gas to trigger and fuel the Mpc-scale radio lobes., Compressed for arXiv, see journal for full resolution images

PHL 6625 is a luminous quasi-stellar object (QSO) at z = 0.3954 located behind the nearby galaxy NGC 247 (z = 0.0005). Hubble Space Telescope (HST) observations revealed an arc structure associated with it. We report on spectroscopic observations with the Very Large Telescope (VLT) and multiwavelength observations from the radio to the X-ray band for the system, suggesting that PHL 6625 and the arc are a close pair of merging galaxies, instead of a strong gravitational lens system. The QSO host galaxy is estimated to be (4-28) x 10^10 M_sun, and the mass of the companion galaxy of is estimated to be M_* = (6.8 +/- 2.4) x 10^9 M_sun, suggesting that this is a minor merger system. The QSO displays typical broad emission lines, from which a black hole mass of about (2-5) x 10^8 M_sun and an Eddington ratio of about 0.01-0.05 can be inferred. The system represents an interesting and rare case where a QSO is associated with an ongoing minor merger, analogous to Arp 142., Comment: ApJ to appear

We use the coadded spectra of 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project observations of 482 quasars with z>1.46 to highlight systematic biases in the SDSS- and BOSS-pipeline redshifts due to the natural diversity of quasar properties. We investigate the characteristics of this bias by comparing the BOSS-pipeline redshifts to an estimate from the centroid of HeII 1640. HeII has a low equivalent width but is often well-defined in high-S/N spectra, does not suffer from self-absorption, and has a narrow component that, when present (the case for about half of our sources), produces a redshift estimate that, on average, is consistent with that determined from [OII] to within 1-sigma of the quadrature sum of the HeII and [OII] centroid measurement uncertainties. The large redshift differences of ~1000 km/s, on average, between the BOSS-pipeline and HeII-centroid redshifts suggest there are significant biases in a portion of BOSS quasar redshift measurements. Adopting the HeII-based redshifts shows that CIV does not exhibit a ubiquitous blueshift for all quasars, given the precision probed by our measurements. Instead, we find a distribution of CIV centroid blueshifts across our sample, with a dynamic range that (i) is wider than that previously reported for this line, and (ii) spans CIV centroids from those consistent with the systemic redshift to those with significant blueshifts of thousands of kilometers per second. These results have significant implications for measurement and use of high-redshift quasar properties and redshifts and studies based thereon., Comment: 13 pages, 7 figures, submitted to ApJ

A. Boehle et. al., We present new, more precise measurements of the mass and distance of our Galaxy's central supermassive black hole, Sgr A∗. These results stem from a new analysis that more than doubles the time baseline for astrometry of faint stars orbiting Sgr A∗, combining 2 decades of speckle imaging and adaptive optics data. Specifically, we improve our analysis of the speckle images by using information about a star's orbit from the deep adaptive optics data (2005-2013) to inform the search for the star in the speckle years (1995-2005). When this new analysis technique is combined with the first complete re-reduction of Keck Galactic Center speckle images using speckle holography, we are able to track the short-period star S0-38 (K-band magnitude = 17, orbital period = 19 yr) through the speckle years. We use the kinematic measurements from speckle holography and adaptive optics to estimate the orbits of S0-38 and S0-2 and thereby improve our constraints of the mass (M ) and distance (R ) of Sgr A∗: M = (4.02 ±0.16 ±0.04) x10 M and 7.86 ±0.14 ±0.04 kpc. The uncertainties in M and R as determined by the combined orbital fit of S0-2 and S0-38 are improved by a factor of 2 and 2.5, respectively, compared to an orbital fit of S0-2 alone and a factor of ∼2.5 compared to previous results from stellar orbits. This analysis also limits the extended dark mass within 0.01 pc to less than 0.13 x10 M at 99.7% confidence, a factor of 3 lower compared to prior work.© 2016. The American Astronomical Society. All rights reserved., RS acknowledges funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no. 614922.

We present an analysis of 43 years (1972 to 2015) of spectroscopic observations of the Seyfert 1 galaxy NGC 5548. This includes 12 years of new unpublished observations (2003 to 2015). We compiled about 1600 H$\beta$ spectra and analyzed the long-term spectral variations of the 5100 \AA\ continuum and the H$\beta$ line. Our analysis is based on standard procedures including the Lomb-Scargle method, which is known to be rather limited to such heterogeneous data sets, and new method developed specifically for this project that is more robust and reveals a $\sim$5700 day periodicity in the continuum light curve, the H$\beta$ light curve and the radial velocity curve of the red wing of the H$\beta$ line. The data are consistent with orbital motion inside the broad emission line region of the source. We discuss several possible mechanisms that can explain this periodicity, including orbiting dusty and dust-free clouds, a binary black hole system, tidal disruption events, and the effect of an orbiting star periodically passing through an accretion disc., Comment: Accepted in ApJS, 64 pages, 10 figures and 4 tables