Search

We present the 2020 version of the Siena Galaxy Atlas (SGA-2020), a multiwavelength optical and infrared imaging atlas of 383,620 nearby galaxies. The SGA-2020 uses optical grz imaging over ≈20,000 deg2 from the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys Data Release 9 and infrared imaging in four bands (spanning 3.4-22 μm) from the 6 year unWISE coadds; it is more than 95% complete for galaxies larger than R(26) ≈ 25″ and r < 18 measured at the 26 mag arcsec−2 isophote in the r band. The atlas delivers precise coordinates, multiwavelength mosaics, azimuthally averaged optical surface-brightness profiles, model images and photometry, and additional ancillary metadata for the full sample. Coupled with existing and forthcoming optical spectroscopy from the DESI, the SGA-2020 will facilitate new detailed studies of the star formation and mass assembly histories of nearby galaxies; enable precise measurements of the local velocity field via the Tully-Fisher and fundamental plane relations; serve as a reference sample of lasting legacy value for time-domain and multimessenger astronomical events; and more.

Deep optical and near-infrared imaging of the entire Galactic plane is essential for understanding our Galaxy’s stars, gas, and dust. The second data release of the Dark Energy Camera (DECam) Plane Survey extends the five-band optical and near-infrared survey of the southern Galactic plane to cover 6.5% of the sky, ∣b∣ ≤ 10°, and 6° > ℓ > −124°, complementary to coverage by Pan-STARRS1. Typical single-exposure effective depths, including crowding effects and other complications, are 23.5, 22.6, 22.1, 21.6, and 20.8 mag in g, r, i, z, and Y bands, respectively, with around 1″ seeing. The survey comprises 3.32 billion objects built from 34 billion detections in 21,400 exposures, totaling 260 hr open shutter time on the DECam at Cerro Tololo. The data reduction pipeline features several improvements, including the addition of synthetic source injection tests to validate photometric solutions across the entire survey footprint. A convenient functional form for the detection bias in the faint limit was derived and leveraged to characterize the photometric pipeline performance. A new postprocessing technique was applied to every detection to debias and improve uncertainty estimates of the flux in the presence of structured backgrounds, specifically targeting nebulosity. The images and source catalogs are publicly available at http://decaps.skymaps.info/.

Modern surveys of gravitational microlensing events have progressed to detecting thousands per year, and surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is the microlensing Einstein radius crossing time (t E) distribution. However, systematics in individual light curves as well as oversimplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian process model. We used light curves for nearly 10,000 OGLE-III and -IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the t E distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and the previous analyses overestimated the number of t E > 100 day events due to their oversimplistic model ignoring parallax effects. We use our fits to identify the hundreds filling a regime in the microlensing parameter space that are 50% pure of black holes. Finally, we have released the largest-ever catalog of Markov Chain Monte Carlo parameter estimates for microlensing events.

The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemodynamical mapping of the Milky Way. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation modes: (i) "Ancillary Science Programs"competitively awarded to Sloan Digital Sky Survey IV PIs through proposal calls in 2015 and 2017 for the pursuit of new scientific avenues outside the main survey, and (ii) an effective 1.5 yr expansion of the survey, known as the Bright Time Extension (BTX), made possible through accrued efficiency gains over the first years of the APOGEE-2N project. For the 23 distinct ancillary programs, we provide descriptions of the scientific aims, target selection, and how to identify these targets within the APOGEE-2 sample. The BTX permitted changes to the main survey strategy, the inclusion of new programs in response to scientific discoveries or to exploit major new data sets not available at the outset of the survey design, and expansions of existing programs to enhance their scientific success and reach. After describing the motivations, implementation, and assessment of these programs, we also leave a summary of lessons learned from nearly a decade of APOGEE-1 and APOGEE-2 survey operations. A companion paper, F. Santana et al. (submitted; AAS29036), provides a complementary presentation of targeting modifications relevant to APOGEE-2 operations in the Southern Hemisphere.

The CatWISE2020 Catalog consists of 1,890,715,640 sources over the entire sky selected from Wide-field Infrared Survey Explorer (WISE) and NEOWISE survey data at 3.4 and 4.6 μm (W1 and W2) collected from 2010 January 7 to 2018 December 13. This data set adds two years to that used for the CatWISE Preliminary Catalog, bringing the total to six times as many exposures spanning over 16 times as large a time baseline as the AllWISE catalog. The other major change from the CatWISE Preliminary Catalog is that the detection list for the CatWISE2020 Catalog was generated using crowdsource from Schlafly et al., while the CatWISE Preliminary Catalog used the detection software used for AllWISE. These two factors result in roughly twice as many sources in the CatWISE2020 Catalog. The scatter with respect to Spitzer photometry at faint magnitudes in the COSMOS field, which is out of the Galactic Plane and at low ecliptic latitude (corresponding to lower WISE coverage depth) is similar to that for the CatWISE Preliminary Catalog. The 90% completeness depth for the CatWISE2020 Catalog is at W1 = 17.7 mag and W2 = 17.5 mag, 1.7 mag deeper than in the CatWISE Preliminary Catalog. In comparison to Gaia, CatWISE2020 motions are accurate at the 20 mas yr-1 level for W1∼15 mag sources and at the ∼100 mas yr-1 level for W1∼17 mag sources. This level of accuracy represents a 12 improvement over AllWISE. The CatWISE catalogs are available in the WISE/NEOWISE Enhanced and Contributed Products area of the NASA/IPAC Infrared Science Archive.

The Dark Energy Spectroscopic Instrument Legacy Surveys, a combination of three ground-based imaging surveys, have mapped 16,000 deg2 in three optical bands (g, r, and z) to a depth 1-2 mag deeper than the Sloan Digital Sky Survey. Our work addresses one of the major challenges of wide-field imaging surveys conducted at ground-based observatories: the varying depth that results from varying observing conditions at Earth-bound sites. To mitigate these effects, the Legacy Surveys (the Dark Energy Camera Legacy Survey, or DECaLS; the Mayall z-band Legacy Survey, or MzLS; and the Beiijing-Arizona Sky Survey, or BASS) employed a unique strategy to dynamically adjust the exposure times as rapidly as possible in response to the changing observing conditions. We present the tiling and observing strategies used by the first two of these surveys. We demonstrate that the tiling and dynamic observing strategies jointly result in a more uniform-depth survey that has higher efficiency for a given total observing time compared with the traditional approach of using fixed exposure times.

We present catalogues of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3, and 3D-HST surveys. These catalogues incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalogue and Moffat profile point spread function fits. We test this catalogue with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SEXTRACTOR catalogues from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogues are publicly available at http://d-scholarship.pitt.edu/36064/. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data.

The DESI Legacy Imaging Surveys (http://legacysurvey.org/) are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ≈14,000 deg2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12, and 22 μm) observed by the Wide-field Infrared Survey Explorer satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We present a new technique to determine distances to major star-forming regions across the Perseus Molecular Cloud, using a combination of stellar photometry, astrometric data, and 12CO spectral-line maps. Incorporating the Gaia DR2 parallax measurements when available, we start by inferring the distance and reddening to stars from their Pan-STARRS1 and Two Micron All Sky Survey photometry, based on a technique presented by Green et al. and implemented in their 3D "Bayestar" dust map of three-quarters of the sky. We then refine their technique by using the velocity slices of a CO spectral cube as dust templates and modeling the cumulative distribution of dust along the line of sight toward these stars as a linear combination of the emission in the slices. Using a nested sampling algorithm, we fit these per-star distance-reddening measurements to find the distances to the CO velocity slices toward each star-forming region. This results in distance estimates explicitly tied to the velocity structure of the molecular gas. We determine distances to the B5, IC 348, B1, NGC 1333, L1448, and L1451 star-forming regions and find that individual clouds are located between ≈275 and 300 pc, with typical combined uncertainties of ≈5%. We find that the velocity gradient across Perseus corresponds to a distance gradient of about 25 pc, with the eastern portion of the cloud farther away than the western portion. We determine an average distance to the complex of 294 ± 17 pc, about 60 pc further than the distance derived to the western portion of the cloud using parallax measurements of water masers associated with young stellar objects. The method we present is not limited to the Perseus Complex, but may be applied anywhere on the sky with adequate CO data in the pursuit of more accurate 3D maps of molecular clouds in the solar neighborhood and beyond.

A precise extinction law is a critical input when interpreting observations of highly reddened sources such as young star clusters and the Galactic Center (GC). We use Hubble Space Telescope observations of a region of moderate extinction and a region of high extinction to measure the optical and near-infrared extinction law (0.8-2.2 μm). The moderate-extinction region is the young massive cluster Westerlund 1 (Wd1; A Ks ∼ 0.6 mag), where 453 proper-motion selected main-sequence stars are used to measure the shape of the extinction law. To quantify the shape, we define the parameter , which behaves similarly to a color-excess ratio, but is continuous as a function of wavelength. The high-extinction region is the GC (A Ks ∼ 2.5 mag), where 819 red clump stars are used to determine the normalization of the law. The best-fit extinction law is able to reproduce the Wd1 main-sequence colors, which previous laws misestimate by 10%-30%. The law is inconsistent with a single power law, even when only the near-infrared filters are considered, and has A F125W/A Ks and A F814W/A Ks values that are 18% and 24% higher than the commonly used Nishiyama et al. law, respectively. Using this law, we recalculate the Wd1 distance to be 3905 ± 422 pc from published observations of the eclipsing binary W13. This new extinction law should be used for highly reddened populations in the Milky Way, such as the Quintuplet cluster and Young Nuclear Cluster. A python code is provided to generate the law for future use.

RR Lyrae stars may be the best practical tracers of Galactic halo (sub-)structure and kinematics. The PanSTARRS1 (PS1) 3π survey offers multi-band, multi-epoch, precise photometry across much of the sky, but a robust identification of RR Lyrae stars in this data set poses a challenge, given PS1's sparse, asynchronous multi-band light curves (≲12 epochs in each of five bands, taken over a 4.5 year period). We present a novel template fitting technique that uses well-defined and physically motivated multi-band light curves of RR Lyrae stars, and demonstrate that we get accurate period estimates, precise to 2 s in >80% of cases. We augment these light-curve fits with other features from photometric time-series and provide them to progressively more detailed machinelearned classification models. From these models, we are able to select the widest (three-fourths of the sky) and deepest (reaching 120 kpc) sample of RR Lyrae stars to date. The PS1 sample of ∼45,000 RRab stars is pure (90%) and complete (80% at 80 kpc) at high galactic latitudes. It also provides distances that are precise to 3%, measured with newly derived period luminosity relations for optical/near-infrared PS1 bands. With the addition of proper motions from Gaia and radial velocity measurements from multi-object spectroscopic surveys, we expect the PS1 sample of RR Lyrae stars to become the premier source for studying the structure, kinematics, and the gravitational potential of the Galactic halo. The techniques presented in this study should translate well to other sparse, multi-band data sets, such as those produced by the Dark Energy Survey and the upcoming Large Synoptic Survey Telescope Galactic plane sub-survey.

A distant, as yet unseen ninth planet has been invoked to explain various observations of the outer solar system. While such a "Planet Nine," if it exists, is most likely to be discovered via reflected light in the optical, it may emit much more strongly at 3-5 μm than simple blackbody predictions would suggest, depending on its atmospheric properties. As a result, Planet Nine may be detectable at 3.4 μm with the Wide-field Infrared Survey Explorer, but single exposures are too shallow except at relatively small distances ( au). We develop a method to search for Planet Nine far beyond the W1 single-exposure sensitivity, to distances as large as 800 au, using inertial coadds of W1 exposures binned into ∼1 day intervals. We apply our methodology to a ∼2000 square degree testbed sky region which overlaps a southern segment of Planet Nine's anticipated orbital path. We do not detect a plausible Planet Nine candidate, but are able to derive a detailed completeness curve, ruling out its presence within the parameter space searched at W1 < 16.66 (90% completeness). Our method uses all publicly available W1 imaging, spanning 2010 January to 2015 December, and will become more sensitive with future NEOWISE-Reactivation releases of additional W1 exposures. We anticipate that our method will be applicable to the entire high Galactic latitude sky, and we will extend our search to that full footprint in the near future.

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ∼220,000 optically variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ∼320 deg2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population based on their Hα emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ∼15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

Using the Pan-STARRS1 survey, we derive limiting magnitude, spatial completeness, and density maps that we use to probe the three-dimensional structure and estimate the stellar mass of the so-called Monoceros Ring. The Monoceros Ring is an enormous and complex stellar sub-structure in the outer Milky Way disk. It is most visible across the large Galactic Anticenter region, 120° < l < 240°, -30° < b < +40°. We estimate its stellar mass density profile along every line of sight in 2° × 2° pixels over the entire 30,000 deg2 Pan-STARRS1 survey using the previously developed MATCH software. By parsing this distribution into a radially smooth component and the Monoceros Ring, we obtain its mass and distance from the Sun along each relevant line of sight. The Monoceros Ring is significantly closer to us in the south (6 kpc) than in the north (9 kpc). We also create 2D cross-sections parallel to the Galactic plane that show 135° of the Monoceros Ring in the south and 170° of the Monoceros Ring in the north. We show that the northern and southern structures are also roughly concentric circles, suggesting that they may be waves rippling from a common origin. Excluding the Galactic plane ∼±4°, we observe an excess mass of 4 × 106M⊙ across 120° < l < 240°. If we interpolate across the Galactic plane, we estimate that this region contains 8×106M⊙. If we assume (somewhat boldly) that the Monoceros Ring is a set of two Galactocentric rings, its total mass is 6×107M⊙. Finally, if we assume that it is a set of two circles centered at a point 4 kpc from the Galactic center in the anti-central direction, as our data suggests, we estimate its mass to be 4 × 107M⊙.

We present deep optical images of the Large and Small Magellanic Clouds (LMC and SMC) using a low cost telephoto lens with a wide field of view to explore stellar substructure in the outskirts of the stellar disk of the LMC (

The spatial structure of stellar populations with different chemical abundances in the Milky Way (MW) contains a wealth of information on Galactic evolution over cosmic time. We use data on 14,699 red-clump stars from the APOGEE survey, covering 4 kpc∼ R≳15 kpc, to determine the structure of mono-abundance populations (MAPs)-stars in narrow bins in [α/Fe] and [Fe H]-accounting for the complex effects of the APOGEE selection function and the spatially variable dust obscuration. We determine that all MAPs with enhanced [α/Fe] are centrally concentrated and are well-described as exponentials with a scale length of 2.2 ± 0.2 kpc over the whole radial range of the disk. We discover that the surface-density profiles of low-[α/Fe] MAPs are complex: they do not monotonically decrease outwards, but rather display a peak radius ranging from ≈5 to ≈13 kpc at low [Fe H]. The extensive radial coverage of the data allows us to measure radial trends in the thickness of each MAP. While high-[α/Fe] MAPs have constant scale heights, low-[α/Fe] MAPs flare. We confirm, now with highprecision abundances, previous results that each MAP contains only a single vertical scale height and that low- [Fe H], low-[α/Fe] and high-[Fe H], high-[α/Fe] MAPs have intermediate (hZ≈300600 pc) scale heights that smoothly bridge the traditional thin- and thick-disk divide. That the high-[α/Fe], thick disk components do not flare is strong evidence against their thickness being caused by radial migration. The correspondence between the radial structure and chemical-enrichment age of stellar populations is clear confirmation of the inside-out growth of galactic disks. The details of these relations will constrain the variety of physical conditions under which stars form throughout the MW disk.

We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1. Using point-spread function (PSF) photometry of 60 million stars with 16 < r < 20, we derive a model of amplifier gain and flat-field corrections with per-run rms residuals of 3 millimagnitudes (mmag) in griz bands and 15 mmag in u band. The new photometric zero points are adjusted to leave the median in the Galactic north unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. The recalibrated stellar PSF photometry of SDSS and PS1 has an rms difference of {9, 7, 7, 8} mmag in griz, respectively, when averaged over 15′ regions.

Inferences about the spatial density or phase-space structure of stellar populations in the Milky Way require a precise determination of the effective survey volume. The volume observed by surveys such as Gaia or near-infrared spectroscopic surveys, which have good coverage of the Galactic midplane region, is highly complex because of the abundant small-scale structure in the three-dimensional interstellar dust extinction. We introduce a novel framework for analyzing the importance of small-scale structure in the extinction. This formalism demonstrates that the spatially complex effect of extinction on the selection function of a pencil-beam or contiguous sky survey is equivalent to a low-pass filtering of the extinction-affected selection function with the smooth density field. We find that the angular resolution of current 3D extinction maps is sufficient for analyzing Gaia sub-samples of millions of stars. However, the current distance resolution is inadequate and needs to be improved by an order of magnitude, especially in the inner Galaxy. We also present a practical and efficient method for properly taking the effect of extinction into account in analyses of Galactic structure through an effective selection function. We illustrate its use with the selection function of red-clump stars in APOGEE using and comparing a variety of current 3D extinction maps.

In area and depth, the Pan-STARRS1 (PS1) 3π survey is unique among many-epoch, multi-band surveys and has enormous potential for the all-sky identification of variable sources. PS1 has observed the sky typically seven times in each of its five bands (grizy) over 3.5 years, but unlike SDSS, not simultaneously across the bands. Here we develop a new approach for quantifying statistical properties of non-simultaneous, sparse, multi-color light curves through light curve structure functions, effectively turning PS1 into a ∼35-epoch survey. We use this approach to estimate variability amplitudes and timescales (ωr, τ) for all point sources brighter than rP1 = 21.5 mag in the survey. With PS1 data on SDSS Stripe 82 as "ground truth," we use a Random Forest Classifier to identify QSOs and RR Lyrae based on their variability and their mean PS1 and WISE colors. We find that, aside from the Galactic plane, QSO and RR Lyrae samples of purity ∼75% and completeness ∼92% can be selected. On this basis we have identified a sample of ∼1,000,000 QSO candidates, as well as an unprecedentedly large and deep sample of ∼150,000 RR Lyrae candidates with distances from ∼10 to ∼120 kpc. Within the Draco dwarf spheroidal, we demonstrate a distance precision of 6% for RR Lyrae candidates. We provide a catalog of all likely variable point sources and likely QSOs in PS1, a total of 25.8 × 106 sources.

We present the discovery of three new Milky Way satellites from our search for compact stellar overdensities in the photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS 1, or PS1) 3π survey. The first satellite, Laevens 3, is located at a heliocentric distance of d = 67 ± 3 kpc. With a total magnitude of MV = -4.4 ± 0.3 and a half-light radius of rh = 7 ± 2 pc, its properties resemble those of outer halo globular clusters. The second system, Draco II/Laevens 4, is a closer and fainter satellite (d ∼ 20 kpc, MV = -2.9 ± 0.8), whose uncertain size (rh = 19-6+8) renders its classification difficult without kinematic information; it could either be a faint and extended globular cluster or a faint and compact dwarf galaxy. The third satellite, Sagittarius II/Laevens 5 (Sgr II), has an ambiguous nature, as it is either the most compact dwarf galaxy or the most extended globular cluster in its luminosity range (rh = 37-8+9 and MV = -5.2 ± 0.4). At a heliocentric distance of 67 ± 5 kpc, this satellite lies intriguingly close to the expected location of the trailing arm of the Sagittarius stellar stream behind the Sagittarius dwarf spheroidal galaxy (Sgr dSph). If confirmed through spectroscopic follow up, this connection would locate this part of the trailing arm of the Sagittarius stellar stream that has so far gone undetected. It would further suggest that Sgr II was brought into the Milky Way halo as a satellite of the Sgr dSph.