We explore in this paper the use of neural networks designed for point-clouds and sets on a new meta-learning task. We present experiments on the astronomical challenge of characterizing the stellar population of stellar streams. Stellar streams are elongated structures of stars in the outskirts of the Milky Way that form when a (small) galaxy breaks up under the Milky Way's gravitational force. We consider that we obtain, for each stream, a small 'support set' of stars that belongs to this stream. We aim to predict if the other stars in that region of the sky are from that stream or not, similar to one-class classification. Each "stream task" could also be transformed into a binary classification problem in a highly imbalanced regime (or supervised anomaly detection) by using the much bigger set of "other" stars and considering them as noisy negative examples. We propose to study the problem in the meta-learning regime: we expect that we can learn general information on characterizing a stream's stellar population by meta-learning across several streams in a fully supervised regime, and transfer it to new streams using only positive supervision. We present a novel use of Deep Sets, a model developed for point-cloud and sets, trained in a meta-learning fully supervised regime, and evaluated in a one-class classification setting. We compare it against Random Forests (with and without self-labeling) in the classic setting of binary classification, retrained for each task. We show that our method outperforms the Random-Forests even though the Deep Sets is not retrained on the new tasks, and accesses only a small part of the data compared to the Random Forest. We also show that the model performs well on a real-life stream when including additional fine-tuning.

We use photometry from the DECam Legacy Survey to detect candidate tidal tails extending ~5 deg on either side of the Palomar 13 globular cluster. The tails are aligned with the proper motion of Palomar 13 and are consistent with its old, metal-poor stellar population. We identify three RR Lyrae stars that are plausibly associated with the tails, in addition to four previously known in the cluster. From these RR Lyrae stars, we find that the mean distance to the cluster and tails is $23.6 \pm 0.2$ kpc and estimate the total (initial) luminosity of the cluster to be $L_V=5.1^{+9.7}_{-3.4}\times 10^3 L_\odot$, consistent with previous claims that its initial luminosity was higher than its current luminosity. Combined with previously-determined proper motion and radial velocity measurements of the cluster, we find that Palomar 13 is on a highly eccentric orbit ($e\sim 0.8$) with a pericenter of ~9 kpc and an apocenter of ~69 kpc, and a recent pericentric passage of the cluster ~75 Myr ago. We note a prominent linear structure in the interstellar dust map that runs parallel to the candidate tidal features, but conclude that reddening due to dust is unlikely to account for the structure that we observe. If confirmed, the Palomar 13 stellar stream would be one of very few streams with a known progenitor system, making it uniquely powerful for studying the disruption of globular clusters, the formation of the stellar halo, and the distribution of matter within our Galaxy.

A product of two Gaussians (or normal distributions) is another Gaussian. That's a valuable and useful fact! Here we use it to derive a refactoring of a common product of multivariate Gaussians: The product of a Gaussian likelihood times a Gaussian prior, where some or all of those parameters enter the likelihood only in the mean and only linearly. That is, a linear, Gaussian, Bayesian model. This product of a likelihood times a prior pdf can be refactored into a product of a marginalized likelihood (or a Bayesian evidence) times a posterior pdf, where (in this case) both of these are also Gaussian. The means and variance tensors of the refactored Gaussians are straightforward to obtain as closed-form expressions; here we deliver these expressions, with discussion. The closed-form expressions can be used to speed up and improve the precision of inferences that contain linear parameters with Gaussian priors. We connect these methods to inferences that arise frequently in physics and astronomy. If all you want is the answer, the question is posed and answered at the beginning of Section 3. We show two toy examples, in the form of worked exercises, in Section 4. The solutions, discussion, and exercises in this Note are aimed at someone who is already familiar with the basic ideas of Bayesian inference and probability.

The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries.

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).

We present new kinematic data for the Ophiuchus stellar stream. Spectra have been taken of member candidates at the MMT telescope using Hectospec, Hectochelle and Binospec, which provide more than 1800 new velocities. Combined with proper motion measurements of stars in the field by the Gaia - DR2 catalog, we have derived stream membership probabilities, resulting in the detection of more than 200 likely members. These data show the stream extends to more than three times the length shown in the discovery data. A spur to the main stream is also detected. The high resolution spectra allow us to resolve the stellar velocity dispersion, found to be $1.6 \pm 0.3 $ km/sec.

The recent detection of circularly polarized, long-duration (>8 hr) low-frequency (~150 MHz) radio emission from the M4.5 dwarf GJ 1151 has been interpreted as arising from a star-planet interaction via the electron cyclotron maser instability. The existence or parameters of the proposed planets have not been determined. Using 20 new HARPS-N observations, we put 99th-percentile upper limits on the mass of any close companion to GJ 1151 at Msini < 5.6 M earth. With no stellar, brown dwarf, or giant planet companion likely in a close orbit, our data are consistent with detected radio emission emerging from a magnetic interaction between a short-period terrestrial-mass planet and GJ 1151.

Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: Its $\sim 20^\circ$ tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and substructure fraction of the dark matter distribution of the Galaxy. However, most details of the phase-space distribution of the stream are not fully explained, and dynamical models that use the stream for other inferences are therefore incomplete. Here we aim to measure distance and kinematic properties along the Pal 5 stream in order to motivate improved models of the system. We use a large catalog of RR Lyrae-type stars (RRLs) with astrometric data from the Gaia mission to probabilistically identify RRLs in the Pal 5 stream. RRLs are useful because they are intrinsically-luminous standard candles and their distances can be inferred with small relative precision ($\sim3\%$). By building a probabilistic model of the Pal 5 cluster and stream in proper motion and distance, we find 27 RRLs consistent with being members of the cluster (10) and stream (17). Using these RRLs, we detect gradients in distance and proper motion along the stream, and provide an updated measurement of the distance to the Pal 5 cluster using the RRLs, $d = 20.6 \pm 0.2~\textrm{kpc}$. We provide a catalog of Pal 5 RRLs with inferred membership probabilities for future modeling work.

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.

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.

The $100^\circ$-long thin stellar stream in the Milky Way halo, GD-1, has an ensemble of features that may be due to dynamical interactions. Using high-resolution MMT/Hectochelle spectroscopy we show that a spur of GD-1-like stars outside of the main stream are kinematically and chemically consistent with the main stream. In the spur, as in the main stream, GD-1 has a low intrinsic radial velocity dispersion, $\sigma_{V_r}\lesssim1\,\rm km\,s^{-1}$, is metal-poor, $\rm [Fe/H]\approx-2.3$, with little $\rm [Fe/H]$ spread and some variation in $\rm [\alpha/Fe]$ abundances, which point to a common globular cluster progenitor. At a fixed location along the stream, the median radial velocity offset between the spur and the main stream is smaller than $0.5\,\rm km\,s^{-1}$, comparable to the measurement uncertainty. A flyby of a massive, compact object can change orbits of stars in a stellar stream and produce features like the spur observed in GD-1. In this scenario, the radial velocity of the GD-1 spur relative to the stream constrains the orbit of the perturber and its current on-sky position to $\approx5,000\,\rm deg^2$. The family of acceptable perturber orbits overlaps the stellar and dark-matter debris of the Sagittarius dwarf galaxy in present-day position and velocity. This suggests that GD-1 may have been perturbed by a globular cluster or an extremely compact dark-matter subhalo formerly associated with Sagittarius.

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.

Theoretical models of stellar evolution predict that most of the lithium inside a star is destroyed as the star becomes a red giant. However, observations reveal that about 1% of red giants are peculiarly rich in lithium, often exceeding the amount in the interstellar medium or predicted from the Big Bang. With only about 150 lithium-rich giants discovered in the past four decades, and no distinguishing properties other than lithium enhancement, the origin of lithium-rich giant stars is one of the oldest problems in stellar astrophysics. Here we report the discovery of 2,330 low-mass (1 to 3$\,M_\odot$) lithium-rich giant stars, which we argue are consistent with internal lithium production that is driven by tidal spin-up by a binary companion. Our sample reveals that most lithium-rich giants have helium-burning cores ($80^{+7}_{-6}\%$), and that the frequency of lithium-rich giants rises with increasing stellar metallicity. We find that while planet accretion may explain some lithium-rich giants, it cannot account for the majority that have helium-burning cores. We rule out most other proposed explanations as the primary mechanism for lithium-rich giants, including all stages related to single star evolution. Our analysis shows that giants remain lithium-rich for only about two million years. A prediction from this lithium depletion timescale is that most lithium-rich giants with a helium-burning core have a binary companion.

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

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.

Stars that escape globular clusters form tidal tails that are predominantly shaped by the global distribution of mass in the Galaxy, but also preserve a historical record of small-scale perturbations. Using deep $grz$ photometry from DECaLS, we present highly probable members of the tidal tails associated with the disrupting globular cluster Palomar 5. These data yield the cleanest view of a stellar stream beyond $\sim20\,\rm kpc$ and reveal: (1) a wide, low surface-brightness extension of the leading tail; (2) significant density variations along the stream; and (3) sharp changes in the direction of both the leading and the trailing tail. In the fiducial Milky Way model, a rotating bar perturbs the Palomar 5 tails and can produce streams with similar width and density profiles to those observed. However, the deviations of the stream track in this simple model do not match those observed in the Palomar 5 trailing tail, indicating the need for an additional source of perturbation. These discoveries open up the possibility of measuring the population of perturbers in the Milky Way, including dark-matter subhalos, with an ensemble of stellar streams and deep photometry alone.

Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic surveys, namely RAVE, APOGEE and LAMOST. In this dataset, we identify several structures using the HDBSCAN clustering algorithm. We discuss their properties and possible links using all the available chemical and dynamical information. Results. In concordance with previous work, we find that stars with [Fe/H] $<-1$ have more retrograde motions than those with [Fe/H] $>-1$. The retrograde halo contains a mixture of debris from objects like Gaia-Enceladus, Sequoia, and even the chemically defined thick-disc. We find that the Sequoia has a smaller range in orbital energies than previously suggested and is confined to high-energy. Sequoia could be a small galaxy in itself, but since it overlaps both in integrals-of-motion space and chemical abundance space with the less bound debris of Gaia-Enceladus, its nature cannot be fully settled yet. In the low-energy part of the halo we find evidence for at least one more distinct structure: Thamnos. Stars in Thamnos are on low inclination, mildly eccentric retrograde orbits, moving at $v_{\phi}\approx-150$ km/s, and are chemically distinct from the other structures. Conclusions. Even with the excellent Gaia DR2 data it remains challenging to piece together all the fragments found in the retrograde halo. At this point, we are very much in need of large datasets with high-quality high-resolution spectra and tailored high-resolution hydrodynamical simulations of galaxy mergers.

We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer HI disk of the Milky Way, using both test particle and Smoothed Particle Hydrodynamics simulations. We find that orbits with low pericenters, $\sim$ 10 kpc, produce disturbances that match the observed outer HI disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of $(\mu_{\alpha}cos\delta, \mu_{\delta}) = (-0.068,0.032) \pm (0.023,-0.031)~\rm mas/yr$, which agrees with results from Torrealba et al. (2019) within the errors, but gives lower mean pericenters, e.g., $\sim$ 15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis (Chakrabarti \& Blitz 2009) of the dwarf that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

The Milky Way is filled with the tidally-disrupted remnants of globular clusters and dwarf galaxies. Determining the properties of these objects -- in particular, initial masses and density profiles -- is relevant to both astronomy and dark matter physics. However, most direct measures of mass cannot be applied to tidal debris, as the systems of interest are no longer in equilibrium. Since phase-space density is conserved during adiabatic phase mixing, Liouville's theorem provides a connection between stellar kinematics as measured by observatories such as Gaia and the original mass of the disrupted system. Accurately recovering the phase-space density is complicated by uncertainties resulting from measurement errors and orbital integration, which both effectively inject entropy into the system, preferentially decreasing the measured density. In this paper, we demonstrate that these two issues can be overcome. First, we measure the phase-space density of the globular cluster M4 in Gaia data, and use Liouville's theorem to derive its mass. We then show that, for tidally disrupted systems, the orbital parameters and thus phase-space density can be inferred by minimizing the phase-space entropy of cold stellar streams. This work is therefore a proof of principle that true phase-space density can be measured and the original properties of the star cluster reconstructed in systems of astrophysical interest.

We present a model for the interaction of the GD-1 stellar stream with a massive perturber that naturally explains many of the observed stream features, including a gap and an off-stream spur of stars. The model involves an impulse by a fast encounter, after which the stream grows a loop of stars at different orbital energies. At specific viewing angles, this loop appears offset from the stream track. A quantitative comparison of the spur and gap features prefers models where the perturber is in the mass range of $10^6\,\rm M_\odot$ to $10^8\,\rm M_\odot$. Orbit integrations back in time show that the stream encounter could not have been caused by any known globular cluster or dwarf galaxy with a determined orbit, and mass, size and impact-parameter arguments show that it could not have been caused by a molecular cloud in the Milky Way disk. The most plausible explanation for the gap-and-spur structure is an encounter with a dark-matter substructure, like those predicted to populate galactic halos in LCDM cosmology. However, the expected densities of LCDM subhalos in this mass range and in this part of the Milky Way are $2-3\,\sigma$ lower than the inferred density of the GD-1 perturber. This observation opens up the possibility that detailed observations of streams could measure the mass spectrum of dark-matter substructures and even identify individual substructures and their orbits in the Galactic halo.

In simple models of the Milky Way, tidally disrupting satellites produce long and thin---nearly one-dimensional---stellar streams. Using astrometric data from the Gaia second data release and photometry from the Dark Energy Survey, we demonstrate that the Jhelum stream, a stellar stream in the inner halo, is a two-dimensional structure. The spatial distribution of highly probable Jhelum members reveals a dense thin component and an associated diffuse, spatially offset component. These two spatial components have indistinguishable proper motions (at $\sigma\sim1\,\rm mas\,yr^{-1}$ level) and a similar ratio of blue straggler to blue horizontal branch stars, which indicates a common origin for the two components. The best-fit orbit to the narrow component (pericenter $8\,\rm kpc$, apocenter $24\,\rm kpc$), however, does not explain the wide component of the Jhelum stream. On the other hand, an older orbital wrap of Jhelum's orbit traces the Indus stream, indicating a possible connection between these two structures and additional complexity in Jhelum's formation. Substructure in the Jhelum progenitor or precession of its tidal debris in the Milky Way potential may explain the observed structure of Jhelum. Future spectroscopic data will enable discrimination between these "nature" and "nurture" formation scenarios. Jhelum adds to the growing list of cold stellar streams that display complex morphologies and promise to reveal the dynamical history of the Milky Way.

In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various parameters of the IMBH-SMBH binary: IMBH mass, semimajor axis and eccentricity. We also calculate the rates of change of the BH binary orbital elements due to those stellar ejections. One of our new findings is that the ejection rate depends (although mildly) on the rotation of the stellar nucleus (its total angular momentum). We also compare the ejection velocity distribution with that produced by the Hills mechanism (stellar binary disruption) and find that the latter produces faster stars on average. Also, the IMBH mechanism produces an ejection velocity distribution which is flattened towards the BH binary plane while the Hills mechanism produces a spherically symmetric one. The results of this paper will allow us in the future to model the ejection of stars by an evolving BH binary and compare both models with \textit{Gaia} observations, for a wide variety of environments (galactic nuclei, globular clusters, the Large Magellanic Clouds, etc.).

In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Our uncertainties about binary star systems (and triples and so on) limit our capabilities in literally every single one of the Thematic Areas identified for Astro2020. We need to understand the population statistics of stellar multiplicity and their variations with stellar type, chemistry, and dynamical environment: Correct interpretation of any exoplanet experiment depends on proper treatment of resolved and unresolved binaries; stellar multiplicity is a direct outcome of star and companion formation; the most precise constraints on stellar structure come from well-characterized binary systems; stellar populations heavily rely on stellar and binary evolution modeling; high-redshift galaxy radiation and reionization is controlled by binary-dependent stellar physics; compact objects are the outcomes of binary evolution; the interpretation of multi-messenger astronomy from gravitational waves, light, and neutrinos relies on understanding the products of binary star evolution; near-Universe constraints on the Hubble constant with Type Ia supernovae and gravitational-wave mergers are subject to systematics related to their binary star progenitors; local measures of dark-matter substructure masses are distorted by binary populations. In order to realize the scientific goals in each of these themes over the next decade, we therefore need to understand how binary stars and stellar multiplets are formed and distributed in the space of masses, composition, age, and orbital properties, and how the distribution evolves with time. This white paper emphasizes the interdisciplinary importance of binary-star science and advocates that coordinated investment from all astrophysical communities will benefit almost all branches of astrophysics.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber.

Tidally-disrupted globular clusters are transformed into thin, dynamically-cold streams of stars that are extremely valuable tracers of the large- and small-scale distribution of mass in the Galaxy. Using data from the Gaia second data release combined with Pan-STARRS photometry, we present a sample of highly-probable members of the longest cold stream in the Milky Way, GD-1. The resulting map of GD-1: (1) extends the apparent length of the stream by $20^\circ$, (2) reveals plausible locations for the progenitor, (3) detects high-contrast gaps along the stream, and (4) indicates the existence of stream members perturbed off the main stream track. These discoveries are only possible because of the exquisite astrometry from Gaia, which permits a clean separation of the stream from Milky Way stars. The additional length and a proper treatment of the progenitor will aid in dynamical modeling of GD-1 for mapping the large-scale dark matter distribution. The complex morphology of the stream points to a turbulent history; detailed phase-space properties of the perturbed stream members could potentially constrain dark matter substructure in the Milky Way.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

We present 1,201 galaxies at $0.05<z<0.45$ that host tidal features, detected from the first $\sim\! 200$ deg$^2$ of imaging from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). All galaxies in the present sample have spectroscopic observations from the Sloan Digital Sky Survey (SDSS) spectroscopic campaigns, generating a sample of 21208 galaxies. Of these galaxies, we identify 214 shell systems and 987 stream systems. For 575 of these systems, we are additionally able to measure the $(g-i)$ colors of the tidal features. We find evidence for star formation in a subset of the streams, with the exception of streams around massive ellipticals, and find that stream host galaxies span the full range of stellar masses in our sample. Galaxies which host shells are predominantly red and massive: we find that observable shells form more frequently around ellipticals than around disc galaxies of the same stellar mass. Although the majority of the shells in our sample are consistent with being formed by minor mergers, $15\% \pm 4.4\%$ of shell host galaxies have $(g-i)$ colors as red as their host galaxy, consistent with being formed by major mergers. These "red shells" are additionally preferentially aligned with the major axis of the host galaxy, as previously predicted from simulations. We suggest that although the bulk of the observable shell population originates from fairly minor mergers, which preferentially form shells that are not aligned with the major axis of the galaxy, major mergers produce a significant number of observable shells.

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.

Short-period binary star systems dissipate orbital energy through tidal interactions that lead to tighter, more circular orbits. When at least one star in a binary has evolved off of the main sequence, orbital circularization occurs for longer-period (~100 day) systems. Past work by (Verbunt & Phinney 1995) has shown that the orbital parameters and the circularization periods of a small sample of binary stars with evolved-star members can be understood within the context of standard tidal circularization theory. Using a sample of binaries with subgiant, giant, and red clump star members that is nearly an order of magnitude larger, we reexamine predictions for tidal circularization of binary stars with evolved members. We confirm that systems predicted by equilibrium-tide theory to have circular orbits generally have negligible measured eccentricities. The circularization period is correlated with the surface gravity (i.e. size) of the evolved member, indicating that the circularization timescale must be shorter than the evolutionary timescale along the giant branch. A few exceptions to the conclusions above are mentioned in the discussion: Some of these exceptions are likely systems in which the spectrum of the secondary biases the radial velocity measurements, but four appear to be genuine, short-period, moderate-eccentricity systems.

Multi-epoch radial velocity measurements of stars can be used to identify stellar, sub-stellar, and planetary-mass companions. Even a small number of observation epochs can be informative about companions, though there can be multiple qualitatively different orbital solutions that fit the data. We have custom-built a Monte Carlo sampler (The Joker) that delivers reliable (and often highly multi-modal) posterior samplings for companion orbital parameters given sparse radial-velocity data. Here we use The Joker to perform a search for companions to 96,231 red-giant stars observed in the APOGEE survey (DR14) with $\geq 3$ spectroscopic epochs. We select stars with probable companions by making a cut on our posterior belief about the amplitude of the stellar radial-velocity variation induced by the orbit. We provide (1) a catalog of 320 companions for which the stellar companion properties can be confidently determined, (2) a catalog of 4,898 stars that likely have companions, but would require more observations to uniquely determine the orbital properties, and (3) posterior samplings for the full orbital parameters for all stars in the parent sample. We show the characteristics of systems with confidently determined companion properties and highlight interesting systems with candidate compact object companions.

Our Galaxy is thought to have undergone an active evolutionary history dominated by star formation, the accretion of cold gas, and, in particular, mergers up to 10 gigayear ago. The stellar halo reveals rich fossil evidence of these interactions in the form of stellar streams, substructures, and chemically distinct stellar components. The impact of dwarf galaxy mergers on the content and morphology of the Galactic disk is still being explored. Recent studies have identified kinematically distinct stellar substructures and moving groups, which may have extragalactic origin. However, there is mounting evidence that stellar overdensities at the outer disk/halo interface could have been caused by the interaction of a dwarf galaxy with the disk. Here we report detailed spectroscopic analysis of 14 stars drawn from two stellar overdensities, each lying about 5 kiloparsecs above and below the Galactic plane - locations suggestive of association with the stellar halo. However, we find that the chemical compositions of these stars are almost identical, both within and between these groups, and closely match the abundance patterns of the Milky Way disk stars. This study hence provides compelling evidence that these stars originate from the disk and the overdensities they are part of were created by tidal interactions of the disk with passing or merging dwarf galaxies.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

The Monoceros Ring (also known as the Galactic Anticenter Stellar Structure) and A13 are stellar overdensities at estimated heliocentric distances of $d \sim 11$ kpc and 15 kpc observed at low Galactic latitudes towards the anticenter of our Galaxy. While these overdensities were initially thought to be remnants of a tidally-disrupted satellite galaxy, an alternate scenario is that they are composed of stars from the Milky Way (MW) disk kicked out to their current location due to interactions between a satellite galaxy and the disk. To test this scenario, we study the stellar populations of the Monoceros Ring and A13 by measuring the number of RR Lyrae and M giant stars associated with these overdensities. We obtain low-resolution spectroscopy for RR Lyrae stars in the two structures and measure radial velocities to compare with previously measured velocities for M giant stars in the regions of the Monoceros Ring and A13, to assess the fraction of RR Lyrae to M giant stars ($f_{RR:MG}$) in A13 and Mon/GASS. We perform velocity modeling on 153 RR Lyrae stars (116 in the Monoceros Ring and 37 in A13) and find that both structures have very low $f_{RR:MG}$. The results support a scenario in which stars in A13 and Mon/GASS formed in the MW disk. We discuss a possible association between Mon/GASS, A13, and the Triangulum-Andromeda overdensity based on their similar velocity distributions and $f_{RR:MG}$.

We present astroplan - an open source, open development, Astropy affiliated package for ground-based observation planning and scheduling in Python. astroplan is designed to provide efficient access to common observational quantities such as celestial rise, set, and meridian transit times and simple transformations from sky coordinates to altitude-azimuth coordinates without requiring a detailed understanding of astropy's implementation of coordinate systems. astroplan provides convenience functions to generate common observational plots such as airmass and parallactic angle as a function of time, along with basic sky (finder) charts. Users can determine whether or not a target is observable given a variety of observing constraints, such as airmass limits, time ranges, Moon illumination/separation ranges, and more. A selection of observation schedulers are included which divide observing time among a list of targets, given observing constraints on those targets. Contributions to the source code from the community are welcome.

Collisions and interactions between gas-rich galaxies are thought to be pivotal stages in their formation and evolution, causing the rapid production of new stars, and possibly serving as a mechanism for fueling supermassive black holes (BH). Harnessing the exquisite spatial resolution (~0.5 arcsec) afforded by the first ~170 deg^2 of the Hyper Suprime-Cam (HSC) Survey, we present our new constraints on the importance of galaxy-galaxy major mergers (1:4) in growing BHs throughout the last ~8 Gyrs. Utilizing mid-infrared observations in the WISE All-Sky survey, we robustly select active galactic nuclei (AGN) and mass-matched control galaxy samples, totaling ~140,000 spectroscopically confirmed systems at i<22 mag. We identify galaxy interaction signatures using a novel machine-learning random forest decision tree technique allowing us to select statistically significant samples of major-mergers, minor-mergers/irregular-systems, and non-interacting galaxies. We use these samples to show that galaxies undergoing mergers are a factor ~2-7 more likely to contain luminous obscured AGN than non-interacting galaxies, and this is independent of both stellar mass and redshift to z < 0.9. Furthermore, based on our comparison of AGN fractions in mass-matched samples, we determine that the most luminous AGN population (L_AGN > 10^45 erg/s) systematically reside in merging systems over non-interacting galaxies. Our findings show that galaxy-galaxy interactions do, on average, trigger luminous AGN activity substantially more often than in secularly evolving non-interacting galaxies, and we further suggest that the BH growth rate may be closely tied to the dynamical time of the merger system.

Stellar tidal streams provide an opportunity to study the motion and structure of the disrupting galaxy as well as the gravitational potential of its host. Streams around the Milky Way are especially promising as phase space positions of individual stars will be measured by ongoing or upcoming surveys. Nevertheless, it remains a challenge to accurately assess distances to stars farther than 10 kpc from the Sun, where we have the poorest knowledge of the Galaxy's mass distribution. To address this we present observations of 32 candidate RR Lyrae stars in the Orphan tidal stream taken as part of the Spitzer Merger History and Shape of the Galactic Halo (SMHASH) program. The extremely tight correlation between the periods, luminosities, and metallicities of RR Lyrae variable stars in the Spitzer IRAC $\mathrm{3.6 \mu m}$ band allows the determination of precise distances to individual stars; the median statistical distance uncertainty to each RR Lyrae star is $2.5\%$. By fitting orbits in an example potential we obtain an upper limit on the mass of the Milky Way interior to 60 kpc of $\mathrm{5.6_{-1.1}^{+1.2}\times 10^{11}\ M_\odot}$, bringing estimates based on the Orphan Stream in line with those using other tracers. The SMHASH data also resolve the stream in line--of--sight depth, allowing a new perspective on the internal structure of the disrupted dwarf galaxy. Comparing with N--body models we find that the progenitor had an initial dark halo mass of approximately $\mathrm{3.2 \times 10^{9}\ M_\odot}$, placing the Orphan Stream's progenitor amongst the classical dwarf spheroidals.

Recent Pan-STARRS data show that the leading arm from the globular cluster Palomar 5 (Pal 5) appears shorter than the trailing arm, while simulations of Pal 5 predict similar angular extents. We demonstrate that including the spinning Galactic bar with appropriate pattern speeds in the dynamical modeling of Pal 5 can reproduce the Pan-STARRS data. As the bar sweeps by, some stream stars experience a difference in net torques near pericenter. This leads to the formation of apparent gaps along Pal 5's tidal streams and these gaps grow due to an energy offset from the rest of the stream members. We conclude that only streams orbiting far from the Galactic center or streams on retrograde orbits (with respect to the bar) can be used to unambiguously constrain dark matter subhalo interactions. Additionally, we expect that the Pal 5 leading arm debris should re-appear south of the Pan-STARRS density truncation.

We report and discuss the discovery of a comoving pair of bright solar-type stars, HD 240430 and HD 240429, with a significant difference in their chemical abundances. The two stars have an estimated 3D separation of $\approx 0.6$ pc ($\approx 0.01$ pc projected) at a distance of $r\approx 100$ pc with nearly identical three-dimensional velocities, as inferred from Gaia TGAS parallaxes and proper motions, and high-precision radial velocity measurements. Stellar parameters determined from high-resolution Keck HIRES spectra indicate that both stars are $\sim 4$ Gyr old. The more metal-rich of the two, HD 240430, shows an enhancement of refractory ($T_C>1200$ K) elements by $\approx 0.2$ dex and a marginal enhancement of (moderately) volatile elements ($T_C<1200$ K, C, N, O, Na, and Mn). This is the largest metallicity difference found in a wide binary pair yet. Additionally, HD 240430 shows an anomalously high surface lithium abundance ($A(\mathrm{Li})=2.75$), higher than its companion by $0.5$ dex. The proximity in phase-space and ages between the two stars suggests that they formed together with the same composition, at odds with the observed differences in metallicity and abundance patterns. We therefore suggest that the star HD~240430, "Kronos", accreted 15 $M_\oplus$ of rocky material after birth, selectively enhancing the refractory elements as well as lithium in its surface and convective envelope.

Members of widely-separated stellar binaries (> 0.1 pc) or loosely bound stellar associations and clusters are born with small spreads in kinematics and chemistry, but disrupt and dissociate on timescales comparable to their orbital time within the Milky Way. The kinematic properties of comoving star pairs---bound wide binaries or unbound but coeval stars---are therefore sensitive to both the smooth tidal field and gravitational perturbation spectrum of the Milky Way. In previous work, we identified >4,000 candidate comoving pairs using only astrometric data from the Tycho-Gaia Astrometric Solution (TGAS) catalog and find a large number of pairs with separations > 1 pc. Without radial velocity measurements for the vast majority of these pairs, the false-positive rate at these separations is significant (~50%). Here we present results from our own low-resolution radial-velocity survey of a random sample of the candidate pairs to identify and validate true comoving star pairs. Of the 311 observed comoving pairs, we confirm 127 comoving pairs with separations as large as ~10 pc, our original search limit. Between separations of 0.001 pc to ~0.5 pc, the number of confirmed pairs per separation decreases with increasing separation (number per log separation is uniform). From ~0.5 pc to 10 pc, the number per separation is approximately uniform (number per log separation increases). We confirm the discovery of a population of comoving star pairs at very large separations, suggesting that disrupted wide binaries and stellar associations remain approximately spatially coherent.

Deep photometric surveys of the Milky Way have revealed diffuse structures encircling our Galaxy far beyond the "classical" limits of the stellar disk. This paper reviews results from our own and other observational programs, which together suggest that, despite their extreme positions, the stars in these structures were formed in our Galactic disk. Mounting evidence from recent observations and simulations implies kinematic connections between several of these distinct structures. This suggests the existence of collective disk oscillations that can plausibly be traced all the way to asymmetries seen in the stellar velocity distribution around the Sun. There are multiple interesting implications of these findings: they promise new perspectives on the process of disk heating, they provide direct evidence for a stellar halo formation mechanism in addition to the accretion and disruption of satellite galaxies, and, they motivate searches of current and near-future surveys to trace these oscillations across the Galaxy. Such maps could be used as dynamical diagnostics in the emerging field of "Galactoseismology", which promises to model the history of interactions between the Milky Way and its entourage of satellites, as well examine the density of our dark matter halo. As sensitivity to very low surface brightness features around external galaxies increases, many more examples of such disk oscillations will likely be identified. Statistical samples of such features not only encode detailed information about interaction rates and mergers, but also about long sought-after dark matter halo densities and shapes. Models for the Milky Way's own Galactoseismic history will therefore serve as a critical foundation for studying the weak dynamical interactions of galaxies across the universe.

Thanks to modern sky surveys, over twenty stellar streams and overdensity structures have been discovered in the halo of the Milky Way. In this paper, we present an analysis of spectroscopic observations of individual stars from one such structure, "A13", first identified as an overdensity using the M giant catalog from the Two Micron All-Sky Survey. Our spectroscopic observations show that stars identified with A13 have a velocity dispersion of $\lesssim$ 40 $\mathrm{km~s^{-1}}$, implying that it is a genuine coherent structure rather than a chance super-position of random halo stars. From its position on the sky, distance ($\sim$15~kpc heliocentric), and kinematical properties, A13 is likely to be an extension of another low Galactic latitude substructure -- the Galactic Anticenter Stellar Structure (also known as the Monoceros Ring) -- towards smaller Galactic longitude and farther distance. Furthermore, the kinematics of A13 also connect it with another structure in the southern Galactic hemisphere -- the Triangulum-Andromeda overdensity. We discuss these three connected structures within the context of a previously proposed scenario that one or all of these features originate from the disk of the Milky Way.

Converting a noisy parallax measurement into a posterior belief over distance requires inference with a prior. Usually this prior represents beliefs about the stellar density distribution of the Milky Way. However, multi-band photometry exists for a large fraction of the \textsl{\small{Gaia}} \textsl{\small{TGAS}} Catalog and is incredibly informative about stellar distances. Here we use \textsl{\small{2MASS}} colors for 1.4 million \textsl{\small{TGAS}} stars to build a noise-deconvolved empirical prior distribution for stars in color--magnitude space. This model contains no knowledge of stellar astrophysics or the Milky Way, but is precise because it accurately generates a large number of noisy parallax measurements under an assumption of stationarity; that is, it is capable of combining the information from many stars. We use the Extreme Deconvolution (\textsl{\small{XD}}) algorithm---an Empirical Bayes approximation to a full hierarchical model of the true parallax and photometry of every star---to construct this prior. The prior is combined with a \textsl{\small{TGAS}} likelihood to infer a precise photometric parallax estimate and uncertainty (and full posterior) for every star. Our parallax estimates are more precise than the \textsl{\small{TGAS}} catalog entries by a median factor of 1.2 (14% are more precise by a factor >2) and are more precise than previous Bayesian distance estimates that use spatial priors. We validate our parallax inferences using members of the Milky Way star cluster M67, which is not visible as a cluster in the \textsl{\small{TGAS}} parallax estimates, but appears as a cluster in our posterior parallax estimates. Our results, including a parallax posterior pdf for each of 1.4 million \textsl{\small{TGAS}} stars, are available in companion electronic tables.

The primary sample of the {\it Gaia} Data Release 1 is the Tycho-Gaia Astrometric Solution (TGAS): $\approx$ 2 million Tycho-2 sources with improved parallaxes and proper motions relative to the initial catalog. This increased astrometric precision presents an opportunity to find new binary stars and moving groups. We search for high-confidence comoving pairs of stars in TGAS by identifying pairs of stars consistent with having the same 3D velocity using a marginalized likelihood ratio test to discriminate candidate comoving pairs from the field population. Although we perform some visualizations using (bias- corrected) inverse parallax as a point estimate of distance, the likelihood ratio is computed with a probabilistic model that includes the covariances of parallax and proper motions and marginalizes the (unknown) true distances and 3D velocities of the stars. We find 13,085 comoving star pairs among 10,606 unique stars with separations as large as 10 pc (our search limit). Some of these pairs form larger groups through mutual comoving neighbors: many of these pair networks correspond to known open clusters and OB associations, but we also report the discovery of several new comoving groups. Most surprisingly, we find a large number of very wide ($>1$ pc) separation comoving star pairs, the number of which increases with increasing separation and cannot be explained purely by false-positive contamination. Our key result is a catalog of high-confidence comoving pairs of stars in TGAS. We discuss the utility of this catalog for making dynamical inferences about the Galaxy, testing stellar atmosphere models, and validating chemical abundance measurements.

We report the discovery of a diffuse stellar cloud with an angular extent $\gtrsim30^{\prime\prime}$, which we term "Sumo Puff", in data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). While we do not have a redshift for this object, it is in close angular proximity to a post-merger galaxy at redshift $z=0.0431$ and is projected within a few virial radii (assuming similar redshifts) of two other ${\sim}L_\star$ galaxies, which we use to bracket a potential redshift range of $0.0055 < z < 0.0431$. The object's light distribution is flat, as characterized by a low Sersic index ($n\sim0.3$). It has a low central $g$-band surface brightness of ${\sim}26.4$ mag arcsec$^{-2}$, large effective radius of ${\sim}13^{\prime\prime}$ (${\sim}11$ kpc at $z=0.0431$ and ${\sim}1.5$ kpc at $z=0.0055$), and an elongated morphology ($b/a\sim0.4$). Its red color ($g-i\sim1$) is consistent with a passively evolving stellar population and similar to the nearby post-merger galaxy, and we may see tidal material connecting Sumo Puff with this galaxy. We offer two possible interpretations for the nature of this object: (1) it is an extreme, galaxy-size tidal feature associated with a recent merger event, or (2) it is a foreground dwarf galaxy with properties consistent with a quenched, disturbed ultra-diffuse galaxy. We present a qualitative comparison with simulations that demonstrates the feasibility of forming a structure similar to this object in a merger event. Follow-up spectroscopy and/or deeper imaging to confirm the presence of the bridge of tidal material will be necessary to reveal the true nature of this object.

The Ophiuchus stellar stream is peculiar: (1) its length is short given the age of its constituent stars, and (2) several probable member stars that lie close in both sky position and velocity have dispersions in these dimensions that far exceed those seen within the stream. The stream's proximity to the Galactic center suggests that the bar must have a significant influence on its dynamical history: The triaxiality and time-dependence of the bar may generate chaotic orbits in the vicinity of the stream that can greatly affect its morphology. We explore this hypothesis with models of stream formation along orbits consistent with Ophiuchus' properties in a Milky Way potential model that includes a rotating bar. We find that in all choices for the rotation parameters of the bar, orbits fit to the stream are strongly chaotic. Mock streams generated along these orbits qualitatively match the observed properties of the stream: because of chaos, stars stripped early generally form low-density, high-dispersion "fans" leaving only the most recently disrupted material detectable as a strong over-density. Our models predict that there should be more low-surface-brightness tidal debris than detected so far, likely with a complex phase-space morphology. The existence of or lack of these features around the Ophiuchus stream would provide an interesting constraint on the properties of the Milky Way bar and would help distinguish between formation scenarios for the stream. This is the first time that chaos has been used to explain the properties of a stellar stream and is the first demonstration of the dynamical importance of chaos in the Galactic halo. The existence of long, thin streams around the Milky Way---presumably formed along non- or weakly-chaotic orbits---may represent only a subset of the total population of disrupted satellites.

Given sparse or low-quality radial-velocity measurements of a star, there are often many qualitatively different stellar or exoplanet companion orbit models that are consistent with the data. The consequent multimodality of the likelihood function leads to extremely challenging search, optimization, and MCMC posterior sampling over the orbital parameters. Here we create a custom Monte Carlo sampler for sparse or noisy radial-velocity measurements of two-body systems that can produce posterior samples for orbital parameters even when the likelihood function is poorly behaved. The six standard orbital parameters for a binary system can be split into four non-linear parameters (period, eccentricity, argument of pericenter, phase) and two linear parameters (velocity amplitude, barycenter velocity). We capitalize on this by building a sampling method in which we densely sample the prior pdf in the non-linear parameters and perform rejection sampling using a likelihood function marginalized over the linear parameters. With sparse or uninformative data, the sampling obtained by this rejection sampling is generally multimodal and dense. With informative data, the sampling becomes effectively unimodal but too sparse: in these cases we follow the rejection sampling with standard MCMC. The method produces correct samplings in orbital parameters for data that include as few as three epochs. The Joker can therefore be used to produce proper samplings of multimodal pdfs, which are still informative and can be used in hierarchical (population) modeling. We give some examples that show how the posterior pdf depends sensitively on the number and time coverage of the observations and their uncertainties.

Pulsating stars, such as Cepheids, Miras, and RR Lyrae stars, are important distance indicators and calibrators of the "cosmic distance ladder", and yet their period-luminosity-metallicity (PLZ) relations are still constrained using simple statistical methods that cannot take full advantage of available data. To enable optimal usage of data provided by the Gaia mission, we present a probabilistic approach that simultaneously constrains parameters of PLZ relations and uncertainties in Gaia parallax measurements. We demonstrate this approach by constraining PLZ relations of type $ab$ RR Lyrae stars in near-infrared W1 and W2 bands, using Tycho-Gaia Astrometric Solution (TGAS) parallax measurements for a sample of $\approx100$ type $ab$ RR Lyrae stars located within 2.5 kpc of the Sun. The fitted PLZ relations are consistent with previous studies, and in combination with other data, deliver distances precise to 6% (once various sources of uncertainty are taken into account). To a precision of 0.05 mas ($1\sigma$), we do not find a statistically significant offset in TGAS parallaxes for this sample of distant RR Lyrae stars (median parallax of 0.8 mas and distance of 1.4 kpc). With only minor modifications, our probabilistic approach can be used to constrain PLZ relations of other pulsating stars, and we intend to apply it to Cepheid and Mira stars in the near future.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Supermassive black hole binaries (SMBHBs) at sub-parsec separations should be common in galactic nuclei, as a result of frequent galaxy mergers. Hydrodynamical simulations of circumbinary discs predict strong periodic modulation of the mass accretion rate on time-scales comparable to the orbital period of the binary. As a result, SMBHBs may be recognized by the periodic modulation of their brightness. We conducted a statistical search for periodic variability in a sample of 35,383 spectroscopically confirmed quasars in the photometric database of the Palomar Transient Factory (PTF). We analysed Lomb-Scargle periodograms and assessed the significance of our findings by modeling each individual quasar's variability as a damped random walk (DRW). We identified 50 quasars with significant periodicity beyond the DRW model, typically with short periods of a few hundred days. We find 33 of these to remain significant after a re-analysis of their periodograms including additional optical data from the intermediate-PTF and the Catalina Real-Time Transient Survey (CRTS). Assuming that the observed periods correspond to the redshifted orbital periods of SMBHBs, we conclude that our findings are consistent with a population of unequal-mass SMBHBs, with a typical mass ratio as low as q = M2/M1 ~ 0.01.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DM*H from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by f{\sigma}8(z), from redshift-space distortions (RSD). We combine measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method. Combined with Planck 2015 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature {\Omega}_K =0.0003+/-0.0026 and a dark energy equation of state parameter w = -1.01+/-0.06, in strong affirmation of the spatially flat cold dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3+/-1.0 km/s/Mpc even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8+/-1.2 km/s/Mpc. Assuming flat {\Lambda}CDM we find {\Omega}_m = 0.310+/-0.005 and H0 = 67.6+/-0.5 km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass sum.

Several long, dynamically cold stellar streams have been observed around the Milky Way Galaxy, presumably formed from the tidal disruption of globular clusters. In integrable potentials---where all orbits are regular---tidal debris phase-mixes close to the orbit of the progenitor system. However, the Milky Way's dark matter halo is expected not to be fully integrable; an appreciable fraction of orbits will be chaotic. This paper examines the influence of chaos on the phase-space morphology of cold tidal streams. Streams even in weakly chaotic regions look very different from those in regular regions. We find that streams can be sensitive to chaos on a much shorter time-scale than any standard prediction (from the Lyapunov or frequency-diffusion times). For example, on a weakly chaotic orbit with a chaotic timescale predicted to be $>$1000 orbital periods ($>$1000 Gyr), the resulting stellar stream is, after just a few 10's of orbits, substantially more diffuse than any formed on a nearby but regular orbit. We find that the enhanced diffusion of the stream stars can be understood by looking at the variance in orbital frequencies of orbit ensembles centered around the parent (progenitor) orbit. Our results suggest that long, cold streams around our Galaxy must exist only on regular (or very nearly regular) orbits; they potentially provide a map of the regular regions of the Milky Way potential. This suggests a promising new direction for the use of tidal streams to constrain the distribution of dark matter around our Galaxy.

This paper presents an example where the morphology of a single stellar stream can be used to rule out a specific galactic potential form without the need for velocity information. We investigate the globular cluster Palomar5 (Pal 5), which is tidally disrupting into a cold, thin stream mapped over 22 degrees on the sky with a typical width of 0.7 degrees. We generate models of this stream by fixing Pal 5's present-day position, distance and radial velocity via observations, while allowing its proper motion to vary. In a spherical dark matter halo we easily find models that fit the observed morphology. However, no plausible Pal 5 model could be found in the triaxial potential of Law & Majewski (2010), which has been proposed to explain the properties of the Sagittarius stream. In this case, the long, thin and curved morphology of the Pal5 stream alone can be used to rule out such a potential configuration. Pal 5 like streams in this potential are either too straight, missing the curvature of the observations, or show an unusual morphology which we dub stream-fanning: a signature sensitive to the triaxiality of a potential. We conclude that the mere existence of other thin tidal streams must provide broad constraints on the orientation and shape of the dark matter halo they inhabit.

The Ophiuchus stellar stream presents a dynamical puzzle: its old stellar populations ($\sim 12$ Gyr) cannot be reconciled with (1) its orbit in a simple model for the Milky Way potential and (2) its short angular extent, both of which imply that the observed stream formed within the last $<1$ Gyr. Recent theoretical work has shown that streams on chaotic orbits may abruptly fan out near their apparent ends; stars in these fans are dispersed in both position and velocity and may be difficult to associate with the stream. Here we present the first evidence of such stream-fanning in the Ophiuchus stream, traced by four blue horizontal branch (BHB) stars beyond the apparent end of the stream. These stars stand out from the background by their high velocities ($v_{\rm los} > 230$ km s$^{-1}$) against $\sim 40$ other stars: their velocities are comparable to those of the stream, but would be exceptional if they were unrelated halo stars. Their positions and velocities are, however, inconsistent with simple extrapolation of the observed cold, high-density portion of the stream. These observations suggest that stream-fanning may be a real, observable effect and, therefore, that Ophiuchus may be on a chaotic orbit. They also show that the Ophiuchus stream is more extended and hence dynamically older than previously thought, easing the stellar population vs. dynamical age tension.

The Triangulum-Andromeda stellar clouds (TriAnd1 and TriAnd2) are a pair of concentric ring- or shell-like over-densities at large $R$ ($\approx$ 30 kpc) and $Z$ ($\approx$ -10 kpc) in the Galactic halo that are thought to have been formed from the accretion and disruption of a satellite galaxy. This paper critically re-examines this formation scenario by comparing the number ratio of RR Lyrae to M giant stars associated with the TriAnd clouds with other structures in the Galaxy. The current data suggest a stellar population for these over-densities ($f_{\rm RR:MG} < 0.38$ at 95% confidence) quite unlike any of the known satellites of the Milky Way ($f_{\rm RR:MG} \approx 0.5$ for the very largest and $f_{\rm RR:MG} >>1$ for the smaller satellites) and more like the population of stars born in the much deeper potential well inhabited by the Galactic disk ($f_{\rm RR:MG} < 0.01$). N-body simulations of a Milky-Way-like galaxy perturbed by the impact of a dwarf galaxy demonstrate that, in the right circumstances, concentric rings propagating outwards from that Galactic disk can plausibly produce similar over-densities. These results provide dramatic support for the recent proposal by Xu et al. (2015) that, rather than stars accreted from other galaxies, the TriAnd clouds could represent stars kicked-out from our own disk. If so, these would be the first populations of disk stars to be found in the Galactic halo and a clear signature of the importance of this second formation mechanism for stellar halos more generally. Moreover, their existence at the very extremities of the disk places strong constraints on the nature of the interaction that formed them.

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT telescopes. We confirm the stream as a fast moving ($v_{los}\sim290$ km s$^{-1}$), kinematically cold group ($\sigma_{v_{los}}\lesssim1$ km s$^{-1}$) of $\alpha$-enhanced and metal-poor stars (${\rm [\alpha/Fe]\sim0.4}$ dex, ${\rm [Fe/H]\sim-2.0}$ dex). Using a probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color-magnitude, coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find that that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than wide, merely appearing highly foreshortened in projection. The analysis of the stellar population contained in the stream suggests that it is $\sim12$ Gyr old, and that its initial stellar mass was $\sim2\times10^4$ $M_{\odot}$ (or at least $\gtrsim7\times10^3$ $M_{\odot}$). Assuming a fiducial Milky Way potential, we fit an orbit to the stream which matches the observed phase-space distribution, except for some tension in the proper motions: the stream has an orbital period of $\sim350$ Myr, and is on a fairly eccentric orbit ($e\sim0.66$) with a pericenter of $\sim3.5$ kpc and an apocenter of $\sim17$ kpc. The phase-space structure and stellar population of the stream show that its progenitor must have been a globular cluster that was disrupted only $\sim240$ Myr ago. We do not detect any significant overdensity of stars along the stream that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the progenitor.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

Graham et al. (2015) discovered a supermassive black hole binary (SMBHB) candidate and identified the detected 5.2 yr period of the optical variability as the orbital period of the binary. Hydrodynamical simulations predict multiple periodic components for the variability of SMBHBs, thus raising the possibility that the true period of the binary is different from 5.2 yr. We analyse the periodogram of PG1302 and find no compelling evidence for additional peaks. We derive upper limits on any additional periodic modulations in the available data, by modeling the light-curve as the sum of red noise and the known 5.2 yr periodic component, and injecting additional sinusoidal signals. We find that, with the current data, we would be able to detect with high significance (false alarm probability <1%) secondary periodic terms with periods in the range predicted by the simulations, if the amplitude of the variability was at least ~0.07 mag (compared to 0.14 mag for the main peak). A three-year follow-up monitoring campaign with weekly observations can increase the sensitivity for detecting secondary peaks ~50%, and would allow a more robust test of predictions from hydrodynamical simulations.

Measuring the masses of companions to single-line spectroscopic binary stars is (in general) not possible because of the unknown orbital plane inclination. Even when the mass of the visible star can be measured, only a lower limit can be placed on the mass of the unseen companion. However, since these inclination angles should be isotropically distributed, for a large enough, unbiased sample, the companion mass distribution can be deconvolved from the distribution of observables. In this work, we construct a hierarchical probabilistic model to infer properties of unseen companion stars given observations of the orbital period and projected radial velocity of the primary star. We apply this model to three mock samples of low-mass white dwarfs (LMWDs, $M\lesssim0.45~\Msun$) and a sample of post-common-envelope binaries. We use a mixture of two Gaussians to model the WD and neutron star (NS) companion mass distributions. Our model successfully recovers the initial parameters of these test data sets. We then apply our model to 55 WDs in the extremely low-mass (ELM) WD Survey. Our maximum a posteriori model for the WD companion population has a mean mass $\mu_{\rm WD} = 0.74~\Msun$, with a standard deviation $\sigma_{\rm WD} = 0.24~\Msun$. Our model constrains the NS companion fraction $f_{\rm NS}$ to be $<$16\% at 68\% confidence. We make samples from the posterior distribution publicly available so that future observational efforts may compute the NS probability for newly discovered LMWDs.

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

The dark matter halo of the Milky Way is expected to be triaxial and filled with substructure. It is hoped that streams or shells of stars produced by tidal disruption of stellar systems will provide precise measures of the gravitational potential to test these predictions. We develop a method for inferring the Galactic potential with tidal streams based on the idea that the stream stars were once close in phase space. Our method can flexibly adapt to any form for the Galactic potential: it works in phase-space rather than action-space and hence relies neither on our ability to derive actions nor on the integrability of the potential. Our model is probabilistic, with a likelihood function and priors on the parameters. The method can properly account for finite observational uncertainties and missing data dimensions. We test our method on synthetic datasets generated from N-body simulations of satellite disruption in a static, multi-component Milky Way including a triaxial dark matter halo with observational uncertainties chosen to mimic current and near-future surveys of various stars. We find that with just four well-measured stream stars, we can infer properties of a triaxial potential with precisions of order 5-7 percent. Without proper motions we obtain 15 percent constraints on potential parameters and precisions around 25 percent for recovering missing phase-space coordinates. These results are encouraging for the eventual goal of using flexible, time-dependent potential models combined with larger data sets to unravel the detailed shape of the dark matter distribution around the Milky Way.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

Many photometric time-domain surveys are driven by specific goals, such as searches for supernovae or transiting exoplanets, which set the cadence with which fields are re-imaged. In the case of the Palomar Transient Factory (PTF), several sub-surveys are conducted in parallel, leading to non-uniform sampling over its $\sim$$20,000 \mathrm{deg}^2$ footprint. While the median $7.26 \mathrm{deg}^2$ PTF field has been imaged $\sim$40 times in \textit{R}-band, $\sim$$2300 \mathrm{deg}^2$ have been observed $>$100 times. We use PTF data to study the trade-off between searching for microlensing events in a survey whose footprint is much larger than that of typical microlensing searches, but with far-from-optimal time sampling. To examine the probability that microlensing events can be recovered in these data, we test statistics used on uniformly sampled data to identify variables and transients. We find that the von Neumann ratio performs best for identifying simulated microlensing events in our data. We develop a selection method using this statistic and apply it to data from fields with $>$10 $R$-band observations, $1.1\times10^9$ light curves, uncovering three candidate microlensing events. We lack simultaneous, multi-color photometry to confirm these as microlensing events. However, their number is consistent with predictions for the event rate in the PTF footprint over the survey's three years of operations, as estimated from near-field microlensing models. This work can help constrain all-sky event rate predictions and tests microlensing signal recovery in large data sets, which will be useful to future time-domain surveys, such as that planned with the Large Synoptic Survey Telescope.

Near-future data from ESA's Gaia mission will provide precise, full phase-space information for hundreds of millions of stars out to heliocentric distances of ~10 kpc. This "horizon" for full phase-space measurements is imposed by the Gaia parallax errors degrading to worse than 10%, and could be significantly extended by an accurate distance indicator. Recent work has demonstrated how Spitzer observations of RR Lyrae stars can be used to make distance estimates accurate to 2%, effectively extending the Gaia, precise-data horizon by a factor of ten in distance and a factor of 1000 in volume. This Letter presents one approach to exploit data of such accuracy to measure the Galactic potential using small samples of stars associated with debris from satellite destruction. The method is tested with synthetic observations of 100 stars from the end point of a simulation of satellite destruction: the shape, orientation, and depth of the potential used in the simulation are recovered to within a few percent. The success of this simple test with such a small sample in a single debris stream suggests that constraints from multiple streams could be combined to examine the Galaxy's dark matter halo in even more detail --- a truly unique opportunity that is enabled by the combination of Spitzer and Gaia with our intimate perspective on the Galaxy.

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We describe a procedure for background subtracting Sloan Digital Sky Survey (SDSS) imaging that improves the resulting detection and photometry of large galaxies on the sky. Within each SDSS drift scan run, we mask out detected sources and then fit a smooth function to the variation of the sky background. This procedure has been applied to all SDSS-III Data Release 8 images, and the results are available as part of that data set. We have tested the effect of our background subtraction on the photometry of large galaxies by inserting fake galaxies into the raw pixels, reanalyzing the data, and measuring them after background subtraction. Our technique results in no size-dependent bias in galaxy fluxes up to half-light radii of 100 arcsec; in contrast, for galaxies of that size the standard SDSS photometric catalog underestimates fluxes by about 1.5 mag. Our results represent a substantial improvement over the standard SDSS catalog results and should form the basis of any analysis of nearby galaxies using the SDSS imaging data.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

Computer representations of real numbers are necessarily discrete, with some finite resolution, discreteness, quantization, or minimum representable difference. We perform astrometric and photometric measurements on stars and co-add multiple observations of faint sources to demonstrate that essentially all of the scientific information in an optical astronomical image can be preserved or transmitted when the minimum representable difference is a factor of two finer than the root-variance of the per-pixel noise. Adopting a representation this coarse reduces bandwidth for data acquisition, transmission, or storage, or permits better use of the system dynamic range, without sacrificing any information for down-stream data analysis, including information on sources fainter than the minimum representable difference itself.