Preparing to Discover the Unknown with Rubin LSST: Time Domain

Perhaps the most exciting promise of the Rubin Observatory Legacy Survey of Space and Time (LSST) is its capability to discover phenomena never before seen or predicted: true astrophysical novelties; but the ability of LSST to make these discoveries will depend on the survey strategy. Evaluating candidate strategies for true novelties is a challenge both practically and conceptually. Unlike traditional astrophysical tracers like supernovae or exoplanets, for anomalous objects, the template signal is by definition unknown. We approach this problem by assessing survey completeness in a phase space defined by object color and flux (and their evolution), and considering the volume explored by integrating metrics within this space with the observation depth, survey footprint, and stellar density. With these metrics, we explore recent simulations of the Rubin LSST observing strategy across the entire observed spatial footprint and in specific Local Volume regions: the Galactic Plane and Magellanic Clouds. Under our metrics, observing strategies with greater diversity of exposures and time gaps tend to be more sensitive to genuinely new transients, particularly over time-gap ranges left relatively unexplored by previous surveys. To assist the community, we have made all of the tools developed publicly available. While here we focus on transients, an extension of the scheme to include proper motions and the detection of associations or populations of interest will be communicated in Paper II of this series. This paper was written with the support of the Vera C. Rubin LSST Transients and Variable Stars and Stars, Milky Way, Local Volume Science Collaborations.

Small Protoplanetary Disks in the Orion Nebula Cluster and OMC1 with ALMA

The Orion Nebula Cluster (ONC) is the nearest dense star-forming region at ∼400 pc away, making it an ideal target to study the impact of high stellar density and proximity to massive stars (the Trapezium) on protoplanetary disk evolution. The OMC1 molecular cloud is a region of high extinction situated behind the Trapezium in which actively forming stars are shielded from the Trapezium’s strong radiation. In this work, we survey disks at high resolution with Atacama Large Millimeter/submillimeter Array at three wavelengths with resolutions of 0.″095 (3 mm; Band 3), 0.″048 (1.3 mm; Band 6), and 0.″030 (0.85 mm; Band 7) centered on radio Source I. We detect 127 sources, including 15 new sources that have not previously been detected at any wavelength. 72 sources are spatially resolved at 3 mm, with sizes from ∼8–100 au. We classify 76 infrared-detected sources as foreground ONC disks and the remainder as embedded OMC1 disks. The two samples have similar disk sizes, but the OMC1 sources have a dense and centrally concentrated spatial distribution, indicating they may constitute a spatially distinct subcluster. We find smaller disk sizes and a lack of large (>75 au) disks in both our samples compared to other nearby star-forming regions, indicating that environmental disk truncation processes are significant. While photoevaporation from nearby massive Trapezium stars may account for the smaller disks in the ONC, the embedded sources in OMC1 are hidden from this radiation and thus must truncated by some other mechanism, possibly dynamical truncation or accretion-driven contraction.

Characterizing the target selection pipeline for the Dark Energy Spectroscopic Instrument Bright Galaxy Survey

ABSTRACT We present the steps taken to produce a reliable and complete input galaxy catalogue for the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) using the photometric Legacy Survey DR8 DECam. We analyse some of the main issues faced in the selection of targets for the DESI BGS, such as star–galaxy separation, contamination by fragmented stars and bright galaxies. Our pipeline utilizes a new way to select BGS galaxies using Gaia photometry and we implement geometrical and photometric masks that reduce the number of spurious objects. The resulting catalogue is cross-matched with the Galaxy And Mass Assembly (GAMA) survey to assess the completeness of the galaxy catalogue and the performance of the target selection. We also validate the clustering of the sources in our BGS catalogue by comparing with mock catalogues and the Sloan Digital Sky Survey (SDSS) data. Finally, the robustness of the BGS selection criteria is assessed by quantifying the dependence of the target galaxy density on imaging and other properties. The largest systematic correlation we find is a 7 per cent suppression of the target density in regions of high stellar density.

Signatures of tidal disruption in the Milky Way globular cluster NGC 6981 (M72)

We study the outer regions of the Milky Way globular cluster NGC 6981 based on publicly available BV photometry and new Dark Energy Camera (DECam) observations, both of which reach nearly 4 mag below the cluster main sequence (MS) turnoff. While the BV data sets reveal the present of extra-tidal features around the cluster, the much larger field of view of the DECam observations allowed us to identify some other tidal features, which extend from the cluster toward the opposite direction to the Milky Way center. Such structural features of clusters arise from stellar density maps built using MS stars, following a cleaning of the cluster color-magnitude diagram to remove the contamination of field stars. We also performed N-body simulations in order to help us to understand the spatial distribution of the extra-tidal debris. The outcomes reveal the presence of long trailing and leading tails that are mostly parallel to the direction of the cluster velocity vector. We find that the cluster loses most of its mass by tidal disruption during its perigalactic passages, each of which lasted nearly 20 Myr. Hence, a decrease in the density of escaping stars near the cluster is expected from our N-body simulations, which, in turn, means that stronger extra-tidal features could be found by exploring much larger areas around NGC 6891.

Evidence of a population of dark subhaloes from Gaia and Pan-STARRS observations of the GD-1 stream

ABSTRACT New data from the Gaia satellite, when combined with accurate photometry from the Pan-STARRS survey, allow us to accurately estimate the properties of the GD-1 stream. Here, we analyse the stellar density variations in the GD-1 stream and show that they cannot be due to known baryonic structures such as giant molecular clouds, globular clusters, or the Milky Way’s bar or spiral arms. A joint analysis of the GD-1 and Pal 5 streams instead requires a population of dark substructures with masses ≈107–$10^9 \ \rm {M}_{\odot }$. We infer a total abundance of dark subhaloes normalized to standard cold dark matter $n_{\rm sub}/n_{\rm sub, CDM} = 0.4 ^{+0.3}_{-0.2}$ (68 per cent), which corresponds to a mass fraction contained in the subhaloes $f_{\rm {sub}} = 0.14 ^{+0.11}_{-0.07} {{\ \rm per\ cent}}$, compatible with the predictions of hydrodynamical simulation of cold dark matter with baryons.

Effects of stellar density on the photoevaporation of circumstellar discs

ABSTRACT Circumstellar discs are the precursors of planetary systems and develop shortly after their host star has formed. In their early stages, these discs are immersed in an environment rich in gas and neighbouring stars, which can be hostile for their survival. There are several environmental processes that affect the evolution of circumstellar discs, and external photoevaporation is arguably one of the most important ones. Theoretical and observational evidence point to circumstellar discs losing mass quickly when in the vicinity of massive, bright stars. In this work, we simulate circumstellar discs in clustered environments in a range of stellar densities, where the photoevaporation mass-loss process is resolved simultaneously with the stellar dynamics, stellar evolution, and the viscous evolution of the discs. Our results indicate that external photoevaporation is efficient in depleting disc masses and that the degree of its effect is related to stellar density. We find that a local stellar density lower than 100 stars pc−2 is necessary for discs massive enough to form planets to survive for 2.0 Myr. There is an order of magnitude difference in the disc masses in regions of projected density 100 versus 104 stars pc−2. We compare our results to observations of the Lupus clouds, the Orion Nebula Cluster, the Orion Molecular Cloud-2, Taurus, and NGC 2024, and find that the trends observed between region density and disc masses are similar to those in our simulations.

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

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’s 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, cleaned, field star 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 mag below the MS turnoff of the cluster 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.5° 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.5° from the cluster center.

Encounters involving planetary systems in birth environments: the significant role of binaries

ABSTRACT Most stars form in a clustered environment. Both single and binary stars will sometimes encounter planetary systems in such crowded environments. Encounter rates for binaries may be larger than for single stars, even for binary fractions as low as 10–20 per cent. In this work, we investigate scatterings between a Sun–Jupiter pair and both binary and single stars as in young clusters. We first perform a set of simulations of encounters involving wide ranges of binaries and single stars, finding that wider binaries have larger cross-sections for the planet’s ejection. Secondly, we consider such scatterings in a realistic population, drawing parameters for the binaries and single stars from the observed population. The scattering outcomes are diverse, including ejection, capture/exchange, and collision. The binaries are more effective than single stars by a factor of several or more in causing the planet’s ejection and collision. Hence, in a cluster, as long as the binary fraction is larger than about 10 per cent, the binaries will dominate the scatterings in terms of these two outcomes. For an open cluster of a stellar density 50 pc−3, a lifetime 100 Myr, and a binary fraction 0.5, we estimate that Jupiters of the order of 1 per cent are ejected, 0.1 per cent collide with a star, 0.1 per cent change ownership, and 10 per cent of the Sun–Jupiter pairs acquire a stellar companion during scatterings. These companions are typically thousands of au distant and in half of the cases (so 5 per cent of all Sun–Jupiter pairs), they can excite the planet’s orbit through Kozai–Lidov mechanism before being stripped by later encounters. Our result suggests that the Solar system may have once had a companion in its birth cluster.

FEDReD

Context. While Gaia enables us to probe the extended local neighbourhood in great detail, the thin disc structure at larger distances remains sparsely explored. Aims. We aim here to build a non-parametric 3D model of the thin disc structures handling both the extinction and the stellar density simultaneously. Methods. We developed a Bayesian deconvolution method in two dimensions: extinction and distance. It uses a reference catalogue whose completeness information defines the selection function. It is designed so that any complementary information from other catalogues can be added. It has also been designed to be robust to outliers, which are frequent in crowded fields, and differential extinction. The prior information is designed to be minimal: only a reference H-R diagram. We derived for this an empirical H-R diagram of the thin disc using Gaia DR2 data, but synthetic isochrone-based H-R diagrams can also be used. Results. We validated the method on simulations and real fields using 2MASS and UKIDSS data complemented by Gaia DR2 photometry and parallaxes. We detail the results of two test fields: a 2MASS field centred around the NGC 4815 open cluster, which shows an over-density of both extinction and stellar density at the cluster distance, and a UKIDSS field at l = 10° where we recover the position of the Galactic bar.

A robust two-parameter description of the stellar profile of elliptical galaxies

Context. The stellar density profile of a galaxy is typically summarised with two numbers: the total stellar mass and half-light radius. The total mass of a galaxy, however, is not a well-defined quantity, due to the finite depth of photometric observations and the arbitrariness of the distinction between galaxy and diffuse intra-group light. This limits our ability to make accurate comparisons between models and observations. Aims. I wish to provide a more robust two-parameter description of the stellar density distribution of elliptical galaxies, in terms of quantities that can be measured unambiguously. Methods. I propose using the stellar mass enclosed within 10 kpc in projection, M*,10, and the mass-weighted stellar density slope within the same aperture, Γ*,10, for this purpose. I measured the distribution in M*,10 and Γ*,10 of a sample of elliptical galaxies from the Sloan Digital Sky Survey and the Galaxy And Mass Assembly survey, using photometry from the Hyper Suprime-Cam survey. I measured, at fixed (M*,10, Γ*,10), what the spread is in the galaxy surface brightness profile and central stellar velocity dispersion within the sample. As a first application, I then compared the observed M*,10 − Γ*,10 relation of elliptical galaxies with that of similarly selected galaxies in the EAGLE REFERENCE simulation. Results. The pair of values of (M*,10, Γ*,10) can be used to predict the stellar density profile in the inner 10 kpc of a galaxy with better than 20% accuracy. Similarly, M*,10 and Γ*,10 can be combined to obtain a proxy for stellar velocity dispersion that is at least as good as the stellar mass fundamental plane. The average stellar density slope of EAGLE elliptical galaxies matches that of observed ones at M*,10 = 1011 M⊙ well, but the EAGLE M*,10 − Γ*,10 relation is shallower and has a larger intrinsic scatter compared to observations. Conclusions. This new parameterisation of the stellar density profile of massive elliptical galaxies provides a more robust way of comparing results from different photometric surveys and from hydrodynamical simulations, with respect to a description based on total stellar mass and half-light radius.