Showing 1–50 of 271 results for author: Catelan, M

Unveiling the nature of two dwarf novae: CRTS J080846.2+313106 and V416 Dra

Authors: Arti Joshi, Márcio Catelan, Simone Scaringi, Axel Schwope, G. C. Anupama, Nikita Rawat, Devendra K. Sahu, Mridweeka Singh, Raya Dastidar, Rama Venkata Subramanian, Srinivas M Rao

Abstract: We present the analysis of optical photometric and spectroscopic observations of two non-magnetic cataclysmic variables, namely CRTS J080846.2+313106 and V416 Dra. CRTS J080846.2+313106 has been found to vary with a period of 4.9116$\pm$0.0003 h, which was not found in earlier studies and is provisionally suggested as the orbital period of the system. In both long-period systems, the observed domi… ▽ More We present the analysis of optical photometric and spectroscopic observations of two non-magnetic cataclysmic variables, namely CRTS J080846.2+313106 and V416 Dra. CRTS J080846.2+313106 has been found to vary with a period of 4.9116$\pm$0.0003 h, which was not found in earlier studies and is provisionally suggested as the orbital period of the system. In both long-period systems, the observed dominant signal at second harmonic of the orbital frequency and the orbital modulation during quiescence are suggestive of ellipsoidal variation from changing aspects of the secondary, with an additional contribution from the accretion stream or hotspot. However, during the outburst, the hotspot itself is overwhelmed by the increased brightness, which is possibly associated with the accretion disc. The mid-eclipse phase for V416 Dra occurs earlier and the width of the eclipse is greater during outbursts compared to quiescence, suggesting an increased accretion disc radius during outbursts. Furthermore, the investigation of accretion disc eclipse in V416 Dra implies that a total disc eclipse is possible during quiescence, whereas the disc seems to be partially obscured during outbursts, which further signifies that the disc may grow in size as the outburst progresses. Optical spectra of CRTS J080846.2+313106 and V416 Dra are typical of dwarf novae during quiescence, and they both show a significant contribution from the M2-4V secondary. The light curve patterns, orbital periods, and spectra observed in both systems look remarkably similar, and seem to resemble the characteristics of U Gem-type dwarf novae. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 14 pages, 11 Figures, and 3 Tables, Accepted for publication in A&A

Digging deeper into the dense Galactic globular cluster Terzan 5 with Electron-Multiplying CCDs. Variable star detection and new discoveries

Authors: R. Figuera Jaimes, M. Catelan, K. Horne, J. Skottfelt, C. Snodgrass, M. Dominik, U. G. Jørgensen, J. Southworth, M. Hundertmark, P. Longa-Peña, S. Sajadian, J. Tregolan-Reed, T. C. Hinse, M. I. Andersen, M. Bonavita, V. Bozza, M. J. Burgdorf, L. Haikala, E. Khalouei, H. Korhonen, N. Peixinho, M. Rabus, S. Rahvar

Abstract: Context. High frame-rate imaging was employed to mitigate the effects of atmospheric turbulence (seeing) in observations of globular cluster Terzan 5. Aims. High-precision time-series photometry has been obtained with the highest angular resolution so far taken in the crowded central region of Terzan 5, with ground-based telescopes, and ways to avoid saturation of the brightest stars in the fiel… ▽ More Context. High frame-rate imaging was employed to mitigate the effects of atmospheric turbulence (seeing) in observations of globular cluster Terzan 5. Aims. High-precision time-series photometry has been obtained with the highest angular resolution so far taken in the crowded central region of Terzan 5, with ground-based telescopes, and ways to avoid saturation of the brightest stars in the field observed. Methods. The Electron-Multiplying Charge Coupled Device (EMCCD) camera installed at the Danish 1.54-m telescope at the ESO La Silla Observatory was employed to produce thousands of short-exposure time images (ten images per second) that were stacked to produce the normal-exposure-time images (minutes). We employed difference image analysis in the stacked images to produce high-precision photometry using the DanDIA pipeline. Results. Light curves of 1670 stars with 242 epochs were analyzed in the crowded central region of Terzan 5 to statistically detect variable stars in the field observed. We present a possible visual counterpart outburst at the position of the pulsar J1748-2446N, and the visual counterpart light curve of the low-mass X-ray binary CX 3. Additionally, we present the discovery of 4 semiregular variables. We also present updated ephemerides and properties of the only RR Lyrae star previously known in the field covered by our observations in Terzan 5. Finally, we report a significant displacement of two sources by ~0.62 and 0.59 arcseconds with respect to their positions in previous images available in the literature. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 22 pages, 18 figures, 7 tables. Accepted for publication in A&A

The VISTA Variables in the Vía Láctea eXtended (VVVX) ESO public survey: Completion of the observations and legacy

Authors: R. K. Saito, M. Hempel, J. Alonso-García, P. W. Lucas, D. Minniti, S. Alonso, L. Baravalle, J. Borissova, C. Caceres, A. N. Chené, N. J. G. Cross, F. Duplancic, E. R. Garro, M. Gómez, V. D. Ivanov, R. Kurtev, A. Luna, D. Majaess, M. G. Navarro, J. B. Pullen, M. Rejkuba, J. L. Sanders, L. C. Smith, P. H. C. Albino, M. V. Alonso , et al. (121 additional authors not shown)

Abstract: The ESO public survey VISTA Variables in the Vía Láctea (VVV) surveyed the inner Galactic bulge and the adjacent southern Galactic disk from $2009-2015$. Upon its conclusion, the complementary VVV eXtended (VVVX) survey has expanded both the temporal as well as spatial coverage of the original VVV area, widening it from $562$ to $1700$ sq. deg., as well as providing additional epochs in… ▽ More The ESO public survey VISTA Variables in the Vía Láctea (VVV) surveyed the inner Galactic bulge and the adjacent southern Galactic disk from $2009-2015$. Upon its conclusion, the complementary VVV eXtended (VVVX) survey has expanded both the temporal as well as spatial coverage of the original VVV area, widening it from $562$ to $1700$ sq. deg., as well as providing additional epochs in $JHK_{\rm s}$ filters from $2016-2023$. With the completion of VVVX observations during the first semester of 2023, we present here the observing strategy, a description of data quality and access, and the legacy of VVVX. VVVX took $\sim 2000$ hours, covering about 4% of the sky in the bulge and southern disk. VVVX covered most of the gaps left between the VVV and the VISTA Hemisphere Survey (VHS) areas and extended the VVV time baseline in the obscured regions affected by high extinction and hence hidden from optical observations. VVVX provides a deep $JHK_{\rm s}$ catalogue of $\gtrsim 1.5\times10^9$ point sources, as well as a $K_{\rm s}$ band catalogue of $\sim 10^7$ variable sources. Within the existing VVV area, we produced a $5D$ map of the surveyed region by combining positions, distances, and proper motions of well-understood distance indicators such as red clump stars, RR Lyrae, and Cepheid variables. In March 2023 we successfully finished the VVVX survey observations that started in 2016, an accomplishment for ESO Paranal Observatory upon 4200 hours of observations for VVV+VVVX. The VVV+VVVX catalogues complement those from the Gaia mission at low Galactic latitudes and provide spectroscopic targets for the forthcoming ESO high-multiplex spectrographs MOONS and 4MOST. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 17 pages, 11 figures (+ appendix). Accepted for publication in Astronomy and Astrophysics in section 14: Catalogs and data

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

A study of nine extreme low mass ratio contact binary systems

Authors: E. Lalounta, P. -E. Christopoulou, A. Papageorgiou, C. E. Ferreira Lopes, M. Catelan

Abstract: Low mass ratio systems (LMR) are a very interesting class of contact eclipsing binaries challenging the theoretical models of stability. These systems are also considered possible progenitors of the rare low-mass optical transients called red novae. In this study, we present the identification of 7 new totally eclipsing LMR systems from Catalina Sky Surveys (CSS) and 77 LMR candidates from the All… ▽ More Low mass ratio systems (LMR) are a very interesting class of contact eclipsing binaries challenging the theoretical models of stability. These systems are also considered possible progenitors of the rare low-mass optical transients called red novae. In this study, we present the identification of 7 new totally eclipsing LMR systems from Catalina Sky Surveys (CSS) and 77 LMR candidates from the All Sky Automated Survey (ASAS$-$3). Using the available CSS light curves and new multiband observations for CSS$\_$J210228.3$-$031048 and CSS$\_$J231513.$+$345335 with the 2.3 m Aristarchos telescope at Helmos Observatory, we estimate their physical and absolute parameters and investigate their stability and their progenitors. The light curves are analyzed by performing a 2$-$dimension scan on the mass ratio $-$ inclination plane with Phoebe$-$0.31 scripter while the errors are estimated using Monte-Carlo simulations and heuristic scanning of the parameter space. Our analysis revealed that all 9 CSS systems have extreme mass ratios from 0.09 to 0.16. Our statistical analysis of well$-$studied LMR contact binaries shows that LMR systems tend to have warmer and more massive primaries. The investigation of the progenitors of both low and higher-mass ratio systems reveals a trend for the former to originate from higher-mass ancestors. Finally, we investigate the stability condition by calculating the ratio of spin angular momentum to orbital angular momentum and other stability indicators in the context of the reliability of the solutions. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures

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.

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

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

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

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

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

The most variable VVV sources: eruptive protostars, dipping giants in the Nuclear Disc and others

Authors: P. W. Lucas, L. C. Smith, Z. Guo, C. Contreras Peña, D. Minniti, N. Miller, J. Alonso-García, M. Catelan, J. Borissova, R. K. Saito, R. Kurtev, M. G. Navarro, C. Morris, H. Muthu, D. Froebrich, V. D. Ivanov, A. Bayo, A. Caratti o Garatti, J. L. Sanders

Abstract: We have performed a comprehensive search of a VISTA Variables in the Via Lactea (VVV) database of 9.5 yr light curves for variable sources with $ΔK_s \ge 4$ mag, aiming to provide a large sample of high amplitude eruptive young stellar objects (YSOs) and detect unusual or new types of infrared variable source. We find 222 variable or transient sources in the Galactic bulge and disc, most of which… ▽ More We have performed a comprehensive search of a VISTA Variables in the Via Lactea (VVV) database of 9.5 yr light curves for variable sources with $ΔK_s \ge 4$ mag, aiming to provide a large sample of high amplitude eruptive young stellar objects (YSOs) and detect unusual or new types of infrared variable source. We find 222 variable or transient sources in the Galactic bulge and disc, most of which are new discoveries. The sample mainly comprises novae, YSOs, microlensing events, Long Period Variable stars (LPVs) and a few rare or unclassified sources. Additionally, we report the discovery of a significant population of aperiodic late-type giant stars suffering deep extinction events, strongly clustered in the Nuclear Disc of the Milky Way. We suggest that these are metal-rich stars in which radiatively driven mass loss has been enhanced by super-solar metallicity. Among the YSOs, 32/40 appear to be undergoing episodic accretion. Long-lasting YSO eruptions have a typical rise time of $\sim$2 yr, somewhat slower than the 6-12 month timescale seen in the few historical events observed on the rise. The outburst durations are usually at least 5 yr, somewhat longer than many lower amplitude VVV events detected previously. The light curves are diverse in nature, suggesting that multiple types of disc instability may occur. Eight long-duration extinction events are seen wherein the YSO dims for a year or more, attributable to inner disc structure. One binary YSO in NGC 6530 displays periodic extinction events (P=59 days) similar to KH 15D. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: Accepted by MNRAS. 38 pages. Full source lists and light curves are available in the Ancillary Files

Spectroscopic confirmation of high-amplitude eruptive YSOs and dipping giants from the VVV survey

Authors: Zhen Guo, P. W. Lucas, R. Kurtev, J. Borissova, C. Contreras Peña, S. N. Yurchenko, L. C. Smith, D. Minniti, R. K. Saito, A. Bayo, M. Catelan, J. Alonso-García, A. Caratti o Garatti, C. Morris, D. Froebrich, J. Tennyson, K. Maucó, A. Aguayo, N. Miller, H. D. S. Muthu

Abstract: During the pre-main-sequence (pre-MS) evolution stage of a star, significant amounts of stellar mass are accreted during episodic accretion events, such as multi-decade FUor-type outbursts. Here, we present a near-infrared spectroscopic follow-up study of 33 high-amplitude (most with $ΔK_s$ > 4 mag) variable sources discovered by the Vista Variables in the Via Lactea (VVV) survey. Based on the spe… ▽ More During the pre-main-sequence (pre-MS) evolution stage of a star, significant amounts of stellar mass are accreted during episodic accretion events, such as multi-decade FUor-type outbursts. Here, we present a near-infrared spectroscopic follow-up study of 33 high-amplitude (most with $ΔK_s$ > 4 mag) variable sources discovered by the Vista Variables in the Via Lactea (VVV) survey. Based on the spectral features, 25 sources are classified as eruptive young stellar objects (YSOs), including 15 newly identified FUors, six with long-lasting but EXor-like bursts of magnetospheric accretion and four displaying outflow-dominated spectra. By examining the photometric behaviours of eruptive YSOs, we found most FUor-type outbursts have higher amplitudes ($ΔK_s$ and $ΔW2$), faster eruptive timescales and bluer infrared colours than the other outburst types. In addition, we identified seven post-main sequence variables apparently associated with deep dipping events and an eruptive star with deep AlO absorption bands resembling those seen in the V838 Mon stellar merger. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 20 pages, 16 figures

Asteroseismological analysis of the polluted ZZ Ceti star G29-38 with TESS

Authors: Murat Uzundag, Francisco C. De Gerónimo, Alejandro H. Córsico, Roberto Silvotti, Paul A. Bradley, Michael H. Montgomery, Márcio Catelan, Odette Toloza, Keaton J. Bell, S. O. Kepler, Leandro G. Althaus, Scot J. Kleinman, Mukremin Kilic, Susan E. Mullally, Boris T. Gänsicke, Karolina Bąkowska, Sam Barber, Atsuko Nitta

Abstract: G\,29$-$38 (TIC~422526868) is one of the brightest ($V=13.1$) and closest ($d = 17.51$\,pc) pulsating white dwarfs with a hydrogen-rich atmosphere (DAV/ZZ Ceti class). It was observed by the {\sl TESS} spacecraft in sectors 42 and 56. The atmosphere of G~29$-$38 is polluted by heavy elements that are expected to sink out of visible layers on short timescales. The photometric {\sl TESS} data set sp… ▽ More G\,29$-$38 (TIC~422526868) is one of the brightest ($V=13.1$) and closest ($d = 17.51$\,pc) pulsating white dwarfs with a hydrogen-rich atmosphere (DAV/ZZ Ceti class). It was observed by the {\sl TESS} spacecraft in sectors 42 and 56. The atmosphere of G~29$-$38 is polluted by heavy elements that are expected to sink out of visible layers on short timescales. The photometric {\sl TESS} data set spans $\sim 51$ days in total, and from this, we identified 56 significant pulsation frequencies, that include rotational frequency multiplets. In addition, we identified 30 combination frequencies in each sector. The oscillation frequencies that we found are associated with $g$-mode pulsations, with periods spanning from $\sim$ 260 s to $\sim$ 1400 s. We identified %three distinct rotational frequency triplets with a mean separation $δν_{\ell=1}$ of 4.67 $μ$Hz and a quintuplet with a mean separation $δν_{\ell=2}$ of 6.67 $μ$Hz, from which we estimated a rotation period of about $1.35 \pm 0.1$ days. We determined a constant period spacing of 41.20~s for $\ell= 1$ modes and 22.58\,s for $\ell= 2$ modes. We performed period-to-period fit analyses and found an asteroseismological model with $M_{\star}/M_{\odot}=0.632 \pm 0.03$, $T_{\rm eff}=11\, 635\pm 178$ K, and $\log{g}=8.048\pm0.005$ (with a hydrogen envelope mass of $M_{\rm H}\sim 5.6\times 10^{-5}M_{\star}$), in good agreement with the values derived from spectroscopy. We obtained an asteroseismic distance of 17.54 pc, which is in excellent agreement with that provided by {\sl Gaia} (17.51 pc). △ Less

Submitted 9 September, 2023; originally announced September 2023.

Comments: 17 pages, Accepted for publication in MNRAS

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

Roman Early-Definition Astrophysics Survey Opportunity: Galactic Roman Infrared Plane Survey (GRIPS)

Authors: Roberta Paladini, Catherine Zucker, Robert Benjamin, David Nataf, Dante Minniti, Gail Zasowski, Joshua Peek, Sean Carey, Lori Allen, Javier Alonso-Garcia, Joao Alves, Friederich Anders, Evangelie Athanassoula, Timothy C. Beers, Jonathan Bird, Joss Bland-Hwathorn, Anthony Brown, Sven Buder, Luca Casagrande, Andrew Casey, Santi Cassisi, Marcio Catelan, Ranga-Ram Chary, Andre-Nicolas Chene, David Ciardi , et al. (45 additional authors not shown)

Abstract: A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yiel… ▽ More A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yield an impressive dataset of ~120 billion sources and map the structure of our Galaxy. The effort would foster subsequent expansions in numerous dimensions (spatial, depth, wavelengths, epochs). Importantly, the survey would benefit from early defintion by the community, namely because the Galactic disk is a complex environment, and different science goals will require trade offs. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: Submitted to the Roman Project on October 22 2021 in response to a call for white papers on early-definition Astrophysics opportunity

A Unique Low-Mass Ratio Contact Eclipsing Binary System under the Period Cut-Off

Authors: A. Papageorgiou, P. -E. Christopoulou, E. Lalounta, C. E. Ferreira Lopes, M. Catelan, A. J. Drake, P. Xantzios, I. Alikakos

Abstract: We present a multi-band photometric analysis of CRTS J163819.6+03485, the first low mass ratio (LMR) contact binary system with a period under the contact binary (CB) period limit. The unprecedented combination of mass ratio and period makes this system unique for eclipsing binary (EB) research. Using new multi-band photometric observations, we explored the parameter space of this unique total EB… ▽ More We present a multi-band photometric analysis of CRTS J163819.6+03485, the first low mass ratio (LMR) contact binary system with a period under the contact binary (CB) period limit. The unprecedented combination of mass ratio and period makes this system unique for eclipsing binary (EB) research. Using new multi-band photometric observations, we explored the parameter space of this unique total EB system through a detailed scan in the mass ratio - inclination plane and using the PIKAIA genetic algorithm optimizer. The best set of relative physical parameters and corresponding uncertainties was adopted through Markov Chain Monte Carlo sampling of the parameter space. The resulting mass ratio of the system is $q = 0.16 \pm 0.01$. The absolute parameters were derived by adopting an empirical mass-luminosity relation. Period changes are also investigated by using new observations and archival photometric light curves from massive astronomical surveys, which revealed in a preliminary solution the presence of a possible low-mass tertiary companion. The origin and evolutionary status of the system are investigated through the detached-binary formation scenario. △ Less

Submitted 20 June, 2023; originally announced June 2023.

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)

Informative regularization for a multi-layer perceptron RR Lyrae classifier under data shift

Authors: Francisco Pérez-Galarce, Karim Pichara, Pablo Huijse, Márcio Catelan, Domingo Mery

Abstract: In recent decades, machine learning has provided valuable models and algorithms for processing and extracting knowledge from time-series surveys. Different classifiers have been proposed and performed to an excellent standard. Nevertheless, few papers have tackled the data shift problem in labeled training sets, which occurs when there is a mismatch between the data distribution in the training se… ▽ More In recent decades, machine learning has provided valuable models and algorithms for processing and extracting knowledge from time-series surveys. Different classifiers have been proposed and performed to an excellent standard. Nevertheless, few papers have tackled the data shift problem in labeled training sets, which occurs when there is a mismatch between the data distribution in the training set and the testing set. This drawback can damage the prediction performance in unseen data. Consequently, we propose a scalable and easily adaptable approach based on an informative regularization and an ad-hoc training procedure to mitigate the shift problem during the training of a multi-layer perceptron for RR Lyrae classification. We collect ranges for characteristic features to construct a symbolic representation of prior knowledge, which was used to model the informative regularizer component. Simultaneously, we design a two-step back-propagation algorithm to integrate this knowledge into the neural network, whereby one step is applied in each epoch to minimize classification error, while another is applied to ensure regularization. Our algorithm defines a subset of parameters (a mask) for each loss function. This approach handles the forgetting effect, which stems from a trade-off between these loss functions (learning from data versus learning expert knowledge) during training. Experiments were conducted using recently proposed shifted benchmark sets for RR Lyrae stars, outperforming baseline models by up to 3\% through a more reliable classifier. Our method provides a new path to incorporate knowledge from characteristic features into artificial neural networks to manage the underlying data shift problem. △ Less

Submitted 11 March, 2023; originally announced March 2023.

Journal ref: Pérez-Galarce, F., Pichara, K., Huijse, P., Catelan, M., & Mery, D. (2023). Informative regularization for a multi-layer perceptron RR Lyrae classifier under data shift. Astronomy and Computing, 100694

Uncertainties in the 12C+12C reaction rate and their impact on the composition of ultra-massive WDs

Authors: F. C. De Gerónimo, M. M. Miller Bertolami, M. Catelan, T. Battich

Abstract: Stars with initial masses 7 Msun . MZAMS . 9 Msun reach temperatures high enough to ignite C under degenerate conditions after the end of He-core burning (Garcia-Berro & Iben 1994). These isolated stars are expected to evolve into the so-called super AGB (SAGB) phase and may end their lives as ultra-massive ONe WDs (see Siess 2006,2007, 2010; Camisassa et al. 2019, and references therein). The exa… ▽ More Stars with initial masses 7 Msun . MZAMS . 9 Msun reach temperatures high enough to ignite C under degenerate conditions after the end of He-core burning (Garcia-Berro & Iben 1994). These isolated stars are expected to evolve into the so-called super AGB (SAGB) phase and may end their lives as ultra-massive ONe WDs (see Siess 2006,2007, 2010; Camisassa et al. 2019, and references therein). The exact proportions of O and Ne found in the core at the end of the SAGB phase will determine the cooling times and pulsational properties of these WDs. Uncertainties affecting the rates of nuclear reactions occurring during the C burning phase should have a measurable impact on the distribution of 16O, 20Ne, 23Na and 24Mg and, consequently, on the evolution of the WD. Here we present a study of the impact of uncertainties in the 12C(12C, α)20Ne and 12C(12C, p)23Na nuclear reaction rates (and their branching ratios) on the chemical structure of intermediate- to high-mass progenitors at the end of the C-burning phase. Using the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA) we computed evolutionary sequences for stars with initial masses 7.25<= MZAMS /Msun <=8.25, from the ZAMS to the SAGB phase, adopting different prescriptions for the 12C+12C burning rates. We found that adopting lower reaction rates for the 12C+12C burning delays C-ignition by at most 2700 yrs, and the ignition takes place in a position further from the center. Our results shows that differences in the 20Ne central abundances remain modest, below 14%. △ Less

Submitted 18 April, 2023; v1 submitted 28 February, 2023; originally announced March 2023.

Comments: Poster presented at the 22nd European Workshop on White Dwarfs. No proceedings were published at the conference

Stellar Variability in Ground-Based Photometric Surveys: An Overview

Authors: M. Catelan

Abstract: Time-resolved ground-based surveys in general, and photometric ones in particular, have played a crucial role in building up our knowledge of the properties, physical nature, and the very existence of the many different classes of variable stars and transient events that are currently known. Here I provide a brief overview of these developments, discussing some of the stumbling blocks that had to… ▽ More Time-resolved ground-based surveys in general, and photometric ones in particular, have played a crucial role in building up our knowledge of the properties, physical nature, and the very existence of the many different classes of variable stars and transient events that are currently known. Here I provide a brief overview of these developments, discussing some of the stumbling blocks that had to be overcome along the way, and others that may still hamper further progress in the area. A compilation of different types of past, present, and future surveys is also provided. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: Invited review, to appear in Mem. Soc. Astr. Italiana

Uranium Abundances and Ages of $R$-process Enhanced Stars with Novel U II Lines

Authors: Shivani P. Shah, Rana Ezzeddine, Alexander P. Ji, Terese Hansen, Ian U. Roederer, Márcio Catelan, Zoe Hackshaw, Erika M. Holmbeck, Timothy C. Beers, Rebecca Surman

Abstract: The ages of the oldest stars shed light on the birth, chemical enrichment, and chemical evolution of the Universe. Nucleocosmochronometry provides an avenue to determining the ages of these stars independent from stellar evolution models. The uranium abundance, which can be determined for metal-poor $r$-process enhanced (RPE) stars, has been known to constitute one of the most robust chronometers… ▽ More The ages of the oldest stars shed light on the birth, chemical enrichment, and chemical evolution of the Universe. Nucleocosmochronometry provides an avenue to determining the ages of these stars independent from stellar evolution models. The uranium abundance, which can be determined for metal-poor $r$-process enhanced (RPE) stars, has been known to constitute one of the most robust chronometers known. So far, U abundance determination has used a $single$ U II line at $\lambda3859$ Å. Consequently, U abundance has been reliably determined for only five RPE stars. Here, we present the first homogeneous U abundance analysis of four RPE stars using two novel U II lines at $\lambda4050$ Å and $\lambda4090$ Å, in addition to the canonical $\lambda3859$ Å line. We find that the U II lines at $\lambda4050$ Å and $\lambda4090$ Å are reliable and render U abundances in agreement with the $\lambda3859$ U abundance, for all the stars. We, thus, determine revised U abundances for RPE stars, 2MASS J09544277+5246414, RAVE J203843.2-002333, HE 1523-0901, and CS 31082-001, using multiple U II lines. We also provide nucleocosmochronometric ages of these stars based on the newly derived U, Th, and Eu abundances. The results of this study open up a new avenue to reliably and homogeneously determine U abundance for a significantly larger number of RPE stars. This will, in turn, enable robust constraints on the nucleocosmochronometric ages of RPE stars, which can be applied to understand the chemical enrichment and evolution in the early Universe, especially of $r$-process elements. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: Resubmitted to ApJ

Three Ultra-short Period Contact Eclipsing Binary Systems Mined from Massive Astronomical Surveys

Authors: A. Papageorgiou, P. -E. Christopoulou, C. E. Ferreira Lopes, E. Lalounta, M. Catelan, A. J. Drake

Abstract: We present the photometric analysis of three ultra-short period total eclipsing binaries in contact configuration, CRTS_J172718.0+431624, OGLE-BLG-ECL-000104, and OGLE-BLG-ECL-000012, mined from massive astronomical surveys. Using the available archival light curves (LCs) from Vista Variables in the Via Lactea (VVV), Optical Gravitational Lensing Experiment (OGLE), Zwicky Transient Facility (ZTF)… ▽ More We present the photometric analysis of three ultra-short period total eclipsing binaries in contact configuration, CRTS_J172718.0+431624, OGLE-BLG-ECL-000104, and OGLE-BLG-ECL-000012, mined from massive astronomical surveys. Using the available archival light curves (LCs) from Vista Variables in the Via Lactea (VVV), Optical Gravitational Lensing Experiment (OGLE), Zwicky Transient Facility (ZTF) and Catalina Sky Survey (CSS) in different passbands and new multi-band photometric observations with the 2.3 m Aristarchos telescope at Helmos Observatory, their relative physical parameters were derived. We explored the parameter space by using PIKAIA genetic algorithm optimizer. The best photometric solution and error budget estimation was adopted for each system through MCMC sampling of the global optimum. The approximate absolute parameters were derived for each contact system adopting an empirical mass-luminosity relation. All three systems have a mass ratio lower than 0.5. The exchange between primary and secondary depths of CRTS_J172718.0+431624 during the years 2016-2022 may be due to spot activity. In addition, we present a detailed analysis of the first well characterized shortest period contact eclipsing binary with total eclipses known so far (OGLE-BLG-ECL-000104). Thanks to VVV and OGLE LCs, new distances were derived for OGLE-BLG-ECL-000104 and OGLE-BLG-ECL-000012 using empirical period-luminosity relations. The origin and evolutionary status of all three ultra-short period contact binaries are thoroughly discussed in the context of the detached-binary formation channel. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Deep Drilling in the Time Domain with DECam: Survey Characterization

Authors: Melissa L. Graham, Robert A. Knop, Thomas Kennedy, Peter E. Nugent, Eric Bellm, Márcio Catelan, Avi Patel, Hayden Smotherman, Monika Soraisam, Steven Stetzler, Lauren N. Aldoroty, Autumn Awbrey, Karina Baeza-Villagra, Pedro H. Bernardinelli, Federica Bianco, Dillon Brout, Riley Clarke, William I. Clarkson, Thomas Collett, James R. A. Davenport, Shenming Fu, John E. Gizis, Ari Heinze, Lei Hu, Saurabh W. Jha , et al. (19 additional authors not shown)

Abstract: This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5… ▽ More This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5 to 6 images per night per filter in $g$, $r$, $i$, and/or $z$, have become publicly available (the proprietary period for this program is waived). We describe the real-time difference-image pipeline and how alerts are distributed to brokers via the same distribution system as the Zwicky Transient Facility (ZTF). In this paper, we focus on the two extragalactic deep fields (COSMOS and ELAIS-S1), characterizing the detected sources and demonstrating that the survey design is effective for probing the discovery space of faint and fast variable and transient sources. We describe and make publicly available 4413 calibrated light curves based on difference-image detection photometry of transients and variables in the extragalactic fields. We also present preliminary scientific analysis regarding Solar System small bodies, stellar flares and variables, Galactic anomaly detection, fast-rising transients and variables, supernovae, and active galactic nuclei. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 22 pages, 17 figures, 2 tables. Accepted to MNRAS

A segmented period-luminosity relation for nearby extragalactic $δ$ Scuti stars

Authors: C. E. Martínez-Vázquez, R. Salinas, A. K. Vivas, M. Catelan

Abstract: The period-luminosity relations (PLR) of Milky Way $δ$ Scuti ($δ$ Sct) stars have been described to the present day by a linear relation. However, when studying extragalactic systems such as the Magellanic Clouds and several dwarf galaxies, we notice for the first time a non-linear behaviour in the PLR of $δ$ Sct stars. Using the largest sample of $\sim 3700$ extragalactic $δ$ Sct stars from data… ▽ More The period-luminosity relations (PLR) of Milky Way $δ$ Scuti ($δ$ Sct) stars have been described to the present day by a linear relation. However, when studying extragalactic systems such as the Magellanic Clouds and several dwarf galaxies, we notice for the first time a non-linear behaviour in the PLR of $δ$ Sct stars. Using the largest sample of $\sim 3700$ extragalactic $δ$ Sct stars from data available in the literature $-$mainly based on OGLE and SuperMACHO survey in the Large Magellanic Cloud (LMC)$-$ we obtain that the best fit to the period-luminosity ($M_V$) plane is given by the following piecewise linear relation with a break at $\log{P} = -1.03 \pm 0.01$ (or $0.093 \pm 0.002$ d) for shorter periods (sp) and longer periods (lp) than the break-point: $$M_V^{sp} = -7.08 (\pm 0.25) \log{P} -5.74 (\pm 0.29) ;\hspace{5pt} \log{P} < -1.03$$ $$M_V^{lp} = M_V^{sp} + 4.38 (\pm 0.32) \cdot (\log{P} + 1.03 (\pm 0.01));\hspace{5pt} \log{P} \geq -1.03$$ Geometric or depth effects in the LMC, metallicity dependence, or different pulsation modes are discarded as possible causes of this segmented PLR seen in extragalactic $δ$ Sct stars. The origin of the segmented relation at $\sim 0.09$ days remains unexplained based on the current data. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: 9 pages, 3 figures, 1 table. Accepted for publication into The Astrophysical Journal Letters

Period Change Rates of Large Magellanic Cloud Cepheids using MESA

Authors: F. Espinoza-Arancibia, M. Catelan, G. Hajdu, N. Rodríguez-Segovia, G. Boggiano, K. Joachimi, C. Muñoz-López, C. Ordenes-Huanca, C. Orquera-Rojas, P. Torres, Á. Valenzuela-Navarro

Abstract: Pulsating stars, such as Cepheids and RR Lyrae, offer us a window to measure and study changes due to stellar evolution. In this work, we study the former by calculating a set of evolutionary tracks of stars with an initial mass of 4 to 7 $M_\odot$, varying the initial rotation rate and metallicity, using the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA). Using Radi… ▽ More Pulsating stars, such as Cepheids and RR Lyrae, offer us a window to measure and study changes due to stellar evolution. In this work, we study the former by calculating a set of evolutionary tracks of stars with an initial mass of 4 to 7 $M_\odot$, varying the initial rotation rate and metallicity, using the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA). Using Radial Stellar Pulsations (RSP), a recently added functionality of MESA, we obtained theoretical instability strip (IS) edges and linear periods for the radial fundamental mode. Period-age, period-age-temperature, period-luminosity, and period-luminosity-temperature relationships were derived for three rotation rates and metallicities, showing a dependence on crossing number, position in the IS, rotation, and metallicity. We calculated period change rates (PCRs) based on the linear periods from RSP. We compared our models with literature results using the Geneva code, and found large differences, as expected due to the different implementations of rotation between codes. In addition, we compared our theoretical PCRs with those measured in our recent work for Large Magellanic Cloud Cepheids. We found good overall agreement, even though our models do not reach the short-period regime exhibited by the empirical data. Implementations of physical processes not yet included in our models, such as pulsation-driven mass loss, an improved treatment of convection that may lead to a better description of the instability strip edges, as well as consideration of a wider initial mass range, could all help improve the agreement with the observed PCRs. △ Less

Submitted 12 October, 2023; v1 submitted 21 September, 2022; originally announced September 2022.

Comments: 19 pages, 17 figures. Accepted by MNRAS

Semi-Supervised Classification and Clustering Analysis for Variable Stars

Authors: R. Pantoja, M. Catelan, K. Pichara, P. Protopapas

Abstract: The immense amount of time series data produced by astronomical surveys has called for the use of machine learning algorithms to discover and classify several million celestial sources. In the case of variable stars, supervised learning approaches have become commonplace. However, this needs a considerable collection of expert-labeled light curves to achieve adequate performance, which is costly t… ▽ More The immense amount of time series data produced by astronomical surveys has called for the use of machine learning algorithms to discover and classify several million celestial sources. In the case of variable stars, supervised learning approaches have become commonplace. However, this needs a considerable collection of expert-labeled light curves to achieve adequate performance, which is costly to construct. To solve this problem, we introduce two approaches. First, a semi-supervised hierarchical method, which requires substantially less trained data than supervised methods. Second, a clustering analysis procedure that finds groups that may correspond to classes or sub-classes of variable stars. Both methods are primarily supported by dimensionality reduction of the data for visualization and to avoid the curse of dimensionality. We tested our methods with catalogs collected from OGLE, CSS, and Gaia surveys. The semi-supervised method reaches a performance of around 90\% for all of our three selected catalogs of variable stars using only $5\%$ of the data in the training. This method is suitable for classifying the main classes of variable stars when there is only a small amount of training data. Our clustering analysis confirms that most of the clusters found have a purity over 90\% with respect to classes and 80\% with respect to sub-classes, suggesting that this type of analysis can be used in large-scale variability surveys as an initial step to identify which classes or sub-classes of variable stars are present in the data and/or to build training sets, among many other possible applications. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 23 pages, 21 figures, 4 tables, submitted to MNRAS

Near-infrared observations of RR Lyrae and Type II Cepheid variables in the metal-rich bulge globular cluster NGC 6441

Authors: A. Bhardwaj, S. M. Kanbur, M. Rejkuba, M. Marconi, M. Catelan, V. Ripepi, H. P. Singh

Abstract: NGC 6441 is a bulge globular cluster with an unusual horizontal branch morphology and a rich population of RR Lyrae (RRL) and Type II Cepheid (T2C) variables that is unexpected for its relatively high metallicity. We present near-infrared (NIR, $JHK_s$) time-series observations of 42 RRL, 8 T2Cs, and 10 eclipsing binary candidate variables in NGC 6441. The multi-epoch observations were obtained us… ▽ More NGC 6441 is a bulge globular cluster with an unusual horizontal branch morphology and a rich population of RR Lyrae (RRL) and Type II Cepheid (T2C) variables that is unexpected for its relatively high metallicity. We present near-infrared (NIR, $JHK_s$) time-series observations of 42 RRL, 8 T2Cs, and 10 eclipsing binary candidate variables in NGC 6441. The multi-epoch observations were obtained using the FLAMINGOS-2 instrument on the 8-m Gemini South telescope. Multi-band data are used to investigate pulsation properties of RRL and T2C variables including their light curves, instability strip, period--amplitude diagrams, and period--luminosity and period--wesenheit relations (PLRs and PWRs) in $JHK_s$ filters. NIR pulsation properties of RRL variables are well fitted with theoretical models that have canonical helium content and mean-metallicity of NGC 6441 ([Fe/H]$=-0.44\pm0.07$~dex). The helium-enhanced RRL models predict brighter NIR magnitudes and bluer colors than the observations of RRL in the cluster. This suggests that these RRL variables in NGC 6441 are either not significantly helium enhanced as previously thought or the impact of such enhancement is smaller in NIR than the predictions of the pulsation models. We also use theoretical calibrations of RRL period--luminosity--metallicity (PLZ) relations to simultaneously estimate the mean reddening, $E(J-K_s)=0.26\pm0.06$~mag, and the distance, $d=12.67\pm0.09$~kpc, to NGC 6441. Our mean reddening value and the distance are consistent with the independent estimates from the bulge reddening map based on red clump stars and the latest Gaia data, respectively. Our distance and reddening values provide a very good agreement between the PLRs of T2Cs in NGC 6441 and those for RRL and T2Cs in Galactic globular clusters that span a broad range of metallicity. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: 17 pages, 14 figures, accepted in A&A, abridged abstract

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

Large amplitude periodic outbursts and long period variables in the VVV VIRAC2-$β$ database

Authors: Zhen Guo, P. W. Lucas, L. C. Smith, C. Clarke, C. Contreras Peña, A. Bayo, C. Briceño, J. Elias, R. G. Kurtev, J. Borissova, J. Alonso-García, D. Minniti, M. Catelan, F. Nikzat, C. Morris, N. Miller

Abstract: The VISTA Variables in the Via Lactea (VVV) survey obtained near-infrared photometry toward the Galactic bulge and the southern disc plane for a decade (2010 - 2019). We designed a modified Lomb-Scargle method to search for large-amplitude ($Δ$Ks > 1.5 mag) mid to long-term periodic variables (P > 10 d) in the 2nd version of VVV Infrared Astrometric Catalogue (VIRAC2-$β$). In total, 1520 periodic… ▽ More The VISTA Variables in the Via Lactea (VVV) survey obtained near-infrared photometry toward the Galactic bulge and the southern disc plane for a decade (2010 - 2019). We designed a modified Lomb-Scargle method to search for large-amplitude ($Δ$Ks > 1.5 mag) mid to long-term periodic variables (P > 10 d) in the 2nd version of VVV Infrared Astrometric Catalogue (VIRAC2-$β$). In total, 1520 periodic sources were discovered, including 59 candidate periodic outbursting young stellar objects (YSOs), based on the unique morphology of the phase-folded light curves, proximity to Galactic HII regions and mid-infrared colours. Five sources are spectroscopically confirmed as accreting YSOs. Both fast-rise/slow-decay and slow-rise/fast-decay periodic outbursts were found, but fast-rise/slow-decay outbursts predominate at the highest amplitudes. The multi-wavelength colour variations are consistent with a variable mass accretion process, as opposed to variable extinction. The cycles are likely to be caused by dynamical perturbations from stellar or planetary companions within the circumstellar disc. An additional search for periodic variability amongst YSO candidates in published Spitzer-based catalogues yielded a further 71 candidate periodic accretors, mostly with lower amplitudes. These resemble cases of pulsed accretion but with unusually long periods and greater regularity. The majority of other long-period variables are pulsating dusty Miras with smooth and symmetric light curves. We find that some Miras have redder $W3 - W4$ colours than previously thought, most likely due to their surface chemical compositions. △ Less

Submitted 21 March, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 29 pages, 25 figures

The VVV survey: Long-period variable stars I. Photometric catalog of ten VVV/OGLE tiles

Authors: F. Nikzat, C. E. Ferreira Lopes, M. Catelan, R. Contreras Ramos, M. Zoccali, A. Rojas-Arriagada, V. F. Braga, D. Minniti, J. Borissova, I. Becker

Abstract: Long-period variable stars (LPVs) are pulsating red giants, primarily in the asymptotic giant branch phase, and they include both Miras and semi-regular variables (SRVs). Their period-age and period-luminosity relations enable us to trace different stellar populations, as they are intrinsically very bright and cover a wide range in distances and ages. The purpose of this study is to establish a ce… ▽ More Long-period variable stars (LPVs) are pulsating red giants, primarily in the asymptotic giant branch phase, and they include both Miras and semi-regular variables (SRVs). Their period-age and period-luminosity relations enable us to trace different stellar populations, as they are intrinsically very bright and cover a wide range in distances and ages. The purpose of this study is to establish a census of LPV stars in a region close to the Galactic center, using the six-year database of the Vista Variables in the Vía Láctea (VVV) ESO Public Survey, as well as to describe the methodology that was employed to search for and characterize LPVs using VVV data. Near-IR surveys such as VVV provide a unique opportunity to probe the high-extinction innermost regions of the Milky Way. The detection and analysis of the intrinsically bright Miras in this region could provide us with an excellent probe of the properties of the Milky Way far behind its bulge. We used point-spread function photometry for all available $K_{s}$-band images in ten VVV tiles, covering $16.4~°^2$ in total, overlapping fields observed in the course of the Optical Gravitational Lensing Experiment (OGLE)-III survey. We designed a method to select LPV candidates, and we used the known variables from OGLE-III and other known variables from the literature to test our approach. The reduced $χ^2$ statistic, along with the flux-independent index $K_{(fi)}$, were used in our analysis. The Lomb-Scargle period search method, Fourier analysis, template fitting, and visual inspection were then performed to refine our sample and characterize the properties of the stars included in our catalog. A final sample of 130 Mira candidates, of which 129 are new discoveries, was thus obtained, with periods in the range between about 80 and 1400~days. △ Less

Submitted 11 March, 2022; originally announced March 2022.

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

New low mass ratio contact binaries in the Catalina Sky Survey

Authors: P. E. Christopoulou, Eleni Lalounta, Athanasios Papageorgiou, C. E. Ferreira Lopes, Marcio Catelan, Andrew J. Drake

Abstract: We present the identification and photometric analysis of 30 new low mass ratio (LMR) totally eclipsing contact binaries found in Catalina Sky Survey data. The LMR candidates are identified using Fourier coefficients and visual inspection. We perform a detailed scan in the parameter plane of mass-ratio (q) versus inclination (i) using Phoebe-0.31 scripter to derive the best (q,i) pair for the init… ▽ More We present the identification and photometric analysis of 30 new low mass ratio (LMR) totally eclipsing contact binaries found in Catalina Sky Survey data. The LMR candidates are identified using Fourier coefficients and visual inspection. We perform a detailed scan in the parameter plane of mass-ratio (q) versus inclination (i) using Phoebe-0.31 scripter to derive the best (q,i) pair for the initial models. The relative physical parameters are determined from the final model of each system. A Monte-Carlo approach was adopted to derive the parameter errors. The resulting parameters confirm the identification. The approximate absolute physical parameters of the systems are estimated based on the light curve solutions and Gaia early Data Release 3 distances. Twelve out of 30 new systems have fill-out factors $f\gt 50 \%$ and $q \leq 0.25$ (deep contact LMR systems), and 8 of them, to within errors, are extreme LMR deep systems with $q\leq0.1$. We discuss the evolutionary status of the 30 LMR systems in comparison with the most updated catalog of LMR systems from the literature. The scenario of the LMR systems as pre-merger candidates forming fast rotating stars is investigated for all systems, new and old, based both on Hut's stability criteria and critical instability mass ratio ($q_{inst}$) relation. CSS$\_$J075848.2+125656, with $\frac{q}{q_{inst}}=1.23\pm 0.23$, and CSS$\_$J093010.1-021624, with $\frac{q}{q_{inst}}=1.25\pm 0.23$, can be considered as merger candidates. △ Less

Submitted 25 February, 2022; originally announced February 2022.

The composition of massive white dwarfs and their dependence on the C-burning modeling

Authors: Francisco C. De Gerónimo, Marcelo M. Miller Bertolami, Francisco Plaza, Márcio Catelan

Abstract: Recent computations of the interior composition of ultra-massive white dwarfs (WD) have suggested that some white dwarfs could be composed of neon (Ne)-dominated cores. This result is at variance with our previous understanding of the chemical structure of massive white dwarfs, where oxygen is the predominant element. In addition, it is not clear whether some hybrid carbon (C) oxygen (O)-Ne white… ▽ More Recent computations of the interior composition of ultra-massive white dwarfs (WD) have suggested that some white dwarfs could be composed of neon (Ne)-dominated cores. This result is at variance with our previous understanding of the chemical structure of massive white dwarfs, where oxygen is the predominant element. In addition, it is not clear whether some hybrid carbon (C) oxygen (O)-Ne white dwarfs might form when convective boundary mixing is accounted for during the propagation of the C-flame, in the C-burning stage. Both Ne-dominated and hybrid CO-Ne core would have measurable consequences for asteroseismological studies based on evolutionary models. In this work we explore in detail to which extent differences in the adopted micro- and macro-physics can explain the different final white dwarf compositions that have been found by different authors. Additionally, we explored the impact of such differences in the cooling times, crystallization and pulsational properties of pulsating WDs. We explore the impact of the intensity of convective boundary mixing during the C-flash, extreme mass-loss rates, and the size of the adopted nuclear networks on the final composition, age, crystallization and pulsational properties of white dwarfs. Based on the insight coming from 3D hydro-dynamical simulations, we expect that the very slow propagation of the carbon flame will be altered by turbulent entrainment affecting the inward propagation of the flame. Also, we find that Ne-dominated chemical profiles of massive WDs recently reported appear in their modeling due to the overlooking of a key nuclear reaction. We find that the inaccuracies in the chemical composition of ultra-massive white dwarfs recently reported lead to differences of 10% in the cooling times and degree of crystallization and about 8% in the period spacing of the models once they reach the ZZ Ceti instability strip. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: 10 pages, 9 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 659, A150 (2022)

Periodic Variable Stars Modulated by Time-Varying Parameters

Authors: Giovanni Motta, Darlin Soto, Márcio Catelan

Abstract: Many astrophysical phenomena are time-varying, in the sense that their brightness change over time. In the case of periodic stars, previous approaches assumed that changes in period, amplitude, and phase are well described by either parametric or piecewise-constant functions. With this paper, we introduce a new mathematical model for the description of the so-called modulated light curves, as foun… ▽ More Many astrophysical phenomena are time-varying, in the sense that their brightness change over time. In the case of periodic stars, previous approaches assumed that changes in period, amplitude, and phase are well described by either parametric or piecewise-constant functions. With this paper, we introduce a new mathematical model for the description of the so-called modulated light curves, as found in periodic variable stars that exhibit smoothly time-varying parameters such as amplitude, frequency, and/or phase. Our model accounts for a smoothly time-varying trend, and a harmonic sum with smoothly time-varying weights. In this sense, our approach is flexible because it avoids restrictive assumptions (parametric or piecewise-constant) about the functional form of trend and amplitudes. We apply our methodology to the light curve of a pulsating RR Lyrae star characterised by the Blazhko effect. To estimate the time-varying parameters of our model, we develop a semi-parametric method for unequally spaced time series. The estimation of our time-varying curves translates into the estimation of time-invariant parameters that can be performed by ordinary least-squares, with the following two advantages: modeling and forecasting can be implemented in a parametric fashion, and we are able to cope with missing observations. To detect serial correlation in the residuals of our fitted model, we derive the mathematical definition of the spectral density for unequally spaced time series. The proposed method is designed to estimate smoothly time-varying trend and amplitudes, as well as the spectral density function of the errors. We provide simulation results and applications to real data. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 26 pages, 6 figures, to be published in The Astrophysical Journal

Period Change Rates in Large Magellanic Cloud Cepheids Revisited

Authors: N. Rodríguez-Segovia, G. Hajdu, M. Catelan, F. Espinoza-Arancibia, G. Boggiano, C. Cenzano, E. Garcés H., K. Joachimi, C. Muñoz-López, C. Ordenes-Huanca, C. Orquera-Rojas, P. Torres, Á. Valenzuela-Navarro

Abstract: The period-change rate (PCR) of pulsating variable stars is a useful probe of changes in their interior structure, and thus of their evolutionary stages. So far, the PCRs of Classical Cepheids in the Large Magellanic Cloud (LMC) have been explored in a limited sample of the total population of these variables. Here we use a template-based method to build observed minus computed (O-C) period diagra… ▽ More The period-change rate (PCR) of pulsating variable stars is a useful probe of changes in their interior structure, and thus of their evolutionary stages. So far, the PCRs of Classical Cepheids in the Large Magellanic Cloud (LMC) have been explored in a limited sample of the total population of these variables. Here we use a template-based method to build observed minus computed (O-C) period diagrams, from which we can derive PCRs for these stars by taking advantage of the long time baseline afforded by the Digital Access to a Sky Century @ Harvard (DASCH) light curves, combined with additional data from the Optical Gravitational Lensing Experiment (OGLE), the MAssive Compact Halo Object (MACHO) project, Gaia's Data Release 2, and in some cases the All-Sky Automated Survey (ASAS). From an initial sample of 2315 sources, our method provides an unprecedented sample of 1303 LMC Classical Cepheids with accurate PCRs, the largest for any single galaxy, including the Milky Way. The derived PCRs are largely compatible with theoretically expected values, as computed by our team using the Modules for Experiments in Stellar Astrophysics (MESA) code, as well as with similar previous computations available in the literature. Additionally, five long-period (P>50 d) sources display a cyclic behavior in their O-C diagrams, which is clearly incompatible with evolutionary changes. Finally, on the basis of their large positive PCR values, two first-crossing Cepheid candidates are identified. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 11 pages, 6 figures. Accepted by MNRAS

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: a pioneering process of community-focused experimental design

Authors: Federica B. Bianco, Željko Ivezić, R. Lynne Jones, Melissa L. Graham, Phil Marshall, Abhijit Saha, Michael A. Strauss, Peter Yoachim, Tiago Ribeiro, Timo Anguita, Franz E. Bauer, Eric C. Bellm, Robert D. Blum, William N. Brandt, Sarah Brough, Màrcio Catelan, William I. Clarkson, Andrew J. Connolly, Eric Gawiser, John Gizis, Renee Hlozek, Sugata Kaviraj, Charles T. Liu, Michelle Lochner, Ashish A. Mahabal , et al. (21 additional authors not shown)

Abstract: Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core scienc… ▽ More Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue. △ Less

Submitted 1 September, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: Submitted as the opening paper of the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

Using classical Cepheids to study the far side of the Milky Way disk. II. The spiral structure in the first and fourth Galactic quadrants

Authors: J. H. Minniti, M. Zoccali, A. Rojas-Arriagada, D. Minniti, L. Sbordone, R. Contreras Ramos, V. F. Braga, M. Catelan, S. Duffau, W. Gieren, M. Marconi, A. A. R. Valcarce

Abstract: In an effort to improve our understanding of the spiral arm structure of the Milky Way, we use Classical Cepheids (CCs) to increase the number of young tracers on the far side of the Galactic disk with accurately determined distances. We use a sample of 30 CCs, discovered using near-infrared photometry from the VISTA Variables in the Vía Láctea survey (VVV) and classified based on their radial vel… ▽ More In an effort to improve our understanding of the spiral arm structure of the Milky Way, we use Classical Cepheids (CCs) to increase the number of young tracers on the far side of the Galactic disk with accurately determined distances. We use a sample of 30 CCs, discovered using near-infrared photometry from the VISTA Variables in the Vía Láctea survey (VVV) and classified based on their radial velocities and metallicities. We combine them with another 20 CCs from the literature for which VVV photometry is available. The compiled sample of CCs with homogeneously computed distances based on VVV infrared photometry was employed as a proof of concept to trace the spiral structure in the poorly explored far side of the disk. Although the use of CCs has some caveats, these variables are currently the only available young tracers in the far side disk for which a numerous sample with accurate distances can be obtained. Therefore, a larger sample could allow us to make a significant step forward in our understanding of the Milky Way disk as a whole. We present preliminary evidence that CCs favor: a spiral arm model with two main arms (Perseus and Scutum-Centaurus) branching out into four arms at galactocentric distances, $R_\mathrm {GC}\gtrsim5-6\,\mathrm{kpc}$; the extension of the Scutum-Centaurus arm behind the Galactic center; a possible connection between the Perseus arm and the Norma tangency direction. The current sample of CCs in the far side of the Galaxy are in the mid-plane, arguing against the presence of a severely warped disk at small Galactocentric distances ($R_\mathrm {GC}\lesssim12\,\mathrm{kpc}$) in the studied area. The discovery and characterization of CCs at near-IR wavelengths appears to be a promising tool to complement studies based on other spiral arm tracers and extend them to the far side of our Galaxy. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: 10 pages, 6 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 654, A138 (2021)

CAPOS: The bulge Cluster APOgee Survey I. Overview and initial ASPCAP results

Authors: Doug Geisler, Sandro Villanova, Julia E. O'Connell, Roger E. Cohen, Christian Moni Bidin, José G. Fernández-Trincado, Cesar Muñoz, Dante Minniti, Manuela Zoccali, Alvaro Rojas-Arriagada, Rodrigo Contreras Ramos, Márcio Catelan, Francesco Mauro, Cristían Cortés, C. E. Ferreira Lopes, Anke Arentsen, Else Starkenburg, Nicolas F. Martin, Baitian Tang, Celeste Parisi, Javier Alonso-García, Felipe Gran, Katia Cunha, Verne Smith, Steven R. Majewski , et al. (17 additional authors not shown)

Abstract: Context. Bulge globular clusters (BGCs) are exceptional tracers of the formation and chemodynamical evolution of this oldest Galactic component. However, until now, observational difficulties have prevented us from taking full advantage of these powerful Galactic archeological tools. Aims. CAPOS, the bulge Cluster APOgee Survey, addresses this key topic by observing a large number of BGCs, most of… ▽ More Context. Bulge globular clusters (BGCs) are exceptional tracers of the formation and chemodynamical evolution of this oldest Galactic component. However, until now, observational difficulties have prevented us from taking full advantage of these powerful Galactic archeological tools. Aims. CAPOS, the bulge Cluster APOgee Survey, addresses this key topic by observing a large number of BGCs, most of which have only been poorly studied previously. Even their most basic parameters, such as metallicity, [α/Fe], and radial velocity, are generally very uncertain. We aim to obtain accurate mean values for these parameters, as well as abundances for a number of other elements, and explore multiple populations. In this first paper, we describe the CAPOS project and present initial results for seven BGCs. Methods. CAPOS uses the APOGEE-2S spectrograph observing in the H band to penetrate obscuring dust toward the bulge. For this initial paper, we use abundances derived from ASPCAP, the APOGEE pipeline. Results. We derive mean [Fe/H] values of $-$0.85$\pm$0.04 (Terzan 2), $-$1.40$\pm$0.05 (Terzan 4), $-$1.20$\pm$0.10 (HP 1), $-$1.40$\pm$0.07 (Terzan 9), $-$1.07$\pm$0.09 (Djorg 2), $-$1.06$\pm$0.06 (NGC 6540), and $-$1.11$\pm$0.04 (NGC 6642) from three to ten stars per cluster. We determine mean abundances for eleven other elements plus the mean [$α$/Fe] and radial velocity. CAPOS clusters significantly increase the sample of well-studied Main Bulge globular clusters (GCs) and also extend them to lower metallicity. We reinforce the finding that Main Bulge and Main Disk GCs, formed in situ, have [Si/Fe] abundances slightly higher than their accreted counterparts at the same metallicity. We investigate multiple populations and find our clusters generally follow the light-element (anti)correlation trends of previous studies of GCs of similar metallicity. We finally explore the abundances ... △ Less

Submitted 31 May, 2021; originally announced June 2021.

Comments: 27 pages, 10 figures, 5 tables. Abridged abstract. Accepted for publication in Astronomy & Astrophysics (A&A)

Informative Bayesian model selection for RR Lyrae star classifiers

Authors: F. Pérez-Galarce, K. Pichara, P. Huijse, M. Catelan, D. Mery

Abstract: Machine learning has achieved an important role in the automatic classification of variable stars, and several classifiers have been proposed over the last decade. These classifiers have achieved impressive performance in several astronomical catalogues. However, some scientific articles have also shown that the training data therein contain multiple sources of bias. Hence, the performance of thos… ▽ More Machine learning has achieved an important role in the automatic classification of variable stars, and several classifiers have been proposed over the last decade. These classifiers have achieved impressive performance in several astronomical catalogues. However, some scientific articles have also shown that the training data therein contain multiple sources of bias. Hence, the performance of those classifiers on objects not belonging to the training data is uncertain, potentially resulting in the selection of incorrect models. Besides, it gives rise to the deployment of misleading classifiers. An example of the latter is the creation of open-source labelled catalogues with biased predictions. In this paper, we develop a method based on an informative marginal likelihood to evaluate variable star classifiers. We collect deterministic rules that are based on physical descriptors of RR Lyrae stars, and then, to mitigate the biases, we introduce those rules into the marginal likelihood estimation. We perform experiments with a set of Bayesian Logistic Regressions, which are trained to classify RR Lyraes, and we found that our method outperforms traditional non-informative cross-validation strategies, even when penalized models are assessed. Our methodology provides a more rigorous alternative to assess machine learning models using astronomical knowledge. From this approach, applications to other classes of variable stars and algorithmic improvements can be developed. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Variable stars in the VVV globular clusters. II. NGC6441, NGC6569, NGC6626 (M28), NGC6656 (M22), 2MASS-GC02, and Terzan10

Authors: Javier Alonso-García, Leigh C. Smith, Márcio Catelan, Dante Minniti, Camila Navarrete, Jura Borissova, Julio A. Carballo-Bello, Rodrigo Contreras Ramos, José G. Fernández-Trincado, Carlos E. Ferreira Lopes, Felipe Gran, Elisa R. Garro, Doug Geisler, Zhen Guo, Maren Hempel, Eamonn Kerins, Philip W. Lucas, Tali Palma, Karla Peña Ramírez, Sebastián Ramírez Alegría, Roberto K. Saito

Abstract: The Galactic globular clusters (GGCs) located in the inner regions of the Milky Way suffer from high extinction that makes their observation challenging. The VVV survey provides a way to explore these GGCs in the near-infrared where extinction effects are highly diminished. We conduct a search for variable stars in several inner GGCs, taking advantage of the unique multi-epoch, wide-field, near-in… ▽ More The Galactic globular clusters (GGCs) located in the inner regions of the Milky Way suffer from high extinction that makes their observation challenging. The VVV survey provides a way to explore these GGCs in the near-infrared where extinction effects are highly diminished. We conduct a search for variable stars in several inner GGCs, taking advantage of the unique multi-epoch, wide-field, near-infrared photometry provided by the VVV survey. We are especially interested in detecting classical pulsators that will help us constrain the physical parameters of these GGCs. In this paper, the second of a series, we focus on NGC6656 (M22), NGC6626 (M28), NGC6569, and NGC6441; these four massive GGCs have known variable sources, but quite different metallicities. We also revisit 2MASS-GC02 and Terzan10, the two GGCs studied in the first paper of this series. We present an improved method and a new parameter that efficiently identify variable candidates in the GGCs. We also use the proper motions of those detected variable candidates and their positions in the sky and in the color-magnitude diagrams to assign membership to the GGCs. We identify and parametrize in the near-infrared numerous variable sources in the studied GGCs, cataloging tens of previously undetected variable stars. We recover many known classical pulsators in these clusters, including the vast majority of their fundamental mode RR Lyrae. We use these pulsators to obtain distances and extinctions toward these objects. Recalibrated period-luminosity-metallicity relations for the RR Lyrae bring the distances to these GGCs to a closer agreement with those reported by Gaia, except for NGC6441. Recovered proper motions for these GGCs also agree with those reported by Gaia, except for 2MASS-GC02, the most reddened GGC in our sample, where the VVV near-infrared measurements provide a more accurate determination of its proper motions. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: Accepted for publication in A&A, 18 pages, 14 Figures, 9 Tables

Journal ref: A&A 651, A47 (2021)

Studies of RR Lyrae Variables in Binary Systems. I.: Evidence of a Trimodal Companion Mass Distribution

Authors: Gergely Hajdu, Grzegorz Pietrzyński, Johanna Jurcsik, Márcio Catelan, Paulina Karczmarek, Bogumił Pilecki, Igor Soszyński, Andrzej Udalski, Ian B. Thompson

Abstract: We present 87 candidates for RR Lyrae variable stars in binary systems, based on our new search using the light-travel time effect (LTTE) and observed - calculated ($O-C$) diagrams in the Galactic bulge time-series photometry of the Optical Gravitational Lensing Experiment. Out of these, 61 are new candidates, while 26 have been announced previously. Furthermore, 12 stars considered as binary cand… ▽ More We present 87 candidates for RR Lyrae variable stars in binary systems, based on our new search using the light-travel time effect (LTTE) and observed - calculated ($O-C$) diagrams in the Galactic bulge time-series photometry of the Optical Gravitational Lensing Experiment. Out of these, 61 are new candidates, while 26 have been announced previously. Furthermore, 12 stars considered as binary candidates in earlier works are discarded from the list, either as they were found to have $O-C$ diagrams incompatible with the LTTE or because their long-term periodicity is definitely caused by the Blazhko effect. This sample of RR Lyrae binary candidates allows us to draw the first firm conclusions about the population of such objects: no candidate has an orbital period below 1000 days, while their occurrence rate steadily increases with increasing period, and peaks between 3000 and 4000 days; however, the decrease in the number of stars toward even longer periods is probably the result of observational biases. The eccentricities show a very significant concentration between 0.25 and 0.3, with a quarter of candidates found in this single bin, overlaid on an otherwise flat distribution between 0.05 and 0.6. Only six stars have higher inferred eccentricities above 0.6. Lastly, the distribution of the mass functions is highly peculiar, exhibiting strong trimodality. We interpret these modes as the presence of three distinct groups of companions, with typical inferred masses of $\sim0.6$, $\sim0.2$, and $\sim0.067\,\mathrm{M}_\odot$, which can be associated with populations of white dwarf and main sequence, red dwarf, and brown dwarf companions, respectively. △ Less

Submitted 16 June, 2021; v1 submitted 8 May, 2021; originally announced May 2021.

Comments: Accepted for publication in the ApJ; 20 pages, 3 figures, 4 tables

VVV survey near-infrared colour catalogue of known variable stars

Authors: Fábio R. Herpich, Carlos E. Ferreira Lopes, Roberto K. Saito, Dante Minniti, Alessandro Ederoclite, Thiago S. Ferreira, Marcio Catelan

Abstract: Context. The Vista Variables in the Via Lactea (VVV) near-infrared variability survey explores some of the most complex regions of the Milky Way bulge and disk in terms of high extinction and high crowding. Aims. We add a new wavelength dimension to the optical information available at the American Association of Variable Star Observers International Variable Star Index (VSX-AAVSO) catalogue to te… ▽ More Context. The Vista Variables in the Via Lactea (VVV) near-infrared variability survey explores some of the most complex regions of the Milky Way bulge and disk in terms of high extinction and high crowding. Aims. We add a new wavelength dimension to the optical information available at the American Association of Variable Star Observers International Variable Star Index (VSX-AAVSO) catalogue to test the VVV survey near-infrared photometry to better characterise these objects. Methods. We cross-matched the VVV and the VSX-AAVSO catalogues along with Gaia Data Release 2 photometry and parallax. Results. We present a catalogue that includes accurate individual coordinates, near-infrared magnitudes (ZY JHKs), extinctions Aks, and distances based on Gaia parallaxes. We also show the near-infrared CMDs and spatial distributions for the different VSX types of variable stars, including important distance indicators, such as RR Lyrae, Cepheids, and Miras. By analysing the photometric flags in our catalogue, we found that about 20% of the stars with measured and verified variability are flagged as non-stellar sources, even when they are outside of the saturation and/or noise regimes. Additionally, we pair-matched our sample with the VIVA catalogue and found that more than half of our sources are missing from the VVV variability list, mostly due to observations with low signal-to-noise ratio or photometric problems with a low percentage due to failures in the selection process. Conclusions. Our results suggest that the current knowledge of the variability in the Galaxy is biased to nearby stars with low extinction. The present catalogue also provides the groundwork for characterising the results of future large variability surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time in the highly crowded and reddened regions of the Galactic plane, as well as follow-up campaigns for △ Less

Submitted 5 February, 2021; originally announced February 2021.

Comments: 14 pages, 9 figures, A&A accepted

Journal ref: A&A 647, A169 (2021)

The elusive tidal tails of the Milky Way globular cluster NGC 7099

Authors: A. E. Piatti, J. A. Carballo-Bello, M. D. Mora, C. Cenzano, C. Navarrete, M. Catelan

Abstract: We present results on the extra-tidal features of the Milky Way globular cluster NGC 7099, using deep gr photometry obtained with the Dark Energy Camera (DECam). We reached nearly 6 mag below the cluster Main Sequence (MS) turnoff, so that we dealt with the most suitable candidates to trace any stellar structure located beyond the cluster tidal radius. From star-by-star reddening corrected color-m… ▽ More We present results on the extra-tidal features of the Milky Way globular cluster NGC 7099, using deep gr photometry obtained with the Dark Energy Camera (DECam). We reached nearly 6 mag below the cluster Main Sequence (MS) turnoff, so that we dealt with the most suitable candidates to trace any stellar structure located beyond the cluster tidal radius. From star-by-star reddening corrected color-magnitude diagrams (CMDs) we defined four adjacent strips along the MS, for which we built the respective stellar density maps, once the contamination by field stars was properly removed. The resulting field star cleaned stellar density maps show a short tidal tail and some scattered debris. Such extra-tidal features are hardly detected when much shallower Gaia DR2 data sets are used and the same CMD field star cleaning procedure is applied. Indeed, by using 2.5 magnitudes below the cluster MS turnoff as the faintest limit (G < 20.5 mag), cluster members turned out to be distributed within the cluster's tidal radius, and some hints for field star density variations are found across a circle of radius 3.5deg centered on the cluster and with similar CMD features as cluster stars. The proper motion distribution of these stars is distinguishable from that of the cluster, with some superposition, which resembles that of stars located beyond 3.5deg from the cluster center. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 7 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 643, A15 (2020)

A revised view of the Canis Major stellar overdensity with DECam and Gaia: new evidence of a stellar warp of blue stars

Authors: Julio A. Carballo-Bello, David Martínez-Delgado, Jesús M. Corral-Santana, Emilio J. Alfaro, Camila Navarrete, A. Katherina Vivas, Márcio Catelan

Abstract: We present DECam imaging combined with Gaia DR2 data to study the Canis Major overdensity. The presence of the so-called Blue Plume stars in a low-pollution area of the color-magnitude diagram allows us to derive the distance and proper motions of this stellar feature along the line of sight of its hypothetical core. The stellar overdensity extends on a large area of the sky at low Galactic latitu… ▽ More We present DECam imaging combined with Gaia DR2 data to study the Canis Major overdensity. The presence of the so-called Blue Plume stars in a low-pollution area of the color-magnitude diagram allows us to derive the distance and proper motions of this stellar feature along the line of sight of its hypothetical core. The stellar overdensity extends on a large area of the sky at low Galactic latitudes, below the plane, and between 230$^\circ < \ell < 255^\circ$. According to the orbit derived for Canis Major, it presents an on-plane rotation around the Milky Way. Moreover, additional overdensities of Blue Plume stars are found around the plane and across the Galaxy, proving that these objects are not only associated with that structure. The spatial distribution of these stars, derived using Gaia astrometric data, confirms that the detection of the Canis Major overdensity results more from the warped structure of the Milky Way disk than from the accretion of a dwarf galaxy. △ Less

Submitted 3 September, 2020; originally announced September 2020.

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

Recovering variable stars in large surveys: EA$_{\rm up}$ Algol-type class in the Catalina Survey

Authors: A. Carmo, C. E. Ferreira Lopes, A. Papageorgiou, F. J. Jablonski, C. V. Rodrigues, A. J. Drake, N. J. G. Cross, M. Catelan

Abstract: The discovery and characterization of Algol eclipsing binaries (EAs) provide an opportunity to contribute for a better picture of the structure and evolution of low-mass stars. However, the cadence of most current photometric surveys hinders the detection of EAs since the separation between observations is usually larger than the eclipse(s) duration and hence few measurements are found at the ecli… ▽ More The discovery and characterization of Algol eclipsing binaries (EAs) provide an opportunity to contribute for a better picture of the structure and evolution of low-mass stars. However, the cadence of most current photometric surveys hinders the detection of EAs since the separation between observations is usually larger than the eclipse(s) duration and hence few measurements are found at the eclipses. Even when those objects are detected as variable, their periods can be missed if an appropriate oversampling factor is not used in the search tools. In this paper, we apply this approach to find the periods of stars cataloged in the Catalina Real-Time Transient Survey (CRTS) as EAs having unknown period (EA$_{\rm up}$). As a result, the periods of $\sim 56\%$ of them were determined. Eight objects were identified as low-mass binary systems and modeled with the Wilson \& Devinney synthesis code combined with a Monte-Carlo Markov Chain optimization procedure. The computed masses and radii are in agreement with theoretical models and show no evidence of inflated radii. This paper is the first of a series aiming to identify suspected binary systems in large surveys. △ Less

Submitted 7 September, 2020; v1 submitted 17 August, 2020; originally announced August 2020.

Comments: 11 pages, 7 figures

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

Using classical Cepheids to study the far side of the Milky Way disk: I. Spectroscopic classification and the metallicity gradient

Authors: J. H. Minniti, L. Sbordone, A. Rojas-Arriagada, M. Zoccali, R. Contreras Ramos, D. Minniti, M. Marconi, V. F. Braga, M. Catelan, S. Duffau, W. Gieren, A. A. R. Valcarce

Abstract: The structure, kinematics, and chemical composition of the far side of the Milky Way disk, beyond the bulge, are still to be revealed. Classical Cepheids (CCs) are young and luminous standard candles. We aim to use a well-characterized sample of these variable stars to study the present time properties of the far side of the Galactic disk. A sample of 45 Cepheid variable star candidates were selec… ▽ More The structure, kinematics, and chemical composition of the far side of the Milky Way disk, beyond the bulge, are still to be revealed. Classical Cepheids (CCs) are young and luminous standard candles. We aim to use a well-characterized sample of these variable stars to study the present time properties of the far side of the Galactic disk. A sample of 45 Cepheid variable star candidates were selected from near infrared time series photometry obtained by the VVV survey. We characterized this sample using high quality near infrared spectra obtained with VLT/X-Shooter, deriving radial velocities and iron abundances for all the sample Cepheids. This allowed us to separate the CCs, which are metal rich and with kinematics consistent with the disk rotation, from type II Cepheids (T2Cs), which are more metal poor and with different kinematics. We estimated individual distances and extinctions using VVV photometry and period-luminosity relations, reporting the characterization of 30 CCs located on the far side of the Galactic disk, plus 8 T2Cs mainly located in the bulge region, of which 10 CCs and 4 T2Cs are new discoveries. This is the first sizeable sample of CCs in this distant region of our Galaxy that has been spectroscopically confirmed. We use their positions, kinematics, and metallicities to confirm that the general properties of the far disk are similar to those of the well-studied disk on the solar side of the Galaxy. In addition, we derive for the first time the radial metallicity gradient on the disk's far side. Considering all the CCs with $R_{\mathrm{GC}} < 17\,\rm{kpc}$, we measure a gradient with a slope of $-0.062 \, \mathrm{dex\, kpc^{-1}}$ and an intercept of $+0.59 \, \rm{dex}$, which is in agreement with previous determinations based on CCs on the near side of the disk. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 16 pages, 10 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 640, A92 (2020)

The VVV Infrared Variability Catalog (VIVA-I)

Authors: C. E. Ferreira Lopes, N. J. G. Cross, M. Catelan, D. Minniti, M. Hempel, P. W. Lucas, R. Angeloni, F. Jablonsky, V. F. Braga, I. C. Leao, F. R. Herpich, J. Alonso-Garcia, A. Papageorgiou, K. Pichara, R. K. Saito, A. Bradley, J. C. Beamin, C. Cortes, J. R. De Medeiros, Christopher. M. P. Russell

Abstract: Thanks to the VISTA Variables in the Via Lactea (VVV) ESO Public Survey it is now possible to explore a large number of objects in those regions. This paper addresses the variability analysis of all VVV point sources having more than 10 observations in VVVDR4 using a novel approach. In total, the near-IR light curves of 288,378,769 sources were analysed using methods developed in the New Insight I… ▽ More Thanks to the VISTA Variables in the Via Lactea (VVV) ESO Public Survey it is now possible to explore a large number of objects in those regions. This paper addresses the variability analysis of all VVV point sources having more than 10 observations in VVVDR4 using a novel approach. In total, the near-IR light curves of 288,378,769 sources were analysed using methods developed in the New Insight Into Time Series Analysis project. As a result, we present a complete sample having 44, 998, 752 variable star candidates (VVV-CVSC), which include accurate individual coordinates, near-IR magnitudes (ZYJHKs), extinctions A(Ks), variability indices, periods, amplitudes, among other parameters to assess the science. Unfortunately, a side effect of having a highly complete sample, is also having a high level of contamination by non-variable (contamination ratio of non-variables to variables is slightly over 10:1). To deal with this, we also provide some flags and parameters that can be used by the community to de-crease the number of variable candidates without heavily decreasing the completeness of the sample. In particular, we cross-identified 339,601 of our sources with Simbad and AAVSO databases, which provide us with information for these objects at other wavelegths. This sub-sample constitutes a unique resource to study the corresponding near-IR variability of known sources as well as to assess the IR variability related with X-ray and Gamma-Ray sources. On the other hand, the other 99.5% sources in our sample constitutes a number of potentially new objects with variability information for the heavily crowded and reddened regions of the Galactic Plane and Bulge. The present results also provide an important queryable resource to perform variability analysis and to characterize ongoing and future surveys like TESS and LSST. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 27 pages, 14 figures

Scalable End-to-end Recurrent Neural Network for Variable star classification

Authors: Ignacio Becker, Karim Pichara, Márcio Catelan, Pavlos Protopapas, Carlos Aguirre, Fatemeh Nikzat

Abstract: During the last decade, considerable effort has been made to perform automatic classification of variable stars using machine learning techniques. Traditionally, light curves are represented as a vector of descriptors or features used as input for many algorithms. Some features are computationally expensive, cannot be updated quickly and hence for large datasets such as the LSST cannot be applied.… ▽ More During the last decade, considerable effort has been made to perform automatic classification of variable stars using machine learning techniques. Traditionally, light curves are represented as a vector of descriptors or features used as input for many algorithms. Some features are computationally expensive, cannot be updated quickly and hence for large datasets such as the LSST cannot be applied. Previous work has been done to develop alternative unsupervised feature extraction algorithms for light curves, but the cost of doing so still remains high. In this work, we propose an end-to-end algorithm that automatically learns the representation of light curves that allows an accurate automatic classification. We study a series of deep learning architectures based on Recurrent Neural Networks and test them in automated classification scenarios. Our method uses minimal data preprocessing, can be updated with a low computational cost for new observations and light curves, and can scale up to massive datasets. We transform each light curve into an input matrix representation whose elements are the differences in time and magnitude, and the outputs are classification probabilities. We test our method in three surveys: OGLE-III, Gaia and WISE. We obtain accuracies of about $95\%$ in the main classes and $75\%$ in the majority of subclasses. We compare our results with the Random Forest classifier and obtain competitive accuracies while being faster and scalable. The analysis shows that the computational complexity of our approach grows up linearly with the light curve size, while the traditional approach cost grows as $N\log{(N)}$. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: 15 pages, 17 figures. To be published in MNRAS

Humps and bumps: The effects of shocks on the optical light curves of fundamental-mode RR Lyrae stars

Authors: Z. Prudil, I. Dékány, R. Smolec, M. Catelan, E. K. Grebel, A. Kunder

Abstract: We present the most extended and homogeneous study carried out so far of the main and early shocks in 1485 RR~Lyrae stars in the Galactic bulge observed by the Optical Gravitational Lensing Experiment (OGLE). We selected non-modulated fundamental-mode RR~Lyrae stars with good-quality photometry. Using a self-developed method, we determined the centers and strengths of main and early shock features… ▽ More We present the most extended and homogeneous study carried out so far of the main and early shocks in 1485 RR~Lyrae stars in the Galactic bulge observed by the Optical Gravitational Lensing Experiment (OGLE). We selected non-modulated fundamental-mode RR~Lyrae stars with good-quality photometry. Using a self-developed method, we determined the centers and strengths of main and early shock features in the phased light curves. We found that the position of both humps and bumps are highly correlated with the pulsation properties of the studied variables. Pulsators with a pronounced main shock are concentrated in the low-amplitude regime of the period-amplitude diagram, while stars with a strong early shock have average and above-average pulsation amplitudes. A connection between the main and early shocks and the Fourier coefficients is also observed. In the color-magnitude diagram (CMD), we see a separation between stars with strong and weak shocks. Variables with a pronounced main shock cluster close to the fundamental red edge of the instability strip (IS), while stars with a strong early shock tend to clump in the center and near the fundamental blue edge of the IS. The appearance of shocks and their properties seem independent of the direction of evolution estimated from the period change rate of the studied stars. In addition, the differences in the period change rate between the two main Oosterhoff groups found in the Galactic bulge suggest that stars of Oosterhoff type I are located close to the zero-age horizontal branch while Oosterhoff type II variables are on their way toward the fundamental red edge of the instability strip, thus having already left the zero-age horizontal branch. △ Less

Submitted 22 January, 2020; originally announced January 2020.

Comments: 14 pages, 10 figures, Accepted for publication in A&A

Journal ref: A&A 635, A66 (2020)

VVV-WIT-01: highly obscured classical nova or protostellar collision?

Authors: P. W. Lucas, D. Minniti, A. Kamble, D. L. Kaplan, N. Cross, I. Dekany, V. D. Ivanov, R. Kurtev, R. K. Saito, L. C. Smith, M. Catelan, N. Masetti, I. Toledo, M. Hempel, M. A. Thompson, C. Contreras Peña, J. Forbrich, M. Krause, J. Dale, J. Borissova, J. Emerson

Abstract: A search of the first Data Release of the VISTA Variables in the Via Lactea (VVV) Survey discovered the exceptionally red transient VVV-WIT-01 ($H-K_s=5.2$). It peaked before March 2010, then faded by $\sim$9.5 mag over the following two years. The 1.6--22 $μ$m spectral energy distribution in March 2010 was well fit by a highly obscured black body with $T \sim 1000$ K and $A_{K_s} \sim 6.6$ mag. T… ▽ More A search of the first Data Release of the VISTA Variables in the Via Lactea (VVV) Survey discovered the exceptionally red transient VVV-WIT-01 ($H-K_s=5.2$). It peaked before March 2010, then faded by $\sim$9.5 mag over the following two years. The 1.6--22 $μ$m spectral energy distribution in March 2010 was well fit by a highly obscured black body with $T \sim 1000$ K and $A_{K_s} \sim 6.6$ mag. The source is projected against the Infrared Dark Cloud (IRDC) SDC G331.062$-$0.294. The chance projection probability is small for any single event ($p \approx 0.01$ to 0.02) which suggests a physical association, e.g. a collision between low mass protostars. However, black body emission at $T \sim 1000$ K is common in classical novae (especially CO novae) at the infrared peak in the light curve, due to condensation of dust $\sim$30--60 days after the explosion. Radio follow up with the Australia Telescope Compact Array (ATCA) detected a fading continuum source with properties consistent with a classical nova but probably inconsistent with colliding protostars. Considering all VVV transients that could have been projected against a catalogued IRDC raises the probability of a chance association to $p=0.13$ to 0.24. After weighing several options, it appears likely that VVV-WIT-01 was a classical nova event located behind an IRDC. △ Less

Submitted 15 January, 2020; originally announced January 2020.

Comments: 13 pages, 7 figures. Accepted by MNRAS

Into the Darkness: Classical and Type II Cepheids in the Zona Galactica Incognita

Authors: István Dékány, Gergely Hajdu, Eva K. Grebel, Márcio Catelan

Abstract: The far side of the Milky Way's disk is one of the most concealed parts of the known Universe due to extremely high interstellar extinction and point source density toward low Galactic latitudes. Large time-domain photometric surveys operating in the near-infrared hold great potential for the exploration of these vast uncharted areas of our Galaxy. We conducted a census of distant classical and ty… ▽ More The far side of the Milky Way's disk is one of the most concealed parts of the known Universe due to extremely high interstellar extinction and point source density toward low Galactic latitudes. Large time-domain photometric surveys operating in the near-infrared hold great potential for the exploration of these vast uncharted areas of our Galaxy. We conducted a census of distant classical and type II Cepheids along the southern Galactic mid-plane using near-infrared photometry from the VISTA Variables in the Vía Láctea survey. We performed a machine-learned classification of the Cepheids based on their infrared light curves using a convolutional neural network. We have discovered 640 distant classical Cepheids with up to ~40 magnitudes of visual extinction, and over 500 type II Cepheids, most of them located in the inner bulge. Intrinsic color indices of individual Cepheids were predicted from sparse photometric data using a neural network, allowing their use as accurate reddening tracers. They revealed a steep, spatially varying near-infrared extinction curve toward the inner bulge. Type II Cepheids in the Galactic bulge were also employed to measure robust mean selective-to-absolute extinction ratios. They trace a centrally concentrated spatial distribution of the old bulge population with a slight elongation, consistent with earlier results from RR Lyrae stars. Likewise, the classical Cepheids were utilized to trace the Galactic warp and various substructures of the Galactic disk, and to uncover significant vertical and radial age gradients of the thin disk population at the far side of the Milky Way. △ Less

Submitted 22 August, 2019; originally announced August 2019.

Comments: Accepted for publication by The Astrophysical Journal