The properties of the AGN torus as revealed from a set of unbiased NuSTAR observations

The obscuration observed in active galactic nuclei (AGN) is mainly caused by dust and gas distributed in a torus-like structure surrounding the supermassive black hole (SMBH). However, properties of the obscuring torus of the AGN in X-ray have not been fully investigated yet due to the lack of high-quality data and proper models. In this work, we perform a broadband X-ray spectral analysis of a large, unbiased sample of obscured AGN (with line-of-sight column density 23$\le$log(NH)$\le$24) in the nearby universe which has high-quality archival NuSTAR data. The source spectra are analyzed using the recently developed borus02 model, which enables us to accurately characterize the physical and geometrical properties of AGN obscuring tori. We also compare our results obtained from the unbiased Compton thin AGN with those of Compton-thick AGN. We find that Compton thin and Compton-thick AGN may possess similar tori, whose average column density is Compton thick (N$\rm _{H,tor,ave}$ $\sim$1.4$\times$10$^{24}$ cm$^{-2}$), but they are observed through different (under-dense or over-dense) regions of the tori. We also find that the obscuring torus medium is significantly inhomogeneous, with the torus average column densities significantly different from their line-of-sight column densities (for most of the sources in the sample). The average torus covering factor of sources in our unbiased sample is c$_f$=0.67, suggesting that the fraction of unobscured AGN is $\sim$33%. We develop a new method to measure the intrinsic line-of-sight column density distribution of AGN in the nearby universe, which we find the result is in good agreement with the constraints from recent population synthesis models.Comment: 16 pages, 14 figures, 7 tables; accepted by A&

Multiwavelength Analysis of Fermi-LAT Blazars with High-Significance Periodicity: Detection of a Long-Term Rising Emission in PG 1553+113

Blazars display variable emission across the entire electromagnetic spectrum, with timescales that can range from a few minutes to several years. Our recent work has shown that a sample of five blazars exhibit hints of periodicity with a global significance $\gtrsim2\,\sigma$ at $\gamma$-ray energies, in the range of 0.1~GeV$<$E$<$800~GeV. In this work, we study their multiwavelength (MWL) emission, covering the X-ray, ultraviolet, optical, and radio bands. We show that three of these blazars present similar periodic patterns in the optical and radio bands. Additionally, fluxes in the different bands of the five blazars are correlated, suggesting a co-spatial origin. Moreover, we detect a long-term ($\approx$10 year) rising trend in the light curves of PG~1553+113, and we use it to infer possible constraints on the binary black hole hypothesis.Comment: 20 pages, 10 figures, 7 table

A Compton-thick nucleus in the dual AGN of Mrk 266

We present results of our analysis of NuSTAR data of the luminous infrared galaxy Mrk 266, which contains two nuclei, SW and NE, resolved in previous Chandra imaging. Combining with the Chandra data, we intepret the hard X-ray spectrum obtained from a NuSTAR observation as resulting from steeply rising flux from a Compton-thick AGN in the SW nucleus which is very faint in the Chandra band, confirming the previous claim of Mazzarella et al. (2012). This hard X-ray component is dominated by reflection, and its intrinsic 2-10 keV luminosity is likely to be ~1e43 erg/s. Although it is bright in soft X-ray, only moderately absorbed NE nucleus has a 2-10 keV luminosity of 4e41 erg/s, placing it in the low-luminosity AGN class. These results have implications for understanding the detectability and duty cycles of emission from dual AGN in heavily obscured mergers.Comment: Accepted for publication in Astronomy and Astrophysics. 7 pages, 5 figures. Language corrections and a minor correction in Fig. 5 include

A Compton-thick nucleus in the dual active galactic nuclei of Mrk 266

We present the results from our analysis of NuSTAR data of the luminous infrared galaxy Mrk 266, which contains two nuclei, south-western (SW) and north-eastern (NE), which were resolved in previous Chandra imaging. Combining this with the Chandra data, we intepret the hard X-ray spectrum obtained from a NuSTAR observation to result from a steeply rising flux from a Compton-thick active galactic nuclei (AGN) in the SW nucleus which is very faint in the Chandra band, confirming the previous claim. This hard X-ray component is dominated by reflection, and its intrinsic 2–10 keV luminosity is likely to be ∼1 × 10⁴³ erg s⁻¹. Although it is bright in soft X-ray, only a moderately absorbed NE nucleus has a 2–10 keV luminosity of 4 × 10⁴¹ erg s⁻¹, placing it in the low-luminosity AGN class. These results have implications for understanding the detectability and duty cycles of emission from dual AGN in heavily obscured mergers

C-GOALS II. Chandra Observations of the Lower Luminosity Sample of Nearby Luminous Infrared Galaxies in GOALS

We analyze Chandra X-ray observatory data for a sample of 63 luminous infrared galaxies (LIRGs), sampling the lower-infrared luminosity range of the Great Observatories All-Sky LIRG survey (GOALS), which includes the most luminous infrared selected galaxies in the local universe. X-rays are detected for 84 individual galaxies within the 63 systems, for which arcsecond resolution X-ray images, fluxes, infrared and X-ray luminosities, spectra and radial profiles are presented. Using X-ray and MIR selection criteria, we find AGN in (31$\pm$5)% of the galaxy sample, compared to the (38$\pm$6)% previously found for GOALS galaxies with higher infrared luminosities (C-GOALS I). Using mid-infrared data, we find that (59$\pm$9)% of the X-ray selected AGN in the full C-GOALS sample do not contribute significantly to the bolometric luminosity of the host galaxy. Dual AGN are detected in two systems, implying a dual AGN fraction in systems that contain at least one AGN of (29$\pm$14)%, compared to the (11$\pm$10)% found for the C-GOALS I sample. Through analysis of radial profiles, we derive that most sources, and almost all AGN, in the sample are compact, with half of the soft X-ray emission generated within the inner $\sim 1$ kpc. For most galaxies, the soft X-ray sizes of the sources are comparable to those of the MIR emission. We also find that the hard X-ray faintness previously reported for the bright C-GOALS I sources is also observed in the brightest LIRGs within the sample, with $L_{\rm FIR}>8\times10^{10}$ L$_{\odot}$.Comment: 24 pages, 13 figures, 11 tables, accepted for publication in A&

MAGIC Observations of the Nearby Short Gamma-Ray Burst GRB 160821B

Acciari, A. V., et al.The coincident detection of GW170817 in gravitational waves and electromagnetic radiation spanning the radio to MeV gamma-ray bands provided the first direct evidence that short gamma-ray bursts (GRBs) can originate from binary neutron star (BNS) mergers. On the other hand, the properties of short GRBs in high-energy gamma-rays are still poorly constrained, with only ∼20 events detected in the GeV band, and none in the TeV band. GRB 160821B is one of the nearest short GRBs known at z = 0.162. Recent analyses of the multiwavelength observational data of its afterglow emission revealed an optical-infrared kilonova component, characteristic of heavy-element nucleosynthesis in a BNS merger. Aiming to better clarify the nature of short GRBs, this burst was automatically followed up with the MAGIC telescopes, starting from 24 s after the burst trigger. Evidence of a gamma-ray signal is found above ∼0.5 TeV at a significance of ∼ 3σ during observations that lasted until 4 hr after the burst. Assuming that the observed excess events correspond to gamma-ray emission from GRB 160821B, in conjunction with data at other wavelengths, we investigate its origin in the framework of GRB afterglow models. The simplest interpretation with one-zone models of synchrotron-self-Compton emission from the external forward shock has difficulty accounting for the putative TeV flux. Alternative scenarios are discussed where the TeV emission can be relatively enhanced. The role of future GeV-TeV observations of short GRBs in advancing our understanding of BNS mergers and related topics is briefly addressed.We 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 and MPG; the Italian INFN and INAF; the Swiss National Fund SNF; the ERDF under the Spanish MINECO (FPA2017-87859-P, FPA2017-85668-P, FPA2017- 82729-C6-2-R, FPA2017-82729-C6-6-R, FPA2017-82729-C6-5- R, AYA2015-71042-P, AYA2016-76012-C3-1-P, ESP2017- 87055-C2-2-P, FPA2017-90566-REDC); 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 No. 320045 is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia “Severo Ochoa” SEV-2016-0588, SEV2015-0548 and SEV-2012-0234, the Unidad de Excelencia “María de Maeztu” MDM-2014-0369 and the “la Caixa” Foundation (fellowship LCF/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. K.N. is thankful for the support by Marie Skłodowska-Curie actions (H2020-MSCA-COFUND-2014, Project P-Sphere GA 665919), and JSPS KAKENHI grant No. JP20KK0067 from MEXT, Japan. L.N. acknowledges funding from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No. 664931. S.I. is supported by JSPS KAKENHI grant No. JP17K05460 from MEXT, Japan, and the RIKEN iTHEMS program

Multiwavelength variability and correlation studies of Mrk 421 during historically low X-ray and γ-ray activity in 2015-2016

Acciari, V. A., et al. (MAGIC Collaboration)We report a characterization of the multiband flux variability and correlations of the nearby (z = 0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, FACT, and other collaborations and instruments from 2014 November till 2016 June. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F>1 TeV 0.1 TeV) γ-rays, which, despite the low activity, show a significant positive correlation with no time lag. The HRkeV and HRTeV show the harder-when-brighter trend observed in many blazars, but the trend flattens at the highest fluxes, which suggests a change in the processes dominating the blazar variability. Enlarging our data set with data from years 2007 to 2014, we measured a positive correlation between the optical and the GeV emission over a range of about 60 d centred at time lag zero, and a positive correlation between the optical/GeV and the radio emission over a range of about 60 d centred at a time lag of 43+9-6 d. This observation is consistent with the radio-bright zone being located about 0.2 parsec downstream from the optical/GeV emission regions of the jet. The flux distributions are better described with a lognormal function in most of the energy bands probed, indicating that the variability in Mrk 421 is likely produced by a multiplicative process.The financial support of the German BMBF and MPG; the Italian INFN and INAF; the Swiss National Fund SNF; the ERDF under the Spanish MINECO (FPA2017-87859-P, FPA2017-85668-P, FPA2017-82729-C6-2-R, FPA2017-82729-C6-6-R, FPA2017-82729-C6-5-R, AYA2015-71042-P, AYA2016-76012-C3-1-P, ESP2017-87055-C2-2-P, FPA2017-90566-REDC); 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 and SEV-2015-0548, the Unidad de Excelencia ‘María de Maeztu’ MDM-2014-0369 and the ‘la Caixa’ Foundation (fellowship LCF/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

Study of the GeV to TeV morphology of the γ Cygni SNR (G 78.2+2.1) with MAGIC and Fermi-LAT: Evidence for cosmic ray escape

MAGIC Collaboration: et al.[Context] Diffusive shock acceleration (DSA) is the most promising mechanism that accelerates Galactic cosmic rays (CRs) in the shocks of supernova remnants (SNRs). It is based on particles scattering caused by turbulence ahead and behind the shock. The turbulence upstream is supposedly generated by the CRs, but this process is not well understood. The dominant mechanism may depend on the evolutionary state of the shock and can be studied via the CRs escaping upstream into the interstellar medium (ISM).[Aims] Previous observations of the γ Cygni SNR showed a difference in morphology between GeV and TeV energies. Since this SNR has the right age and is at the evolutionary stage for a significant fraction of CRs to escape, our aim is to understand γ-ray emission in the vicinity of the γ Cygni SNR.[Methods] We observed the region of the γ Cygni SNR with the MAGIC Imaging Atmospheric Cherenkov telescopes between 2015 May and 2017 September recording 87 h of good-quality data. Additionally, we analysed Fermi-LAT data to study the energy dependence of the morphology as well as the energy spectrum in the GeV to TeV range. The energy spectra and morphology were compared against theoretical predictions, which include a detailed derivation of the CR escape process and their γ-ray generation.[Results] The MAGIC and Fermi-LAT data allowed us to identify three emission regions that can be associated with the SNR and that dominate at different energies. Our hadronic emission model accounts well for the morphology and energy spectrum of all source components. It constrains the time-dependence of the maximum energy of the CRs at the shock, the time-dependence of the level of turbulence, and the diffusion coefficient immediately outside the SNR shock. While in agreement with the standard picture of DSA, the time-dependence of the maximum energy was found to be steeper than predicted, and the level of turbulence was found to change over the lifetime of the SNR.The financial support of the German BMBF and MPG; the Italian INFN and INAF; the Swiss National Fund SNF; the ERDF under the Spanish MINECO (FPA2017-87859-P, FPA2017-85668-P, FPA2017-82729-C6-2-R, FPA2017-82729-C6-6-R, FPA2017-82729-C6-5-R, AYA2015-71042-P, AYA2016-76012-C3-1-P, ESP2017-87055-C2-2-P, FPA2017- 90566-REDC); 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 no. 320045 is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia “Severo Ochoa” SEV-2016-0588 and SEV-2015-0548, the Unidad de Excelencia “María de Maeztu” MDM2014-0369 and the “la Caixa” F oundation (fellowship LCF/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, 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.Peer reviewe