Search

Whiting 1 is a faint and young globular cluster in the halo of the Milky Way, and was suggested to have originated in the Sagittarius spherical dwarf galaxy (Sgr dSph). In this paper, we use the deep DESI Legacy Imaging Surveys to explore tentative spatial connection between Whiting 1 and the Sgr dSph. We redetermine the fundamental parameters of Whiting 1 and use the best-fitting isochrone (age $\tau$=6.5 Gyr, metalicity Z=0.005 and $\rm d_{\odot}$=26.9 kpc) to construct a theoretical matched filter for the extra-tidal features searching. Without any smooth technique to the matched filter density map, we detect a round-shape feature with possible leading and trailing tails on either side of the cluster. This raw image is not totally new compared to old discoveries, but confirms that no more large-scale features can be detected under a depth of r<=22.5 mag. In our results, the whole feature stretches 0.1-0.2 degree along the orbit of Whiting 1, which gives a much larger area than the cluster core. The tails on both sides of the cluster align along the orbital direction of the Sgr dSph as well as the cluster itself, which implies that these debris are probably stripped remnants of Whiting 1 by the Milky Way., Comment: 10 pages, 9 figures, accepted for the publication of ApJ

We collect a sample of stars observed both in LAMOST and Gaia which have colors implying a temperature hotter than 7000 K. We train a machine learning algorithm on LAMOST spectroscopic data which has been tagged with stellar classifications and metallicities, and use this machine to construct a catalog of Blue Horizontal Branch stars (BHBs) with metallicity information. Another machine is trained using Gaia parallaxes to predict absolute magnitudes for these stars. The final catalog of 13,693 BHBs is thought to be about 86\% pure, with $\sigma_{[Fe/H]}\sim$0.35 dex and $\sigma_{G}\sim$0.31 mag. These values are confirmed via comparison to globular clusters, although a covariance error seems to affect our magnitude and abundance estimates. We analyze a subset of this catalog in the Galactic Halo. We find that BHB populations in the outer halo appear redder, which could imply a younger population, and that the metallicity gradient is relatively flat around [Fe/H] = -1.9 dex over our sample footprint. We find that our metal rich BHB stars are on more radial velocity dispersion dominated orbits ($\beta \sim 0.70$) at all radii than our metal poor BHB stars ($\beta \sim 0.62$)., Comment: 13 Pages, 7 Figures, accepted to apj

Although the stellar halo accounts for just ~1% of the total stellar mass of the Milky Way, the kinematics of halo stars encode valuable information about the origins and evolution of our Galaxy. It has been shown that the high transverse velocity stars in Gaia DR2 reveal a double sequence in the Hertzsprung-Russell (HR) diagram, indicating a bifurcation in the local stellar halo within 1 kpc. We fit these stars by updating the popular Besan\c{c}on/Galaxia model, incorporating the latest observational results for the stellar halo and an improved kinematic description for the thick-disk from Sch\"onrich & Binney (2012). We are able to obtain a good match to the Gaia data and provide new constraints on the properties of the Galactic disc and stellar halo. In particular, we show that the kinematically defined thick disc contribution to this high velocity tail is ~13%. We look in greater detail using chemistry from LAMOST DR5, identifying a population of retrograde stars with thick-disc chemistry. Our thick disc kinematic model cannot account for this population and so we conclude there is likely to be a contribution from heated or accreted stars in the Solar Neighbourhood. We also investigate proposed dynamical substructures in this sample, concluding that they are probably due to resonant orbits rather than accreted populations. Finally we provide new insights on the nature of the two sequences and their relation with past accretion events and the primordial Galactic disc., Comment: Accepted for publication in MNRAS

The Milky Way's progenitor experienced several merger events which left their imprints on the stellar halo, including the Gaia-Sausage/Enceladus. Recently, it has been proposed that this event perturbed the proto-disk and gave rise to a metal rich ([Fe/H] $>-1$) low angular momentum ($v_{\phi} < 100$ km/s) stellar population. These stars have dynamical and chemical properties different from the accreted stellar halo, but are continuous with the canonical thick disk. In this letter, we use a hydrodynamical simulation of an isolated galaxy which develops clumps that produce a bimodal thin$+$thick disk chemistry to explore whether it forms such a population. We demonstrate clump scattering forms a metal-rich, low angular momentum population, without the need for a major merger. We show that, in the simulation, these stars have chemistry, kinematics and density distribution in good agreement with those in the Milky Way., Comment: Accepted for publication in ApJL

Identifying the structure of our Galaxy has always been fraught with difficulties, and while modern surveys continue to make progress building a map of the Milky Way, there is still much to understand. The arm and bar features are important drivers in shaping the interstellar medium, but their exact nature and influence still require attention. We present results of smoothed particle hydrodynamic simulations of gas in the Milky Way including star formation, stellar feedback, and ISM cooling, when exposed to different arm and bar features, with the aim of better understanding how well newly formed stars trace out the underlying structure of the Galaxy. The bar is given a faster pattern speed than the arms, resulting in a complex, time-dependent morphology and star formation. Inter-arm branches and spurs are easily influenced by the bar, especially in the two-armed spiral models where there is a wide region of resonance overlap in the disc. As the bar over-takes the spiral arms it induces small boosts in star formation and enhances spiral features, which occur at regularly spaced beat-like intervals. The locations of star formation events are similar to those seen in observational data, and do not show a perfect 1:1 correspondence with the underlying spiral potential, though arm tangencies are generally well traced by young stars. Stellar velocity fields from the newly formed stars are compared to data from Gaia DR2, showing that the spiral and bar features can reproduce many of the non-axisymmetric features seen in the data. A simple analytical model is used to show many of these feature are a natural response of gas to rigidly rotating spiral and bar potentials., Comment: 21 pages, 26 figures, accepted for publication in MNRAS. Movies of simulations included in this work can be found at: https://zenodo.org/record/3239994

Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the "thin" and "thick" discs and the halo. We show that in the Solar neighborhood, there exists a large population of metal-rich ([Fe/H]>-0.7) stars on highly eccentric orbits. By studying the evolution of elemental abundances, kinematics and stellar ages in the plane of azimuthal velocity v_phi and metallicity [Fe/H], we demonstrate that this metal-rich halo-like component, which we dub the Splash, is linked to the alpha-rich (or "thick") disc. Splash stars have little to no angular momentum and many are on retrograde orbits. They are predominantly old, but not as old as the stars deposited into the Milky Way in the last major merger. We argue, in agreement with several recent studies, that the Splash stars may have been born in the Milky Way's proto-disc prior to the massive ancient accretion event which drastically altered their orbits. We can not, however, rule out other (alternative) formation channels. Taking advantage of the causal connection between the merger and the Splash, we put constraints of the epoch of the last massive accretion event to have finished 9.5 Gyr ago. The link between the local metal-rich and metal-poor retrograde stars is confirmed using a large suite of cutting-edge numerical simulations of the Milky Way's formation., Comment: Submitted to MNRAS. Comments are welcome!

Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly populated of old stars ($\sim 10\,\rm{Gyr}$), whereas Monoceros has a pronounced tail of younger stars ($6-10\, \rm{Gyr}$) as revealed by their cumulative age distributions. Put togehter, all of this evidence support predictions from simulations of the interaction of the Sagittarius dwarf with the Milky Way, which argue that the Anticenter Stream (ACS) is the remains of a tidal tail of the Galaxy excited during Sgr's first pericentric passage after it crossed the virial radius, whereas Monoceros consists of the composite stellar populations excited during the more extended phases of the interaction. We suggest that the ACS can be used to constrain the Galactic potential, particularly its flattening, setting strong limits on the existence of a dark disc. Importantly, the ACS can be viewed as a stand-alone fossil of the chemical enrichment history of the Galactic disc., Comment: 6 pages, 6 Figure, submitted to MNRAS, comments welcome

Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces Overdensity, a diffuse sub-structure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces's elongated shape is aligned with the Magellanic Stream. Using follow-up VLT FORS2 spectroscopy, we have measured the velocity distribution of the Pisces candidate member stars and have shown it to be as broad as that of the Galactic halo but offset to negative velocities. Using a suite of numerical simulations, we demonstrate that the structure has many properties in common with the predicted behaviour of the Magellanic wake, i.e. the Galactic halo overdensity induced by the in-fall of the Magellanic Clouds., Comment: Submitted to MNRAS Letters. Comments welcome!

In this paper we investigate how the chemical and kinematic properties of stars vary as a function of age. Using data from a variety of photometric, astrometric and spectroscopic surveys, we calculate the ages, phase space information and orbits for $\sim$125,000 stars covering a wide range of stellar parameters. We find indications that the inner regions of the disk reached high levels of enrichment early, while the outer regions were more substantially enriched in intermediate and recent epochs. We consider these enrichment histories through comparison of the ages of stars, their metallicities, and kinematic properties, such as their angular momentum in the solar neighborhood (which is a proxy for orbital radius). We calculate rates at which the velocity dispersions evolve, investigate the Oort constants for different aged populations (finding a slightly negative $\partial V_{C} / \partial R$ and $\partial V_{R} / \partial R$ for all ages, being most negative for the oldest stars), as well as examine the behavior of the velocity vertex deviation angle as a function of age (which we find to fall from $\sim$15 degrees for the 2 Gyr aged population to $\sim$6 degrees at around 6.5 Gyr of age, after which it remains unchanged). We find evidence for stellar churning, and find that the churned stars have a slightly younger age distribution than the rest of the data., Comment: 18 Pages, 14 Figures, Accepted ApJ

The stellar velocity ellipsoid of the solar neighbour is re-examined using intermediate-old mono-abundances stellar groups with high quality chemistry data together with parallaxes and proper motions from Gaia DR1. We find the average velocity dispersion values for the three space velocity components for the thin and thick disc of (\sigma_{U},\sigma_{V},\sigma_{W})_{thin} = (33 \pm 4, 28 \pm 2, 23 \pm 2) and (\sigma_{U},\sigma_{V},\sigma_{W})_{thick} = (57 \pm 6, 38 \pm 5, 37 \pm 4) km s^{-1}, respectively. The mean values of the ratio between the semi-axes of the velocity ellipsoid for the thin disc are found to be, \sigma_{V}/\sigma_{U} = 0.70 \pm 0.13 and \sigma_{W}/\sigma_{U} is 0.64 \pm 0.08, while for the thick disc \sigma_{V}/\sigma_{U} = 0.67 \pm 0.11 and \sigma_{W}/\sigma_{U} is 0.66 \pm 0.11. Inputting these dispersions into the linear Str\"omberg relation for the thin disc groups, we find the Sun's velocity with respect to the LSR in Galactic rotation to be V_{\sun} = 13.9 \pm 3.4 km s^{-1}. A relation is found between the vertex deviation and the chemical abundances for the thin disc, ranging from -5 to +40^{\circ} as iron-abundance increases. For the thick disc we find a vertex deviation of l_{uv} \sim -15^{\circ}. The tilt angle (l_{uw}) in the U-W plane for the thin disc groups ranges from -10 to +15^\circ, but there is no evident relation between l_{uw} and the mean abundances. However we find a weak relation for l_{uw} as a function of iron abundances and \alpha-elements for most of the groups in the thick disc, where the tilt angle decreases from -5 to -20^\circ when [Fe/H] decreases and [\alpha/Fe] increases. The velocity anisotropy parameter is independent of the chemical group abundances and its value is nearly constant for both discs (\beta \sim 0.5), suggesting that the combined disc is dynamically relaxed., Comment: 13 pages, 13 figures, accepted for publication in MNRAS

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., Comment: 13 pages, accepted for publication in ApJ

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling., Comment: 23 pages; accepted by ApJ

Clues to the origins and evolution of our Galaxy can be found in the kinematics of stars around us. Remnants of accreted satellite galaxies produce over- densities in velocity-space, which can remain coherent for much longer than spatial over-densities. This chapter reviews a number of studies that have hunted for these accretion relics, both in the nearby solar-neighborhood and the more-distant stellar halo. Many observational surveys have driven this field forwards, from early work with the Hipparcos mission, to contemporary surveys like RAVE & SDSS. This active field continues to flourish, providing many new discoveries, and will be revolutionised as the Gaia mission delivers precise proper motions for a billion stars in our Galaxy., Comment: 27 pages, 10 figures. Chapter from Springer ASSL Volume entitled "Tidal Streams in the Local Group and Beyond". Affluent readers may wish to purchase the full volume here: http://link.springer.com/book/10.1007/978-3-319-19336-6

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., Comment: 15 pages, 13 figures, accepted by ApJ. Minor issue fixed with Fig 12 and Section 5.3

Using LAMOST spectroscopic data, we find a strong signal of a comoving group of stars in the constellation of Draco. The group, observed near the apocenter of its orbit, is 2.6 kpc from the Sun with a metallicity of -0.64 dex. The system is observed as a streaming population of unknown provenance with mass of about 2.1E4 solar masses and an absolute V band magnitude of about -3.6. Its high metallicity, diffuse physical structure, and eccentric orbit may indicate that the progenitor satellite was a globular cluster rather than a dwarf galaxy or an open cluster., Comment: 6 pages, 4 Figures, 1 Table, Accepted to ApJL

In this paper we analyze a sample of metal-rich (>-0.8 dex) main sequence stars in the extended solar neighborhood, investigating kinematic outliers from the background population. The data, which are taken from the Sloan Digital Sky Survey, are kinematically profiled as a function of distance from the Galactic plane using full six dimensional phase space information. Each star is examined in the context of these kinematic profiles and likelihoods are assigned to quantify whether a star matches the underlying profile. Since some of these stars are likely to have been ejected from the disc, we trace back their orbits in order to determine potential ejection radii. We find that objects with low probability (i.e. `outliers') are typically more metal poor, faster and, most importantly, have a tendency to originate from the inner Galaxy compared to the underlying population. We also compose a sample of stars with velocities exceeding the local escape velocity. Although we do not discount that our sample could be contaminated by objects with spurious proper motions, a number of stars appear to have been ejected from the disc with exceptionally high velocities. Some of these are consistent with being ejected from the spiral arms and hence are a rich resource for further study. Finally we look at objects whose orbits are consistent with them being ejected at high speeds from the Galactic center. Of these objects we find that one, J135855.65+552538.19, is inconsistent with halo, bulge and disk kinematics and could plausibly have been ejected from the Galactic nucleus via a Hills mechanism., Comment: 17 Pages, 12 Figures, Accepted to AJ

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., Comment: 11 pages, 12 figures; accepted for publication in AJ

We extract the resonant orbits from an N-body bar that is a good representation of the Milky Way, using the method recently introduced by Molloy et al. (2015). By decomposing the bar into its constituent orbit families, we show that they are intimately connected to the boxy-peanut shape of the density. We highlight the imprint due solely to resonant orbits on the kinematic landscape towards the Galactic centre. The resonant orbits are shown to have distinct kinematic features and may be used to explain the cold velocity peak seen in the APOGEE commissioning data (Nidever at al., 2012). We show that high velocity peaks are a natural consequence of the motions of stars in the 2:1 orbit family and that stars on other higher order resonances can contribute to the peaks. The locations of the peaks vary with bar angle and, with the tacit assumption that the observed peaks are due to the 2:1 family, we find that the locations of the high velocity peaks correspond to bar angles in the range 10 < theta_bar < 25 (deg). However, some important questions about the nature of the peaks remain, such as their apparent absence in other surveys of the Bulge and the deviations from symmetry between equivalent fields in the north and south. We show that the absence of a peak in surveys at higher latitudes is likely due to the combination of a less prominent peak and a lower number density of bar supporting orbits at these latitudes., Comment: 7 Figures, 1 Table, Now includes figures & discussion of higher order resonances, Minor revisions to text throughout, Conclusions unchanged

We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature in the Milky Way with an N-body bar/bulge model. From the solar perspective, the distance distribution of particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between the two peaks qualitatively match the observed trends towards the Galactic bulge. We confirm clear signatures of cylindrical rotation in the pattern of mean radial velocity across the bar/bulge region. We also find possible imprints of coherent orbital motion inside the bar structure in the radial velocity distribution along l=0 degree, where the near and far sides of the bar/bulge show excesses of approaching and receding particles. The coherent orbital motion is also reflected in the slight displacement of the zero-velocity-line in the mean radial velocity, and the displacement of the maximum/minimum in the mean longitudinal proper motion across the bulge region. We find some degree of anisotropy in the stellar velocity within the X-shape, but the underlying orbital family of the X-shape cannot be clearly distinguished. Two potential applications of the X-shape in previous literature are tested, i.e., bulge rotation and Galactic center measurements. We find that the proper motion difference between the two sides of the X-shape can be used to estimate the mean azimuthal streaming motion of the bulge, but not the pattern speed of the bar. We also demonstrate that the Galactic center can be located with the X-shape, but the accuracy depends on the fitting scheme, the number of fields, and their latitudinal coverage., Comment: Minor changes to match the ApJ accepted version; 17 pages; emulateapj format. The electronic tables of our model result are available upon request

Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-dependent systematic errors in the derived absolute proper motions. The spectroscopically confirmed quasars were used to test the internal systematic and random error of the proper motions. The systematic errors of the overall proper motions in the SCUSS catalog are estimated as -0.08 and -0.06 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}, respectively. The random errors of the proper motions in the SCUSS catalog are estimated independently as 4.2 and 4.4 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}. There are no obvious position-, magnitude-, and color-dependent systematic errors of the SCUSS proper motions. The random error of the proper motions goes up with the magnitude from about 3 mas/yr at u < 18.0 mag to about 7 mas/yr at u = 22.0 mag. The proper motions of stars in SCUSS catalog are compared with those in the SDSS catalog, and they are highly consistent., Comment: 8 pages, 10 figures ,PASP accepted

We investigate properties of Galactic microlensing events in which a stellar object is lensed by a neutron star. For an all-sky photometric microlensing survey, we determine the number of lensing events caused by $\sim10^{5}$ potentially-observable radio pulsars to be $\sim0.2\ \rm{yr^{-1}}$ for $10^{10}$ background stellar sources. We expect a few detectable events per year for the same number of background sources from an astrometric microlensing survey. We show that such a study could lead to precise measurements of radio pulsar masses. For instance, if a pulsar distance could be constrained through radio observations, then its mass would be determined with a precision of $\sim10\%$. We also investigate the time-scale distributions for neutron star events, finding that they are much shorter than had been previously thought. For photometric events towards the Galactic centre that last $\sim15$ days, around $7\%$ will have a neutron star lens. This fraction drops rapidly for longer time-scales. Away from the bulge region we find that neutron stars will contribute $\sim40\%$ of the events that last less than $\sim10$ days. These results are in contrast to earlier work which found that the maximum fraction of neutron star events would occur on time-scales of hundreds of days., Comment: 10 pages, accepted for publication in ApJ. v2 updated to reflect change of title in proof stage

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., Comment: 5 pages; Contribution to proceedings for "The Milky Way Unravelled by Gaia" conference, Barcelona, Dec 2015

We describe a method to extract resonant orbits from N-body simulations exploiting the fact that they close in a frame rotating with a constant pattern speed. Our method is applied to the N-body simulation of the Milky Way by Shen et al. (2010). This simulation hosts a massive bar, which drives strong resonances and persistent angular momentum exchange. Resonant orbits are found throughout the disc, both close to the bar itself and out to the very edges of the disc. Using Fourier spectrograms, we demonstrate that the bar is driving kinematic substructure even in the very outer parts of the disc. We identify two major orbit families in the outskirts of the disc that make significant contributions to the kinematic landscape, namely the m:l = 3:-2 and 1:-1 families resonating with the pattern speed of the bar. A mechanism is described that produces bimodal distributions of Galactocentric radial velocities at selected azimuths in the outer disc. It occurs as a result of the temporal coherence of particles on the 3:-2 resonant orbits, which causes them to arrive simultaneously at pericentre or apocentre. This resonant clumping, due to the in-phase motion of the particles through their epicycle, leads to both inward and outward moving groups which belong to the same orbital family and consequently produce bimodal radial velocity distributions. This is a possible explanation of the bimodal velocity distributions observed towards the Galactic anti-Centre by Liu et al. (2012). Another consequence is that transient overdensities appear and dissipate (in a symmetric fashion) on timescales equal to the their epicyclic period resulting in a periodic pulsing of the disc's surface density., Comment: 11 Figures, 1 Table. Accepted for publication in ApJ. Version 2 reflects minor changes to the text. Animation referenced in Figure 7 is available at http://hubble.shao.ac.cn/~shen/resonantclumping/DensMovie.mp4

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., Comment: Accepted for publication in RAA; 12 pages, 6 figures (Fig. 1 at reduced resolution); full table of corrections available in online journal, with arxiv ancillary files (as ASCII table), or by request

We present a catalog of RR Lyrae stars (RRLs) observed by the Xuyi Schmidt Telescope Photometric Survey (XDSS). The area we consider is located in the North Galactic Cap, covering 376.75 sq deg at RA $\approx$ 150 deg and Dec $\approx$ 27 deg down to a magnitude limit of i $\approx$ 19. Using the variability information afforded by the multi-epoch nature of our XDSS data, combined with colors from the Sloan Digital Sky Survey, we are able to identify candidate RRLs. We find 318 candidates, derive distances to them and estimate the detection efficiency. The majority of our candidates have more than 12 observations and for these we are able to calculate periods. These also allows us to estimate our contamination level, which we predict is between 30% to 40%. Finally we use the sample to probe the halo density profile in the 9-49 kpc range and find that it can be well fitted by a double power law. We find good agreement between this model and the models derived for the South Galactic Cap using the Watkins et al. (2009) and Sesar et al. (2010) RRL data-sets, after accounting for possible contamination in our data-set from Sagittarius stream members. We consider non-spherical double power law models of the halo density profile and again find agreement with literature data-sets, although we have limited power to constrain the flattening due to our small survey area. Much tighter constraints will be placed by current and future wide-area surveys, most notably ESA's astrometric Gaia mission. Our analysis demonstrates that surveys with a limited number of epochs can effectively be mined for RRLs. Our complete sample is provided as accompanying online material., Comment: 14 pages, ApJ (in press)

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., Comment: 24 pages, 20 figures, submitted to ApJ

We investigate the extinction together with the radial velocity dispersion and distribution of red clump stars in the anti-center direction using spectra obtained with Hectospec on the MMT. We find that extinction peaks at Galactocentric radii of about 9.5 and 12.5 kpc, right in front of the locations of the Perseus and Outer arms and in line with the relative position of dust and stars in external spiral galaxies. The radial velocity dispersion peaks around 10kpc, which coincides with the location of the Perseus arm, yields an estimated arm-interarm density contrast of 1.3-1.5 and is in agreement with previous studies. Finally, we discover that the radial velocity distribution bifurcates around 10-11 kpc into two peaks at +27 km/s and -4 km/s. This seems to be naturally explained by the presence of the outer Lindblad resonance of the Galactic bar, but further observations will be needed to understand if the corotation resonance of the spirals arms also plays a role., Comment: 8 pages, 2 figures, accepted for publication in ApJL

The present-day state of the Milky Way disk can tell us much about the history of our Galaxy and provide insights into its formation. We have constructed a high-precision catalogue of disk stars using data from the Sloan Digital Sky Survey (SDSS) and use these stars to probe the heating history as well as investigating the detailed phase-space distribution. We also show how this sample can be used to probe the global properties of the Milky Way disk, employing the Jeans equations to provide a simple model of the potential close to the disk. Our model is in excellent agreement with others in the literature and provides an indication that the disk, rather than the halo, dominates the circular speed at the solar neighborhood. The work presented in these proceedings has been published as "Slicing and dicing the Milky Way disc in SDSS" (Smith et al. 2012)., Comment: Contribution to the Third Subaru Conference, Galactic Archaeology: Near Field Cosmology and the Formation of the Milky Way ed. Wako Aoki. See arXiv:1111.6920 for corresponding ApJ publication

We use the Stripe 82 proper motion catalogue of Bramich et al. (2008) to study the kinematics of Galactic disk stars in the solar neighborhood. We select samples of dwarf stars with reliable spectra and proper motions. They have cylindrical polar radius between 7 < R < 9 kpc, heights from the Galactic plane satisfying z < 2 kpc and span a range of metallicities -1.5 < [Fe/H] < 0. We develop a method for calculating and correcting for the halo contamination in our sample using the distribution of rotational velocities. Two Gaussians representing disk and halo populations are used to fit the radial and vertical velocity distributions via maximum likelihood methods. For the azimuthal velocities the same technique is used, except that a skewed non-Gaussian functional form now represents the disk velocity distribution. This enables us to compute the dispersions and cross-terms (the tilt and the vertex deviation) of the velocity ellipsoid as a function of height and metallicity. We also investigate the rotation lag of the disk, finding that the more metal-poor stars rotate significantly slower than the metal-rich stars. These samples provide important constraints on heating mechanisms in the Galactic disk and can be used for a variety of applications. We present one such application, employing the Jeans equations to provide a simple model of the potential close to the disk. Our model is in excellent agreement with others in the literature and provides an indication the disk, rather than the halo, dominates the circular speed at the solar neighborhood. We obtain a surface mass density within 1.1 kpc of around 66 M/pc^2 and estimate a local halo density of 0.015 M/pc^3 = 0.57 GeV/cm^3., Comment: 16 pages. ApJ, in press. Further minor typos fixed to match final proof version

This paper presents an investigation into the gravitational microlensing of quasars by stars and stellar remnants in the Milky Way. We present predictions for the all-sky microlensing optical depth, time-scale distributions and event rates for future large-area sky surveys. As expected, the total event rate increases rapidly with increasing magnitude limit, reflecting the fact that the number density of quasars is a steep function of magnitude. Surveys such as Pan-STARRS and LSST should be able to detect more than ten events per year, with typical event durations of around one month. Since microlensing of quasar sources suffers from fewer degeneracies than lensing of Milky Way sources, they could be used as a powerful tool for recovering the mass of the lensing object in a robust, often model-independent, manner. As a consequence, for a subset of these events it will be possible to directly `weigh' the star (or stellar remnant) that is causing the lensing signal, either through higher order microlensing effects and/or high-precision astrometric observations of the lens star (using, for example, Gaia or SIM-lite). This means that such events could play a crucial role in stellar astronomy. Given the current operational timelines for Pan-STARRS and LSST, by the end of the decade they could potentially detect up to 100 events. Although this is still too few events to place detailed constraints on Galactic models, consistency checks can be carried out and such samples could lead to exciting and unexpected discoveries., Comment: 11 pages, 8 figures. MNRAS (in press). Minor revisions according to referee's report; mainly presentational issues and clarification of a few items in the discussion; results and conclusions remain unchanged

Gravitational microlensing surveys target very dense stellar fields in the local group. As a consequence the microlensed source stars are often blended with nearby unresolved stars. The presence of `blending' is a cause of major uncertainty when determining the lensing properties of events towards the Galactic centre. After demonstrating empirical cases of blending we utilize Monte Carlo simulations to probe the effects of blending. We generate artificial microlensing events using an HST luminosity function convolved to typical ground-based seeing, adopting a range of values for the stellar density and seeing. We find that a significant fraction of bright events are blended, contrary to the oft-quoted assumption that bright events should be free from blending. We probe the effect that this erroneous assumption has on both the observed event timescale distribution and the optical depth, using realistic detection criteria relevent to the different surveys. Importantly, under this assumption the latter quantity appears to be reasonably unaffected across our adopted values for seeing and density. The timescale distribution is however biased towards smaller values, even for the least dense fields. The dominant source of blending is from lensing of faint source stars, rather than lensing of bright source stars blended with nearby fainter stars. We also explore other issues, such as the centroid motion of blended events and the phenomena of `negative' blending. Furthermore, we breifly note that blending can affect the determination of the centre of the red clump giant region from an observed luminosity function. This has implications for a variety of studies, e.g. mapping extinction towards the bulge and attempts to constrain the parameters of the Galactic bar through red clump giant number counts. (Abridged), Comment: 18 pages, 10 figures. MNRAS (in press)

Red clump giant stars can be used as distance indicators to trace the mass distribution of the Galactic bar. We use RCG stars from 44 bulge fields from the OGLE-II microlensing collaboration database to constrain analytic tri-axial models for the Galactic bar. We find the bar major axis is oriented at an angle of 24 - 27 degrees to the Sun-Galactic centre line-of-sight. The ratio of semi-major and semi-minor bar axis scale lengths in the Galactic plane x_0, y_0, and vertical bar scale length z_0, is x_0 : y_0 : z_0 = 10 : 3.5 : 2.6, suggesting a slightly more prolate bar structure than the working model of Gerhard (2002) which gives the scale length ratios as x_0 : y_0 : z_0 = 10 : 4 : 3 ., Comment: 15 pages, 10 figures, accepted for publication in MNRAS. Supplementary material available online: 10 pages, 10 figures

We describe the discovery of the longest microlensing event ever observed, OGLE-1999-BUL-32, also independently identified by the MACHO collaboration as MACHO-99-BLG-22. This unique event has an Einstein radius crossing time of 641 days. The high quality data obtained with difference image analysis shows a small but significant parallax signature. This parallax effect allows one to determine the Einstein radius projected onto the observer plane as rE^hat ~ 29.2AU. The transverse velocity projected onto the observer plane is about 79km/s. We argue that the lens is likely to be have a mass of at least a few solar masses, i.e., it could be a stellar black hole. The black hole hypothesis can be tested using the astrometric microlensing signature with the soon-to-be installed Advanced Camera for Surveys on board the Hubble Space Telescope. Deep X-ray and radio images may also be useful for revealing the nature of the object., Comment: submitted to Monthly Notices

We show that the perturbation series for quarkonia energies diverges at large orders. This results in a perturbative ambiguity in the energy that scales as e^(-1/a*Lambda) where a is the Bohr radius of quarkonium and Lambda is the QCD scale parameter. This ambiguity is associated with a nonperturbative contribution to the energy from distances of order 1/Lambda and greater. This contribution is separate from that of the gluon condensate., Comment: 6 pages, 2 figures

The top quark decays more quickly than the strong-interaction time scale, $\lqcd^{-1}$, and might be expected to escape the effects of nonperturbative QCD. Nevertheless, the top-quark pole mass, like the mass of a stable heavy quark, is ambiguous by an amount proportional to $\lqcd$., Comment: 9 pages, LaTeX

We calculate the O(alpha_s) and O(alpha_W m_t^2/M_W^2) corrections to the production of a single top quark via the weak process q qbar -> t bbar at the Fermilab Tevatron and the CERN Large Hadron Collider. An accurate calculation of the cross section is necessary in order to extract |V_tb| from experiment., Comment: LaTeX, 13 pages, replaced with version to appear in Phys. Rev. D