Search

We obtain a quantitative star formation history (SFH) of a shell-like structure ('shell') located in the northeastern part of the Small Magellanic Cloud (SMC). We use the Survey of the MAgellanic Stellar History (SMASH) to derive colour-magnitude diagrams (CMDs), reaching below the oldest main-sequence turnoff, from which we compute the SFHs with CMD fitting techniques. We present, for the first time, a novel technique that uses red clump (RC) stars from the CMDs to assess and account for the SMC's line-of-sight depth effect present during the SFH derivation. We find that accounting for this effect recovers a more accurate SFH. We quantify a 7 kpc line-of-sight depth present in the CMDs, in good agreement with depth estimates from RC stars in the northeastern SMC. By isolating the stellar content of the northeastern shell and incorporating the line-of-sight depth into our calculations, we obtain an unprecedentedly detailed SFH. We find that the northeastern shell is primarily composed of stars younger than 500 Myrs, with significant star formation enhancements around 250 Myr and 450 Myr. These young stars are the main contributors to the shell's structure. We show synchronicity between the northeastern shell's SFH with the Large Magellanic Cloud's (LMC) northern arm, which we attribute to the interaction history of the SMC with the LMC and the Milky Way (MW) over the past 500 Myr. Our results highlight the complex interplay of ram pressure stripping and the influence of the MW's circumgalactic medium in shaping the SMC's northeastern shell., Comment: 17 pages, 13 figures. Accepted to MNRAS for publication

We show with Gaia XP spectroscopy that extremely metal-rich stars in the Milky Way (EMR; $[M/H]_{XP} > 0.5$) - but only those - are largely confined to a tight "knot" at the center of the Galaxy. This EMR knot is round in projection, has a fairly abrupt edge near $\sim 1.5$kpc, and is a dynamically hot system. This central knot also contains very metal-rich (VMR; $+0.2\le [M/H]_{XP} \le +0.4$) stars. However, in contrast to EMR stars, the bulk of VMR stars form an extended, highly flattened distribution in the inner Galaxy ($R_{\mathrm{GC}}\lesssim 5$ kpc). We draw on TNG50 simulations of Milky Way analogs for context and find that compact, metal-rich knots confined to $<1.5$kpc are a universal feature. In typical simulated analogs, the top 5-10% most metal-rich stars are confined to a central knot; however, in our Milky Way data this fraction is only 0.1%. Dust-penetrating wide-area near-infrared spectroscopy, such as SDSS-V, will be needed for a rigorous estimate of the fraction of stars in the Galactic EMR knot. Why in our Milky Way only EMR giants are confined to such a central knot remains to be explained. Remarkably, the central few kiloparsecs of the Milky Way harbor both the highest concentration of metal-poor stars (the `poor old heart') and almost all EMR stars. This highlights the stellar population diversity at the bottom of galactic potential wells., Comment: 11 pages, 7 figures, submitted to ApJ

We determine radial- and age-abundance gradients of the Small Magellanic Cloud (SMC) using spectra of 2,062 red giant branch (RGB) field stars observed by SDSS-IV / APOGEE-2S. With coverage out to $\sim$9 kpc in the SMC, these data taken with the high resolution ($R \sim 22,500$) APOGEE $H$-band spectrograph afford the opportunity to measure extensive radial gradients for as many as 24 abundance ratios. The SMC is found to have an overall metallicity gradient of $-$0.0546 $\pm$ 0.0043 dex/kpc. Ages are calculated for every star to explore the evolution of the different abundance gradients. As a function of age, many of the gradients show a feature 3.66--5.58 Gyr ago, which is especially prominent in the [X/H] gradients. Initially many gradients flatten until about $\sim$5.58 Gyr ago, but then steepen in more recent times. We previously detected similar evolutionary patterns in the Large Magellanic Cloud (LMC) which are attributed to a recent interaction between the LMC and SMC. It is inferred that the feature in the SMC gradients was caused by the same interaction. The age-[X/Fe] trends, which track average [X/Fe] over time, are flat, demonstrating a slow enrichment history for the SMC. When comparing the SMC gradients to the LMC and MW, normalized to disk scale length ($R_\text{d}$), the [X/Fe] and [X/Mg] gradients are similar, but there is a dichotomy between the dwarfs and the Milky Way (MW) for the [X/H] gradients. The median MW [X/H] gradient around $-$0.125 dex/$R_\text{d}$ whilst the Clouds have gradients of about $-$0.075 dex/$R_\text{d}$., Comment: 27 pages, 22 figures, and 11 tables

We present the abundance gradients of the Large Magellanic Cloud (LMC) for 25 elemental abundance ratios and their respective temporal evolution as well as age-[X/Fe] trends using 6130 LMC field red giant branch (RGB) stars observed by SDSS-IV / APOGEE-2S. APOGEE is a high resolution ($R$ $\sim$22,500) $H$-band spectroscopic survey that gathered data on the LMC with broad radial and azimuthal coverage out to $\sim$10\degr. The calculated overall metallicity gradient of the LMC with no age binning is $-$0.0380 $\pm$ 0.0022 dex/kpc. We also find that many of the abundance gradients show a U-shaped trend as functions of age. This trend is marked by a flattening of the gradient but then a general steepening at more recent times. The extreme point at which all these gradients (with the U-shaped trend) begin to steepen is $\gtrsim$2 Gyr ago. In addition, some of the age-[X/Fe] trends show an increase starting a few Gyr before the extreme point in the gradient evolutions. A subset of the age-[X/Fe] trends also show maxima concurrent with the gradients' extreme points, further pinpointing a major event in the history of the LMC $\sim$2 Gyr ago. This time frame is consistent with a previously proposed interaction between the Magellanic Clouds suggesting that this is most likely the cause of the distinct trend in the gradients and age-[X/Fe] trends., Comment: 25 pages, 19 figures, and 10 tables

The Magellanic Cloud system represents a unique laboratory for study of both interacting dwarf galaxies and the ongoing process of the formation of the Milky Way and its halo. We focus on one aspect of this complex, 3 body interaction - the dynamical perturbation of the Small Magellanic Cloud (SMC) by the Large Magellanic Cloud (LMC), and specifically potential tidal effects on the SMC's eastern side. Using Gaia astrometry and the precise radial velocities and multi-element chemical abundances from APOGEE-2 DR17, we explore the well-known distance bimodality on the eastern side of the SMC. Through estimated stellar distances, proper motions, and radial velocities, we characterize the kinematics of the two populations in the bimodality and compare their properties with those of SMC populations elsewhere. Moreover, while all regions explored by APOGEE seem to show a single chemical enrichment history, the metallicity distribution function (MDF), of the "far" stars on the eastern periphery of the SMC is found to resemble that for the more metal-poor fields of the western periphery, whereas the MDF for the "near" stars on the eastern periphery resembles that for stars in the SMC center. The closer eastern periphery stars also show radial velocities (corrected for SMC rotation and bulk motion) that are, on average, approaching us relative to all other SMC populations sampled. We interpret these trends as evidence that the near stars on the eastern side of the SMC represent material pulled out of the central SMC as part of its tidal interaction with the LMC.

I report the discovery of a stellar stream (Sutlej) using Gaia DR3 proper motions and XP metallicities located ~15 degrees north of the Small Magellanic Cloud (SMC). The stream is composed of two parallel linear components ("branches") approximately ~8 x 0.6 degrees in size and separated by 2.5 degrees. The stars have a mean proper motion of (pmra,pmdec)=(+0.08 mas/yr,-1.41 mas/yr) which is quite similar to the proper motion of stars on the western side of the SMC. The color magnitude diagram of the stream stars has a clear red giant branch, horizontal branch, and main sequence turnoff that is well-matched by a PARSEC isochrone of 10 Gyr, [Fe/H]=-1.8 at 32 kpc and a total stellar mass of ~33,000 Msun. The stream is spread out over an area of 9.6 square degrees and has a surface brightness of 32.5 mag/arcsec^2. The metallicity of the stream stars from Gaia XP spectra extend over -2.5 < [M/H] < -1.0 with a median of [M/H]=-1.8. The tangential velocity of the stream stars is 214 km/s compared to the values of 448 km/s for the Large Magellanic Cloud and 428 km/s for the SMC. While the radial velocity of the stream is not yet known, a comparison of the space velocities using a range of assumed radial velocities, shows that the stream is unlikely to be associated with the Magellanic Clouds. The tangential velocity vector is misaligned with the stream by ~25 degrees which might indicate an important gravitational influence from the nearby Magellanic Clouds., Comment: 8 pages, 10 figures, 1 table. Submitted to MNRAS

In this paper, we analyse the metallicity structure of the Magellanic Clouds using parameters derived from the Gaia DR3 low-resolution XP spectra, astrometry and photometry. We find that the qualitative behavior of the radial metallicity gradients in the LMC and SMC are quite similar, with both of them having a metallicity plateau at intermediate radii and a second at larger radii. The LMC has a first metallicity plateau at [Fe/H]$\approx$-0.8 for 3$-$7$^{\circ}$, while the SMC has one at [Fe/H]$\approx$-1.1 at 3$-$5$^{\circ}$. The outer LMC periphery has a fairly constant metallicity of [Fe/H]$\approx$-1.0 (10$-$18$^{\circ}$), while the outer SMC periphery has a value of [Fe/H]$\approx$-1.3 (6$-$10$^{\circ}$). The sharp drop in metallicity in the LMC at $\sim$8$^{\circ}$ and the marked difference in age distributions in these two regions suggests that there were two important evolutionary phases in the LMC. In addition, we find that the Magellanic periphery substructures, likely Magellanic debris, are mostly dominated by LMC material stripped off in old interactions with the SMC. This presents a new picture in contrast with the popular belief that the debris around the Clouds had been mostly stripped off from the SMC due to having a lower mass. We perform a detailed analysis for each known substructure and identify its potential origin based on metallicities and motions with respect to each galaxy., Comment: 11 pages, 7 figures. Accepted to MNRAS for publication. Comments welcome!

We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), $\alpha$-element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low-$\alpha$ disk, high-$\alpha$ disk, and total population independently. The low-$\alpha$ disk exhibits a negative radial metallicity gradient of $-0.06 \pm 0.001$ dex kpc$^{-1}$, which flattens with distance from the midplane. The high-$\alpha$ disk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]-[Fe/H] plane and in the [Mg/Fe]-age relation persist across the entire disk. The age estimates have typical uncertainties of $\sim0.15$ in $\log$(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy., Comment: 41 pages, 32 figures, accepted to ApJ

The Milky Way stellar disk has both a thin and a thick component. The thin disk is composed mostly of younger stars ($\lesssim$8 Gyr) with a lower abundance of $\alpha$ elements, while the thick disk contains predominantly older stars ($\gtrsim$8--12 Gyr) with a higher $\alpha$ abundance, giving rise to an $\alpha$-bimodality most prominent at intermediate metallicities. A proposed explanation for the bimodality is an episode of clumpy star formation, where high-$\alpha$ stars form in massive clumps that appear in the first few Gyrs of the Milky Way's evolution, while low-$\alpha$ stars form throughout the disk and over a longer time span. To better understand the evolution of clumps, we track them and their constituent stars in two clumpy Milky Way simulations that reproduce the $\alpha$-abundance bimodality, one with 10% and the other with 20% supernova feedback efficiency. We investigate the paths that these clumps take in the chemical space ([O/Fe]--[Fe/H]) as well as their mass, star formation rate (SFR), formation location, lifetime, and merger history. The clumps in the simulation with lower feedback last longer on average, with several lasting hundreds of Myr. Some of the clumps do not reach high-$\alpha$, but the ones that do on average had a higher SFR, longer lifetime, greater mass, and form closer to the galactic center than the ones that do not. Most clumps that reach high-$\alpha$ merge with others and eventually spiral into the galactic center, but shed stars along the way to form most of the thick disk component., Comment: 14 pages, 14 figures, accepted for publication in ApJ

As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($\Delta t\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$\sigma$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own., Comment: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

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., Comment: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

Stellar ages are critical for understanding the temporal evolution of a galaxy. We calculate the ages of over 6000 red giant branch stars in the Large Magellanic Cloud (LMC) observed with SDSS-IV / APOGEE-S. Ages are derived using multi-band photometry, spectroscopic parameters (T$_\text{eff}$, $\log{g}$, [Fe/H], and [$\alpha$/Fe]) and stellar isochrones and the assumption that the stars lie in a thin inclined plane to get accurate distances. The isochrone age and extinction are varied until a best match is found for the observed photometry. We perform validation using the APOKASC sample, which has asteroseismic masses and accurate ages, and find that our uncertainties are $\sim$20% and range from $\sim$1$-$3 Gyr for the calculated age values. Here we present the LMC age map as well as the age-radius relation and an accurate age-metallicity relation (AMR). The age map and age-radius relation reveal that recent star formation in the galaxy was more centrally located and that there is a slight dichotomy between the north and south with the northern fields being slightly younger. The northern fields that cover a known spiral arm have median ages of $\gtrsim$ 2 Gyr, which is the time when an interaction with the SMC is suggested to have happened. The AMR is mostly flat especially for older ages although recently (about 2.0-2.5 Gyr ago) there is an increase in the median [Fe/H]. Based on the time frame, this might also be attributed to the close interaction between the LMC and SMC., Comment: 21 pages, 22 figures, submitted to MNRAS

We present initial results from our JWST NIRSpec program to study the $\alpha$-abundances in the M31 disk. The Milky Way has two chemically-defined disks, the low-$\alpha$ and high-$\alpha$ disks, which are closely related to the thin and thick disks, respectively. The origin of the two populations and the $\alpha$-bimodality between them is not entirely clear, although there are now several models that can reproduce the observed features. To help constrain the models and discern the origin, we have undertaken a study of the chemical abundances of the M31 disk using JWST NIRSpec, in order to determine whether stars in M31's disk also show an $\alpha$-abundance bimodality. Approximately 100 stars were observed in our single NIRSpec field at a projected distance of 18 kpc from the M31 center. The 1-D extracted spectra have an average signal-to-noise ratio of 85 leading to statistical metallicity precision of 0.016 dex, $\alpha$-abundance precision of 0.012 dex, and a radial velocity precision 8 km/s. The initial results indicate that, in contrast to the Milky Way, there is no $\alpha$-bimodality in the M31 disk, and no low-$\alpha$ sequence. The entire stellar population falls along a single chemical sequence very similar to the MW's high-alpha component which had a high star formation rate. While this is somewhat unexpected, the result is not that surprising based on other studies that found the M31 disk has a larger velocity dispersion than the MW and is dominated by a thick component. M31 has had a more active accretion and merger history than the MW which might explain the chemical differences., Comment: 8 pages, 4 figures, IAU Symposium 377, Early Disk-Galaxy Formation: From JWST to the Milky Way

We present the first detailed chemical analysis from APOGEE-2S observations of stars in six regions of recently discovered substructures in the outskirts of the Magellanic Clouds extending to 20 degrees from the LMC center. We also present, for the first time, the metallicity and alpha-abundance radial gradients of the LMC and SMC out to 11 degrees and 6 degrees, respectively. Our chemical tagging includes 13 species including light, alpha, and Fe-peak elements. We find that the abundances of all of these chemical elements in stars populating two regions in the northern periphery - along the northern "stream"-like feature - show good agreement with the chemical patterns of the LMC, and thus likely have an LMC origin. For substructures located in the southern periphery of the LMC, we find more complex chemical and kinematical signatures, indicative of a mix of LMC-like and SMC-like populations. However, the southern region closest to the LMC shows better agreement with the LMC, whereas that closest to the SMC shows a much better agreement with the SMC chemical pattern. When combining this information with 3-D kinematical information for these stars, we conclude that the southern region closest to the LMC has likely an LMC origin, whereas that closest to the SMC has an SMC origin, and the other two southern regions have a mix of LMC and SMC origins. Our results add to the evidence that the southern substructures of the LMC periphery are the product of close interactions between the LMC and SMC, and thus likely hold important clues that can constrain models of their detailed dynamical histories., Comment: 21 pages, 15 figures, 9 tables; Accepted for publication in A&A

Despite the universe containing primordial thorium (Th) of sufficient abundance to appear in stellar spectra, detection of Th has to date been tentative and based on just a few weak and blended lines. Here, we present convincing evidence not only for the first Th detection in a magnetic chemically peculiar Ap star but also for the first detection of Th III in a stellar spectrum. CPD-62 2717 was initially recognized as a highly-magnetized Ap star thanks to resolved magnetically split lines captured in $H$-band spectra from the SDSS/APOGEE survey. The star was subsequently pinpointed as extraordinarily peculiar when careful inspection of the $H$-band line content revealed the presence of five lines of Th III, none of which are detected in the other $\sim1500$ APOGEE-observed Ap stars. Follow-up with the VLT+UVES confirmed a similarly peculiar optical spectrum featuring dozens of Th III lines, among other peculiarities. Unlike past claims of Th detection, and owing to high-resolution observations of the strong ($\sim$8$-$12$\,$kG) magnetic field of CPD-62 2717, the detection of Th III can in this case be supported by matches between the observed and theoretical magnetic splitting patterns. Comparison of CPD-62 2717 to stars for which Th overabundances have been previously reported (e.g., Przybylski's Star) indicate that only for CPD-62 2717 is the Th detection certain. Along with the focus on Th III, we use time series measurements of the magnetic field modulus to constrain the rotation period of CPD-62 2717 to $\sim$4.8 years, thus establishing it as a new example of a super-slowly-rotating Ap star.

In Paper I we showed that clumps in high-redshift galaxies, having a high star formation rate density (\Sigma_SFR), produce disks with two tracks in the [Fe/H]-[\alpha/Fe] chemical space, similar to that of the Milky Way's (MW's) thin + thick disks. Here we investigate the effect of clumps on the bulge's chemistry. The chemistry of the MW's bulge is comprised of a single track with two density peaks separated by a trough. We show that the bulge chemistry of an N-body + smoothed particle hydrodynamics clumpy simulation also has a single track. Star formation within the bulge is itself in the high-\Sigma_SFR clumpy mode, which ensures that the bulge's chemical track follows that of the thick disk at low [Fe/H] and then extends to high [Fe/H], where it peaks. The peak at low metallicity instead is comprised of a mixture of in-situ stars and stars accreted via clumps. As a result, the trough between the peaks occurs at the end of the thick disk track. We find that the high-metallicity peak dominates near the mid-plane and declines in relative importance with height, as in the MW. The bulge is already rapidly rotating by the end of the clump epoch, with higher rotation at low [\alpha/Fe]. Thus clumpy star formation is able to simultaneously explain the chemodynamic trends of the MW's bulge, thin + thick disks and the Splash., Comment: 25 pages, in press at ApJ

The Astro Data Lab is preparing to host large spectroscopic datasets such as a copy of the Dark Energy Spectroscopic Instrument (DESI) survey, which is projected to include approximately 40 million spectra of galaxies and quasars as well as over 10 million spectra of stars by 2026. Currently, we serve DR16 spectra from the Sloan Digital Sky Survey (SDSS), including Baryon Oscillation Spectroscopic Survey (BOSS), and Extended BOSS (eBOSS) spectra. A spectral viewer tool allows users to visually and interactively inspect spectra. Given the large size of these spectroscopic datasets, a typical use case might consist of a selection or query for a subset of objects of interest (e.g., a subsample of stars or galaxies or quasars), followed by visual inspection of the selected spectra. It is anticipated that in some cases, users will want to go through a long list of spectra (e.g., thousands) quickly while looking for specific features. This document contains a description of the requirements for such a spectral viewer tool to be incorporated within the Astro Data Lab environment at NSF's NOIRLab. For each object, the spectral viewer will display the observed spectrum and, if available, the noise spectrum, sky spectrum, and best-fit template spectrum. Users will be able to control the display interactively after they launch the tool as part of their Data Lab workflow. The primary objective will be to support the visualization of spectroscopic datasets hosted at the Astro Data Lab but this requirements document could be a useful reference or inspiration for other applications and/or other datasets in the astronomy community., Comment: Spectral Viewer requirements document for astrophysical software, 27 pages

The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM), including input catalogs and selection functions for their numerous scientific objectives. We describe the production of the targeting databases and their calibration- and scientifically-focused components. DR18 also includes ~25,000 new SDSS spectra and supplemental information for X-ray sources identified by eROSITA in its eFEDS field. We present updates to some of the SDSS software pipelines and preview changes anticipated for DR19. We also describe three value-added catalogs (VACs) based on SDSS-IV data that have been published since DR17, and one VAC based on the SDSS-V data in the eFEDS field., Comment: Accepted to ApJS

We report the discovery of a new example of the rare class of highly magnetized, rapidly rotating, helium enhanced, early B stars that produce anomalously wide hydrogen emission due to a centrifugal magnetosphere (CM). The star is Trumpler 16-26, a B1.5 V member of the Trumpler 16 open cluster. A CM was initially suspected based on hydrogen Brackett series emission observed in SDSS/APOGEE $H$-band spectra. Similar to the other stars of this type, the emission was highly variable and at all times remarkable due to the extreme velocity separations of the double peaks (up to 1300 km s$^{-1}$.) Another clue lay in the TESS lightcurve, which shows two irregular eclipses per cycle when phased with the likely 0.9718115 day rotation period, similar to the behavior of the well known CM host star $\sigma$ Ori E. To confirm a strong magnetic field and rotation-phase-locked variability, we initiated a follow-up campaign consisting of optical spectropolarimetry and spectroscopy. The associated data revealed a longitudinal magnetic field varying between $-3.1$ and $+1.6$ kG with the period found from photometry. The optical spectra confirmed rapid rotation ($v \sin i=195$ km s$^{-1}$), surface helium enhancement, and wide, variable hydrogen emission. Tr16-26 is thus confirmed as the 20$^{\rm th}$ known, the fourth most rapidly rotating, and the faintest CM host star yet discovered. With a projected dipole magnetic field strength of $B_{\rm d}>11$ kG, Tr16-26 is also among the most magnetic CM stars.

We present a map of the total intrinsic reddening across ~90 deg$^{2}$ of the Large Magellanic Cloud (LMC) derived using optical (ugriz) and near-infrared (IR; YJKs) spectral energy distributions (SEDs) of background galaxies. The reddening map is created from a sample of 222,752 early-type galaxies based on the LEPHARE $\chi^{2}$ minimisation SED-fitting routine. We find excellent agreement between the regions of enhanced intrinsic reddening across the central (4x4 deg$^2$) region of the LMC and the morphology of the low-level pervasive dust emission as traced by far-IR emission. In addition, we are able to distinguish smaller, isolated enhancements that are coincident with known star-forming regions and the clustering of young stars observed in morphology maps. The level of reddening associated with the molecular ridge south of 30 Doradus is, however, smaller than in the literature reddening maps. The reduced number of galaxies detected in this region, due to high extinction and crowding, may bias our results towards lower reddening values. Our map is consistent with maps derived from red clump stars and from the analysis of the star formation history across the LMC. This study represents one of the first large-scale categorisations of extragalactic sources behind the LMC and as such we provide the LEPHARE outputs for our full sample of ~2.5 million sources., Comment: 18 pages, 16 figures, accepted for publication in MNRAS

We report the first 3-D kinematical measurements of 88 stars in the direction of several recently discovered substructures in the southern periphery of the Large Magellanic Cloud (LMC) using a combination of Gaia proper motions and radial velocities from the APOGEE-2 survey. More specifically, we explore stars lie in assorted APOGEE-2 pointings in a region of the LMC periphery where various overdensities of stars have previously been identified in maps of stars from Gaia and DECam. By using a model of the LMC disk rotation, we find that a sizeable fraction of the APOGEE-2 stars have extreme space velocities that are distinct from, and not a simple extension of, the LMC disk. Using N-body hydrodynamical simulations of the past dynamical evolution and interaction of the LMC and Small Magellanic Cloud (SMC), we explore whether the extreme velocity stars may be accounted for as tidal debris created in the course of that interaction. We conclude that the combination of LMC and SMC debris produced from their interaction is a promising explanation, although we cannot rule out other possible origins, and that these new data should be used to constrain future simulations of the LMC-SMC interaction. We also conclude that many of the stars in the southern periphery of the LMC lie out of the LMC plane by several kpc. Given that the metallicity of these stars suggest they are likely of Magellanic origin, our results suggest that a wider exploration of the past interaction history of the Magellanic Clouds is needed.

We present analysis of the proper-motion (PM) field of the red clump stars in the Large Magellanic Cloud (LMC) disk using the Gaia Early Data Release 3 catalog. Using a kinematic model based on old stars with 3D velocity measurements, we construct the residual PM field by subtracting the center-of-mass motion and internal rotation motion components. The residual PM field reveals asymmetric patterns, including larger residual PMs in the southern disk. Comparisons between the observed residual PM field with those of five numerical simulations of an LMC analog that is subject to the tidal fields of the Milky Way and the Small Magellanic Cloud (SMC) show that the present-day LMC is not in dynamical equilibrium. We find that both the observed level of disk heating (PM residual root-mean-square of 0.057$\pm$0.002 mas yr$^{-1}$) and kinematic asymmetry are not reproduced by Milky Way tides or if the SMC impact parameter is larger than the size of the LMC disk. This measured level of disk heating provides a novel and important method to validate numerical simulations of the LMC-SMC interaction history. Our results alone put constraints on an impact parameter $\lesssim$10 kpc and impact timing $<$250 Myr. When adopting the impact timing constraint of $\sim$140--160 Myr ago from previous studies, our results suggest that the most recent SMC encounter must have occurred with an impact parameter of $\sim$5 kpc. We also find consistent radial trends in the kinematically- and geometrically-derived disk inclination and line-of-node position angles, indicating a common origin., Comment: 17 pages, 10 figures, resubmitted to AAS Journal

This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys., Comment: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey has obtained high-resolution spectra for thousands of red giant stars distributed among the massive satellite galaxies of the Milky Way (MW): the Large and Small Magellanic Clouds (LMC/SMC), the Sagittarius Dwarf (Sgr), Fornax (Fnx), and the now fully disrupted \emph{Gaia} Sausage/Enceladus (GSE) system. We present and analyze the APOGEE chemical abundance patterns of each galaxy to draw robust conclusions about their star formation histories, by quantifying the relative abundance trends of multiple elements (C, N, O, Mg, Al, Si, Ca, Fe, Ni, and Ce), as well as by fitting chemical evolution models to the [$\alpha$/Fe]-[Fe/H] abundance plane for each galaxy. Results show that the chemical signatures of the starburst in the MCs observed by Nidever et al. in the $\alpha$-element abundances extend to C+N, Al, and Ni, with the major burst in the SMC occurring some 3-4 Gyr before the burst in the LMC. We find that Sgr and Fnx also exhibit chemical abundance patterns suggestive of secondary star formation epochs, but these events were weaker and earlier ($\sim$~5-7 Gyr ago) than those observed in the MCs. There is no chemical evidence of a second starburst in GSE, but this galaxy shows the strongest initial star formation as compared to the other four galaxies. All dwarf galaxies had greater relative contributions of AGB stars to their enrichment than the MW. Comparing and contrasting these chemical patterns highlight the importance of galaxy environment on its chemical evolution., Comment: 37 pages, 21 figures, accepted for publication in ApJ

Despite extensive searches and the relative proximity of solar system objects (SSOS) to Earth, many remain undiscovered and there is still much to learn about their properties and interactions. This work is the first in a series dedicated to detecting and analyzing SSOs in the all-sky NOIRLab Source Catalog (NSC). We search the first data release of the NSC with CANFind, a Computationally Automated NSC tracklet Finder. NSC DR1 contains 34 billion measurements of 2.9 billion unique objects, which CANFind categorizes as belonging to "stationary" (distant stars, galaxies) or moving (SSOs) objects via an iterative clustering method. Detections of stationary bodies for proper motion (mu) less than 2.5"/hr (0.017 degrees/day) are identified and analyzed separately. Remaining detections belonging to hi-mu objects are clustered together over single nights to form "tracklets". Each tracklet contains detections of an individual moving object, and is validated based on spatial linearity and motion through time. Proper motions are then calculated and used to connect tracklets and other unassociated measurements over multiple nights by predicting their locations at common times forming "tracks". This method extracted 527,055 tracklets from NSC DR1 in an area covering 29,971 square degrees of the sky. The data show distinct groups of objects with similar observed mu in ecliptic coordinates, namely Main Belt Asteroids, Jupiter Trojans, and Kuiper Belt Objects. Apparent magnitudes range from 10-25 mag in the ugrizY and VR bands. Color-color diagrams show a bimodality of tracklets between primarily carbonaceous and siliceous groups, supporting prior studies., Comment: 18 pages, 14 figures, 4 tables

APOGEE is a high-resolution (R sim 22,000), near-infrared, multi-epoch, spectroscopic survey of the Milky Way. The second generation of the APOGEE project, APOGEE-2, includes an expansion of the survey to the Southern Hemisphere called APOGEE-2S. This expansion enabled APOGEE to perform a fully panoramic mapping of all the main regions of the Milky Way; in particular, by operating in the H-band, APOGEE is uniquely able to probe the dust-hidden inner regions of the Milky Way that are best accessed from the Southern Hemisphere. In this paper we present the targeting strategy of APOGEE-2S, with special attention to documenting modifications to the original, previously published plan. The motivation for these changes is explained as well as an assessment of their effectiveness in achieving their intended scientific objective. In anticipation of this being the last paper detailing APOGEE targeting, we present an accounting of all such information complete through the end of the APOGEE-2S project; this includes several main survey programs dedicated to exploration of major stellar populations and regions of the Milky Way, as well as a full list of programs contributing to the APOGEE database through allocations of observing time by the Chilean National Time Allocation Committee (CNTAC) and the Carnegie Institution for Science (CIS). This work was presented along with a companion article, R. Beaton et al. (submitted; AAS29028), presenting the final target selection strategy adopted for APOGEE-2 in the Northern Hemisphere., Comment: 43 pages; 5 figures; 6 Tables; 1 Appendix; Submitted to Journal and Under Review; Posting to accompany papers using the SDSS-IV/APOGEE-2 Data Release 17 scheduled for December 2021

LIN 358 and SMC N73 are two symbiotic binaries in the halo of the Small Magellanic Cloud, each composed of a hot white dwarf accreting from a cool giant companion. In this work, we characterize these systems using a combination of SED-fitting to the extant photometric data spanning a broad wavelength range (X-ray/ultraviolet to near-infrared), detailed analysis of the APOGEE spectra for the giant stars, and orbit fitting to high quality radial velocities from the APOGEE database. Using the calculated Roche lobe radius for the giant component and the mass ratio for each system, it is found that LIN 358 is likely undergoing mass transfer via wind Roche lobe overflow while the accretion mechanism for SMC N73 remains uncertain. This work presents the first orbital characterization for both of these systems (yielding periods of >270 and >980 days, respectively, for SMC N73 and LIN 358) and the first global SED fitting for SMC N73. In addition, variability was identified in APOGEE spectra of LIN 358 spanning 17 epochs over two years that may point to a time variable accretion rate as the product of an eccentric orbit., Comment: 19 pages, 12 figures

Complex A is a high-velocity cloud that is traversing through the Galactic halo toward the Milky Way's disk. We combine both new and archival Green Bank Telescope observations to construct a spectroscopically resolved HI~21-cm map of this entire complex at a $17.1\lesssim\log{\left({N_{\rm HI},\,1\sigma}/{\rm cm}^{-2}\right)}\lesssim17.9$ sensitivity for a ${\rm FWHM}=20~{\rm km}\,{\rm s}^{-1}$ line and $\Delta\theta=9.1\,{\rm arcmins}$ or $17\lesssim\Delta d_{\theta}\lesssim30~\rm pc$ spatial resolution. We find that that Complex A is has a Galactic standard of rest frame velocity gradient of $\Delta\rm v_{GSR}/\Delta L=25~{\rm km}\,{\rm s}^{-1}/{\rm kpc}$ along its length, that it is decelerating at a rate of $\langle a\rangle_{\rm GSR}=55~{\rm km}/{\rm yr}^2$, and that it will reach the Galactic plane in $\Delta t\lesssim70~{\rm Myrs}$ if it can survive the journey. We have identify numerous signatures of gas disruption. The elongated and multi-core structure of Complex A indicates that either thermodynamic instabilities or shock-cascade processes have fragmented this stream. We find Rayleigh-Taylor fingers on the low-latitude edge of this HVC; many have been pushed backward by ram-pressure stripping. On the high-latitude side of the complex, Kelvin-Helmholtz instabilities have generated two large wings that extend tangentially off Complex A. The tips of these wings curve slightly forward in the direction of motion and have an elevated \hi\ column density, indicating that these wings are forming Rayleigh-Taylor globules at their tips and that this gas is becoming entangled with unseen vortices in the surrounding coronal gas. These observations provide new insights on the survivability of low-metallicity gas streams that are accreting onto $L_\star$ galaxies., Comment: Figure 5 of the paper is a video in the published that can be found here: https://www.youtube.com/watch?v=YSbJetacspo

The Large and Small Magellanic Clouds (LMC and SMC) are the largest satellite galaxies of the Milky Way and close enough to allow for a detailed exploration of their structure and formation history. The Survey of the MAgellanic Stellar History (SMASH) is a community Dark Energy Camera (DECam) survey of the Magellanic Clouds using $\sim$50 nights to sample over $\sim$2400 deg$^2$ centered on the Clouds at $\sim$20% filling factor (but with contiguous coverage in the central regions) and to depths of $\sim$24th mag in $ugriz$. The primary goals of SMASH are to map out the extended stellar peripheries of the Clouds and uncover their complicated interaction and accretion history as well as to derive spatially-resolved star formation histories of the central regions and create a "movie" of their past star formation. Here we announce the second SMASH public data release (DR2), which contains all 197 fully-calibrated DECam fields including the main body fields in the central regions. The DR2 data are available through the Astro Data Lab hosted by the NSF's National Optical-Infrared Astronomy Research Laboratory. We highlight three science cases that make use of the SMASH DR2 data and will be published in the future: (1) preliminary star formation histories of the LMC; (2) the search for Magellanic star clusters using citizen scientists; and, (3) photometric metallicities of Magellanic Cloud stars using the DECam $u$-band., Comment: 18 pages, 9 figures. Accepted for publication in The Astronomical Journal

We announce the second data release (DR2) of the NOIRLab Source Catalog (NSC), using 412,116 public images from CTIO-4m+DECam, the KPNO-4m+Mosaic3 and the Bok-2.3m+90Prime. NSC DR2 contains over 3.9 billion unique objects, 68 billion individual source measurements, covers $\approx$35,000 square degrees of the sky, has depths of $\approx$23rd magnitude in most broadband filters with $\approx$1-2% photometric precision, and astrometric accuracy of $\approx$7 mas. Approximately 1.9 billion objects within $\approx$30,000 square degrees of sky have photometry in three or more bands. There are several improvements over NSC DR1. DR2 includes 156,662 (61%) more exposures extending over 2 more years than in DR1. The southern photometric zeropoints in $griz$ are more accurate by using the Skymapper DR1 and ATLAS-Ref2 catalogs, and improved extinction corrections were used for high-extinction regions. In addition, the astrometric accuracy is improved by taking advantage of Gaia DR2 proper motions when calibrating the WCS of individual images. This improves the NSC proper motions to $\sim$2.5 mas/yr (precision) and $\sim$0.2 mas/yr (accuracy). The combination of sources into unique objects is performed using a DBSCAN algorithm and mean parameters per object (such as mean magnitudes, proper motion, etc.) are calculated more robustly with outlier rejection. Finally, eight multi-band photometric variability indices are calculated for each object and variable objects are flagged (23 million objects). NSC DR2 will be useful for exploring solar system objects, stellar streams, dwarf satellite galaxies, QSOs, variable stars, high-proper motion stars, and transients. Several examples of these science use cases are presented. The NSC DR2 catalog is publicly available via the NOIRLab's Astro Data Lab science platform., Comment: 20 pages, 16 figures, submitted to AJ

The DESI Milky Way Survey (MWS) will observe $\ge$8 million stars between $16 < r < 19$ mag, supplemented by observations of brighter targets under poor observing conditions. The survey will permit an accurate determination of stellar kinematics and population gradients; characterize diffuse substructure in the thick disk and stellar halo; enable the discovery of extremely metal-poor stars and other rare stellar types; and improve constraints on the Galaxy's 3D dark matter distribution from halo star kinematics. MWS will also enable a detailed characterization of the stellar populations within 100 pc of the Sun, including a complete census of white dwarfs. The target catalog from the preliminary selection described here is public., Comment: 4 pages, 1 figure; published in Res. Notes AAS

We present a map of the total intrinsic reddening across ~34 deg$^{2}$ of the Small Magellanic Cloud (SMC) derived using optical ($ugriz$) and near-infrared (IR; $YJK_{\mathrm{s}}$) spectral energy distributions (SEDs) of background galaxies. The reddening map is created using a subsample of 29,274 galaxies with low levels of intrinsic reddening based on the LePhare $\chi^{2}$ minimisation SED-fitting routine. We find statistically significant enhanced levels of reddening associated with the main body of the SMC compared with regions in the outskirts [$\Delta E(B-V)\simeq 0.3$ mag]. A comparison with literature reddening maps of the SMC shows that, after correcting for differences in the volume of the SMC sampled, there is good agreement between our results and maps created using young stars. In contrast, we find significant discrepancies between our results and maps created using old stars or based on longer wavelength far-IR dust emission that could stem from biased samples in the former and uncertainties in the far-IR emissivity and the optical properties of the dust grains in the latter. This study represents one of the first large-scale categorisations of extragalactic sources behind the SMC and as such we provide the LePhare outputs for our full sample of ~500,000 sources., Comment: Accepted for publication in MNRAS, 13 pages, 6 figures and 4 tables. Tables 1, 3 and 4 will be available in full as Supporting Information in the online version of the article

Draco C1 is a known symbiotic binary star system composed of a carbon red giant and a hot, compact companion -- likely a white dwarf -- belonging to the Draco dwarf spheroidal galaxy. From near-infrared spectroscopic observations taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2), part of Sloan Digital Sky Survey IV, we provide updated stellar parameters for the cool, giant component, and constrain the temperature and mass of the hot, compact companion. Prior measurements of the periodicity of the system, based on only a few epochs of radial velocity data or relatively short baseline photometric observations, were sufficient only to place lower limits on the orbital period ($P > 300$ days). For the first time, we report precise orbital parameters for the binary system: With 43 radial velocity measurements from APOGEE spanning an observational baseline of more than 3 years, we definitively derive the period of the system to be $1220.0^{+3.7}_{-3.5}$ days. Based on the newly derived orbital period and separation of the system, together with estimates of the radius of the red giant star, we find that the hot companion must be accreting matter from the dense wind of its evolved companion., Comment: 8 pages, 4 figures, 1 table. Accepted for publication in ApJL

The periphery of the Small Magellanic Cloud (SMC) can unlock important information regarding galaxy formation and evolution in interacting systems. Here, we present a detailed study of the extended stellar structure of the SMC using deep colour-magnitude diagrams (CMDs), obtained as part of the Survey of the MAgellanic Stellar History (SMASH). Special care was taken in the decontamination of our data from MW foreground stars, including from foreground globular clusters NGC 362 and 47 Tuc. We derived the SMC surface brightness using a ``conservative'' approach from which we calculated the general parameters of the SMC, finding a staggered surface brightness profile. We also traced the fainter outskirts by constructing a stellar density profile. This approach, based on stellar counts of the oldest main sequence turn-off (MSTO) stars, uncovered a tidally disrupted stellar feature that reaches as far out as 12 degrees from the SMC centre. We also serendipitously found a faint feature of unknown origin located at $\sim 14$ degrees from the centre of the SMC and that we tentatively associated to a more distant structure. We compared our results to in-house simulations of a $1\times10^{9} M_\odot$ SMC, finding that its elliptical shape can be explained by its tidal disruption under the combined presence of the MW and the LMC. Finally, we found that the older stellar populations show a smooth profile while the younger component presents a jump in the density followed by a flat profile, confirming the heavily disturbed nature of the SMC., Comment: 18 pages, 15 figures, submitted to MNRAS

Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe], [Ti/Fe], [Mn/Fe], [Fe/H], and [Ni/Fe], abundances ratios which together span multiple nucleosynthetic channels. In a high quality sample of 182,538 giant stars, we detect twenty-one candidate clusters with more than fifteen members. Our candidate clusters are more chemically homogeneous than a population of non-member stars with similar [Mg/Fe] and [Fe/H], even in abundances not used for tagging. Group members are consistent with having the same age and fall along a single stellar-population track in logg vs. Teff space. Each group's members are distributed over multiple kpc, and the spread in their radial and azimuthal actions increases with age. We qualitatively reproduce this increase using N-body simulations of cluster dissolution in Galactic potentials that include transient winding spiral arms. Observing our candidate birth clusters with high-resolution spectroscopy in other wavebands to investigate their chemical homogeneity in other nucleosynthetic groups will be essential to confirming the efficacy of strong chemical tagging. Our initially spatially-compact but now widely dispersed candidate clusters will provide novel limits on chemical evolution and orbital diffusion in the Galactic disc, and constraints on star formation in loosely-bound groups., Comment: 15 pages, 9 figures, accepted by MNRAS

Many problems in contemporary astrophysics---from understanding the formation of black holes to untangling the chemical evolution of galaxies---rely on knowledge about binary stars. This, in turn, depends on discovery and characterization of binary companions for large numbers of different kinds of stars in different chemical and dynamical environments. Current stellar spectroscopic surveys observe hundreds of thousands to millions of stars with (typically) few observational epochs, which allows binary discovery but makes orbital characterization challenging. We use a custom Monte Carlo sampler (The Joker) to perform discovery and characterization of binary systems through radial-velocities, in the regime of sparse, noisy, and poorly sampled multi-epoch data. We use it to generate posterior samplings in Keplerian parameters for 232,531 sources released in APOGEE Data Release 16. Our final catalog contains 19,635 high-confidence close-binary (P < few years, a < few AU) systems that show interesting relationships between binary occurrence rate and location in the color-magnitude diagram. We find notable faint companions at high masses (black-hole candidates), at low masses (substellar candidates), and at very close separations (mass-transfer candidates). We also use the posterior samplings in a (toy) hierarchical inference to measure the long-period binary-star eccentricity distribution. We release the full set of posterior samplings for the entire parent sample of 232,531 stars. This set of samplings involves no heuristic "discovery" threshold and therefore can be used for myriad statistical purposes, including hierarchical inferences about binary-star populations and sub-threshold searches., Comment: 35 pages, 12 figures, but check out Figure 5. Data tables available at http://adrian.pw/apogee-dr16.html Submitted to AAS journals

Using 3D positions and kinematics of stars relative to the Sagittarius (Sgr) orbital plane and angular momentum, we identify 166 Sgr stream members observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) that also have Gaia DR2 astrometry. This sample of 63/103 stars in the Sgr trailing/leading arm are combined with an APOGEE sample of 710 members of the Sgr dwarf spheroidal core (385 of them newly presented here) to establish differences of 0.6 dex in median metallicity and 0.1 dex in [$\alpha$/Fe] between our Sgr core and dynamically older stream samples. Mild chemical gradients are found internally along each arm, but these steepen when anchored by core stars. With a model of Sgr tidal disruption providing estimated dynamical ages (i.e., stripping times) for each stream star, we find a mean metallicity gradient of 0.12 +/- 0.03 dex/Gyr for stars stripped from Sgr over time. For the first time, an [$\alpha$/Fe] gradient is also measured within the stream, at 0.02 +/- 0.01 dex/Gyr using magnesium abundances and 0.04 +/- 0.01 dex/Gyr using silicon, which imply that the Sgr progenitor had significant radial abundance gradients. We discuss the magnitude of those inferred gradients and their implication for the nature of the Sgr progenitor within the context of the current family of Milky Way satellite galaxies, and suggest that more sophisticated Sgr models are needed to properly interpret the growing chemodynamical detail we have on the Sgr system., Comment: 27 pages, 11 figures, 2 tables. Accepted for publishing in The Astrophysical Journal

We investigate the Fe, C, N, O, Mg, Al, Si, K, Ca, Ce and Nd abundances of 2283 red giant stars in 31 globular clusters from high-resolution spectra observed in both the northern and southern hemisphere by the SDSS-IV APOGEE-2 survey. This unprecedented homogeneous dataset, largest to date, allows us to discuss the intrinsic Fe spread, the shape and statistics of Al-Mg and N-C anticorrelations as a function of cluster mass, luminosity, age and metallicity for all 31 clusters. We find that the Fe spread does not depend on these parameters within our uncertainties including cluster metallicity, contradicting earlier observations. We do not confirm the metallicity variations previously observed in M22 and NGC 1851. Some clusters show a bimodal Al distribution, while others exhibit a continuous distribution as has been previously reported in the literature. We confirm more than 2 populations in $\omega$ Cen and NGC 6752, and find new ones in M79. We discuss the scatter of Al by implementing a correction to the standard chemical evolution of Al in the Milky Way. After correction, its dependence on cluster mass is increased suggesting that the extent of Al enrichment as a function of mass was suppressed before the correction. We observe a turnover in the Mg-Al anticorrelation at very low Mg in $\omega$ Cen, similar to the pattern previously reported in M15 and M92. $\omega$ Cen may also have a weak K-Mg anticorrelation, and if confirmed, it would be only the third cluster known to show such a pattern., Comment: Published in MNRAS, 31 pages, 24 figures, 7 tables. Further typos are corrected, now identical in text, figures and tables to the published paper

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17)., Comment: DR16 release: Monday Dec 9th 2019. This is the alphabetical order SDSS-IV collaboration data release paper. 25 pages, 6 figures, accepted by ApJS on 11th May 2020. Minor changes clarify or improve text and figures relative to v1

Time series observations of a single dithered field centered on the diffuse dwarf satellite galaxy Crater II were obtained with the Dark Energy Camera (DECam) at the 4m Blanco Telescope at Cerro Tololo Inter-American Observatory, Chile, uniformly covering up to two half-light radii. Analysis of the $g$ and $i$ time series results in the identification and characterization of 130 periodic variable stars, including 98 RR Lyrae stars, 7 anomalous Cepheids, and 1 SX Phoenicis star belonging to the Crater II population, and 24 foreground variables of different types. Using the large number of ab-type RR Lyrae stars present in the galaxy, we obtained a distance modulus to Crater II of $(m-M)_0=20.333\pm 0.004$ (stat) $\pm 0.07$ (sys). The distribution of the RR Lyrae stars suggests an elliptical shape for Crater II, with an ellipticity of 0.24 and a position angle of $153^\circ$. From the RR Lyrae stars we infer a small metallicity dispersion for the old population of Crater II of only 0.17 dex. There are hints that the most metal-poor stars in that narrow distribution have a wider distribution across the galaxy, while the slightly more metal rich part of the population is more centrally concentrated. Given the features in the color-magnitude diagram of Crater II, the anomalous Cepheids in this galaxy must have formed through a binary evolution channel of an old population., Comment: Accepted for publication in MNRAS

We report spectroscopic measurements of stars in the recently discovered young stellar association Price-Whelan 1 (PW 1), which was found in the vicinity of the Leading Arm (LA) of the Magellanic Stream. We obtained Magellan+MIKE high-resolution spectra of the 28 brightest stars in PW 1 and used The Cannon to determine their stellar parameters. We find that the mean metallicity of PW 1 is [Fe/H]=-1.23 with a small scatter of 0.06 dex and the mean radial velocity is Vhelio=276.7 km/s with a dispersion of 11.0 km/s. Our results are consistent in Teff, logg, and [Fe/H] with the young and metal-poor characteristics (116 Myr and [Fe/H]=-1.1) determined for PW 1 from our discovery paper. We find a strong correlation between the spatial pattern of the PW 1 stars and the LA II gas with an offset of -10.15 deg in L_MS and +1.55 deg in B_MS. The similarity in metallicity, velocity, and spatial patterns indicates that PW 1 likely originated in LA II. We find that the spatial and kinematic separation between LA II and PW 1 can be explained by ram pressure from Milky Way gas. Using orbit integrations that account for the LMC and MW halo and outer disk gas, we constrain the halo gas density at the orbital pericenter of PW 1 to be n_halo (17 kpc) = 2.7 (3) x 10^-3 atoms/cm^3 and the disk gas density at the midplane at 20 kpc to be n_disk (20 kpc,0) = 6.0 (1.8) x 10^-2 atoms/cm^3. We, therefore, conclude that PW 1 formed from the LA II of the Magellanic Stream, making it a powerful constraint on the Milky Way-Magellanic interaction., Comment: 18 pages, 13 figures, 1 table, submitted to ApJ

We present a method to map the total intrinsic reddening of a foreground extinguishing medium via the analysis of spectral energy distributions (SEDs) of background galaxies. In this pilot study, we implement this technique in two distinct regions of the Small Magellanic Cloud (SMC) - the bar and the southern outskirts - using a combination of optical and near-infrared $ugrizYJK_{\mathrm{s}}$ broadband imaging. We adopt the LePhare $\chi^{2}$-minimisation SED-fitting routine and various samples of galaxies and/or quasi-stellar objects to investigate the intrinsic reddening. We find that only when we construct reddening maps using objects classified as galaxies with low levels of intrinsic reddening (i.e. ellipticals/lenticulars and early-type spirals), the resultant maps are consistent with previous literature determinations i.e. the intrinsic reddening of the SMC bar is higher than that in the outer environs. We employ two sets of galaxy templates - one theoretical and one empirical - to test for template dependencies in the resulting reddening maps and find that the theoretical templates imply systematically higher reddening values by up to 0.20 mag in $E(B-V)$. A comparison with previous reddening maps, based on the stellar components of the SMC, typically shows reasonable agreement. There is, however, significant variation amongst the literature reddening maps as to the level of intrinsic reddening associated with the bar. Thus, it is difficult to unambiguously state that instances of significant discrepancies are the result of appreciable levels of dust not accounted for in some literature reddening maps or whether they reflect issues with our adopted methodology., Comment: Accepted for publication in MNRAS, 21 pages, 13 figures and 4 tables

Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from APOGEE-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of $\sim$0.05--1.5 \msun, we find 399 binaries, including both RV variables and SB2s. The mass ratio distribution of SB2s is consistent with a uniform for $q<0.95$ with an excess of twins with $q>0.95$. The period distribution is consistent with what has been observed in close binaries ($<10$ AU) in the evolved populations. Three systems are found to have $q\sim$0.1, with a companion located within the brown dwarf desert. There are not any strong trends in the multiplicity fraction (MF) as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at $\Sigma_*\sim30$ stars/pc$^{-2}$, and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of $\sim$2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal mass systems that is less effective in binaries with lower $q$., Comment: 25 pages, 20 figures. Accepted to AJ

The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey provides precise chemical abundances of 18 chemical elements for $\sim$ 176,000 red giant stars distributed over much of the Milky Way Galaxy (MW), and includes observations of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). The APOGEE chemical abundance patterns of Sgr have revealed that it is chemically distinct from the MW in most chemical elements. We employ a \emph{k}-means clustering algorithm to 6-dimensional chemical space defined by [(C+N)/Fe], [O/Fe], [Mg/Fe], [Al/Fe], [Mn/Fe], and [Ni/Fe] to identify 62 MW stars in the APOGEE sample that have Sgr-like chemical abundances. Of the 62 stars, 35 have \emph{Gaia} kinematics and positions consistent with those predicted by \emph{N}-body simulations of the Sgr stream, and are likely stars that have been stripped from Sgr during the last two pericenter passages ($<$ 2 Gyr ago). Another 20 of the 62 stars exhibit chemical abundances indistinguishable from the Sgr stream stars, but are on highly eccentric orbits with median $r_{\rm apo} \sim $ 25 kpc. These stars are likely the `accreted' halo population thought to be the result of a separate merger with the MW 8-11 Gyr ago. We also find one hypervelocity star candidate. We conclude that Sgr was enriched to [Fe/H] $\sim$ -0.2 before its most recent pericenter passage. If the `accreted halo' population is from one major accretion event, then this progenitor galaxy was enriched to at least [Fe/H] $\sim$ -0.6, and had a similar star formation history to Sgr before merging., Comment: 15 pages, 10 figures, accepted for publication in ApJ

We report the first APOGEE metallicities and alpha-element abundances measured for 3600 red giant stars spanning a large radial range of both the Large (LMC) and Small Magellanic Clouds (SMC), the largest Milky Way dwarf galaxies. Our sample is an order of magnitude larger than that of previous studies, and extends to much larger radial distances. These are the first results presented that make use of the newly installed Southern APOGEE instrument on the du Pont telescope at Las Campanas Observatory. Our unbiased sample of the LMC spans a large range in metallicity, from [Fe/H]=-0.2 to very metal-poor stars with [Fe/H]=-2.5, the most metal-poor Magellanic Clouds (MCs) stars detected to date. The LMC [alpha/Fe]-[Fe/H] distribution is very flat over a large metallicity range, but rises by ~0.1 dex at -1.0<[Fe/H]<-0.5. We interpret this as a sign of the known recent increase in MC star-formation activity, and are able to reproduce the pattern with a chemical evolution model that includes a recent "starburst". At the metal-poor end, we capture the increase of [alpha/Fe] with decreasing [Fe/H], and constrain the "alpha-knee" to [Fe/H]<-2.2 in both MCs, implying a low star-formation efficiency of ~0.01 Gyr^-1. The MC knees are more metal poor than those of less massive Milky Way (MW) dwarf galaxies such as Fornax, Sculptor, or Sagittarius. One possible interpretation is that the MCs formed in a lower-density environment than the MW, a hypothesis that is consistent with the paradigm that the MCs fell into the MW's gravitational potential only recently., Comment: 22 pages, 20 figues, 1 table, accepted for publication in ApJ

The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical structure, including showing that an identification of chemical and geometric thick discs is not exact, but the origin of the chemical dichotomy has remained unclear. Here we demonstrate that a dichotomy arises naturally if the early gas-rich disc fragments, leading to some fraction of the star formation occurring in clumps of the type observed in high-redshift galaxies. These clumps have high star formation rate density. They, therefore, enrich rapidly, moving from the low-alpha to the high-alpha sequence, while more distributed star formation produces the low-alpha sequence. We demonstrate that this model produces a chemically-defined thick disc that has many of the properties of the Milky Way's thick disc. Because clump formation is common in high redshift galaxies, we predict that chemical bimodalities are common in massive galaxies., Comment: Accepted to MNRAS

We report the discovery of a young (${\rm age} \sim 130~{\rm Myr}$), low-mass ($M \sim 1200~{\rm M}_\odot$), metal-poor ($[{\rm Fe}/{\rm H}] \sim -1.1$) stellar association at a heliocentric distance $D \approx 29~{\rm kpc}$, placing it far into the Milky Way halo. At its present Galactocentric position $(R, z) \sim (23, 15)~{\rm kpc}$, the association is (on the sky) near the leading arm of the gas stream emanating from the Magellanic cloud system, but is located $\approx 60^\circ$ from the Large Magellanic Cloud (LMC) center on the other side of the Milky Way disk. If we assume that the cluster is co-located with HI gas in the stream, we directly measure the distance to the leading arm of the Magellanic stream. The measured distance is inconsistent with stream predictions from models of the LMC/SMC interaction and infall into the Milky Way that do not account for ram pressure and gas interaction with Milky Way disk. The estimated age of the cluster is consistent with the time of last passage of the leading arm gas through the Galactic midplane, and we therefore speculate that this star-formation event was triggered by its the last disk midplane passage. Most details of this idea remain a puzzle: the Magellanic stream has low column density, the Milky Way disk at this large radius has low gas density, and the relative velocity of the leading arm gas and Milky Way gas is large. However it formed, the discovery of a young stellar cluster in the Milky Way halo presents an interesting opportunity for study. This cluster was discovered with Gaia astrometry and photometry alone, but folow-up DECam photometry was crucial for measuring its properties., Comment: 18 pages, 9 figures, 6 tables; published in ApJ

We map the trends of elemental abundance ratios across the Galactic disk, spanning R = 3-15 kpc and midplane distance |Z|= 0-2 kpc, for 15 elements in a sample of 20,485 stars measured by the SDSS/APOGEE survey (O, Na, Mg, Al, Si, P, S, K, Ca, V, Cr, Mn, Fe, Co, Ni). Adopting Mg rather than Fe as our reference element, and separating stars into two populations based on [Fe/Mg], we find that the median trends of [X/Mg] vs. [Mg/H] in each population are nearly independent of location in the Galaxy. The full multi-element cartography can be summarized by combining these nearly universal median sequences with our measured metallicity distribution functions and the relative proportions of the low-[Fe/Mg] (high-alpha) and high-[Fe/Mg] (low-alpha) populations, which depend strongly on R and |Z|. We interpret the median sequences with a semi-empirical "2-process" model that describes both the ratio of core collapse and Type Ia supernova contributions to each element and the metallicity dependence of the supernova yields. These observationally inferred trends can provide strong tests of supernova nucleosynthesis calculations. Our results lead to a relatively simple picture of abundance ratio variations in the Milky Way, in which the trends at any location can be described as the sum of two components with relative contributions that change systematically and smoothly across the Galaxy. Deviations from this picture and future extensions to other elements can provide further insights into the physics of stellar nucleosynthesis and unusual events in the Galaxy's history., Comment: 29 pp, 22 figs, submitted to AAS journals, comments welcome

The nature of the Triangulum-Andromeda (TriAnd) system has been debated since the discovery of this distant, low-latitude Milky Way (MW) overdensity more than a decade ago. Explanations for its origin are either as a halo substructure from the disruption of a dwarf galaxy or a distant extension of the Galactic disk. We test these hypotheses using chemical abundances of a dozen TriAnd members from the Sloan Digital Sky Survey's 14th Data Release of Apache Point Observatory Galactic Evolution Experiment (APOGEE) data to compare to APOGEE abundances of stars with similar metallicity from both the Sagittarius (Sgr) dSph, and the outer MW disk. We find that TriAnd stars are chemically distinct from Sgr across a variety of elements, (C+N), Mg, K, Ca, Mn, and Ni, with a separation in [X/Fe] of about 0.1 to 0.4 dex depending on the element. Instead, the TriAnd stars, with a median metallicity of about -0.8, exhibit chemical abundance ratios similar to those of the lowest metallicity ([Fe/H] ~ -0.7) stars in the outer Galactic disk, and are consistent with expectations of extrapolated chemical gradients in the outer disk of the MW. These results suggest that TriAnd is associated with the MW disk, and, therefore, that the disk extends to this overdensity --- i.e., past a Galactocentric radius of 24 kpc --- albeit vertically perturbed about 7 kpc below the nominal disk midplane in this region of the Galaxy., Comment: 8 pages, 4 figures, 1 table. Accepted for publishing in The Astrophysical Journal Letters

We present the detection of very extended stellar populations around the Large Magellanic Cloud (LMC) out to R~21 degrees, or ~18.5 kpc at the LMC distance of 50 kpc, as detected in the Survey of the MAgellanic Stellar History (SMASH) performed with the Dark Energy Camera on the NOAO Blanco 4m Telescope. The deep (g~24) SMASH color magnitude diagrams (CMDs) clearly reveal old (~9 Gyr), metal-poor ([Fe/H]=-0.8 dex) main-sequence stars at a distance of 50 kpc. The surface brightness of these detections is extremely low with our most distant detection having 34 mag per arcsec squared in g-band. The SMASH radial density profile breaks from the inner LMC exponential decline at ~13-15 degrees and a second component at larger radii has a shallower slope with power-law index of -2.2 that contributes ~0.4% of the LMC's total stellar mass. In addition, the SMASH densities exhibit large scatter around our best-fit model of ~70% indicating that the envelope of stellar material in the LMC periphery is highly disturbed. We also use data from the NOAO Source catalog to map the LMC main-sequence populations at intermediate radii and detect a steep dropoff in density on the eastern side of the LMC (at R~8 deg) as well as an extended structure to the far northeast. These combined results confirm the existence of a very extended, low-density envelope of stellar material with disturbed shape around the LMC. The exact origin of this structure remains unclear but the leading options include a classical accreted halo or tidally stripped outer disk material., Comment: 9 pages, 5 figues, 1 table, submitted to ApJ