Search

We present the first application of data-driven techniques for dynamical system analysis based on Koopman theory to variable stars. We focus on light curves of RRLyrae type variables, in the Galactic globular cluster $\omega$ Centauri. Light curves are thus summarized by a handful of complex eigenvalues, corresponding to oscillatory or fading dynamical modes. We find that variable stars of the RRc subclass can be summarized in terms of fewer ($\approx 8$) eigenvalues, while RRab need comparatively more ($\approx 12$). This result can be leveraged for classification and reflects the simpler structure of RRc light curves. We then consider variable stars displaying secular variations due to the Tseraskaya-Blazhko effect and find a change in relevant eigenvalues with time, with possible implications for the physical interpretation of the effect., Comment: Accepted (spotlight) at the Scaling AI for Science workshop at ICML 2024

Stellar mergers and accretion events have been crucial in shaping the evolution of the Milky Way (MW). These events have been dynamically identified and chemically characterised using red giants and main-sequence stars. RR Lyrae (RRL) variables can play a crucial role in tracing the early formation of the MW since they are ubiquitous, old (t$\ge$10 Gyr) low-mass stars and accurate distance indicators. We exploited Data Release 3 of the GALAH survey to identify 78 field RRLs suitable for chemical analysis. Using synthetic spectra calculations, we determined atmospheric parameters and abundances of Fe, Mg, Ca, Y, and Ba. Most of our stars exhibit halo-like chemical compositions, with an iron peak around [Fe/H]$\approx -$1.40, and enhanced Ca and Mg content. Notably, we discovered a metal-rich tail, with [Fe/H] values ranging from $-$1 to approximately solar metallicity. This sub-group includes almost ~1/4 of the sample, it is characterised by thin disc kinematics and displays sub-solar $\alpha$-element abundances, marginally consistent with the majority of the MW stars. Surprisingly, they differ distinctly from typical MW disc stars in terms of the s-process elements Y and Ba. We took advantage of similar data available in the literature and built a total sample of 535 field RRLs for which we estimated kinematical and dynamical properties. We found that metal-rich RRLs (1/3 of the sample) likely represent an old component of the MW thin disc. We also detected RRLs with retrograde orbits and provided preliminary associations with the Gaia-Sausage-Enceladus, Helmi, Sequoia, Sagittarius, and Thamnos stellar streams., Comment: Accepted for publication in MNRAS. 29 pages, 20 figures

Binary Cepheids with giant companions are crucial for studying the physical properties of Cepheid variables, providing the best means to measure their masses. Systems composed of two Cepheids are even more important but to date, only one such system in the Large Magellanic Cloud (LMC) was known. Our current aim is to increase the number of these systems tenfold and provide their basic characteristics. The final goal is to obtain the physical properties of the component Cepheids, including their masses and radii, and to learn about their evolution in the multiple systems, also revealing their origin. We started a spectroscopic monitoring of nine unresolved pairs of Cepheids from the OGLE catalog, to check if they are gravitationally bound. Two of these so-called double Cepheids are located in the LMC, five in the Small Magellanic Cloud (SMC), and two in the Milky Way (MW). We report the spectroscopic detection of binarity of all 9 of these double Cepheids with orbital periods from 2 to 18 years. This increases the number of known binary double (BIND) Cepheids from 1 to 10 and triples the number of all confirmed double-lined binary (SB2) Cepheids. For five BIND Cepheids disentangled pulsational light curves of the components show anti-correlated phase shifts due to orbital motion. We show the first empirical evidence that typical period-luminosity relations (PLRs) are rather binary Cepheid PLRs that include the companion's light. The statistics of pulsation period ratios of BIND Cepheids do not agree with those expected for pairs of the same-age Cepheids. These ratios together with the mass ratios far from unity suggest merger-origin of at least one component for about half of the systems. The SMC and MW objects are the first found in SB2 systems composed of giants in their host galaxies. The Milky Way BIND Cepheids are also the closest such systems, being located at about 11 and 26 kpc., Comment: 13 pages, 7 figures, 4 tables, accepted for publication in Astronomy & Astrophysics

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, Comment: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

Outflows in the Active Galactic Nucleus (AGN) are considered to play a key role in the host galaxy evolution through transfer of a large amount of energy. A Narrow Line Region (NLR) in the AGN is composed of ionized gas extending from pc-scales to kpc-scales. It has been suggested that shocks are required for ionization of the NLR gas. If AGN outflows generate such shocks, they will sweep through the NLR and the outflow energy will be transferred into a galaxy-scale region. In order to study contribution of the AGN outflow to the NLR-scale shock, we measure the [\ion{Fe}{2}]$\lambda12570$/[\ion{P}{2}]$\lambda11886$ line ratio, which is a good tracer of shocks, using near-infrared spectroscopic observations with WINERED (Warm INfrared Echelle spectrograph to Realize Extreme Dispersion and sensitivity) mounted on the New Technology Telescope. Among 13 Seyfert galaxies we observed, the [\ion{Fe}{2}] and [\ion{P}{2}] lines were detected in 12 and 6 targets, respectively. The [\ion{Fe}{2}]/[\ion{P}{2}] ratios in 4 targets were found to be higher than 10, which implies the existence of shocks. We also found that the shock is likely to exist where an ionized outflow, i.e., a blue wing in [\ion{S}{3}]$\lambda9533$, is present. Our result implies that the ionized outflow present over a NLR-scale region sweeps through the interstellar medium and generates a shock., Comment: Accepted for ApJ, 20 pages, 11 figures

The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value., Comment: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

Wide-field and deep DECam multi-band photometry, combined with HST data for the core of the Galactic globular cluster NGC 2808, allowed us to study the distribution of various stellar sub-populations and stars in different evolutionary phases out to the cluster tidal radius. We used the C_ugi = (u-g)-(g-i) index to identify three chemically distinct sub-populations along the red giant branch and compared their spatial distributions. The most light-element enriched sub-population (P3) is more centrally concentrated; however, it shows a more extended distribution in the external regions of the cluster compared to the primordial (P1) and intermediate (P2) composition populations. Furthermore, the P3 sub-population centroid is off-center relative to those of the P1 and P2 groups. We also analyzed the spatial distribution of horizontal branch stars and found that the relative fraction of red horizontal branch stars increases for radial distances larger than ~ 1.5' while that of the blue and hotter stars decreases. These new observations, combined with literature spectroscopic measurements, suggest that the red horizontal branch stars are the progeny of all the stellar sub-populations in NGC 2808, i.e. primordial and light-element enhanced, while the blue stars are possibly the result of a combination of the "hot-flasher" and the "helium-enhanced" scenarios. A similar distribution of different red giant branch sub-populations and horizontal branch stars was also found for the most massive Galactic globular cluster, omega Cen, based on combined DECam and HST data, which suggests the two may share a similar origin., Comment: Accepted for publication in AJ; 22 pages; 16 figures; the full version of table 2 will be available in electronic form with the published version

RR Lyrae (RRL) are old, low-mass radially pulsating variable stars in their core helium burning phase. They are popular stellar tracers and primary distance indicators, since they obey to well defined period-luminosity relations in the near-infrared regime. Their photometric identification is not trivial, indeed, RRL samples can be contaminated by eclipsing binaries, especially in large datasets produced by fully automatic pipelines. Interpretable machine-learning approaches for separating eclipsing binaries from RRL are thus needed. Ideally, they should be able to achieve high precision in identifying RRL while generalizing to new data from different instruments. In this paper, we train a simple logistic regression classifier on Catalina Sky Survey (CSS) light curves. It achieves a precision of 87% at 78% recall for the RRL class on unseen CSS light curves. It generalizes on out-of-sample data (ASAS/ASAS-SN light curves) with a precision of 85% at 96% recall. We also considered a L1-regularized version of our classifier, which reaches 90% sparsity in the light-curve features with a limited trade-off in accuracy on our CSS validation set and -- remarkably -- also on the ASAS/ASAS-SN light curve test set. Logistic regression is natively interpretable, and regularization allows us to point out the parts of the light curves that matter the most in classification. We thus achieved both good generalization and full interpretability., Comment: To appear on The Astrophysical Journal. 13 pages, 8 figures, 3 tables

RR Lyrae variables are excellent population II distance indicators thanks to their well-defined period-luminosity relations (PLRs) at infrared wavelengths. We present results of near-infrared (NIR) monitoring of Galactic globular clusters to empirically quantify the metallicity dependence of NIR PLRs for RR Lyrae variables. Our sample includes homogeneous, accurate, and precise photometric data for 964 RR Lyrae variables in 11 globular clusters covering a large metallicity range ($\Delta\textrm{[Fe/H]}\sim2$~dex). We derive $JHK_s$ band period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations anchored using 346 Milky Way field RR Lyrae stars with {\it Gaia} parallaxes, and simultaneously solved for independent distances to globular clusters. We find a significant metallicity dependence of $\sim0.2$~mag/dex in $JHK_s$ band PLZ and PWZ relations for RR Lyrae stars independent of the adopted metallicity scale. The metallicity coefficients and the zero-points of the empirical PLZ and PWZ relations are in excellent agreement with the predictions from the horizontal branch evolution and pulsation models. Furthermore, RR Lyrae based distances to our sample of globular clusters are also statistically consistent with other independent measurements in the literature. Our recommended empirical $JHK_s$ band PLZ relations are also provided for RR Lyrae based distance measurements., Comment: 10 pages, 4 figures, abridged abstract, accepted in the ApJ Letters

We present new empirical infrared Period-Luminosity-Metallicity (PLZ) and Period-Wesenheit-Metallicity (PWZ) relations for RR Lyrae based on the latest Gaia EDR3 parallaxes. The relations are provided in the WISE $W1$ and $W2$ bands, as well as in the $W(W1, V - W1)$ and $W(W2, V - W2)$ Wesenheit magnitudes. The relations are calibrated using a very large sample of Galactic halo field RR Lyrae stars with homogeneous spectroscopic [Fe/H] abundances (over 1,000 stars in the $W1$ band), covering a broad range of metallicities ($-2.5 \lesssim \textrm{[Fe/H]} \lesssim 0.0$). We test the performance of our PLZ and PWZ relations by determining the distance moduli of both galactic and extragalactic stellar associations: the Sculptor dwarf spheroidal galaxy in the Local Group (finding $\bar{\mu}_{0}=19.47 \pm 0.06$), the Galactic globular clusters M4 ($\bar{\mu}_{0}=11.16 \pm 0.05$) and the Reticulum globular cluster in the Large Magellanic Cloud ($\bar{\mu}_{0}=18.23 \pm 0.06$). The distance moduli determined through all our relations are internally self-consistent (within $\lesssim$ 0.05 mag) but are systematically smaller (by $\sim$ 2-3$\sigma$) than previous literature measurements taken from a variety of methods/anchors. However, a comparison with similar recent RR Lyrae empirical relations anchored with EDR3 likewise shows to varying extents a systematically smaller distance modulus for PLZ/PWZ RR Lyrae relations., Comment: Accepted by ApJ, 14 pages, 5 Figures, 2 Tables

We report the discovery of a surprising binary configuration of the double-mode Cepheid OGLE-LMC-CEP-1347 pulsating in the first (P_1=0.690d) and second overtone (P_2=0.556d) modes. The orbital period (P_orb=59d) of the system is five times shorter than the shortest known to date (310d) for a binary Cepheid. The Cepheid itself is also the shortest-period one ever found in a binary system and the first double-mode Cepheid in a spectroscopically double-lined binary. OGLE-LMC-CEP-1347 is most probably on its first crossing through the instability strip, as inferred from both its short period and fast period increase, consistent with evolutionary models, and from the short orbital period (not expected for binary Cepheids whose components have passed through the red giant phase). Our evolutionary analysis yielded a first-crossing Cepheid with a mass in a range of 2.9-3.4 Msun (lower than any measured Cepheid mass), consistent with observations. The companion is a stable star, at least two times fainter and less massive than the Cepheid (preliminary mass ratio q=0.55), while also redder and thus at the subgiant or more advanced evolutionary stage. To match these characteristics, the Cepheid has to be a product of binary interaction, most likely a merger of two less massive stars, which makes it the second known classical Cepheid of binary origin. Moreover, further evolution of the components may lead to another binary interaction., Comment: 6 pages, 4 figures, 1 table, published in The Astrophysical Journal Letters

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

The revolutionary power of future Rubin-LSST observations will allow us to significantly improve the physics of pulsating stars, including RR Lyrae. In this context, an updated theoretical scenario predicting all the relevant pulsation observables in the corresponding photometric filters is mandatory. The bolometric light curves based on a recently computed extensive set of nonlinear convective pulsation models for RR Lyrae stars, covering a broad range in metal content and transformed into the Rubin-LSST photometric system. Predicted Rubin-LSST mean magnitudes and pulsation amplitudes have been adopted to built the Bailey diagrams (luminosity amplitude vs period) and the color-color diagrams in these bands. The current findings indicate that the gLSST-rLSST, rLSST-iLSST colors obey to a well defined linear relation with the metal content. Moreover, the Period Luminosity relations display in the reddest filters (rLSST,iLSST,zLSST,yLSST) a significant dependence on the assumed metal abundance. In particular, more metal-rich RR Lyrae are predicted to be fainter at fixed period. Metal-dependent Period-Wesenheit relations for different combinations of optical and NIR filters are also provided. These represent powerful tools to infer individual distances independently of reddening uncertainties, once the metal abundance is known and no relevant deviations from the adopted extinction law occur. Finally, we also derived new linear and quadratic absolute magnitude metallicity relations (gLSST vs [Fe/H]) and the metallicity coefficient is consistent with previous findings concerning the B and the V band., Comment: Accepted for publication on ApJ

We searched for diffuse interstellar bands (DIBs) in the 0.91$<\lambda<$1.33 $\mu$m region by analyzing the near-infrared (NIR) high-resolution ($R=20,000$ and 28,000) spectra of 31 reddened early-type stars ($0.04

We present accurate and deep multi-band ($g,r,i$) photometry of the Local Group dwarf irregular galaxy NGC 6822. The images were collected with wide field cameras at 2m/4m- (INT,CTIO,CFHT) and 8m-class telescopes (SUBARU) covering a 2 square degrees FoV across the center of the galaxy. We performed PSF photometry of $\approx$7,000 CCD images and the final catalog includes more than 1 million objects. We developed a new approach to identify candidate field and galaxy stars, and performed a new estimate of the galaxy center by using old stellar tracers finding that it differs by 1.15 (RA) and 1.53 (DEC) arcmin from previous estimates. We also found that young (Main Sequence, Red Supergiants), intermediate (Red Clump, Asymptotic Giant Branch [AGB]) and old (Red Giant Branch [RGB]) stars display different radial distributions. Old stellar population is spherically distributed and extends to radial distances larger than previously estimated ($\sim$1 degree). The young population shows a well defined bar and a disk-like distribution, as suggested by radio measurements, that is off-center compared with old population. We discuss pros and cons of the different diagnostics adopted to identify AGB stars and develop new ones based on optical-NIR-MIR color-color diagrams (CCDs) to characterize Oxygen and Carbon (C) rich stars. We found a mean population ratio between Carbon and M-type (C/M) stars of 0.67$\pm$0.08 (optical/NIR/MIR) and we used the observed C/M ratio with empirical C/M-metallicity relations to estimate a mean iron abundance of [Fe/H]$\sim$-1.25 ($\sigma$=0.04 dex) that agrees quite well with literature estimates., Comment: Accepted for publication in ApJ, 34 pages, 22 figures, 6 tables

We present new period-$\phi_{31}$-[Fe/H] relations for first overtone RRL stars (RRc), calibrated over a broad range of metallicities ($-2.5 < \textrm{[Fe/H]}< 0.0$) utilizing the largest currently available set of Galactic halo field RRL with homogeneous spectroscopic metallicities. Our relations are defined in the optical (ASAS-SN $V$-band) and, inaugurally, in the infrared (WISE $W1$ and $W2$ bands). Our $V$-band relation can reproduce individual RRc spectroscopic metallicities with a dispersion of 0.30 dex over the entire metallicity range of our calibrator sample (an RMS smaller than what we found for other relations in literature including non-linear terms). Our infrared relation has a similar dispersion in the low and intermediate metallicity range ($\textrm{[Fe/H]} < -0.5$) but tends to underestimate the [Fe/H] abundance around solar metallicity. We tested our relations by measuring both the metallicity of the Sculptor dSph and a sample of Galactic globular clusters, rich in both RRc and RRab stars. The average metallicity we obtain for the combined RRL sample in each cluster is within $\pm 0.08$ dex of their spectroscopic metallicities. The infrared and optical relations presented in this work will enable deriving reliable photometric RRL metallicities in conditions where spectroscopic measurements are not feasible; e.g., in distant galaxies or reddened regions (observed with upcoming Extremely Large Telescopes and the James Webb Space Telescope), or in the large sample of new RRL that will be discovered in large-area time-domain photometric surveys (such as LSST and the Roman space telescope)., Comment: Accepted by ApJ, 22 pages, 9 Figures, 3 Tables

The Extremely Large Telescope and the Thirty Meter Telescope will use state of the art multiconjugate adaptive optics (MCAO) systems to obtain the full D4 advantage that their apertures can provide. However, to reach the full astrometric potential of these facilities for on-sky science requires understanding any residual astrometric distortions from these systems and find ways to measure and eliminate them. In this work, we use Gemini multiconjugate adaptive optic system (GeMS) observations of the core of NGC 6723 to better understand the on-sky astrometric performance of MCAO. We develop new methods to measure the astrometric distortion fields of the observing system, which probe the distortion at the highest possible spatial resolution. We also describe methods for examining the time-variable and static components of the astrometric distortion. When applied to the GeMS Gemini South Adaptive Optics Imager (GSAOI) data, we are able to see the effect of the field rotator at the subpixel level, and we are able to empirically derive the distortion due to the optical design of GeMS-GSAOI. We argue that the resulting distortion maps are a valuable tool to measure and monitor the on-sky astrometric performance of future instrumentation. Our overall astrometry pipeline produces high-quality proper motions with an uncertainty floor of 45 uas per year. We measure the proper motion dispersion profile of NGC 6723 from a radius of 10 arcsec out to 1 arcmin based on 12000 stars. We also produce a high-quality optical-near infrared color-magnitude diagram, which clearly shows the extreme horizontal branch and main-sequence knee of this cluster.

Metal-deficient stars are important tracers for understanding the early formation of the Galaxy. Recent large-scale surveys with both photometric and spectroscopic data have reported an increasing number of metal-deficient stars whose kinematic features are consistent with those of the disk stellar populations. We report the discovery of an RR~Lyrae variable (hereafter RRL) that is located within the thick disk and has an orbit consistent with the thick-disk kinematics. Our target RRL (HD 331986) is located at around 1 kpc from the Sun and, with V=11.3, is among the 130 brightest RRLs known so far. However, this object was scarcely studied because it is in the midplane of the Galaxy, the Galactic latitude around -1 deg. Its near-infrared spectrum (0.91-1.32 micron) shows no absorption line except hydrogen lines of the Paschen series, suggesting [Fe/H] less than -2.5. It is the most metal-deficient RRL, at least, among the RRLs whose orbits are consistent with the disk kinematics, although we cannot determine to which of the disk and the halo it belongs. This unique RRL would provide us with essential clues for studying the early formation of stars in the inner Galaxy with further investigations, including high-resolution optical spectroscopy., Comment: 13 pages, 7 figures, 6 tables, accepted for publication in ApJ