Transient low-mass X-ray binaries (LMXBs) are discovered largely by X-ray and gamma-ray all-sky monitors. The X-ray outburst is also accompanied by an optical brightening, which empirically can precede the detection of X-rays. Newly sensitive optical synoptic surveys may offer a complementary pathway for discovery and potential for insight into the initial onset and propagation of the thermal instability that leads to the ionization of the accretion disk. We use the Zwicky Transient Facility (ZTF) alert stream to perform a comprehensive search at optical wavelengths for previously undiscovered outbursting LMXBs. Our pipeline first crossmatches the positions of the alerts to cataloged X-ray sources, and then analyzes the 30 day lightcurve of matched alerts by thresholding on differences with an 8 day exponentially weighted moving average. In addition to a 19 month long live search, we ran our pipeline over 4 yr of ZTF archival data, recovering four known LMXBs. We also independently detected an outburst of MAXI J1957+032 in the live search and found the first outburst of Swift J1943.4+0228, an unclassified X-ray transient, in 10 yr. Using Monte Carlo simulations of the Galactic LMXB population, we estimate that 29% of outbursting LMXBs are detectable by ZTF and that 4.4% of LMXBs would be present in the crossmatched X-ray catalogs, giving an estimated Galactic population of ${3390}_{-1930}^{+3980}$ . We estimate that our current pipeline can detect 1.3% of all outbursting LMXBs, including those previously unknown, but that Rubin Observatory’s Legacy Survey of Space and Time will be able to detect 43% of outbursting LMXBs.

The MOOSE (Monitoring Observations of SMC X-1 Excursions) program uses the Neutron Star Interior Composition Explorer Mission (NICER) to monitor the high mass X-ray binary SMC X-1 during its superorbital period excursions. Here we perform X-ray spectral analyses of 26 NICER observations of SMC X-1, taken at the tail-end of the excursion between 2021-04-01 and 2022-01-05. We use a single spectral model to fit spectra observed in high, intermediate and low states, using a combination of a partial covering fraction model, a black-body disc, and a power-law component. We find that the partial covering fraction varies significantly with the superorbital state during superorbital excursion. Our findings suggest that the low/high state in SMC X-1 is caused by a very high obscuration of the accretion disk.

Context. The inner Galaxy is a complex environment, and the relative contributions of different formation scenarios to its observed morphology and stellar properties are still debated. The different components are expected to have different spatial, kinematic, and metallicity distributions, and a combination of photometric, spectroscopic, and astrometric large-scale surveys is needed to study the formation and evolution of the Galactic bulge. Aims. The Blanco DECam Bulge Survey (BDBS) provides near-ultraviolet to near-infrared photometry for approximately 250 million unique stars over more than 200 square degrees of the southern Galactic bulge. By combining BDBS photometry with the latest Gaia astrometry, we aim to characterize the chemodynamics of red clump stars across the BDBS footprint using an unprecedented sample size and sky coverage. Methods. Our field of view of interest is |ℓ| ≤ 10°, −10° ≤b ≤ −3°. We constructed a sample of approximately 2.3 million red clump giants in the bulge with photometric metallicities, BDBS photometric distances, and proper motions. Photometric metallicities are derived from a (u − i)0 versus [Fe/H] relation; astrometry, including precise proper motions, is from the third data release (DR3) of the ESA satellite Gaia. We studied the kinematics of the red clump stars as a function of sky position and metallicity by investigating proper-motion rotation curves, velocity dispersions, and proper-motion correlations across the southern Galactic bulge. Results. By binning our sample into eight metallicity bins in the range of −1.5 dex < [Fe/H] < +1 dex, we find that metal-poor red clump stars exhibit lower rotation amplitudes, at ∼29 km s−1 kpc−1. The peak of the angular velocity is ∼39 km s−1 kpc−1 for [Fe/H] ∼ −0.2 dex, exhibiting declining rotation at higher [Fe/H]. The velocity dispersion is higher for metal-poor stars, while metal-rich stars show a steeper gradient with Galactic latitude, with a maximum dispersion at low latitudes along the bulge minor axis. Only metal-rich stars ([Fe/H] ≳ −0.5 dex) show clear signatures of the bar in their kinematics, while the metal-poor population exhibits isotropic motions with an axisymmetric pattern around Galactic longitude ℓ = 0. Conclusions. This work describes the largest sample of bulge stars with distance, metallicity, and astrometry reported to date, and shows clear kinematic differences with metallicity. The global kinematics over the bulge agrees with earlier studies. However, we see striking changes with increasing metallicity, and, for the first time, kinematic differences for stars with [Fe/H]> − 0.5, suggesting that the bar itself may have kinematics that depends on metallicity.

The Milky Way Bulge extra-tidal star survey is a spectroscopic survey with the goal of identifying stripped globular cluster stars from inner Galaxy clusters. In this way, an indication of the fraction of metal-poor bulge stars that originated from globular clusters can be determined. We observed and analyzed stars in and around BH 261, an understudied globular cluster in the bulge. From seven giants within the tidal radius of the cluster, we measured an average heliocentric radial velocity of 〈RV〉 = −61 ± 2.6 km s ^−1 with a radial velocity dispersion of 〈 σ 〉 = 6.1 ± 1.9 km s ^−1 . The large velocity dispersion may have arisen from tidal heating in the cluster’s orbit about the Galactic center, or because BH 261 has a high dynamical mass as well as a high mass-to-light ratio. From spectra of five giants, we measure an average metallicity of 〈[Fe/H]〉 = −1.1 ± 0.2 dex. We also spectroscopically confirm an RR Lyrae star in BH 261, which yields a distance to the cluster of 7.1 ± 0.4 kpc. Stars with 3D velocities and metallicities consistent with BH 261 reaching to ∼0.°5 from the cluster are identified. A handful of these stars are also consistent with the spatial distribution of potential debris from models focusing on the most recent disruption of the cluster.

Context. The inner Galaxy is a complex environment, and the relative contributions of di fferent formation scenarios to its observed morphology and stellar properties are still debated. The di fferent components are expected to have di fferent spatial, kinematic, and metallicity distributions, and a combination of photometric, spectroscopic, and astrometric large-scale surveys is needed to study the formation and evolution of the Galactic bulge. Aims. The Blanco DECam Bulge Survey (BDBS) provides near-ultraviolet to near-infrared photometry for approximately 250 million unique stars over more than 200 square degrees of the southern Galactic bulge. By combining BDBS photometry with the latest Gaia astrometry, we aim to characterize the chemodynamics of red clump stars across the BDBS footprint using an unprecedented sample size and sky coverage. Methods. Our field of view of interest is |l| <= 10 degrees, -10 degrees <= b -3 degrees. We constructed a sample of approximately 2 :3 million red clump giants in the bulge with photometric metallicities, BDBS photometric distances, and proper motions. Photometric metallicities are derived from a (u- i)(0) versus [Fe=H] relation; astrometry, including precise proper motions, is from the third data release (DR3) of the ESA satellite Gaia. We studied the kinematics of the red clump stars as a function of sky position and metallicity by investigating proper-motion rotation curves, velocity dispersions, and proper-motion correlations across the southern Galactic bulge. Results. By binning our sample into eight metallicity bins in the range of 1.5 dex < [Fe/H] < +1 dex, we find that metal-poor red clump stars exhibit lower rotation amplitudes, at similar to 29 km s 1 kpc 1. The peak of the angular velocity is similar to 39 km s 1 kpc 1 for [Fe/H] similar to 0.2 dex, exhibiting declining rotation at higher [Fe/H]. The velocity dispersion is higher for metal-poor stars, while metalrich stars show a steeper gradient with Galactic latitude, with a maximum dispersion at low latitudes along the bulge minor axis. Only metal-rich stars ([Fe/H] >= 0.5 dex) show clear signatures of the bar in their kinematics, while the metal-poor population exhibits isotropic motions with an axisymmetric pattern around Galactic longitude l = 0. Conclusions. This work describes the largest sample of bulge stars with distance, metallicity, and astrometry reported to date, and shows clear kinematic di fferences with metallicity. The global kinematics over the bulge agrees with earlier studies. However, we see striking changes with increasing metallicity, and, for the first time, kinematic di fferences for stars with [Fe/H] > -0.5, suggesting that the bar itself may have kinematics that depends on metallicity.

Patchick 99 is a candidate globular cluster located in the direction of the Galactic bulge, with a proper motion almost identical to the field and extreme field star contamination. A recent analysis suggests it is a low-luminosity globular cluster with a population of RR Lyrae stars. We present new spectra of stars in and around Patchick 99, targeting specifically the 3 RR Lyrae stars associated with the cluster as well as the other RR Lyrae stars in the field. A sample of 53 giant stars selected from proper motions and a position on CMD are also observed. The three RR Lyrae stars associated with the cluster have similar radial velocities and distances, and two of the targeted giants also have radial velocities in this velocity regime and [Fe/H] metallicities that are slightly more metal-poor than the field. Therefore, if Patchick 99 is a bonafide globular cluster, it would have a radial velocity of -92+/-10 km s-1, a distance of 6.7+/-0.4 kpc (as determined from the RR Lyrae stars), and an orbit that confines it to the inner bulge.

A complete map of the youngest stellar populations of the Milky Way in the era of all-sky surveys is one of the most challenging goals in modern astrophysics. The characterization of the youngest stellar components is crucial not only for a global overview of the Milky Way’s structure, of the Galactic thin disk, and its spiral arms, but also for local studies. In fact, the identification of star-forming regions (SFRs) and the comparison with the environment in which they form are also fundamental to put SFRs in the context of the surrounding giant molecular clouds and to understand still unknown physical mechanisms related to star and planet formation processes. In 10 yr of observations, the Vera C. Rubin Legacy Survey of Space and Time (Rubin LSST) will achieve an exquisite photometric depth that will allow us to significantly extend the volume within which we will be able to discover new SFRs and to enlarge the region of our own Galaxy we have detailed knowledge about. We describe here a metric that estimates the total number of young stars with ages t < 10 Myr and masses >0.3 M _⊙ that will be detected with the Rubin LSST observations in the gri bands at a 5 σ magnitude significance. We examine the results of our metric adopting the most recent simulated Rubin LSST survey strategies in order to evaluate the impact that different observing strategies might have on our science case.

The Blanco DECam Bulge Survey (BDBS) provides near-ultraviolet to near-infrared photometry for ~250 million unique stars. By combining BDBS photometry with the latest Gaia astrometry, we characterize the chemo-dynamics of red clump stars across the BDBS footprint, using an unprecedented sample size and sky coverage. We construct a sample of ~2.3 million red clump giants in the bulge with photometric metallicities, BDBS photometric distances, and proper motions. We study the kinematics of the red clump stars as a function of sky position and metallicity, by investigating proper motion rotation curves, velocity dispersions, and proper motion correlations across the southern Galactic bulge. We find that metal-poor red clump stars exhibit lower rotation amplitudes, at ~29 km s$^{-1}$ kpc^{-1}. The peak of the angular velocity is ~39 km s^{-1} kpc^{-1} for [Fe/H] ~ -0.2 dex, exhibiting declining rotation at higher [Fe/H]. The velocity dispersion is higher for metal-poor stars, while metal-rich stars show a steeper gradient with Galactic latitude, with a maximum dispersion at low latitudes along the bulge minor axis. Only metal-rich stars ([Fe/H] >~ -0.5 dex) show clear signatures of the bar in their kinematics, while the metal-poor population exhibits isotropic motions with an axisymmetric pattern around Galactic longitude l = 0. This work reports the largest sample of bulge stars with distance, metallicity, and astrometry and shows clear kinematic differences with metallicity. The global kinematics over the bulge agrees with earlier studies. However, we see striking changes with increasing metallicity and for the first time, see kinematic differences for stars with [Fe/H]>-0.5, suggesting that the bar itself may have kinematics that depends on metallicity.

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.

The Nancy Grace Roman Telescope's High Latitude Wide Area Survey will have a number of synergies with the Vera Rubin Observatory's Legacy Survey of Space and Time (LSST), particularly for extragalactic star clusters. Understanding the nature of star clusters and star cluster systems are key topics in many areas of astronomy, chief among them stellar evolution, high energy astrophysics, galaxy assembly/dark matter, the extragalactic distance scale, and cosmology. One of the challenges will be disentangling the age/metallicity degeneracy because young ($\sim$Myr) metal-rich clusters have similar SEDs to old ($\sim$Gyr) metal-poor clusters. Rubin will provide homogeneous, $ugrizy$ photometric coverage, and measurements in the red Roman filters will help break the age-metallicity and age-extinction degeneracies, providing the first globular cluster samples that cover wide areas while essentially free of contamination from Milky Way stars. Roman's excellent spatial resolution will also allow measurements of cluster sizes. We advocate for observations of a large sample of galaxies with a range of properties and morphologies in the Rubin/LSST footprint matching the depth of the LSST Wide-Fast-Deep field $i$ band limit (26.3 mag), and recommend adding the F213 filter to the survey.

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution ( 0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe.

ABSTRACT This is the third paper in a series that attempts to observe a clear signature of the Galactic bar/bulge using kinematic observations of the bulge stellar populations in low foreground extinction windows. We report on the detection of ∼100 000 new proper motions in four fields covering the far side of the Galactic bar/bulge, at negative longitudes. Our proper motions have been obtained using observations from the Advance Camera for Surveys (ACS), on board of the Hubble Space Telescope (HST), with a time-baseline of 8–9 years, which has produced accuracies better than 0.5 mas yr−1 for a significant fraction of the stellar populations with F814W < 23 mag. Interestingly, as shown in previous works, the Hess diagrams show a strikingly similar proper motion distribution to fields closer to the Galactic center and consistent with an old stellar population. The observed kinematics point to a significant bulge rotation, which seems to predominate even in fields as far as l ≃ −8°, and is also reflected in the changes of the velocity ellipsoid in the l, b plane as a function of distance.

Galactic science encompasses a wide range of subjects in the study of the Milky Way and Magellanic Clouds, from young stellar objects to X-ray binaries. Mapping these populations, and exploring transient phenomena within them, are among the primary science goals of the Vera C. Rubin Observatory's Legacy Survey of Space and Time. While early versions of the survey strategy dedicated relatively few visits to the Galactic Plane region, more recent strategies under consideration envision a higher cadence within selected regions of high scientific interest. The range of galactic science presents a challenge in evaluating which strategies deliver the highest scientific returns. Here we present metrics designed to evaluate Rubin survey strategy simulations, based on the cadence of observations they deliver within regions of interest to different topics in galactic science, using variability categories defined by timescale. We also compare the fractions of exposures obtained in each filter with those recommended for the different science goals. We find that the baseline_v2.x simulations deliver observations of the high-priority regions at sufficiently high cadence to reliably detect variability on timescales >10 days or more. Follow-up observations may be necessary to properly characterize variability, especially transients, on shorter timescales. Combining the regions of interest for all the science cases considered, we identify those areas of the Galactic Plane and Magellanic Clouds of highest priority. We recommend that these refined survey footprints be used in future simulations to explore rolling cadence scenarios, and to optimize the sequence of observations in different bandpasses.

We present wide-field, high resolution maps of the color excess for 14 globular clusters toward the Southern Galactic bulge. The maps were derived using Gaia EDR3 astrometry and stellar photometry from the Blanco DECam Bulge Survey, which is a deep, wide-field ugriz Y photometric survey of the southern Galactic bulge. Comparisons with WISE 12 μ m images of thermal continuum emission demonstrate that the maps presented here trace interstellar extinction by dust down to 5″ scales. We use the reddening-corrected photometry of proper motion-selected cluster stars to build color–magnitude diagrams for the target globular clusters, which show residual broadening in excess of that expected from the photometric errors alone. This residual broadening is likely to be driven by star-to-star elemental abundance variations.

We advocate for a Galactic center (GC) field to be added to the Galactic Bulge Time Domain Survey (GBTDS). The new field would yield high-cadence photometric and astrometric measurements of an unprecedented ${\sim}$3.3 million stars toward the GC. This would enable a wide range of science cases, such as finding star-compact object binaries that may ultimately merge as LISA-detectable gravitational wave sources, constraining the mass function of stars and compact objects in different environments, detecting populations of microlensing and transiting exoplanets, studying stellar flares and variability in young and old stars, and monitoring accretion onto the central supermassive black hole. In addition, high-precision proper motions and parallaxes would open a new window into the large-scale dynamics of stellar populations at the GC, yielding insights into the formation and evolution of galactic nuclei and their co-evolution with the growth of the supermassive black hole. We discuss the possible trade-offs between the notional GBTDS and the addition of a GC field with either an optimal or minimal cadence. Ultimately, the addition of a GC field to the GBTDS would dramatically increase the science return of Roman and provide a legacy dataset to study the mid-plane and innermost regions of our Galaxy.

SMC X-1 is a high-mass X-ray binary showing superorbital modulation with an unstable period. Previous monitoring shows three excursion events in 1996--1998, 2005--2007, and 2014--2016. The superorbital period drifts from >60 days to <40 days and then evolves back during an excursion. Here we report a new excursion event of SMC X-1 in 2020--2021, indicating that the superorbital modulation has an unpredictable, chaotic nature. We trace the spin-period evolution and find that the spin-up rate accelerated one year before the onset of this new excursion, which suggests a possible inside-out process connecting the spin-up acceleration and the superorbital excursion. This results in a deviation of the spin period residual, similar to the behaviour of the first excursion in 1996--1998. In further analysis of the pulse profile evolution, we find that the pulsed fraction shows a long-term evolution and may be connected to the superorbital excursion. These discoveries deepen the mystery of SMC X-1 because they cannot be solely interpreted by the warped disc model. Upcoming pointed observations and theoretical studies may improve our understanding of the detailed accretion mechanisms taking place.

The Vera C. Rubin Observatory will undertake the Legacy Survey of Space and Time, providing an unprecedented, volume-limited catalog of star clusters in the Southern Sky, including Galactic and extragalactic star clusters. The Star Clusters subgroup of the Stars, Milky Way and Local Volume Working Group has identified key areas where Rubin Observatory will enable significant progress in star cluster research. This roadmap represents our science cases and preparation for studies of all kinds of star clusters from the Milky Way out to distances of tens of megaparsecs.

Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, can be explored by taking advantage of the multiband filters option available in Rubin LSST, in particular the u, g, r, i filters that enable us to discriminate between photospheric stellar properties and accretion signatures. The homogeneity and depth of sky coverage that will be achieved with LSST will provide us with a unique opportunity to characterize the time evolution of disk accretion as a function of age and varying environmental conditions (e.g., field crowdedness, massive neighbors, metallicity) by targeting different star-forming regions. In this contribution to the Rubin LSST Survey Strategy Optimization Focus Issue, we discuss how implementing a dense observing cadence to explore short-term variability in YSOs represents a key complementary effort to the Wide–Fast–Deep observing mode that will be used to survey the sky over the full duration of the main survey (≈10 yr). The combination of these two modes will be vital to investigate the connection between the inner-disk dynamics and longer-term eruptive variability behaviors, such as those observed on EX Lupi–type objects.

With their emission mainly coming from a relativistic jet pointing toward us, blazars are fundamental sources for studying extragalactic jets and their central engines, consisting of supermassive black holes fed by accretion disks. They are also candidate sources of high-energy neutrinos and cosmic rays. Because of the jet orientation, the nonthermal blazar emission is Doppler beamed; its variability is unpredictable, and it occurs on timescales from less than 1 hr to years. Comprehension of the diverse mechanisms producing the flux and spectral changes requires well-sampled multiband light curves over long time periods. In particular, outbursts are the best test bench for shedding light on the underlying physics, especially when studied in a multiwavelength context. The Vera C. Rubin Legacy Survey of Space and Time (Rubin-LSST) will monitor the southern sky for 10 yr in six photometric bands, offering a formidable tool for studying blazar variability features in a statistical way. The alert system will allow us to trigger follow-up observations of outstanding events, especially at high (keV-to-GeV) and very high (TeV) energies. We here examine the simulated Rubin-LSST survey strategies with the aim of understanding which cadences are more suitable for blazar variability science. Our metrics include light curve and color sampling. We also investigate the problem of saturation, which will affect the brightest and many flaring sources, and will have a detrimental impact on follow-up observations.

ABSTRACT We use an N-body+smoothed particle hydrodynamics simulation of an isolated barred galaxy to study the age dependence of bulge longitudinal proper motion (μl) rotation curves. We show that close to the minor axis (|l| ∼ 0○) the relatively young stars rotate more rapidly than the old stars, as found by Hubble Space Telescope in the Milky Way’s (MW's) bulge. This behaviour would be expected also if the MW were unbarred. At larger |l| a different behaviour emerges. Because younger stars trace a strong bar, their galactocentric radial motions dominate their μl at |l| ∼ 6○, leading to a reversal in the sign of $\langle {\rm{\mu _l}}\rangle $. This results in a rotation curve with forbidden velocities (negative $\langle {\rm{\mu _l}}\rangle $ at positive longitudes, and positive $\langle {\rm{\mu _l}}\rangle $ at negative longitudes). The old stars, instead, trace a much weaker bar and thus their kinematics are more axisymmetric, resulting in no forbidden velocities. We develop metrics of the difference in the $\langle {\rm{\mu _l}}\rangle $ rotation curves of young and old stars, and forbidden velocities. We use these to predict the locations where rotation curve reversals can be observed by HST and the Vera Rubin Telescope. Such measurements would represent support for the amplitude of the bar being a continuous function of age, as predicted by kinematic fractionation, in which the bar strength variations are produced purely by differences in the random motions of stellar populations at bar formation.