Search

The question of whether and how the properties of radio galaxies (RGs) are connected with the large-scale environment is still an open issue. For this work we measured the large-scale galaxies' density around RGs present in the revised Third Cambridge Catalog of radio sources (3CR) with 0.02 < z < 0.3. The goal is to determine whether the accretion mode and morphology of RGs are related to the richness of the environment. We considered RGs at 0.05 < z < 0.3 for a comparison between optical spectroscopic classes, and those within 0.02 < z < 0.1 to study the differences between the radio morphological types. Photometric data from the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) survey were used to search for "red sequences" within an area of 500 kpc of radius around each RG. We find that 1) RGs span over a large range of local galaxies' density, from isolated sources to those in rich environments, 2) the richness distributions of the various classes are not statistically different, and 3) the radio luminosity is not connected with the source environment. Our results suggest that the RG properties are independent of the local galaxies density, which is in agreement with some previous analyses, but contrasting with other studies. We discuss the possible origin of this discrepancy. An analysis of a larger sample is needed to put out results on a stronger statistical basis., Comment: Pre-proofs version - Accepted for publication in A&A

We present a long-term and intraday variability study on optical multiwaveband ($U\!BV\!RI$) data from the blazar AO 0235+164 collected by various telescopes for $\sim$44 years (1975--2019). The blazar was found to be significantly variable over the years in all wavebands with a variation of about six magnitudes between its low and active states. The variations in the different wavebands are highly correlated without any time lag. We did not observe any significant trend in color variation with time, but we observed a bluer-when-brighter trend between the $B-I$ color index and the $R$-magnitude. Optical $BV\!R$-band spectral energy distributions always show a convex shape. Significant intraday variability was frequently seen in the quasi-simultaneous observations of AO\,0235+164 made on 22 nights in $R$ and $V$-bands by the CASLEO and CAHA telescopes during 1999--2019. We also estimated the central supermassive black-hole mass of $7.9\times10^7 M_{\odot}$ by analyzing the broad Mg II emission line in AO 0235+164's spectrum. We briefly explore the probable physical scenarios responsible for the observed variability., Comment: 27 pages, 10 figures. Accepted for publication in The Astrophysical Journal Supplement Series (ApJS)

Gaia DR3 contains 1.8 billion sources with G-band photometry, 1.5 billion of which with BP and RP photometry, complemented by positions on the sky, parallax, and proper motion. The median number of field-of-view transits in the three photometric bands is between 40 and 44 measurements per source and covers 34 months of data collection. We pursue a classification of Galactic and extra-galactic objects that are detected as variable by Gaia across the whole sky. Supervised machine learning (eXtreme Gradient Boosting and Random Forest) was employed to generate multi-class, binary, and meta-classifiers that classified variable objects with photometric time series in the G, BP, and RP bands. Classification results comprise 12.4 million sources (selected from a much larger set of potential variable objects) and include about 9 million variable stars classified into 22 variability types in the Milky Way and nearby galaxies such as the Magellanic Clouds and Andromeda, plus thousands of supernova explosions in distant galaxies, 1 million active galactic nuclei, and almost 2.5 million galaxies. The identification of galaxies was made possible by the artificial variability of extended objects as detected by Gaia, so they were published in the galaxy_candidates table of the Gaia DR3 archive, separate from the classifications of genuine variability (in the vari_classifier_result table). The latter contains 24 variability classes or class groups of periodic and non-periodic variables (pulsating, eclipsing, rotating, eruptive, cataclysmic, stochastic, and microlensing), with amplitudes from a few milli-magnitudes to several magnitudes., Comment: 106 pages, 513 figure panels, 4 tables. Accepted by Astronomy & Astrophysics (in press)

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will detect an unprecedentedly large sample of actively accreting supermassive black holes with typical accretion disk (AD) sizes of a few light days. This brings us to face challenges in the reverberation mapping (RM) measurement of AD sizes in active galactic nuclei (AGNs) using interband continuum delays. We examine the effect of LSST cadence strategies on AD RM using our metric AGNTimeLagMetric. It accounts for redshift, cadence, the magnitude limit, and magnitude corrections for dust extinction. Running our metric on different LSST cadence strategies, we produce an atlas of the performance estimations for LSST photometric RM measurements. We provide an upper limit on the estimated number of quasars for which the AD time lag can be computed within 01000 sources in each Deep Drilling field (DDF, 10 sq. deg) in any filter, with the redshift distribution of these sources peaking at z~1. We find the LSST observation strategies with a good cadence (~ 5 days) and a long cumulative season (~9 yr), as proposed for LSST DDF, are favored for the AD size measurement. We create synthetic LSST light curves for the most suitable DDF cadences and determine RM time lags to demonstrate the impact of the best cadences based on the proposed metric., Comment: accepted for publication in the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey's first light., Comment: 202 pages (in book format) 34 figures plus chapter heading figures (13)

One of the novelties of the Gaia-DR3 with respect to the previous data releases is the publication of the multiband light curves of about 1 million AGN. The goal of this work was the creation of a catalogue of variable AGN, whose selection was based on Gaia data only. We first present the implementation of the methods to estimate the variability parameters into a specific object study module for AGN. Then we describe the selection procedure that led to the definition of the high-purity variable AGN sample and analyse the properties of the selected sources. We started from a sample of millions of sources, which were identified as AGN candidates by 11 different classifiers based on variability processing. Because the focus was on the variability properties, we first defined some pre-requisites in terms of number of data points and mandatory variability parameters. Then a series of filters was applied using only Gaia data and the Gaia Celestial Reference Frame 3 (Gaia-CRF3) sample as a reference.The resulting Gaia AGN variable sample, named GLEAN, contains about 872000 objects, more than 21000 of which are new identifications. We checked the presence of contaminants by cross-matching the selected sources with a variety of galaxies and stellar catalogues. The completeness of GLEAN with respect to the variable AGN in the last Sloan Digital Sky Survey quasar catalogue is about 47%, while that based on the variable AGN of the Gaia-CRF3 sample is around 51%. From both a comparison with other AGN catalogues and an investigation of possible contaminants, we conclude that purity can be expected to be above 95%. Multiwavelength properties of these sources are investigated. In particular, we estimate that about 4% of them are radio-loud. We finally explore the possibility to evaluate the time lags between the flux variations of the multiple images of strongly lensed quasars, and show one case., Comment: 19 pages, 31 figures, 2 table. This paper is part of Gaia Data Release 3 (DR3). In press for A&A

The study of short-term variability properties in AGN jets has the potential to shed light on their particle acceleration and emission mechanisms. We report results from a four-day coordinated multi-wavelength campaign on the highly-peaked blazar (HBL) Mkn 421 in 2019 January. We obtained X-ray data from AstroSAT, BVRI photometry with the Whole Earth Blazar Telescope (WEBT), and TeV data from FACT to explore short-term multi-wavelength variability in this HBL. The X-ray continuum is rapidly variable on time-scales of tens of ks. Fractional variability amplitude increases with energy across the synchrotron hump, consistent with previous studies; we interpret this observation in the context of a model with multiple cells whose emission spectra contain cutoffs that follow a power-law distribution. We also performed time-averaged and time-resolved (time-scales of 6 ks) spectral fits; a broken power-law model fits all spectra well; time-resolved spectral fitting reveals the usual hardening when brightening behaviour. Intra-X-ray cross correlations yield evidence for the 0.6-0.8 keV band to likely lead the other bands by an average of 4.6 +- 2.6 ks, but only during the first half of the observation. The source displayed minimal night-to-night variability at all wavebands thus precluding significant interband correlations during our campaign. The broadband SED is modeled well with a standard one-zone leptonic model, yielding jet parameters consistent with those obtained from previous SEDs of this source., Comment: This article has been accepted for publication in The Monthly Notices of the Royal Astronomical Society (2022), and is published in Volume 513, Issue 2, pp.1662-1679. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved. 19 pages; 12 figures

We present a long term optical $R$ band light curve analysis of the gravitationally lensed blazar AO 0235+164 in the time span 1982 - 2019. Several methods of analysis lead to the result that there is a periodicity of ~8.13 years present in these data. In addition, each of these five major flares are apparently double-peaked, with the secondary peak following the primary one by ~2 years. Along with the well known system, OJ 287, our finding constitutes one of the most secure cases of long term quasi-periodic optical behaviour in a blazar ever found. A binary supermassive black hole system appears to provide a good explanation for these results., Comment: 7 pages, 3 figures, 1 table, Accepted for publication in MNRAS

With their emission mainly coming from a relativistic jet pointing towards us, blazars are fundamental sources to study extragalactic jets and their central engines, consisting of supermassive black holes (SMBHs) fed by accretion discs. They are also candidate sources of high-energy neutrinos and cosmic rays. Because of the jet orientation, the non-thermal blazar emission is Doppler beamed; its variability is unpredictable and occurs on time-scales from less than one hour to years. The comprehension of the diverse mechanisms producing the flux and spectral changes requires well-sampled multiband light curves on long time periods. In particular, outbursts are the best test bench to shed light on the underlying physics, especially when studied in a multiwavelength context. The Vera C. Rubin Legacy Survey of Space and Time (Rubin-LSST) will monitor the southern sky for ten years in six photometric bands, offering a formidable tool to study blazar variability features in a statistical way. The alert system will allow us to trigger follow-up observations of outstanding events, especially at high (keV-to-GeV) and very high (TeV) energies. We here examine the simulated Rubin-LSST survey strategies with the aim of understanding which cadences are more suitable for the blazar variability science. Our metrics include light curve and colour sampling. We also investigate the problem of saturation, which will affect the brightest and many flaring sources, and will have a detrimental impact on follow-up observations., Comment: 18 pages, 12 figures, accepted for publication in the ApJS Rubin Cadence Focus Issue

The limiting temporal resolution of a time-domain survey in detecting transient behavior is set by the time between observations of the same sky area. We analyze the distribution of visit separations for a range of Vera C. Rubin Observatory cadence simulations. Current simulations are strongly peaked at the 22 minute visit pair separation and provide effectively no constraint on temporal evolution within the night. This choice will necessarily prevent Rubin from discovering a wide range of astrophysical phenomena in time to trigger rapid followup. We present a science-agnostic metric to supplement detailed simulations of fast-evolving transients and variables and suggest potential approaches for improving the range of timescales explored., Comment: 7 pages, 2 figures. Submitted to the ApJS Rubin Cadence Focus Issue

Active galactic nuclei come in many varieties. A minority of them are radio-loud, and exhibit two opposite prominent plasma jets extending from the proximity of the supermassive black hole up to megaparsec distances. When one of the relativistic jets is oriented closely to the line of sight, its emission is Doppler beamed and these objects show extreme variability properties at all wavelengths. These are called "blazars". The unpredictable blazar variability, occurring on a continuous range of time-scales, from minutes to years, is most effectively investigated in a multi-wavelength context. Ground-based and space observations together contribute to give us a comprehensive picture of the blazar emission properties from the radio to the gamma-ray band. Moreover, in recent years, a lot of effort has been devoted to the observation and analysis of the blazar polarimetric radio and optical behaviour, showing strong variability of both the polarisation degree and angle. The Whole Earth Blazar Telescope (WEBT) Collaboration, involving many tens of astronomers all around the globe, has been monitoring several blazars since 1997. The results of the corresponding data analysis have contributed to the understanding of the blazar phenomenon, particularly stressing the viability of a geometrical interpretation of the blazar variability. We review here the most significant polarimetric results achieved in the WEBT studies., Comment: Review published in "Galaxies" as part of the Special Issue "Polarimetry as a Probe of Magnetic Fields in AGN Jets", Academic Editors: Margo Aller, Jose L. G\'omez and Eric Perlman

We present a unique method that allows the LSST to scan the sky for stellar variability on short timescales. The operational component of the strategy requires LSST to take star trail images. The image processing component uses deep learning to sift for transient events on timescales down to 10 ms. We advocate for enabling this observing mode with LSST, as coupling this capability with the LSST's tremendous 319.5 m$^2$deg$^2$ etendue will produce the first wide area optical survey of the universe on these timescales. We explain how these data will advance both planned lines of investigation and enable new research in the areas of stellar flares, cataclysmic variables, active galactic nuclei, Kuiper Belt objects, gamma-ray bursts, and fast radio bursts.

The detection of gamma-ray emission from narrow-line Seyfert 1 galaxies (NLSy1) has challenged the idea that large black hole (BH) masses ($\ge$10$^8$ M$_{\odot}$) are needed to launch relativistic jets. We present near-infrared imaging data of the gamma-ray-emitting NLSy1 PKS 1502+036 obtained with the Very Large Telescope. Its surface brightness profile, extending to $\sim$ 20 kpc, is well described by the combination of a nuclear component and a bulge with a Sersic index $n$ = 3.5, which is indicative of an elliptical galaxy. A circumnuclear structure observed near PKS 1502+036 may be the result of galaxy interactions. A BH mass of about $\sim 7 \times 10^{8}$ M$_{\odot}$ has been estimated by the bulge luminosity. The presence of an additional faint disc component cannot be ruled out with the present data, but this would reduce the BH mass estimate by only $\sim$ 30%. These results, together with analogous findings obtained for FBQS J1644+2619, indicate that the relativistic jets in gamma-ray-emitting NLSy1 are likely produced by massive black holes at the center of elliptical galaxies., Comment: 5 pages, 2 figures, 1 table, and 2 supplementary figures. Monthly Notices of the Royal Astronomical Society Letters, Volume 478, Issue 1, L66-L71

The discovery of $\gamma$-ray emission from radio-loud narrow-line Seyfert 1 (NLSy1) galaxies has questioned the need for large black hole masses (> 10$^8$ M$_{\odot}$) to launch relativistic jets. We present near-infrared data of the $\gamma$-ray-emitting NLSy1 FBQS J1644+2619 that were collected using the camera CIRCE (Canarias InfraRed Camera Experiment) at the 10.4-m Gran Telescopio Canarias to investigate the structural properties of its host galaxy and to infer the black hole mass. The 2D surface brightness profile is modelled by the combination of a nuclear and a bulge component with a S\'ersic profile with index $n$ = 3.7, indicative of an elliptical galaxy. The structural parameters of the host are consistent with the correlations of effective radius and surface brightness against absolute magnitude measured for elliptical galaxies. From the bulge luminosity, we estimated a black hole mass of (2.1$\pm$0.2) $\times$10$^8$ M$_{\odot}$, consistent with the values characterizing radio-loud active galactic nuclei., Comment: 5 pages, 3 figures, 1 table. Monthly Notices of the Royal Astronomical Society Letter, Vol. 469, L11-L15

In a previous paper, we proposed a new method to select low-power BL Lacs (LPBLs) based on mid-infrared emission and flux contrast through the Ca II spectral break; that study led to the selection of a complete sample formed by 34 LPBLs with 0.05

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings., Comment: 9 pages, 4 Figures, ApJ Letter accepted

Among active galactic nuclei, BL Lac objects show extreme properties that have been interpreted as the effect of relativistic beaming on the emission from a plasma jet oriented close to the line of sight. The Doppler amplification of the jet emission makes them ideal targets for studying jet physics. In particular, low-power BL Lacs (LPBL) are very interesting because they probe the jet formation and emission processes at the lowest levels of accretion. However, they are difficult to identify since their emission is swamped by the radiation from the host galaxy in most observing bands. In this paper we propose a new LPBL selection method based on the mid-infrared emission, in addition to the traditional optical indices. We considered the radio-selected sample of Best & Heckman (2012, MNRAS, 421, 1569) and cross-matched it with the WISE all-sky survey. In a new diagnostic plane including the W2-W3 color and the Dn(4000) index, LPBL are located in a region scarcely populated by other sources. By filtering objects with small emission line equivalent width, we isolated 36 LPBL candidates up to redshift 0.15. Their radio luminosity at 1.4 GHz spans the range log L_r = 39.2-41.5 [erg/s]. Considering the completeness of our sample, we analyzed the BL Lac luminosity function (RLF), finding a dramatic paucity of LPBL with respect to the extrapolation of the RLF toward low power. This requires a break in the RLF located at log L_r~40.6 [erg/s]. The consequent peak in the BL Lacs number density is possibly the manifestation of a minimum power required to launch a relativistic jet., Comment: Accepted for publication in A&A

We constrain the behavior of the radio luminosity function (RLF) of two classes of active galactic nuclei (AGN) namely AGN of low radio power (LRP) and BL Lac objects. The extrapolation of the observed steep RLFs to low power predicts a space density of such objects that exceeds that of the sources that can harbor them and this requires a break to a shallower slope. For LRP AGN we obtain P_br,LRP > 10^20.5 W/Hz at 1.4 GHz to limit their density to be smaller than that of elliptical galaxies with black hole masses M_BH > 10^7.5 solar masses. By combining this value with the limit derived by the observations the break must occur at P_br,LRP~10^20.5-10^21.5 W/Hz. For BL Lacs we find P_br,BLLAC > 10^23.3 W/Hz otherwise they would outnumber the density of weak-lined and compact radio sources, while the observations indicate P_br,BLLAC < 10^24.5 W/Hz. In the framework of the AGN unified model a low luminosity break in the RLF of LRP AGN must correspond to a break in the RLF of BL Lacs. The ratio between P_br,LRP and P_br,BLLAC is ~10^3, as expected for a jet Doppler factor of ~10., Comment: Accepted for publication in MNRAS letters

We present UBVRI light curves of BL Lacertae from May 2000 to January 2001, obtained by 24 telescopes in 11 countries. More than 15000 observations were performed in that period, which was the extension of the Whole Earth Blazar Telescope (WEBT) campaign originally planned for July-August 2000. Rapid flux oscillations are present all the time, involving variations up to a few tenths of mag on hour time scales, and witnessing an intense intraday activity of this source. Colour indexes have been derived by coupling the highest precision B and R data taken by the same instrument within 20 min and after subtracting the host galaxy contribution from the fluxes. The 620 indexes obtained show that the optical spectrum is weakly sensitive to the long-term trend, while it strictly follows the short-term flux behaviour, becoming bluer when the brightness increases. Thus, spectral changes are not related to the host galaxy contribution, but they are an intrinsic feature of fast flares. We suggest that the achromatic mechanism causing the long-term flux base-level modulation can be envisaged in a variation of the relativistic Doppler beaming factor, and that this variation is likely due to a change of the viewing angle. Discrete correlation function (DCF) analysis reveals the existence of a characteristic time scale of variability of about 7 h in the light curve of the core WEBT campaign, while no measurable time delay between variations in the B and R bands is found., Comment: 14 pages, 8 PostScript figures, 5 JPEG figures, in press for A&A