Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

By combining LAMOST DR4 and Gaia DR2 common red clump stars with age and proper motion, we analyze the amplitude evolution of the stellar warp independently of any assumption with a simple model. The greatest height of the warp disk increases with Galactocentric distance in different populations and it is dependent on the age: the younger stellar populations exhibit stronger warp features than the old ones, accompanied with the warp amplitude $\gamma ({\rm age})$ decreasing with age and its first derivative $\dot{\gamma} ({\rm age})$ is different from zero. The azimuth of line of nodes $\phi _w$ is stable at $-$5 degree without clear time evolution, which perfectly confirms some of previous works. All these self-consistent evidences support that our Galactic warp should most likely be a long lived, but non-steady structure and not a transient one, which is supporting the warp is originated from gas infall onto the disk or other hypotheses that suppose that the warp mainly affects to the gas, and consequently younger populations tracing the gas are stronger than the older ones. In other words, the Galactic warp is induced by the non-gravitational interaction over the disk models.

We report the discovery of two ultra-faint stellar systems found in early data from the DECam Local Volume Exploration survey (DELVE). The first system, Centaurus I (DELVE J1238-4054), is identified as a resolved overdensity of old and metal-poor stars with a heliocentric distance of ${\rm D}_{\odot} = 116.3_{-0.6}^{+0.6}$ kpc, a half-light radius of $r_h = 2.3_{-0.3}^{+0.4}$ arcmin, an age of $\tau > 12.85$ Gyr, a metallicity of $Z = 0.0002_{-0.0002}^{+0.0001}$, and an absolute magnitude of $M_V = -5.55_{-0.11}^{+0.11}$ mag. This characterization is consistent with the population of ultra-faint satellites, and confirmation of this system would make Centaurus I one of the brightest recently discovered ultra-faint dwarf galaxies. Centaurus I is detected in Gaia DR2 with a clear and distinct proper motion signal, confirming that it is a real association of stars distinct from the Milky Way foreground; this is further supported by the clustering of blue horizontal branch stars near the centroid of the system. The second system, DELVE 1 (DELVE J1630-0058), is identified as a resolved overdensity of stars with a heliocentric distance of ${\rm D}_{\odot} = 19.0_{-0.6}^{+0.5} kpc$, a half-light radius of $r_h = 0.97_{-0.17}^{+0.24}$ arcmin, an age of $\tau = 12.5_{-0.7}^{+1.0}$ Gyr, a metallicity of $Z = 0.0005_{-0.0001}^{+0.0002}$, and an absolute magnitude of $M_V = -0.2_{-0.6}^{+0.8}$ mag, consistent with the known population of faint halo star clusters. Given the low number of probable member stars at magnitudes accessible with Gaia DR2, a proper motion signal for DELVE 1 is only marginally detected. We compare the spatial position and proper motion of both Centaurus I and DELVE 1 with simulations of the accreted satellite population of the Large Magellanic Cloud (LMC) and find that neither is likely to be associated with the LMC.

A microphysics-agnostic meta-model of rotational glitches in rotation-powered pulsars is developed, wherein the globally averaged internal stress accumulates as a Brownian process between glitches, and a glitch is triggered once a critical threshold is surmounted. Precise, falsifiable predictions are made regarding long-term event statistics in individual pulsars. For example, the Spearman cross-correlation coefficient between the size of a glitch and the waiting time until the next glitch should exceed 0.25 in all pulsars. Among the six pulsars with the most recorded glitches, PSR J0537$-$6910 and PSR J0835$-$4510 are consistent with the predictions of the meta-model, while PSR J1740$-$3015 and PSR J0631$+$1036 are not. PSR J0534$+$2200 and PSR J1341$-$6220 are only consistent with the meta-model, if there exists an undetected population of small glitches with small waiting times, which we do not resolve. The results are compared with a state-dependent Poisson process, another microphysics-agnostic meta-model in the literature. The results are also applied briefly to recent pulse-to-pulse observations of PSRJ0835$-$4510, which appear to reveal evidence for a negative fluctuation in rotation frequency just prior to the 2016 glitch.

We perform analysis of the three-dimensional kinematics of Milky Way disk stars in mono-age populations. We focus on stars between Galactocentric distances of $R=6$ and 14 \,kpc, selected from the combined LAMOST DR4 red clump giant stars and Gaia DR2 proper motion catalogue. We confirm the 3D asymmetrical motions of recent works and provide time tagging of the Galactic outer disk asymmetrical motions near the anticenter direction out to Galactocentric distances of 14\,kpc. Radial Galactocentric motions reach values up to 10 km s$^{-1}$, depending on the age of the population, and present a north-south asymmetry in the region corresponding to density and velocity substructures that were sensitive to the perturbations in the early 6 \,Gyr. After that time, the disk stars in this asymmetrical structure have become kinematically hotter, and are thus not sensitive to perturbations, and we find the structure is a relatively younger population. With quantitative analysis, we find stars both above and below the plane at $R\gtrsim 9$ kpc that exhibit bending mode motions of which the sensitive duration is around 8 \,Gyr. We speculate that the in-plane asymmetries might not be mainly caused by a fast rotating bar, intrinsically elliptical outer disk, secular expansion of the disk, or streams. Spiral arm dynamics, out-of-equilibrium models, minor mergers or others are important contributors. Vertical motions might be dominated by bending and breathing modes induced by complicated inner or external perturbers. It is likely that many of these mechanisms are coupled together.

We report the discovery of a $>1^\circ$ ($\sim50$ kpc) long stellar tidal stream emanating from the dwarf galaxy DDO 44, a likely satellite of Local Volume galaxy NGC 2403 located $\sim70$ kpc in projection from its companion. NGC 2403 is a roughly Large Magellanic Cloud stellar-mass galaxy 3 Mpc away, residing at the outer limits of the M 81 group. We are mapping a large region around NGC 2403 as part of our MADCASH (Magellanic Analogs' Dwarf Companions and Stellar Halos) survey, reaching point source depths (90% completeness) of ($g, i$) = (26.5, 26.2). Density maps of old, metal-poor RGB stars reveal tidal streams extending on two sides of DDO 44, with the streams directed toward NGC 2403. We estimate total luminosities of the original DDO 44 system (dwarf and streams combined) to be $M_{i, \rm{tot}} = -13.4$ and $M_{g, \rm{tot}} = -12.6$, with $\sim25-30\%$ of the luminosity in the streams. Analogs of $\sim$LMC-mass hosts with massive tidally disrupting satellites are rare in the Illustris simulations, especially at large separations such as that of DDO 44. The few analogs that are present in the models suggest that even low-mass hosts can efficiently quench their massive satellites.

We characterize the kinematic and chemical properties of $\sim$3,000 Sagittarius (Sgr) stream stars, including K-giants, M-giants, and BHBs, select from SEGUE-2, LAMOST, and SDSS separately in Integrals-of-Motion space. The orbit of Sgr stream is quite clear from the velocity vector in $X$-$Z$ plane. Stars traced by K-giants and M-giants present the apogalacticon of trailing steam is $\sim$ 100 kpc. The metallicity distributions of Sgr K-, M-giants, and BHBs present that the M-giants are on average the most metal-rich population, followed by K-giants and BHBs. All of the K-, M-giants, and BHBs indicate that the trailing arm is on average more metal-rich than leading arm, and the K-giants show that the Sgr debris is the most metal-poor part. The $\alpha$-abundance of Sgr stars exhibits a similar trend with the Galactic halo stars at lower metallicity ([Fe/H] $<\sim$ $-$1.0 dex), and then evolve down to lower [$\alpha$/Fe] than disk stars at higher metallicity, which is close to the evolution pattern of $\alpha$-element of Milky Way dwarf galaxies. We find $V_Y$ and metallicity of K-giants have gradients along the direction of line-of-sight from the Galactic center in $X$-$Z$ plane, and the K-giants show that $V_Y$ increases with metallicity at [Fe/H] $>\sim-$1.5 dex. After dividing the Sgr stream into bright and faint stream according to their locations in equatorial coordinate, the K-giants and BHBs show that the bright and faint stream present different $V_Y$ and metallicities, the bright stream is on average higher in $V_Y$ and metallicity than the faint stream.

We investigate the three-dimensional asymmetrical kinematics and present time stamps of the Milky Way disk between Galactocentric distances of $R=12$ and 15 \,kpc, using red clump stars selected from the LAMOST Galactic survey, also with proper motion measurements provided by the Gaia DR2. We discover velocity substructure above the Galactic plane corresponding to a density dip found recently ("south-middle opposite" density structure[R $\sim$ 12-15 \,kpc, Z $\sim$ 1.5 \,kpc] discovered in \citet{Wang2018b, Wang2018c}) in the radial and azimuthal velocity. For the vertical velocity, we detect clear vertical bulk motions or bending mode motions, which has no clear north-south asymmetry corresponding to the in-plane asymmetrical features. In the subsample of stars with different ages, we find that there is little temporal evolution of the in-plane asymmetry from 0$-$14 \,Gyr, which means the structure is sensitive to the perturbations in almost cosmic time possibly. We propose that the possible scenario of this asymmetric velocity structure is caused by the mechanisms generated in-plane, rather than vertical perturbations.

In this work, we present new catalogs of M giant and M dwarf stars from the LAMOST DR5. In total, 39,796 M giants and 501,152 M dwarfs are identified from the classification pipeline. The template-fitting results contain M giants with 7 temperature subtypes from M0 to M6, M dwarfs with 18 temperature subtypes from K7.0 to M8.5 and 12 metallicity subclasses from dMr to usdMp. We cross-matched our M-type catalog with the 2MASS and WISE catalog to obtain infrared magnitude and colors. Adopting the distances derived from the parallaxes in \gaia{} DR2, the M_G vs. (G_bp-G_rp)_0 diagram shows that there are also early-type stars and white dwarf-M dwarf binaries included in our M type stars sample, with a contamination rate of about 4.6% for M giants and 0.48% for M dwarfs. We found that CaH spectral indices are an efficient selection criteria for carbon stars. A total of 289 carbon stars were identified from the M giants sample, and further confirmed by LAMOST spectra.

Stellar halos offer fossil evidence for hierarchical structure formation. Since halo assembly is predicted to be scale-free, stellar halos around low-mass galaxies constrain properties such as star formation in the accreted subhalos and the formation of dwarf galaxies. However, few observational searches for stellar halos in dwarfs exist. Here we present gi photometry of resolved stars in isolated Local Group dwarf irregular galaxy IC 1613 ($M_{\star} \sim 10^8 M_{\odot})$. These Subaru/Hyper Suprime-Cam observations are the widest and deepest of IC 1613 to date. We measure surface density profiles of young main-sequence, intermediate to old red giant branch, and ancient horizontal branch stars outside of 12' ($\sim 2.6$ kpc; 2.5 half-light radii) from the IC 1613 center. All of the populations extend to ~24' ($\sim 5.2$ kpc; 5 half-light radii), with the older populations best fit by a broken exponential in these outer regions. Comparison with earlier studies sensitive to IC 1613's inner regions shows that the density of old stellar populations steepens substantially with distance from the center; we trace the $g$-band effective surface brightness to an extremely faint limit of $\sim 33.7$ mag arcsec$^{-2}$. Conversely, the distribution of younger stars follows a single, shallow exponential profile in the outer regions, demonstrating different formation channels for the younger and older components of IC 1613. The outermost, intermediate-age and old stars have properties consistent with those expected for accreted stellar halos, though future observational and theoretical work is needed to definitively distinguish this scenario from other possibilities.

Among the five pulsars with the most recorded rotational glitches, only PSR J0534$+$2200 is found to have an autocorrelation between consecutive glitch sizes which differs significantly from zero (Spearman correlation coefficient $\rho=-0.46$, p-value $=0.046$). No statistically compelling autocorrelations between consecutive waiting times are found. The autocorrelation observations are interpreted within the framework of a predictive meta-model describing stress-release in terms of a state-dependent Poisson process. Specific combinations of size and waiting time autocorrelations are identified, alongside combinations of cross-correlations and size and waiting time distributions, that are allowed or excluded within the meta-model. For example, future observations of any "quasiperiodic" glitching pulsar, such as PSR J0537$-$6910, should not reveal a positive waiting time autocorrelation. The implications for microphysical models of the stress-release process driving pulsar glitches are discussed briefly.

The ultra-faint dwarf galaxy Leo V has shown both photometric overdensities and kinematic members at large radii, along with a tentative kinematic gradient, suggesting that it may have undergone a close encounter with the Milky Way. We investigate these signs of disruption through a combination of i) high-precision photometry obtained with the Hubble Space Telescope (HST), ii) two epochs of stellar spectra obtained with the Hectochelle Spectrograph on the MMT, and iii) measurements from the Gaia mission. Using the HST data, we examine one of the reported stream-like overdensities at large radii, and conclude that it is not a true stellar stream, but instead a clump of foreground stars and background galaxies. Our spectroscopic analysis shows that one known member star is likely a binary, and challenges the membership status of three others, including two distant candidates that had formerly provided evidence for overall stellar mass loss. We also find evidence that the proposed kinematic gradient across Leo V might be due to small number statistics. We update the systemic proper motion of Leo V, finding $(\mu_\alpha \cos\delta, \mu_\delta)= (0.009\pm0.560$, $-0.777\pm0.314)$ mas yr$^{-1}$, which is consistent with its reported orbit that did not put Leo V at risk of being disturbed by the Milky Way. These findings remove most of the observational clues that suggested Leo V was disrupting, however, we also find new plausible member stars, two of which are located >5 half-light radii from the main body. These stars require further investigation. Therefore, the nature of Leo V still remains an open question.

A census of the satellite population around dwarf galaxy primary hosts in environments outside the Local Group is essential to understanding $\Lambda$CDM galaxy formation and evolution on the smallest scales. We present deep optical Hubble Space Telescope imaging of the gas-rich, faint dwarf galaxy Antlia B ($M_V = -9.4$) -- a likely satellite of NGC 3109 ($D = 1.3$ Mpc) -- discovered as part of our ongoing survey of primary host galaxies similar to the Magellanic Clouds. We derive a new tip of the red giant branch (TRGB) distance of $D = 1.35 \pm 0.06$ Mpc ($m-M = 25.65 \pm 0.10$), consistent with membership in the nearby NGC 3109 dwarf association. The color-magnitude diagram shows both a prominent old, metal-poor stellar component and confirms a small population of young, blue stars with ages $\lesssim 1$ Gyr. We use the color-magnitude diagram fitting algorithm MATCH to derive the star formation history and find that it is consistent with the typical dwarf irregular or transitional dwarf galaxy (dTrans) in the Local Group. Antlia B shows relatively constant stellar mass growth for the first $\sim 10-11$ Gyr and almost no growth in the last $\sim 2-3$ Gyr. Despite being gas-rich, Antlia B shows no evidence of active star formation (i.e., no H$\alpha$ emission) and should therefore be classified as a dTrans dwarf. Both Antlia B and the Antlia dwarf (dTrans) are likely satellites of NGC 3109 suggesting that the cessation of ongoing star formation in these galaxies may be environmentally driven. Future work studying the gas kinematics and distribution in Antlia B will explore this scenario in greater detail. Our work highlights the fact that detailed studies of nearby dwarf galaxies in a variety of environments may continue to shed light on the processes that drive the star formation history and evolution of dwarf galaxies more generally.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

Studying our Galaxy, the Milky Way (MW), gives us a close-up view of the interplay between cosmology, dark matter, and galaxy formation. In the next decade our understanding of the MW's dynamics, stellar populations, and structure will undergo a revolution thanks to planned and proposed astrometric, spectroscopic and photometric surveys, building on recent advances by the Gaia astrometric survey. Together, these new efforts will measure three-dimensional positions and velocities and numerous chemical abundances for stars to the MW's edge and well into the Local Group, leading to a complete multidimensional view of our Galaxy. Studies of the multidimensional Milky Way beyond the Gaia frontier---from the edge of the Galactic disk to the edge of our Galaxy's dark matter halo---will unlock new scientific advances across astrophysics, from constraints on dark matter to insights into galaxy formation.

We selected blue horizontal branch (BHB) stars within the expected distance range and sky position of the Hermus Stream from Data Release 10 of the Sloan Digital Sky Survey. We identify a moving group of $19$ BHB stars that are concentrated within two degrees of the Hermus Stream, between $10$ and $14$ kpc from the Sun. The concentration in velocity is inconsistent with a Gaussian distribution with 98% confidence (2.33 sigma).The stars in the moving group have line-of-sight velocities of v$_{\rm gsr} \sim 50$ km s$^{-1}$, a velocity dispersion of $\sigma_v\lesssim11$ km s$^{-1}$, a line-of-sight depth of $\sim 1$ kpc, and a metallicity of [Fe/H] $=-2.1 \pm 0.4$. The best-fit orbit has a perigalacticon of $\sim4$ kpc, apogalacticon of $\sim 17$ kpc, orbital period of $\sim247$ Myr, eccentricity $e = 0.62$, and inclination $i \sim75^\circ$ from $b=90^\circ$. The BHB stars in the stream are estimated to be $12$ Gyr old. An N-body simulation of a mass-follows-light ultrafaint dwarf galaxy with mass $10^6M_{\odot}$ and radius $40$ pc is consistent with the observed properties. The properties of the identified moving group of 19 BHB stars are close enough to those of the Hermus Stream (which is traced predominantly in turnoff stars) that we find it likely that they are associated. If that is the case, then our orbit fit would imply that there is no relationship between the Hermus and Phoenix streams, as previously proposed.

We use LAMOST DR4 M giants combined with Gaia DR2 proper motions and ALLWISE photometry to obtain an extremely pure sample of Sagittarius (Sgr) stream stars. Using TiO5 and CaH spectral indices as an indicator, we selected out a large sample of M giant stars from M dwarf stars in LAMOST DR4 spectra. Considering the position, distance, proper motion and the angular momentum distribution, we obtained 164 pure Sgr stream stars. We find that the trailing arm has higher energy than the leading arm in same angular momentum. The trailing arm we detected extends to a heliocentric distance of $\sim 130$ kpc at $\tilde\Lambda_{\odot}\sim 170^{\circ}$, which is consistent with the feature found in RR Lyrae in \citet{sgr2017}. Both of these detections of Sgr, in M giants and in RR Lyrae, imply that the Sgr stream may contain multiple stellar populations with a broad metallicity range.

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

Glitching pulsars fall broadly into two statistical classes: those with Poisson-like waiting times and power-law sizes, and those with unimodal waiting times and sizes. Previous glitch modeling based on a state-dependent Poisson process readily generates Poisson-like behaviour but struggles to produce unimodal waiting times or sizes. Here it is shown that, when some of the inputs to the model are modified, both classes of statistical behaviour can be reproduced by varying a single control parameter related to the spin-down rate. The implications for past and future glitch observations and the underlying microphysical mechanism are explored briefly.

It is an open question whether glitch activity in individual pulsars varies on decadal time-scales. The Crab pulsar has experienced 23 spin-up glitches in the last 36 years, interrupting an otherwise monotonic deceleration. A homogeneous Poisson process, i.e. a process with constant rate, is not sufficient to describe the time-ordered distribution of glitch epochs in the Crab pulsar. There are signs of clustering at the $2\sigma$ level when testing with Ripley's $K$ function. Two alternative, inhomogeneous models with one and two step-wise rate changes are found to have higher relative evidence (Bayes factors of 1.74 and 2.86 respectively) than the homogeneous Poisson process. The distinction between clustering, where events are correlated, and rate variation is discussed. The implications for glitch microphysics, in particular trigger mechanisms based on avalanche processes, are briefly discussed.

We present proper motion measurements of Bo\"{o}tes III, an enigmatic stellar satellite of the Milky Way, utilizing data from the second data release of the Gaia mission. By selecting 15 radial velocity confirmed members of Bo\"{o}tes III, along with a likely RR Lyrae member in the vicinity, we measure an error weighted mean proper motion of $(\mu_{\alpha} \cos{\delta}, \mu_\delta) = (-1.14, -0.98)\pm(0.18, 0.20)$ mas/yr. We select and present further stars that may be Bo\"{o}tes III members based on their combined proper motion and position in the color magnitude diagram. We caution against assigning membership to stars that are not confirmed spectroscopically, as we demonstrate that there are contaminating stars from the disrupting globular cluster NGC 5466 in the vicinity of the main body of Bo\"{o}tes III, but we note that our results are consistent with previous Bo\"{o}tes III proper motion estimates that did not include spectroscopic members. Based on the measured proper motion and other known properties of Bo\"{o}tes III, we derive its Galactocentric velocity and compute its orbit given canonical Milky Way potentials with halo masses of both 0.8$\times$10$^{12}$ M$_{\odot}$ and 1.6$\times$10$^{12}$ M$_{\odot}$. These orbits robustly show that Bo\"{o}tes III passed within $\sim$12 kpc of the Galactic center on an eccentric orbit roughly 140 Myr ago. Additionally, the proper motion of Bo\"{o}tes III is in excellent agreement with predictions for the retrograde motion of the coincident Styx stellar stream. Given this, along with the small pericenter and metallicity spread of Bo\"{o}tes III itself, we suggest that it is a disrupting dwarf galaxy giving rise to the Styx stellar stream.

Triangulum (M33) is the most massive satellite galaxy of Andromeda (M31), with a stellar mass of about $3\times10^9\; M_{\odot}$. Based on abundance matching techniques, M33's total mass at infall is estimated to be of order $10^{11}\; M_{\odot}$. $\Lambda$CDM simulations predict that M33-mass halos host several of their own low mass satellite companions, yet only one candidate M33 satellite has been detected in deep photometric surveys to date. This `satellites of satellites' hierarchy has recently been explored in the context of the dwarf galaxies discovered around the Milky Way's Magellanic Clouds in the Dark Energy Survey. Here we discuss the number of satellite galaxies predicted to reside within the virial radius ($\sim$160 kpc) of M33 based on $\Lambda$CDM simulations. We also calculate the expected number of satellite detections in N fields of data using various ground--based imagers. Finally, we discuss how these satellite population predictions change as a function of M33's past orbital history. If M33 is on its first infall into M31's halo, its proposed satellites are expected to remain bound to M33 today. However, if M33 experienced a recent tidal interaction with M31, the number of surviving satellites depends strongly on the distance achieved at pericenter due to the effects of tidal stripping. We conclude that a survey extending to $\sim$100 kpc around M33 would be sufficient to constrain its orbital history and a majority of its satellite population. In the era of WFIRST, surveys of this size will be considered small observing programs.

We present deep Magellan/Megacam stellar photometry of four recently discovered faint Milky Way satellites: Sagittarius II (Sgr II), Reticulum II (Ret II), Phoenix II (Phe II), and Tucana III (Tuc III). Our photometry reaches ~2-3 magnitudes deeper than the discovery data, allowing us to revisit the properties of these new objects (e.g., distance, structural properties, luminosity measurements, and signs of tidal disturbance). The satellite color-magnitude diagrams show that they are all old (~13.5 Gyr) and metal-poor ([Fe/H]$\lesssim-2.2$). Sgr II is particularly interesting as it sits in an intermediate position between the loci of dwarf galaxies and globular clusters in the size-luminosity plane. The ensemble of its structural parameters is more consistent with a globular cluster classification, indicating that Sgr II is the most extended globular cluster in its luminosity range. The other three satellites land directly on the locus defined by Milky Way ultra-faint dwarf galaxies of similar luminosity. Ret II is the most elongated nearby dwarf galaxy currently known for its luminosity range. Our structural parameters for Phe II and Tuc III suggest that they are both dwarf galaxies. Tuc III is known to be associated with a stellar stream, which is clearly visible in our matched-filter stellar density map. The other satellites do not show any clear evidence of tidal stripping in the form of extensions or distortions. Finally, we also use archival HI data to place limits on the gas content of each object.

We analyze chemical abundances of stars in the Sagittarius (Sgr) tidal stream using high-resolution Gemini+GRACES spectra of 42 members of the highest surface brightness portions of both the trailing and leading arms. Targets were chosen using a 2MASS+WISE color-color selection, combined with LAMOST radial velocities. In this study, we analyze [Fe/H] and alpha-elements produced by both hydrostatic (O, Mg) and explosive (Si, Ca, Ti) nucleosynthetic processes. The average [Fe/H] for our Sgr stream stars is lower than that for stars in the Sgr core, and stars in the trailing and leading arms show systematic differences in [Fe/H]. Both hydrostatic and explosive elements are depleted relative to Milky Way (MW) disk and halo stars, with a larger gap between the MW trend and Sgr stars for the hydrostatic elements. Chemical abundances of Sgr stream stars show similar patterns to those measured in the core of the Sgr dSph. We explore the ratio of hydrostatic to explosive alpha-elements [$\alpha_{\rm h/ex}$] (which we refer to as the "HEx ratio"). Our observed HEx ratio trends for Sgr debris are deficient relative to MW stars. Via simple chemical evolution modeling, we show that these HEx ratio patterns are consistent with a Sgr IMF that lacks the most massive stars. This study provides a link between the chemical properties in the intact Sgr core and the significant portion of the Sgr system's luminosity that is estimated to currently reside in the streams.

We present a three dimensional velocity analysis of Milky Way disk kinematics using LAMOST K giant stars and the GPS1 proper motion catalogue. We find that Galactic disk stars near the anticenter direction (in the range of Galactocentric distance between $R=8$ and $13$ kpc and vertical position between $Z=-2$ and $2$ kpc) exhibit asymmetrical motions in the Galactocentric radial, azimuthal, and vertical components. Radial motions are not zero, thus departing from circularity in the orbits; they increase outwards within $R\lesssim 12$ kpc, show some oscillation in the northern ($0 < Z < 2$ kpc) stars, and have north-south asymmetry in the region corresponding to a well-known nearby northern structure in the velocity field. There is a clear vertical gradient in azimuthal velocity, and also an asymmetry that shifts from a larger azimuthal velocity above the plane near the solar radius to faster rotation below the plane at radii of 11-12 kpc. Stars both above and below the plane at $R\gtrsim 9$ kpc exhibit net upward vertical motions. We discuss some possible mechanisms that might create the asymmetrical motions, such as external perturbations due to dwarf galaxy minor mergers or dark matter sub-halos, warp dynamics, internal processes due to spiral arms or the Galactic bar, and (most likely) a combination of some or all of these components.

We report the discovery of distant RR Lyrae stars, including the most distant known in the Milky Way, using data taken in the $g-$band with the Dark Energy Camera as part of the High cadence Transient Survey (HiTS; 2014 campaign). We detect a total of 173 RR Lyrae stars over a ~120 deg^2 area, including both known RR Lyrae and new detections. The heliocentric distances d_H of the full sample range from 9 to >200 kpc, with 18 of them beyond 90 kpc. We identify three sub-groups of RR Lyrae as members of known systems: the Sextans dwarf spheroidal galaxy, for which we report 46 new discoveries, and the ultra-faint dwarf galaxies Leo IV and Leo V. Following an MCMC methodology, we fit spherical and ellipsoidal profiles of the form rho(R) ~ R^n to the radial density distribution of RR Lyrae in the Galactic halo. The best fit corresponds to the spherical case, for which we obtain a simple power-law index of n=-4.17^{+0.18}_{-0.20}, consistent with recent studies made with samples covering shorter distances. The pulsational properties of the outermost RR Lyrae in the sample (d_H>90 kpc) differ from the ones in the halo population at closer distances. The distribution of the stars in a Period-Amplitude diagram suggest they belong to Oosterhoff-intermediate or Oosterhoff II groups, similar to what is found in the ultra-faint dwarf satellites around the Milky Way. The new distant stars discovered represent an important addition to the few existing tracers of the Milky Way potential in the outer halo.

We investigate the hypothesized tidal disruption of the Hercules ultra-faint dwarf galaxy (UFD). Previous tidal disruption studies of the Hercules UFD have been hindered by the high degree of foreground contamination in the direction of the dwarf. We bypass this issue by using RR Lyrae stars, which are standard candles with a very low field-volume density at the distance of Hercules. We use wide-field imaging from the Dark Energy Camera on CTIO to identify candidate RR Lyrae stars, supplemented with observations taken in coordination with the Beijing-Arizona Sky Survey on the Bok Telescope. Combining color, magnitude, and light-curve information, we identify three new RR Lyrae stars associated with Hercules. All three of these new RR Lyrae stars lie outside its published tidal radius. When considered with the nine RR Lyrae stars already known within the tidal radius, these results suggest that a substantial fraction of Hercules' stellar content has been stripped. With this degree of tidal disruption, Hercules is an interesting case between a visibly disrupted dwarf (such as the Sagittarius dwarf spheroidal galaxy) and one in dynamic equilibrium. The degree of disruption also shows that we must be more careful with the ways we determine object membership when estimating dwarf masses in the future. One of the three discovered RR Lyrae stars sits along the minor axis of Hercules, but over two tidal radii away. This type of debris is consistent with recent models that suggest Hercules' orbit is aligned with its minor axis.

We present deep, wide-field Subaru Hyper Suprime-Cam photometry of two recently discovered satellites of the Milky Way (MW): Columba I and Triangulum II. The color magnitude diagrams of both objects point to exclusively old and metal-poor stellar populations. We re-derive structural parameters and luminosities of these satellites, and find $M_{\rm V, Col~I} = -4.2\pm0.2$ for Col I and $M_{\rm V, Tri~II} = -1.2\pm0.4$ for Tri II, with corresponding half-light radii of $r_{\rm h, Col~I} = 117\pm17$ pc and $r_{\rm h, Tri~II} = 21\pm4$ pc. The properties of both systems are consistent with observed scaling relations for MW dwarf galaxies. Based on archival data, we derive upper limits on the neutral gas content of these dwarfs, and find that they lack HI, as do the majority of observed satellites within the MW virial radius. Neither satellite shows evidence of tidal stripping in the form of extensions or distortions in matched-filter stellar density maps or surface density profiles. However, the smaller Tri II system is relatively metal-rich for its luminosity (compared to other MW satellites), possibly because it has been tidally stripped. Through a suite of orbit simulations, we show that Tri II is approaching pericenter of its eccentric orbit, a stage at which tidal debris is unlikely to be seen. In addition, we find that Tri II may be on its first infall into the MW, which helps explain its unique properties among MW dwarfs. Further evidence that Tri II is likely an ultra-faint dwarf comes from its stellar mass function, which is similar to those of other MW dwarfs.

We present UV, optical, and NIR photometry of the first electromagnetic counterpart to a gravitational wave source from Advanced LIGO/Virgo, the binary neutron star merger GW170817. Our data set extends from the discovery of the optical counterpart at $0.47$ days to $18.5$ days post-merger, and includes observations with the Dark Energy Camera (DECam), Gemini-South/FLAMINGOS-2 (GS/F2), and the {\it Hubble Space Telescope} ({\it HST}). The spectral energy distribution (SED) inferred from this photometry at $0.6$ days is well described by a blackbody model with $T\approx 8300$ K, a radius of $R\approx 4.5\times 10^{14}$ cm (corresponding to an expansion velocity of $v\approx 0.3c$), and a bolometric luminosity of $L_{\rm bol}\approx 5\times10^{41}$ erg s$^{-1}$. At $1.5$ days we find a multi-component SED across the optical and NIR, and subsequently we observe rapid fading in the UV and blue optical bands and significant reddening of the optical/NIR colors. Modeling the entire data set we find that models with heating from radioactive decay of $^{56}$Ni, or those with only a single component of opacity from $r$-process elements, fail to capture the rapid optical decline and red optical/NIR colors. Instead, models with two components consistent with lanthanide-poor and lanthanide-rich ejecta provide a good fit to the data, the resulting "blue" component has $M_\mathrm{ej}^\mathrm{blue}\approx 0.01$ M$_\odot$ and $v_\mathrm{ej}^\mathrm{blue}\approx 0.3$c, and the "red" component has $M_\mathrm{ej}^\mathrm{red}\approx 0.04$ M$_\odot$ and $v_\mathrm{ej}^\mathrm{red}\approx 0.1$c. These ejecta masses are broadly consistent with the estimated $r$-process production rate required to explain the Milky Way $r$-process abundances, providing the first evidence that BNS mergers can be a dominant site of $r$-process enrichment.

We present the Dark Energy Camera (DECam) discovery of the optical counterpart of the first binary neutron star merger detected through gravitational wave emission, GW170817. Our observations commenced 10.5 hours post-merger, as soon as the localization region became accessible from Chile. We imaged 70 deg$^2$ in the $i$ and $z$ bands, covering 93\% of the initial integrated localization probability, to a depth necessary to identify likely optical counterparts (e.g., a kilonova). At 11.4 hours post-merger we detected a bright optical transient located $10.6''$ from the nucleus of NGC\,4993 at redshift $z=0.0098$, consistent (for $H_0 = 70$\, km s$^{-1}$ Mpc$^{-1}$) with the distance of $40 \pm 8$\, Mpc reported by the LIGO Scientific Collaboration and the Virgo Collaboration (LVC). At detection the transient had magnitudes $i\approx 17.30$ and $z\approx 17.45$, and thus an absolute magnitude of $M_i = -15.7$, in the luminosity range expected for a kilonova. We identified 1,500 potential transient candidates. Applying simple selection criteria aimed at rejecting background events such as supernovae, we find the transient associated with NGC\,4993 as the only remaining plausible counterpart, and reject chance coincidence at the 99.5\% confidence level. We therefore conclude that the optical counterpart we have identified near NGC\,4993 is associated with GW170817. This discovery ushers in the era of multi-messenger astronomy with gravitational waves, and demonstrates the power of DECam to identify the optical counterparts of gravitational-wave sources.

The radial number density and flattening of the Milky Way's stellar halo is measured with $\mathrm{5351}$ metal-poor ([Fe/H]$<-1$) K giants from LAMOST DR3, using a nonparametric method which is model independent and largely avoids the influence of halo substucture. The number density profile is well described by a single power law with index $5.03^{+0.64}_{-0.64}$, and flattening that varies with radius. The stellar halo traced by LAMOST K giants is more flattened at smaller radii, and becomes nearly spherical at larger radii. The flattening, $q$, is about 0.64, 0.8, 0.96 at $r=15$, 20 and 30 kpc (where $r=\sqrt{R^2+\left[Z/q\left(r\right)\right]^2}$), respectively. Moreover, the leading arm of the Sagittarius dwarf galaxy tidal stream in the north, and the trailing arm in the south, are significant in the residual map of density distribution. In addition, an unknown overdensity is identified in the residual map at (R,Z)=(30,15) kpc.

Recent discovery of many dwarf satellite galaxies in the direction of the Small and Large Magellanic Clouds (SMC and LMC) provokes questions of their origins, and what they can reveal about galaxy evolution theory. Here, we predict the satellite stellar mass function of Magellanic Cloud-mass host galaxies using abundance matching and reionization models applied to the \textit{Caterpillar} simulations. Specifically focusing on the volume within $50$~kpc of the LMC, we predict a mean of 4-8 satellites with stellar mass $M_* > 10^4\, \mathrm{M_\odot}$, and 3-4 satellites with $80 < M_* \leq 3000\, \mathrm{M_\odot}$. Surprisingly, all $12$ currently known satellite candidates have stellar masses of $80 < M_* \leq 3000\, \mathrm{M_\odot}$. Reconciling the dearth of large satellites and profusion of small satellites is challenging and may require a combination of a major modification of the \mstarmhalo relationship (steep, but with an abrupt flattening at $10^3\, \mathrm{M_\odot}$), late reionization for the Local Group ($z_{\rm{reion}} \lesssim 9$ preferred), and/or strong tidal stripping. We can more robustly predict that $\sim 53\%$ of satellites within this volume were accreted together with the LMC and SMC, and $\sim 47\%$ were only ever Milky Way satellites. Observing satellites of isolated LMC-sized field galaxies is essential to placing the LMC in context, and to better constrain the \mstarmhalo relationship. Modeling known LMC-sized galaxies within $8$ Mpc, we predict 1-6 (2-12) satellites with $M_* > 10^5\, \mathrm{M_\odot}$ ($M_* > 10^4\, \mathrm{M_\odot}$) within the virial volume of each, and 1-3 (1-7) within a single $1.5^{\circ}$ diameter field of view, making their discovery likely.

We confirm, quantify, and provide a table of the coherent velocity substructure of the Milky Way disk within 2 kpc of the Sun towards the Galactic anticenter, with 0.2 kpc resolution. We use the radial velocities of ~340,000 F-type stars obtained with the Guoshoujing Telescope (also known as the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST), and proper motions derived from the PPMXL catalog. The PPMXL proper motions have been corrected to remove systematic errors by subtracting the average proper motions of galaxies and QSOs that have been confirmed in the LAMOST spectroscopic survey, and that are within 2.5 degrees of the star's position. We provide the resulting table of systematic offsets derived from the PPMXL proper motion measurements of extragalactic objects identified in the LAMOST spectroscopic survey. Using the corrected phase- space stellar sample, we find statistically significant deviations in the bulk disk velocity of 20 km/s or more in the three dimensional velocities of Galactic disk stars. The bulk velocity varies significantly over length scales of half a kpc or less. The rotation velocity of the disk increases by 20 km/s from the Sun's position to 1.5 kpc outside the solar circle. Disk stars in the second quadrant, within 1 kpc of the Sun, are moving radially towards the Galactic center and vertically towards a point a few tenths of a kpc above the Galactic plane; looking down on the disk, the stars appear to move in a circular streaming motion with a radius of order 1 kpc.

During the analysis of RR Lyrae stars discovered in the High cadence Transient Survey (HiTS) taken with the Dark Energy Camera at the 4-m telescope at Cerro Tololo Inter-American Observatory, we found a group of three very distant, fundamental mode pulsator RR Lyrae (type ab). The location of these stars agrees with them belonging to the Leo V ultra-faint satellite galaxy, for which no variable stars have been reported to date. The heliocentric distance derived for Leo V based on these stars is 173 +/- 5 kpc. The pulsational properties (amplitudes and periods) of these stars locate them within the locus of the Oosterhoff II group, similar to most other ultra-faint galaxies with known RR Lyrae stars. This serendipitous discovery shows that distant RR Lyrae stars may be used to search for unknown faint stellar systems in the outskirts of the Milky Way.

We present a statistical method to derive the stellar density profiles of the Milky Way from spectroscopic survey data, taking into account selection effects. We assume that the selection function of the spectroscopic survey is based on photometric colors and magnitudes and possibly altered during observations and data reductions. Then the underlying selection function for a line-of-sight can be well recovered by comparing the distribution of the spectroscopic stars in a color-magnitude plane with that of the photometric dataset. Subsequently, the stellar density profile along a line-of-sight can be derived from the spectroscopically measured stellar density profile multiplied by the selection function. The method is validated using Galaxia mock data with two different selection functions. We demonstrate that the derived stellar density profiles well reconstruct the true ones not only for the full targets, but also for the sub-populations selected from the full dataset. Finally, the method is applied to map the density profiles for the Galactic disk and halo, respectively, using the LAMOST RGB stars. The Galactic disk extends to about R=19 kpc, where the disk still contributes about 10% to the total stellar surface density. Beyond this radius, the disk smoothly transitions to the halo without any truncation, bending, or broken. Moreover, no over-density corresponding to the Monoceros ring is found in the Galactic anti-center direction. The disk shows moderate north-south asymmetry at radii larger than 12 kpc. On the other hand, the R-Z tomographic map directly shows that the stellar halo is substantially oblate within a Galactocentric radius of 20 kpc and gradually becomes nearly spherical beyond 30 kpc.

Single-exposure spectra in large spectral surveys are valuable for time domain studies such as stellar variability, but there is no available method to eliminate cosmic rays for single-exposure, multi-fiber spectral images. In this paper, we describe a new method to detect and remove cosmic rays in multi-fiber spectroscopic single exposures. Through the use of two-dimensional profile fitting and a noise model that considers the position-dependent errors, we successfully detect as many as 80% of the cosmic rays and correct the cosmic ray polluted pixels to an average accuracy of 97.8%. Multiple tests and comparisons with both simulated data and real LAMOST data show that the method works properly in detection rate, false detection rate, and validity of cosmic ray correction.

We present a multi-epoch Hubble Space Telescope (HST) study of stellar proper motions (PMs) for four fields along the Orphan Stream. We determine absolute PMs of several individual stars per target field using established techniques that utilize distant background galaxies to define a stationary reference frame. Five Orphan Stream stars are identified in one of the four fields based on combined color-magnitude and PM information. The average PM is consistent with the existing model of the Orphan stream by Newberg et al. In addition to the Orphan stream stars, we detect stars that likely belong to other stellar streams. To identify which stellar streams these stars belong to, we examine the 2-d bulk motion of each group of stars on the sky by subtracting the PM contribution of the solar motion (which is a function of position on the sky and distance) from the observed PMs, and comparing the vector of net motion with the spatial extent of known stellar streams. By doing this, we identify candidate stars in the Sagittarius and Lethe streams, and a newly-found stellar stream at a distance of ~17 kpc, which we tentatively name the "Parallel stream". Together with our Sagittarius stream study (Sohn et al., 2015, ApJ, 803, 56), this work demonstrates that even in the Gaia era, HST will continue to be advantageous in measuring PMs of old stellar populations on a star-by-star basis, especially for distances beyond ~10 kpc.

Stars ejected from the Galactic centre can be used to place important constraints on the Milky Way potential. Since existing hypervelocity stars are too distant to accurately determine orbits, we have conducted a search for nearby candidates using full three-dimensional velocities. Since the efficacy of such studies are often hampered by deficiencies in proper motion catalogs, we have chosen to utilize the reliable, high-precision SDSS Stripe 82 proper motion catalog. Although we do not find any candidates which have velocities in excess of the escape speed, we identify 226 stars on orbits that are consistent with Galactic centre ejection. This number is significantly larger than what we would expect for halo stars on radial orbits and cannot be explained by disk or bulge contamination. If we restrict ourselves to metal-rich stars, we find 29 candidates with [Fe/H] > -0.8 dex and 10 with [Fe/H] > -0.6 dex. Their metallicities are more consistent with what we expect for bulge ejecta, and so we believe these candidates are especially deserving of further study. We have supplemented this sample using our own radial velocities, developing an algorithm to use proper motions for optimizing candidate selection. This technique provides considerable improvement on the blind spectroscopic sample of SDSS, being able to identify candidates with an efficiency around 20 times better than a blind search.

We report the discovery of the faintest known dwarf galaxy satellite of an LMC stellar-mass host beyond the Local Group, based on deep imaging with Subaru/Hyper Suprime-Cam. MADCASH J074238+652501-dw lies $\sim$35 kpc in projection from NGC 2403, a dwarf spiral galaxy at $D$$\approx$3.2 Mpc. This new dwarf has $M_{g} = -7.4\pm0.4$ and a half-light radius of $168\pm70$ pc, at the calculated distance of $3.39\pm0.41$ Mpc. The color-magnitude diagram reveals no evidence of young stellar populations, suggesting that MADCASH J074238+652501-dw is an old, metal-poor dwarf similar to low luminosity dwarfs in the Local Group. The lack of either detected HI gas ($M_{\rm HI}/L_{V} < 0.69 M_\odot/L_\odot$, based on Green Bank Telescope observations) or $GALEX$ NUV/FUV flux enhancement is consistent with a lack of young stars. This is the first result from the MADCASH (Magellanic Analog Dwarf Companions And Stellar Halos) survey, which is conducting a census of the stellar substructure and faint satellites in the halos of Local Volume LMC analogs via resolved stellar populations. Models predict a total of $\sim$4-10 satellites at least as massive as MADCASH J074238+652501-dw around a host with the mass of NGC 2403, with 2-3 within our field of view, slightly more than the one such satellite observed in our footprint.

We report two new tidal debris nearby the Sagittarius (Sgr) tidal stream in the north Galactic cap identified from the M giant stars in LAMOST DR2 data. The M giant stars with sky area of $210^\circ<$\Lambda$<290^\circ$, distance of 10--20kpc, and [Fe/H]$<-0.75$ show clear bimodality in velocity distribution. We denote the two peaks as Vel-3+83 for the one within mean velocity of -3kms$^{-1}$ with respect to that of the well observed Sgr leading tail at the same $\Lambda$ and Vel+162+26 for the other one with mean velocity of 162kms$^{-1}$ with respect to the Sgr leading tail. Although the projected $\Lambda$--$V_{gar}$ relation of Vel-3+83 is very similar to the Sgr leading tail, the opposite trend in $\Lambda$--distance relation against the Sgr leading tail suggests Vel-3+83 has a different 3D direction of motion with any branch of the simulated Sgr tidal stream from Law & Majewski. Therefore, we propose it to be a new tidal debris not related to the Sgr stream. Similarly, the another substructure Vel+162+26, which is the same one as the NGC group discovered by Chou et al., also moves toward a different direction with the Sgr stream, implying that it may have different origin with the Sgr tidal stream.

Using a spectroscopically confirmed sample of M-giants, M-dwarfs and quasars from the LAMOST survey, we assess how well WISE $\&$ 2MASS color-cuts can be used to select M-giant stars. The WISE bands are very efficient at separating M-giants from M-dwarfs and we present a simple classification that can produce a clean and relatively complete sample of M-giants. We derive a new photometric relation to estimate the metallicity for M-giants, calibrated using data from the APOGEE survey. We find a strong correlation between the $(W1-W2)$ color and $\rm [M/H]$, where almost all of the scatter is due to photometric uncertainties. We show that previous photometric distance relations, which are mostly based on stellar models, may be biased and devise a new empirical distance relation, investigating trends with metallicity and star formation history. Given these relations, we investigate the properties of M-giants in the Sagittarius stream. The offset in the orbital plane between the leading and trailing tails is reproduced and, by identifying distant M-giants in the direction of the Galactic anti-center, we confirm that the previously detected debris in the outer halo is the apocenter of the trailing tail. We also find tentative evidence supporting an existing overdensity near the leading tail in the Northern Galactic hemisphere, possibly an extension to the trailing tail (so-called Branch C). We have measured the metallicity distribution along the stream, finding a clear metallicity offset between the leading and trailing tails, in agreement with models for the stream formation. We include an online table of M-giants to facilitate further studies.

Recent years have seen the discovery of an ever growing number of stellar debris streams and clouds. These structures are typically detected as extended and often curvilinear overdensities of metal-poor stars that stand out from the foreground disk population. The streams typically stretch tens of degrees or more across the sky, even encircling the Galaxy, and range in heliocentric distance from 3 to 100 kpc. This chapter summarizes the techniques used for finding such streams and provides tables giving positions, distances, velocities, and metallicities, where available, for all major streams and clouds that have been detected as of January 2015. Sky maps of the streams are also provided. Properties of individual tidal debris structures are discussed.

To place the highly substructured stellar halos of the Milky Way and M31 in a larger context of hierarchical galaxy formation, it is necessary to understand the prevalence and properties of tidal substructure around external galaxies. This chapter details the current state of our observational knowledge of streams in galaxies in and beyond the Local Group, which are studied both in resolved stellar populations and in integrated light. Modeling of individual streams in extragalactic systems is hampered by our inability to obtain resolved stellar kinematics in the streams, though many streams contain alternate luminous kinematic tracers, such as globular clusters or planetary nebulae. We compare the observed structures to the predictions of models of galactic halo formation, which provide insight in the number and properties of streams expected around Milky Way like galaxies. More specifically, we discuss the inferences that can be made about stream progenitors based only on observed morphologies. We expand our discussion to consider hierarchical accretion at lower mass scales, in particular the observational evidence that substructure exists on smaller mass scales and the effects accretion events may have on the evolution of dwarf galaxies (satellite or isolated). Lastly, we discuss potential correlations between the presence of substructure in the halo and the structural properties of the disk. While many exciting discoveries have been made of tidal substructures around external galaxies, the "global" questions of galaxy formation and evolution via hierarchical accretion await a more complete census of the low surface brightness outskirts of galaxies in and beyond the Local Group.

We derive the fraction of substructure in the Galactic halo using a sample of over 10,000 spectroscopically-confirmed halo giant stars from the LAMOST spectroscopic survey. By observing 100 synthetic models along each line of sight with the LAMOST selection function in that sky area, we statistically characterize the expected halo populations. We define as SHARDS (Stellar Halo Accretion Related Debris Structures) any stars in >3-sigma excesses above the model predictions. We find that at least 10% of the Milky Way halo stars from LAMOST are part of SHARDS. By running our algorithm on smooth halos observed with the LAMOST selection function, we show that the LAMOST data contain excess substructure over all Galactocentric radii R_GC < 40 kpc, beyond what is expected due to statistical fluctuations and incomplete sampling of a smooth halo. The level of substructure is consistent with the fraction of stars in SHARDS in model halos created entirely from accreted satellites. This work illustrates the potential of vast spectroscopic surveys with high filling factors over large sky areas to recreate the merging history of the Milky Way.

We use test-particle integrations to show that epicyclic motions excited by a pericentre passage of a dwarf galaxy could account for bulk vertical velocity streaming motions recently observed in the Galactic stellar disc near the Sun. We use fixed potential test-particle integrations to isolate the role of phase wrapping of epicyclic perturbations from bending and breathing waves or modes, which require self-gravity to oscillate. Perturbations from a fairly massive Sagittarius dwarf galaxy, $M_d \sim 2.5 \times 10^{10} M_\odot$, are required to account for the size of the observed streaming motions from its orbital pericentre approximately a Gyr ago. A previous passage of the dwarf through the Galactic disc approximately 2.2 Gyr ago (with a then more massive dwarf galaxy) is less effective. If phase wrapping of epicyclic perturbations is responsible for stellar streaming motions in the Galactic disc, then there should be variations in velocity gradients on scales of a few kpc in the vicinity of the Sun.

We use about 200,000 FGK type main-sequence stars from the LAMOST DR1 data to map the local stellar kinematics. With the velocity de-projection technique, we are able to derive the averaged 3 dimensional velocity and velocity ellipsoids using only the line-of-sight velocity for the stars with various effective temperatures within $100 < |z| < 500$ pc. Using the mean velocities of the cool stars, we derive the solar motion of ($U_{\!\odot}$, $V_{\!\odot}$, $W_{\!\odot}$)=(9.58$\pm2.39$, 10.52$\pm1.96$, 7.01$\pm1.67$)$km\,s^{-1}$ with respect to the local standard of rest. Moreover, we find that the stars with ${T}_{\rm eff}>6000$K show a net asymmetric motion of $\sim3 km\,s^{-1}$ in $\langle W\rangle$ compared to the stars with ${T}_{\rm eff}<6000$K. And their azimuthal velocity increases when $|z|$ increases. This peculiar motion in the warmer stars is likely because they are young and not completely relaxed, although other reasons, such as the resonance induced by the central rotating bar or the spiral structures, and the perturbation of the merging dwarf galaxies, can not be ruled out. The derived velocity dispersions and cross terms for the data are approximately consistent with previous studies. We also find that the vertical gradients of $\sigma_{U}$ and $\sigma_V$ are larger than that of $\sigma_W$ . And the vertical gradient of $\sigma_U$ shows clear correlation with ${T}_{\rm eff}$, while the other two do not. Finally, our sample shows vertex deviation of about 11$^\circ$, at $300 < |z| < 500$pc, but roughly zero at $100 < |z| < 300$pc.

We present a method to estimate distances to stars with spectroscopically derived stellar parameters. The technique is a Bayesian approach with likelihood estimated via comparison of measured parameters to a grid of stellar isochrones, and returns a posterior probability density function for each star's absolute magnitude. This technique is tailored specifically to data from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) survey. Because LAMOST obtains roughly 3000 stellar spectra simultaneously within each ~5-degree diameter "plate" that is observed, we can use the stellar parameters of the observed stars to account for the stellar luminosity function and target selection effects. This removes biasing assumptions about the underlying populations, both due to predictions of the luminosity function from stellar evolution modeling, and from Galactic models of stellar populations along each line of sight. Using calibration data of stars with known distances and stellar parameters, we show that our method recovers distances for most stars within ~20%, but with some systematic overestimation of distances to halo giants. We apply our code to the LAMOST database, and show that the current precision of LAMOST stellar parameters permits measurements of distances with ~40% error bars. This precision should improve as the LAMOST data pipelines continue to be refined.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

We present a multi-epoch Hubble Space Telescope (HST) study of stellar proper motions (PMs) for four fields spanning 200 degrees along the Sagittarius (Sgr) stream: one trailing arm field, one field near the Sgr dwarf spheroidal tidal radius, and two leading arm fields. We determine absolute PMs of dozens of individual stars per field, using established techniques that use distant background galaxies as stationary reference frame. Stream stars are identified based on combined color-magnitude diagram and PM information. The results are broadly consistent with the few existing PM measurements for the Sgr galaxy and the trailing arm. However, our new results provide the highest PM accuracy for the stream to date, the first PM measurements for the leading arm, and the first PM measurements for individual stream stars; we also serendipitously determine the PM of the globular cluster NGC~6652. In the trailing-arm field, the individual PMs allow us to kinematically separate trailing-arm stars from leading-arm stars that are 360 degrees further ahead in their orbit. Also, in three of our fields we find indications that two distinct kinematical components may exist within the same arm and wrap of the stream. Qualitative comparison of the HST data to the predictions of the Law & Majewski and Penarrubia et al. N-body models show that the PM measurements closely follow the predicted trend with Sgr longitude. This provides a successful consistency check on the PM measurements, as well as on these N-body approaches (which were not tailored to fit any PM data).

The origin of the Galactic halo stellar structure known as the Monoceros ring is still under debate. In this work, we study that halo substructure using deep CFHT wide-field photometry obtained for the globular clusters NGC2419 and Koposov2, where the presence of Monoceros becomes significant because of their coincident projected position. Using Sloan Digital Sky Survey photometry and spectroscopy in the area surrounding these globulars and beyond, where the same Monoceros population is detected, we conclude that a second feature, not likely to be associated with Milky Way disk stars along the line-of-sight, is present as foreground population. Our analysis suggests that the Monoceros ring might be composed of an old stellar population of age t ~ 9Gyr and a new component ~ 4Gyr younger at the same heliocentric distance. Alternatively, this detection might be associated with a second wrap of Monoceros in that direction of the sky and also indicate a metallicity spread in the ring. The detection of such a low-density feature in other sections of this halo substructure will shed light on its nature.

We show that in the anticenter region, between Galactic longitudes of $110^\circ<l<229^\circ$, there is an oscillating asymmetry in the main sequence star counts on either side of the Galactic plane using data from the Sloan Digital Sky Survey. This asymmetry oscillates from more stars in the north at distances of about 2 kpc from the Sun to more stars in the south at 4-6 kpc from the Sun to more stars in the north at distances of 8-10 kpc from the Sun. We also see evidence that there are more stars in the south at distances of 12-16 kpc from the Sun. The three more distant asymmetries form roughly concentric rings around the Galactic center, opening in the direction of the Milky Way's spiral arms. The northern ring, 9 kpc from the Sun, is easily identified with the previously discovered Monoceros Ring. Parts of the southern ring at 14 kpc from the Sun (which we call the TriAnd Ring) have previously been identified as related to the Monoceros Ring and others have been called the Triangulum Andromeda Overdensity. The two nearer oscillations are approximated by a toy model in which the disk plane is offset by of the order 100 pc up and then down at different radii. We also show that the disk is not azimuthally symmetric around the Galactic anticenter and that there could be a correspondence between our observed oscillations and the spiral structure of the Galaxy. Our observations suggest that the TriAnd and Monoceros Rings (which extend to at least 25 kpc from the Galactic center) are primarily the result of disk oscillations.

Hypervelocity stars are those that have speeds exceeding the escape speed and are hence unbound from the Milky Way. We investigate a sample of low-mass hypervelocity candidates obtained using data from the high-precision SDSS Stripe 82 catalogue, which we have combined with spectroscopy from the 200-inch Hale Telescope at Palomar Observatory. We find four good candidates, but without metallicities it is difficult to pin-down their distances and therefore total velocities. Our best candidate has a significant likelihood that it is escaping the Milky Way for a wide-range of metallicities.

We present a census of the 12,060 spectra of blue objects ($(g-r)_0<-0.25$) in the Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8). As part of the data release, all of the spectra were cross-correlated with 48 template spectra of stars, galaxies and QSOs to determine the best match. We compared the blue spectra by eye to the templates assigned in SDSS DR8. 10,856 of the objects matched their assigned template, 170 could not be classified due to low signal-to-noise (S/N), and 1034 were given new classifications. We identify 7458 DA white dwarfs, 1145 DB white dwarfs, 273 rarer white dwarfs (including carbon, DZ, DQ, and magnetic), 294 subdwarf O stars, 648 subdwarf B stars, 679 blue horizontal branch stars, 1026 blue stragglers, 13 cataclysmic variables, 129 white dwarf - M dwarf binaries, 36 objects with spectra similar to DO white dwarfs, 179 QSOs, and 10 galaxies. We provide two tables of these objects, sample spectra that match the templates, figures showing all of the spectra that were grouped by eye, and diagnostic plots that show the positions, colors, apparent magnitudes, proper motions, etc. for each classification. Future surveys will be able to use templates similar to stars in each of the classes we identify to classify blue stars, including rare types, automatically.

We quantify and correct systematic errors in PPMXL proper motions using extragalactic sources from the first two LAMOST data releases and the Veron-Cetty & Veron Catalog of Quasars. Although the majority of the sources are from the Veron catalog, LAMOST makes important contributions in regions that are not well-sampled by previous catalogs, particularly at low Galactic latitudes and in the south Galactic cap. We show that quasars in PPMXL have measureable and significant proper motions, which reflect the systematic zero-point offsets present in the catalog. We confirm the global proper motion shifts seen by Wu, Ma, & Zhou (2011), and additionally find smaller-scale fluctuations of the QSO-derived corrections to an absolute frame. We average the proper motions of 158,106 extragalactic objects in bins of 3x3 degrees and present a table of proper motion corrections.

We present a support vector machine classifier to identify the K giant stars from the LAMOST survey directly using their spectral line features. The completeness of the identification is about 75% for tests based on LAMOST stellar parameters. The contamination in the identified K giant sample is lower than 2.5%. Applying the classification method to about 2 million LAMOST spectra observed during the pilot survey and the first year survey, we select 298,036 K giant candidates. The metallicities of the sample are also estimated with uncertainty of $0.13\sim0.29$\,dex based on the equivalent widths of Mg$_{\rm b}$ and iron lines. A Bayesian method is then developed to estimate the posterior probability of the distance for the K giant stars, based on the estimated metallicity and 2MASS photometry. The synthetic isochrone-based distance estimates have been calibrated using 7 globular clusters with a wide range of metallicities. The uncertainty of the estimated distance modulus at $K=11$\,mag, which is the median brightness of the K giant sample, is about 0.6\,mag, corresponding to $\sim30$% in distance. As a scientific verification case, the trailing arm of the Sagittarius stream is clearly identified with the selected K giant sample. Moreover, at about 80\,kpc from the Sun, we use our K giant stars to confirm a detection of stream members near the apo-center of the trailing tail. These rediscoveries of the features of the Sagittarius stream illustrate the potential of the LAMOST survey for detecting substructures in the halo of the Milky Way.

We report the first hypervelocity star (HVS) discovered from the LAMOST spectroscopic survey. It is a B-type star with a heliocentric radial velocity about 620 km/s, which projects to a Galactocentric radial velocity component of ~477 km/s. With a heliocentric distance of ~13 kpc and an apparent magnitude of ~13 mag, it is the nearest bright HVS currently known. With a mass of ~9Msun, it is one of the three most massive HVSs discovered so far. The star is clustered on the sky with many other known HVSs, with the position suggesting a possible connection to Galactic center structures. With the current poorly-determined proper motion, a Galactic center origin of this HVS remains consistent with the data at the 1sigma level, while a disk run-away origin cannot be excluded. We discuss the potential of the LAMOST survey to discover a large statistical sample of HVSs of different types.

We searched the light curves of over 40000 stars in the 2MASS calibration database, spanning approximately 4 years, for objects that have significant day long dimming events. We also searched the multi-color light curves for red dimming events that could be due to transient extinction. In the color independent sigma-limited search, we found 46 previously unknown eclipsing binaries, 6 periodic variable stars likely to be intrinsic pulsators and 21 young stellar objects in the $\rho$ Ophiucus star formation region previously studied by Parks et al. 2013. An additional 11 objects exhibited dimming events and most of these are unclassified. The search for red dimming events primarily reveals a population of low luminosity active galaxies that become bluer when they are brighter, and variable young stellar objects exhibiting high cross-correlation coefficients between color and brightness. The young stellar objects exhibit brightness and color variations in the direction of interstellar extinction whereas the active galaxies can have a bowed distribution in color and magnitude with reduced variation in color when the object is brightest. Among the objects that are usually quiescent (not strongly variable), we failed to find any dimming events deeper than 0.2 magnitude and lasting longer than a day. Two of the least embedded young stellar objects, however, dimmed by 0.2 mag for longer than a day without strong color variation. Having failed to find new exotic objects, we explored ways to eliminate commonly found objects so that a larger number of objects may be searched. We find that all but one eclipsing binary are excluded by requiring moderate color variation during a dimming event and most of the active galaxies and young stellar objects are excluded by placing a limit on the standard deviation of the magnitude distribution in the light curve.

We find that Galactic disk stars near the anticenter exhibit velocity asymmetries in both the Galactocentric radial and vertical components across the mid-plane as well as azimuthally. These findings are based on LAMOST spectroscopic velocities for a sample of ~400,000 F-type stars, combined with proper motions from the PPMXL catalog for which we have derived corrections to the zero points based in part on spectroscopically discovered galaxies and QSOs from LAMOST. In the region within 2 kpc outside the Sun's radius and +/-2 kpc from the Galactic midplane, we show that stars above the plane exhibit net outward radial motions with downward vertical velocities, while stars below the plane have roughly the opposite behavior. We discuss this in the context of other recent findings, and conclude that we are likely seeing the signature of vertical disturbances to the disk due to an external perturbation.

We trace the Cetus Polar Stream (CPS) with blue horizontal branch (BHB) and red giant stars (RGBs) from Data Release 8 of the Sloan Digital Sky Survey (SDSS DR8). Using a larger dataset than was available previously, we are able to refine the measured distance and velocity to this tidal debris star stream in the south Galactic cap. Assuming the tidal debris traces the progenitor's orbit, we fit an orbit to the CPS and find that the stream is confined between ~24-36 kpc on a rather polar orbit inclined 87 degrees to the Galactic plane. The eccentricity of the orbit is 0.20, and the period ~700 Myr. If we instead matched N-body simulations to the observed tidal debris, these orbital parameters would change by 10% or less. The CPS stars travel in the opposite direction to those from the Sagittarius tidal stream in the same region of the sky. Through N-body models of satellites on the best-fitting orbit, and assuming that mass follows light, we show that the stream width, line-of-sight depth, and velocity dispersion imply a progenitor of at least 10^8 solar masses. However, the density of stars along the stream requires either a disruption time on the order of one orbit, or a stellar population that is more centrally concentrated than the dark matter. We suggest that an ultra-faint dwarf galaxy progenitor could reproduce a large stream width and velocity dispersion without requiring a very recent deflection of the progenitor into its current orbit. We find that most Cetus stars have metallicities of -2.5 < [Fe/H] < -2.0, similar to the observed metallicities of the ultra-faint dwarfs. Our simulations suggest that the parameters of the dwarf galaxy progenitors, including their dark matter content, could be constrained by observations of their tidal tails through comparison of the debris with N-body simulations.

A total of $\sim640,000$ objects from LAMOST pilot survey have been publicly released. In this work, we present a catalog of DA white dwarfs from the entire pilot survey. We outline a new algorithm for the selection of white dwarfs by fitting S\'ersic profiles to the Balmer H$\beta$, H$\gamma$ and H$\delta$ lines of the spectra, and calculating the equivalent width of the CaII K line. 2964 candidates are selected by constraining the fitting parameters and the equivalent width of CaII K line. All the spectra of candidates are visually inspected. We identify 230 (59 of them are already in Villanova and SDSS WD catalog) DA white dwarfs, 20 of which are DA white dwarfs with non-degenerate companions. In addition, 128 candidates are classified as DA white dwarf/subdwarfs, which means the classifications are ambiguous. The result is consistent with the expected DA white dwarf number estimated based on the LEGUE target selection algorithm.

The LAMOST spectroscopic survey for Galactic structure and evolution has been in operation since October 2012, following a one-year pilot survey. The pilot survey produced a data release containing over 600,000 stellar spectra. By cross-checking with a large time series photometric database of RR Lyrae stars in high Galactic latitude regions, we found a total number of 157 stars that have been observed with LAMOST. In this sample, we successfully captured three RR Lyrae stars in the fast expansion phase, all of them showing significant hypersonic shock wave features in the spectra. By fitting the H{\alpha} line shape, we determine that the emission features seen within the broader H{\alpha} absorption line suggest hypersonic relative motions in the atmospheres of these three objects. With a further LAMOST survey of millions of stars, we will capture a large sample of RR Lyrae stars in their hypersonic expansion phase, and therefore provide a large database for the study of the internal structure and pulsation mechanism of RR Lyrae stars.

We report the kinematic discovery of the Pisces Stellar Stream (PSS), at Galactic longitude of roughly l=135 degrees and longitude between -39 < b < -36 deg. We originally identified this halo substructure from velocities of red giant branch stars in Sloan Digital Sky Survey (SDSS) Data Release 8, and confirmed its presence in turnoff stars from SDSS photometric data. The PSS is a narrow, kinematically cold tidal stream, with a velocity dispersion of 8 km/s. Its metallicity is [Fe/H]=-2.2, with ~0.3 dex dispersion. The color-magnitude signature of the stream turnoff, combined with our measured metallicity, places the PSS at a distance of 35+/-3 kpc. The Pisces Stellar Stream is the same as the previously announced "Triangulum stream" and part of the proposed "stream a". We rule out an association of the PSS with other previously known Milky Way substructures in the same region of the sky.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

We outline the design of the dark nights portion of the LAMOST Pilot Survey, which began observations in October 2011. In particular, we focus on Milky Way stellar candidates that are targeted for the LEGUE (LAMOST Experiment for Galactic Understanding and Exploration) survey. We discuss the regions of sky in which spectroscopic candidates were selected, and the motivations for selecting each of these sky areas. Some limitations due to the unique design of the telescope are discussed, including the requirement that a bright (V < 8) star be placed at the center of each plate for wavefront sensing and active optics corrections. The target selection categories and scientific goals motivating them are briefly discussed, followed by a detailed overview of how these selection functions were realized. We illustrate the difference between the overall input catalog - Sloan Digital Sky Survey (SDSS) photometry - and the final targets selected for LAMOST observation.

The weather at Xinglong Observing Station, where the Guo Shou Jing Telescope (GSJT) is located, is strongly affected by the monsoon climate in north-east China. The LAMOST survey strategy is constrained by these weather patterns. In this paper, we present a statistics on observing hours from 2004 to 2007, and the sky brightness, seeing, and sky transparency from 1995 to 2011 at the site. We investigate effects of the site conditions on the survey plan. Operable hours each month shows strong correlation with season: on average there are 8 operable hours per night available in December, but only 1-2 hours in July and August. The seeing and the sky transparency also vary with seasons. Although the seeing is worse in windy winters, and the atmospheric extinction is worse in the spring and summer, the site is adequate for the proposed scientific program of LAMOST survey. With a Monte Carlo simulation using historical data on the site condition, we find that the available observation hours constrain the survey footprint from 22h to 16h in right ascension; the sky brightness allows LAMOST to obtain the limit magnitude of V = 19.5mag with S/N = 10.

We describe the footprint and input catalog for bright nights in the LAMOST Pilot Survey, which began in October 2011. Targets are selected from two stripes in the north and south Galactic Cap regions, centered at $\alpha$= 29$^\circ$, with 10$^\circ$ width in declination, covering right ascension of 135$^\circ-290^\circ$ and -30$^\circ$ to 30$^\circ$ respectively. We selected spectroscopic targets from a combination of the SDSS and 2MASS point source catalogs. The catalog of stars defining the field centers (as required by the Shack-Hartmann wavefront sensor at the center of the LAMOST field) consists of all V < 8m stars from the Hipparcos catalog. We employ a statistical selection algorithm that assigns priorities to targets based on their positions in multidimensional color/magnitude space. This scheme overemphasizes rare objects and de-emphasizes more populated regions of magnitude and color phase space, while ensuring a smooth, well-understood selection function. A demonstration of plate design is presented based on the Shack-Hartmann star catalog and an input catalog that was generated by our target selection routines.

We describe the current plans for a spectroscopic survey of millions of stars in the Milky Way galaxy using the Guo Shou Jing Telescope (GSJT, formerly the Large Area Multi-Object Spectroscopic Telescope - LAMOST). The survey will obtain spectra for 2.5 million stars brighter than $r<19$ during dark/grey time, and 5 million stars brighter than $r<17$ or $J<16$ on nights that are moonlit or have low transparency. The survey will begin in fall of 2012, and will run for at least four years. The telescope design constrains the optimal declination range for observations to $10^\circ<\delta<50^\circ$, and site conditions lead to an emphasis on stars in the direction of the Galactic anticenter. The survey is divided into three parts with different target selection strategies: disk, anticenter, and spheroid. The resulting dataset will be used to study the merger history of the Milky Way, the substructure and evolution of the disks, the nature of the first generation of stars through identification of the lowest metallicity stars, and star formation through study of open clusters and the OB associations. Detailed design of the LEGUE survey will be completed after a review of the results of the pilot survey in summer 2012.

We describe the target selection algorithm for the low latitude disk portion of the LAMOST Pilot Survey, which aims to test systems in preparation for the LAMOST spectroscopic survey. We use the PPMXL (Roeser et al. 2010) astrometric catalog, which provides positions, proper motions, B/R/I magnitudes (mostly) from USNO-B (Monet et al. 2003) and J/H/Ks from The Two Micron All Sky Survey (2MASS, see Skrutskie et al. 2006) as well. We chose 8 plates along the Galactic plane, in the region $0^\circ<\alpha<67^\circ$ and $42^\circ<\delta<59^\circ$, that cover 22 known open clusters with a range of ages. Adjacent plates may have small overlapping. Each plate covers an area $2.5^\circ$ in radius,with central star (for Shack-Hartmann guider) brighter than $\sim8^{\rm th}$ magnitude. For each plate, we create an input catalog in the magnitude range $11.3<Imag<16.3$ and $Bmag$ available from PPMXL. The stars are selected to satisfy the requirements of the fiber positioning system and have a uniform distribution in the $I$ vs. $B-I$ color-magnitude diagram. Our final input catalog consists of 12,000 objects on each of 8 plates that are observable during the winter observing season in Xinglong Station of the National Astronomical Observatory of China.

We confirm that there are at least three separate low-latitude over-densities of blue F turnoff stars near the Milky Way anti-center: the Monoceros Ring, the Anti-Center Stream (ACS), and the Eastern Banded Structure (EBS). There might also be a small number of normal thick disk stars at the same location. The ACS is a tilted component that extends to higher Galactic latitude at lower Galactic longitude, 10 kpc from the Sun towards the anti-center. It has a sharp cutoff on the high latitude side. Distance, velocity, and proper motion measurements are consistent with previous orbit fits. The mean metallicity is [Fe/H]$=-0.96 \pm 0.03$, which is lower than the thick disk and Monoceros Ring. The Monoceros Ring is a higher density substructure that is present at $15\arcdeg<b<22\arcdeg$ at all longitudes probed in this survey. The structure likely continues towards lower latitudes. The distances are consistent with a constant distance from the Galactic Center of 17.6 kpc. The mean line-of-sight velocity of the structure is consistent with a thick disk rotation. However, the velocity dispersion of these stars is $\sim 15$ km s$^{-1}$, and the metallicity is [Fe/H]$=-0.80 \pm 0.01$. Both of these quantities are lower than the canonical thick disk. We suggest that this ring structure is likely different from the thick disk, though its association with the disk cannot be definitively ruled out. The Eastern Banded Structure (EBS) is detected primarily photometrically, near $(l,b)=(225\arcdeg,30\arcdeg)$, at a distance of 10.9 kpc from the Sun.

We describe a general target selection algorithm that is applicable to any survey in which the number of available candidates is much larger than the number of objects to be observed. This routine aims to achieve a balance between a smoothly-varying, well-understood selection function and the desire to preferentially select certain types of targets. Some target-selection examples are shown that illustrate different possibilities of emphasis functions. Although it is generally applicable, the algorithm was developed specifically for the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) survey that will be carried out using the Chinese Guo Shou Jing Telescope. In particular, this algorithm was designed for the portion of LEGUE targeting the Galactic halo, in which we attempt to balance a variety of science goals that require stars at fainter magnitudes than can be completely sampled by LAMOST. This algorithm has been implemented for the halo portion of the LAMOST pilot survey, which began in October 2011.

We present three-dimensional space velocities of stars selected to be consistent with membership in the Virgo stellar substructure. Candidates were selected from SA 103, a single 40x40 arcmin field from our proper motion (PM) survey in Kapteyn's Selected Areas (SAs), based on the PMs, SDSS photometry, and follow-up spectroscopy of 215 stars. The signature of the Virgo substructure is clear in the SDSS color-magnitude diagram (CMD) centered on SA 103, and 16 stars are identified that have high Galactocentric-frame radial velocities (V_GSR > 50 km/s) and lie near the CMD locus of Virgo. The implied distance to the Virgo substructure from the candidates is 14+/-3 kpc. We derive mean kinematics from these 16 stars, finding a radial velocity V_GSR = 153+/-22 km/s and proper motions (mu_alpha*cos(delta), mu_delta) = (-5.24, -0.91)+/-(0.43, 0.46) mas/yr. From the mean kinematics of these members, we determine that the Virgo progenitor was on an eccentric (e ~ 0.8) orbit that recently passed near the Galactic center (pericentric distance R_p ~ 6 kpc). This destructive orbit is consistent with the idea that the substructure(s) in Virgo originated in the tidal disruption of a Milky Way satellite. N-body simulations suggest that the entire cloud-like Virgo substructure (encompassing the "Virgo Overdensity" and the "Virgo Stellar Stream") is likely the tidal debris remnant from a recently-disrupted massive (~10^9 M_sun) dwarf galaxy. The model also suggests that some other known stellar overdensities in the Milky Way halo (e.g., the Pisces Overdensity and debris near NGC 2419 and SEGUE 1) are explained by the disruption of the Virgo progenitor.

We study the shapes of Milky Way satellites in the context of the tidal stirring scenario for the formation of dwarf spheroidal galaxies. The standard procedures used to measure shapes involve smoothing and binning of data and thus may not be sufficient to detect structural properties like bars, which are usually subtle in low surface brightness systems. Taking advantage of the fact that in nearby dwarfs photometry of individual stars is available we introduce discrete measures of shape based on the two-dimensional inertia tensor and the Fourier bar mode. We apply these measures of shape first to a variety of simulated dwarf galaxies formed via tidal stirring of disks embedded in dark matter halos and orbiting the Milky Way. In addition to strong mass loss and randomization of stellar orbits, the disks undergo morphological transformation that typically involves the formation of a triaxial bar after the first pericenter passage. These tidally induced bars persist for a few Gyr before being shortened towards a more spherical shape if the tidal force is strong enough. We test this prediction by measuring in a similar way the shape of nearby dwarf galaxies, satellites of the Milky Way. We detect inner bars in Ursa Minor, Sagittarius, LMC and possibly Carina. In addition, six out of eleven studied dwarfs show elongated stellar distributions in the outer parts that may signify transition to tidal tails. We thus find the shapes of Milky Way satellites to be consistent with the predictions of the tidal stirring model.

We present three-dimensional kinematics of Sagittarius (Sgr) trailing tidal debris in six fields located 70-130 degrees along the stream from the Sgr dwarf galaxy core. The data are from our proper-motion (PM) survey of Kapteyn's Selected Areas, in which we have measured accurate PMs to faint magnitudes in 40x40 arcmin fields evenly spaced across the sky. The radial velocity (RV) signature of Sgr has been identified among our follow-up spectroscopic data in four of the six fields and combined with mean PMs of spectroscopically-confirmed members to derive space motions of Sgr debris based on 15-64 confirmed stream members per field. These kinematics are compared to predictions of the Law & Majewski (2010) model of Sgr disruption; we find reasonable agreement with model predictions in RVs and PMs along Galactic latitude. However, an upward adjustment of the Local Standard of Rest velocity Theta_LSR from its standard 220 km/s to at least $232\pm14$ km/s (and possibly as high as $264\pm23$ km/s) is necessary to bring 3-D model debris kinematics and our measurements into agreement. Satisfactory model fits that simultaneously reproduce known position, distance, and radial velocity trends of the Sgr tidal streams, while significantly increasing Theta_LSR}, could only be achieved by increasing the Galactic bulge and disk mass while leaving the dark matter halo fixed to the best-fit values from Law & Majewski (2010). We derive low-resolution spectroscopic abundances along this stretch of the Sgr stream and find a constant [Fe/H] = -1.15 (with ~0.5 dex scatter in each field -- typical for dwarf galaxy populations) among the four fields with reliable measurements. A constant metallicity suggests that debris along the ~60-degree span of this study was all stripped from Sgr on the same orbital passage.

We have measured the mean three-dimensional kinematics of stars in Kapteyn's Selected Area (SA) 76 (l=209.3, b=26.4 degrees) that were selected to be Anticenter Stream (ACS) members on the basis of their radial velocities, proper motions, and location in the color-magnitude diagram. From a total of 31 stars ascertained to be ACS members primarily from its main sequence turnoff, a mean ACS radial velocity (derived from spectra obtained with the Hydra multi-object spectrograph on the WIYN 3.5m telescope) of V_helio = 97.0 +/- 2.8 km/s was determined, with an intrinsic velocity dispersion sigma_0 = 12.8 \pm 2.1 km/s. The mean absolute proper motions of these 31 ACS members are mu_alpha cos (delta) = -1.20 +/- 0.34 mas/yr and mu_delta = -0.78 \pm 0.36 mas/yr. At a distance to the ACS of 10 \pm 3 kpc, these measured kinematical quantities produce an orbit that deviates by ~30 degrees from the well-defined swath of stellar overdensity constituting the Anticenter Stream in the western portion of the Sloan Digital Sky Survey footprint. We explore possible explanations for this, and suggest that our data in SA 76 are measuring the motion of a kinematically cold sub-stream among the ACS debris that was likely a fragment of the same infalling structure that created the larger ACS system. The ACS is clearly separated spatially from the majority of claimed Monoceros ring detections in this region of the sky; however, with the data in hand, we are unable to either confirm or rule out an association between the ACS and the poorly-understood Monoceros structure.

We report the results of a spectroscopic study of the Bootes III (BooIII) stellar overdensity carried out with the Hectospec multifiber spectrograph on the MMT telescope. Radial velocities have been measured for 193 BooIII candidate stars selected to have magnitudes and colors consistent with its upper main sequence and lower red giant branch, as well as a number of horizontal branch candidates. From 20 identified candidate BooIII members, we measure a systemic velocity of V_sun=197.5+-3.8 km/s and a velocity dispersion of sigma_o=14.0+-3.2 km/s. We use the somewhat large velocity dispersion and the implied highly radial orbit, along with morphological evidence from Grillmair (2009) and stellar abundances, to argue that BooIII is likely the first known object observed in a transitional state between being a bound dwarf galaxy and a completely unbound tidal stream.

We present absolute proper motions in Kapteyn Selected Area (SA) 103. This field is located 7 degrees west of the center of the Virgo Stellar Stream (VSS, Duffau et al. 2006), and has a well-defined main sequence representing the stream. In SA 103 we identify one RR Lyrae star as a member of the VSS according to its metallicity, radial velocity and distance. VSS candidate turnoff stars and subgiant stars have proper motions consistent with that of the RR Lyrae star. The 3D velocity data imply an orbit with a pericenter of 11 kpc and an apocenter of ~90 kpc. Thus, the VSS comprises tidal debris found near the pericenter of a highly destructive orbit. Examining the six globular clusters at distances larger than 50 kpc from the Galactic center, and the proposed orbit of the VSS, we find one tentative association, NGC 2419. We speculate that NGC 2419 is possibly the nucleus of a disrupted system of which the VSS is a part.

We have obtained radial velocity measurements for stars in two, widely-separated fields in the Anticenter Stream. Combined with SDSS/USNO-B proper motions, the new measurements allow us to establish that the stream is on a nearly circular, somewhat inclined, prograde orbit around the Galaxy. While the orbital eccentricity is similar to that previously determined for the Monoceros stream, the sizes, inclinations, and positions of the orbits for the two systems differ significantly. Integrating our best fitting Anticenter Stream orbit forward, we find that it is closely aligned along and lies almost on top of a stream-like feature previously designated the "Eastern Banded Structure". The position of this feature coincides with the apogalacticon of the orbit. We tentatively conclude that this feature is the next wrap of the Anticenter Stream.

We present a 3D kinematical analysis of stars located in Kapteyn Selected Area 71 (l = 167.1, b = -34.7), where previously a stellar excess was found (Dinescu et al. 2002, Newberg et al. 2002). Previous findings indicated that the stellar excess has a cold kinematical signature as inferred from proper motions, and was initially associated with debris from the Sagittarius dwarf galaxy (Sgr) -- namely the southern trailing tail. We have obtained radial velocities using the Hydra multiobject spectrograph on the WIYN 3.5 m telescope. Results for 183 proper-motion selected stars indicate that the dominant population in this stellar excess is not debris from Sgr, but rather a population that kinematically belongs to the ring-like stream that is now known as the Monoceros stream. The orbit determined for this population agrees very well with the predictions for the Monoceros stream from Penarrubia et al. (2005). The radial-velocity dispersion of this population is between 20 and 30 km/s, lower than that of the Galactic field. Also, the shape of the radial-velocity distribution shows a sharp cut-off on one side, which is more in line with model predictions of the disruption of a satellite rather than with the distribution of the Galactic field. Despite the fact that we now believe most of the stars in the stellar excess to be part of Monoceros, about ten stars in this stellar excess have highly negative radial velocities, which is a clear indication of their membership to the Sgr trailing tidal tail.

We report the results of a spectroscopic study of the Bootes (Boo) dwarf spheroidal (dSph) galaxy carried out with the WIYN telescope and the Hydra multifiber spectrograph. Radial velocities have been measured for 58 Boo candidate stars selected to have magnitudes and colors consistent with its red and asymptotic giant branches. Within the 13' half-light radius, seven members of Boo yield a systemic velocity of V_r=95.6+-3.4 km/s and a velocity dispersion of 6.6+-2.3 km/s. This implies a mass on the order of 1 x 10^7 M_sun, similar to the inferred masses of other Galactic dSphs. Adopting a total Boo luminosity of L=1.8 x 10^4 L_sun to 8.6 x 10^4 L_sun implies M/L ~ 610 to 130, making Boo, the most distorted known Milky Way dwarf galaxy, potentially also the darkest. From the spectra of Boo member stars we estimate its metallicity to be [Fe/H] ~ -2.5, which would make it the most metal poor dSph known to date.

We describe a high-precision, deep (to V ~ 19-21) absolute proper-motion survey that samples ~50 lines of sight in the Kapteyn Selected Areas along declination zones -15, 0 and 15 degrees. In many fields the astrometric baseline reaches nearly a century. We demonstrate that these data provide typical per star precisions between ~ 1 and 3 mas/yr to the above magnitude limits, with the absolute reference frame established by numerous extragalactic sources in each survey field. Combined with existing and ongoing photometric and radial velocity surveys in these fields, these astrometric data will enable, among other things, accurate, detailed dynamical modeling of satellite interactions with our Galaxy. In this contribution we describe the astrometric part of our survey and show preliminary results along the trailing tail of the Sagittarius dwarf galaxy, and in the Monoceros ring region.

(Abridged) We observed the X-ray bright E3 galaxy NGC 1600 and nearby group members with the Chandra X-ray Observatory ACIS S3. Unresolved emission dominates the observation; however, we resolved some emission into 71 sources, most of which are low-mass X-ray binaries (LMXBs) associated with NGC 1600. Twenty-one sources have L_X>2E39 ergs/s (0.3-10.0 keV; assuming they are at the distance of NGC 1600) marking them as ultraluminous X-ray point source (ULX) candidates; we expect only 11+/-2 are unrelated sources. NGC 1600 may have the largest number of ULX candidates in an early-type galaxy to date; however, cosmic variance in the number of background AGN cannot be ruled out. There are indications that there are also a large number of unresolved point sources. We propose that these are associated with globular clusters (GCs), and that NGC 1600 has a large GC specific frequency. The unresolved emission appears to have two gas components. We interpret the soft inner component (a<~25", kT~0.85 keV) as the interstellar medium of NGC 1600 and the hotter outer component (a>~25", kT~1.5 keV) as the intragroup medium of the NGC 1600 group. The X-ray image shows several interesting structures, excess emission roughly cospatial with H$\alpha$ and dust filaments, possible holes in the X-ray emission roughly coincident with radio lobes, excess emission to the northeast possibly indicating the center of the group potential, and a tail of X-ray emission from the group galaxy NGC 1603 probably due to ram-pressure stripping