Showing 1–50 of 77 results for author: Förster, F

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

Temporal Stamp Classifier: Classifying Short Sequences of Astronomical Alerts

Authors: Daniel Neira O., Pablo A. Estévez, Francisco Förster

Abstract: In this work, we propose a deep learning-based classification model of astronomical objects using alerts reported by the Zwicky Transient Facility (ZTF) survey. The model takes as inputs sequences of stamp images and metadata contained in each alert, as well as features from the All-WISE catalog. The proposed model, called temporal stamp classifier, is able to discriminate between three classes of… ▽ More In this work, we propose a deep learning-based classification model of astronomical objects using alerts reported by the Zwicky Transient Facility (ZTF) survey. The model takes as inputs sequences of stamp images and metadata contained in each alert, as well as features from the All-WISE catalog. The proposed model, called temporal stamp classifier, is able to discriminate between three classes of astronomical objects: Active Galactic Nuclei (AGN), Super-Novae (SNe) and Variable Stars (VS), with an accuracy of approximately 98% in the test set, when using 2 to 5 detections. The results show that the model performance improves with the addition of more detections. Simple recurrence models obtain competitive results with those of more complex models such as LSTM.We also propose changes to the original stamp classifier model, which only uses the first detection. The performance of the latter model improves with changes in the architecture and the addition of random rotations, achieving a 1.46% increase in test accuracy. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: Accepted in International Joint Conference on Neural Networks 2024

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.

AstroInformatics: Recommendations for Global Cooperation

Authors: Ashish Mahabal, Pranav Sharma, Rana Adhikari, Mark Allen, Stefano Andreon, Varun Bhalerao, Federica Bianco, Anthony Brown, S. Bradley Cenko, Paula Coehlo, Jeffery Cooke, Daniel Crichton, Chenzhou Cui, Reinaldo de Carvalho, Richard Doyle, Laurent Eyer, Bernard Fanaroff, Christopher Fluke, Francisco Forster, Kevin Govender, Matthew J. Graham, Renée Hložek, Puji Irawati, Ajit Kembhavi, Juna Kollmeier , et al. (23 additional authors not shown)

Abstract: Policy Brief on "AstroInformatics, Recommendations for Global Collaboration", distilled from panel discussions during S20 Policy Webinar on Astroinformatics for Sustainable Development held on 6-7 July 2023. The deliberations encompassed a wide array of topics, including broad astroinformatics, sky surveys, large-scale international initiatives, global data repositories, space-related data, regi… ▽ More Policy Brief on "AstroInformatics, Recommendations for Global Collaboration", distilled from panel discussions during S20 Policy Webinar on Astroinformatics for Sustainable Development held on 6-7 July 2023. The deliberations encompassed a wide array of topics, including broad astroinformatics, sky surveys, large-scale international initiatives, global data repositories, space-related data, regional and international collaborative efforts, as well as workforce development within the field. These discussions comprehensively addressed the current status, notable achievements, and the manifold challenges that the field of astroinformatics currently confronts. The G20 nations present a unique opportunity due to their abundant human and technological capabilities, coupled with their widespread geographical representation. Leveraging these strengths, significant strides can be made in various domains. These include, but are not limited to, the advancement of STEM education and workforce development, the promotion of equitable resource utilization, and contributions to fields such as Earth Science and Climate Science. We present a concise overview, followed by specific recommendations that pertain to both ground-based and space data initiatives. Our team remains readily available to furnish further elaboration on any of these proposals as required. Furthermore, we anticipate further engagement during the upcoming G20 presidencies in Brazil (2024) and South Africa (2025) to ensure the continued discussion and realization of these objectives. The policy webinar took place during the G20 presidency in India (2023). Notes based on the seven panels will be separately published. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 7 pages

Domain Adaptation via Minimax Entropy for Real/Bogus Classification of Astronomical Alerts

Authors: Guillermo Cabrera-Vives, César Bolivar, Francisco Förster, Alejandra M. Muñoz Arancibia, Manuel Pérez-Carrasco, Esteban Reyes

Abstract: Time domain astronomy is advancing towards the analysis of multiple massive datasets in real time, prompting the development of multi-stream machine learning models. In this work, we study Domain Adaptation (DA) for real/bogus classification of astronomical alerts using four different datasets: HiTS, DES, ATLAS, and ZTF. We study the domain shift between these datasets, and improve a naive deep le… ▽ More Time domain astronomy is advancing towards the analysis of multiple massive datasets in real time, prompting the development of multi-stream machine learning models. In this work, we study Domain Adaptation (DA) for real/bogus classification of astronomical alerts using four different datasets: HiTS, DES, ATLAS, and ZTF. We study the domain shift between these datasets, and improve a naive deep learning classification model by using a fine tuning approach and semi-supervised deep DA via Minimax Entropy (MME). We compare the balanced accuracy of these models for different source-target scenarios. We find that both the fine tuning and MME models improve significantly the base model with as few as one labeled item per class coming from the target dataset, but that the MME does not compromise its performance on the source dataset. △ Less

Submitted 14 August, 2023; originally announced August 2023.

MSC Class: 85-08 ACM Class: I.2.1

Journal ref: ICML 2023 Workshop on Machine Learning for Astrophysics

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

A characterization of ASAS-SN core-collapse supernova environments with VLT+MUSE: I. Sample selection, analysis of local environments, and correlations with light curve properties

Authors: Thallis Pessi, Jose L. Prieto, Joseph P. Anderson, Lluís Galbany, Joseph D. Lyman, Christopher Kochanek, Subo Dong, Francisco Forster, Raul González-Díaz, Santiago Gonzalez-Gaitan, Claudia P. Gutiérrez, Thomas W. -S. Holoien, Philip A. James, Cristina Jiménez-Palau, Evelyn J. Johnston, Hanindyo Kuncarayakti, Fabián Rosales-Ortega, Sebastian F. Sánchez, Steve Schulze, Benjamin Shappee

Abstract: The analysis of core-collapse supernova (CCSN) environments can provide important information on the life cycle of massive stars and constrain the progenitor properties of these powerful explosions. The MUSE instrument at the VLT enables detailed local environment constraints of the progenitors of large samples of CCSNe. Using a homogeneous SN sample from the ASAS-SN survey has enabled us to perfo… ▽ More The analysis of core-collapse supernova (CCSN) environments can provide important information on the life cycle of massive stars and constrain the progenitor properties of these powerful explosions. The MUSE instrument at the VLT enables detailed local environment constraints of the progenitors of large samples of CCSNe. Using a homogeneous SN sample from the ASAS-SN survey has enabled us to perform a minimally biased statistical analysis of CCSN environments. We analyze 111 galaxies observed by MUSE that hosted 112 CCSNe detected or discovered by the ASAS-SN survey between 2014 and 2018. The majority of the galaxies were observed by the the AMUSING survey. Here we analyze the immediate environment around the SN locations and compare the properties between the different CCSN types and their light curves. We used stellar population synthesis and spectral fitting techniques to derive physical parameters for all HII regions detected within each galaxy, including the star formation rate (SFR), H$α$ equivalent width (EW), oxygen abundance, and extinction. We found that stripped-envelope (SE) SNe occur in environments with a higher median SFR, H$α$ EW, and oxygen abundances than SNe II and SNe IIn/Ibn. The distributions of SNe II and IIn are very similar, indicating that these events explode in similar environments. For the SESNe, SNe Ic have higher median SFRs, H$α$ EWs, and oxygen abundances than SNe Ib. SNe IIb have environments with similar SFRs and H$α$ EWs to SNe Ib, and similar oxygen abundances to SNe Ic. We also show that the postmaximum decline rate, $s$, of SNe II correlates with the H$α$ EW, and that the luminosity and the $Δm_{15}$ parameter of SESNe correlate with the oxygen abundance, H$α$ EW, and SFR at their environments. This suggests a connection between the explosion mechanisms of these events to their environment properties. △ Less

Submitted 28 September, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: 27 pages, 16 figures, accepted at A&A

Journal ref: A&A 677, A28 (2023)

Multi-scale stamps for real-time classification of alert streams

Authors: Ignacio Reyes-Jainaga, Francisco Förster, Alejandra M. Muñoz Arancibia, Guillermo Cabrera-Vives, Amelia Bayo, Franz E. Bauer, Javier Arredondo, Esteban Reyes, Giuliano Pignata, A. M. Mourão, Javier Silva-Farfán, Lluís Galbany, Alex Álvarez, Nicolás Astorga, Pablo Castellanos, Pedro Gallardo, Alberto Moya, Diego Rodríguez

Abstract: In recent years, automatic classifiers of image cutouts (also called "stamps") have shown to be key for fast supernova discovery. The Vera C. Rubin Observatory will distribute about ten million alerts with their respective stamps each night, enabling the discovery of approximately one million supernovae each year. A growing source of confusion for these classifiers is the presence of satellite gli… ▽ More In recent years, automatic classifiers of image cutouts (also called "stamps") have shown to be key for fast supernova discovery. The Vera C. Rubin Observatory will distribute about ten million alerts with their respective stamps each night, enabling the discovery of approximately one million supernovae each year. A growing source of confusion for these classifiers is the presence of satellite glints, sequences of point-like sources produced by rotating satellites or debris. The currently planned Rubin stamps will have a size smaller than the typical separation between these point sources. Thus, a larger field of view stamp could enable the automatic identification of these sources. However, the distribution of larger stamps would be limited by network bandwidth restrictions. We evaluate the impact of using image stamps of different angular sizes and resolutions for the fast classification of events (AGNs, asteroids, bogus, satellites, SNe, and variable stars), using data from the Zwicky Transient Facility. We compare four scenarios: three with the same number of pixels (small field of view with high resolution, large field of view with low resolution, and a multi-scale proposal) and a scenario with the full stamp that has a larger field of view and higher resolution. Compared to small field of view stamps, our multi-scale strategy reduces misclassifications of satellites as asteroids or supernovae, performing on par with high-resolution stamps that are 15 times heavier. We encourage Rubin and its Science Collaborations to consider the benefits of implementing multi-scale stamps as a possible update to the alert specification. △ Less

Submitted 14 July, 2023; v1 submitted 25 April, 2023; originally announced April 2023.

Comments: Accepted for publication in ApJL

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

Carnegie Supernova Project-II: Near-infrared spectral diversity and template of Type Ia Supernovae

Authors: Jing Lu, Eric Y. Hsiao, Mark M. Phillips, Christopher R. Burns, Chris Ashall, Nidia Morrell, Lawrence Ng, Sahana Kumar, Melissa Shahbandeh, Peter Hoeflich, E. Baron, Syed Uddin, Maximilian D. Stritzinger, Nicholas B. Suntzeff, Charles Baltay, Scott Davis, Tiara R. Diamond, Gaston Folatelli, Francisco Förster, Jonathan Gagné, Lluís Galbany, Christa Gall, Santiago González-Gaitán, Simon Holmbo, Robert P. Kirshner , et al. (8 additional authors not shown)

Abstract: We present the largest and most homogeneous collection of near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia): 339 spectra of 98 individual SNe obtained as part of the Carnegie Supernova Project-II. These spectra, obtained with the FIRE spectrograph on the 6.5 m Magellan Baade telescope, have a spectral range of 0.8--2.5 $μ$m. Using this sample, we explore the NIR spectral diversity of SNe… ▽ More We present the largest and most homogeneous collection of near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia): 339 spectra of 98 individual SNe obtained as part of the Carnegie Supernova Project-II. These spectra, obtained with the FIRE spectrograph on the 6.5 m Magellan Baade telescope, have a spectral range of 0.8--2.5 $μ$m. Using this sample, we explore the NIR spectral diversity of SNe Ia and construct a template of spectral time series as a function of the light-curve-shape parameter, color stretch $s_{BV}$. Principal component analysis is applied to characterize the diversity of the spectral features and reduce data dimensionality to a smaller subspace. Gaussian process regression is then used to model the subspace dependence on phase and light-curve shape and the associated uncertainty. Our template is able to predict spectral variations that are correlated with $s_{BV}$, such as the hallmark NIR features: Mg II at early times and the $H$-band break after peak. Using this template reduces the systematic uncertainties in K-corrections by ~90% compared to those from the Hsiao template. These uncertainties, defined as the mean K-correction differences computed with the color-matched template and observed spectra, are on the level of $4\times10^{-4}$ mag on average. This template can serve as the baseline spectral energy distribution for light-curve fitters and can identify peculiar spectral features that might point to compelling physics. The results presented here will substantially improve future SN~Ia cosmological experiments, for both nearby and distant samples. △ Less

Submitted 17 March, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: 38 pages, 21 figures, accepted to APJ

Near-infrared and Optical Nebular-phase Spectra of Type Ia Supernovae SN 2013aa and SN 2017cbv in NGC 5643

Authors: Sahana Kumar, Eric Y. Hsiao, Chris Ashall, Mark M. Phillips, Nidia Morrell, Peter Hoeflich, Chris R. Burns, Lluis Galbany, Eddie Baron, Carlos Contreras, Scott Davis, Tiara Diamond, Francisco Forster, Melissa L. Graham, Emir Karamehmetoglu, Robert P. Kirshner, Baerbel Koribalski, Kevin Krisciunas, Jing Lu, G. H. Marion, Priscila J. Pessi, Anthony L. Piro, Melissa Shahbandeh, Maximillian D Stritzinger, Nicholas B. Suntzeff , et al. (1 additional authors not shown)

Abstract: We present multi-wavelength time-series spectroscopy of SN 2013aa and SN 2017cbv, two Type Ia supernovae (SNe Ia) on the outskirts of the same host galaxy, NGC 5643. This work utilizes new nebular-phase near-infrared (NIR) spectra obtained by the Carnegie Supernova Project-II, in addition to previously published optical and NIR spectra. By measuring nebular-phase [Fe II] lines in both the optical… ▽ More We present multi-wavelength time-series spectroscopy of SN 2013aa and SN 2017cbv, two Type Ia supernovae (SNe Ia) on the outskirts of the same host galaxy, NGC 5643. This work utilizes new nebular-phase near-infrared (NIR) spectra obtained by the Carnegie Supernova Project-II, in addition to previously published optical and NIR spectra. By measuring nebular-phase [Fe II] lines in both the optical and NIR, we examine the explosion kinematics and test the efficacy of several emission line fitting techniques commonly used in the literature. The NIR [Fe II] 1.644 $μ$m line provides the most robust velocity measurements against variations due to the choice of the fit method and line blending. The resulting effects on velocity measurements due to choosing different fit methods, initial fit parameters, continuum and line profile functions, and fit region boundaries were also investigated. The NIR [Fe II] velocities yield the same radial shift direction as velocities measured using the optical [Fe II] 7155 A line, but the sizes of the shifts are consistently and substantially lower, pointing to a potential issue in optical studies. The NIR [Fe II] 1.644 $μ$m emission profile shows a lack of significant asymmetry in both SNe Ia, and the observed low velocities elevate the importance for correcting for any radial velocity contribution from the host galaxy's rotation. The low [Fe II] velocities measured in the NIR at nebular phases disfavors most progenitor scenarios in close double-degenerate systems for both SN 2013aa and SN 2017cbv. The time evolution of the NIR [Fe II] 1.644 $μ$m line also indicates moderately high progenitor white dwarf central density and potentially high magnetic fields. These sibling SNe Ia were well observed at both early and late times, providing an excellent opportunity to study the intrinsic diversity of SNe Ia. △ Less

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

Comments: submitted to ApJ on Oct 4, 2022 accepted for publication on Dec 19, 2022

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

Type II supernovae from the Carnegie Supernova Project-I. III. Understanding SN II diversity through correlations between physical and observed properties

Authors: L. Martinez, J. P. Anderson, M. C. Bersten, M. Hamuy, S. González-Gaitán, M. Orellana, M. Stritzinger, M. M. Phillips, C. P. Gutiérrez, C. Burns, T. de Jaeger, K. Ertini, G. Folatelli, F. Förster, L. Galbany, P. Hoeflich, E. Y. Hsiao, N. Morrell, P. J. Pessi, N. B. Suntzeff

Abstract: SNe II show great photometric and spectroscopic diversity which is attributed to the varied physical characteristics of their progenitor and explosion properties. In this study, the third of a series of papers where we analyse a sample of SNe II observed by the Carnegie Supernova Project-I, we present correlations between their observed and physical properties. Our analysis shows that explosion en… ▽ More SNe II show great photometric and spectroscopic diversity which is attributed to the varied physical characteristics of their progenitor and explosion properties. In this study, the third of a series of papers where we analyse a sample of SNe II observed by the Carnegie Supernova Project-I, we present correlations between their observed and physical properties. Our analysis shows that explosion energy is the physical property that correlates with the highest number of parameters. We recover previously suggested relationships between the hydrogen-rich envelope mass and the plateau duration, and find that more luminous SNe II with higher expansion velocities, faster declining light curves, and higher Ni masses are consistent with higher energy explosions. In addition, faster declining SNe II are also compatible with more concentrated Ni in the inner regions of the ejecta. Positive trends are found between the initial mass, explosion energy, and Ni mass. While the explosion energy spans the full range explored with our models, the initial mass generally arises from a relatively narrow range. Observable properties were measured from our grid of models to determine the effect of each physical parameter on the observed SN II diversity. We argue that explosion energy is the physical parameter causing the greatest impact on SN II diversity, when assuming standard single-star evolution as in the models used in this study. The inclusion of pre-SN models assuming higher mass loss produces a significant increase in the strength of some correlations, particularly those between the progenitor hydrogen-rich envelope mass and the plateau and optically thick phase durations. These differences clearly show the impact of having different treatments of stellar evolution, implying that changes in the assumption of standard single-star evolution are necessary for a complete understanding of SN II diversity. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Astronomy & Astrophysics. Paper I: arXiv:2111.06519; Paper II: arXiv:2111.06529

Journal ref: A&A 660, A42 (2022)

Deep Attention-Based Supernovae Classification of Multi-Band Light-Curves

Authors: Óscar Pimentel, Pablo A. Estévez, Francisco Förster

Abstract: In astronomical surveys, such as the Zwicky Transient Facility, supernovae (SNe) are relatively uncommon objects compared to other classes of variable events. Along with this scarcity, the processing of multi-band light-curves is a challenging task due to the highly irregular cadence, long time gaps, missing-values, few observations, etc. These issues are particularly detrimental to the analysis o… ▽ More In astronomical surveys, such as the Zwicky Transient Facility, supernovae (SNe) are relatively uncommon objects compared to other classes of variable events. Along with this scarcity, the processing of multi-band light-curves is a challenging task due to the highly irregular cadence, long time gaps, missing-values, few observations, etc. These issues are particularly detrimental to the analysis of transient events: SN-like light-curves. We offer three main contributions: 1) Based on temporal modulation and attention mechanisms, we propose a Deep attention model (TimeModAttn) to classify multi-band light-curves of different SN types, avoiding photometric or hand-crafted feature computations, missing-value assumptions, and explicit imputation/interpolation methods. 2) We propose a model for the synthetic generation of SN multi-band light-curves based on the Supernova Parametric Model, allowing us to increase the number of samples and the diversity of cadence. Thus, the TimeModAttn model is first pre-trained using synthetic light-curves. Then, a fine-tuning process is performed. The TimeModAttn model outperformed other Deep Learning models, based on Recurrent Neural Networks, in two scenarios: late-classification and early-classification. Also, the TimeModAttn model outperformed a Balanced Random Forest (BRF) classifier (trained with real data), increasing the balanced-$F_1$score from $\approx.525$ to $\approx.596$. When training the BRF with synthetic data, this model achieved similar performance to the TimeModAttn model proposed while still maintaining extra advantages. 3) We conducted interpretability experiments. High attention scores were obtained for observations earlier than and close to the SN brightness peaks. This also correlated with an early highly variability of the learned temporal modulation. △ Less

Submitted 25 November, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: Submitted to AJ on 14-Jan-2022

Type II supernovae from the Carnegie Supernova Project-I. II. Physical parameter distributions from hydrodynamical modelling

Authors: L. Martinez, M. C. Bersten, J. P. Anderson, M. Hamuy, S. González-Gaitán, F. Förster, M. Orellana, M. Stritzinger, M. M. Phillips, C. P. Gutiérrez, C. Burns, C. Contreras, T. de Jaeger, K. Ertini, G. Folatelli, L. Galbany, P. Hoeflich, E. Y. Hsiao, N. Morrell, P. J. Pessi, N. B. Suntzeff

Abstract: Linking supernovae to their progenitors is a powerful method for furthering our understanding of the physical origin of their observed differences, while at the same time testing stellar evolution theory. In this second study of a series of three papers where we characterise SNe II to understand their diversity, we derive progenitor properties (initial and ejecta masses, and radius), explosion ene… ▽ More Linking supernovae to their progenitors is a powerful method for furthering our understanding of the physical origin of their observed differences, while at the same time testing stellar evolution theory. In this second study of a series of three papers where we characterise SNe II to understand their diversity, we derive progenitor properties (initial and ejecta masses, and radius), explosion energy, $^{56}$Ni mass, and its degree of mixing within the ejecta for a large sample of SNe II. This dataset was obtained by the Carnegie Supernova Project-I and is characterised by a high cadence of their optical and NIR light curves and optical spectra that were homogeneously observed and processed. A large grid of hydrodynamical models and a fitting procedure based on MCMC methods were used to fit the bolometric light curve and the evolution of the photospheric velocity of 53 SNe II. We infer ejecta masses between 7.9 and 14.8 $M_{\odot}$, explosion energies between 0.15 and 1.40 foe, and $^{56}$Ni masses between 0.006 and 0.069 $M_{\odot}$. We define a subset of 24~SNe (the `gold sample') with well-sampled bolometric light curves and expansion velocities for which we consider the results more robust. Most SNe~II in the gold sample ($\sim$88%) are found with ejecta masses in the range of $\sim$8-10 $M_{\odot}$, coming from low zero-age main-sequence masses (9-12 $M_{\odot}$). The modelling of the initial-mass distribution of the gold sample gives an upper mass limit of 21.3$^{+3.8}_{-0.4}$ $M_{\odot}$ and a much steeper distribution than that for a Salpeter massive-star IMF. This IMF incompatibility is due to the large number of low-mass progenitors found -- when assuming standard stellar evolution. This may imply that high-mass progenitors lose more mass during their lives than predicted. However, a deeper analysis of all stellar evolution assumptions is required to test this hypothesis. △ Less

Submitted 8 February, 2022; v1 submitted 11 November, 2021; originally announced November 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 660, A41 (2022)

Type II supernovae from the Carnegie Supernova Project-I. I. Bolometric light curves of 74 SNe II using uBgVriYJH photometry

Authors: L. Martinez, M. C. Bersten, J. P. Anderson, M. Hamuy, S. González-Gaitán, M. Stritzinger, M. M. Phillips, C. P. Gutiérrez, C. Burns, C. Contreras, T. de Jaeger, K. Ertini, G. Folatelli, F. Förster, L. Galbany, P. Hoeflich, E. Y. Hsiao, N. Morrell, M. Orellana, P. J. Pessi, N. B. Suntzeff

Abstract: The present study is the first of a series of three papers where we characterise the type II supernovae (SNe~II) from the Carnegie Supernova Project-I to understand their diversity in terms of progenitor and explosion properties. In this first paper, we present bolometric light curves of 74 SNe~II. We outline our methodology to calculate the bolometric luminosity, which consists of the integration… ▽ More The present study is the first of a series of three papers where we characterise the type II supernovae (SNe~II) from the Carnegie Supernova Project-I to understand their diversity in terms of progenitor and explosion properties. In this first paper, we present bolometric light curves of 74 SNe~II. We outline our methodology to calculate the bolometric luminosity, which consists of the integration of the observed fluxes in numerous photometric bands ($uBgVriYJH$) and black-body (BB) extrapolations to account for the unobserved flux at shorter and longer wavelengths. BB fits were performed using all available broadband data except when line blanketing effects appeared. Photometric bands bluer than $r$ that are affected by line blanketing were removed from the fit, which makes near-infrared (NIR) observations highly important to estimate reliable BB extrapolations to the infrared. BB fits without NIR data produce notably different bolometric light curves, and therefore different estimates of SN~II progenitor and explosion properties when data are modelled. We present two methods to address the absence of NIR observations: (a) colour-colour relationships from which NIR magnitudes can be estimated using optical colours, and (b) new prescriptions for bolometric corrections as a function of observed SN~II colours. Using our 74 SN~II bolometric light curves, we provide a full characterisation of their properties based on several observed parameters. We measured magnitudes at different epochs, as well as durations and decline rates of different phases of the evolution. An analysis of the light-curve parameter distributions was performed, finding a wide range and a continuous sequence of observed parameters which is consistent with previous analyses using optical light curves. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 660, A40 (2022)

Carnegie Supernova Project: The First Homogeneous Sample of "Super-Chandrasekhar Mass"/2003fg-like Type Ia Supernova

Authors: C. Ashall, J. Lu, E. Y. Hsiao, P. Hoeflich, M. M. Phillips, L. Galbany, C. R. Burns, C. Contreras, K. Krisciunas, N. Morrell, M. D. Stritzinger, N. B. Suntzeff, F. Taddia, J. Anais, E. Baron, P. J. Brown, L. Busta, A. Campillay, S. Castellón, C. Corco, S. Davis, G. Folatelli, F. Forster, W. L. Freedman, C. Gonzaléz , et al. (16 additional authors not shown)

Abstract: We present a multi-wavelength photometric and spectroscopic analysis of thirteen "Super-Chandrasekhar Mass"/2003fg-like type Ia Supernova (SNe~Ia). Nine of these objects were observed by the Carnegie Supernova Project. 2003fg-like have slowly declining light curves ($Δm_{15}$(B) $<$1.3 mag), and peak absolute $B$-band magnitudes between $-19<M_{B}<-21$~mag. Many 2003fg-like are located in the same… ▽ More We present a multi-wavelength photometric and spectroscopic analysis of thirteen "Super-Chandrasekhar Mass"/2003fg-like type Ia Supernova (SNe~Ia). Nine of these objects were observed by the Carnegie Supernova Project. 2003fg-like have slowly declining light curves ($Δm_{15}$(B) $<$1.3 mag), and peak absolute $B$-band magnitudes between $-19<M_{B}<-21$~mag. Many 2003fg-like are located in the same part of the luminosity width relation as normal SNe~Ia. In the optical $B$ and $V$ bands, 2003fg-like look like normal SNe~Ia, but at redder wavelengths they diverge. Unlike other luminous SNe~Ia, 2003fg-like generally have only one $i$-band maximum which peaks after the epoch of $B$-band maximum, while their NIR light curve rise times can be $\gtrsim$40 days longer than those of normal SNe~Ia. They are also at least one magnitude brighter in the NIR bands than normal SNe~Ia, peaking above $M_H < -19$~mag, and generally have negative Hubble residuals, which may be the cause of some systematics in dark energy experiments. Spectroscopically, 2003fg-like exhibit peculiarities such as unburnt carbon well past maximum light, a large spread (8000--12000~km/s) in SiII $λ$6355 velocities at maximum light with no rapid early velocity decline, and no clear $H$-band break at +10~d, e. We find that SNe with a larger pseudo equivalent width of CII at maximum light have lower SiII $λ$6355 velocities and slower declining light curves. There are also multiple factors that contribute to the peak luminosity of 2003fg-like. The explosion of a C-O degenerate core inside a carbon-rich envelope is consistent with these observations. Such a configuration may come from the core degenerate scenario. △ Less

Submitted 31 August, 2021; v1 submitted 22 June, 2021; originally announced June 2021.

Comments: Accepted for publication in ApJ

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)

The G 305 Star-forming Region: II. Irregular variable stars

Authors: N. Medina, J. Borissova, R. Kurtev, J. Alonso-García, Carlos G. Román-Zúñiga, A. Bayo, Marina Kounkel, Alexandre Roman-Lopes, P. W. Lucas, K. R. Covey, Francisco Förster, Dante Minniti, Lucia Adame, Jesús Hernández

Abstract: We present a catalog of 167 newly discovered, irregular variables spanning a $\sim$7 deg${^2}$ area that encompasses the G 305 star-forming complex, one of the most luminous giant H II regions in the Galaxy. We aim to unveil and characterize the young stellar object (YSO) population of the region by analyzing the $K_{\rm s}$-band variability and $JHK_{\rm s}$ infrared colors from the {\it VISTA Va… ▽ More We present a catalog of 167 newly discovered, irregular variables spanning a $\sim$7 deg${^2}$ area that encompasses the G 305 star-forming complex, one of the most luminous giant H II regions in the Galaxy. We aim to unveil and characterize the young stellar object (YSO) population of the region by analyzing the $K_{\rm s}$-band variability and $JHK_{\rm s}$ infrared colors from the {\it VISTA Variables in the Vía Láctea} (VVV) survey. Additionally, SDSS-IV APOGEE-2 infrared spectra of selected objects are analyzed. The sample show relatively high amplitudes ($0.661<ΔK_{\rm S} <3.521$ mag). Most of them resemble sources with outbursts with amplitude $>1$ mag and duration longer than a few days, typically at least a year, known as {\it Eruptive Variables}. About 60% are likely to be Class II/Flat/I objects. This is also confirmed by the spectral index $α$ when available. From the analysis of APOGEE-2 near-infrared spectra of sources in the region, another 122 stars are classified as YSOs, and displays some infrared variability. The measured effective temperature $T_{\rm eff}$ peak is around 4000K and they are slightly super-solar in metal abundance. The modal radial velocity is approximately $-$41 km/s. Combining available catalogs of YSOs in the region with our data, we investigate the spatial distributions of 700 YSOs. They are clearly concentrated within the central cavity formed by the massive clusters Danks 1 and 2. The calculated surface density for the entire catalog is 0.025 YSOs/pc$^{-2}$, while the central cavity contains 10 times more objects per area (0.238 YSOs/pc$^{-2}$). △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 22 pages, 13 figures

Report number: AAS26957R1

Optical and Near-Infrared Observations of the Nearby Type Ia Supernova 2017cbv

Authors: Lingzhi Wang, Carlos Contreras, Maokai Hu, Mario A. Hamuy, Eric Y. Hsiao, David J. Sand, Joseph P. Anderson, Chris Ashall, Christopher R. Burns, Juncheng Chen, Tiara R. Diamond, Scott Davis, Francisco Förster, Lluís Galbany, Santiago González-Gaitán, Mariusz Gromadzki, Peter Hoeflich, Wenxiong Li, G. H. Marion, Nidia Morrell, Giuliano Pignata, Jose L. Prieto, Mark M. Phillips, Melissa Shahbandeh, Nicholas Suntzeff , et al. (6 additional authors not shown)

Abstract: Supernova (SN) 2017cbv in NGC 5643 is one of a handful of type Ia supernovae (SNe~Ia) reported to have excess blue emission at early times. This paper presents extensive $BVRIYJHK_s$-band light curves of SN 2017cbv, covering the phase from $-16$ to $+125$ days relative to $B$-band maximum light. SN 2017cbv reached a $B$-band maximum of 11.710$\pm$0.006~mag, with a post-maximum magnitude decline… ▽ More Supernova (SN) 2017cbv in NGC 5643 is one of a handful of type Ia supernovae (SNe~Ia) reported to have excess blue emission at early times. This paper presents extensive $BVRIYJHK_s$-band light curves of SN 2017cbv, covering the phase from $-16$ to $+125$ days relative to $B$-band maximum light. SN 2017cbv reached a $B$-band maximum of 11.710$\pm$0.006~mag, with a post-maximum magnitude decline $Δm_{15}(B)$=0.990$\pm$0.013 mag. The supernova suffered no host reddening based on Phillips intrinsic color, Lira-Phillips relation, and the CMAGIC diagram. By employing the CMAGIC distance modulus $μ=30.58\pm0.05$~mag and assuming $H_0$=72~$\rm km \ s^{-1} \ Mpc^{-1}$, we found that 0.73~\msun $^{56}$Ni was synthesized during the explosion of SN 2017cbv, which is consistent with estimates using reddening-free and distance-free methods via the phases of the secondary maximum of the NIR-band light curves. We also present 14 near-infrared spectra from $-18$ to $+49$~days relative to the $B$-band maximum light, providing constraints on the amount of swept-up hydrogen from the companion star in the context of the single degenerate progenitor scenario. No $Paβ$ emission feature was detected from our post-maximum NIR spectra, placing a hydrogen mass upper limit of 0.1 $M_{\odot}$. The overall optical/NIR photometric and NIR spectral evolution of SN 2017cbv is similar to that of a normal SN~Ia, even though its early evolution is marked by a flux excess no seen in most other well-observed normal SNe~Ia. We also compare the exquisite light curves of SN 2017cbv with some $M_{ch}$ DDT models and sub-$M_{ch}$ double detonation models. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: 20 figures, Accepted for publication in The Astrophysical Journal

The early discovery of SN 2017ahn: signatures of persistent interaction in a fast declining Type II supernova

Authors: L. Tartaglia, D. J. Sand, J. H. Groh, S. Valenti, S. D. Wyatt, K. A. Bostroem, P. J. Brown, S. Yang, J. Burke, T. -W. Chen, S. Davis, F. Förster, L. Galbany, J. Haislip, D. Hiramatsu, G. Hosseinzadeh, D. A. Howell, E. Y. Hsiao, S. W. Jha, V. Kouprianov, H. Kuncarayakti, J. D. Lyman, C. McCully, M. M. Phillips, A. Rau , et al. (3 additional authors not shown)

Abstract: We present high-cadence, comprehensive data on the nearby ($D\simeq33\,\rm{Mpc}$) Type II SN 2017ahn, discovered within $\sim$1 day of explosion, from the very early phases after explosion to the nebular phase. The observables of SN 2017ahn show a significant evolution over the $\simeq470\,\rm{d}$ of our follow-up campaign, first showing prominent, narrow Balmer lines and other high-ionization fea… ▽ More We present high-cadence, comprehensive data on the nearby ($D\simeq33\,\rm{Mpc}$) Type II SN 2017ahn, discovered within $\sim$1 day of explosion, from the very early phases after explosion to the nebular phase. The observables of SN 2017ahn show a significant evolution over the $\simeq470\,\rm{d}$ of our follow-up campaign, first showing prominent, narrow Balmer lines and other high-ionization features purely in emission (i.e. flash spectroscopy features), which progressively fade and lead to a spectroscopic evolution similar to that of more canonical Type II supernovae. Over the same period, the decline of the light curves in all bands is fast, resembling the photometric evolution of linearly declining H-rich core-collapse supernovae. The modeling of the light curves and early flash spectra suggest a complex circumstellar medium surrounding the progenitor star at the time of explosion, with a first dense shell produced during the very late stages of its evolution being swept up by the rapidly expanding ejecta within the first $\sim6\,\rm{d}$ of the supernova evolution, while signatures of interaction are observed also at later phases. Hydrodynamical models support the scenario in which linearly declining Type II supernovae are predicted to arise from massive yellow super/hyper giants depleted of most of their hydrogen layers. △ Less

Submitted 12 November, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: 26 pages (21+Appendices), 19 figures, 4 tables, accepted for publication on ApJ

Progenitor properties of type II supernovae: fitting to hydrodynamical models using Markov chain Monte Carlo methods

Authors: L. Martinez, M. C. Bersten, J. P. Anderson, S. González-Gaitán, F. Förster, G. Folatelli

Abstract: The progenitor and explosion properties of type II supernovae (SNe II) are fundamental to understand the evolution of massive stars. Special interest has been given to the range of initial masses of their progenitors, but despite the efforts made, it is still uncertain. Direct imaging of progenitors in pre-explosion images point out an upper initial mass cutoff of $\sim$18$M_{\odot}$. However, thi… ▽ More The progenitor and explosion properties of type II supernovae (SNe II) are fundamental to understand the evolution of massive stars. Special interest has been given to the range of initial masses of their progenitors, but despite the efforts made, it is still uncertain. Direct imaging of progenitors in pre-explosion images point out an upper initial mass cutoff of $\sim$18$M_{\odot}$. However, this is in tension with previous studies in which progenitor masses inferred by light curve modelling tend to favour high-mass solutions. Moreover, it has been argued that light curve modelling alone cannot provide a unique solution for the progenitor and explosion properties of SNe II. We develop a robust method which helps us to constrain the physical parameters of SNe II by fitting simultaneously their bolometric light curve and the evolution of the photospheric velocity to hydrodynamical models using statistical inference techniques. Pre-supernova red supergiant models were created using the stellar evolution code MESA, varying the initial progenitor mass. The explosion of these progenitors was then processed through hydrodynamical simulations, where the explosion energy, synthesised nickel mass, and the latter's spatial distribution within the ejecta were changed. We compare to observations via Markov chain Monte Carlo methods. We apply this method to a well-studied set of SNe with an observed progenitor in pre-explosion images and compare with results in the literature. Progenitor mass constraints are found to be consistent between our results and those derived by pre-SN imaging and the analysis of late-time spectral modelling. We have developed a robust method to infer progenitor and explosion properties of SN II progenitors which is consistent with other methods in the literature, which suggests that hydrodynamical modelling is able to accurately constrain physical properties of SNe II. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: Accepted for publication on Astronomy and Astrophysics

Journal ref: A&A 642, A143 (2020)

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

Studying Type II supernovae as cosmological standard candles using the Dark Energy Survey

Authors: T. de Jaeger, L. Galbany, S. González-Gaitán, R. Kessler, A. V. Filippenko, F. Förster, M. Hamuy, P. J. Brown, T. M. Davis, C. P. Gutiérrez, C. Inserra, G F. Lewis, A. Möller, D. Scolnic, M. Smith, D. Brout, D. Carollo, R. J. Foley, K. Glazebrook, S. R. Hinton, E. Macaulay, B. Nichol, M. Sako, N. E. Sommer, B. E. Tucker , et al. (63 additional authors not shown)

Abstract: Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lemaître law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility… ▽ More Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lemaître law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility of SNe II for deriving accurate extragalactic distances, which will be an asset for the next generation of telescopes where more-distant SNe II will be discovered. More specifically, we present a sample from the Dark Energy Survey Supernova Program (DES-SN) consisting of 15 SNe II with photometric and spectroscopic information spanning a redshift range up to 0.35. Combining our DES SNe with publicly available samples, and using the standard candle method (SCM), we construct the largest available Hubble diagram with SNe II in the Hubble flow (70 SNe II) and find an observed dispersion of 0.27 mag. We demonstrate that adding a colour term to the SN II standardisation does not reduce the scatter in the Hubble diagram. Although SNe II are viable as distance indicators, this work points out important issues for improving their utility as independent extragalactic beacons: find new correlations, define a more standard subclass of SNe II, construct new SN II templates, and dedicate more observing time to high-redshift SNe II. Finally, for the first time, we perform simulations to estimate the redshift-dependent distance-modulus bias due to selection effects. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 39 pages, 22 figures, 10 tables, Accepted for publication in MNRAS

Design and operation of the ATLAS Transient Science Server

Authors: K. W. Smith, S. J. Smartt, D. R. Young, J. L. Tonry, L. Denneau, H. Flewelling, A. N. Heinze, H. J. Weiland, B. Stalder, A. Rest, C. W. Stubbs, J. P. Anderson, T. -W. Chen, P. Clark, A. Do, F. Förster, M. Fulton, J. Gillanders, O. R. McBrien, D. O'Neill, S. Srivastav, D. E. Wright

Abstract: The Asteroid Terrestrial impact Last Alert System (ATLAS) system consists of two 0.5m Schmidt telescopes with cameras covering 29 square degrees at plate scale of 1.86 arcsec per pixel. Working in tandem, the telescopes routinely survey the whole sky visible from Hawaii (above $δ> -50^{\circ}$) every two nights, exposing four times per night, typically reaching $o < 19$ magnitude per exposure when… ▽ More The Asteroid Terrestrial impact Last Alert System (ATLAS) system consists of two 0.5m Schmidt telescopes with cameras covering 29 square degrees at plate scale of 1.86 arcsec per pixel. Working in tandem, the telescopes routinely survey the whole sky visible from Hawaii (above $δ> -50^{\circ}$) every two nights, exposing four times per night, typically reaching $o < 19$ magnitude per exposure when the moon is illuminated and $c < 19.5$ per exposure in dark skies. Construction is underway of two further units to be sited in Chile and South Africa which will result in an all-sky daily cadence from 2021. Initially designed for detecting potentially hazardous near earth objects, the ATLAS data enable a range of astrophysical time domain science. To extract transients from the data stream requires a computing system to process the data, assimilate detections in time and space and associate them with known astrophysical sources. Here we describe the hardware and software infrastructure to produce a stream of clean, real, astrophysical transients in real time. This involves machine learning and boosted decision tree algorithms to identify extragalactic and Galactic transients. Typically we detect 10-15 supernova candidates per night which we immediately announce publicly. The ATLAS discoveries not only enable rapid follow-up of interesting sources but will provide complete statistical samples within the local volume of 100 Mpc. A simple comparison of the detected supernova rate within 100 Mpc, with no corrections for completeness, is already significantly higher (factor 1.5 to 2) than the current accepted rates. △ Less

Submitted 2 June, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: 27 pages, 12 figures. Accepted for publication in PASP on 2020 May 15

Asteroids' Size Distribution and Colors from HiTS

Authors: J. Peña, C. Fuentes, F. Förster, J. Martínez-Palomera, G. Cabrera-Vives, J. C. Maureira, P. Huijse, P. A. Estévez, L. Galbany, S. González-Gaitán, Th. de Jaeger

Abstract: We report the observations of solar system objects during the 2015 campaign of the High cadence Transient Survey (HiTS). We found 5740 bodies (mostly Main Belt asteroids), 1203 of which were detected in different nights and in $g'$ and $r'$. Objects were linked in the barycenter system and their orbital parameters were computed assuming Keplerian motion. We identified 6 near Earth objects, 1738 Ma… ▽ More We report the observations of solar system objects during the 2015 campaign of the High cadence Transient Survey (HiTS). We found 5740 bodies (mostly Main Belt asteroids), 1203 of which were detected in different nights and in $g'$ and $r'$. Objects were linked in the barycenter system and their orbital parameters were computed assuming Keplerian motion. We identified 6 near Earth objects, 1738 Main Belt asteroids and 4 Trans-Neptunian objects. We did not find a $g'-r'$ color-size correlation for $14<H_{g'}<18$ ($1<D<10$ km) asteroids. We show asteroids' colors are disturbed by HiTS' 1.6 hour cadence and estimate that observations should be separated by at most 14 minutes to avoid confusion in future wide-field surveys like LSST. The size distribution for the Main Belt objects can be characterized as a simple power law with slope $\sim0.9$, steeper than in any other survey, while data from HiTS 2014's campaign is consistent with previous ones (slopes $\sim0.68$ at the bright end and $\sim0.34$ at the faint end). This difference is likely due to the ecliptic distribution of the Main Belt since 2015's campaign surveyed farther from the ecliptic than did 2014's and most previous surveys. △ Less

Submitted 13 March, 2020; v1 submitted 11 March, 2020; originally announced March 2020.

Comments: 17 pages, 18 figures

Journal ref: The Astronomical Journal, Volume 159, Number 4, Page 148, Year 2020

Introducing the Search for Intermediate-mass Black-hole In Nearby Galaxy (SIBLING) Survey

Authors: Jorge Martínez-Palomera, Paulina Lira, India Bhalla-Ladd, Francisco Förster, Richard M. Plotkin

Abstract: Intermediate mass black holes (IMBHs) have masses between the $10^2\!-\!10^6$ M$_\odot$ and are key to our understanding of the formation of massive black holes. The known population of IMBH remains small, with a few hundred candidates and only a handful of them confirmed as bona-fide IMBHs. Until now, the most widely used selection method is based on spectral analysis. Here we present a methodolo… ▽ More Intermediate mass black holes (IMBHs) have masses between the $10^2\!-\!10^6$ M$_\odot$ and are key to our understanding of the formation of massive black holes. The known population of IMBH remains small, with a few hundred candidates and only a handful of them confirmed as bona-fide IMBHs. Until now, the most widely used selection method is based on spectral analysis. Here we present a methodology to select IMBH candidates via optical variability analysis of the nuclear region of local galaxies ($z \leqslant 0.35$). Active IMBH accreting at low rates show small amplitude variability with time scales of hours, as it is seen in one of the known IMBH NGC4395. We found a sample of $\sim \!500$ galaxies evidencing fast and small amplitude variation in their weekly based light curves. We estimate an average occupancy fraction of 4\% and a surface density of $\sim \!3$ deg$^{-2}$, which represent an increase by a factor of $\sim\!40$ compared to previous searches. A large fraction ($78\%$) of the candidates are in spiral galaxies. We preliminary confirm the AGN nature of 22 sources via BPT diagrams using SDSS legacy spectra. Further confirmation of these candidates will require multiwavelength observations, especially in X-ray and radio bands. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: 16 pages (including Appendix), 15 figures

Enabling real-time multi-messenger astrophysics discoveries with deep learning

Authors: E. A. Huerta, Gabrielle Allen, Igor Andreoni, Javier M. Antelis, Etienne Bachelet, Bruce Berriman, Federica Bianco, Rahul Biswas, Matias Carrasco, Kyle Chard, Minsik Cho, Philip S. Cowperthwaite, Zachariah B. Etienne, Maya Fishbach, Francisco Förster, Daniel George, Tom Gibbs, Matthew Graham, William Gropp, Robert Gruendl, Anushri Gupta, Roland Haas, Sarah Habib, Elise Jennings, Margaret W. G. Johnson , et al. (35 additional authors not shown)

Abstract: Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravit… ▽ More Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: Invited Expert Recommendation for Nature Reviews Physics. The art work produced by E. A. Huerta and Shawn Rosofsky for this article was used by Carl Conway to design the cover of the October 2019 issue of Nature Reviews Physics

Journal ref: Nature Reviews Physics volume 1, pages 600-608 (2019)

Luminous Type II supernovae for their low expansion velocities

Authors: Ó. Rodríguez, G. Pignata, J. P. Anderson, T. J. Moriya, A. Clocchiatti, F. Förster, J. L. Prieto, M. M. Phillips, C. R. Burns, C. Contreras, G. Folatelli, C. P. Gutiérrez, M. Hamuy, N. I. Morrell, M. D. Stritzinger, N. B. Suntzeff, S. Benetti, E. Cappellaro, N. Elias-Rosa, A. Pastorello, M. Turatto, J. Maza, R. Antezana, R. Cartier, L. González , et al. (5 additional authors not shown)

Abstract: We present optical and near-IR data of three Type II supernovae (SNe II), SN 2008bm, SN 2009aj, and SN 2009au. These SNe display the following common characteristics: signs of early interaction of the ejecta with circumstellar material (CSM), blue $B\!-\!V$ colours, weakness of metal lines, low expansion velocities, and $V$-band absolute magnitudes 2-3 mag brighter than those expected for normal S… ▽ More We present optical and near-IR data of three Type II supernovae (SNe II), SN 2008bm, SN 2009aj, and SN 2009au. These SNe display the following common characteristics: signs of early interaction of the ejecta with circumstellar material (CSM), blue $B\!-\!V$ colours, weakness of metal lines, low expansion velocities, and $V$-band absolute magnitudes 2-3 mag brighter than those expected for normal SNe II based on their expansion velocities. Two more SNe reported in the literature (SN 1983K and LSQ13fn) share properties similar to our sample. Analysing this set of five SNe II, which are luminous for their low expansion velocities (LLEV), we find that their properties can be reproduced assuming ejecta-CSM interaction that lasts between 4-11 weeks post explosion. The contribution of this interaction to the radiation field seems to be the dominant component determining the observed weakness of metal lines in the spectra rather than the progenitor metallicity. Based on hydrodynamic simulations, we find that the interaction of the ejecta with a CSM of ~3.6 M$_\odot$ can reproduce the light curves and expansion velocities of SN 2009aj. Using data collected by the Chilean Automatic Supernova Search, we estimate an upper limit for the LLEV SNe II fraction to be 2-4 per cent of all normal SNe II. With the current data-set, it is not clear whether the LLEV events are a separated class of SNe II with a different progenitor system, or whether they are the extreme of a continuum mediated by CSM interaction with the rest of the normal SN II population. △ Less

Submitted 22 April, 2020; v1 submitted 9 October, 2019; originally announced October 2019.

Comments: 20 pages, 15 figures, Accepted for publication to MNRAS

Fast algorithms for slow moving asteroids: constraints on the distribution of Kuiper Belt Objects

Authors: Peter J. Whidden, J. Bryce Kalmbach, Andrew J. Connolly, R. Lynne Jones, Hayden Smotherman, Dino Bektesevic, Colin Slater, Andrew C. Becker, Željko Ivezić, Mario Jurić, Bryce Bolin, Joachim Moeyens, Francisco Förster, V. Zach Golkhou

Abstract: We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over… ▽ More We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over 10^10 possible asteroid trajectories in stacks of the order 10-15 4K x 4K images in under a minute using a single consumer grade GPU. We apply this algorithm to data from the 2015 campaign of the High Cadence Transient Survey (HiTS) obtained with the Dark Energy Camera (DECam). We find 39 previously unknown Kuiper Belt Objects in the 150 square degrees of the survey. Comparing these asteroids to an existing model for the inclination distribution of the Kuiper Belt we demonstrate that we recover a KBO population above our detection limit consistent with previous studies. Software used in this analysis is made available as an open source package. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: Accepted to AJ

Carnegie Supernova Project-II: The Near-infrared Spectroscopy Program

Authors: E. Y. Hsiao, M. M. Phillips, G. H. Marion, R. P. Kirshner, N. Morrell, D. J. Sand, C. R. Burns, C. Contreras, P. Hoeflich, M. D. Stritzinger, S. Valenti, J. P. Anderson, C. Ashall, C. Baltay, E. Baron, D. P. K. Banerjee, S. Davis, T. R. Diamond, G. Folatelli, Wendy L. Freedman, F. Förster, L. Galbany, C. Gall, S. González-Gaitán, A. Goobar , et al. (20 additional authors not shown)

Abstract: Shifting the focus of Type Ia supernova (SN Ia) cosmology to the near-infrared (NIR) is a promising way to significantly reduce the systematic errors, as the strategy minimizes our reliance on the empirical width-luminosity relation and uncertain dust laws. Observations in the NIR are also crucial for our understanding of the origins and evolution of these events, further improving their cosmologi… ▽ More Shifting the focus of Type Ia supernova (SN Ia) cosmology to the near-infrared (NIR) is a promising way to significantly reduce the systematic errors, as the strategy minimizes our reliance on the empirical width-luminosity relation and uncertain dust laws. Observations in the NIR are also crucial for our understanding of the origins and evolution of these events, further improving their cosmological utility. Any future experiments in the rest-frame NIR will require knowledge of the SN Ia NIR spectroscopic diversity, which is currently based on a small sample of observed spectra. Along with the accompanying paper, Phillips et al. (2018), we introduce the Carnegie Supernova Project-II (CSP-II), to follow up nearby SNe Ia in both the optical and the NIR. In particular, this paper focuses on the CSP-II NIR spectroscopy program, describing the survey strategy, instrumental setups, data reduction, sample characteristics, and future analyses on the data set. In collaboration with the Harvard-Smithsonian Center for Astrophysics (CfA) Supernova Group, we obtained 661 NIR spectra of 157 SNe Ia. Within this sample, 451 NIR spectra of 90 SNe Ia have corresponding CSP-II follow-up light curves. Such a sample will allow detailed studies of the NIR spectroscopic properties of SNe Ia, providing a different perspective on the properties of the unburned material, radioactive and stable nickel produced, progenitor magnetic fields, and searches for possible signatures of companion stars. △ Less

Submitted 3 December, 2018; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: 20 pages, 7 figures, accepted for publication in PASP

The delay of shock breakout due to circumstellar material seen in most Type II Supernovae

Authors: F. Förster, T. J. Moriya, J. C. Maureira, J. P. Anderson, S. Blinnikov, F. Bufano, G. Cabrera-Vives, A. Clocchiatti, Th. de Jaeger, P. A. Estévez, L. Galbany, S. González-Gaitán, G. Gräfener, M. Hamuy, E. Hsiao, P. Huentelemu, P. Huijse, H. Kuncarayakti, J. Martínez-Palomera, G. Medina, F. Olivares E., G. Pignata, A. Razza, I. Reyes, J. San Martín , et al. (13 additional authors not shown)

Abstract: Type II supernovae (SNe) originate from the explosion of hydrogen-rich supergiant massive stars. Their first electromagnetic signature is the shock breakout, a short-lived phenomenon which can last from hours to days depending on the density at shock emergence. We present 26 rising optical light curves of SN II candidates discovered shortly after explosion by the High cadence Transient Survey (HiT… ▽ More Type II supernovae (SNe) originate from the explosion of hydrogen-rich supergiant massive stars. Their first electromagnetic signature is the shock breakout, a short-lived phenomenon which can last from hours to days depending on the density at shock emergence. We present 26 rising optical light curves of SN II candidates discovered shortly after explosion by the High cadence Transient Survey (HiTS) and derive physical parameters based on hydrodynamical models using a Bayesian approach. We observe a steep rise of a few days in 24 out of 26 SN II candidates, indicating the systematic detection of shock breakouts in a dense circumstellar matter consistent with a mass loss rate $\dot{M} > 10^{-4} M_\odot yr^{-1}$ or a dense atmosphere. This implies that the characteristic hour timescale signature of stellar envelope SBOs may be rare in nature and could be delayed into longer-lived circumstellar material shock breakouts in most Type II SNe. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: Published in Nature Astronomy (https://www.nature.com/articles/s41550-018-0563-4). 41 pages including methods. 5 figures in main text) + 8 figures in methods

Journal ref: Nature Astronomy, 2018

The High Cadence Transient Survey (HITS): Compilation and characterization of light-curve catalogs

Authors: Jorge Martínez-Palomera, Francisco Förster, Pavlos Protopapas, Juan Carlos Maureira, Paulina Lira, Guillermo Cabrera-Vives, Pablo Huijse, Lluis Galbany, Thomas de Jaeger, Santiago González-Gaitán, Gustavo Medina, Giuliano Pignata, Jaime San Martín, Mario Hamuy, Ricardo R. Muñoz

Abstract: The High Cadence Transient Survey (HiTS) aims to discover and study transient objects with characteristic timescales between hours and days, such as pulsating, eclipsing and exploding stars. This survey represents a unique laboratory to explore large etendue observations from cadences of about 0.1 days and to test new computational tools for the analysis of large data. This work follows a fully \t… ▽ More The High Cadence Transient Survey (HiTS) aims to discover and study transient objects with characteristic timescales between hours and days, such as pulsating, eclipsing and exploding stars. This survey represents a unique laboratory to explore large etendue observations from cadences of about 0.1 days and to test new computational tools for the analysis of large data. This work follows a fully \textit{Data Science} approach: from the raw data to the analysis and classification of variable sources. We compile a catalog of ${\sim}15$ million object detections and a catalog of ${\sim}2.5$ million light-curves classified by variability. The typical depth of the survey is $24.2$, $24.3$, $24.1$ and $23.8$ in $u$, $g$, $r$ and $i$ bands, respectively. We classified all point-like non-moving sources by first extracting features from their light-curves and then applying a Random Forest classifier. For the classification, we used a training set constructed using a combination of cross-matched catalogs, visual inspection, transfer/active learning and data augmentation. The classification model consists of several Random Forest classifiers organized in a hierarchical scheme. The classifier accuracy estimated on a test set is approximately $97\%$. In the unlabeled data, $3\,485$ sources were classified as variables, of which $1\,321$ were classified as periodic. Among the periodic classes we discovered with high confidence, 1 $δ$-scutti, 39 eclipsing binaries, 48 rotational variables and 90 RR-Lyrae and for the non-periodic classes we discovered 1 cataclysmic variables, 630 QSO, and 1 supernova candidates. The first data release can be accessed in the project archive of HiTS. △ Less

Submitted 7 September, 2018; v1 submitted 3 September, 2018; originally announced September 2018.

Comments: 22 pages including 10 figures and 9 tables. Accepted for publication in AJ. For associated files, see http://astro.cmm.uchile.cl/HiTS/

Enhanced Rotational Invariant Convolutional Neural Network for Supernovae Detection

Authors: Esteban Reyes, Pablo A. Estévez, Ignacio Reyes, Guillermo Cabrera-Vives, Pablo Huijse, Rodrigo Carrasco-Davis, Francisco Förster

Abstract: In this paper, we propose an enhanced CNN model for detecting supernovae (SNe). This is done by applying a new method for obtaining rotational invariance that exploits cyclic symmetry. In addition, we use a visualization approach, the layer-wise relevance propagation (LRP) method, which allows finding the relevant pixels in each image that contribute to discriminate between SN candidates and artif… ▽ More In this paper, we propose an enhanced CNN model for detecting supernovae (SNe). This is done by applying a new method for obtaining rotational invariance that exploits cyclic symmetry. In addition, we use a visualization approach, the layer-wise relevance propagation (LRP) method, which allows finding the relevant pixels in each image that contribute to discriminate between SN candidates and artifacts. We introduce a measure to assess quantitatively the effect of the rotational invariant methods on the LRP relevance heatmaps. This allows comparing the proposed method, CAP, with the original Deep-HiTS model. The results show that the enhanced method presents an augmented capacity for achieving rotational invariance with respect to the original model. An ensemble of CAP models obtained the best results so far on the HiTS dataset, reaching an average accuracy of 99.53%. The improvement over Deep-HiTS is significant both statistically and in practice. △ Less

Submitted 10 August, 2018; originally announced August 2018.

Comments: 8 pages, 5 figures. Accepted for publication in proceedings of the IEEE World Congress on Computational Intelligence (IEEE WCCI), Rio de Janeiro, Brazil, 8-13 July, 2018

Seeing Double: ASASSN-18bt Exhibits a Two-Component Rise in the Early-Time K2 Light Curve

Authors: B. J. Shappee, T. W. -s. Holoien, M. R. Drout, K. Auchettl, M. D. Stritzinger, C. S. Kochanek, K. Z. Stanek, E. Shaya, G. Narayan, J. S. Brown, S. Bose, D. Bersier, J. Brimacombe, Ping Chen, Subo Dong, S. Holmbo, B. Katz, J. A. Munnoz, R. L. Mutel, R. S. Post, J. L. Prieto, J. Shields, D. Tallon, T. A. Thompson, P. J. Vallely , et al. (88 additional authors not shown)

Abstract: On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the A… ▽ More On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems. △ Less

Submitted 23 November, 2018; v1 submitted 30 July, 2018; originally announced July 2018.

Comments: 17 pages, 8 figures, 3 Tables. Accepted to ApJ. This work is part of a number of papers analyzing ASASSN-18bt, with coordinated papers from Dimitriadis et al. (2018) and Li et al. (2018)

Deep Learning for Image Sequence Classification of Astronomical Events

Authors: Rodrigo Carrasco-Davis, Guillermo Cabrera-Vives, Francisco Förster, Pablo A. Estévez, Pablo Huijse, Pavlos Protopapas, Ignacio Reyes, Jorge Martínez-Palomera, Cristóbal Donoso

Abstract: We propose a new sequential classification model for astronomical objects based on a recurrent convolutional neural network (RCNN) which uses sequences of images as inputs. This approach avoids the computation of light curves or difference images. This is the first time that sequences of images are used directly for the classification of variable objects in astronomy. The second contribution of th… ▽ More We propose a new sequential classification model for astronomical objects based on a recurrent convolutional neural network (RCNN) which uses sequences of images as inputs. This approach avoids the computation of light curves or difference images. This is the first time that sequences of images are used directly for the classification of variable objects in astronomy. The second contribution of this work is the image simulation process. We generate synthetic image sequences that take into account the instrumental and observing conditions, obtaining a realistic, set of movies for each astronomical object. The simulated dataset is used to train our RCNN classifier. This approach allows us to generate datasets to train and test our RCNN model for different astronomical surveys and telescopes. We aim at building a simulated dataset whose distribution is close enough to the real dataset, so that a fine tuning could match the distributions between real and simulated dataset. To test the RCNN classifier trained with the synthetic dataset, we used real-world data from the High cadence Transient Survey (HiTS) obtaining an average recall of 85%, improved to 94% after performing fine tuning with 10 real samples per class. We compare the results of our model with those of a light curve random forest classifier. The proposed RCNN with fine tuning has a similar performance on the HiTS dataset compared to the light curve classifier, trained on an augmented training set with 10 real samples per class. The RCNN approach presents several advantages in an alert stream classification scenario, such as a reduction of the data pre-processing, faster online evaluation and easier performance improvement using a few real data samples. These results encourage us to use this method for alert brokers systems that will process alert streams generated by new telescopes such as the Large Synoptic Survey Telescope. △ Less

Submitted 7 November, 2018; v1 submitted 10 July, 2018; originally announced July 2018.

Comments: 20 pages, 20 figures (corrected compilation errors). This is an Accepted Manuscript version of an article accepted for publication in Publications of the Astronomical Society of the Pacific. Nether the Astronomical Society of the Pacific nor IOP Publishing Ltd is responsible for any errors or omissions in this version of the manuscript or any version derived from it

Asteroids in the High cadence Transient Survey

Authors: J. Peña, C. Fuentes, F. Förster, J. C. Maureira, J. San Martín, J. Littín, P. Huijse, G. Cabrera-Vives, P. A. Estévez, L. Galbany, S. González-Gaitán, J. Martínez, Th. de Jaeger, M. Hamuy

Abstract: We report on the serendipitous observations of Solar System objects imaged during the High cadence Transient Survey (HiTS) 2014 observation campaign. Data from this high cadence, wide field survey was originally analyzed for finding variable static sources using Machine Learning to select the most-likely candidates. In this work we search for moving transients consistent with Solar System objects… ▽ More We report on the serendipitous observations of Solar System objects imaged during the High cadence Transient Survey (HiTS) 2014 observation campaign. Data from this high cadence, wide field survey was originally analyzed for finding variable static sources using Machine Learning to select the most-likely candidates. In this work we search for moving transients consistent with Solar System objects and derive their orbital parameters. We use a simple, custom detection algorithm to link trajectories and assume Keplerian motion to derive the asteroid's orbital parameters. We use known asteroids from the Minor Planet Center (MPC) database to assess the detection efficiency of the survey and our search algorithm. Trajectories have an average of nine detections spread over 2 days, and our fit yields typical errors of $σ_a\sim 0.07 ~{\rm AU}$, $σ_{\rm e} \sim 0.07 $ and $σ_i\sim 0.^{\circ}5~ {\rm deg}$ in semi-major axis, eccentricity, and inclination respectively for known asteroids in our sample. We extract 7,700 orbits from our trajectories, identifying 19 near Earth objects, 6,687 asteroids, 14 Centaurs, and 15 trans-Neptunian objects. This highlights the complementarity of supernova wide field surveys for Solar System research and the significance of machine learning to clean data of false detections. It is a good example of the data--driven science that LSST will deliver. △ Less

Submitted 8 June, 2018; originally announced June 2018.

Comments: 9 pages, 7 figures

Journal ref: The Astronomical Journal, Volume 155, Year 2018, Page 135

A Fast-Evolving, Luminous Transient Discovered by K2/Kepler

Authors: A. Rest, P. M. Garnavich, D. Khatami, D. Kasen, B. E. Tucker, E. J. Shaya, R. P. Olling, R. Mushotzky, A. Zenteno, S. Margheim, G. Strampelli, D. James, R. C. Smith, F. Förster, V. A. Villar

Abstract: For decades optical time-domain searches have been tuned to find ordinary supernovae, which rise and fall in brightness over a period of weeks. Recently, supernova searches have improved their cadences and a handful of fast-evolving luminous transients (FELTs) have been identified. FELTs have peak luminosities comparable to Type Ia supernovae, but rise to maximum in $<10$ days and fade from view i… ▽ More For decades optical time-domain searches have been tuned to find ordinary supernovae, which rise and fall in brightness over a period of weeks. Recently, supernova searches have improved their cadences and a handful of fast-evolving luminous transients (FELTs) have been identified. FELTs have peak luminosities comparable to Type Ia supernovae, but rise to maximum in $<10$ days and fade from view in $<$month. Here we present the most extreme example of this class thus far, KSN2015K, with a rise time of only 2.2 days and a time above half-maximum ($t_{1/2}$) of only 6.8 days. Here we show that, unlike Type Ia supernovae, the light curve of KSN2015K was not powered by the decay of radioactive elements. We further argue that it is unlikely that it was powered by continuing energy deposition from a central remnant (a magnetar or black hole). Using numerical radiation hydrodynamical models, we show that the light curve of KSN2015K is well fit by a model where the supernova runs into external material presumably expelled in a pre-supernova mass loss episode. The rapid rise of KSN2015K therefore probes the venting of photons when a hypersonic shock wave breaks out of a dense extended medium. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Comments: 29 pages, 4 figures, 3 supplemental figures, Nature Astronomy, in press

Type IIP supernova light curves affected by the acceleration of red supergiant winds

Authors: Takashi J. Moriya, Francisco Förster, Sung-Chul Yoon, Götz Gräfener, Sergei I. Blinnikov

Abstract: We introduce the first synthetic light-curve model set of Type IIP supernovae exploded within circumstellar media in which the acceleration of the red supergiant winds is taken into account. Because wind acceleration makes the wind velocities near the progenitors low, the density of the immediate vicinity of the red supergiant supernova progenitors can be higher than that extrapolated by using a c… ▽ More We introduce the first synthetic light-curve model set of Type IIP supernovae exploded within circumstellar media in which the acceleration of the red supergiant winds is taken into account. Because wind acceleration makes the wind velocities near the progenitors low, the density of the immediate vicinity of the red supergiant supernova progenitors can be higher than that extrapolated by using a constant terminal wind velocity. Therefore, even if the mass-loss rate of the progenitor is relatively low, it can have a dense circumstellar medium at the immediate stellar vicinity and the early light curves of Type IIP supernovae are significantly affected by it. We adopt a simple beta velocity law to formulate the wind acceleration. We provide bolometric and multicolor light curves of Type IIP supernovae exploding within such accelerated winds from the combinations of three progenitors, 12 - 16 Msun; five beta, 1-5; seven mass-loss rates, 1e-5 - 1e-2 Msun/yr; and four explosion energies, 0.5e51 - 2e51 erg. All the light curve models are available at https://goo.gl/o5phYb. When the circumstellar density is sufficiently high, our models do not show a classical shock breakout as a consequence of the interaction with the dense and optically-thick circumstellar media. Instead, they show a delayed 'wind breakout', substantially affecting early light curves of Type IIP supernovae. We find that the mass-loss rates of the progenitors need to be 1e-3 - 1e-2 Msun/yr to explain typical rise times of 5 - 10 days in Type IIP supernovae assuming a dense circumstellar radius of 1e15 cm. △ Less

Submitted 9 March, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: 12 pages, 9 figures, 2 tables, accepted by Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 476, Issue 2, p.2840-2851 (2018)

Discovery of distant RR Lyrae stars in the Milky Way using DECam

Authors: Gustavo Medina, Ricardo Muñoz, Katherina Vivas, Jeffrey Carlin, Francisco Förster, Jorge Martínez, Lluis Galbany, Santiago González-Gaitán, Mario Hamuy, Thomas de Jaeger, Juan Carlos Maureira, Jaime San Martín

Abstract: We report the discovery of distant RR Lyrae stars, including the most distant known in the Milky Way, using data taken in the $g-$band with the Dark Energy Camera as part of the High cadence Transient Survey (HiTS; 2014 campaign). We detect a total of 173 RR Lyrae stars over a ~120 deg^2 area, including both known RR Lyrae and new detections. The heliocentric distances d_H of the full sample range… ▽ More We report the discovery of distant RR Lyrae stars, including the most distant known in the Milky Way, using data taken in the $g-$band with the Dark Energy Camera as part of the High cadence Transient Survey (HiTS; 2014 campaign). We detect a total of 173 RR Lyrae stars over a ~120 deg^2 area, including both known RR Lyrae and new detections. The heliocentric distances d_H of the full sample range from 9 to >200 kpc, with 18 of them beyond 90 kpc. We identify three sub-groups of RR Lyrae as members of known systems: the Sextans dwarf spheroidal galaxy, for which we report 46 new discoveries, and the ultra-faint dwarf galaxies Leo IV and Leo V. Following an MCMC methodology, we fit spherical and ellipsoidal profiles of the form rho(R) ~ R^n to the radial density distribution of RR Lyrae in the Galactic halo. The best fit corresponds to the spherical case, for which we obtain a simple power-law index of n=-4.17^{+0.18}_{-0.20}, consistent with recent studies made with samples covering shorter distances. The pulsational properties of the outermost RR Lyrae in the sample (d_H>90 kpc) differ from the ones in the halo population at closer distances. The distribution of the stars in a Period-Amplitude diagram suggest they belong to Oosterhoff-intermediate or Oosterhoff II groups, similar to what is found in the ultra-faint dwarf satellites around the Milky Way. The new distant stars discovered represent an important addition to the few existing tracers of the Milky Way potential in the outer halo. △ Less

Submitted 5 February, 2018; originally announced February 2018.

Comments: Accepted for publication in The Astrophysical Journal

Robust period estimation using mutual information for multi-band light curves in the synoptic survey era

Authors: Pablo Huijse, Pablo A. Estevez, Francisco Forster, Scott F. Daniel, Andrew J. Connolly, Pavlos Protopapas, Rodrigo Carrasco, Jose C. Principe

Abstract: The Large Synoptic Survey Telescope (LSST) will produce an unprecedented amount of light curves using six optical bands. Robust and efficient methods that can aggregate data from multidimensional sparsely-sampled time series are needed. In this paper we present a new method for light curve period estimation based on the quadratic mutual information (QMI). The proposed method does not assume a part… ▽ More The Large Synoptic Survey Telescope (LSST) will produce an unprecedented amount of light curves using six optical bands. Robust and efficient methods that can aggregate data from multidimensional sparsely-sampled time series are needed. In this paper we present a new method for light curve period estimation based on the quadratic mutual information (QMI). The proposed method does not assume a particular model for the light curve nor its underlying probability density and it is robust to non-Gaussian noise and outliers. By combining the QMI from several bands the true period can be estimated even when no single-band QMI yields the period. Period recovery performance as a function of average magnitude and sample size is measured using 30,000 synthetic multi-band light curves of RR Lyrae and Cepheid variables generated by the LSST Operations and Catalog simulators. The results show that aggregating information from several bands is highly beneficial in LSST sparsely-sampled time series, obtaining an absolute increase in period recovery rate up to 50%. We also show that the QMI is more robust to noise and light curve length (sample size) than the multiband generalizations of the Lomb Scargle and Analysis of Variance periodograms, recovering the true period in 10-30% more cases than its competitors. A python package containing efficient Cython implementations of the QMI and other methods is provided. △ Less

Submitted 11 September, 2017; originally announced September 2017.

Comments: Accepted for publication ApJ Supplement Series: Special Issue on Solar/Stellar Astronomy Big Data

Serendipitous discovery of RR Lyrae stars in the Leo V ultra-faint galaxy

Authors: Gustavo E. Medina, Ricardo R. Muñoz, A. Katherina Vivas, Francisco Förster, Jeffrey L. Carlin, Jorge Martinez, Lluis Galbany, Santiago González-Gaitán, Mario Hamuy, Thomas de Jaeger, Juan Carlos Maureira, Jaime San Martín

Abstract: During the analysis of RR Lyrae stars discovered in the High cadence Transient Survey (HiTS) taken with the Dark Energy Camera at the 4-m telescope at Cerro Tololo Inter-American Observatory, we found a group of three very distant, fundamental mode pulsator RR Lyrae (type ab). The location of these stars agrees with them belonging to the Leo V ultra-faint satellite galaxy, for which no variable st… ▽ More During the analysis of RR Lyrae stars discovered in the High cadence Transient Survey (HiTS) taken with the Dark Energy Camera at the 4-m telescope at Cerro Tololo Inter-American Observatory, we found a group of three very distant, fundamental mode pulsator RR Lyrae (type ab). The location of these stars agrees with them belonging to the Leo V ultra-faint satellite galaxy, for which no variable stars have been reported to date. The heliocentric distance derived for Leo V based on these stars is 173 +/- 5 kpc. The pulsational properties (amplitudes and periods) of these stars locate them within the locus of the Oosterhoff II group, similar to most other ultra-faint galaxies with known RR Lyrae stars. This serendipitous discovery shows that distant RR Lyrae stars may be used to search for unknown faint stellar systems in the outskirts of the Milky Way. △ Less

Submitted 31 July, 2017; originally announced August 2017.

Comments: Accepted in ApJ Letters

A New, Faint Population of X-ray Transients

Authors: F. E. Bauer, E. Treister, K. Schawinski, S. Schulze, B. Luo, D. M. Alexander, W. N. Brandt, A. Comastri, F. Forster, R. Gilli, D. A. Kann, K. Maeda, K. Nomoto, M. Paolillo, P. Ranalli, D. P. Schneider, O. Shemmer, M. Tanaka, A. Tolstov, N. Tominaga, P. Tozzi, C. Vignali, J. Wang, Y. Xue, G. Yang

Abstract: We report on the detection of a remarkable new fast high-energy transient found in the Chandra Deep Field-South, robustly associated with a faint ($m_{\rm R}=27.5$ mag, $z_{\rm ph}$$\sim$2.2) host in the CANDELS survey. The X-ray event is comprised of 115$^{+12}_{-11}$ net 0.3-7.0 keV counts, with a light curve characterised by a $\approx$100 s rise time, a peak 0.3-10 keV flux of $\approx$5… ▽ More We report on the detection of a remarkable new fast high-energy transient found in the Chandra Deep Field-South, robustly associated with a faint ($m_{\rm R}=27.5$ mag, $z_{\rm ph}$$\sim$2.2) host in the CANDELS survey. The X-ray event is comprised of 115$^{+12}_{-11}$ net 0.3-7.0 keV counts, with a light curve characterised by a $\approx$100 s rise time, a peak 0.3-10 keV flux of $\approx$5$\times$10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, and a power-law decay time slope of $-1.53\pm0.27$. The average spectral slope is $Γ=1.43^{+0.23}_{-0.13}$, with no clear spectral variations. The \hbox{X-ray} and multi-wavelength properties effectively rule out the vast majority of previously observed high-energy transients. A few theoretical possibilities remain: an "orphan" X-ray afterglow from an off-axis short-duration Gamma-ray Burst (GRB) with weak optical emission; a low-luminosity GRB at high redshift with no prompt emission below $\sim$20 keV rest-frame; or a highly beamed Tidal Disruption Event (TDE) involving an intermediate-mass black hole and a white dwarf with little variability. However, none of the above scenarios can completely explain all observed properties. Although large uncertainties exist, the implied rate of such events is comparable to those of orphan and low-luminosity GRBs as well as rare TDEs, implying the discovery of an untapped regime for a known transient class, or a new type of variable phenomena whose nature remains to be determined. △ Less

Submitted 14 February, 2017; originally announced February 2017.

Comments: 18 pages, 7 figures, 4 tables. Accepted for publication in MNRAS

Continuum Foreground Polarization and Na~I Absorption in Type Ia SNe

Authors: P. Zelaya, A. Clocchiatti, D. Baade, P. Höflich, J. Maund, F. Patat, J. R. Quinn, E. Reilly, L. Wang, J. C. Wheele, F. Förster, S. González-Gaitán

Abstract: We present a study of the continuum polarization over the 400--600 nm range of 19 Type Ia SNe obtained with FORS at the VLT. We separate them in those that show Na I D lines at the velocity of their hosts and those that do not. Continuum polarization of the sodium sample near maximum light displays a broad range of values, from extremely polarized cases like SN 2006X to almost unpolarized ones lik… ▽ More We present a study of the continuum polarization over the 400--600 nm range of 19 Type Ia SNe obtained with FORS at the VLT. We separate them in those that show Na I D lines at the velocity of their hosts and those that do not. Continuum polarization of the sodium sample near maximum light displays a broad range of values, from extremely polarized cases like SN 2006X to almost unpolarized ones like SN 2011ae. The non--sodium sample shows, typically, smaller polarization values. The continuum polarization of the sodium sample in the 400--600 nm range is linear with wavelength and can be characterized by the mean polarization (P$_{\rm{mean}}$). Its values span a wide range and show a linear correlation with color, color excess, and extinction in the visual band. Larger dispersion correlations were found with the equivalent width of the Na I D and Ca II H & K lines, and also a noisy relation between P$_{\rm{mean}}$ and $R_{V}$, the ratio of total to selective extinction. Redder SNe show stronger continuum polarization, with larger color excesses and extinctions. We also confirm that high continuum polarization is associated with small values of $R_{V}$. The correlation between extinction and polarization -- and polarization angles -- suggest that the dominant fraction of dust polarization is imprinted in interstellar regions of the host galaxies. We show that Na I D lines from foreground matter in the SN host are usually associated with non-galactic ISM, challenging the typical assumptions in foreground interstellar polarization models. △ Less

Submitted 18 January, 2017; originally announced January 2017.

Comments: 46 pages, 17 figures, Accepted by ApJ

Deep-HiTS: Rotation Invariant Convolutional Neural Network for Transient Detection

Authors: Guillermo Cabrera-Vives, Ignacio Reyes, Francisco Förster, Pablo A. Estévez, Juan-Carlos Maureira

Abstract: We introduce Deep-HiTS, a rotation invariant convolutional neural network (CNN) model for classifying images of transients candidates into artifacts or real sources for the High cadence Transient Survey (HiTS). CNNs have the advantage of learning the features automatically from the data while achieving high performance. We compare our CNN model against a feature engineering approach using random f… ▽ More We introduce Deep-HiTS, a rotation invariant convolutional neural network (CNN) model for classifying images of transients candidates into artifacts or real sources for the High cadence Transient Survey (HiTS). CNNs have the advantage of learning the features automatically from the data while achieving high performance. We compare our CNN model against a feature engineering approach using random forests (RF). We show that our CNN significantly outperforms the RF model reducing the error by almost half. Furthermore, for a fixed number of approximately 2,000 allowed false transient candidates per night we are able to reduce the miss-classified real transients by approximately 1/5. To the best of our knowledge, this is the first time CNNs have been used to detect astronomical transient events. Our approach will be very useful when processing images from next generation instruments such as the Large Synoptic Survey Telescope (LSST). We have made all our code and data available to the community for the sake of allowing further developments and comparisons at https://github.com/guille-c/Deep-HiTS. △ Less

Submitted 2 January, 2017; originally announced January 2017.

Journal ref: The Astrophysical Journal, 2017

The High Cadence Transient Survey (HiTS) - I. Survey design and supernova shock breakout constraints

Authors: Francisco Förster, Juan C. Maureira, Jaime San Martín, Mario Hamuy, Jorge Martínez, Pablo Huijse, Guillermo Cabrera, Lluís Galbany, Thomas de Jaeger, Santiago González-Gaitán, Joseph P. Anderson, Hanindyo Kuncarayakti, Giuliano Pignata, Filomena Bufano, Jorge Littín, Felipe Olivares, Gustavo Medina, R. Chris Smith, A. Katherina Vivas, Pablo A. Estévez, Ricardo Muñoz, Eduardo Vera

Abstract: We present the first results of the High cadence Transient Survey (HiTS), a survey whose objective is to detect and follow up optical transients with characteristic timescales from hours to days, especially the earliest hours of supernova (SN) explosions. HiTS uses the Dark Energy Camera (DECam) and a custom made pipeline for image subtraction, candidate filtering and candidate visualization, whic… ▽ More We present the first results of the High cadence Transient Survey (HiTS), a survey whose objective is to detect and follow up optical transients with characteristic timescales from hours to days, especially the earliest hours of supernova (SN) explosions. HiTS uses the Dark Energy Camera (DECam) and a custom made pipeline for image subtraction, candidate filtering and candidate visualization, which runs in real-time to be able to react rapidly to the new transients. We discuss the survey design, the technical challenges associated with the real-time analysis of these large volumes of data and our first results. In our 2013, 2014 and 2015 campaigns we have detected more than 120 young SN candidates, but we did not find a clear signature from the short-lived SN shock breakouts (SBOs) originating after the core collapse of red supergiant stars, which was the initial science aim of this survey. Using the empirical distribution of limiting-magnitudes from our observational campaigns we measured the expected recovery fraction of randomly injected SN light curves which included SBO optical peaks produced with models from Tominaga et al. (2011) and Nakar & Sari (2010). From this analysis we cannot rule out the models from Tominaga et al. (2011) under any reasonable distributions of progenitor masses, but we can marginally rule out the brighter and longer-lived SBO models from Nakar & Sari (2010) under our best-guess distribution of progenitor masses. Finally, we highlight the implications of this work for future massive datasets produced by astronomical observatories such as LSST. △ Less

Submitted 12 September, 2016; originally announced September 2016.

Comments: 30 pages, 14 figures, accepted for publication in ApJ