Search

This is an Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Context. Active galactic nuclei (AGN) are thought to be responsible for the suppression of star formation in massive ∼1010 M⊙ galaxies. While this process is a key feature in numerical simulations of galaxy formation, it has not been unambiguously confirmed in observational studies yet. Aims. The characterization of the star formation rate (SFR) in AGN host galaxies is challenging as AGN light contaminates most SFR tracers. Furthermore, the various SFR tracers are sensitive to different timescales of star formation from approximately a few to 100 Myr. We aim to obtain and compare SFR estimates from different tracers for AGN host galaxies in the Close AGN Reference Survey (CARS) to provide new observational insights into the recent SFR history of those systems. Methods. We constructed integrated panchromatic spectral energy distributions to measure the far infrared (FIR) luminosity as a tracer for the recent (< 100 Myr) SFR. In addition we used the integral-field unit observation of the CARS targets to employ the Hα luminosity decontaminated by AGN excitation as a proxy for the current (< 5 Myr) SFR. Results. We find that significant differences in specific SFR of the AGN host galaxies as compared with the larger galaxy population disappear once cold gas mass, in addition to stellar mass, is used to predict the SFR for a specific AGN host. Only a tentative trend with the inclination of the host galaxy remains, such that SFR appears slightly lower than expected when the galaxies of unobscured AGN appear more edge-on along our line-of-sight, particular for dust-insensitive FIR-based SFRs. We identify individual galaxies with a significant difference in their SFR which can be related to a recent enhancement or decline in their SFR history that might be related to various processes including interactions, gas consumption, outflows, and AGN feedback. Conclusions. AGN can be present in various stages of galaxy evolution which makes it difficult to relate the SFR solely to the impact of the AGN. Our study shows that stellar mass alone is an insufficient parameter to estimate the expected SFR of an AGN host galaxy compared to the underlying non-AGN galaxy population. We do not find any strong evidence for a global positive or negative AGN feedback in the CARS sample. However, there is tentative evidence that (1) the relative orientation of the AGN engine with respect to the host galaxies might alter the efficiency of AGN feedback and that (2) the recent SFH is an additional tool to identify rapid changes in galaxy growth driven by the AGN or other processes. © I. Smirnova-Pinchukova et al. 2022., I.S.P. greatly appreciates financial support from the DLR through grant 50OR2006. B.H. is financially supported through DFG grant GE625/17-1. I.S.P. and B.H. also acknowledge travel support from the DAAD via grant 57509925. T.A.D. acknowledges support from the UK Science and Technology Facilities Council through grant ST/S00033X/1. M.G. acknowledges partial support by NASA Chandra GO8-19104X/GO9-20114X and HST GO-15890.020/023-A, and the BlackHoleWeather program. S.B., C.O., and M.S. acknowledge support from the Natural Sciences and Engineering Research Council (NSERC) of Canada. M.P.T. 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 grants PGC2018-098915-B-C21 and PID2020-117404GB-C21 (MCI/AEI/FEDER, UE). The work of T.C. was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA. Based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere under ESO programme 094.B-0345(A) and 095.B-0015(A). Based on observations collected at the Centro Astronómico Hispano-Alemán (CAHA) at Calar Alto, operated jointly by Junta de Andalucía and Consejo Superior de Investigaciones Científicas (IAA-CSIC). This project used data obtained with the Dark Energy Camera (DECam), which was constructed by the Dark Energy Survey (DES) Collaboration. Funding for the DES Projects has been provided by the US Department of Energy, the US 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 for Cosmological Physics at the University of Chicago, 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, Financiadora de Estudos e Projetos, Fundaçã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 Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Enérgeticas, Medioambientales y Tecnológicas–Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciències de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, NSF’s NOIRLab, the University of Nottingham, the Ohio State University, the OzDES Membership Consortium, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University. Based on observations at Cerro Tololo Inter-American Observatory, NSF’s NOIRLab (NOIRLab Prop. ID 2017A-0914; PI: Grant Tremblay), which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation. The James Clerk Maxwell Telescope is operated by the East Asian Observatory on behalf of The National Astronomical Observatory of Japan; Academia Sinica Institute of Astronomy and Astrophysics; the Korea Astronomy and Space Science Institute; Center for Astronomical Mega-Science (as well as the National Key R&D Program of China with No. 2017YFA0402700). Additional funding support is provided by the Science and Technology Facilities Council of the United Kingdom and participating universities and organizations in the United Kingdom and Canada. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research has made use of the NASA/IPAC Infrared Science Archive, which is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology. This research has made use of the APASS database, located at the AAVSO website. Funding for APASS has been provided by the Robert Martin Ayers Sciences Fund. Supported by the international Gemini Observatory, a program of NSF’s NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation, on behalf of the Gemini partnership of Argentina, Brazil, Canada, Chile, the Republic of Korea, and the United States of America. The Science, Technology and Facilities Council is acknowledged by JN for support through the Consolidated Grant Cosmology and Astrophysics at Portsmouth, ST/S000550/1.

HE 0435-5304 from Hamburg European Southern Observatory survey is a quasar that appears in the literature with two conflicting redshift values: $\sim 1.2$ and $\sim 0.4$. It was used in the studies of the intergalactic medium through fitting of the narrow absorption lines in its ultraviolet (UV) spectrum. This source is also known historically as a luminous infrared galaxy. We present optical spectra of HE 0435-5304, aiming to precisely measure its redshift and to study its physical properties. In particular, properties of its active nucleus, which is studied in the context of the source being identified here as an ultra-luminous infrared galaxy, allow us to place this quasar in the context of the general population. Fitting the spectra, we focused on modeling H$\beta$ and [O III] lines. Based on these, we derived the virial black hole mass, bolometric luminosity, and Eddington ratio of the active galactic nucleus (AGN). Additionally, we performed broad band photometry fitting which allows us to quantify host galaxy parameters. The improved redshift value of HE 0435-5304 is estimated to $0.42788 \pm 0.00027$ based on the [O II] line, which is mostly consistent with the narrowest components of the other emission lines. The source was found to be a relatively massive and luminous AGN whose host galaxy is actively forming stars. Although its stellar population seems to be heavily obscured, we did not find evidence for significant obscuration of the nucleus. We conclude that the AGN HE 0435-5304 is a rather prominent iron emitter from the extreme type-A population very close to the narrow-line Seyfert 1 group. The fact that the width of the H$\beta$ line appears to be systematically growing in its broadest component with time may suggest that this AGN is changing its broad line region. Due to the influence of atmospheric effects, this finding is uncertain., Comment: 19 pages, 11 figures, accepted for publication in A&A

Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilise the availability of radio and $\gamma$-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculate the multiband time lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process regression. We find that the jet particle and magnetic field energy densities are in equipartition ($k_\textrm{r} = 1.08\pm0.18$). The jet apex is located $z_\textrm{91.5 GHz}= 22 - 645$ $R_\textrm{s}$ ($2 - 20 \times 10^{-3}$ pc) upstream of the 3 mm radio core; at that position, the magnetic field amplitude is $B_\textrm{core}^\textrm{91.5 GHz}= 3 - 10$ G. Our results are in good agreement with earlier studies, which utilised very-long-baseline interferometry. Furthermore, we investigate the temporal relation between the ejection of radio and $\gamma$-ray flares. Our results are in favour of the $\gamma$-ray emission being associated with the radio emission. We are able to tentatively connect the ejection of features identified at 43 and 86 GHz to prominent $\gamma$-ray flares. Finally, we compute the multiplicity parameter $\lambda$ and the Michel magnetisation $\sigma_\textrm{M}$ and find that they are consistent with a jet launched by the Blandford & Znajek 1977 mechanism., Comment: 13 pages, 6 figures, accepted for publication in A&A

The REX (Radio-Emitting X-ray sources) is a catalogue produced by cross-matching X-ray data from the ROSAT-PSPC archive of pointed observations and radio data from the NRAO VLA Sky Survey, aimed at the selection of blazars. From the REX catalogue, we select a well defined and statistically complete sample of high-energy peaked BL Lac (HBL). HBL are expected to be the most numerous class of extragalactic TeV emitting sources. Specifically, we have considered only the REX sources in the currently planned CTA extragalactic survey area satisfying specific criteria and with an optical spectroscopic confirmation. We obtain 46 HBL candidates that we called Te-REX (TeV-emitting REX). We estimate the very high-energy gamma-ray emission, in the TeV domain, using an empirical approach i.e. using specific statistical relations between gamma-rays (at GeV energies) and radio/X-rays properties observed in bright HBL from the literature. We compare the spectral energy distributions (SEDs) with the sensitivities of current and upcoming Cherenkov telescopes and we predict that 14 Te-REX could be detectable with 50 hours of observations of CTA and 7 of them also with current Cherenkov facilities in 50 hours. By extrapolating these numbers on the total extragalactic sky, we predict that about 800 HBL could be visible in pointed CTA observations and about 400 with current Cherenkov telescopes in 50 hours. Interestingly, our predictions show that a non-negligible fraction (about 30%) of the HBL that will be detectable by CTA is composed of relatively weak objects whose optical nuclear emission is swamped by the host-galaxy light and not (yet) detected by Fermi-LAT., Comment: Accepted in MNRAS

Full list of authors: Acciari, V. A.; Ansoldi, S.; Antonelli, L. A.; Arbet Engels, A.; Artero, M.; Asano, K.; Babić, A.; Baquero, A.; Barres de Almeida, U.; Barrio, J. A.; Batković, I.; Becerra González, J.; Bednarek, W.; Bellizzi, L.; Bernardini, E.; Bernardos, M.; Berti, A.; Besenrieder, J.; Bhattacharyya, W.; Bigongiari, C.; Blanch, O.; Bošnjak, Ž.; Busetto, G.; Carosi, R.; Ceribella, G.; Cerruti, M.; Chai, Y.; Chilingarian, A.; Cikota, S.; Colak, S. M.; Colombo, E.; Contreras, J. L.; Cortina, J.; Covino, S.; D'Amico, G.; D'Elia, V.; da Vela, P.; Dazzi, F.; de Angelis, A.; de Lotto, B.; Delfino, M.; Delgado, J.; Delgado Mendez, C.; Depaoli, D.; di Pierro, F.; di Venere, L.; Do Souto Espiñeira, E.; Dominis Prester, D.; Donini, A.; Doro, M.; Fallah Ramazani, V.; Fattorini, A.; Ferrara, G.; Fonseca, M. V.; Font, L.; Fruck, C.; Fukami, S.; García López, R. J.; Garczarczyk, M.; Gasparyan, S.; Gaug, M.; Giglietto, N.; Giordano, F.; Gliwny, P.; Godinović, N.; Green, J. G.; Green, D.; Hadasch, D.; Hahn, A.; Heckmann, L.; Herrera, J.; Hoang, J.; Hrupec, D.; Hütten, M.; Inada, T.; Inoue, S.; Ishio, K.; Iwamura, Y.; Jiménez, I.; Jormanainen, J.; Jouvin, L.; Kajiwara, Y.; Karjalainen, M.; Kerszberg, D.; Kobayashi, Y.; Kubo, H.; Kushida, J.; Lamastra, A.; Lelas, D.; Leone, F.; Lindfors, E.; Lombardi, S.; Longo, F.; López-Coto, R.; López-Moya, M.; López-Oramas, A.; Loporchio, S.; Machado de Oliveira Fraga, B.; Maggio, C.; Majumdar, P.; Makariev, M.; Mallamaci, M.; Maneva, G.; Manganaro, M.; Maraschi, L.; Mariotti, M.; Martínez, M.; Mazin, D.; Menchiari, S.; Mender, S.; Mićanović, S.; Miceli, D.; Miener, T.; Minev, M.; Miranda, J. M.; Mirzoyan, R.; Molina, E.; Moralejo, A.; Morcuende, D.; Moreno, V.; Moretti, E.; Neustroev, V.; Nigro, C.; Nilsson, K.; Nishijima, K.; Noda, K.; Nozaki, S.; Ohtani, Y.; Oka, T.; Otero-Santos, J.; Paiano, S.; Palatiello, M.; Paneque, D.; Paoletti, R.; Paredes, J. M.; Pavletić, L.; Peñil, P.; Perennes, C.; Persic, M.; Prada Moroni, P. G.; Prandini, E.; Priyadarshi, C.; Puljak, I.; Ribó, M.; Rico, J.; Righi, C.; Rugliancich, A.; Saha, L.; Sahakyan, N.; Saito, T.; Sakurai, S.; Satalecka, K.; Saturni, F. G.; Schmidt, K.; Schweizer, T.; Sitarek, J.; Šnidarić, I.; Sobczynska, D.; Spolon, A.; Stamerra, A.; Strom, D.; Strzys, M.; Suda, Y.; Surić, T.; Takahashi, M.; Tavecchio, F.; Temnikov, P.; Terzić, T.; Teshima, M.; Tosti, L.; Truzzi, S.; Tutone, A.; Ubach, S.; van Scherpenberg, J.; Vanzo, G.; Vazquez Acosta, M.; Ventura, S.; Verguilov, V.; Vigorito, C. F.; Vitale, V.; Vovk, I.; Will, M.; Wunderlich, C.; Zarić, D.; Baack, D.; Balbo, M.; Biederbeck, N.; Biland, A.; Bretz, T.; Buss, J.; Dorner, D.; Eisenberger, L.; Elsaesser, D.; Hildebrand, D.; Iotov, R.; Mannheim, K.; Neise, D.; Noethe, M.; Paravac, A.; Rhode, W.; Schleicher, B.; Sliusar, V.; Walter, R.; D'Ammando, F.; Horan, D.; Lien, A. Y.; Baloković, M.; Madejski, G. M.; Perri, M.; Verrecchia, F.; Leto, C.; Lähteenmäki, A.; Tornikoski, M.; Ramakrishnan, V.; Järvelä, E.; Vera, R. J. C.; Chamani, W.; Villata, M.; Raiteri, C. M.; Gupta, A. C.; Pandey, A.; Fuentes, A.; Agudo, I.; Casadio, C.; Semkov, E.; Ibryamov, S.; Marchini, A.; Bachev, R.; Strigachev, A.; Ovcharov, E.; Bozhilov, V.; Valcheva, A.; Zaharieva, E.; Damljanovic, G.; Vince, O.; Larionov, V. M.; Borman, G. A.; Grishina, T. S.; Hagen-Thorn, V. A.; Kopatskaya, E. N.; Larionova, E. G.; Larionova, L. V.; Morozova, D. A.; Nikiforova, A. A.; Savchenko, S. S.; Troitskiy, I. S.; Troitskaya, Y. V.; Vasilyev, A. A.; Merkulova, O. A.; Chen, W. P; Samal, M.; Lin, H. C.; Moody, J. W.; Sadun, A. C.; Jorstad, S. G.; Marscher, A. P.; Weaver, Z. R.; Feige, M.; Kania, J.; Kopp, M.; Kunkel, L.; Reinhart, D.; Scherbantin, A.; Schneider, L.; Lorey, C.; Acosta-Pulido, J. A.; Carnerero, M. I.; Carosati, D.; Kurtanidze, S. O.; Kurtanidze, O. M.; Nikolashvili, M. G.; Chigladze, R. A.; Ivanidze, R. Z.; Kimeridze, G. N.; Sigua, L. A.; Joner, M. D.; Spencer, M.; Giroletti, M.; Marchili, N.; Righini, S.; Rizzi, N.; Bonnoli, G.; MAGIC Collaboration; Fact Collaboration.-- This is an Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Aims. We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with a special focus on the multi-band flux correlations. Methods. The dataset has been collected through an extensive multi-wavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina, and Metsahovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and NuSTAR observations allowed a precise measurement of the falling segments of the two spectral components. Results. The very-high-energy (VHE; E 100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow dierent patterns, partly due to substantial changes in the Compton dominance for a few days without a simultaneous increase in the X-ray flux (i.e., orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including NuSTAR observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. The peak frequency of the synchrotron bump varies by two orders of magnitude throughout the campaign. Our multi-band correlation study also hints at an anti-correlation between UV-optical and X-ray at a significance higher than 3. A VHE flare observed on MJD 57788 (2017 February 4) shows gamma-ray variability on multi-hour timescales, with a factor ten increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario.We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors, from 0 min = 2104 to 0 max = 6105. © 2021 Georg Thieme Verlag. All rights reserved., The MAGIC Collaboration would like to thank the Instituto de Astrofísica de Canarias for the excellent working conditions at the Observatorio del Roque de los Muchachos in La Palma. The financial support of the German BMBF, MPG and HGF; the Italian INFN and INAF; the Swiss National Fund SNF; the ERDF under the Spanish Ministerio de Ciencia e Innovación (MICINN) (FPA2017-87859-P, FPA2017-85668-P, FPA2017-82729-C6-5-R, FPA2017-90566-REDC, PID2019-104114RB-C31, PID2019-104114RB-C32, PID2019-105510GB-C31,PID2019-107847RB-C41, PID2019-107847RB-C42, PID2019-107847RB-C44, PID2019-107988GB-C22); the Indian Department of Atomic Energy; the Japanese ICRR, the University of Tokyo, JSPS, and MEXT; the Bulgarian Ministry of Education and Science, National RI Roadmap Project DO1-268/16.12.2019 and the Academy of Finland grant nr. 320045 is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia “Severo Ochoa” SEV-2016-0588, SEV-2017-0709 and CEX2019-000920-S, and “María de Maeztu” CEX2019-000918-M, the Unidad de Excelencia “María de Maeztu” MDM-2015-0509-18-2 and the “la Caixa” Foundation (fellowship LCF/BQ/PI18/11630012) and by the CERCA program of the Generalitat de Catalunya; by the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project 13.12.1.3.02; by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3; the Polish National Research Centre grant UMO-2016/22/M/ST9/00382; and by the Brazilian MCTIC, CNPq and FAPERJ. The important contributions from ETH Zurich grants ETH-10.08-2 and ETH-27.12-1 as well as the funding by the Swiss SNF and the German BMBF (Verbundforschung Astro- und Astroteilchenphysik) and HAP (Helmoltz Alliance for Astroparticle Physics) are gratefully acknowledged. Part of this work is supported by Deutsche Forschungsgemeinschaft (DFG) within the Collaborative Research Center SFB 876 “Providing Information by Resource-Constrained Analysis”, project C3. We are thankful for the very valuable contributions from E. Lorenz, D. Renker and G. Viertel during the early phase of the project. We thank the Instituto de Astrofísica de Canarias for allowing us to operate the telescope at the Observatorio del Roque de los Muchachos in La Palma, the Max-Planck-Institut für Physik for providing us with the mount of the former HEGRA CT3 telescope, and the MAGIC collaboration for their support. The Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat à l’Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucléaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d’Études Spatiales in France. This work performed in part under DOE Contract DE-AC02-76SF00515. This work made use of data from the NuSTAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by the National Aeronautics and Space Administration. We thank the NuSTAR Operations, Software, and Calibration teams for support with the execution and analysis of these observations. This research has made use of the NuSTAR Data Analysis Software (NuSTARDAS) jointly developed by the ASI Science Data Center (ASDC; Italy) and the California Institute of Technology (USA). This research has also made use of the XRT Data Analysis Software (XRTDAS) developed under the responsibility of the ASI Science Data Center (ASDC), Italy. A.A.E and D.P acknowledge support from the Deutsche Forschungs gemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2094 – 390783311. M. B. acknowledges support from the YCAA Prize Postdoctoral Fellowship and from the Black Hole Initiative at Harvard University, which is funded in part by the Gordon and Betty Moore Foundation (grant GBMF8273) and in part by the John Templeton Foundation. This publication makes use of data obtained at the Metsähovi Radio Observatory, operated by Aalto University in Finland. This research has made use of data from the OVRO 40-m monitoring program (Richards et al. 2011) which is supported in part by NASA grants NNX08AW31G, NNX11A043G, and NNX14AQ89G and NSF grants AST-0808050 and AST-1109911. I.A. acknowledges financial support from the Spanish “Ministerio de Ciencia e Innovación” (MCINN) through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía-CSIC (SEV-2017-0709). Acquisition and reduction of the MAPCAT data was supported in part by MICINN through grants AYA2016-80889-P and PID2019-107847RB-C44. The MAPCAT observations were carried out at the German-Spanish Calar Alto Observatory, which is jointly operated by Junta de Andalucía and Consejo Superior de Investigaciones Científicas. C.C. acknowledges support from the European Research Council (ERC) under the European Union Horizon 2020 research and innovation program under the grant agreement No 771282. This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grants KP-06-H28/3 (2018), KP-06-H38/4 (2019) and KP-06-KITAJ/2 (2020). We acknowledge support by Bulgarian National Science Fund under grant DN18-10/2017 and National RI Roadmap Projects DO1-277/16.12.2019 and DO1-268/16.12.2019 of the Ministry of Education and Science of the Republic of Bulgaria. This research was supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia (contract No 451-03-68/2020-14/200002). G.D. acknowledges observing grant support from the Institute of Astronomy and Rozhen NAO BAS through the bilateral joint research project “Gaia Celestial Reference Frame (CRF) and fast variable astronomical objects” (2020–2022, head – G. Damljanovic). The BU group was supported in part by NASA Fermi guest investigator program grants 80NSSC19K1505 and 80NSSC20K1566. This study was based in part on observations conducted using the 1.8 m Perkins Telescope Observatory (PTO) in Arizona, which is owned and operated by Boston University. This article is partly based on observations made with the LCOGT Telescopes, one of whose nodes is located at the Observatorios de Canarias del IAC on the island of Tenerife in the Observatorio del Teide. This article is also based partly on data obtained with the STELLA robotic telescopes in Tenerife, an AIP facility jointly operated by AIP and IAC. The Abastumani team acknowledges financial support by the Shota Rustaveli National Science Foundation under contract FR-19-6174. Based on observations with the Medicina telescope operated by INAF – Istituto di Radioastronomia.

At optical wavelengths, blazar Electric Vector Position Angle (EVPA) rotations linked with gamma-ray activity have been the subject of intense interest and systematic investigation for over a decade. One difficulty in the interpretation of EVPA rotations is the inherent 180{\deg} ambiguity in the measurements. It is therefore essential, when studying EVPA rotations, to ensure that the typical time-interval between successive observations -- i.e. the cadence -- is short enough to ensure that the correct modulo 180{\deg} value is selected. This optimal cadence depends on the maximum intrinsic EVPA rotation speed in blazars, which is currently not known. In this paper we address the following questions for the RoboPol sample: What range of rotation speeds for rotations greater than 90{\deg} can we expect? What observation cadence is required to detect such rotations? Have rapid rotations been missed in EVPA rotation studies thus far? What fraction of data is affected by the ambiguity? And how likely are detected rotations affected by the ambiguity? We answer these questions with three seasons of optical polarimetric observations of a statistical sample of blazars sampled weekly with the RoboPol instrument and an additional season with daily observations. We model the distribution of EVPA changes on time scales from 1-30 days and estimate the fraction of changes exceeding 90{\deg}. We show that at least daily observations are necessary to measure >96% of optical EVPA variability in the RoboPol sample of blazars correctly and that intra-day observations are needed to measure the fastest rotations that have been seen thus far., Comment: accepted to MNRAS, 20 pages, 15 figures

Full list of authors: Acciari, V. A.; Ansoldi, S.; Antonelli, L. A.; Arbet Engels, A. ; Artero, M.; Asano, K.; Baack, D.; Babić, A.; Baquero, A.; Barres de Almeida, U.; Barrio, J. A.; Batković, I.; Becerra González, J.; Bednarek, W.; Bellizzi, L.; Bernardini, E.; Bernardos, M.; Berti, A.; Besenrieder, J.; Bhattacharyya, W.; Bigongiari, C.; Biland, A.; Blanch, O.; Bošnjak, Ž.; Busetto, G.; Carosi, R.; Ceribella, G.; Cerruti, M.; Chai, Y.; Chilingarian, A.; Cikota, S.; Colak, S. M.; Colombo, E.; Contreras, J. L.; Cortina, J.; Covino, S.; D'Amico, G.; D'Elia, V.; da Vela, P.; Dazzi, F.; de Angelis, A.; de Lotto, B.; Delfino, M.; Delgado, J.; Delgado Mendez, C.; Depaoli, D.; Pierro, F. Di; Venere, L. Di; Do Souto Espiñeira, E.; Dominis Prester, D.; Donini, A.; Dorner, D.; Doro, M.; Elsaesser, D.; Fallah Ramazani, V.; Fattorini, A.; Ferrara, G.; Fonseca, M. V.; Font, L.; Fruck, C.; Fukami, S.; García López, R. J.; Garczarczyk, M.; Gasparyan, S.; Gaug, M.; Giglietto, N.; Giordano, F.; Gliwny, P.; Godinović, N.; Green, J. G.; Green, D.; Hadasch, D.; Hahn, A.; Heckmann, L.; Herrera, J.; Hoang, J.; Hrupec, D.; Hütten, M.; Inada, T.; Inoue, S.; Ishio, K.; Iwamura, Y.; Jiménez, I.; Jormanainen, J.; Jouvin, L.; Kajiwara, Y.; Karjalainen, M.; Kerszberg, D.; Kobayashi, Y.; Kubo, H.; Kushida, J.; Lamastra, A.; Lelas, D.; Leone, F.; Lindfors, E.; Lombardi, S.; Longo, F.; López-Coto, R.; López-Moya, M.; López-Oramas, A.; Loporchio, S.; Machado de Oliveira Fraga, B.; Maggio, C.; Majumdar, P.; Makariev, M.; Mallamaci, M.; Maneva, G.; Manganaro, M.; Mannheim, K.; Maraschi, L.; Mariotti, M.; Martínez, M.; Mazin, D.; Menchiari, S.; Mender, S.; Mićanović, S.; Miceli, D.; Miener, T.; Minev, M.; Miranda, J. M.; Mirzoyan, R.; Molina, E.; Moralejo, A.; Morcuende, D.; Moreno, V.; Moretti, E.; Neustroev, V.; Nigro, C.; Nilsson, K.; Nishijima, K.; Noda, K.; Nozaki, S.; Ohtani, Y.; Oka, T.; Otero-Santos, J.; Paiano, S.; Palatiello, M.; Paneque, D.; Paoletti, R.; Paredes, J. M.; Pavletić, L.; Peñil, P.; Perennes, C.; Persic, M.; Prada Moroni, P. G.; Prandini, E.; Priyadarshi, C.; Puljak, I.; Rhode, W.; Ribó, M.; Rico, J.; Righi, C.; Rugliancich, A.; Saha, L.; Sahakyan, N.; Saito, T.; Sakurai, S.; Satalecka, K.; Saturni, F. G.; Schleicher, B.; Schmidt, K.; Schweizer, T.; Sitarek, J.; Šnidarić, I.; Sobczynska, D.; Spolon, A.; Stamerra, A.; Strom, D.; Strzys, M.; Suda, Y.; Surić, T.; Takahashi, M.; Tavecchio, F.; Temnikov, P.; Terzić, T.; Teshima, M.; Tosti, L.; Truzzi, S.; Tutone, A.; Ubach, S.; van Scherpenberg, J.; Vanzo, G.; Vazquez Acosta, M.; Ventura, S.; Verguilov, V.; Vigorito, C. F.; Vitale, V.; Vovk, I.; Will, M.; Wunderlich, C.; Zarić, D.; Bissaldi, E.; Bonnoli, G.; Cutini, S.; D'Ammando, F.; Nabizadeh, A.; Marchini, A.; Orienti, M.; MAGIC Collaboration., We report here on the first multiwavelength (MWL) campaign on the blazar TXS 1515-273, undertaken in 2019 and extending from radio to very-high-energy gamma-rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotron-peaked source. MAGIC observations resulted in the first-time detection of the source in VHE with a statistical significance of 7.6σ. The average integral VHE flux of the source is 6 ± 1 per cent of the Crab nebula flux above 400 GeV. X-ray coverage was provided by Swift-XRT, XMM-Newton, and NuSTAR. The long continuous X-ray observations were separated by ∼9 h, both showing clear hour scale flares. In the XMM-Newton data, both the rise and decay time-scales are longer in the soft X-ray than in the hard X-ray band, indicating the presence of a particle cooling regime. The X-ray variability time-scales were used to constrain the size of the emission region and the strength of the magnetic field. The data allowed us to determine the synchrotron peak frequency and classify the source as a flaring high, but not extreme synchrotron-peaked object. Considering the constraints and variability patterns from the X-ray data, we model the broad-band spectral energy distribution. We applied a simple one-zone model, which could not reproduce the radio emission and the shape of the optical emission, and a two-component leptonic model with two interacting components, enabling us to reproduce the emission from radio to VHE band. © 2021 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., The financial support of the German BMBF, MPG, and HGF; the Italian INFN and INAF; the Swiss National Fund SNF; the ERDF under the Spanish Ministerio de Ciencia e Innovacion (MICINN) (FPA201787859-P, FPA2017-85668-P, FPA2017-82729-C6-5-R, FPA201790566-REDC, PID2019-104114RB-C31, PID2019-104114RB-C32, PID2019-105510GB-C31, PID2019-107847RB-C41, PID2019107847RB-C42, PID2019-107847RB-C44, PID2019-107988GBC22); the Indian Department of Atomic Energy; the Japanese ICRR, the University of Tokyo, JSPS, and MEXT; the Bulgarian Ministry of Education and Science, National RI Roadmap Project DO1268/16.12.2019 and the Academy of Finland grant nr. 317637 and 320045 are gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia `Severo Ochoa' SEV-20160588, SEV-2017-0709 and CEX2019-000920-S, and `Maria de Maeztu' CEX2019-000918-M, the Unidad de Excelencia `Maria de Maeztu' MDM-2015-0509-18-2 and the `la Caixa' Foundation (fellowshipLCF/BQ/PI18/11630012), by the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project 13.12.1.3.02, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3, the Polish National Research Centre grant UMO-2016/22/M/ST9/00382 and by the Brazilian MCTIC, CNPq, and FAPERJ. The Fermi-LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, theCommissariat à l'Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucleaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase from the following agencies is also gratefully acknowledged: the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Etudes Spatiales in France.

We study the connection between the X-ray and UV properties of the broad absorption line (BAL) wind in the highly X-ray variable quasar PG 2112+059 by comparing Chandra-ACIS data with contemporaneous UV HST/STIS spectra in three different epochs. We observe a correlation whereby an increase in the equivalent-widths (EWs) of the BALs is accompanied by a redder UV spectrum. The growth in the BALs EWs is also accompanied by a significant dimming in soft X-ray emission (2 keV) are only accompanied by minor spectral variations of the UV-BALs and do not show significant changes in the EW of BALs. These trends suggest a wind-shield scenario where the outflow inclination with respect to the line of sight is decreasing and/or the wind mass is increasing. These changes elevate the covering fraction and/or column densities of the BALs and are likely accompanied by a nearly contemporaneous increase in the column density of the shield., 15 pages, 9 figures, accepted in MNRAS

We present multi-wavelength spectral and temporal variability analysis of PKS 0027-426 using optical griz observations from DES (Dark Energy Survey) between 2013-2018 and VOILETTE (VEILS Optical Light curves of Extragalactic TransienT Events) between 2018-2019 and near infrared (NIR) JKs observations from VEILS (VISTAExtragalactic Infrared Legacy Survey) between 2017-2019. Multiple methods of cross-correlation of each combination of light curve provides measurements of possible lags between optical-optical, optical-NIR, and NIR-NIR emission, for each observation season and for the entire observational period. Inter-band time lag measurements consistently suggest either simultaneous emission or delays between emission regions on timescales smaller than the cadences of observations. The colour-magnitude relation between each combination of filters was also studied to determine the spectral behaviour of PKS 0027-426. Our results demonstrate complex colour behaviour that changes between bluer when brighter (BWB), stable when brighter (SWB) and redder when brighter (RWB) trends over different timescales and using different combinations of optical filters. Additional analysis of the optical spectra is performed to provide further understanding of this complex spectral behaviour., Accepted for publication in MNRAS. 33 pages, 28 figures, 11 tables

Supermassive black holes launch highly relativistic jets with velocities reaching Lorentz factors as high as $\Gamma>50$. How the jets accelerate to such high velocities and where along the jet do they reach terminal velocity are open questions that are tightly linked to their structure, launching and dissipation mechanisms. Changes in the beaming factor along the jets could potentially reveal jet acceleration, deceleration, or bending. We aim to (1) quantify the relativistic effects in multiple radio frequencies and (2) study possible jet velocity--viewing angle variations at parsec scales. We used the state-of-the-art code Magnetron to model light curves from the University of Michigan Radio Observatory and the Mets\"{a}hovi Radio Observatory's monitoring programs in five frequencies covering about 25 years of observations in the 4.8-37~GHz range for 61 sources. We supplement our data set with high-frequency radio observations in the 100-340~GHz range from ALMA, CARMA, and SMA. For each frequency we estimate the Doppler factor which we use to quantify possible changes in the relativistic effects along the jets. The majority of our sources do not show any statistically significant difference in their Doppler factor across frequencies. This is consistent with constant velocity in a conical jet, as expected at parsec scales. However, our analysis reveals 17 sources where relativistic beaming changes as a function of frequency. In the majority of cases the Doppler factor increases towards lower frequencies. Only 1253-053 shows the opposite behavior. By exploring their jet properties we find that the jet of 0420-014 is likely bent across the 4.8-340~GHz range. For 0212+735 the jet is likely parabolic, and still accelerating in the 4.8-37~GHz range. We discuss possible interpretations for the trends found in the remaining sources., Comment: 15 pages, 5 figure, 3 Tables, accepted for publication in A&A

For understanding the diversity of jetted active galactic nuclei (AGN) and especially the puzzling wide range in their radio-loudness, it is important to understand what role the magnetic fields play in setting the power of relativistic jets in AGN. We have performed multi-frequency (4-24 GHz) VLBA phase-referencing observations of the radio-intermediate quasar III Zw 2 using three nearby calibrators as reference sources to estimate jet magnetic flux by measuring the core-shift effect. By combining the self-referencing core-shift of each calibrator with the phase-referencing core-shifts, we obtained an upper limit of 0.16 mas for the core-shift between 4 and 24 GHz in III Zw 2. By assuming equipartition between magnetic and particle energy densities and adopting the flux-freezing approximation, we further estimated the upper limit for both magnetic field strength and poloidal magnetic flux threading the black hole. We find that the upper limit to the measured magnetic flux is smaller by at least a factor of five compared to the value predicted by the magnetically arrested disk (MAD) model. An alternative way to derive the jet magnetic field strength from the turnover of the synchrotron spectrum leads to an even smaller upper limit. Hence, the central engine of III Zw 2 has not reached the MAD state, which could explain why it has failed to develop a powerful jet, even though the source harbours a fast-spinning black hole. However, it generates an intermittent jet, which is possibly triggered by small scale magnetic field fluctuations as predicted by the magnetic flux paradigm of Sikora & Begelman (2013). We propose here that combining black hole spin measurements with magnetic field measurements from the VLBI core-shift observations of AGN over a range of jet powers could provide a strong test for the dominant factor setting the jet power relative to the accretion power available., 29 pages, 22 figures, 3 Tables, accepted for publication with A&A

Ricci, C., et al., The merger of two or more galaxies can enhance the inflow of material from galactic scales into the close environments of active galactic nuclei (AGNs), obscuring and feeding the supermassive black hole (SMBH). Both recent simulations and observations of AGN in mergers have confirmed that mergers are related to strong nuclear obscuration. However, it is still unclear how AGN obscuration evolves in the last phases of the merger process. We study a sample of 60 luminous and ultra-luminous IR galaxies (U/LIRGs) from the GOALS sample observed by NuSTAR. We find that the fraction of AGNs that are Compton thick (CT; NH ≥ 1024 cm-2) peaks at 74+14 −19 per cent at a late merger stage, prior to coalescence, when the nuclei have projected separations (dsep) of 0.4-6 kpc. A similar peak is also observed in the median NH [(1.6 ± 0.5) × 1024 cm−2]. The vast majority (85+7 −9 per cent) of the AGNs in the final merger stages (dsep ≤ 10 kpc) are heavily obscured (NH ≥ 1023 cm−2), and the median NH of the accreting SMBHs in our sample is systematically higher than that of local hard X-ray-selected AGN, regardless of the merger stage. This implies that these objects have very obscured nuclear environments, with the NH ≥ 1023 cm−2 gas almost completely covering the AGN in late mergers. CT AGNs tend to have systematically higher absorption-corrected X-ray luminosities than less obscured sources. This could either be due to an evolutionary effect, with more obscured sources accreting more rapidly because they have more gas available in their surroundings, or to a selection bias. The latter scenario would imply that we are still missing a large fraction of heavily obscured, lower luminosity (L2−10 ≤ 1043 erg s−1) AGNs in U/LIRGs., LCH was supported by the National Key Research and Development Program of China (2016YFA0400702) and the National Science Foundation of China (11721303 and 11991052). CR acknowledges support from the Fondecyt Iniciacion grant 11190831. ET acknowledges support from CATA-Basal AFB-170002, FONDECYT Regular grant 1190818, ANID Anillo ACT172033, and Millennium Nucleus NCN19_058 (TITANs). FEB acknowledges support from ANID – Millennium Science Initiative Program – ICN12_009, CATA-Basal – AFB-170002, and FONDECYT Regular – 1190818 and 1200495. SA gratefully acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 789410). VU acknowledges funding support from NASA Astrophysics Data Analysis Program (ADAP) Grant 80NSSC20K0450. AMM acknowledges support from the National Science Foundation under grant number 2009416. KI acknowledges support by the Spanish MICINN under grant Proyecto/AEI/10.13039/501100011033 and ‘Unit of excellence María de Maeztu 2020-2023’ awarded to ICCUB (CEX2019-000918-M). PA acknowledges financial support from ANID Millennium Nucleus NCN19-058 (TITANS) and the Max Planck Society through a Partner Group. HI acknowledges support from JSPS KAKENHI Grant Number JP19K23462.

We analyse optical datacubes of the inner kiloparsec of 30 local ($z\le0.02$) active galactic nuclei (AGN) hosts that our research group, AGNIFS, has collected over the past decade via observations with the integral field units of the Gemini Multi-Object Spectrographs. Spatial resolutions range between $50~{\rm pc}$ and $300~{\rm pc}$ and spectral coverage is from $4800~\mathring{A}$ or $5600~\mathring{A}$ to $7000~\mathring{A}$, at velocity resolutions of $\approx 50~{\rm km~s^{-1}}$. We derive maps of the gas excitation and kinematics, determine the AGN ionisation axis -- which has random orientation relative to the galaxy, and the kinematic major axes of the emitting gas. We find that rotation dominates the gas kinematics in most cases, but is disturbed by the presence of inflows and outflows. Outflows have been found in 21 nuclei, usually along the ionisation axis. The gas velocity dispersion is traced by $W_{80}$ (velocity width encompassing 80 per cent of the line flux), adopted as a tracer of outflows. In 7 sources $W_{80}$ is enhanced perpendicularly to the ionisation axis, indicating lateral expansion of the outflow. We have estimated mass-outflow rates $\dot{M}$ and powers $\dot{E}$, finding median values of $\log\,[\dot{M}/({\rm\,M_\odot\,yr^{-1}})]=-2.1_{-1.0}^{+1.6}$ and $\log\,[\dot{E}/({\rm\,erg\,s^{-1}})]=38.5_{-0.9}^{+1.8}$, respectively. Both quantities show a mild correlation with the AGN luminosity ($L_{\rm AGN}$). $\dot{E}$ is of the order of 0.01 $L_{\rm AGN}$ for 4 sources, but much lower for the majority (9) of the sources, with a median value of $\log\,[\dot{E}/L_{\rm AGN}]=-5.34_{-0.9}^{+3.2}$ indicating that typical outflows in the local Universe are unlikely to significantly impact their host galaxy evolution.

We present 3-mm observations of the quasar J0015+1842 at z=2.63 with the NOrthern Extended Millimeter Array (NOEMA). Our data reveals molecular gas, traced via a Gaussian CO(3-2) line, with a remarkably large velocity dispersion (FWHM=1010+/-120 km/s) and corresponding to a total molecular mass MH2~(3.4-17)x10^10 Msun, depending on the adopted CO-to-H2 conversion factor alphaCO=(0.8-4.0) Msun (km/s pc^2)^-1. Assuming the 3-mm continuum emission is thermal, we derive a dust mass of the order of Mdust ~5x10^8 Msun. J0015+1842 is located in the molecular gas-rich region in the IR vs CO line luminosity diagram, in-between the main locus of main-sequence and sub-millimetre galaxies and that of most other AGNs targeted so far for CO measurements. While the large velocity dispersion of the CO line suggests a merging system, J0015+1842 is observed to be a regular, only very moderately dust-reddened (Av~0.3-0.4) type-I quasar from its UV-optical spectrum, from which we infer a mass of the super-massive black hole be around MBH~6x10^8 Msun. We suggest that J0015+1842 is observed at a galaxy evolutionary stage where a massive merger has brought significant amounts of gas towards an actively accreting super-massive black hole (quasar). While the host still contains a large amount of dust and molecular gas with high velocity dispersion, the quasar has already cleared the way towards the observer, likely through powerful outflows as recently revealed by optical observations of the same object. High angular resolution observations of this and similar systems} should help determining better the respective importance of evolution and orientation in the appearance of quasars and their host galaxies and have the potential to investigate early feedback and star-formation processes in galaxies in their quasar phases., Comment: 5 pages, 5 figures, accepted for publication in A\&A

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

Nuclear Activity in Galaxies Across Cosmic Time, Proceedings of the conference held 7-11 October 2019 in Addis Ababa, Ethiopia. Edited by Mirjana Pović et al. Proceedings of the International Astronomical Union, Volume 356, pp. 12-16, The J-PAS survey will soon start observing thousands of square degrees of the Northern Sky with its unique set of 56 narrow band filters covering the entire optical wavelength range, providing, effectively, a low resolution spectra for every object detected. Active galaxies and quasars, thanks to their strong emission lines, can be easily identified and characterized with J-PAS data. A variety of studies can be performed, from IFU-like analysis of local AGN, to clustering of high-z quasars. We also expect to be able to extract intrinsic physical quasar properties from the J-PAS pseudo-spectra, including continuum slope and emission line luminosities. Here we show the first attempts of using the QSFit software package to derive the properties for 22 quasars at 0.8 < z < 2 observed by the miniJPAS survey, the first deg2 of J-PAS data obtained with an interim camera. Results are compared with the ones obtained by applying the same software to SDSS quasar spectra., 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)

In order to understand the possible mechanisms of recurrent jet activity in radio galaxies and quasars, which are still unclear, we have identified such sources with a large range of linear sizes (220 $-$ 917 kpc), and hence time scales of episodic activity. Here we present high-sensitivity 607-MHz Giant Metrewave Radio Telescope (GMRT) images of 21 possible double-double radio galaxies (DDRGs) identified from the FIRST survey to confirm their episodic nature. These GMRT observations show that none of the inner compact components suspected to be hot-spots of the inner doubles are cores having a flat radio spectrum, confirming the episodic nature of these radio sources. We have indentified a new DDRG with a candidate quasar, and have estimated the upper spectral age limits for eight sources which showed marginal evidence of steepening at higher frequencies. The estimated age limits (11 $-$ 52 Myr) are smaller than those of the large-sized ($\sim$ 1 Mpc) DDRGs., Comment: Accepted for publication in MNRAS. 14 pages, 7 figures

We present the molecular gas morphology and kinematics of seven nearby Seyfert galaxies obtained from our 230~GHz ALMA observations. The CO J=2-1 kinematics within the inner $\sim30$" ($\lesssim9$~kpc) reveals rotation patterns that have been explored using the Bertola rotation model and a modified version of the Kinemetry package. The latter algorithm reveals various deviations from pure circular rotation in the inner kiloparsec of all seven galaxies, including kinematic twists, decoupled and counter-rotating cores. A comparison of the global molecular gas and stellar kinematics show overall agreement in the position angle of the major axis and the systemic velocity, but larger discrepancies in the disc inclination. The residual maps obtained with both the methods shows the presence of non-circular motions in most of the galaxies. Despite its importance, a detailed interpretation of the physics responsible for non-circular motions will be discussed in a forthcoming work., Comment: 28 pages in total, 20 figures, 4 tables. Accepted for publication in MNRAS

We present a pilot study aimed to identify large-scale galaxy properties that could play a role in activating a quiescent nucleus. To do so, we compare the properties of two isolated nearby active galaxies and their non-active twins selected from the Calar Alto Legacy Integral Field Area (CALIFA) survey. This pilot sample includes two barred and two unbarred galaxies. We characterize the stellar and ionized gas kinematics and also their stellar content. We obtain simple kinematic models by fitting the full stellar and ionized gas velocity fields and just the approaching or receding sides. We find that the analysed active galaxies present lopsided discs and higher values of the global stellar angular momentum (λR) than their non-active twins. This could be indicating that the stellar discs of the AGN gained angular momentum from the inflowing gas that triggered the nuclear activity. The inflow of gas could have been produced by a twisted disc instability in the case of the unbarred AGN, and by the bar in the case of the barred AGN. In addition, we find that the central regions of the studied active galaxies show older stellar populations than their non-active twins. The next step is to statistically explore these galaxy properties in a larger sample of twin galaxies. © 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society., This study uses data provided by the Calar Alto Legacy Integral Field Area (CALIFA) survey (http://califa.caha.es/).Based on observations collected at the Centro Astronomico Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max-Planck-Institut fur Astronomie and the Instituto de Astrofisica de Andalucia (CSIC). IMC acknowledges the support of the Instituto de Astrofisica de Canarias via an Astrophysicist Resident fellowship. BGL acknowledges support from the Spanish Ministerio de Economia y Competitividad (MINECO) by the grant AYA2015-68217-P. CRA acknowledges the Ramon y Cajal Program of the Spanish Ministry of Economy and Competitiveness through project RYC-2014-15779 and the Spanish Plan Nacional de Astronomia y Astrofisica under grant AYA2016-76682-C3-2-P. TRL acknowledges support via grantsAYA2014-56795-P and AYA2016-77237-C3-1-P from the Spanish Government. JFB acknowledges support from the Spanish Ministerio de Economia y Competitividad (MINECO) by the grant AYA2013-48226-C3-1P. SFS thanks the CONACyT programs CB-285080 and DGAPA IA101217 grants for their support. IMP and JM acknowledges support from the Spanish Ministerio de Economia y Competitividad(MINECO) by the grant AYA2016-76682-C3-1-P.

We map the co-eval growth of galaxies and their central supermassive black holes in detail by measuring the incidence of Active Galactic Nuclei (AGN) in galaxies as a function of star formation rate (SFR) and redshift (to z~4). We combine large galaxy samples with deep Chandra X-ray imaging to measure the probability distribution of specific black hole accretion rates (LX relative to stellar mass) and derive robust AGN fractions and average specific accretion rates. First, we consider galaxies along the main sequence of star formation. We find a linear correlation between the average SFR and both the AGN fraction and average specific accretion rate across a wide range in stellar mass ($M_* \sim 10^{8.5-11.5}M_\odot$) and to at least z~2.5, indicating that AGN in main-sequence galaxies are driven by the stochastic accretion of cold gas. We also consider quiescent galaxies and find significantly higher AGN fractions than predicted, given their low SFRs, indicating that AGN in quiescent galaxies are fuelled by additional mechanisms (e.g. stellar winds). Next, we bin galaxies according to their SFRs relative to the main sequence. We find that the AGN fraction is significantly elevated for galaxies that are still star-forming but with SFRs below the main sequence, indicating further triggering mechanisms enhance AGN activity within these sub-main-sequence galaxies. We also find that the incidence of high-accretion-rate AGN is enhanced in starburst galaxies and evolves more mildly with redshift than within the rest of the galaxy population, suggesting mergers play a role in driving AGN activity in such high-SFR galaxies., Comment: 19 pages, 15 figures. This is a pre-copyedited, author-produced version of an article accepted for publication in MNRAS following peer review

We quantify the luminosity contribution of active galactic nuclei (AGN) to the 12 $\mu$m, mid-infrared (MIR; 5-38 $\mu$m), and the total IR (5-1000 $\mu$m) emission in the local AGN detected in the all-sky 70-month Swift/Burst Alert Telescope (BAT) ultra hard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into AGN and starburst components using AGN torus and star-forming galaxy templates. This enables us to recover the AGN torus emission also for low-luminosity end, down to $\log (L_{14-150}/{\rm erg}~{\rm s}^{-1}) \simeq 41$, which typically have significant host galaxy contamination. We find that the luminosity contribution of the AGN to the 12 $\mu$m, the MIR, and the total IR band is an increasing function of the 14-150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 $\mu$m. The obtained total IR AGN luminosity through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median of the dust covering factor is always smaller than that of the X-ray obscuration fraction above the AGN bolometric luminosity of $\log (L_{\rm bol}/{\rm erg}~{\rm s}^{-1}) \simeq 42.5$. Considering that X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, it indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission., Comment: 21 pages, 15 figures, accepted for publication in ApJ. The full list of Table 1 is available at http://www.kusastro.kyoto-u.ac.jp/~ichikawa/Table1_MR_20181107.txt

Fundamental differences in the radio properties of red quasars (QSOs), as compared to blue QSOs, have been recently discovered, positioning them as a potential key population in the evolution of galaxies and black holes across cosmic time. To elucidate their nature, we exploited a rich compilation of photometry and spectroscopic data to model their spectral energy distributions (SEDs) from the UV to the FIR and characterise their emission-line properties. Following a systematic comparison approach, we infer the AGN accretion, obscuration, and host galaxy properties in a sample of ~1800 QSOs at 0.21000 km/s) in red QSOs. We find that red QSOs that exhibit evidence for high-velocity winds present a stronger signature of the infrared excess, suggesting a causal connection between reddening and the presence of hot dust in QSO winds. We propose that dusty winds at nuclear scales are potentially the physical ingredient responsible for the colours in red QSOs, as well as a key parameter for the regulation of accretion material in the nucleus., Comment: Accepted for publication in Astronomy & Astrophysics. 22 pages, 16 figures. Abstract abridged

Context. Many low-redshift active galactic nuclei harbor a supermassive black hole accreting matter at low or medium rates. At such rates, the accretion flow usually consists of a cold optically thick disk, plus a hot, low density, collisionless corona. In the latter component, charged particles can be accelerated to high energies by various mechanisms. Aims. We aim to investigate, in detail, nonthermal processes in hot accretion flows onto supermassive black holes, covering a wide range of accretion rates and luminosities. Methods. We developed a model consisting of a thin Shakura-Sunyaev disk plus an inner hot accretion flow or corona, modeled as a radiatively inefficient accretion flow, where nonthermal processes take place. We solved the transport equations for relativistic particles and estimated the spectral energy distributions resulting from nonthermal interactions between the various particle species and the fields in the source. Results. We covered a variety of scenarios, from low accretion rates up to 10% of the Eddington limit, and identified the relevant cooling mechanisms in each case. The presence of hadrons in the hot flow is decisive for the spectral shape, giving rise to secondary particles and gamma-ray cascades. We applied our model to the source IC 4329A, confirming earlier results that showed evidence of nonthermal particles in the corona., 16 pages, 13 figures, matching published version in A&A

arXiv:2103.03584v1, We use the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) project as a template to demonstrate how deep spectrophotometric surveys covering large cosmological volumes over extended fields (1–15 deg2) with a mid-IR imaging spectrometer (17–36 µm) in conjunction with deep photometry with a far-IR camera, at wavelengths which are not affected by dust extinction can answer the most crucial questions in current galaxy evolution studies. A SPICA-like mission will be able for the first time to provide an unobscured three-dimensional (3D, i.e. x, y, and redshift z) view of galaxy evolution back to an age of the universe of less than ∼ 2 Gyrs, in the mid-IR rest frame. This survey strategy will produce a full census of the Star Formation Rate (SFR) in the universe, using polycyclic aromatic hydrocarbons (PAH) bands and fine-structure ionic lines, reaching the characteristic knee of the galaxy luminosity function, where the bulk of the population is distributed, at any redshift up to z ∼ 3.5. Deep follow-up pointed spectroscopic observations with grating spectrometers onboard the satellite, across the full IR spectral range (17–210 µm), would simultaneously measure Black Hole Accretion Rate (BHAR), from high-ionisation fine-structure lines, and SFR, from PAH and low- to mid-ionisation lines in thousands of galaxies from solar to low metallicities, down to the knee of their luminosity functions. The analysis of the resulting atlas of IR spectra will reveal the physical processes at play in evolving galaxies across cosmic time, especially its heavily dust-embedded phase during the activity peak at the cosmic noon (z ∼ 1–3), through IR emission lines and features that are insensitive to the dust obscuration., LS and JAFO acknowledge financial support by the Agenzia Spaziale Italiana (ASI) under the research contract 2018-31-HH.0. AAH acknowledges support from grants PGC2018-094671-B-I00 (MCIU/AEI/FEDER,UE) and No. MDM-2017-0737 at Unidad de Excelencia María de Maeztu-Centro de Astrobiología (INTACSIC). FJC acknowledges financial support from the Spanish Ministry MCIU under project RTI2018-096686-B-C21 (MCIU/AEI/FEDER/UE), cofunded by FEDER funds and from the Agencia Estatal de Investigación, Unidad de Excelencia María de Maeztu, ref. MDM-2017-0765. HD acknowledges financial support from the Spanish Ministry of Science, Innovation and Universities (MICIU) under the 2014 Ramón y Cajal program RYC-2014-15686 and AYA2017-84061-P, the later one co-financed by FEDER (European Regional Development Funds). AL acknowledges the support from Comunidad de Madrid through the Atracción de Talento grant 2017-T1/TIC-5213. MPS acknowledges support from the Comunidad de Madrid, Spain, through Atracción de Talento Investigador Grant 2018-T1/TIC-11035.

Full list of authors: Cui, Yu-Zhu; Hada, Kazuhiro; Kino, Motoki; Sohn, Bong-Won; Park, Jongho; Ro, Hyun-Wook; Sawada-Satoh, Satoko; Jiang, Wu; Cui, Lang; Honma, Mareki; Shen, Zhi-Qiang; Tazaki, Fumie; An, Tao; Cho, Ilje; Zhao, Guang-Yao; Cheng, Xiao-Peng; Niinuma, Kotaro; Wajima, Kiyoaki; Zhang, Ying-Kang; Kawaguchi, Noriyuki; Algaba, Juan-Carlos; Koyama, Shoko; Hirota, Tomoya; Yonekura, Yoshinori; Sakai, Nobuyuki; Xia, Bo; Jiang, Yong-Bin; Yu, Lin-Feng; Gou, Wei; Hwang, Ju-Yeon; Jiang, Yong-Chen; Sun, Yun-Xia; Jung, Dong-Kyu; Kim, Hyo-Ryoung; Kim, Jeong-Sook; Kobayashi, Hideyuki; Lee, Jee-Won; Lee, Jeong-Ae; Zhang, Hua; Li, Guang-Hui; Xu, Zhi-Qiang; Li, Peng; Oh, Jung-Hwan; Oh, Se-Jin; Oh, Chung-Sik; Oyama, Tomoaki; Roh, Duk-Gyoo; Shibata, Katsunori-M.; Guo, Wen; Zhao, Rong-Bing; Zhong, Wei-Ye; Wang, Jin-Qing; Yang, Wen-Jun; Yan, Hao; Yeom, Jae-Hwan; Li, Bin; Li, Xiao-Fei; Yuan, Jian-Ping; Dong, Jian; Chen, Zhong; Akiyama, Kazunori; Asada, Keiichi; Byun, Do-Young; Hagiwara, Yoshiaki; Hodgson, Jeffrey; Jung, Tae-Hyun; Kim, Kee-Tae; Lee, Sang-Sung; Yi, Kunwoo; Liu, Qing-Hui; Liu, Xiang; Lu, Ru-Sen; Nakamura, Masanori; Trippe, Sascha; Wang, Na; Wang, Xue-Zheng; Zhang, Bo, The East Asian Very Long Baseline Interferometry (VLBI) Network (EAVN) is a rapidly evolving international VLBI array that is currently promoted under joint efforts among China, Japan and Korea. EAVN aims at forming a joint VLBI Network by combining a large number of radio telescopes distributed over East Asian regions. After the combination of the Korean VLBI Network (KVN) and the VLBI Exploration of Radio Astrometry (VERA) into KaVA, further expansion with the joint array in East Asia is actively promoted. Here we report the first imaging results (at 22 and 43 GHz) of bright radio sources obtained with KaVA connected to Tianma 65-m and Nanshan 26-m Radio Telescopes in China. To test the EAVN imaging performance for different sources, we observed four active galactic nuclei (AGN) having different brightness and morphology. As a result, we confirmed that the Tianma 65-m Radio Telescope (TMRT) significantly enhances the overall array sensitivity, a factor of 4 improvement in baseline sensitivity and 2 in image dynamic range compared to the case of KaVA only. The addition of the Nanshan 26-m Radio Telescope (NSRT) further doubled the east-west angular resolution. With the resulting high-dynamic-range, high-resolution images with EAVN (KaVA+TMRT+NSRT), various fine-scale structures in our targets, such as the counter-jet in M87, a kink-like morphology of the 3C 273 jet and the weak emission in other sources are successfully detected. This demonstrates the powerful capability of EAVN to study AGN jets and to achieve other science goals in general. Ongoing expansion of EAVN will further enhance the angular resolution, detection sensitivity and frequency coverage of the network. © 2021 National Astronomical Observatories, CAS and IOP Publishing Ltd., TMRT is operated by Shanghai Astronomical Observatory. NSRT is operated by Xinjiang Astronomical Observatory. KVN is a facility operated at by the Korea Astronomy and Space Science Institute. VERA is a facility operated at National Astronomical Observatory of Japan in collaboration with associated universities in Japan. This work is supported by The Graduate University for Advanced Studies (SOKENDAI). Y.C. is supported by the Japanese Government (MEXT) Scholarship. This work is supported by JSPS KAKENHI Grant Numbers JP18K03656 (M.K.), JP18H03721 (K.N., K.H. and M.K.), JP19H01943 (K.H., F.T. and Y.H.) and JP18KK0090 (K.H. and F.T.). K.H. is supported by the Mitsubishi Foundation (grant number 201911019). J.P. is supported by an EACOA Fellowship awarded by the East Asia Core Observatories Association, which consists of the Academia Sinica Institute of Astronomy and Astrophysics, the National Astronomical Observatory of Japan, the Center for Astronomical Mega-Science, the Chinese Academy of Sciences, and the Korea Astronomy and Space Science Institute. J.P. and I.C. acknowledge the financial support from the National Research Foundation (NRF) of Korea via Global Ph.D. Fellowship Grant 2014H1A2A1018695 and 2015H1A2A1033752, respectively. S.T. acknowledges support from the NRF via Grant 2019R1F1A1059721. This work is supported by the Major Program of the National Natural Science Foundation of China (NSFC, Grant Nos. 11590780 and 11590784), the Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KJCX1-YW-18), the Scientific Program of Shanghai Municipality (08DZ1160100) and Key Laboratory for Radio Astronomy, CAS. W.J. acknowledges support from NSFC (Grant No. 11803071). L.C. is supported by the National Key R&D Program of China (Grant No. 2018YFA0404602), the CAS 'Light of West China' Program (Grant No. 2018-XBQNXZ-B-021), the NSFC (Grant Nos. U2031212 and 61931002) and the Youth Innovation Promotion Association of the CAS (Grant No. 2017084). J.C.A. acknowledges support from Fundamental Research Grant Scheme (FRGS) FRGS/1/2019/STG02/UM/02/6. R.S. Lu is supported by the Max Planck Partner Group of the MPG and the CAS and acknowledges the support by the Key Program of the NSFC (Grant No. 11933007) and the Research Program of Fundamental and Frontier Sciences, CAS (Grant No. ZDBS-LY-SLH011)., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

The role of radio mode feedback in non radio-loud quasars needs to be explored in depth to determine its true importance. Its effects can be identified based on the evidence of interactions between the radio structures and the ambient ionised gas. We investigate this in a sample of 13 optically selected type-2 quasars (QSO2) at z, 32 pages, 39 figures. Accepted for publication in Astronomy & Astrophysics

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon., 19 pages, 9 figures, ApJ in press

This is an open access article under the CC BY license https://creativecommons.org/licenses/by/4.0, Fanaroff-Riley (FR) 0 radio galaxies compose a new class of radio galaxies, which are usually weaker but much more numerous than the well-established class of FR 1 and FR 2 galaxies. The latter classes have been proposed as sources of the ultra-high-energy cosmic rays (UHECRs) with energies reaching up to ∼1020 eV. Based on this conjecture, the possibility of UHECR acceleration and survival in an FR 0 source environment is examined in this work. In doing so, an average spectral energy distribution (SED) based on data from the FR 0 catalog (FR0CAT) is compiled. The resulting photon fields are used as targets for UHECRs, which suffer from electromagnetic pair production, photo-disintegration, photo-meson production losses, and synchrotron radiation. Multiple mechanisms are discussed to assess the UHECR acceleration probability, including Fermi-I order and gradual shear accelerations, and particle escape from the source region. This work shows that in a hybrid scenario, combining Fermi and shear accelerations, FR 0 galaxies can contribute to the observed UHECR flux, as long as Γj≳1.6, where shear acceleration starts to dominate over escape. Even in less optimistic scenarios, FR 0s can be expected to contribute to the cosmic-ray flux between the knee and the ankle. Our results are relatively robust with respect to the realized magnetic turbulence model and the speed of the accelerating shocks. © 2021 The Authors., This work acknowledges financial support from the Austrian Science Fund (FWF) under grant agreement number I 4144-N27 and from the Slovenian Research Agency -ARRS (project no. N10111). MB has for this project received funding from the European Union's Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement No 847476. The views and opinions expressed herein do not necessarily reflect those of the European Commission. G.B. 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). This research has made use of the NASA/IPAC Extragalactic Database (NED), which is funded by the National Aeronautics and Space Administration (NASA) and operated by the California Institute of Technology (CIT). Furthermore, part of this work is based on archival data, software, or online services provided by the Space Science Data Center (SSDC) via ASI.

Penton, A. et al (DES Collaboration), We present the statistical methods that have been developed to analyse the OzDES reverberation mapping sample. To perform this statistical analysis we have created a suite of customizable simulations that mimic the characteristics of each source in the OzDES sample. These characteristics include: the variability in the photometric and spectroscopic light curves, the measurement uncertainties, and the observational cadence. By simulating the sources in the OzDES sample that contain the C iv emission line, we developed a set of criteria that rank the reliability of a recovered time-lag depending on the agreement between different recovery methods, the magnitude of the uncertainties, and the rate at which false positives were found in the simulations. These criteria were applied to simulated light curves and these results used to estimate the quality of the resulting Radius-Luminosity relation. We grade the results using three quality levels (gold, silver, and bronze). The input slope of the R-L relation was recovered within 1σ for each of the three quality samples, with the gold standard having the lowest dispersion with a recovered a R-L relation slope of 0.454 ± 0.016 with an input slope of 0.47. Future work will apply these methods to the entire OzDES sample of 771 AGN., 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 MICINN under grants ESP2017-89838, PGC2018-094773, PGC2018-102021, 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) do e-Universo (CNPq grant 465376/2014-2).

Several evidences indicate that Lyman Break Galaxies (LBGs) in the Epoch of Reionization (redshift z > 6) might host massive black holes (MBHs). We address this question by using a merger-tree model combined with tight constraints from the 7 Ms Chandra survey and the known high-z super-MBH population. We find that a typical LBG with MUV = −22 residing in an Mh ≈ 1012 M⊙ halo at z = 6 host an MBH with mass M• ≈ 2 × 108 M⊙. Depending on the fraction, fseed, of early haloes planted with a direct collapse black hole seed (Mseed = 105M⊙), the model suggests two possible scenarios: (i) if fseed = 1, MBHs in LBGs mostly grow by merging and must accrete at a low (λE ≃ 10−3) Eddington ratio not to exceed the experimental X-ray luminosity upper bound $L_\mathrm{ X}^* = 10^{42.5} {\rm erg\, s}^{-1}$; (ii) if fseed = 0.05, accretion dominates (λE ≃ 0.22) and MBH emission in LBGs must be heavily obscured. In both scenarios the UV luminosity function is largely dominated by stellar emission up to very bright mag, $M_{\rm UV} \lower.5ex\hbox{$\,\, \buildrel\gt \over \sim \,\,$}-23$, with BH emission playing a subdominant role. Scenario (i) poses extremely challenging, and possibly unphysical, requirements on DCBH formation. Scenario (ii) entails testable implications on the physical properties of LBGs involving the FIR luminosity, emission lines, and the presence of outflows.

Feedback from an active galactic nucleus (AGN) is often implicated as a mechanism that leads to the quenching of galactic star formation. However, AGN-driven quenching is challenging to reconcile with observations that AGN hosts tend to harbour equal (or even excess) amounts of gas compared with inactive galaxies of similar stellar mass. In this paper, we investigate whether AGN feedback happens on sub-galactic (kpc) scales, an effect that might be difficult to detect with global gas measurements. Using kpc-scale measurements of molecular gas (ΣH2) and stellar mass (Σ∗) surface densities, taken from the Extragalactic Data base for Galaxy Evolution-Calar Alto Legacy Integral Field Area survey, we show that the gas fractions of central AGN regions are typically a factor of ∼2 lower than those in star-forming regions. Based on four galaxies with the best spaxel statistics, the difference between AGN and star-forming gas fractions is seen even within a given galaxy, indicating that AGN feedback is able to deplete the molecular gas reservoir in the central few kpc. © 2021 The Author(s)., The authors acknowledge an NSERC Discovery Grant (SLE), NSF AST-1616199 (TW), support by the Deutsche Forschungsgemeinschaft project number SFB956A (DC), NSF AST-1615960 (AB, SNV), PID2019-109067GB-I00, P18-FRJ-2595, SEV-2017-0709 (RGB), IA-100420, IN100519, CF19-39578, CB-285080, and FC-2016-01-1916 (SFS and JBB),and Fondecyt grant number 1190684 (MR). Support for CARMA construction was derived from the Gordon and Betty Moore Foundation, the Eileen and Kenneth Norris Foundation, the Caltech Associates, the states of California, Illinois, and Maryland, and the NSF. Funding for CARMA development and operations was supported by NSF and the CARMA partner universities. We acknowledge the usage of the HyperLeda data base (http://leda.univ-lyon1.fr).

We report the relationship between the $Gaia$--VLBI position differences and the magnitudes of source structure effects in VLBI observations. Because the $Gaia$--VLBI position differences are statistically significant for a considerable number of common sources, we attempt to discuss and explain these position differences based on VLBI observations and available source images at cm-wavelengths. Based on the derived closure amplitude root-mean-square (CARMS), which quantifies the magnitudes of source structure effects in the VLBI observations used for building the third realization of the International Celestial Reference Frame, the arc lengths and normalized arc lengths of the position differences are examined in detail. The radio jet directions and the directions of the $Gaia$--VLBI position differences are investigated for a small sample of sources. Both the arc lengths and normalized arc lengths of the $Gaia$ and VLBI positions are found to increase with the CARMS values. The majority of the sources with statistically significant position differences are associated with the sources having extended structure. Radio source structure is the one of the major factors of these position differences, and it can be the dominate factor for a number of sources. The vectors of the $Gaia$ and VLBI position differences are parallel to the radio-jet directions, which is confirmed with stronger evidence., Comment: 13 pages, 7 figures, and 5 tables. Accepted to be published in Astronomy & Astrophysics

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

Relatively few X-ray sources are known that have low-mass galaxies as hosts. This is an important restriction on studies of active galactic nuclei (AGNs), hence black holes, and of X-ray binaries (XRBs) in low-mass galaxies; addressing it requires very large samples of both galaxies and X-ray sources. Here, we have matched the X-ray point sources found in the XXL-N field of the XXL survey (with an X-ray flux limit of ∼6 × 10−15 erg s−1 cm−2 in the [0.5–2] keV band) to galaxies with redshifts from the Galaxy And Mass Assembly (GAMA) G02 survey field (down to a magnitude limit r = 19.8) in order to search for AGNs and XRBs in GAMA galaxies, particularly those of low optical luminosity or stellar mass (fainter than Mr = −19 or $M_* \lesssim 10^{9.5}\, \mathrm{M}_{\odot }$). Out of a total of 1200 low-mass galaxies in the overlap region, we find a total of 28 potential X-ray source hosts, though this includes possible background contaminants. From a combination of photometry (optical and infrared colours), positional information, and optical spectra, we deduce that most of the ≃20 X-ray sources genuinely in low-mass galaxies are high-mass X-ray binaries in star-forming galaxies. None of the matched sources in a low-mass galaxy has a BPT classification as an AGN, and even ignoring this requirement, none passes both criteria of close match between the X-ray source position and optical galaxy centre (separation ≤3 arcsec) and high [O iii] line luminosity (above 1040.3 erg s−1).

Nuclear Activity in Galaxies Across Cosmic Time, Proceedings of the conference held 7-11 October 2019 in Addis Ababa, Ethiopia. Edited by Mirjana Pović et al. Proceedings of the International Astronomical Union, Volume 356, pp. 158-162, We tested how the AGN contribution (5%-75% of the total flux) may affect different morphological parameters commonly used in galaxy classification. We carried out all analysis at z ∼ 0 and at higher redshifts that correspond to the COSMOS field. Using a local training sample of >2000 visually classified galaxies, we carried out all measurements with and without the central source, and quantified how the contribution of a bright nuclear point source could affect different morphological parameters, such as: Abraham and Concelice-Bershady indices, Gini, Asymmetry, M20 moment of light, and Smoothness. We found that concentration indexes are less sensitive to both redshift and brightness in comparison to the other parameters. We also found that all parameters change significantly with AGN contribution. At z ∼ 0, up to a 10% of AGN contribution the morphological classification will not be significantly affect, but for ⩾ 25% of AGN contribution late-type spirals follow the range of parameters of elliptical galaxies and can therefore be misclassified early types., TGW acknowledges the support from Bule Hora University. TGW, MP, and ZBA acknowledge financial support from the Ethiopian Space Science and Technology Institute (ESSTI) under the Ethiopian Ministry of Innovation and Technology (MoIT). TGW and JM acknowledge support by the grant CSIC I-COOP 2017, COOPA20168. MP, JM, JP, and IMP acknowledge the support from the Spanish Ministry of Science, Innovation and Universities (MICIU) through project AYA2016-76682C3-1-P and 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)

We present an analysis of the mid-infrared properties of the youngest radio galaxies, based on low-resolution data provided by the WISE satellite. We restrict our analysis to sources with available X-ray data that constitute the earliest phase of radio galaxy evolution, that is, those classified as gigahertz peaked spectrum and/or compact symmetric objects. We compare the MIR color distribution in our sample to that in the general population of local active galactic nuclei (AGN), in the population of evolved FR II radio galaxies, and also in the population of radio galaxies with recurrent jet activity. We conclude that the triggering of radio jets in AGN does not differentiate between elliptical hosts with substantially different fractions of young stars; instead there is a relationship between the jet duty cycle and the ongoing star formation.

The classification of the optical spectra of active galactic nuclei (AGN) into different types is well founded on AGN physics, but it involves some degree of human oversight and cannot be reliably scaled to large data sets. Machine learning (ML) tackles such a classification problem in a fast and reproducible way, but is often perceived as a black box. However, ML interpretability and explainability are active research areas in computer science, increasingly providing us with tools to alleviate this issue. We applied ML interpretability tools to a classifier trained to predict AGN type from spectra, to demonstrate the use of such tools in this context. We trained a support-vector machine on 3346 high-quality, low redshift AGN spectra from SDSS DR15 with an existing reliable classification as type 1, type 2, or intermediate type. On a selection of test-set spectra, we computed the gradient of the predicted class probability and we built saliency maps. We also visualized the high-dimensional space of AGN spectra using t-distributed stochastic neighbor embedding (t-SNE), showing where the spectra for which we computed a saliency map are located. Regions that affect the predicted AGN type often coincide with physically relevant features, such as spectral lines. t-SNE visualization shows good separability of type 1 and type 2 spectra, while intermediate-type spectra either lie in-between as expected or appear mixed with type 2 spectra. Saliency maps show why a given AGN type was predicted by our classifier, resulting in a physical interpretation in terms of regions of the spectrum that affected its decision, making it no longer a black box. These regions coincide with those used by human experts such as relevant spectral lines, and are even used in a similar way, with the classifier e.g. effectively measuring the width of a line by weighing its center and its tails oppositely., Comment: Under review by A&A

Blazar variability appears to be stochastic in nature. However, a possibility of low-dimensional chaos was considered in the past, but with no unambiguous detection so far. If present, it would constrain the emission mechanism by suggesting an underlying dynamical system. We rigorously searched for signatures of chaos in Fermi-Large Area Telescope light curves of 11 blazars. The data were comprehensively investigated using the methods of nonlinear time-series analysis: phase-space reconstruction, fractal dimension, and maximal Lyapunov exponent (mLE). We tested several possible parameters affecting the outcomes, in particular the mLE, in order to verify the spuriousness of the outcomes. We found no signs of chaos in any of the analysed blazars. Blazar variability is either truly stochastic in nature or governed by high-dimensional chaos that can often resemble randomness. © 2021 The Author(s)., OO thanks the Astronomical Observatory of the Jagiellonian University for the summer internship during which this research began. MT acknowledges support by the Polish National Science Center through the OPUS grant No. 2017/25/B/ST9/01208. The work of N. Z is supported by the South African Research Chairs Initiative (grant no. 64789) of the Department of Science and Innovation and the National Research Foundation of South Africa. J.P.-G. acknowledges funding support from Spanish public funds for research under project ESP2017-87676-C5-5-R and 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 ' isica de Andaluc ' ia (SEV-2017-0709).

We present Multi-Unit Spectroscopic Explorer (MUSE) integral-field spectroscopy of ESO 253-G003, which hosts a known active galactic nucleus (AGN) and the periodic nuclear transient ASASSN-14ko, observed as part of the All-weather MUse Supernova Integral-field of Nearby Galaxies survey. The MUSE observations reveal that the inner region hosts two AGN separated by 1.4±0.1 kpc (≍1 ′′. 7). The brighter nucleus has asymmetric broad permitted emission-line profiles and is associated with the archival AGN designation. The fainter nucleus does not have a broad emission-line component but exhibits other AGN characteristics, including vFWHM≈700 km~s−1 forbidden line emission, log10([OIII]/Hβ)≈1.1 , and high-excitation potential emission lines, such as [Fe VII] λ6086 and He II λ4686. The host galaxy exhibits a disturbed morphology with large kpc-scale tidal features, potential outflows from both nuclei, and a likely superbubble. A circular relativistic disc model cannot reproduce the asymmetric broad emission-line profiles in the brighter nucleus, but two non-axisymmetric disc models provide good fits to the broad emission-line profiles: an elliptical disc model and a circular disc + spiral arm model. Implications for the periodic nuclear transient ASASSN-14ko are discussed., MAT acknowledges support from the DOE CSGF through grant no. DE-SC0019323. BJS and CSK are supported by NSF grant no. AST-1907570. BJS is also supported by NASA grant no. 80NSSC19K1717 and NSF grants AST-1920392 and AST-1911074. CSK is supported by NSF grant no. AST-181440. KAA is supported by the Danish National Research Foundation (DNRF132). Support for JLP is provided in part by FONDECYT through grant n.1191038 and by the Ministry for the Economy, Development, and Tourism’s Millennium Science Initiative through grant no. IC120009, awarded to The Millennium Institute of Astrophysics, MAS. LG acknowledges financial support from the Spanish Ministry of Science, Innovation and Universities (MICIU) under the 2019 Ramón y Cajal program RYC2019-027683 and from the Spanish MICIU project PID2020-115253GA-I00. Parts of this research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. LG was funded by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 839090, and partially supported by the Spanish grant no. PGC2018-095317-B-C21 within the European Funds for Regional Development (FEDER). Support for TW-SH was provided by NASA through the NASA Hubble Fellowship grant no. #HST-HF2-51458.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. Based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere under ESO programme 096.D-0296(A).

AGN are a key ingredient for understanding galactic evolution. AGN-driven outflows are one of the manifestations of feedback. The AO mode for MUSE at the VLT permits to study the innermost tens of parsecs of nearby AGN in the optical. We present a detailed analysis of the ionised gas in the central regions of NGC 7130, an archetypical composite Seyfert and nuclear starburst galaxy. We achieve an angular resolution of 0.17$^{\prime\prime}$ (50 pc). We performed a multi-component analysis of the main ISM lines and identified nine kinematic components, six of which correspond to the outflow. The outflow is biconic and has velocities of a few $100\,{\rm km\,s^{-1}}$ with respect to the disc. We decompose the approaching side of the outflow into a broad and a narrow component with typical velocity dispersions below and above $\sim200\,{\rm km\,s^{-1}}$, respectively. The blueshifted narrow component has substructure, in particular a collimated plume aligned with the radio jet, indicating that it may be jet-powered. The redshifted lobe is composed of two Narrow Components and a Broad Component. An additional redshifted component is seen outside the main outflow axis. Line ratio diagnostics indicate that the outflow gas in the main axis is AGN-powered whereas the off-axis component has LINER properties. The ionised gas mass outflow rate is $\dot{M}=1.2\pm0.7\,M_{\odot}\,{\rm yr^{-1}}$ and the kinetic power is $\dot{E}_{\rm kin}=(2.7\pm2.0)\times10^{41}\,{\rm erg\,s^{-1}}$, which corresponds to $F_{\rm kin}=0.12\pm0.09\%$ of the bolometric AGN power. The combination of high angular resolution integral field spectroscopy and a careful multi-component decomposition allows a uniquely detailed view of the outflow in NGC 7130, illustrating that AGN kinematics are more complex than traditionally derived from less sophisticated data and analyses. (abridged), Published in A&A. This version of the preprint has been corrected with the modifications appearing in the A&A corrigendum

The catalogue is 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/645/A74, Since its launch in 1999, the XMM-Newton mission has compiled the largest catalogue of serendipitous X-ray sources, with the 3XMM being the third version of this catalogue. This was possible thanks to the combination of a large effective area (5000 cm2 at 1 keV) and a wide field of view (30 arcmin). The 3XMM-DR6 catalogue contains about 470 000 unique X-ray sources over an area of 982 deg2. A significant fraction of these (100 178 sources) have reliable optical, near-(NIR), and mid-infrared (MIR) counterparts in the SDSS, PANSTARRS, VIDEO, UKIDSS, and WISE surveys. In a previous paper we presented photometric redshifts for these sources using the TPZ machine-learning algorithm. About a quarter of these (22 677) have adequate photon statistics, meaning that a reliable X-ray spectrum can be extracted. Owing to both the X-ray counts selection and the optical counterpart constraint, the sample above is biased towards the bright sources. Here, we present XMMFITCAT-Z: a spectral fit catalogue for these sources using the Bayesian X-ray Analysis technique. To demonstrate the potential of the present catalogue, we comment on the optical and MIR colours of the 765 X-ray absorbed sources with NH > 1022 cm−2. We show that a considerable fraction of X-ray-selected AGNs would not be classified as AGNs following the MIR W1–W2 versus W2 selection criterion. These are AGNs with lower luminosities, where the contribution of the host galaxy to the MIR emission is non-negligible. Only one-third of obscured AGNs in X-rays present red colours or r–W2 > 6. Also, it appears that the r–W2 criterion, often used in the literature for the selection of obscured AGNs, produces very different X-ray absorbed AGN samples compared to the standard X-ray selection criteria., This work is part of the Enhanced XMM-Newton Spectralfit Database project, funded by the European Space Agency (ESA) under the PRODEX program. AR acknowledges support of this work by the PROTEAS II project (MIS 5002515), which is implemented under the “Reinforcement of the Research and Innovation Infrastructure” action, funded by the “Competitiveness, Entrepreneurship and Innovation” operational programme (NSRF 2014-2020) and co-financed by Greece and the European Union (European Regional Development Fund). AC acknowledges financial support from the Spanish Ministry MCIU under project RTI2018-096686-B-C21 (MCIU/AEI/FEDER/UE), cofunded by FEDER funds and from the Agencia Estatal de Investigación, Unidad de Excelencia María de Maeztu, ref. MDM-2017-0765.

The regular monitoring of flat-spectrum radio quasars (FSRQs) in $\gamma$-rays by Fermi-LAT since past 12 years indicated six sources who exhibited extreme $\gamma$-ray outbursts crossing daily flux of $10^{-5}$ photons/cm$^{2}$/s. We obtained nearly-simultaneous multi-wavelength data of these sources in radio to $\gamma$-ray waveband from OVRO, Steward Observatory, SMARTS, Swift-UVOT, Swift-XRT, and Fermi-LAT. The time-averaged broadband Spectral Energy Distributions (SEDs) of these sources in quiescent states were studied to get an idea about the underlying baseline radiation processes. We modeled the SEDs using one-zone leptonic synchrotron and inverse-Compton emission scenario from broken power-law electron energy distribution inside a spherical plasma blob, relativistically moving down a conical jet. The model takes into account inverse-Compton scattering of externally and locally originated seed photons in the jet. The big blue bumps visible in quiescent state SEDs helped to estimate the accretion disk luminosities and central black hole masses. We found a correlation between the magnetic field inside the emission region and the ratio of emission region distance to disk luminosity, which implies that the magnetic field decreases with an increase in emission region distance and decrease in disk luminosity, suggesting a disk-jet connection. The high-energy index of the electron distribution was also found to be correlated with observed $\gamma$-ray luminosity as $\gamma$-rays are produced by high-energy particles. In most cases, kinetic power carried by electrons can account for jet radiation power as jets become radiatively inefficient during quiescent states., Comment: 13 pages, 7 figures. Accepted in MNRAS

Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $\nu_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $\gamma$-ray telescope FACT during a high $\gamma$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $\gamma$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $\gamma$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $\Gamma$ index of the intrinsic spectrum in the VHE $\gamma$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz, Comment: Accepted in MNRAS. Corresponding authors: M. Manganaro, A. Arbet Engels, D. Dorner

We use the Gemini NIFS instrument to map the H22.1218μm and Brγ flux distributions in the inner 0.04–2 kpc of a sample of 36 nearby active galaxies (0.001 ≲ z ≲ 0.056) at spatial resolutions from 4 to 250 pc. We find extended emission in 34 galaxies. In ∼55 per cent of them, the emission in both lines is most extended along the galaxy major axis, while in the other 45 per cent the extent follows a distinct orientation. The emission of H2 is less concentrated than that of Brγ, presenting a radius that contains half of the flux 60 per cent greater, on average. The H2 emission is driven by thermal processes – X-ray heating and shocks – at most locations for all galaxies, where 0.4 6 (seen in 40 per cent of the galaxies), shocks are the main H2 excitation mechanism, while in regions with H2/Brγ < 0.4 (25 per cent of the sample) the H2 emission is produced by fluorescence. The only difference we found between type 1 and type 2 active galactic nucleus (AGN) was in the nuclear emission-line equivalent widths that are smaller in type 1 than in type 2 due to a larger contribution to the continuum from the hot dusty torus in the former. The gas masses in the inner 125 pc radius are in the range 101−104 M⊙ for the hot H2 and 103−106 M⊙ for the ionized gas and would be enough to power the AGN in our sample for 105−108 yr at their current accretion rates.

We study the impact of AGN ionised outflows on star formation in high-redshift AGN hosts, by combining NIR IFS observations, mapping the H$\alpha$ emission and [OIII] outflows, with matched-resolution observations of the rest-frame FIR emission. We present high-resolution ALMA Band 7 observations of eight X-ray selected AGN at z~2 from the SUPER sample, targeting the rest-frame ~260 um continuum at ~2 kpc (0.2'') resolution. We detected 6 out of 8 targets with S/N>10 in the ALMA maps, with continuum flux densities F = 0.27-2.58 mJy and FIR half-light radii Re = 0.8-2.1 kpc. The FIR Re of our sample are comparable to other AGN and star-forming galaxies at a similar redshift from the literature. However, we find that the mean FIR size in X-ray AGN (Re = 1.16+/- 0.11 kpc) is slightly smaller than in non-AGN (Re = 1.69+/-0.13 kpc). From SED fitting, we find that the main contribution to the 260 um flux density is dust heated by star formation, with < 4% contribution from AGN-heated dust and < 1% from synchrotron emission. The majority of our sample show different morphologies for the FIR (mostly due to reprocessed stellar emission) and the ionised gas emission (H$\alpha$ and [OIII], mostly due to AGN emission). This could be due to the different locations of dust and ionised gas, the different sources of the emission (stars and AGN), or the effect of dust obscuration. We are unable to identify any residual H$\alpha$ emission, above that dominated by AGN, that could be attributed to star formation. Under the assumption that the FIR emission is a reliable tracer of obscured star formation, we find that the obscured star formation activity in these AGN host galaxies is not clearly affected by the ionised outflows. However, we cannot rule out that star formation suppression is happening on smaller spatial scales than the ones we probe with our observations (< 2 kpc) or on different timescales., Comment: 32 pages, 13 figures. Accepted for publication in A&A

We performed a detailed spectral and timing analysis of a Seyfert 1 galaxy Mrk 509 using data from the Neil Gehrels Swift observatory that spanned over ~13 years between 2006 and 2019. To study the variability properties from the optical/UV to X-ray emission, we used a total of 275 pointed observations in this work. The average spectrum over the entire duration exhibits a strong soft X-ray excess above the power-law continuum. The soft X-ray excess is well described by two thermal components with temperatures of kT_BB1 ~120 eV and kT_BB2 ~460 eV. The warm thermal component is likely due to the presence of an optically thick and warm Comptonizing plasma in the inner accretion disk. The fractional variability amplitude is found to be decreasing with increasing wavelength, i.e. from the soft X-ray to UV/optical emission. However, the hard X-ray (2-8 keV) emission shows very low variability. The strength of the correlation within the UV and the optical bands (0.95-0.99) is found to be stronger than the correlation between the UV/Optical and X-ray bands (0.40-0.53). These results clearly suggest that the emitting regions of the X-ray and UV/optical emission are likely distinct or partly interacting. Having removed the slow variations in the light curves, we find that the lag spectrum is well described by the 4/3 rule for the standard Shakura-Sunyaev accretion disk when we omit X-ray lags. All these results suggest that the real disk is complex, and the UV emission is likely reprocessed in the accretion disk to give X-ray and optical emission., This article has been accepted to the Publications of the Astronomical Society of Australia (PASA). This article contains 17 pages and 13 figures

Nuclear Activity in Galaxies Across Cosmic Time, Proceedings of the conference held 7-11 October 2019 in Addis Ababa, Ethiopia. Edited by Mirjana Pović et al. Proceedings of the International Astronomical Union, Volume 356, pp. 3-8, In this paper we summarise the research that is currently going on in Ethiopia and East-Africa in extragalactic astronomy and physics of active galaxies and active galactic nuclei (AGN). The study is focused on some of the still open questions such as: what are the stellar ages and populations of ultra hard X-ray detected AGN and connection between AGN and their host galaxies?, what are the properties of AGN in galaxy clusters and the role that environment has in triggering nuclear activity?, what are the morphological properties of AGN and how precisely we can deal with morphological classification of active galaxies?, what are the properties of galaxies in the green valley and the role of AGN in galaxy evolution?, and what are the properties of radio-loud and radio-quiet quasars (QSO) and dichotomy between the two?. Each of these questions has been developed under one specific project that will be briefly introduced. These projects involve 6 PhD and 3 MSc students and collaborations between Ethiopia, Rwanda, South Africa, Uganda, Tanzania, Spain, Italy, and Chile. With all projects we aim: first, to contribute to our general knowledge about AGN, and second, to contribute to the development in astronomy and science in Ethiopia and East-Africa., Support of the Ethiopian Space Science and Technology Institute (ESSTI) under the Ethiopian Ministry of Innovation and Technology (MInT) and Spanish MEC under AYA2016-76682-C3-1- P are gratefully acknowledged.

We present a dynamical study of the narrow-line regions in two nearby type 2 quasars (QSO2s). We construct dynamical models based on detailed photoionization models of the emission-line gas, including the effects of internal dust, to apply to observations of large-scale outflows from these active galactic nuclei (AGNs). We use Mrk 477 and Mrk 34 in order to test our models against recent Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) observations of [O iii] emission-line kinematics, since these AGNs possess more energetic outflows than found in Seyfert galaxies. We find that the outflows within 500 pc are consistent with radiative acceleration of dusty gas, however the outflows in Mrk 34 are significantly more extended and may not be directly accelerated by radiation. We characterize the properties of X-ray winds found from the expansion of [O iii]-emitting gas close to the black hole. We show that such winds possess the kinetic energy density to disturb [O iii] gas at ∼1.8 kpc, and have sufficient energy to entrain the [O iii] clouds at ∼1.2 kpc. Assuming that the X-ray wind possesses the same radial mass distribution as the [O iii] gas, we find that the peak kinetic luminosity for this wind is 2 per cent of Mrk 34’s bolometric luminosity, which is in the 0.5–5 per cent range required by some models for efficient feedback. Our work shows that, although the kinetic luminosity as measured from [O iii]-emitting gas is frequently low, X-ray winds may provide more than one order of magnitude higher kinetic power.

Nuclear Activity in Galaxies Across Cosmic Time, Proceedings of the conference held 7-11 October 2019 in Addis Ababa, Ethiopia. Edited by Mirjana Pović et al. Proceedings of the International Astronomical Union, Volume 356, pp. 61-65, We present a new tool for the analysis of the optical emission lines of the gas in the Narrow Line Region (NLR) around Active Galactic Nuclei (AGNs). This new tool can be used in large samples of objects in a consistent way using different sets of optical emission-lines taking into the account possible variations from the O/H - N/O relation. The code compares certain observed emission-line ratios with the predictions from a large grid of photoionization models calculated under the most usual conditions in the NLR of AGNs to calculate the total oxygen abundance, nitrogen-to-oxygen ratio and ionization parameter. We applied our method to a sample of Seyfert 2 galaxies with optical emission-line fluxes from the literature. Our results confirm the high metallicity of the objects of the sample and provide consistent values with the direct method. The usage of models to calculate precise ICFs is mandatory when only optical emission lines are available to derive chemical abundances using the direct method in NLRs of AGN., 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)

Full list of authors: Raiteri, C. M.; Villata, M.; Larionov, V. M.; Jorstad, S. G.; Marscher, A. P.; Weaver, Z. R.; Acosta-Pulido, J. A.; Agudo, I.; Andreeva, T.; Arkharov, A.; Bachev, R.; Benítez, E.; Berton, M.; Björklund, I.; Borman, G. A.; Bozhilov, V.; Carnerero, M. I.; Carosati, D.; Casadio, C.; Chen, W. P. Damljanovic, G.; D'Ammando, F.; Escudero, J.; Fuentes, A.; Giroletti, M.; Grishina, T. S.; Gupta, A. C.; Hagen-Thorn, V. A.; Hart, M.; Hiriart, D.; Hou, W. -J.; Ivanov, D.; Kim, J. -Y.; Kimeridze, G. N.; Konstantopoulou, C.; Kopatskaya, E. N.; Kurtanidze, O. M.; Kurtanidze, S. O.; Lähteenmäki, A.; Larionova, E. G.; Larionova, L. V.; Marchili, N.; Markovic, G.; Minev, M.; Morozova, D. A.; Myserlis, I.; Nakamura, M.; Nikiforova, A. A.; Nikolashvili, M. G.; Otero-Santos, J.; Ovcharov, E.; Pursimo, T.; Rahimov, I.; Righini, S.; Sakamoto, T.; Savchenko, S. S.; Semkov, E. H.; Shakhovskoy, D.; Sigua, L. A.; Stojanovic, M.; Strigachev, A.; Thum, C.; Tornikoski, M.; Traianou, E.; Troitskaya, Y. V.; Troitskiy, I. S.; Tsai, A.; Valcheva, A.; Vasilyev, A. A.; Vince, O.; Zaharieva, E., Among active galactic nuclei, blazars show extreme variability properties. We here investigate the case of the BL Lac object S4 0954+65 with data acquired in 2019-2020 by the Transiting Exoplanet Survey Satellite (TESS) and by the Whole Earth Blazar Telescope (WEBT) Collaboration. The 2-min cadence optical light curves provided by TESS during three observing sectors of nearly 1 month each allow us to study the fast variability in great detail. We identify several characteristic short-term time-scales, ranging from a few hours to a few days. However, these are not persistent, as they differ in the various TESS sectors. The long-term photometric and polarimetric optical and radio monitoring undertaken by the WEBT brings significant additional information, revealing that (i) in the optical, long-term flux changes are almost achromatic, while the short-term ones are strongly chromatic; (ii) the radio flux variations at 37 GHz follow those in the optical with a delay of about 3 weeks; (iii) the range of variation of the polarization degree and angle is much larger in the optical than in the radio band, but the mean polarization angles are similar; (iv) the optical long-term variability is characterized by a quasi-periodicity of about 1 month. We explain the source behaviour in terms of a rotating inhomogeneous helical jet, whose pitch angle can change in time. © 2021 The Author(s)., This paper includes data collected by the TESS mission. Funding for the TESS mission is provided by the NASA Explorer Program. Partly based on observations made with the Nordic Optical Telescope, owned in collaboration by the University of Turku and Aarhus University, and operated jointly by Aarhus University, the University of Turku and the University of Oslo, representing Denmark, Finland and Norway, the University of Iceland and Stockholm University at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias. This article is partly based on observations made in the Observatorios de Canarias del IAC with the Liverpool telescope operated on the island of La Palma by the Liverpool John Moores University in the Observatorio del Roque de Los Muchachos. This article is partly based on observations made with the LCOGT Telescopes, one of whose nodes is located at the Observatorios de Canarias del IAC on the island of Tenerife in the Observatorio del Teide. This article is partly based on observations made with the IAC-80 operated on the island of Tenerife by the Instituto de Astrofisica de Canarias in the Spanish Observatorio del Teide. Many thanks are due to the IAC support astronomers and telescope operators for supporting the observations at the IAC-80 telescope. This publication makes use of data obtained at Metsähovi Radio Observatory, operated by Aalto University in Finland. This research has made use of NASA’s Astrophysics Data System and of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. The research at Boston University was supported by NASA grants 80NSSC20K1566 (Fermi Guest Investigator Program) and 80NSSC21K0243 (TESS Guest Investigator Program). This study was based (in part) on observations conducted using the 1.8 m Perkins Telescope Observatory (PTO) in Arizona (USA), which is owned and operated by Boston University. GD, MS, GM, and OV acknowledge the observing grant support from the Institute of Astronomy and Rozhen NAO BAS through the bilateral joint research project ‘Gaia Celestial Reference Frame (CRF) and fast variable astronomical objects’ (2020–2022, leader is G.Damljanovic), and support by the Ministry of Education, Science and Technological Development of the Republic of Serbia (contract No 451-03-68/2020-14/200002) This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grants DN 18-13/2017, KP-06-H28/3 (2018), KP-06-H38/4 (2019), and KP-06-KITAJ/2 (2020). SOK acknowledges financial support by Shota Rustaveli National Science Foundation of Georgia under contract PHDF-18-354 EB acknowledges support from DGAPA-PAPIIT GRANT IN113320. This work is partly based upon observations carried out at the Observatorio Astronómico Nacional on the Sierra San Pedro Mártir (OAN- SPM), Baja California, Mexico. We acknowledge support by Bulgarian National Science Fund under grant DN18-10/2017 and National RI Roadmap Project DO1-383/18.12.2020 of the Ministry of Education and Science of the Republic of Bulgaria. IA acknowledges financial support from the Spanish ‘Ministerio de Ciencia e Innovación’ (MCINN) through the ‘Center of Excellence Severo Ochoa’ award for the Instituto de Astrofísica de Andalucía-CSIC (SEV-2017-0709). Acquisition and reduction of the POLAMI data was supported in part by MICINN through grants AYA2016-80889-P and PID2019-107847RB-C44. The POLAMI observations were carried out at the IRAM 30m Telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain).