Showing 1–50 of 57 results for author: Carnerero, M I

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Authors: Gaia Collaboration, P. Panuzzo, T. Mazeh, F. Arenou, B. Holl, E. Caffau, A. Jorissen, C. Babusiaux, P. Gavras, J. Sahlmann, U. Bastian, Ł. Wyrzykowski, L. Eyer, N. Leclerc, N. Bauchet, A. Bombrun, N. Mowlavi, G. M. Seabroke, D. Teyssier, E. Balbinot, A. Helmi, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne , et al. (390 additional authors not shown)

Abstract: Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is exp… ▽ More Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way. △ Less

Submitted 19 April, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 23 pages, accepted fro publication in A&A Letters. New version with small fixes

Optical variability of the blazar 3C 371: from minute to year timescales

Authors: J. Otero-Santos, C. M. Raiteri, J. A. Acosta-Pulido, M. I. Carnerero, M. Villata, S. S. Savchenko, D. Carosati, W. P. Chen, S. O. Kurtanidze, M. D. Joner, E. Semkov, T. Pursimo, E. Benítez, G. Damljanovic, G. Apolonio, G. A. Borman, V. Bozhilov, F. J. Galindo-Guil, T. S. Grishina, V. A. Hagen-Thorn, D. Hiriart, H. Y. Hsiao, S. Ibryamov, R. Z. Ivanidze, G. N. Kimeridze , et al. (19 additional authors not shown)

Abstract: The BL Lac object 3C 371 has been observed by the Transiting Exoplanet Survey Satellite (\textit{TESS}) for approximately a year, between July 2019 and July 2020, with an unmatched 2-minute observing cadence. In parallel, the Whole Earth Blazar Telescope (WEBT) Collaboration organized an extensive observing campaign, providing three years of continuous optical monitoring between 2018 and 2020. The… ▽ More The BL Lac object 3C 371 has been observed by the Transiting Exoplanet Survey Satellite (\textit{TESS}) for approximately a year, between July 2019 and July 2020, with an unmatched 2-minute observing cadence. In parallel, the Whole Earth Blazar Telescope (WEBT) Collaboration organized an extensive observing campaign, providing three years of continuous optical monitoring between 2018 and 2020. These datasets allow for a thorough investigation of the variability of the source. The goal of this study is evaluating the optical variability of 3C 371. Taking advantage of the remarkable cadence of \textit{TESS} data, we aim to characterize the intra-day variability (IDV) displayed by the source and identify its shortest variability timescale. With this estimate, constraints on the size of the emitting region and black hole mass can be calculated. Moreover, WEBT data is used to investigate long-term variability (LTV), including understanding spectral behaviour of the source and the polarization variability. Based on the derived characteristics, information on the origin of the variability on different timescales is extracted. We evaluated the variability applying the variability amplitude tool that quantifies how variable the emission is. Moreover, we employed common tools like ANOVA (ANalysis Of VAariance) tests, wavelet and power spectral density (PSD) analyses to characterize the shortest variability timescales present in the emission and the underlying noise affecting the data. Short- and long-term colour behaviours have been evaluated to understand the spectral behaviour of the source. The polarized emission was analyzed, studying its variability and possible rotation patterns of the electric vector position angle (EVPA). Flux distributions of IDV and LTV were also studied with the aim of linking the flux variations to turbulent and/or accretion disc related processes. △ Less

Submitted 4 April, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: Accepted in A&A, 20 pages, 19 figures

Journal ref: A&A 686, A228 (2024)

The variability patterns of the TeV blazar PG 1553+113 from a decade of MAGIC and multi-band observations

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (242 additional authors not shown)

Abstract: PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m a… ▽ More PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

Extreme photometric and polarimetric variability of blazar S4 0954+65 at its maximum optical and $γ$-ray brightness levels

Authors: C. M. Raiteri, M. Villata, M. I. Carnerero, S. S. Savchenko, S. O. Kurtanidze, V. V. Vlasyuk, A. Marchini, K. Matsumoto, C. Lorey, M. D. Joner, K. Gazeas, D. Carosati, D. O. Mirzaqulov, J. A. Acosta Pulido, I. Agudo, R. Bachev, E. Benítez, G. A. Borman, P. Calcidese, W. P. Chen, G. Damljanovic, S. A. Ehgamberdiev, D. Elsässer, M. Feige, A. Frasca , et al. (42 additional authors not shown)

Abstract: In 2022 the BL Lac object S4 0954+65 underwent a major variability phase, reaching its historical maximum brightness in the optical and $γ$-ray bands. We present optical photometric and polarimetric data acquired by the Whole Earth Blazar Telescope (WEBT) Collaboration from 2022 April 6 to July 6. Many episodes of unprecedented fast variability were detected, implying an upper limit to the size of… ▽ More In 2022 the BL Lac object S4 0954+65 underwent a major variability phase, reaching its historical maximum brightness in the optical and $γ$-ray bands. We present optical photometric and polarimetric data acquired by the Whole Earth Blazar Telescope (WEBT) Collaboration from 2022 April 6 to July 6. Many episodes of unprecedented fast variability were detected, implying an upper limit to the size of the emitting region as low as $10^{-4}$ parsec. The WEBT data show rapid variability in both the degree and angle of polarization. We analyse different models to explain the polarization behaviour in the framework of a twisting jet model, which assumes that the long-term trend of the flux is produced by variations in the emitting region viewing angle. All the models can reproduce the average trend of the polarization degree, and can account for its general anticorrelation with the flux, but the dispersion of the data requires the presence of intrinsic mechanisms, such as turbulence, shocks, or magnetic reconnection. The WEBT optical data are compared to $γ$-ray data from the Fermi satellite. These are analysed with both fixed and adaptive binning procedures. We show that the strong correlation between optical and $γ$-ray data without measurable delay assumes different slopes in faint and high brightness states, and this is compatible with a scenario where in faint states we mainly see the imprint of the geometrical effects, while in bright states the synchrotron self-Compton process dominates. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 12 pages, 13 figures, in press for MNRAS

Gaia Focused Product Release: Sources from Service Interface Function image analysis -- Half a million new sources in omega Centauri

Authors: Gaia Collaboration, K. Weingrill, A. Mints, J. Castañeda, Z. Kostrzewa-Rutkowska, M. Davidson, F. De Angeli, J. Hernández, F. Torra, M. Ramos-Lerate, C. Babusiaux, M. Biermann, C. Crowley, D. W. Evans, L. Lindegren, J. M. Martín-Fleitas, L. Palaversa, D. Ruz Mieres, K. Tisanić, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, A. Barbier , et al. (378 additional authors not shown)

Abstract: Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This ne… ▽ More Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This new pipeline produced half a million additional Gaia sources in the region of the omega Centauri ($ω$ Cen) cluster, which are published with this Focused Product Release. We discuss the dedicated SIF CF data reduction pipeline, validate its data products, and introduce their Gaia archive table. Our aim is to improve the completeness of the {\it Gaia} source inventory in a very dense region in the sky, $ω$ Cen. An adapted version of {\it Gaia}'s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. We validated the results by comparing them to the public {\it Gaia} DR3 catalogue and external Hubble Space Telescope data. With this Focused Product Release, 526\,587 new sources have been added to the {\it Gaia} catalogue in $ω$ Cen. Apart from positions and brightnesses, the additional catalogue contains parallaxes and proper motions, but no meaningful colour information. While SIF CF source parameters generally have a lower precision than nominal {\it Gaia} sources, in the cluster centre they increase the depth of the combined catalogue by three magnitudes and improve the source density by a factor of ten. This first SIF CF data publication already adds great value to the {\it Gaia} catalogue. It demonstrates what to expect for the fourth {\it Gaia} catalogue, which will contain additional sources for all nine SIF CF regions. △ Less

Submitted 8 November, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Journal ref: A&A 680, A35 (2023)

Gaia Focused Product Release: A catalogue of sources around quasars to search for strongly lensed quasars

Authors: Gaia Collaboration, A. Krone-Martins, C. Ducourant, L. Galluccio, L. Delchambre, I. Oreshina-Slezak, R. Teixeira, J. Braine, J. -F. Le Campion, F. Mignard, W. Roux, A. Blazere, L. Pegoraro, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, A. Barbier, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, R. Guerra , et al. (376 additional authors not shown)

Abstract: Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those ex… ▽ More Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 35 pages, 60 figures, accepted for publication by Astronomy and Astrophysics

Journal ref: A&A 685, A130 (2024)

Gaia Focused Product Release: Radial velocity time series of long-period variables

Authors: Gaia Collaboration, Gaia Collaboration, M. Trabucchi, N. Mowlavi, T. Lebzelter, I. Lecoeur-Taibi, M. Audard, L. Eyer, P. García-Lario, P. Gavras, B. Holl, G. Jevardat de Fombelle, K. Nienartowicz, L. Rimoldini, P. Sartoretti, R. Blomme, Y. Frémat, O. Marchal, Y. Damerdji, A. G. A. Brown, A. Guerrier, P. Panuzzo, D. Katz, G. M. Seabroke, K. Benson , et al. (382 additional authors not shown)

Abstract: The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the… ▽ More The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged) △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 36 pages, 38 figures

Multi-year characterisation of the broad-band emission from the intermittent extreme BL Lac 1ES~2344+514

Authors: H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland, O. Blanch , et al. (210 additional authors not shown)

Abstract: The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, f… ▽ More The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $ν_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $ν_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 682, A114 (2024)

Variability and evolution of the optical polarization of a sample of gamma-ray blazars

Authors: J. Otero-Santos, J. A. Acosta-Pulido, J. Becerra González, C. M. Raiteri, M. I. Carnerero, N. Castro Segura, O. González-Martín, A. Luashvili

Abstract: We present a polarization variability analysis of a sample of 26 $γ$-ray blazars monitored by the Steward Observatory between 2008 and 2018 in the optical band. We investigate the properties and long-term variability of their optical polarization, searching for differences between blazar types. We observe that BL Lac objects are typically less polarized and less variable than flat spectrum radio q… ▽ More We present a polarization variability analysis of a sample of 26 $γ$-ray blazars monitored by the Steward Observatory between 2008 and 2018 in the optical band. We investigate the properties and long-term variability of their optical polarization, searching for differences between blazar types. We observe that BL Lac objects are typically less polarized and less variable than flat spectrum radio quasars (FSRQs). Moreover, BL Lacs display a distribution of their polarization angle typically oriented in a preferential direction, contrary to the rather random distribution of FSRQs. For the latter blazar type, as well as those sources showing a bright stellar emission, we take into account the depolarizing effect introduced by the broad line region and the host galaxy on the measured polarization degree. In this sample we also observe that BL Lacs present an uncorrelated evolution of the flux and the polarization. Contrary, FSRQs show a correlation before the depolarization correction, that is lost however after considering this effect. In addition, we study the behaviour of the polarization angle, searching for angle rotations in its long-term evolution. We derive that the FSRQs studied here show rotations more frequently than BL Lac objects by a factor $\sim$1.5. During these periods we also observe a systematic decrease of the polarization fraction, as well as a marginal flux increase, not significant however to connect rotations with optical flares. We interpret these results within the extended shock-in-jet scenario, able to explain the overall features observed here for the polarization of the blazar sample. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS

Gaia Data Release 3: All-sky classification of 12.4 million variable sources into 25 classes

Authors: Lorenzo Rimoldini, Berry Holl, Panagiotis Gavras, Marc Audard, Joris De Ridder, Nami Mowlavi, Krzysztof Nienartowicz, Grégory Jevardat de Fombelle, Isabelle Lecoeur-Taïbi, Lea Karbevska, Dafydd W. Evans, Péter Ábrahám, Maria I. Carnerero, Gisella Clementini, Elisa Distefano, Alessia Garofalo, Pedro García-Lario, Roy Gomel, Sergei A. Klioner, Katarzyna Kruszyńska, Alessandro C. Lanzafame, Thomas Lebzelter, Gábor Marton, Tsevi Mazeh, Roberto Molinaro , et al. (9 additional authors not shown)

Abstract: 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 obje… ▽ More 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. △ Less

Submitted 7 March, 2023; v1 submitted 30 November, 2022; originally announced November 2022.

Comments: 106 pages, 513 figure panels, 4 tables. Accepted by Astronomy & Astrophysics (in press)

Journal ref: A&A 674, A14 (2023)

Multiwavelength periodicity search in a sample of $γ$-ray bright blazars

Authors: J. Otero-Santos, P. Peñil, J. A. Acosta-Pulido, J. Becerra González, C. M. Raiteri, M. I. Carnerero, M. Villata

Abstract: We present the results of a long-term periodicity search in a sample of $γ$-ray blazars within a multiwavelength context. These blazars have been selected from the Steward Observatory sample as part of its optical monitoring program between 2008 and 2018. We study 15 sources with a temporal coverage in their optical total and polarized emission sufficiently large ($>9$ years) to perform a reliable… ▽ More We present the results of a long-term periodicity search in a sample of $γ$-ray blazars within a multiwavelength context. These blazars have been selected from the Steward Observatory sample as part of its optical monitoring program between 2008 and 2018. We study 15 sources with a temporal coverage in their optical total and polarized emission sufficiently large ($>9$ years) to perform a reliable long-term periodicity analysis. We collect data from several observatories to extend the coverage, enabling the search of longer periods. In addition, data are also gathered in the high-energy ($E>100$ MeV) $γ$-ray band from the $\textit{Fermi}$ Large Area Telescope; and in the 15-GHz radio band from the Owens Valley Radio Observatory. We identify 5 promising candidates to host quasi-periodic emission, AO 0235+164, PKS 1222+216, Mrk 501, BL Lacertae and 1ES 2344+514 with periods in one or more bands and statistical significances $\sim$3$σ$ after trial factor correction. AO 0235+164 shows a period of $\sim$8.2 years in the R band; PKS 1222+216 has a quasi-periodic modulation in its total and polarized optical emission of $\sim$1.6 years; Mrk 501 displays a $\sim$5-year quasi-periodicity in optical and radio wavelengths; BL Lacertae presents a period of $\sim$1.8 years in its polarized emission; and 1ES 2344+514 shows a hint of a $\sim$5.5-year period in its optical R band. We interpret these results in the framework of the most common models and scenarios, namely the presence of a binary supermassive black hole system; or geometrical effects like helical or precessing jets. △ Less

Submitted 11 November, 2022; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS

Multi-messenger characterization of Mrk 501 during historically low X-ray and $γ$-ray activity

Authors: MAGIC collaboration, H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder , et al. (300 additional authors not shown)

Abstract: We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the… ▽ More We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the highest occurring at X-rays and very-high-energy (VHE) $γ$-rays. A significant correlation ($>$3$σ$) between X-rays and VHE $γ$-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between Swift-XRT and Fermi-LAT. We additionally find correlations between high-energy $γ$-rays and radio, with the radio lagging by more than 100 days, placing the $γ$-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE $γ$-rays from mid-2017 to mid-2019 with a stable VHE flux ($>$0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2-year-long low-state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock. △ Less

Submitted 5 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

Journal ref: ApJS 266 37 (2023)

Rubin Observatory LSST Transients and Variable Stars Roadmap

Authors: Kelly M. Hambleton, Federica B. Bianco, Rachel Street, Keaton Bell, David Buckley, Melissa Graham, Nina Hernitschek, Michael B. Lund, Elena Mason, Joshua Pepper, Andrej Prsa, Markus Rabus, Claudia M. Raiteri, Robert Szabo, Paula Szkody, Igor Andreoni, Simone Antoniucci, Barbara Balmaverde, Eric Bellm, Rosaria Bonito, Giuseppe Bono, Maria Teresa Botticella, Enzo Brocato, Katja Bucar Bricman, Enrico Cappellaro , et al. (57 additional authors not shown)

Abstract: 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 T… ▽ More 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. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 202 pages (in book format) 34 figures plus chapter heading figures (13)

Gaia Data Release 3: Summary of the content and survey properties

Authors: Gaia Collaboration, A. Vallenari, A. G. A. Brown, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, L. Eyer, R. Guerra, A. Hutton, C. Jordi, S. A. Klioner, U. L. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, C. Soubiran , et al. (431 additional authors not shown)

Abstract: We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photom… ▽ More We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photometry in the G, G$_{BP}$, and G$_{RP}$ pass-bands already present in the Early Third Data Release. GDR3 introduces an impressive wealth of new data products. More than 33 million objects in the ranges $G_{rvs} < 14$ and $3100 <T_{eff} <14500 $, have new determinations of their mean radial velocities based on data collected by Gaia. We provide G$_{rvs}$ magnitudes for most sources with radial velocities, and a line broadening parameter is listed for a subset of these. Mean Gaia spectra are made available to the community. The GDR3 catalogue includes about 1 million mean spectra from the radial velocity spectrometer, and about 220 million low-resolution blue and red prism photometer BPRP mean spectra. The results of the analysis of epoch photometry are provided for some 10 million sources across 24 variability types. GDR3 includes astrophysical parameters and source class probabilities for about 470 million and 1500 million sources, respectively, including stars, galaxies, and quasars. Orbital elements and trend parameters are provided for some $800\,000$ astrometric, spectroscopic and eclipsing binaries. More than $150\,000$ Solar System objects, including new discoveries, with preliminary orbital solutions and individual epoch observations are part of this release. Reflectance spectra derived from the epoch BPRP spectral data are published for about 60\,000 asteroids. Finally, an additional data set is provided, namely the Gaia Andromeda Photometric Survey (abridged) △ Less

Submitted 30 July, 2022; originally announced August 2022.

Comments: 23 pages, 2 figures

Gaia Data Release 3: The first Gaia catalogue of variable AGN

Authors: Maria I. Carnerero, Claudia M. Raiteri, Lorenzo Rimoldini, Deborah Busonero, Enrico Licata, Nami Mowlavi, Isabelle Lecoeur-Taïbi, Marc Audard, Berry Holl, Panagiotis Gavras, Krzysztof Nienartowicz, Grégory Jevardat de Fombelle, Ruth Carballo, Gisella Clementini, Ludovic Delchambre, Sergei Klioner, Mario G. Lattanzi, Laurent Eyer

Abstract: 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 A… ▽ More 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. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 19 pages, 31 figures, 2 table. This paper is part of Gaia Data Release 3 (DR3). In press for A&A

Journal ref: A&A 674, A24 (2023)

Gaia Data Release 3: Reflectance spectra of Solar System small bodies

Authors: Gaia Collaboration, L. Galluccio, M. Delbo, F. De Angeli, T. Pauwels, P. Tanga, F. Mignard, A. Cellino, A. G. A. Brown, K. Muinonen, A. Penttila, S. Jordan, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, L. Eyer, R. Guerra, A. Hutton, C. Jordi , et al. (422 additional authors not shown)

Abstract: The Gaia mission of the European Space Agency (ESA) has been routinely observing Solar System objects (SSOs) since the beginning of its operations in August 2014. The Gaia data release three (DR3) includes, for the first time, the mean reflectance spectra of a selected sample of 60 518 SSOs, primarily asteroids, observed between August 5, 2014, and May 28, 2017. Each reflectance spectrum was deriv… ▽ More The Gaia mission of the European Space Agency (ESA) has been routinely observing Solar System objects (SSOs) since the beginning of its operations in August 2014. The Gaia data release three (DR3) includes, for the first time, the mean reflectance spectra of a selected sample of 60 518 SSOs, primarily asteroids, observed between August 5, 2014, and May 28, 2017. Each reflectance spectrum was derived from measurements obtained by means of the Blue and Red photometers (BP/RP), which were binned in 16 discrete wavelength bands. We describe the processing of the Gaia spectral data of SSOs, explaining both the criteria used to select the subset of asteroid spectra published in Gaia DR3, and the different steps of our internal validation procedures. In order to further assess the quality of Gaia SSO reflectance spectra, we carried out external validation against SSO reflectance spectra obtained from ground-based and space-borne telescopes and available in the literature. For each selected SSO, an epoch reflectance was computed by dividing the calibrated spectrum observed by the BP/RP at each transit on the focal plane by the mean spectrum of a solar analogue. The latter was obtained by averaging the Gaia spectral measurements of a selected sample of stars known to have very similar spectra to that of the Sun. Finally, a mean of the epoch reflectance spectra was calculated in 16 spectral bands for each SSO. The agreement between Gaia mean reflectance spectra and those available in the literature is good for bright SSOs, regardless of their taxonomic spectral class. We identify an increase in the spectral slope of S-type SSOs with increasing phase angle. Moreover, we show that the spectral slope increases and the depth of the 1 um absorption band decreases for increasing ages of S-type asteroid families. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 30 pages, 26 figures

Gaia Data Release 3. Summary of the variability processing and analysis

Authors: L. Eyer, M. Audard, B. Holl, L. Rimoldini, M. I. Carnerero, G. Clementini, J. De Ridder, E. Distefano, D. W. Evans, P. Gavras, R. Gomel, T. Lebzelter, G. Marton, N. Mowlavi, A. Panahi, V. Ripepi, L. Wyrzykowski, K. Nienartowicz, G. Jevardat de Fombelle, I. Lecoeur-Taibi, L. Rohrbasser, M. Riello, P. Garcia-Lario, A. C. Lanzafame, T. Mazeh , et al. (33 additional authors not shown)

Abstract: Context. Gaia has been in operations since 2014. The third Gaia data release expands from the early data release (EDR3) in 2020 by providing 34 months of multi-epoch observations that allowed us to probe, characterise and classify systematically celestial variable phenomena. Aims. We present a summary of the variability processing and analysis of the photometric and spectroscopic time series of… ▽ More Context. Gaia has been in operations since 2014. The third Gaia data release expands from the early data release (EDR3) in 2020 by providing 34 months of multi-epoch observations that allowed us to probe, characterise and classify systematically celestial variable phenomena. Aims. We present a summary of the variability processing and analysis of the photometric and spectroscopic time series of 1.8 billion sources done for Gaia DR3. Methods. We used statistical and Machine Learning methods to characterise and classify the variable sources. Training sets were built from a global revision of major published variable star catalogues. For a subset of classes, specific detailed studies were conducted to confirm their class membership and to derive parameters that are adapted to the peculiarity of the considered class. Results. In total, 10.5 million objects are identified as variable in Gaia DR3 and have associated time series in G, GBP, and GRP and, in some cases, radial velocity time series. The DR3 variable sources subdivide into 9.5 million variable stars and 1 million Active Galactic Nuclei/Quasars. In addition, supervised classification identified 2.5 million galaxies thanks to spurious variability induced by the extent of these objects. The variability analysis output in the DR3 archive amounts to 17 tables containing a total of 365 parameters. We publish 35 types and sub-types of variable objects. For 11 variable types, additional specific object parameters are published. An overview of the estimated completeness and contamination of most variability classes is provided. Conclusions. Thanks to Gaia we present the largest whole-sky variability analysis based on coherent photometric, astrometric, and spectroscopic data. Later Gaia data releases will more than double the span of time series and the number of observations, thus allowing for an even richer catalogue in the future. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 18 pages, 12 figures, submitted to Astronomy & Astrophysics

Gaia Data Release 3: Mapping the asymmetric disc of the Milky Way

Authors: Gaia Collaboration, R. Drimmel, M. Romero-Gomez, L. Chemin, P. Ramos, E. Poggio, V. Ripepi, R. Andrae, R. Blomme, T. Cantat-Gaudin, A. Castro-Ginard, G. Clementini, F. Figueras, M. Fouesneau, Y. Fremat, K. Jardine, S. Khanna, A. Lobel, D. J. Marshall, T. Muraveva, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou , et al. (431 additional authors not shown)

Abstract: With the most recent Gaia data release the number of sources with complete 6D phase space information (position and velocity) has increased to well over 33 million stars, while stellar astrophysical parameters are provided for more than 470 million sources, in addition to the identification of over 11 million variable stars. Using the astrophysical parameters and variability classifications provid… ▽ More With the most recent Gaia data release the number of sources with complete 6D phase space information (position and velocity) has increased to well over 33 million stars, while stellar astrophysical parameters are provided for more than 470 million sources, in addition to the identification of over 11 million variable stars. Using the astrophysical parameters and variability classifications provided in Gaia DR3, we select various stellar populations to explore and identify non-axisymmetric features in the disc of the Milky Way in both configuration and velocity space. Using more about 580 thousand sources identified as hot OB stars, together with 988 known open clusters younger than 100 million years, we map the spiral structure associated with star formation 4-5 kpc from the Sun. We select over 2800 Classical Cepheids younger than 200 million years, which show spiral features extending as far as 10 kpc from the Sun in the outer disc. We also identify more than 8.7 million sources on the red giant branch (RGB), of which 5.7 million have line-of-sight velocities, allowing the velocity field of the Milky Way to be mapped as far as 8 kpc from the Sun, including the inner disc. The spiral structure revealed by the young populations is consistent with recent results using Gaia EDR3 astrometry and source lists based on near infrared photometry, showing the Local (Orion) arm to be at least 8 kpc long, and an outer arm consistent with what is seen in HI surveys, which seems to be a continuation of the Perseus arm into the third quadrant. Meanwhile, the subset of RGB stars with velocities clearly reveals the large scale kinematic signature of the bar in the inner disc, as well as evidence of streaming motions in the outer disc that might be associated with spiral arms or bar resonances. (abridged) △ Less

Submitted 5 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 35 pages, 27 figures, accepted for publication in A&A special Gaia DR3 issue. V2: abstract completed. V3: complete author list and link to data: https://drive.google.com/drive/u/1/folders/1yOJPjYmM7QK5XVsqaiSOTuwDQNti2LlZ

Journal ref: A&A 674, A37 (2023)

Gaia Data Release 3: Pulsations in main sequence OBAF-type stars

Authors: Gaia Collaboration, J. De Ridder, V. Ripepi, C. Aerts, L. Palaversa, L. Eyer, B. Holl, M. Audard, L. Rimoldini, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, R. Guerra, A. Hutton, C. Jordi, S. A. Klioner, U. L. Lammers, L. Lindegren , et al. (423 additional authors not shown)

Abstract: The third Gaia data release provides photometric time series covering 34 months for about 10 million stars. For many of those stars, a characterisation in Fourier space and their variability classification are also provided. This paper focuses on intermediate- to high-mass (IHM) main sequence pulsators M >= 1.3 Msun) of spectral types O, B, A, or F, known as beta Cep, slowly pulsating B (SPB), del… ▽ More The third Gaia data release provides photometric time series covering 34 months for about 10 million stars. For many of those stars, a characterisation in Fourier space and their variability classification are also provided. This paper focuses on intermediate- to high-mass (IHM) main sequence pulsators M >= 1.3 Msun) of spectral types O, B, A, or F, known as beta Cep, slowly pulsating B (SPB), delta Sct, and gamma Dor stars. These stars are often multi-periodic and display low amplitudes, making them challenging targets to analyse with sparse time series. All datasets used in this analysis are part of the Gaia DR3 data release. The photometric time series were used to perform a Fourier analysis, while the global astrophysical parameters necessary for the empirical instability strips were taken from the Gaia DR3 gspphot tables, and the vsini data were taken from the Gaia DR3 esphs tables. We show that for nearby OBAF-type pulsators, the Gaia DR3 data are precise and accurate enough to pinpoint them in the Hertzsprung-Russell diagram. We find empirical instability strips covering broader regions than theoretically predicted. In particular, our study reveals the presence of fast rotating gravity-mode pulsators outside the strips, as well as the co-existence of rotationally modulated variables inside the strips as reported before in the literature. We derive an extensive period-luminosity relation for delta Sct stars and provide evidence that the relation features different regimes depending on the oscillation period. Finally, we demonstrate how stellar rotation attenuates the amplitude of the dominant oscillation mode of delta Sct stars. △ Less

Submitted 16 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Journal ref: A&A 674, A36 (2023)

Gaia Data Release 3: A Golden Sample of Astrophysical Parameters

Authors: Gaia Collaboration, O. L. Creevey, L. M. Sarro, A. Lobel, E. Pancino, R. Andrae, R. L. Smart, G. Clementini, U. Heiter, A. J. Korn, M. Fouesneau, Y. Frémat, F. De Angeli, A. Vallenari, D. L. Harrison, F. Thévenin, C. Reylé, R. Sordo, A. Garofalo, A. G. A. Brown, L. Eyer, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux , et al. (423 additional authors not shown)

Abstract: Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for a half billion stars. In this work we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. We query the astrophysical parameter tables along with other tables in Gaia DR3 to derive the samples… ▽ More Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for a half billion stars. In this work we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. We query the astrophysical parameter tables along with other tables in Gaia DR3 to derive the samples of the stars of interest. We validate our results by using the Gaia catalogue itself and by comparison with external data. We have produced six homogeneous samples of stars with high quality astrophysical parameters across the HR diagram for the community to exploit. We first focus on three samples that span a large parameter space: young massive disk stars (~3M), FGKM spectral type stars (~3M), and UCDs (~20K). We provide these sources along with additional information (either a flag or complementary parameters) as tables that are made available in the Gaia archive. We furthermore identify 15740 bone fide carbon stars, 5863 solar-analogues, and provide the first homogeneous set of stellar parameters of the Spectro Photometric Standard Stars. We use a subset of the OBA sample to illustrate its usefulness to analyse the Milky Way rotation curve. We then use the properties of the FGKM stars to analyse known exoplanet systems. We also analyse the ages of some unseen UCD-companions to the FGKM stars. We additionally predict the colours of the Sun in various passbands (Gaia, 2MASS, WISE) using the solar-analogue sample. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 35 pages, (incl 6 pages references, acknowledgements, affiliations), 37 figures, A&A accepted

Journal ref: A&A 674, A39 (2023)

Gaia Data Release 3: The extragalactic content

Authors: Gaia Collaboration, C. A. L. Bailer-Jones, D. Teyssier, L. Delchambre, C. Ducourant, D. Garabato, D. Hatzidimitriou, S. A. Klioner, L. Rimoldini, I. Bellas-Velidis, R. Carballo, M. I. Carnerero, C. Diener, M. Fouesneau, L. Galluccio, P. Gavras, A. Krone-Martins, C. M. Raiteri, R. Teixeira, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux , et al. (422 additional authors not shown)

Abstract: The Gaia Galactic survey mission is designed and optimized to obtain astrometry, photometry, and spectroscopy of nearly two billion stars in our Galaxy. Yet as an all-sky multi-epoch survey, Gaia also observes several million extragalactic objects down to a magnitude of G~21 mag. Due to the nature of the Gaia onboard selection algorithms, these are mostly point-source-like objects. Using data prov… ▽ More The Gaia Galactic survey mission is designed and optimized to obtain astrometry, photometry, and spectroscopy of nearly two billion stars in our Galaxy. Yet as an all-sky multi-epoch survey, Gaia also observes several million extragalactic objects down to a magnitude of G~21 mag. Due to the nature of the Gaia onboard selection algorithms, these are mostly point-source-like objects. Using data provided by the satellite, we have identified quasar and galaxy candidates via supervised machine learning methods, and estimate their redshifts using the low resolution BP/RP spectra. We further characterise the surface brightness profiles of host galaxies of quasars and of galaxies from pre-defined input lists. Here we give an overview of the processing of extragalactic objects, describe the data products in Gaia DR3, and analyse their properties. Two integrated tables contain the main results for a high completeness, but low purity (50-70%), set of 6.6 million candidate quasars and 4.8 million candidate galaxies. We provide queries that select purer sub-samples of these containing 1.9 million probable quasars and 2.9 million probable galaxies (both 95% purity). We also use high quality BP/RP spectra of 43 thousand high probability quasars over the redshift range 0.05-4.36 to construct a composite quasar spectrum spanning restframe wavelengths from 72-100 nm. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: Accepted to A&A

Gaia Data Release 3: Stellar multiplicity, a teaser for the hidden treasure

Authors: Gaia Collaboration, F. Arenou, C. Babusiaux, M. A. Barstow, S. Faigler, A. Jorissen, P. Kervella, T. Mazeh, N. Mowlavi, P. Panuzzo, J. Sahlmann, S. Shahaf, A. Sozzetti, N. Bauchet, Y. Damerdji, P. Gavras, P. Giacobbe, E. Gosset, J. -L. Halbwachs, B. Holl, M. G. Lattanzi, N. Leclerc, T. Morel, D. Pourbaix, P. Re Fiorentin , et al. (425 additional authors not shown)

Abstract: The Gaia DR3 Catalogue contains for the first time about eight hundred thousand solutions with either orbital elements or trend parameters for astrometric, spectroscopic and eclipsing binaries, and combinations of them. This paper aims to illustrate the huge potential of this large non-single star catalogue. Using the orbital solutions together with models of the binaries, a catalogue of tens of t… ▽ More The Gaia DR3 Catalogue contains for the first time about eight hundred thousand solutions with either orbital elements or trend parameters for astrometric, spectroscopic and eclipsing binaries, and combinations of them. This paper aims to illustrate the huge potential of this large non-single star catalogue. Using the orbital solutions together with models of the binaries, a catalogue of tens of thousands of stellar masses, or lower limits, partly together with consistent flux ratios, has been built. Properties concerning the completeness of the binary catalogues are discussed, statistical features of the orbital elements are explained and a comparison with other catalogues is performed. Illustrative applications are proposed for binaries across the H-R diagram. The binarity is studied in the RGB/AGB and a search for genuine SB1 among long-period variables is performed. The discovery of new EL CVn systems illustrates the potential of combining variability and binarity catalogues. Potential compact object companions are presented, mainly white dwarf companions or double degenerates, but one candidate neutron star is also presented. Towards the bottom of the main sequence, the orbits of previously-suspected binary ultracool dwarfs are determined and new candidate binaries are discovered. The long awaited contribution of Gaia to the analysis of the substellar regime shows the brown dwarf desert around solar-type stars using true, rather than minimum, masses, and provides new important constraints on the occurrence rates of substellar companions to M dwarfs. Several dozen new exoplanets are proposed, including two with validated orbital solutions and one super-Jupiter orbiting a white dwarf, all being candidates requiring confirmation. Beside binarity, higher order multiple systems are also found. △ Less

Submitted 11 June, 2022; originally announced June 2022.

Comments: 60 pages, 60 figures. Accepted for publication in Astronomy & Astrophysics (2022-06-09). The catalogue of binary masses is available for download from the ESA Gaia DR3 Archive and will be available from the CDS/VizieR service

Journal ref: A&A 674, A34 (2023)

Gaia Data Release 3: Chemical cartography of the Milky Way

Authors: Gaia Collaboration, A. Recio-Blanco, G. Kordopatis, P. de Laverny, P. A. Palicio, A. Spagna, L. Spina, D. Katz, P. Re Fiorentin, E. Poggio, P. J. McMillan, A. Vallenari, M. G. Lattanzi, G. M. Seabroke, L. Casamiquela, A. Bragaglia, T. Antoja, C. A. L. Bailer-Jones, R. Andrae, M. Fouesneau, M. Cropper, T. Cantat-Gaudin, U. Heiter, A. Bijaoui, A. G. A. Brown , et al. (425 additional authors not shown)

Abstract: Gaia DR3 opens a new era of all-sky spectral analysis of stellar populations thanks to the nearly 5.6 million stars observed by the RVS and parametrised by the GSP-spec module. The all-sky Gaia chemical cartography allows a powerful and precise chemo-dynamical view of the Milky Way with unprecedented spatial coverage and statistical robustness. First, it reveals the strong vertical symmetry of the… ▽ More Gaia DR3 opens a new era of all-sky spectral analysis of stellar populations thanks to the nearly 5.6 million stars observed by the RVS and parametrised by the GSP-spec module. The all-sky Gaia chemical cartography allows a powerful and precise chemo-dynamical view of the Milky Way with unprecedented spatial coverage and statistical robustness. First, it reveals the strong vertical symmetry of the Galaxy and the flared structure of the disc. Second, the observed kinematic disturbances of the disc -- seen as phase space correlations -- and kinematic or orbital substructures are associated with chemical patterns that favour stars with enhanced metallicities and lower [alpha/Fe] abundance ratios compared to the median values in the radial distributions. This is detected both for young objects that trace the spiral arms and older populations. Several alpha, iron-peak elements and at least one heavy element trace the thin and thick disc properties in the solar cylinder. Third, young disc stars show a recent chemical impoverishment in several elements. Fourth, the largest chemo-dynamical sample of open clusters analysed so far shows a steepening of the radial metallicity gradient with age, which is also observed in the young field population. Finally, the Gaia chemical data have the required coverage and precision to unveil galaxy accretion debris and heated disc stars on halo orbits through their [alpha/Fe] ratio, and to allow the study of the chemo-dynamical properties of globular clusters. Gaia DR3 chemo-dynamical diagnostics open new horizons before the era of ground-based wide-field spectroscopic surveys. They unveil a complex Milky Way that is the outcome of an eventful evolution, shaping it to the present day (abridged). △ Less

Submitted 11 June, 2022; originally announced June 2022.

Comments: Astronomy and Astrophysics (accepted, in press)

Journal ref: A&A 674, A38 (2023)

Gaia Early Data Release 3: The celestial reference frame (Gaia-CRF3)

Authors: Gaia Collaboration, S. A. Klioner, L. Lindegren, F. Mignard, J. Hernández, M. Ramos-Lerate, U. Bastian, M. Biermann, A. Bombrun, A. de Torres, E. Gerlach, R. Geyer, T. Hilger, D. Hobbs, U. L. Lammers, P. J. McMillan, H. Steidelmüller, D. Teyssier, C. M. Raiteri, S. Bartolomé, M. Bernet, J. Castañeda, M. Clotet, M. Davidson, C. Fabricius , et al. (426 additional authors not shown)

Abstract: Gaia-CRF3 is the celestial reference frame for positions and proper motions in the third release of data from the Gaia mission, Gaia DR3 (and for the early third release, Gaia EDR3, which contains identical astrometric results). The reference frame is defined by the positions and proper motions at epoch 2016.0 for a specific set of extragalactic sources in the (E)DR3 catalogue. We describe the c… ▽ More Gaia-CRF3 is the celestial reference frame for positions and proper motions in the third release of data from the Gaia mission, Gaia DR3 (and for the early third release, Gaia EDR3, which contains identical astrometric results). The reference frame is defined by the positions and proper motions at epoch 2016.0 for a specific set of extragalactic sources in the (E)DR3 catalogue. We describe the construction of Gaia-CRF3, and its properties in terms of the distributions in magnitude, colour, and astrometric quality. Compact extragalactic sources in Gaia DR3 were identified by positional cross-matching with 17 external catalogues of quasars (QSO) and active galactic nuclei (AGN), followed by astrometric filtering designed to remove stellar contaminants. Selecting a clean sample was favoured over including a higher number of extragalactic sources. For the final sample, the random and systematic errors in the proper motions are analysed, as well as the radio-optical offsets in position for sources in the third realisation of the International Celestial Reference Frame (ICRF3). The Gaia-CRF3 comprises about 1.6 million QSO-like sources, of which 1.2 million have five-parameter astrometric solutions in Gaia DR3 and 0.4 million have six-parameter solutions. The sources span the magnitude range G = 13 to 21 with a peak density at 20.6 mag, at which the typical positional uncertainty is about 1 mas. The proper motions show systematic errors on the level of 12 $μ$as yr${}^{-1}$ on angular scales greater than 15 deg. For the 3142 optical counterparts of ICRF3 sources in the S/X frequency bands, the median offset from the radio positions is about 0.5 mas, but exceeds 4 mas in either coordinate for 127 sources. We outline the future of the Gaia-CRF in the next Gaia data releases. △ Less

Submitted 30 October, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Journal ref: A&A 667, A148 (2022)

A statistical study of the optical spectral variability in gamma-ray blazars

Authors: J. Otero-Santos, J. A. Acosta-Pulido, J. Becerra González, A. Luashvili, N. Castro Segura, O. González-Martín, C. M. Raiteri, M. I. Carnerero

Abstract: Blazars optical emission is generally dominated by relativistic jets, although the host galaxy, accretion disk and broad line region (BLR) may also contribute significantly. Disentangling their contributions has been challenging for years due to the dominance of the jet. To quantify the contributions to the spectral variability, we use the statistical technique for dimensionality reduction Non-Neg… ▽ More Blazars optical emission is generally dominated by relativistic jets, although the host galaxy, accretion disk and broad line region (BLR) may also contribute significantly. Disentangling their contributions has been challenging for years due to the dominance of the jet. To quantify the contributions to the spectral variability, we use the statistical technique for dimensionality reduction Non-Negative Matrix Factorization on a spectroscopic data set of 26 $γ$-ray blazars. This technique allows to model large numbers of spectra in terms of a reduced number of components.We use a priori knowledge to obtain components associated to meaningful physical processes. The sources are classified according to their optical spectrum as host-galaxy dominated BL Lac objects (BL Lacs), BL Lacs, or Flat Spectrum Radio Quasars (FSRQs). Host-galaxy sources show less variability, as expected, and bluer-when-brighter trends, as the other BL Lacs. For FSRQs, more complicated colour-flux behaviours are observed: redder-when-brighter for low states saturating above a certain level and, in some cases, turning to bluer-when-brighter. We are able to reproduce the variability observed during 10 years using only 2 to 4 components, depending on the type. The simplest scenario corresponds to host-galaxy blazars, whose spectra are reconstructed using the stellar population and a power law for the jet. BL Lac spectra are reproduced using from 2 to 4 power laws. Different components can be associated to acceleration/cooling processes taking place in the jet. The reconstruction of FSRQs also incorporates a QSO-like component to account for the BLR, plus a very steep power law, associated to the accretion disk. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS

Blazar variability with the Vera C. Rubin Legacy Survey of Space and Time (LSST)

Authors: Claudia M. Raiteri, Maria I. Carnerero, Barbara Balmaverde, Eric C. Bellm, William Clarkson, Filippo D'Ammando, Maurizio Paolillo, Gordon T. Richards, Massimo Villata, Peter Yoachim, Ilsang Yoon

Abstract: 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 va… ▽ More 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. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 18 pages, 12 figures, accepted for publication in the ApJS Rubin Cadence Focus Issue

Investigation of the correlation patterns and the Compton dominance variability of Mrk 421 in 2017

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, O. Blanch, Ž. Bošnjak, G. Busetto, R. Carosi , et al. (263 additional authors not shown)

Abstract: We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicin… ▽ More We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina and Metsähovi. 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. 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 different patterns, partly due to substantial changes in the Compton dominance during 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. Our multi-band correlation study also hints at an anti-correlation between UV/optical and X-ray at a significance higher than 3 sigmas. A VHE flare observed on 2017 February 4 shows gamma-ray variability on multi-hour timescales, with a factor 10 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. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 38 pages, 23 figures, accepted for publication in Astronomy and Astrophysics. Corresponding authors: Axel Arbet-Engels, David Paneque

Journal ref: A&A 655, A89 (2021)

Gaia Early Data Release 3: The Galactic anticentre

Authors: Gaia Collaboration, T. Antoja, P. McMillan, G. Kordopatis, P. Ramos, A. Helmi, E. Balbinot, T. Cantat-Gaudin, L. Chemin, F. Figueras, C. Jordi, S. Khanna, M. Romero-Gomez, G. Seabroke, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, A. Hutton, F. Jansen , et al. (395 additional authors not shown)

Abstract: We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of the Milky Way structure and evolution. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. We explore the disturbances of the current d… ▽ More We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of the Milky Way structure and evolution. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. We explore the disturbances of the current disc, the spatial and kinematical distributions of early accreted versus in-situ stars, the structures in the outer parts of the disc, and the orbits of open clusters Berkeley 29 and Saurer 1. We find that: i) the dynamics of the Galactic disc are very complex with vertical asymmetries, and new correlations, including a bimodality with disc stars with large angular momentum moving vertically upwards from below the plane, and disc stars with slightly lower angular momentum moving preferentially downwards; ii) we resolve the kinematic substructure (diagonal ridges) in the outer parts of the disc for the first time; iii) the red sequence that has been associated with the proto-Galactic disc that was present at the time of the merger with Gaia-Enceladus-Sausage is currently radially concentrated up to around 14 kpc, while the blue sequence that has been associated with debris of the satellite extends beyond that; iv) there are density structures in the outer disc, both above and below the plane, most probably related to Monoceros, the Anticentre Stream, and TriAnd, for which the Gaia data allow an exhaustive selection of candidate member stars and dynamical study; and v) the open clusters Berkeley~29 and Saurer~1, despite being located at large distances from the Galactic centre, are on nearly circular disc-like orbits. We demonstrate how, once again, the Gaia are crucial for our understanding of the different pieces of our Galaxy and their connection to its global structure and history. △ Less

Submitted 26 April, 2021; v1 submitted 14 January, 2021; originally announced January 2021.

Comments: Gaia EDR3 performance verification paper, version 2 closer to published version in A&A, complete list of authors

Journal ref: A&A 649, A8 (2021)

Gaia Early Data Release 3: The Gaia Catalogue of Nearby Stars

Authors: Gaia Collaboration, R. L. Smart, L. M. Sarro, J. Rybizki, C. Reylé, A. C. Robin, N. C. Hambly, U. Abbas, M. A. Barstow, J. H. J. de Bruijne, B. Bucciarelli, J. M. Carrasco, W. J. Cooper, S. T. Hodgkin, E. Masana, D. Michalik, J. Sahlmann, A. Sozzetti, A. G. A. Brown, A. Vallenari, T. Prusti, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans , et al. (398 additional authors not shown)

Abstract: We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of obj… ▽ More We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of objects within 100\,pc from the full catalogue used selected training sets, machine-learning procedures, astrometric quantities, and solution quality indicators to determine a probability that the astrometric solution is reliable. The training set construction exploited the astrometric data, quality flags, and external photometry. For all candidates we calculated distance posterior probability densities using Bayesian procedures and mock catalogues to define priors. Any object with reliable astrometry and a non-zero probability of being within 100\,pc is included in the catalogue. We have produced a catalogue of \NFINAL\ objects that we estimate contains at least 92\% of stars of stellar type M9 within 100\,pc of the Sun. We estimate that 9\% of the stars in this catalogue probably lie outside 100\,pc, but when the distance probability function is used, a correct treatment of this contamination is possible. We produced luminosity functions with a high signal-to-noise ratio for the main-sequence stars, giants, and white dwarfs. We examined in detail the Hyades cluster, the white dwarf population, and wide-binary systems and produced candidate lists for all three samples. We detected local manifestations of several streams, superclusters, and halo objects, in which we identified 12 members of \G\ Enceladus. We present the first direct parallaxes of five objects in multiple systems within 10\,pc of the Sun. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 45 Pages, 39 figures in main part and 18 in appendix, tables on CDS

Journal ref: A&A 649, A6 (2021)

Gaia Early Data Release 3: Acceleration of the solar system from Gaia astrometry

Authors: Gaia Collaboration, S. A. Klioner, F. Mignard, L. Lindegren, U. Bastian, P. J. McMillan, J. Hernández, D. Hobbs, M. Ramos-Lerate, M. Biermann, A. Bombrun, A. de Torres, E. Gerlach, R. Geyer, T. Hilger, U. Lammers, H. Steidelmüller, C. A. Stephenson, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, O. L. Creevey, D. W. Evans , et al. (392 additional authors not shown)

Abstract: Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar system barycentre with respect to the rest frame of the… ▽ More Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar system barycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of the Gaia astrometric solution. Methods. The effect of the acceleration is obtained as a part of the general expansion of the vector field of proper motions in Vector Spherical Harmonics (VSH). Various versions of the VSH fit and various subsets of the sources are tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with the Gaia astrometric solution are used to get a better idea on possible systematic errors in the estimate. Results. Our best estimate of the acceleration based on Gaia EDR3 is $(2.32 \pm 0.16) \times 10^{-10}$ m s${}^{-2}$ (or $7.33 \pm 0.51$ km s$^{-1}$ Myr${}^{-1}$) towards $α= 269.1^\circ \pm 5.4^\circ$, $δ= -31.6^\circ \pm 4.1^\circ$, corresponding to a proper motion amplitude of $5.05 \pm 0.35$ $μ$as yr${}^{-1}$. This is in good agreement with the acceleration expected from current models of the Galactic gravitational potential. We expect that future Gaia data releases will provide estimates of the acceleration with uncertainties substantially below 0.1 $μ$as yr${}^{-1}$. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: A&A, accepted

Journal ref: A&A 649, A9 (2021)

Gaia Early Data Release 3: Structure and properties of the Magellanic Clouds

Authors: Gaia Collaboration, X. Luri, L. Chemin, G. Clementini, H. E. Delgado, P. J. McMillan, M. Romero-Gómez, E. Balbinot, A. Castro-Ginard, R. Mor, V. Ripepi, L. M. Sarro, M. -R. L. Cioni, C. Fabricius, A. Garofalo, A. Helmi, T. Muraveva, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans , et al. (395 additional authors not shown)

Abstract: We compare the Gaia DR2 and Gaia EDR3 performances in the study of the Magellanic Clouds and show the clear improvements in precision and accuracy in the new release. We also show that the systematics still present in the data make the determination of the 3D geometry of the LMC a difficult endeavour; this is at the very limit of the usefulness of the Gaia EDR3 astrometry, but it may become feasib… ▽ More We compare the Gaia DR2 and Gaia EDR3 performances in the study of the Magellanic Clouds and show the clear improvements in precision and accuracy in the new release. We also show that the systematics still present in the data make the determination of the 3D geometry of the LMC a difficult endeavour; this is at the very limit of the usefulness of the Gaia EDR3 astrometry, but it may become feasible with the use of additional external data. We derive radial and tangential velocity maps and global profiles for the LMC for the several subsamples we defined. To our knowledge, this is the first time that the two planar components of the ordered and random motions are derived for multiple stellar evolutionary phases in a galactic disc outside the Milky Way, showing the differences between younger and older phases. We also analyse the spatial structure and motions in the central region, the bar, and the disc, providing new insights into features and kinematics. Finally, we show that the Gaia EDR3 data allows clearly resolving the Magellanic Bridge, and we trace the density and velocity flow of the stars from the SMC towards the LMC not only globally, but also separately for young and evolved populations. This allows us to confirm an evolved population in the Bridge that is slightly shift from the younger population. Additionally, we were able to study the outskirts of both Magellanic Clouds, in which we detected some well-known features and indications of new ones. △ Less

Submitted 4 January, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: This paper is part of the "demonstration papers" released with Gaia EDR3: https://www.cosmos.esa.int/web/gaia/earlydr3

Journal ref: A&A 649, A7 (2021)

Gaia Early Data Release 3: Summary of the contents and survey properties

Authors: Gaia Collaboration, A. G. A Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, A. Hutton, F. Jansen, C. Jordi, S. A. Klioner, U. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, C. Soubiran, N. A. Walton, F. Arenou , et al. (401 additional authors not shown)

Abstract: We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motio… ▽ More We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motions, and the (G_BP-G_RP) colour are also available. The passbands for G, G_BP, and G_RP are provided as part of the release. For ease of use, the 7 million radial velocities from Gaia DR2 are included in this release, after the removal of a small number of spurious values. New radial velocities will appear as part of Gaia DR3. Finally, Gaia EDR3 represents an updated materialisation of the celestial reference frame (CRF) in the optical, the Gaia-CRF3, which is based solely on extragalactic sources. The creation of the source list for Gaia EDR3 includes enhancements that make it more robust with respect to high proper motion stars, and the disturbing effects of spurious and partially resolved sources. The source list is largely the same as that for Gaia DR2, but it does feature new sources and there are some notable changes. The source list will not change for Gaia DR3. Gaia EDR3 represents a significant advance over Gaia DR2, with parallax precisions increased by 30 percent, proper motion precisions increased by a factor of 2, and the systematic errors in the astrometry suppressed by 30--40 percent for the parallaxes and by a factor ~2.5 for the proper motions. The photometry also features increased precision, but above all much better homogeneity across colour, magnitude, and celestial position. A single passband for G, G_BP, and G_RP is valid over the entire magnitude and colour range, with no systematics above the 1 percent level. △ Less

Submitted 9 June, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: Accepted for A&A Special Issue on Gaia EDR3, 21 pages, 2 figures. This version includes the updates in the erratum (https://doi.org/10.1051/0004-6361/202039657e)

Journal ref: A&A 650, C3 (2021)

The dual nature of blazar fast variability. Space and ground observations of S5 0716+714

Authors: C. M. Raiteri, M. Villata, D. Carosati, E. Benítez, S. O. Kurtanidze, A. C. Gupta, D. O. Mirzaqulov, F. D'Ammando, V. M. Larionov, T. Pursimo, J. A. Acosta-Pulido, G. V. Baida, B. Balmaverde, G. Bonnoli, G. A. Borman, M. I. Carnerero, W. -P. Chen, V. Dhiman, A. Di Maggio, S. A. Ehgamberdiev, D. Hiriart, G. N. Kimeridze, O. M. Kurtanidze, C. S. Lin, J. M. Lopez , et al. (17 additional authors not shown)

Abstract: Blazar S5 0716+714 is well-known for its short-term variability, down to intra-day time-scales. We here present the 2-min cadence optical light curve obtained by the TESS space telescope in 2019 December - 2020 January and analyse the object fast variability with unprecedented sampling. Supporting observations by the Whole Earth Blazar Telescope Collaboration in B, V , R, and I bands allow us to i… ▽ More Blazar S5 0716+714 is well-known for its short-term variability, down to intra-day time-scales. We here present the 2-min cadence optical light curve obtained by the TESS space telescope in 2019 December - 2020 January and analyse the object fast variability with unprecedented sampling. Supporting observations by the Whole Earth Blazar Telescope Collaboration in B, V , R, and I bands allow us to investigate the spectral variability during the TESS pointing. The spectral analysis is further extended in frequency to the UV and X-ray bands with data from the Neil Gehrels Swift Observatory. We develop a new method to unveil the shortest optical variability time-scales. This is based on progressive de-trending of the TESS light curve by means of cubic spline interpolations through the binned fluxes, with decreasing time bins. The de-trended light curves are then analysed with classical tools for time-series analysis (periodogram, auto-correlation and structure functions). The results show that below 3 d there are significant characteristic variability time-scales of about 1.7, 0.5, and 0.2 d. Variability on time-scales < 0.2 d is strongly chromatic and must be ascribed to intrinsic energetic processes involving emitting regions, likely jet sub-structures, with dimension less than about milliparsec. In contrast, flux changes on time-scales > 0.5 d are quasi-achromatic and are probably due to Doppler factor changes of geometric origin. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: 17 pages, 18 figures, accepted by MNRAS

Multi-Wavelength Variability of BL Lacertae Measured with High Time Resolution

Authors: Zachary R. Weaver, K. E. Williamson, S. G. Jorstad, A. P. Marscher, V. M. Larionov, C. M. Raiteri, M. Villata, J. A. Acosta-Pulido, R. Bachev, G. V. Baida, T. J. Balonek, E. Benitez, G. A. Borman, V. Bozhilov, M. I. Carnerero, D. Carosati, W. P. Chen, G. Damljanovic, V. Dhiman, D. J. Dougherty, S. A. Ehgamberdiev, T. S. Grishina, A. C. Gupta, M. Hart, D. Hiriart , et al. (32 additional authors not shown)

Abstract: In an effort to locate the sites of emission at different frequencies and physical processes causing variability in blazar jets, we have obtained high time-resolution observations of BL Lacertae over a wide wavelength range: with the \emph{Transiting Exoplanet Survey Satellite} (TESS) at 6,000-10,000 Å with 2-minute cadence; with the Neil Gehrels \emph{Swift} satellite at optical, UV, and X-ray ba… ▽ More In an effort to locate the sites of emission at different frequencies and physical processes causing variability in blazar jets, we have obtained high time-resolution observations of BL Lacertae over a wide wavelength range: with the \emph{Transiting Exoplanet Survey Satellite} (TESS) at 6,000-10,000 Å with 2-minute cadence; with the Neil Gehrels \emph{Swift} satellite at optical, UV, and X-ray bands; with the Nuclear Spectroscopic Telescope Array at hard X-ray bands; with the \emph{Fermi} Large Area Telescope at $γ$-ray energies; and with the Whole Earth Blazar Telescope for measurement of the optical flux density and polarization. All light curves are correlated, with similar structure on timescales from hours to days. The shortest timescale of variability at optical frequencies observed with TESS is $\sim 0.5$ hr. The most common timescale is $13\pm1$~hr, comparable with the minimum timescale of X-ray variability, 14.5 hr. The multi-wavelength variability properties cannot be explained by a change solely in the Doppler factor of the emitting plasma. The polarization behavior implies that there are both ordered and turbulent components to the magnetic field in the jet. Correlation analysis indicates that the X-ray variations lag behind the $γ$-ray and optical light curves by up to $\sim 0.4$ days. The timescales of variability, cross-frequency lags, and polarization properties can be explained by turbulent plasma that is energized by a shock in the jet and subsequently loses energy to synchrotron and inverse Compton radiation in a magnetic field of strength $\sim3$ G △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 33 pages, 25 figures, 14 tables. Accepted to ApJ

A New Multi-Wavelength Census of Blazars

Authors: A. Paggi, M. Bonato, C. M. Raiteri, M. Villata, G. De Zotti, M. I. Carnerero

Abstract: Context:Blazars are the rarest and most powerful active galactic nuclei, playing a crucial and growing role in today multi-frequency and multi-messenger astrophysics. Current blazar catalogs, however, are incomplete and particularly depleted at low Galactic latitudes. Aims: We aim at augmenting the current blazar census to build a catalog of blazar candidates with homogeneous sky coverage that can… ▽ More Context:Blazars are the rarest and most powerful active galactic nuclei, playing a crucial and growing role in today multi-frequency and multi-messenger astrophysics. Current blazar catalogs, however, are incomplete and particularly depleted at low Galactic latitudes. Aims: We aim at augmenting the current blazar census to build a catalog of blazar candidates with homogeneous sky coverage that can provide candidate counterparts to unassociated gamma-ray sources, sources of high-energy neutrino emission, and ultra-high energy cosmic rays. Methods: Starting from the ALMA Calibrator Catalog we built a catalog of 1580 blazar candidates (ALMA Blazar Candidates, ABC) for which we collect multi-wavelength information. We also compared ABC sources with existing blazar catalogs. Results: The ABC catalogue fills the lack of low Galactic latitude sources in current blazar catalogues. ABC sources are significantly dimmer than known blazars in Gaia g band, and they appear bluer in SDSS and WISE colors. The majority of ABC sources (~ 90%) have optical spectra that classify them as QSO, while the remaining sources resulted galactic objects. ABC sources are similar in X-rays to known blazar, while in gamma-rays they are on average dimmer and softer, indicating a significant contribution of FSRQ sources. Making use of WISE colours, we classified 715 ABC sources as candidate gamma-ray blazar of different classes. Conclusions: We built a new catalogue of 1580 candidate blazars with a rich multi-wavelength data-set, filling the lack of low Galactic latitude sources in current blazar catalogues. This will be particularly important to identify the source population of high energy neutrinos or ultra-high energy cosmic rays. The data collected by the upcoming LSST surveys will provide a key tool to investigate the possible blazar nature of these sources. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 53 pages, 32 figures, 6 tables. Accepted for publication on A&A

Journal ref: A&A 641, A62 (2020)

Unveiling the monster heart: unbeamed properties of blazar 4C 71.07

Authors: C. M. Raiteri, J. A. Acosta Pulido, M. Villata, M. I. Carnerero, P. Romano, S. Vercellone

Abstract: 4C 71.07 is a high-redshift blazar whose optical radiation is dominated by quasar-like nuclear emission. We here present the results of a spectroscopic monitoring of the source to study its unbeamed properties. We obtained 24 optical spectra at the Nordic Optical Telescope (NOT) and William Herschel Telescope (WHT) and 3 near-infrared spectra at the Telescopio Nazionale Galileo (TNG). They show no… ▽ More 4C 71.07 is a high-redshift blazar whose optical radiation is dominated by quasar-like nuclear emission. We here present the results of a spectroscopic monitoring of the source to study its unbeamed properties. We obtained 24 optical spectra at the Nordic Optical Telescope (NOT) and William Herschel Telescope (WHT) and 3 near-infrared spectra at the Telescopio Nazionale Galileo (TNG). They show no evidence of narrow emission lines. The estimate of the systemic redshift from the H$β$ and H$α$ broad emission lines leads to $z_{\rm sys}=2.2130 \pm 0.0004$. Notwithstanding the nearly face-on orientation of the accretion disc, the high-ionization emission lines present large broadening as well as noticeable blueshifts, which increase with the ionizing energy of the corresponding species. This is a clear indication of strong ionized outflows. Line broadening and blueshift appear correlated. We applied scaling relationships to estimate the mass of the supermassive black hole from the Balmer and C IV lines, taking into account the prescriptions to correct for outflow. They give $M_{\rm BH} \sim 2 \times 10^9 \, M_\odot$. We derived an Eddington luminosity $L_{\rm Edd} \sim 2.5 \times 10^{47} \rm \, erg \, s^{-1}$ $\sim L_{\rm disc}$, and a broad line region luminosity $L_{\rm BLR} \sim 1.5 \times 10^{46} \rm \, erg \, s^{-1}$. The line fluxes do not show significant variability in time. In particular, there is no line reaction to the jet flaring activity detected in 2015 October and November. This implies that the jet gives no contribution to the photoionization of the broad line region in the considered period. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: 13 pages, 13 figures, in press for MNRAS

Multiwavelength behaviour of the blazar 3C279: decade-long study from $γ$-ray to radio

Authors: V. M. Larionov, S. G. Jorstad, A. P. Marscher, M. Villata, C. M. Raiteri, P. S. Smith, I. Agudo, S. S. Savchenko, D. A. Morozova, J. A. Acosta-Pulido, M. F. Aller, H. D. Aller, T. S. Andreeva, A. A. Arkharov, R. Bachev, G. Bonnoli, G. A. Borman, V. Bozhilov, P. Calcidese, M. I. Carnerero, D. Carosati, C. Casadio, W. -P. Chen, G. Damljanovic, A. V. Dementyev , et al. (62 additional authors not shown)

Abstract: We report the results of decade-long (2008-2018) $γ$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $γ$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $γ$-ray-optical flux-flux relation changes with activi… ▽ More We report the results of decade-long (2008-2018) $γ$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $γ$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $γ$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $γ$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $τ=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 21 pages, 18 figures, accepted for publication in MNRAS

Journal ref: MNRAS 492 (2020) 3829

Quasi-periodic behaviour in the optical and γ-ray light curves of blazars 3C 66A and B2 1633+38

Authors: J. Otero-Santos, J. A. Acosta-Pulido, J. Becerra González, C. M. Raiteri, V. M. Larionov, P. Peñil, P. S. Smith, C. Ballester Niebla, G. A. Borman, M. I. Carnerero, N. Castro Segura, T. S. Grishina, E. N. Kopatskaya, E. G. Larionova, D. A. Morozova, A. A. Nikiforova, S. S. Savchenko, Yu. V. Troitskaya, I. S. Troitsky, A. A. Vasilyev, M. Villata

Abstract: We report on quasi-periodic variability found in two blazars included in the Steward Observatory Blazar Monitoring data sample: the BL Lac object 3C 66A and the Flat Spectrum Radio Quasar B2 1633+38. We collect optical photometric and polarimetric data in V and R bands of these sources from different observatories: St. Petersburg University, Crimean Astrophysical Observatory, WEBT-GASP, Catalina R… ▽ More We report on quasi-periodic variability found in two blazars included in the Steward Observatory Blazar Monitoring data sample: the BL Lac object 3C 66A and the Flat Spectrum Radio Quasar B2 1633+38. We collect optical photometric and polarimetric data in V and R bands of these sources from different observatories: St. Petersburg University, Crimean Astrophysical Observatory, WEBT-GASP, Catalina Real-Time Transient Survey, Steward Observatory, STELLA Robotic Observatory and Katzman Automatic Imaging Telescope. In addition, an analysis of the $γ$-ray light curves from $\textit{Fermi}$-LAT is included. Three methods are used to search for any periodic behaviour in the data: the Z-transform Discrete Correlation Function, the Lomb-Scargle periodogram and the Weighted Wavelet Z-transform. We find evidences of possible quasi-periodic variability in the optical photometric data of both sources with periods of $\sim$3 years for 3C 66A and $\sim$1.9 years for B2 1633+38, with significances between 3$σ$ and 5$σ$. Only B2 1633+38 shows evidence of this behaviour in the optical polarized data set at a confidence level of 2$σ$-4$σ$. This is the first reported evidence of quasi-periodic behaviour in the optical light curve of B2 1633+38. Also a hint of quasi-periodic behaviour is found in the $γ$-ray light curve of B2 1633+38 with a confidence level $\geqslant$2$σ$, while no periodicity is observed for 3C 66A in this energy range. We propose different jet emission models that could explain the quasi-periodic variability and the differences found between these two sources. △ Less

Submitted 5 February, 2020; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS

Investigating the multiwavelength behaviour of the flat spectrum radio quasar CTA 102 during 2013-2017

Authors: F. D'Ammando, C. M. Raiteri, M. Villata, J. A. Acosta-Pulido, I. Agudo, A. A. Arkharov, R. Bachev, G. V. Baida, E. Benitez, G. A. Borman, W. Boschin, V. Bozhilov, M. S. Butuzova, P. Calcidese, M. I. Carnerero, D. Carosati, C. Casadio, N. Castro-Segura, W. -P. Chen, G. Damljanovic, A. Di Paola, J. Echevarria, N. V. Efimova, Sh. A. Ehgamberdiev, C. Espinosa , et al. (72 additional authors not shown)

Abstract: We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and… ▽ More We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($γ$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $γ$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $γ$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $γ$-ray activity is found. The $γ$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $γ$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission. △ Less

Submitted 11 November, 2019; v1 submitted 8 October, 2019; originally announced October 2019.

Comments: 24 pages, 18 figures, and 4 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The beamed jet and quasar core of the distant blazar 4C 71.07

Authors: C. M. Raiteri, M. Villata, M. I. Carnerero, J. A. Acosta-Pulido, D. O. Mirzaqulov, V. M. Larionov, P. Romano, S. Vercellone

Abstract: The object 4C 71.07 is a high-redshift blazar whose spectral energy distribution shows a prominent big blue bump and a strong Compton dominance. We present the results of a two-year multiwavelength campaign led by the Whole Earth Blazar Telescope (WEBT) to study both the quasar core and the beamed jet of this source. The WEBT data are complemented by ultraviolet and X-ray data from Swift, and by g… ▽ More The object 4C 71.07 is a high-redshift blazar whose spectral energy distribution shows a prominent big blue bump and a strong Compton dominance. We present the results of a two-year multiwavelength campaign led by the Whole Earth Blazar Telescope (WEBT) to study both the quasar core and the beamed jet of this source. The WEBT data are complemented by ultraviolet and X-ray data from Swift, and by gamma-ray data by Fermi. The big blue bump is modelled by using optical and near-infrared mean spectra obtained during the campaign, together with optical and ultraviolet quasar templates. We give prescriptions to correct the source photometry in the various bands for the thermal contribution, in order to derive the non-thermal jet flux. The role of the intergalactic medium absorption is analysed in both the ultraviolet and X-ray bands. We provide opacity values to deabsorb ultraviolet data, and derive a best-guess value for the hydrogen column density through the analysis of X-ray spectra. We estimate the disc and jet bolometric luminosities, accretion rate, and black hole mass. Light curves do not show persistent correlations among flux changes at different frequencies. We study the polarimetric behaviour and find no correlation between polarisation degree and flux, even when correcting for the dilution effect of the big blue bump. Similarly, wide rotations of the electric vector polarisation angle do not seem to be connected with the source activity. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: 15 pages, 14 figures, in press for MNRAS

Target of Opportunity Observations of Gravitational Wave Events with LSST

Authors: R. Margutti, P. Cowperthwaite, Z. Doctor, K. Mortensen, C. P. Pankow, O. Salafia, V. A. Villar, K. Alexander, J. Annis, I. Andreoni, A. Baldeschi, B. Balmaverde, E. Berger, M. G. Bernardini, C. P. L. Berry, F. Bianco, P. K. Blanchard, E. Brocato, M. I. Carnerero, R. Cartier, S. B. Cenko, R. Chornock, L. Chomiuk, C. M. Copperwheat, M. W. Coughlin , et al. (57 additional authors not shown)

Abstract: The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s… ▽ More The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars ($\sim$10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: White paper for LSST cadence optimization- ToOs

Blazars and Fast Radio Bursts with LSST

Authors: C. M. Raiteri, M. I. Carnerero, B. Balmaverde, F. D'Ammando, C. Righi, A. Possenti, E. Pian, A. Capetti, M. Villata, M. Giroletti, P. Romano, A. Stamerra, F. Tavecchio, S. Vercellone

Abstract: The aim of this white paper is to discuss the observing strategies for the LSST Wide-Fast-Deep that would improve the study of blazars (emission variability, census, environment) and Fast Radio Bursts (FRBs). For blazars, these include the adoption of: i) a reference filter to allow reconstruction of a well-sampled light curve not affected by colour changes effects; ii) two snapshots/visit with di… ▽ More The aim of this white paper is to discuss the observing strategies for the LSST Wide-Fast-Deep that would improve the study of blazars (emission variability, census, environment) and Fast Radio Bursts (FRBs). For blazars, these include the adoption of: i) a reference filter to allow reconstruction of a well-sampled light curve not affected by colour changes effects; ii) two snapshots/visit with different exposure times to avoid saturation during flaring states; iii) a rolling cadence to get better-sampled light curves at least in some time intervals. We also address the potential importance of Target of Opportunity (ToO) observations of blazar neutrino sources, and the advantages of a Minisurvey with a star trail cadence (see white paper by David Thomas et al.) for both the blazar science and the detection of possible very fast optical counterparts of FRBs. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: Contribution in response of the Call for White Papers on LSST Cadence Optimization

AGILE, Fermi, Swift, and GASP-WEBT multi-wavelength observations of the high-redshift blazar 4C $+$71.07 in outburst

Authors: S. Vercellone, P. Romano, G. Piano, V. Vittorini, I. Donnarumma, P. Munar-Adrover, C. M. Raiteri, M. Villata, F. Verrecchia, F. Lucarelli, C. Pittori, A. Bulgarelli, V. Fioretti, M. Tavani, J. A. Acosta-Pulido, I. Agudo, A. A. Arkharov, U. Bach, R. Bachev, G. A. Borman, M. S. Butuzova, M. I. Carnerero, C. Casadio, G. Damljanovic, F. D'Ammando , et al. (34 additional authors not shown)

Abstract: The flat-spectrum radio quasar 4C $+$71.07 is a high-redshift ($z=2.172$), $γ$-loud blazar whose optical emission is dominated by the thermal radiation from accretion disc. 4C $+$71.07 has been detected in outburst twice by the AGILE $γ$-ray satellite during the period end of October - mid November 2015, when it reached a $γ$-ray flux of the order of… ▽ More The flat-spectrum radio quasar 4C $+$71.07 is a high-redshift ($z=2.172$), $γ$-loud blazar whose optical emission is dominated by the thermal radiation from accretion disc. 4C $+$71.07 has been detected in outburst twice by the AGILE $γ$-ray satellite during the period end of October - mid November 2015, when it reached a $γ$-ray flux of the order of $F_{\rm E>100\,MeV} = (1.2 \pm 0.3)\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ and $F_{\rm E>100\,MeV} = (3.1 \pm 0.6)\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, respectively, allowing us to investigate the properties of the jet and of the emission region. We investigated its spectral energy distribution by means of almost simultaneous observations covering the cm, mm, near-infrared, optical, ultra-violet, X-ray and $γ$-ray energy bands obtained by the GASP-WEBT Consortium, the Swift and the AGILE and Fermi satellites. The spectral energy distribution of the second $γ$-ray flare (the one whose energy coverage is more dense) can be modelled by means of a one-zone leptonic model, yielding a total jet power of about $4\times10^{47}$ erg s$^{-1}$. During the most prominent $γ$-ray flaring period our model is consistent with a dissipation region within the broad-line region. Moreover, this class of high-redshift, large-mass black-hole flat-spectrum radio quasars might be good targets for future $γ$-ray satellites such as e-ASTROGAM. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: Accepted for publication in A&A. 9 pages, 4 Figures, 3 Tables

Journal ref: A&A 621, A82 (2019)

The extreme HBL behaviour of Markarian 501 during 2012

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin, E. Colombo, J. L. Contreras, J. Cortina, S. Covino , et al. (254 additional authors not shown)

Abstract: A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0… ▽ More A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays. △ Less

Submitted 14 August, 2018; v1 submitted 13 August, 2018; originally announced August 2018.

Comments: 25 pages, 9 figures, accepted for publication in A&A. Corresponding authors: Gareth Hughes (gareth.hughes@cfa.harvard.edu), David Paneque (dpaneque@mppmu.mpg.de), Amit Shukla (amit.shukla@astro.uni-wuerzburg.de)

Journal ref: A&A 620, A181 (2018)

The Bright $γ$-ray Flare of 3C 279 in June 2015: AGILE Detection and Multifrequency Follow-up Observations

Authors: C. Pittori, F. Lucarelli, F. Verrecchia, C. M. Raiteri, M. Villata, V. Vittorini, M. Tavani, S. Puccetti, M. Perri, I. Donnarumma, S. Vercellone, J. A. Acosta-Pulido, R. Bachev, E. Benitez, G. A. Borman, M. I. Carnerero, D. Carosati, W. P. Chen, Sh. A. Ehgamberdiev, A. Goded, T. S. Grishina, D. Hiriart, H. Y. Hsiao, S. G. Jorstad, G. N. Kimeridze , et al. (22 additional authors not shown)

Abstract: We report the AGILE detection and the results of the multifrequency follow-up observations of a bright $γ$-ray flare of the blazar 3C 279 in June 2015. We use AGILE-GRID and Fermi-LAT $γ$-ray data, together with Swift-XRT, Swift-UVOT, and ground-based GASP-WEBT optical observations, including polarization information, to study the source variability and the overall spectral energy distribution dur… ▽ More We report the AGILE detection and the results of the multifrequency follow-up observations of a bright $γ$-ray flare of the blazar 3C 279 in June 2015. We use AGILE-GRID and Fermi-LAT $γ$-ray data, together with Swift-XRT, Swift-UVOT, and ground-based GASP-WEBT optical observations, including polarization information, to study the source variability and the overall spectral energy distribution during the $γ$-ray flare. The $γ$-ray flaring data, compared with as yet unpublished simultaneous optical data which allow to set constraints on the big blue bump disk luminosity, show very high Compton dominance values of $\sim 100$, with a ratio of $γ$-ray to optical emission rising by a factor of three in a few hours. The multi-wavelength behavior of the source during the flare challenges one-zone leptonic theoretical models. The new observations during the June 2015 flare are also compared with already published data and non-simultaneous historical 3C 279 archival data. △ Less

Submitted 20 March, 2018; originally announced March 2018.

Comments: 13 pages, 5 figures,5 Tables. Accepted for publication in ApJ

The detection of the blazar S4 0954+65 at very-high-energy with the MAGIC telescopes during an exceptionally high optical state

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak, P. Colin , et al. (147 additional authors not shown)

Abstract: The very-high-energy (VHE, $\gtrsim 100$ GeV) $γ$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physi… ▽ More The very-high-energy (VHE, $\gtrsim 100$ GeV) $γ$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. We investigate the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and put it in context with possible emission scenarios. We collect photometric data in radio, optical, X-ray, and $γ$ ray. We study both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy $γ$ ray\ (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100$^\circ$. (...) The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars, i.e. inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to flat spectrum radio quasars (FSRQs) and low-frequency peaked BL Lac objects (LBL). △ Less

Submitted 12 January, 2018; originally announced January 2018.

Comments: 15 pages, 7 figures, submitted to A&A

Journal ref: A&A 617, A30 (2018)

An UXor among FUors: extinction related brightness variations of the young eruptive star V582 Aur

Authors: P. Ábrahám, Á. Kóspál, M. Kun, O. Fehér, G. Zsidi, J. A. Acosta-Pulido, M. I. Carnerero, D. Garcia-Álvarez, A. Moór, B. Cseh, G. Hajdu, O. Hanyecz, J. Kelemen, L. Kriskovics, G. Marton, Gy. Mező, L. Molnár, A. Ordasi, G. Rodriguez-Coira, K. Sárneczky, Á. Sódor, R. Szakáts, E. Szegedi-Elek, A. Szing, A. Farkas-Takács , et al. (2 additional authors not shown)

Abstract: V582 Aur is an FU Ori-type young eruptive star in outburst since $\sim$1985. The eruption is currently in a relatively constant plateau phase, with photometric and spectroscopic variability superimposed. Here we will characterize the progenitor of the outbursting object, explore its environment, and analyse the temporal evolution of the eruption. We are particularly interested in the physical orig… ▽ More V582 Aur is an FU Ori-type young eruptive star in outburst since $\sim$1985. The eruption is currently in a relatively constant plateau phase, with photometric and spectroscopic variability superimposed. Here we will characterize the progenitor of the outbursting object, explore its environment, and analyse the temporal evolution of the eruption. We are particularly interested in the physical origin of the two deep photometric dips, one occurred in 2012, and one is ongoing since 2016. We collected archival photographic plates, and carried out new optical, infrared, and millimeter wave photometric and spectroscopic observations between 2010 and 2017, with high sampling rate during the current minimum. Beside analysing the color changes during fading, we compiled multiepoch spectral energy distributions, and fitted them with a simple accretion disk model. Based on pre-outburst data and a millimeter continuum measurement, we suggest that the progenitor of the V582 Aur outburst is a low-mass T Tauri star with average properties. The mass of an unresolved circumstellar structure, probably a disk, is 0.04 M$_{\odot}$. The optical and near-infrared spectra demonstrate the presence of hydrogen and metallic lines, show the CO bandhead in absorption, and exhibit a variable H$α$ profile. The color variations strongly indicate that both the $\sim$1 year long brightness dip in 2012, and the current minimum since 2016 are caused by increased extinction along the line of sight. According to our accretion disk models, the reddening changed from $A_V$=4.5 mag to 12.5 mag, while the accretion rate remained practically constant. Similarly to the models of the UXor phenomenon of intermediate and low-mass young stars, orbiting disk structures could be responsible for the eclipses. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: 20 pages, 11 figures, accepted for publication in ApJ

Dissecting the long-term emission behaviour of the BL Lac object Mrk 421

Authors: M. I. Carnerero, C. M. Raiteri, M. Villata, J. A. Acosta-Pulido, V. M. Larionov, P. S. Smith, F. D'Ammando

Abstract: We report on long-term multi-wavelength monitoring of blazar Mrk 421 by the GASP-WEBT collaboration and Steward Observatory, and by the Swift and Fermi satellites. We study the source behaviour in the period 2007-2015, characterized by several extreme flares. The ratio between the optical, X-ray, and γ-ray fluxes is very variable. The γ-ray flux variations show fair correlation with the optical on… ▽ More We report on long-term multi-wavelength monitoring of blazar Mrk 421 by the GASP-WEBT collaboration and Steward Observatory, and by the Swift and Fermi satellites. We study the source behaviour in the period 2007-2015, characterized by several extreme flares. The ratio between the optical, X-ray, and γ-ray fluxes is very variable. The γ-ray flux variations show fair correlation with the optical ones starting from 2012. We analyse spectropolarimetric data and find wavelength-dependence of the polarisation degree (P), which is compatible with the presence of the host galaxy, and no wavelength-dependence of the electric vector polarisation angle (EVPA). Optical polarimetry shows a lack of simple correlation between P and flux and wide rotations of the EVPA. We build broad-band spectral energy distributions with simultaneous near-infrared and optical data from the GASP-WEBT and UV and X-ray data from the Swift satellite. They show strong variability in both flux and X-ray spectral shape and suggest a shift of the synchrotron peak up to a factor $\sim 50$ in frequency. The interpretation of the flux and spectral variability is compatible with jet models including at least two emitting regions that can change their orientation with respect to the line of sight. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Comments: 18 pages, 27 figures, 1 table, in press for MNRAS

Multiwavelength observations of a VHE gamma-ray flare from PKS 1510-089 in 2015

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, R. Carosi, A. Carosi, A. Chatterjee, P. Colin , et al. (151 additional authors not shown)

Abstract: Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In… ▽ More Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models. △ Less

Submitted 28 February, 2017; v1 submitted 28 October, 2016; originally announced October 2016.

Comments: 12 pages, 10 figures, accepted for a publication in A&A

Journal ref: A&A 603, A29 (2017)

Exceptional outburst of the blazar CTA 102 in 2012: The GASP-WEBT campaign and its extension

Authors: V. M. Larionov, M. Villata, C. M. Raiteri, S. G. Jorstad, A. P. Marscher, I. Agudo, P. S. Smith, J. A. Acosta-Pulido, M. J. Arévalo, A. A. Arkharov, R. Bachev, D. A. Blinov, G. Borisov, G. A. Borman, V. Bozhilov, A. Bueno, M. I. Carnerero, D. Carosati, C. Casadio, W. P. Chen, D. P. Clemens, A. Di Paola, Sh. A. Ehgamberdiev, J. L. Gómez, P. A González-Morales , et al. (27 additional authors not shown)

Abstract: After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished… ▽ More After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: 11 pages, 14 figures, accepted in MNRAS