Showing 1–28 of 28 results for author: Sanchez-Saez, P

SDSS1335+0728: The awakening of a $\sim 10^6 M_{\odot}$ black hole

Authors: P. Sánchez-Sáez, L. Hernández-García, S. Bernal, A. Bayo, G. Calistro Rivera, F. E. Bauer, C. Ricci, A. Merloni, M. J. Graham, R. Cartier, P. Arévalo, R. J. Assef, A. Concas, D. Homan, M. Krumpe, P. Lira, A. Malyali, M. L. Martínez-Aldama, A. M. Muñoz Arancibia, A. Rau, G. Bruni, F. Förster, M. Pavez-Herrera, D. Tubín-Arenas, M. Brightman

Abstract: The galaxy SDSS1335+0728, which had exhibited no prior optical variations during the preceding two decades, began showing significant nuclear variability in the Zwicky Transient Facility (ZTF) alert stream from December 2019 (as ZTF19acnskyy). Its behaviour suggests that SDSS1335+0728 hosts a $\sim 10^6 M_{\odot}$ black hole (BH) that is currently in the process of `turning on'. We present a multi… ▽ More The galaxy SDSS1335+0728, which had exhibited no prior optical variations during the preceding two decades, began showing significant nuclear variability in the Zwicky Transient Facility (ZTF) alert stream from December 2019 (as ZTF19acnskyy). Its behaviour suggests that SDSS1335+0728 hosts a $\sim 10^6 M_{\odot}$ black hole (BH) that is currently in the process of `turning on'. We present a multi-wavelength photometric analysis and spectroscopic follow-up performed with the aim of better understanding the origin of the nuclear variations detected in SDSS1335+0728. We used archival photometry and spectroscopic data to study the state of SDSS1335+0728 prior to December 2019, and new observations from Swift, SOAR/Goodman, VLT/X-shooter, and Keck/LRIS taken after its turn-on to characterise its current state. We find that: (a) since 2021, the UV flux is four times brighter than the flux reported by GALEX in 2004; (b) since June 2022, the mid-infrared flux has risen more than two times, and the W1-W2 WISE colour has become redder; (c) since February 2024, the source has begun showing X-ray emission; (d) the narrow emission line ratios are now consistent with a more energetic ionising continuum; (e) broad emission lines are not detected; and (f) the [OIII] line increased its flux $\sim 3.6$ years after the first ZTF alert, which implies a relatively compact narrow-line-emitting region. We conclude that the variations observed in SDSS1335+0728 could be either explained by an AGN that is just turning on or by an exotic tidal disruption event (TDE). If the former is true, SDSS1335+0728 is one of the strongest cases of an AGN observed in the process of activating. If the latter, it would correspond to the longest and faintest TDE ever observed (or another class of still unknown nuclear transient). Future observations of SDSS1335+0728 are crucial to further understand its behaviour. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. ESO press release available at https://www.eso.org/public/news/eso2409/. Abstract shortened for arXiv

ATAT: Astronomical Transformer for time series And Tabular data

Authors: G. Cabrera-Vives, D. Moreno-Cartagena, N. Astorga, I. Reyes-Jainaga, F. Förster, P. Huijse, J. Arredondo, A. M. Muñoz Arancibia, A. Bayo, M. Catelan, P. A. Estévez, P. Sánchez-Sáez, A. Álvarez, P. Castellanos, P. Gallardo, A. Moya, D. Rodriguez-Mancini

Abstract: The advent of next-generation survey instruments, such as the Vera C. Rubin Observatory and its Legacy Survey of Space and Time (LSST), is opening a window for new research in time-domain astronomy. The Extended LSST Astronomical Time-Series Classification Challenge (ELAsTiCC) was created to test the capacity of brokers to deal with a simulated LSST stream. We describe ATAT, the Astronomical Trans… ▽ More The advent of next-generation survey instruments, such as the Vera C. Rubin Observatory and its Legacy Survey of Space and Time (LSST), is opening a window for new research in time-domain astronomy. The Extended LSST Astronomical Time-Series Classification Challenge (ELAsTiCC) was created to test the capacity of brokers to deal with a simulated LSST stream. We describe ATAT, the Astronomical Transformer for time series And Tabular data, a classification model conceived by the ALeRCE alert broker to classify light-curves from next-generation alert streams. ATAT was tested in production during the first round of the ELAsTiCC campaigns. ATAT consists of two Transformer models that encode light curves and features using novel time modulation and quantile feature tokenizer mechanisms, respectively. ATAT was trained on different combinations of light curves, metadata, and features calculated over the light curves. We compare ATAT against the current ALeRCE classifier, a Balanced Hierarchical Random Forest (BHRF) trained on human-engineered features derived from light curves and metadata. When trained on light curves and metadata, ATAT achieves a macro F1-score of 82.9 +- 0.4 in 20 classes, outperforming the BHRF model trained on 429 features, which achieves a macro F1-score of 79.4 +- 0.1. The use of Transformer multimodal architectures, combining light curves and tabular data, opens new possibilities for classifying alerts from a new generation of large etendue telescopes, such as the Vera C. Rubin Observatory, in real-world brokering scenarios. △ Less

Submitted 16 May, 2024; v1 submitted 5 May, 2024; originally announced May 2024.

Physical Properties of Type II Supernovae Inferred from ZTF and ATLAS Photometric Data

Authors: Javier Silva-Farfán, Francisco Förster, Takashi J. Moriya, L. Hernández-García, A. M. Muñoz Arancibia, P. Sánchez-Sáez, Joseph P. Anderson, John L. Tonry, Alejandro Clocchiatti

Abstract: We report an analysis of a sample of 186 spectroscopically confirmed Type II supernova (SN) light curves (LCs) obtained from a combination of Zwicky Transient Facility (ZTF) and Asteroid Terrestrial-impact Last Alert System (ATLAS) observations. We implement a method to infer physical parameters from these LCs using hydrodynamic models that take into account the progenitor mass, the explosion ener… ▽ More We report an analysis of a sample of 186 spectroscopically confirmed Type II supernova (SN) light curves (LCs) obtained from a combination of Zwicky Transient Facility (ZTF) and Asteroid Terrestrial-impact Last Alert System (ATLAS) observations. We implement a method to infer physical parameters from these LCs using hydrodynamic models that take into account the progenitor mass, the explosion energy, and the presence of circumstellar matter (CSM). The CSM is modelled via the mass loss rate, wind acceleration at the surface of the progenitor star with a $β$ velocity law, and the CSM radius. We also infer the time of explosion, attenuation (A$_V$), and the redshift for each SN. Our results favor low-mass progenitor stars (M$_{ZAMS}$\,$<$14\,$M_\odot$) with a dense CSM ($\dot{M}$ $>$ 10$^{-3}$ [M$_\odot$ yr$^{-1}$], a CSM radius of $\sim$ 10$^{15}$ cm, and $β$ $>$ 2). Additionally, we find that the redshift inferred from the supernova LCs is significantly more accurate than that inferred using the host galaxy photometric redshift, suggesting that this method could be used to infer more accurate host galaxy redshifts from large samples of SNe II in the LSST era. Lastly, we compare our results with similar works from the literature. △ Less

Submitted 3 June, 2024; v1 submitted 19 April, 2024; originally announced April 2024.

The Wide-field Spectroscopic Telescope (WST) Science White Paper

Authors: Vincenzo Mainieri, Richard I. Anderson, Jarle Brinchmann, Andrea Cimatti, Richard S. Ellis, Vanessa Hill, Jean-Paul Kneib, Anna F. McLeod, Cyrielle Opitom, Martin M. Roth, Paula Sanchez-Saez, Rodolfo Smiljanic, Eline Tolstoy, Roland Bacon, Sofia Randich, Angela Adamo, Francesca Annibali, Patricia Arevalo, Marc Audard, Stefania Barsanti, Giuseppina Battaglia, Amelia M. Bayo Aran, Francesco Belfiore, Michele Bellazzini, Emilio Bellini , et al. (192 additional authors not shown)

Abstract: The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integ… ▽ More The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate △ Less

Submitted 12 April, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

A Newborn AGN in a Starforming Galaxy

Authors: P. Arévalo, E. López-Navas, M. L. Martínez-Aldama, P. Lira, S. Bernal, P. Sánchez-Sáez, M. Salvato, L. Hernández-García, C. Ricci, A. Merloni, M. Krumpe

Abstract: We report on the finding of a newborn AGN, i.e. current AGN activity in a galaxy previously classified as non-active, and characterize its evolution. Black hole ignition event candidates were selected from a parent sample of spectrally classified non-active galaxies (2.394.312 objects), that currently show optical flux variability indicative of a type I AGN, according to the ALeRCE light curve cla… ▽ More We report on the finding of a newborn AGN, i.e. current AGN activity in a galaxy previously classified as non-active, and characterize its evolution. Black hole ignition event candidates were selected from a parent sample of spectrally classified non-active galaxies (2.394.312 objects), that currently show optical flux variability indicative of a type I AGN, according to the ALeRCE light curve classifier. A second epoch spectrum for a sample of candidate newborn AGN were obtained with the SOAR telescope to search for new AGN features. We present spectral results for the most convincing case of new AGN activity, for a galaxy with a previous star-forming optical classification, where the second epoch spectrum shows the appearance of prominent, broad Balmer lines without significant changes in the narrow line flux ratios. Long term optical lightcurves show a steady increase in luminosity starting 1.5 years after the SDSS spectrum was taken and continuing for at least 7 years. MIR colors from the WISE catalog have also evolved from typical non active galaxy colors to AGN-like colors and recent X-ray flux detections confirm its AGN nature. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A Letters

Journal ref: A&A 683, L8 (2024)

Selection of powerful radio galaxies with machine learning

Authors: R. Carvajal, I. Matute, J. Afonso, R. P. Norris, K. J. Luken, P. Sánchez-Sáez, P. A. C. Cunha, A. Humphrey, H. Messias, S. Amarantidis, D. Barbosa, H. A. Cruz, H. Miranda, A. Paulino-Afonso, C. Pappalardo

Abstract: We developed and trained a pipeline of three machine learning (ML) models than can predict which sources are more likely to be an AGN and to be detected in specific radio surveys. Also, it can estimate redshift values for predicted radio-detectable AGNs. These models, which combine predictions from tree-based and gradient-boosting algorithms, have been trained with multi-wavelength data from near-… ▽ More We developed and trained a pipeline of three machine learning (ML) models than can predict which sources are more likely to be an AGN and to be detected in specific radio surveys. Also, it can estimate redshift values for predicted radio-detectable AGNs. These models, which combine predictions from tree-based and gradient-boosting algorithms, have been trained with multi-wavelength data from near-infrared-selected sources in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring field. Training, testing, calibration, and validation were carried out in the HETDEX field. Further validation was performed on near-infrared-selected sources in the Stripe 82 field. In the HETDEX validation subset, our pipeline recovers 96% of the initially labelled AGNs and, from AGNs candidates, we recover 50% of previously detected radio sources. For Stripe 82, these numbers are 94% and 55%. Compared to random selection, these rates are two and four times better for HETDEX, and 1.2 and 12 times better for Stripe 82. The pipeline can also recover the redshift distribution of these sources with $σ_{\mathrm{NMAD}}$ = 0.07 for HETDEX ($σ_{\mathrm{NMAD}}$ = 0.09 for Stripe 82) and an outlier fraction of 19% (25% for Stripe 82), compatible with previous results based on broad-band photometry. Feature importance analysis stresses the relevance of near- and mid-infrared colours to select AGNs and identify their radio and redshift nature. Combining different algorithms in ML models shows an improvement in the prediction power of our pipeline over a random selection of sources. Tree-based ML models (in contrast to deep learning techniques) facilitate the analysis of the impact that features have on the predictions. This prediction can give insight into the potential physical interplay between the properties of radio AGNs (e.g. mass of black hole and accretion rate). △ Less

Submitted 1 December, 2023; v1 submitted 20 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A, 24 pages and 21 figures. Updated to include information on retrieval of data and models, which can be obtained from https://zenodo.org/doi/10.5281/zenodo.10220008

Journal ref: A&A 679, A101 (2023)

Multi-Class Deep SVDD: Anomaly Detection Approach in Astronomy with Distinct Inlier Categories

Authors: Manuel Pérez-Carrasco, Guillermo Cabrera-Vives, Lorena Hernández-García, Francisco Forster, Paula Sánchez-Sáez, Alejandra Muñoz Arancibia, Nicolás Astorga, Franz Bauer, Amelia Bayo, Martina Cádiz-Leyton, Marcio Catelan

Abstract: With the increasing volume of astronomical data generated by modern survey telescopes, automated pipelines and machine learning techniques have become crucial for analyzing and extracting knowledge from these datasets. Anomaly detection, i.e. the task of identifying irregular or unexpected patterns in the data, is a complex challenge in astronomy. In this paper, we propose Multi-Class Deep Support… ▽ More With the increasing volume of astronomical data generated by modern survey telescopes, automated pipelines and machine learning techniques have become crucial for analyzing and extracting knowledge from these datasets. Anomaly detection, i.e. the task of identifying irregular or unexpected patterns in the data, is a complex challenge in astronomy. In this paper, we propose Multi-Class Deep Support Vector Data Description (MCDSVDD), an extension of the state-of-the-art anomaly detection algorithm One-Class Deep SVDD, specifically designed to handle different inlier categories with distinct data distributions. MCDSVDD uses a neural network to map the data into hyperspheres, where each hypersphere represents a specific inlier category. The distance of each sample from the centers of these hyperspheres determines the anomaly score. We evaluate the effectiveness of MCDSVDD by comparing its performance with several anomaly detection algorithms on a large dataset of astronomical light-curves obtained from the Zwicky Transient Facility. Our results demonstrate the efficacy of MCDSVDD in detecting anomalous sources while leveraging the presence of different inlier categories. The code and the data needed to reproduce our results are publicly available at https://github.com/mperezcarrasco/AnomalyALeRCE. △ Less

Submitted 10 August, 2023; v1 submitted 9 August, 2023; originally announced August 2023.

Comments: Accepted to ICML 2023 Workshop on Machine Learning for Astrophysics

Improving the selection of changing-look AGNs through multi-wavelength photometric variability

Authors: E. López-Navas, P. Sánchez-Sáez, P. Arévalo, S. Bernal, M. J. Graham, L. Hernández-García, D. Homan, M. Krumpe, G. Lamer, P. Lira, M. L. Martínez-Aldama, A. Merloni, S. Ríos, M. Salvato, D. Stern, D. Tubín-Arenas

Abstract: We present second epoch optical spectra for 30 changing-look (CL) candidates found by searching for Type-1 optical variability in a sample of active galactic nuclei (AGNs) spectroscopically classified as Type 2. We use a random-forest-based light curve classifier and spectroscopic follow-up, confirming 50 per cent of candidates as turning-on CLs. In order to improve this selection method and to be… ▽ More We present second epoch optical spectra for 30 changing-look (CL) candidates found by searching for Type-1 optical variability in a sample of active galactic nuclei (AGNs) spectroscopically classified as Type 2. We use a random-forest-based light curve classifier and spectroscopic follow-up, confirming 50 per cent of candidates as turning-on CLs. In order to improve this selection method and to better understand the nature of the not-confirmed CL candidates, we perform a multi-wavelength variability analysis including optical, mid-infrared (MIR) and X-ray data, and compare the results from the confirmed and not-confirmed CLs identified in this work. We find that most of the not-confirmed CLs are consistent with weak Type 1s dominated by host-galaxy contributions, showing weaker optical and MIR variability. On the contrary, the confirmed CLs present stronger optical fluctuations and experience a long (from five to ten years) increase in their MIR fluxes and the colour W1-W2 over time. In the 0.2-2.3 keV band, at least four out of 11 CLs with available SRG/eROSITA detections have increased their flux in comparison with archival upper limits. These common features allow us to select the most promising CLs from our list of candidates, leading to nine sources with similar multi-wavelength photometric properties to our CL sample. The use of machine learning algorithms with optical and MIR light curves will be very useful to identify CLs in future large-scale surveys. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Journal ref: Monthly Notices of the Royal Astronomical Society, 2023; stad1893

The universal power spectrum of Quasars in optical wavelengths: Break timescale scales directly with both black hole mass and accretion rate

Authors: P. Arevalo, E. Churazov, P. Lira, P. Sanchez-Saez, S. Bernal, L. Hernandez-Garcia, E. Lopez-Navas, P. Patel

Abstract: Aims: Establish the dependence of variability properties, such as characteristic timescales and variability amplitude, on basic quasar parameters such as black hole mass and accretion rate, controlling for the rest-frame wavelength of emission. Methods: Using large catalogs of quasars, we selected the g-band light curves for 4770 objects from the Zwicky Transient Facility archive. All selected obj… ▽ More Aims: Establish the dependence of variability properties, such as characteristic timescales and variability amplitude, on basic quasar parameters such as black hole mass and accretion rate, controlling for the rest-frame wavelength of emission. Methods: Using large catalogs of quasars, we selected the g-band light curves for 4770 objects from the Zwicky Transient Facility archive. All selected objects fall into a narrow redshift bin, $0.6<z<0.7$, but cover a wide range of accretion rates in Eddington units (REdd) and black hole masses ($M$). We grouped these objects into 26 independent bins according to these parameters, calculated low-resolution $g$-band variability power spectra for each of these bins, and approximated the power spectra with a simple analytic model that features a break at a timescale $t_b$. Results: We found a clear dependence of the break timescale $t_b$ on REdd, on top of the known dependence of $t_b$ on the black hole mass $M$. In our fits, $t_b\propto M^{0.65 - 0.55}$ REdd $^{0.35 - 0.3}$, where the ranges in the exponents correspond to the best-fitting parameters of different power spectrum models. Scaling $t_b$ to the orbital timescale of the innermost stable circular orbit (ISCO), $t_{\rm ISCO}$, results approximately in $t_{b}/t_{\rm ISCO} \propto ($REdd$/M)^{0.35}$. The observed values of $t_b$ are $\sim 10$ longer than the orbital timescale at the light-weighted average radius of the disc region emitting in the (observer frame) $g$-band. The different scaling of the break frequency with $M$ and REdd shows that the shape of the variability power spectrum cannot be solely a function of the quasar luminosity, even for a single rest-frame wavelength. Finally, the best-fitting models have slopes above the break in the range -2.5 and -3. A slope of -2, as in the damped random walk models, fits the data significantly worse. △ Less

Submitted 19 November, 2023; v1 submitted 19 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A

Optical variability in Quasars: Scaling with black hole mass and Eddington ratio depend on the observed timescales

Authors: Patricia Arévalo, Paulina Lira, Paula Sánchez-Sáez, Priyanjali Patel, Elena López-Navas, Eugene Churazov, Lorena Hernández-García

Abstract: Quasars emission is highly variable, and this variability gives us clues to understand the accretion process onto supermassive black holes. We can expect variability properties to correlate with the main physical properties of the accreting black hole, i.e., its mass and accretion rate. It has been established that the relative amplitude of variability anti-correlates with the accretion rate.The d… ▽ More Quasars emission is highly variable, and this variability gives us clues to understand the accretion process onto supermassive black holes. We can expect variability properties to correlate with the main physical properties of the accreting black hole, i.e., its mass and accretion rate. It has been established that the relative amplitude of variability anti-correlates with the accretion rate.The dependence of the variance on black hole mass has remained elusive, and contradicting results, including positive, negative, or no correlation, have been reported. In this work, we show that the key to these contradictions lies in the timescales of variability studied (e.g., the length of the light curves available). By isolating the variance on different timescales as well as mass and accretion rate bins we show that there is indeed a negative correlation between black hole mass and variance and that this anti-correlation is stronger for shorter timescale fluctuations. The behavior can be explained in terms of a universal variability power spectrum for all quasars, resembling a broken power law where the variance is constant at low temporal frequencies and then drops continuously for frequencies higher than a characteristic frequency $f_b$, where $f_b$ correlates with the black hole mass. Furthermore, to explain all the variance results presented here, not only the normalization of this power spectrum must anti-correlate with the accretion rate, but also the shape of the power spectra at short timescales must depend on this parameter as well. △ Less

Submitted 19 October, 2023; v1 submitted 27 April, 2023; originally announced April 2023.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Persistent and occasional: searching for the variable population of the ZTF/4MOST sky using ZTF data release 11

Authors: P. Sánchez-Sáez, J. Arredondo, A. Bayo, P. Arévalo, F. E. Bauer, G. Cabrera-Vives, M. Catelan, P. Coppi, P. A. Estévez, F. Förster, L. Hernández-García, P. Huijse, R. Kurtev, P. Lira, A. M. Muñoz Arancibia, G. Pignata

Abstract: We present a variability, color and morphology based classifier, designed to identify transients, persistently variable, and non-variable sources, from the Zwicky Transient Facility (ZTF) Data Release 11 (DR11) light curves of extended and point sources. The main motivation to develop this model was to identify active galactic nuclei (AGN) at different redshift ranges to be observed by the 4MOST C… ▽ More We present a variability, color and morphology based classifier, designed to identify transients, persistently variable, and non-variable sources, from the Zwicky Transient Facility (ZTF) Data Release 11 (DR11) light curves of extended and point sources. The main motivation to develop this model was to identify active galactic nuclei (AGN) at different redshift ranges to be observed by the 4MOST ChANGES project. Still, it serves as a more general time-domain astronomy study. The model uses nine colors computed from CatWISE and PS1, a morphology score from PS1, and 61 single-band variability features computed from the ZTF DR11 g and r light curves. We trained two versions of the model, one for each ZTF band. We used a hierarchical local classifier per parent node approach, where each node was composed of a balanced random forest model. We adopted a 17-class taxonomy, including non-variable stars and galaxies, three transient classes, five classes of stochastic variables, and seven classes of periodic variables. The macro averaged precision, recall and F1-score are 0.61, 0.75, and 0.62 for the g-band model, and 0.60, 0.74, and 0.61, for the r-band model. When grouping the four AGN classes into one single class, its precision, recall, and F1-score are 1.00, 0.95, and 0.97, respectively, for both the g and r bands. We applied the model to all the sources in the ZTF/4MOST overlapping sky, avoiding ZTF fields covering the Galactic bulge, including 86,576,577 light curves in the g-band and 140,409,824 in the r-band. Only 0.73\% of the g-band light curves and 2.62\% of the r-band light curves were classified as stochastic, periodic, or transient with high probability ($P_{init}\geq0.9$). We found that, in general, more reliable results are obtained when using the g-band model. Using the latter, we identified 384,242 AGN candidates, 287,156 of which have $P_{init}\geq0.9$. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. Abstract shortened for arXiv. Tables containing the classifications and features for the ZTF g and r bands, and the labeled sets will be available at CDS. Individual catalogs per class and band, as well as the labeled set catalogs, can be downloaded at Zenodo DOI:10.5281/zenodo.7826045

Journal ref: A&A 675, A195 (2023)

Multiwavelength monitoring of the nucleus in PBC J2333.9-2343: the giant radio galaxy with a blazar-like core

Authors: L. Hernández-García, F. Panessa, G. Bruni, L. Bassani, P. Arévalo, V. M. Patiño-Alvarez, A. Tramacere, P. Lira, P. Sánchez-Sáez, F. E. Bauer, V. Chavushyan, R. Carraro, F. Förster, A. M. Muñoz Arancibia, P. Ubertini

Abstract: PBC J2333.9-2343 is a giant radio galaxy at z = 0.047 with a bright central core associated to a blazar nucleus. If the nuclear blazar jet is a new phase of the jet activity, then the small orientation angle suggest a dramatic change of the jet direction. We present observations obtained between September 2018 and January 2019 (cadence larger than three days) with Effeslberg, SMARTS-1.3m, ZTF, ATL… ▽ More PBC J2333.9-2343 is a giant radio galaxy at z = 0.047 with a bright central core associated to a blazar nucleus. If the nuclear blazar jet is a new phase of the jet activity, then the small orientation angle suggest a dramatic change of the jet direction. We present observations obtained between September 2018 and January 2019 (cadence larger than three days) with Effeslberg, SMARTS-1.3m, ZTF, ATLAS, Swift, and Fermi-LAT, and between April-July 2019 (daily cadence) with SMARTS-1.3m and ATLAS. Large (>2x) flux increases are observed on timescales shorter than a month, which are interpreted as flaring events. The cross correlation between the SMARTS-1.3m monitoring in the NIR and optical shows that these data do not show significant time lag within the measured errors. A comparison of the optical variability properties between non-blazars and blazars AGN shows that PBC J2333.9-2343 has properties more comparable to the latter. The SED of the nucleus shows two peaks, that were fitted with a one zone leptonic model. Our data and modelling shows that the high energy peak is dominated by External Compton from the dusty torus with mild contribution from Inverse Compton from the jet. The derived jet angle of 3 degrees is also typical of a blazar. Therefore, we confirm the presence of a blazar-like core in the center of this giant radio galaxy, likely a Flat Spectrum Radio Quasar with peculiar properties. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS. 15 pages, 6 figures, 2 appendix including tables and figures

The Type 1 and Type 2 AGN dichotomy according to their ZTF optical variability

Authors: E. López-Navas, P. Arévalo, S. Bernal, Matthew J. Graham, L. Hernández-García, P. Lira, P. Sánchez-Sáez

Abstract: The scarce optical variability studies in spectrally classified Type 2 active galactic nuclei (AGNs) have led to the discovery of anomalous objects that are incompatible with the simplest unified models (UM). This paper focuses on the exploration of different variability features that allows to separate between obscured, Type 2 AGNs, and the variable, unobscured Type 1s. We analyse systematically… ▽ More The scarce optical variability studies in spectrally classified Type 2 active galactic nuclei (AGNs) have led to the discovery of anomalous objects that are incompatible with the simplest unified models (UM). This paper focuses on the exploration of different variability features that allows to separate between obscured, Type 2 AGNs, and the variable, unobscured Type 1s. We analyse systematically the Zwicky Transient Facility, 2.5 years long light curves of ~ 15000 AGNs from the Sloan Digital Sky Survey Data Release 16, which are generally considered Type 2s due to the absence of strong broad emission lines (BELs). Consistently with the expectations from the UM, the variability features are distributed differently for distinct populations, with spectrally classified weak Type 1s showing 1 order of magnitude larger variances than the Type 2s. We find that the parameters given by the damped random walk model leads to broader Hα equivalent width for objects with τ_g > 16 d and long term structure function SF{\infty},g> 0.07 mag. By limiting the variability features, we find that ~ 11 per cent of Type 2 sources show evidence for optical variations. A detailed spectral analysis of the most variable sources (~1 per cent of the Type 2 sample) leads to the discovery of misclassified Type 1s with weak BELs and changing-state candidates. This work presents one of the largest systematic investigations of Type 2 AGN optical variability to date, in preparation for future large photometric surveys. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS. 13 pages, 9 figures

The LSST era of supermassive black holes accretion-disk reverberation mapping

Authors: Andjelka B. Kovacevic, Viktor Radovic, Dragana Ilic, Luka C. Popovic, Roberto J. Assef, Paula Sanchez-Saez, Robert Nikutta, Claudia M. Raiteri, Ilsang Yoon, Yasaman Homayouni, Yan-Rong Li, Neven Caplar, Bozena Czerny, Swayamtrupta Panda, Claudio Ricci, Isidora Jankov, Hermine Landt, Christian Wolf, Jelena Kovacevic-Dojcinovic, Masa Lakicevic, Djorđe Savic, Oliver Vince, Sasa Simic, Iva Cvorovic-Hajdinjak, Sladjana Marceta-Mandic

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

Submitted 12 August, 2022; originally announced August 2022.

Comments: accepted for publication in the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

DELIGHT: Deep Learning Identification of Galaxy Hosts of Transients using Multi-resolution Images

Authors: Francisco Förster, Alejandra M. Muñoz Arancibia, Ignacio Reyes, Alexander Gagliano, Dylan Britt, Sara Cuellar-Carrillo, Felipe Figueroa-Tapia, Ava Polzin, Yara Yousef, Javier Arredondo, Diego Rodríguez-Mancini, Javier Correa-Orellana, Amelia Bayo, Franz E. Bauer, Márcio Catelan, Guillermo Cabrera-Vives, Raya Dastidar, Pablo A. Estévez, Giuliano Pignata, Lorena Hernandez-Garcia, Pablo Huijse, Esteban Reyes, Paula Sánchez-Sáez, Mauricio Ramirez, Daniela Grandón , et al. (3 additional authors not shown)

Abstract: We present DELIGHT, or Deep Learning Identification of Galaxy Hosts of Transients, a new algorithm designed to automatically and in real-time identify the host galaxies of extragalactic transients. The proposed algorithm receives as input compact, multi-resolution images centered at the position of a transient candidate and outputs two-dimensional offset vectors that connect the transient with the… ▽ More We present DELIGHT, or Deep Learning Identification of Galaxy Hosts of Transients, a new algorithm designed to automatically and in real-time identify the host galaxies of extragalactic transients. The proposed algorithm receives as input compact, multi-resolution images centered at the position of a transient candidate and outputs two-dimensional offset vectors that connect the transient with the center of its predicted host. The multi-resolution input consists of a set of images with the same number of pixels, but with progressively larger pixel sizes and fields of view. A sample of \nSample galaxies visually identified by the ALeRCE broker team was used to train a convolutional neural network regression model. We show that this method is able to correctly identify both relatively large ($10\arcsec < r < 60\arcsec$) and small ($r \le 10\arcsec$) apparent size host galaxies using much less information (32 kB) than with a large, single-resolution image (920 kB). The proposed method has fewer catastrophic errors in recovering the position and is more complete and has less contamination ($< 0.86\%$) recovering the cross-matched redshift than other state-of-the-art methods. The more efficient representation provided by multi-resolution input images could allow for the identification of transient host galaxies in real-time, if adopted in alert streams from new generation of large etendue telescopes such as the Vera C. Rubin Observatory. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: Submitted to The Astronomical Journal on Aug 5th, 2022. Comments and suggestions are welcome

From Data to Software to Science with the Rubin Observatory LSST

Authors: Katelyn Breivik, Andrew J. Connolly, K. E. Saavik Ford, Mario Jurić, Rachel Mandelbaum, Adam A. Miller, Dara Norman, Knut Olsen, William O'Mullane, Adrian Price-Whelan, Timothy Sacco, J. L. Sokoloski, Ashley Villar, Viviana Acquaviva, Tomas Ahumada, Yusra AlSayyad, Catarina S. Alves, Igor Andreoni, Timo Anguita, Henry J. Best, Federica B. Bianco, Rosaria Bonito, Andrew Bradshaw, Colin J. Burke, Andresa Rodrigues de Campos , et al. (75 additional authors not shown)

Abstract: The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the po… ▽ More The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: White paper from "From Data to Software to Science with the Rubin Observatory LSST" workshop

A structure function analysis of VST-COSMOS AGN

Authors: D. De Cicco, F. E. Bauer, M. Paolillo, P. Sánchez-Sáez, W. N. Brandt, F. Vagnetti, G. Pignata, M. Radovich, M. Vaccari

Abstract: We present our sixth work in a series dedicated to variability studies of active galactic nuclei (AGN) based on the survey of the COSMOS field by the VLT Survey Telescope (VST). Its 54 r-band visits over 3.3 yr and single-visit depth of 24.6 r-band mag make this dataset a valuable scaled-down version that can help forecast the performance of the Rubin Observatory Legacy Survey of Space and Time (L… ▽ More We present our sixth work in a series dedicated to variability studies of active galactic nuclei (AGN) based on the survey of the COSMOS field by the VLT Survey Telescope (VST). Its 54 r-band visits over 3.3 yr and single-visit depth of 24.6 r-band mag make this dataset a valuable scaled-down version that can help forecast the performance of the Rubin Observatory Legacy Survey of Space and Time (LSST). This work is centered on the analysis of the structure function (SF) of VST-COSMOS AGN, investigating possible differences in its shape and slope related to how the AGN were selected, and explores possible connections between the ensemble variability of AGN and black-hole mass, accretion rate, bolometric luminosity, redshift, and obscuration of the source. Given its features, our dataset opens up the exploration of samples ~2 mag fainter than most of the literature to date. We identify several samples of AGN - 677 in total - obtained by a variety of selection techniques which partly overlap. Our analysis compares results for the various samples. We split each sample in two based on the median of the physical property of interest, and analyze differences in the shape and slope of the SF, and possible causes. While the shape of the SF does not change with depth, it is highly affected by the type of AGN (unobscured/obscured) included in the sample. Where a linear region can be identified, we find that the variability amplitude anticorrelates with accretion rate and bolometric luminosity, consistent with previous literature on the topic, while no dependence on black-hole mass emerges from this study. With its longer baseline and denser and more regular sampling, the LSST will allow an improved characterization of the SF and its dependencies on the mentioned physical properties over much larger AGN samples. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 18 pages, 10 figures; accepted for publication in A&A

Two years of optical and NIR observations of the superluminous supernova UID 30901 discovered by the UltraVISTA SN survey

Authors: Emilio D. Hueichapán, Carlos Contreras V., Regis Cartier, Paulina Lira, Paula Sanchez-Saez, Bo Milvang-Jensen, Johan P. U. Fynbo, Joseph P. Anderson, Mario Hamuy

Abstract: We present deep optical and near-infrared photometry of UID 30901, a superluminous supernova (SLSN) discovered during the UltraVISTA survey. The observations were obtained with VIRCAM ($YJHK_{s}$) mounted on the VISTA telescope, DECam ($griz$) on the Blanco telescope, and SUBARU Hyper Suprime-Cam (HSC; $grizy$). These multi-band observations comprise +700 days making UID 30901 one of the best phot… ▽ More We present deep optical and near-infrared photometry of UID 30901, a superluminous supernova (SLSN) discovered during the UltraVISTA survey. The observations were obtained with VIRCAM ($YJHK_{s}$) mounted on the VISTA telescope, DECam ($griz$) on the Blanco telescope, and SUBARU Hyper Suprime-Cam (HSC; $grizy$). These multi-band observations comprise +700 days making UID 30901 one of the best photometrically followed SLSNe to date. The host galaxy of UID 30901 is detected in a deep HST F814W image with an AB magnitude of $27.3 \pm 0.2$. While no spectra exist for the SN or its host galaxy, we perform our analysis assuming $z = 0.37$, based on the photometric redshift of a possible host galaxy found at a projected distance of 7 kpc. Fitting a blackbody to the observations, the radius, temperature, and bolometric light curve are computed. We find a maximum bolometric luminosity of $5.4 \pm 0.34 \times 10^{43}$ erg s$^{-1}$. A flattening in the light curve beyond 600 days is observed and several possible causes are discussed. We find the observations to clearly favour a SLSN type I, and plausible power sources such as the radioactive decay of $^{56}$Ni and the spin-down of a magnetar are compared to the data. We find that the magnetar model yields a good fit to the observations with the following parameters: a magnetic field $B = 1.4 \pm 0.3 \times 10^{14} \ G$, spin period of $P = 6.0 \pm 0.1 \ ms$ and ejecta mass $M_{ej} = 11.9^{+4.8}_{-6.4} M_{\odot}$. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 22 pages, 15 figures, Accepted for publication in MNRAs

Confirming new changing--look AGNs discovered through optical variability using a random-forest based light curve classifier

Authors: E. López-Navas, M. L. Martínez-Aldama, S. Bernal, P. Sánchez-Sáez, P. Arévalo, Matthew J. Graham, L. Hernández-García, P. Lira, P. A. Rojas Lobos

Abstract: Determining the frequency and duration of changing--look (CL) active galactic nuclei (AGNs) phenomena, where the optical broad emission lines appear or disappear, is crucial to understand the evolution of the accretion flow around supermassive black holes. We present a strategy to select new CL candidates starting from a spectroscopic type 2 AGNs sample and searching for current type 1 photometric… ▽ More Determining the frequency and duration of changing--look (CL) active galactic nuclei (AGNs) phenomena, where the optical broad emission lines appear or disappear, is crucial to understand the evolution of the accretion flow around supermassive black holes. We present a strategy to select new CL candidates starting from a spectroscopic type 2 AGNs sample and searching for current type 1 photometric variability. We use the publicly available Zwicky Transient Facility (ZTF) alert stream and the Automatic Learning for the Rapid Classification of Events (ALeRCE) light curve classifier to produce a list of CL candidates with a highly automated algorithm, resulting in 60 candidates. Visual inspection reduced the sample to 30. We performed new spectroscopic observations of six candidates of our clean sample, without further refinement, finding the appearance of clear broad Balmer lines in four of them and tentative evidence of type changes in the remaining two, which suggests a promising success rate of $\geq66$ per cent for this CL selection method. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 5 pages

Photoreverberation mapping of quasars in the context of LSST observing strategies

Authors: Isidora Jankov, Andjelka B. Kovacevic, Dragana Ilic, Luka C. Popovic, Viktor Radovic, Iva Cvorovic-Hajdinjak, Robert Nikutta, Paula Sanchez-Saez

Abstract: The upcoming photometric surveys, such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) will monitor unprecedented number of active galactic nuclei (AGN) in a decade long campaign. Motivated by the science goals of LSST, which includes the harnessing of broadband light curves of AGN for photometric reverberation mapping (PhotoRM), we implement the existing formalism to estimate th… ▽ More The upcoming photometric surveys, such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) will monitor unprecedented number of active galactic nuclei (AGN) in a decade long campaign. Motivated by the science goals of LSST, which includes the harnessing of broadband light curves of AGN for photometric reverberation mapping (PhotoRM), we implement the existing formalism to estimate the lagged response of the emission line flux to the continuum variability using only mutli-band photometric light curves. We test the PhotoRM method on a set of 19 artificial light curves simulated using a stochastic model based on the Damped Random Walk process. These light curves are sampled using different observing strategies, including the two proposed by the LSST, in order to compare the accuracy of time-lag retrieval based on different observing cadences. Additionally, we apply the same procedure for time-lag retrieval to the observed photometric light curves of NGC 4395, and compare our results to the existing literature. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: The manuscript was accepted for publication in the Astronomische Nachrichten Special Issue

Conditional Neural Process for non-parametric modeling of AGN light curve

Authors: Iva Cvorovic-Hajdinjak, Andjelka B. Kovacevic, Dragana Ilic, Luka C. Popovic, Xinyu Dai, Isidora Jankov, Viktor Radovic, Paula Sanchez-Saez, Robert Nikutta

Abstract: The consequences of complex disturbed environments in the vicinity of a supermassive black hole are not well represented by standard statistical models of optical variability in active galactic nuclei (AGN). Thus, developing new methodologies for investigating and modeling AGN light curves is crucial. Conditional Neural Processes (CNPs) are nonlinear function models that forecast stochastic time-s… ▽ More The consequences of complex disturbed environments in the vicinity of a supermassive black hole are not well represented by standard statistical models of optical variability in active galactic nuclei (AGN). Thus, developing new methodologies for investigating and modeling AGN light curves is crucial. Conditional Neural Processes (CNPs) are nonlinear function models that forecast stochastic time-series based on a finite amount of known data without the use of any additional parameters or prior knowledge (kernels). We provide a CNP algorithm that is specifically designed for simulating AGN light curves. It was trained using data from the All-Sky Automated Survey for Supernovae, which included 153 AGN. We present CNP modeling performance for a subsample of five AGNs with distinctive difficult-to-model properties. The performance of CNP in predicting temporal flux fluctuation was assessed using a minimizing loss function, and the results demonstrated the algorithm's usefulness. Our preliminary parallelization experiments show that CNP can efficiently handle large amounts of data. These results imply that CNP can be more effective than standard tools in modeling large volumes of AGN data (as anticipated from time-domain surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time). △ Less

Submitted 18 November, 2021; originally announced November 2021.

Comments: The manuscript was accepted for publication in the Astronomische Nachrichten Special Issue

Searching for changing-state AGNs in massive datasets -- I: applying deep learning and anomaly detection techniques to find AGNs with anomalous variability behaviours

Authors: P. Sánchez-Sáez, H. Lira, L. Martí, N. Sánchez-Pi, J. Arredondo, F. E. Bauer, A. Bayo, G. Cabrera-Vives, C. Donoso-Oliva, P. A. Estévez, S. Eyheramendy, F. Förster, L. Hernández-García, A. M. Muñoz Arancibia, M. Pérez-Carrasco, M. Sepúlveda, J. R. Vergara

Abstract: The classic classification scheme for Active Galactic Nuclei (AGNs) was recently challenged by the discovery of the so-called changing-state (changing-look) AGNs (CSAGNs). The physical mechanism behind this phenomenon is still a matter of open debate and the samples are too small and of serendipitous nature to provide robust answers. In order to tackle this problem, we need to design methods that… ▽ More The classic classification scheme for Active Galactic Nuclei (AGNs) was recently challenged by the discovery of the so-called changing-state (changing-look) AGNs (CSAGNs). The physical mechanism behind this phenomenon is still a matter of open debate and the samples are too small and of serendipitous nature to provide robust answers. In order to tackle this problem, we need to design methods that are able to detect AGN right in the act of changing-state. Here we present an anomaly detection (AD) technique designed to identify AGN light curves with anomalous behaviors in massive datasets. The main aim of this technique is to identify CSAGN at different stages of the transition, but it can also be used for more general purposes, such as cleaning massive datasets for AGN variability analyses. We used light curves from the Zwicky Transient Facility data release 5 (ZTF DR5), containing a sample of 230,451 AGNs of different classes. The ZTF DR5 light curves were modeled with a Variational Recurrent Autoencoder (VRAE) architecture, that allowed us to obtain a set of attributes from the VRAE latent space that describes the general behaviour of our sample. These attributes were then used as features for an Isolation Forest (IF) algorithm, that is an anomaly detector for a "one class" kind of problem. We used the VRAE reconstruction errors and the IF anomaly score to select a sample of 8,809 anomalies. These anomalies are dominated by bogus candidates, but we were able to identify 75 promising CSAGN candidates. △ Less

Submitted 12 July, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: Accepted for publication in the Astronomical Journal (AJ)

Journal ref: AJ 162 206 (2021)

A random forest-based selection of optically variable AGN in the VST-COSMOS field

Authors: D. De Cicco, F. E. Bauer, M. Paolillo, S. Cavuoti, P. Sánchez-Sáez, W. N. Brandt, G. Pignata, M. Vaccari, M. Radovich

Abstract: The survey of the COSMOS field by the VLT Survey Telescope is an appealing testing ground for variability studies of active galactic nuclei (AGN). With 54 r-band visits over 3.3 yr and a single-visit depth of 24.6 r-band mag, the dataset is also particularly interesting in the context of performance forecasting for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This work is… ▽ More The survey of the COSMOS field by the VLT Survey Telescope is an appealing testing ground for variability studies of active galactic nuclei (AGN). With 54 r-band visits over 3.3 yr and a single-visit depth of 24.6 r-band mag, the dataset is also particularly interesting in the context of performance forecasting for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This work is the fifth in a series dedicated to the development of an automated, robust, and efficient methodology to identify optically variable AGN, aimed at deploying it on future LSST data. We test the performance of a random forest (RF) algorithm in selecting optically variable AGN candidates, investigating how the use of different AGN labeled sets (LSs) and features sets affects this performance. We define a heterogeneous AGN LS and choose a set of variability features and optical and near-infrared colors based on what can be extracted from LSST data. We find that an AGN LS that includes only Type I sources allows for the selection of a highly pure (91%) sample of AGN candidates, obtaining a completeness with respect to spectroscopically confirmed AGN of 69% (vs. 59% in our previous work). The addition of colors to variability features mildly improves the performance of the RF classifier, while colors alone prove less effective than variability in selecting AGN as they return contaminated samples of candidates and fail to identify most host-dominated AGN. We observe that a bright (r < 21 mag) AGN LS is able to retrieve candidate samples not affected by the magnitude cut, which is of great importance as faint AGN LSs for LSST-related studies will be hard to find and likely imbalanced. We estimate a sky density of 6.2 million AGN for the LSST main survey down to our current magnitude limit. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 25 pages, 14 figures; accepted for publication in A&A

Journal ref: A&A 645, A103 (2021)

Alert Classification for the ALeRCE Broker System: The Light Curve Classifier

Authors: P. Sánchez-Sáez, I. Reyes, C. Valenzuela, F. Förster, S. Eyheramendy, F. Elorrieta, F. E. Bauer, G. Cabrera-Vives, P. A. Estévez, M. Catelan, G. Pignata, P. Huijse, D. De Cicco, P. Arévalo, R. Carrasco-Davis, J. Abril, R. Kurtev, J. Borissova, J. Arredondo, E. Castillo-Navarrete, D. Rodriguez, D. Ruz-Mieres, A. Moya, L. Sabatini-Gacitúa, C. Sepúlveda-Cobo , et al. (1 additional authors not shown)

Abstract: We present the first version of the ALeRCE (Automatic Learning for the Rapid Classification of Events) broker light curve classifier. ALeRCE is currently processing the Zwicky Transient Facility (ZTF) alert stream, in preparation for the Vera C. Rubin Observatory. The ALeRCE light curve classifier uses variability features computed from the ZTF alert stream, and colors obtained from AllWISE and ZT… ▽ More We present the first version of the ALeRCE (Automatic Learning for the Rapid Classification of Events) broker light curve classifier. ALeRCE is currently processing the Zwicky Transient Facility (ZTF) alert stream, in preparation for the Vera C. Rubin Observatory. The ALeRCE light curve classifier uses variability features computed from the ZTF alert stream, and colors obtained from AllWISE and ZTF photometry. We apply a Balanced Random Forest algorithm with a two-level scheme, where the top level classifies each source as periodic, stochastic, or transient, and the bottom level further resolves each of these hierarchical classes, amongst 15 total classes. This classifier corresponds to the first attempt to classify multiple classes of stochastic variables (including core- and host-dominated active galactic nuclei, blazars, young stellar objects, and cataclysmic variables) in addition to different classes of periodic and transient sources, using real data. We created a labeled set using various public catalogs (such as the Catalina Surveys and {\em Gaia} DR2 variable stars catalogs, and the Million Quasars catalog), and we classify all objects with $\geq6$ $g$-band or $\geq6$ $r$-band detections in ZTF (868,371 sources as of 2020/06/09), providing updated classifications for sources with new alerts every day. For the top level we obtain macro-averaged precision and recall scores of 0.96 and 0.99, respectively, and for the bottom level we obtain macro-averaged precision and recall scores of 0.57 and 0.76, respectively. Updated classifications from the light curve classifier can be found at the \href{http://alerce.online}{ALeRCE Explorer website}. △ Less

Submitted 19 November, 2020; v1 submitted 7 August, 2020; originally announced August 2020.

Comments: 39 pages, 24 figures, 5 tables, 4 apendices. Accepted for publication in the Astronomical Journal (AJ)

Alert Classification for the ALeRCE Broker System: The Real-time Stamp Classifier

Authors: Rodrigo Carrasco-Davis, Esteban Reyes, Camilo Valenzuela, Francisco Förster, Pablo A. Estévez, Giuliano Pignata, Franz E. Bauer, Ignacio Reyes, Paula Sánchez-Sáez, Guillermo Cabrera-Vives, Susana Eyheramendy, Márcio Catelan, Javier Arredondo, Ernesto Castillo-Navarrete, Diego Rodríguez-Mancini, Daniela Ruz-Mieres, Alberto Moya, Luis Sabatini-Gacitúa, Cristóbal Sepúlveda-Cobo, Ashish A. Mahabal, Javier Silva-Farfán, Ernesto Camacho-Iñiquez, Lluís Galbany

Abstract: We present a real-time stamp classifier of astronomical events for the ALeRCE (Automatic Learning for the Rapid Classification of Events) broker. The classifier is based on a convolutional neural network, trained on alerts ingested from the Zwicky Transient Facility (ZTF). Using only the \textit{science, reference} and \textit{difference} images of the first detection as inputs, along with the met… ▽ More We present a real-time stamp classifier of astronomical events for the ALeRCE (Automatic Learning for the Rapid Classification of Events) broker. The classifier is based on a convolutional neural network, trained on alerts ingested from the Zwicky Transient Facility (ZTF). Using only the \textit{science, reference} and \textit{difference} images of the first detection as inputs, along with the metadata of the alert as features, the classifier is able to correctly classify alerts from active galactic nuclei, supernovae (SNe), variable stars, asteroids and bogus classes, with high accuracy ($\sim$94\%) in a balanced test set. In order to find and analyze SN candidates selected by our classifier from the ZTF alert stream, we designed and deployed a visualization tool called SN Hunter, where relevant information about each possible SN is displayed for the experts to choose among candidates to report to the Transient Name Server database. From June 26th 2019 to February 28th 2021, we have reported 6846 SN candidates to date (11.8 candidates per day on average), of which 971 have been confirmed spectroscopically. Our ability to report objects using only a single detection means that 70\% of the reported SNe occurred within one day after the first detection. ALeRCE has only reported candidates not otherwise detected or selected by other groups, therefore adding new early transients to the bulk of objects available for early follow-up. Our work represents an important milestone toward rapid alert classifications with the next generation of large etendue telescopes, such as the Vera C. Rubin Observatory. △ Less

Submitted 3 June, 2021; v1 submitted 7 August, 2020; originally announced August 2020.

Comments: Submitted to AAS on Jun 30th. Comments welcome

The Automatic Learning for the Rapid Classification of Events (ALeRCE) Alert Broker

Authors: F. Förster, G. Cabrera-Vives, E. Castillo-Navarrete, P. A. Estévez, P. Sánchez-Sáez, J. Arredondo, F. E. Bauer, R. Carrasco-Davis, M. Catelan, F. Elorrieta, S. Eyheramendy, P. Huijse, G. Pignata, E. Reyes, I. Reyes, D. Rodríguez-Mancini, D. Ruz-Mieres, C. Valenzuela, I. Alvarez-Maldonado, N. Astorga, J. Borissova, A. Clocchiatti, D. De Cicco, C. Donoso-Oliva, M. J. Graham , et al. (15 additional authors not shown)

Abstract: We introduce the Automatic Learning for the Rapid Classification of Events (ALeRCE) broker, an astronomical alert broker designed to provide a rapid and self--consistent classification of large etendue telescope alert streams, such as that provided by the Zwicky Transient Facility (ZTF) and, in the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). ALeRCE is a Chilean--l… ▽ More We introduce the Automatic Learning for the Rapid Classification of Events (ALeRCE) broker, an astronomical alert broker designed to provide a rapid and self--consistent classification of large etendue telescope alert streams, such as that provided by the Zwicky Transient Facility (ZTF) and, in the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). ALeRCE is a Chilean--led broker run by an interdisciplinary team of astronomers and engineers, working to become intermediaries between survey and follow--up facilities. ALeRCE uses a pipeline which includes the real--time ingestion, aggregation, cross--matching, machine learning (ML) classification, and visualization of the ZTF alert stream. We use two classifiers: a stamp--based classifier, designed for rapid classification, and a light--curve--based classifier, which uses the multi--band flux evolution to achieve a more refined classification. We describe in detail our pipeline, data products, tools and services, which are made public for the community (see \url{https://alerce.science}). Since we began operating our real--time ML classification of the ZTF alert stream in early 2019, we have grown a large community of active users around the globe. We describe our results to date, including the real--time processing of $9.7\times10^7$ alerts, the stamp classification of $1.9\times10^7$ objects, the light curve classification of $8.5\times10^5$ objects, the report of 3088 supernova candidates, and different experiments using LSST-like alert streams. Finally, we discuss the challenges ahead to go from a single-stream of alerts such as ZTF to a multi--stream ecosystem dominated by LSST. △ Less

Submitted 7 August, 2020; originally announced August 2020.

Comments: Submitted to AAS on Jun 29th. Preview for LSST PCW 2020. Comments welcome

The QUEST-La Silla AGN Variability Survey: selection of AGN candidates through optical variability

Authors: P. Sánchez-Sáez, P. Lira, R. Cartier, N. Miranda, L. C. Ho, P. Arévalo, F. E. Bauer, P. Coppi, C. Yovaniniz

Abstract: We used data from the QUEST-La Silla Active Galactic Nuclei (AGN) variability survey to construct light curves for 208,583 sources over $\sim 70$ deg$^2$, with a a limiting magnitude $r \sim 21$. Each light curve has at least 40 epochs and a length of $\geq 200$ days. We implemented a Random Forest algorithm to classify our objects as either AGN or non-AGN according to their variability features a… ▽ More We used data from the QUEST-La Silla Active Galactic Nuclei (AGN) variability survey to construct light curves for 208,583 sources over $\sim 70$ deg$^2$, with a a limiting magnitude $r \sim 21$. Each light curve has at least 40 epochs and a length of $\geq 200$ days. We implemented a Random Forest algorithm to classify our objects as either AGN or non-AGN according to their variability features and optical colors, excluding morphology cuts. We tested three classifiers, one that only includes variability features (RF1), one that includes variability features and also $r-i$ and $i-z$ colors (RF2), and one that includes variability features and also $g-r$, $r-i$, and $i-z$ colors (RF3). We obtained a sample of high probability candidates (hp-AGN) for each classifier, with 5,941 candidates for RF1, 5,252 candidates for RF2, and 4,482 candidates for RF3. We divided each sample according to their $g-r$ colors, defining blue ($g-r\leq 0.6$) and red sub-samples ($g-r>0.6$). We find that most of the candidates known from the literature belong to the blue sub-samples, which is not necessarily surprising given that, unlike for many literature studies, we do not cut our sample to point-like objects. This means that we can select AGN that have a significant contribution from redshifted starlight in their host galaxies. In order to test the efficiency of our technique we performed spectroscopic follow-up, confirming the AGN nature of 44 among 54 observed sources (81.5\% of efficiency). From the campaign we concluded that RF2 provides the purest sample of AGN candidates. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series

The QUEST-La Silla AGN variability survey: connection between AGN variability and black hole physical properties

Authors: P. Sánchez-Sáez, P. Lira, J. Mejía-Restrepo, L. C. Ho, P. Arévalo, M. Kim, R. Cartier, P. Coppi

Abstract: We present our statistical analysis of the connection between active galactic nuclei (AGN) variability and physical properties of the central supermassive black hole (SMBH). We constructed optical light curves using data from the QUEST-La Silla AGN variability survey. To model the variability, we used the structure function, among the excess variance and the amplitude from Damp Random Walk (DRW) m… ▽ More We present our statistical analysis of the connection between active galactic nuclei (AGN) variability and physical properties of the central supermassive black hole (SMBH). We constructed optical light curves using data from the QUEST-La Silla AGN variability survey. To model the variability, we used the structure function, among the excess variance and the amplitude from Damp Random Walk (DRW) modeling. For the measurement of SMBH physical properties, we used public spectra from the Sloan Digital Sky Survey (SDSS). Our analysis is based on an original sample of 2345 sources detected in both SDSS and QUEST-La Silla. For 1473 of these sources we could perform a proper measurement of the spectral and variability properties, and 1348 of these sources were classified as variable ($91.5\%$). We found that the amplitude of the variability ($A$) depends solely on the rest frame emission wavelength and the Eddington ratio, where $A$ anti-correlates with both $λ_{rest}$ and $L/L_{\text{Edd}}$. This suggests that AGN variability does not evolve over cosmic time, and its amplitude is inversely related to the accretion rate. We found that the logarithmic gradient of the variability ($γ$) does not correlate significantly with any SMBH physical parameter, since there is no statistically significant linear regression model with an absolute value of the slope higher than 0.1. Finally, we found that the general distribution of $γ$ measured for our sample differs from the distribution of $γ$ obtained for light curves simulated from a DRW process. For 20.6\% of the variable sources in our sample, a DRW model is not appropriate to describe the variability, since $γ$ differs considerably from the expected value of 0.5. △ Less

Submitted 27 August, 2018; v1 submitted 2 August, 2018; originally announced August 2018.

Comments: 26 pages, 17 figures, 9 tables. Minor formatting changes. Accepted for publication in ApJ