We use K-giant stars selected from the LAMOST DR5 to study the variation of the rotational velocity of the galactic halo at different space positions. Modelling the rotational velocity distribution with both the halo and disk components, we find that the rotational velocity of the halo population decreases almost linearly with increasing vertical distance to the galactic disk plane, $Z$, at fixed galactocentric radius, $R$. The samples are separated into two parts with $6<R<12$ kpc and $12<R<20$ kpc. We derive that the decreasing rates along $Z$ for the two subsamples are $-3.07\pm0.63$ and $-1.89\pm0.37$ km s$^{-1}$ kpc$^{-1}$, respectively. Compared with the TNG simulations, we suggest that this trend is probably caused by the interaction between the disk and halo. The results from the simulations show that only the oblate halo can provide a decreasing rotational velocity with an increasing $Z$. This indicates that the Galactic halo is oblate with galactocentric radius $R<20$ kpc. On the other hand, the flaring of the disk component (mainly the thick disk) is clearly traced by this study, with $R$ between 12 and 20 kpc, the disk can vertically extend to $6\sim10$ kpc above the disk plane. What is more interesting is that, we find the Gaia-Enceladus-Sausage (GES) component has a significant contribution only in the halo with $R<12$ kpc, i.e. a fraction of 23$-$47\%. While in the outer subsample, the contribution is too low to be well constrained.

We apply the vertical Jeans equation to the kinematics of Milky Way stars in the solar neighbourhood to measure the local dark matter density. More than 90,000 G- and K-type dwarf stars are selected from the cross-matched sample of LAMOST DR5 and Gaia DR2 for our analyses. The mass models applied consist of a single exponential stellar disc, a razor thin gas disc and a constant dark matter density. We first consider the simplified vertical Jeans equation which ignores the tilt term and assumes a flat rotation curve. Under a Gaussian prior on the total stellar surface density, the local dark matter density inferred from Markov Chain Monte Carlo simulations is $0.0133_{-0.0022}^{+0.0024}\ {\rm M}_{\odot}\,{\rm pc}^{-3}$. The local dark matter densities for subsamples in an azimuthal angle range of $-10^{\circ} < \phi < 5^{\circ}$ are consistent within their 1$\sigma$ errors. However, the northern and southern subsamples show a large discrepancy due to plateaux in the northern and southern vertical velocity dispersion profiles. These plateaux may be the cause of the different estimates of the dark matter density between the north and south. Taking the tilt term into account has little effect on the parameter estimations and does not explain the north and south asymmetry. Taking half of the difference of $\sigma_{z}$ profiles as unknown systematic errors, we then obtain consistent measurements for the northern and southern subsamples. We discuss the influence of the vertical data range, the scale height of the tracer population, the vertical distribution of stars and the sample size on the uncertainty of the determination of the local dark matter density.

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

We analyze the anisotropy profile of the Milky Way's smooth, diffuse stellar halo using SDSS/SEGUE blue horizontal branch stars and SDSS/SEGUE and LAMOST K giants. These intrinsically luminous stars allow us to probe the halo to approximately 100 kpc from the Galactic center. Line-of-sight velocities, distances, metallicities, and proper motions are available for all stars via SDSS/SEGUE, LAMOST, and Gaia, and we use these data to construct a full 7D set consisting of positions, space motions, and metallicity. We remove substructure from our samples using integrals of motion based on the method of Xue et al. We find radially dominated kinematic profiles with nearly constant anisotropy within 20 kpc, beyond which the anisotropy profile gently declines although remains radially dominated to the furthest extents of our sample. Independent of star type or substructure removal, the anisotropy depends on metallicity, such that the orbits of the stars become less radial with decreasing metallicity. For $-1.7<$ [Fe/H] $<-1$, the smooth, diffuse halo anisotropy profile begins to decline at Galactocentric distances $\sim20$ kpc, from $\beta\sim0.9$ to 0.7 for K giants and from $\beta\sim0.8$ to 0.1 for blue horizontal branch stars. For [Fe/H] $<-1.7$, the smooth, diffuse halo anisotropy remains constant along all distances with $0.2<\beta<0.7$ depending on the metallicity range probed, although independent on star type. These samples are ideal for estimating the total Galactic mass as they represent the virialized stellar halo system.

The effective vacuum energy density contributed by the non-trivial contortion distribution and the bare vacuum energy density can be viewed as the energy density of the auxiliary quintessence field potential. We find that the negative bare vacuum energy density from string landscape leads to a monotonically decreasing quintessence potential while the positive one from swampland leads to the meta stable or stable de Sitter like potential. Moreover, the non-trivial Brans-Dicke like coupling between quintessence field and gravitation field is necessary in the latter case.

A large number of new members ($\sim$150) of the Cetus Stream (CS) were identified from their clustering features in dynamical space using 6D kinematic data by combining LAMOST DR5 and Gaia DR2 surveys. They map a diffuse structure that extends over at least 100 degrees in the northern and southern Galactic hemispheres, at heliocentric distances between 20 to 50 kpc. Taking advantage of this expanded dataset, we model the stream with a suite of tailored N-body simulations. Our findings exclude the possibility that the NGC 5824 globular cluster is the core of the progenitor of the stream, as postulated by previous studies. Our best models, which successfully reproduce the features of the CS indicate that the progenitor is likely a dwarf galaxy of $\sim$ 2$\times$10$^9$M$_{\odot}$, with a diffuse disc morphology. The merger occured $\sim$ 5 Gyr ago and since then it has experienced approximately eight apo-center passages. Our results suggest that NGC 5824 was either a globular cluster situated off-centre in the dwarf progenitor or, alternatively, it was the nuclear star cluster of another dwarf galaxy that has very similar orbit as the progenitor of the CS. In both scenarios, the progenitor systems would leave streams around NGC 5824, but with distinct distance distributions. To discriminate between these scenarios, the detection and accurate distance measurements of the predicted stream around the GC are crucial, which will be possible in the upcoming LSST era. Our simulations also predict that part of the Southern Cetus stream is very likely the newly discovered Palca stream, and possibly related to another, more diffuse Southern substructure, the Eridanus-Pheonix overdensity.

With high spatial resolution, polarimetric imaging of a supermassive black hole, like M87$^\star$ or Sgr A$^\star$, by the Event Horizon Telescope can be used to probe the existence of ultralight bosonic particles, such as axions. Such particles can accumulate around a rotating black hole through the superradiance mechanism, forming an axion cloud. When linearly polarized photons are emitted from an accretion disk near the horizon, their position angles oscillate due to the birefringent effect when traveling through the axion background. In particular, the observations of supermassive black holes M87$^\star$ (Sgr A$^\star$) can probe the dimensionless axion-photon coupling $c = 2 \pi g_{a \gamma} f_a$ for axions with mass around $O(10^{-20})$~eV ($O( 10^{-17}$)~eV) and decay constant $f_a < O(10^{16})$ GeV, which is complimentary to other axion measurements.

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.

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

We use a novel cluster identification tool StarGO to explore the metal poor ([Fe/H] $<$ -1.5) outer stellar halo (d $>$ 15 kpc) of the Milky Way using data from Gaia, LAMOST and SDSS. Our method is built using an unsupervised learning algorithm, a self-organizing map, which trains a 2-D neural network to learn the topological structures of a data set from an n-D input space. Using a 4-D space of angular momentum and orbital energy, we identify three distinct groups corresponding to the Sagittarius, Orphan, and Cetus Streams. For the first time we are able to discover a northern counterpart to the Cetus stream. We test the robustness of this new detection using mock data and find that the significance is more than 5-sigma. We also find that the existing southern counterpart bifurcates into two clumps with different radial velocities. By exploiting the visualization power of StarGO, we attach MW globular clusters to the same trained neural network. The Sagittarius stream is found to have five related clusters, confirming recent literature studies, and the Cetus stream has one associated cluster, NGC 5824. This latter association has previously been postulated, but can only now be truly confirmed thanks to the high-precision Gaia proper motions and large numbers of stellar spectra from LAMOST. The large metallicity dispersion of the stream indicates that the progenitor cannot be a globular cluster. Given the mean metallicity of the stream, we propose that the stream is the result of a merger of a low-mass dwarf galaxy that hosted a large nuclear star cluster (NGC 5824).

We construct a large halo K-giant sample by combining the positions, distances, radial velocities, and metallicities of over 13,000 LAMOST DR5 halo K giants with the Gaia DR2 proper motions, which covers a Galactocentric distance range of 5-120 kpc. Using a position-velocity clustering estimator (the 6Distance), we statistically quantify the presence of position-velocity substructure at high significance: K giants have more close pairs in position-velocity space than a smooth stellar halo. We find that the amount of substructure in the halo increases with increasing distance and metallicity. With a percolation algorithm named friends-of-friends (FoF) to identify groups, we identify members belonging to Sagittarius (Sgr) Streams, Monoceros Ring, Virgo overdensity, Hercules-Aquila Cloud, Orphan Streams and other unknown substructures and find that the Sgr streams account for a large part of grouped stars beyond 20 kpc and enhance the increase of substructure with distance and metallicity. For the first time, we identify spectroscopic members of Monoceros Ring in the south and north Galactic hemisphere, which presents a rotation of about 185 km s^{-1} and mean metallicity is -0.66 dex.

A class of dark photon dark matter models with ultralight masses would lead to oscillation of a test body through a coupling with baryons or $B-L$ charge. This periodical oscillation of an observer results in swing of a star's apparent position due to the effect of aberration of light, which could be probed with high-precision astrometry observations of stars in the Milky Way. We propose to use the observations of stellar positions of a number of stars by Gaia to search for this kind of dark photon dark matter. We show that this astrometry method is able to give promising sensitivities to search for the dark photon dark matter in the mass range of $10^{-23}\sim10^{-21}$ eV.

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.

Based on the [Fe/H] versus [Mg/Fe] diagram and distances from APOGEE data release 14, we compare the spatial distributions, the l-Vlos diagram and the abundance gradients between high-[Mg/Fe] and low-[Mg/Fe] sequences. The two sequences are clearly shown at 5<|Z|<10 kpc in the metallicity range of -1.6 <[Fe/H] <-0.7, where the halo at |Z| > 10 kpc consists of low-[Mg/Fe] stars only. In the intermediate-metallicity range of -1.1 <[Fe/H]<-0.7, a [Mg/Fe] gradient is detected for stars at |Z|=10-30 kpc and it flattens out at |Z|>30 kpc. The l-Vlos diagram is adopted to separate halo stars from the disk by defining the transition metallicity, which is of [Fe/H]~ -1.1 dex for the high-[Mg/Fe] sequence and of [Fe/H]~-0.7 dex for the low-[Mg/Fe] sequence. The R and |Z| distributions for the high-[Mg/Fe] sequence, the thick disk at -1.1<[Fe/H]<-0.7 and the in situ halo at -1.6<[Fe/H]<-1.1, have a cutoff at R~15 kpc and |Z|~10 kpc, beyond which low-[Mg/Fe] halo stars are the main contributions. In the metallicity range of -1.6<[Fe/H]<-0.7, there is a negative metallicity gradient for the high-[Mg/Fe] halo at |Z|<8-10 kpc, while only a marginal or no slope in the [Fe/H] versus |Z| diagram for the low-[Mg/Fe] halo at |Z|<8-10 kpc, beyond which both the high-[Mg/Fe] halo and low-[Mg/Fe] halo flatten out toward |Z| > 20 kpc. These results indicate a complicated formation history of the Galaxy and we may see a hint of a three-section halo, i.e. the inner in situ halo within $|Z|~8-10$ kpc, the intermediately outer dual-mode halo at |Z|~10-30 kpc, and the extremely outer accreted halo with |Z|>30 kpc.

With the help of Gaia DR2, we are able to obtain the full 6-D phase space information for stars from LAMOST DR5. With high precision of position, velocity, and metallicity, the rotation of the local stellar halo is presented using the K giant stars with [Fe/H]$<-1$ dex within 4 kpc from the Sun. By fitting the rotational velocity distribution with three Gaussian components, stellar halo, disk, and a counter-rotating hot population, we find that the local halo progradely rotates with $V_T=+27^{+4}_{-5}$ km s$^{-1}$ providing the local standard of rest velocity of $V_{LSR}=232$ km s$^{-1}$. Meanwhile, we obtain the dispersion of rotational velocity is $\sigma_{T}=72^{+4}_{-4}$ km s$^{-1}$. Although the rotational velocity strongly depends on the choice of $V_{LSR}$, the trend of prograde rotation is substantial even when $V_{LSR}$ is set at as low as 220 km s$^{-1}$. Moreover, we derive the rotation for subsamples with different metallicities and find that the rotational velocity is essentially not correlated with [Fe/H]. This may hint a secular evolution origin of the prograde rotation. It shows that the metallicity of the progradely rotating halo is peaked within -1.9$<$[Fe/H]$<$-1.6 without considering the selection effect. We also find a small fraction of counter-rotating stars with larger dispersion and lower metallicity. Finally, the disk component rotates with $V_T=+182^{+6}_{-6}$ km s$^{-1}$ and $\sigma_T=45^{+3}_{-3}$ km s$^{-1}$, which is quite consistent with the metal-weak thick disk population.

The accelerating expansion of universe can be described by the non-zero cosmological constant or the dark energy. However, the origin of the dark energy remains a mystery of modern physics. The local Lorentz invariance is the most exact symmetry of the Nature on the one hand, but all quantum gravity theories predict Lorentz violation on the other hand. The local Lorentz violation induced by the quantum gravity at the very early universe may be transformed into large scale by inflation. Combining the low-$l$ anomalies of the CMB spectrum, we propose that the local Lorentz invariance may be broken at the large scale. We construct the effective gravity at the cosmic scale with a local $SO(3)$ symmetry. The theory exhibits non-trivial contortion distribution even with scalar matter source. The FRW like solution of the theory is analyzed and the contortion distribution contributes a dark energy like effect which is responsible for the accelerating expansion of the universe. It reveals that the dark energy may be the remnants of quantum gravity in this sense.

We study the radial metallicity gradient $\Delta[M/H]/\Delta R_g$ as a function of [Mg/Fe] and $|Z|$ with the help of a guiding radius based on the Apache Point Observatory Galactic Evolution Experiment and Gaia and then analyze the radial migration effect on the radial metallicity gradient and metallicity-rotation gradient between the Galactic thin and thick disks. The derived trend of gradient $\Delta[M/H]/\Delta R_g$ versus [Mg/Fe] shows a transition at [Mg/Fe] $\sim 0.18$ dex, below which the gradient is negative and varies a little as [Mg/Fe] increases; however, it changes sharply in [Mg/Fe] ranges of 0.16-0.18, above which the gradient increases linearly with increasing [Mg/Fe], being a positive value at [Mg/Fe]$\gtrsim 0.22$ dex. These positive gradients in the high-[Mg/Fe] populations are found at $|Z| < 0.8$ kpc, and there are nearly no gradients toward higher $|Z|$. By comparing the metallicity distributions, the radial metallicity gradients $\Delta[M/H]/\Delta R$ and the metallicity-rotation gradients between the total sample and $|R-R_g|<2$ kpc subsample (or $|R-R_g|>2$ kpc subsample), we find that, for the thick disk, blurring flattens the gradient $\Delta[M/H]/\Delta R$ and favors metal-poor high-eccentricity stars. These stars are responsible for the measured positive metallicity-rotation gradient of the thick disk.

The anisotropy parameter $\beta$ characterizes the extent to which orbits in stellar systems are predominantly radial or tangential, and is likely to constrain, for the stellar halo of the Milky Way, scenarios for its formation and evolution. We have measured the anisotropy $\beta$ as a function of Galactocentric radius from $5-100$ kpc for over 8600 metal poor ([Fe/H] $<-1.3$) halo K giants from the LAMOST catalog with line-of-sight velocities and distances, matched to proper motions from the second $Gaia$ data release. We construct full 6-D positions and velocities for the K giants to directly measure the 3 components of the velocity dispersion $(\sigma_r, \sigma_\theta, \sigma_\phi)$ (in spherical coordinates). We find that the orbits in the halo are radial over our full Galactocentric distance range reaching over 100 kpc. The anisotropy remains remarkably unchanged with Galactocentric radius from approximately 5 to 25 kpc, with an amplitude that depends on the metallicity of the stars, dropping from $\beta \approx 0.9$ for $-1.8 \leq$ [Fe/H] $< -1.3$ (for the bulk of the stars) to $\beta \approx 0.6$ for the lowest metallicities ([Fe/H] $< -1.8$). Considering our sample as a whole, $\beta\approx0.8$ and, beyond 25 kpc, the orbits gradually become less radial and anisotropy decreases to $\beta<0.3$ past 100 kpc. Within 8 kpc, $\beta<0.8$. The measurement of anisotropy is affected by substructure and streams, particularly beyond a Galactocentric distance of approximately 25 kpc, where the Sagittarius stream is prominent in the data. These results are complimentary to recent analysis of simulations by Loebman et al. and of SDSS/$Gaia$ DR1 data by Belokurov et al.

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.

We statistically quantify the amount of substructure in the Milky Way stellar halo using a sample of 4568 halo K giant stars at Galactocentric distances ranging over 5-125 kpc. These stars have been selected photometrically and confirmed spectroscopically as K giants from the Sloan Digital Sky Survey's SEGUE project. Using a position-velocity clustering estimator (the 4distance) and a model of a smooth stellar halo, we quantify the amount of substructure in the halo, divided by distance and metallicity. Overall, we find that the halo as a whole is highly structured. We also confirm earlier work using BHB stars which showed that there is an increasing amount of substructure with increasing Galactocentric radius, and additionally find that the amount of substructure in the halo increases with increasing metallicity. Comparing to resampled BHB stars, we find that K giants and BHBs have similar amounts of substructure over equivalent ranges of Galactocentric radius. Using a friends-of-friends algorithm to identify members of individual groups, we find that a large fraction (~33%) of grouped stars are associated with Sgr, and identify stars belonging to other halo star streams: the Orphan Stream, the Cetus Polar Stream, and others, including previously unknown substructures. A large fraction of sample K giants (more than 50%) are not grouped into any substructure. We find also that the Sgr stream strongly dominates groups in the outer halo for all except the most metal-poor stars, and suggest that this is the source of the increase of substructure with Galactocentric radius and metallicity.

The dipole anomaly in the power spectrum of CMB may indicate that the Lorentz boost invarianc is violated at cosmic scale. We assume that the Lorentz symmetry is violated partly from the scale of galaxy. We employ the symmetry of very special relativity as an example to illustrate the Lorentz violation effect by constructing the corresponding gauge theories as the effective gravitational theory at the large scale. We find the common feather of these gravitation models is the non-triviality of spacetime torsion and contorsion even if the matter source is of only scalar matter. The presence of non-trivial contorsion contributes an effective enenrgy-momentum distribution which may account for part of dark matter effect.

We characterise the radial density, metallicity and flattening profile of the Milky Way's stellar halo, based on the large sample of 1757 spectroscopically confirmed giant stars from SDSS/SEGUE-2 after excising stars that were algorithmically attributed to apparent halo substructure (including the Sagittarius stream). Compared to BHB stars or RR Lyrae, giants are more readily understood tracers of the overall halo star population, with less bias in age or metallicity. The well-characterized selection function of the sample enables forward modelling of those data, based on ellipsoidal stellar density models, $\nu_* (R,z)$, with Einasto profiles and (broken) power laws for their radial dependence, combined with a model for the metallicity gradient and the flattening profile. Among models with constant flattening, these data are reasonably well fit by an Einasto profile of $n=3.1\pm 0.5$ with an effective radius $\rm r_{eff} = 15\pm2~$kpc and a flattening of $q=0.7\pm 0.02$; or comparably well by an equally flattened broken power-law, with radial slopes of $\alpha_{in}=2.1\pm 0.3$ and $\alpha_{out}=3.8\pm 0.1$, with a break-radius of $r_{break}=18\pm1$~kpc; this is largely consistent with earlier work. We find a modest, but significant metallicity gradient within the "outer" stellar halo, $\rm [Fe/H]$ decreasing outward. If we allow for a variable flattening $q = f(r_{GC} )$, we find the distribution of halo giants to be considerably more flattened at small radii, $q({\rm 10~kpc})\sim 0.57$, compared to $q(>30{\rm kpc})\sim 0.8$. Remarkably, the data are then very well fit by a single power-law of index $\rm \sim 4.2\pm0.1$ of the variable $r_q\equiv\sqrt{R^2+(z/q(r))^2}$. In this simple and better fitting model, there is a break in flattening at $\sim 20$~kpc, instead of a break in the radial density function.

We present an analysis of a new, large sample of field blue-straggler stars (BSSs) in the thick disk and halo system of the Galaxy, based on stellar spectra obtained during the Sloan Digital Sky Survey (SDSS) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE). Using estimates of stellar atmospheric parameters obtained from application of the SEGUE Stellar Parameter Pipeline, we obtain a sample of some 8000 BSSs, which are considered along with a previously selected sample of some 4800 blue horizontal-branch (BHB) stars. We derive the ratio of BSSs to BHB stars, F$_{\rm BSS/BHB}$, as a function of Galactocentric distance and distance from the Galactic plane. The maximum value found for F$_{\rm BSS/BHB}$ is $\sim~$4.0 in the thick disk (at 3 kpc $<$ $|$Z$|$ $<$ 4 kpc), declining to on the order of $\sim~1.5-2.0$ in the inner-halo region; this ratio continues to decline to $\sim~$1.0 in the outer-halo region. We associate a minority of field BSSs with a likely extragalactic origin; at least 5$\%$ of the BSS sample exhibit radial velocities, positions, and distances commensurate with membership in the Sagittarius Stream.

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 present an online catalog of distance determinations for $\rm 6036$ K giants, most of which are members of the Milky Way's stellar halo. Their medium-resolution spectra from SDSS/SEGUE are used to derive metallicities and rough gravity estimates, along with radial velocities. Distance moduli are derived from a comparison of each star's apparent magnitude with the absolute magnitude of empirically calibrated color-luminosity fiducials, at the observed $(g-r)_0$ color and spectroscopic [Fe/H]. We employ a probabilistic approach that makes it straightforward to properly propagate the errors in metallicities, magnitudes, and colors into distance uncertainties. We also fold in ${\it prior}$ information about the giant-branch luminosity function and the different metallicity distributions of the SEGUE K-giant targeting sub-categories. We show that the metallicity prior plays a small role in the distance estimates, but that neglecting the luminosity prior could lead to a systematic distance modulus bias of up to 0.25 mag, compared to the case of using the luminosity prior. We find a median distance precision of $16\%$, with distance estimates most precise for the least metal-poor stars near the tip of the red-giant branch. The precision and accuracy of our distance estimates are validated with observations of globular and open clusters. The stars in our catalog are up to 125 kpc distant from the Galactic center, with 283 stars beyond 50 kpc, forming the largest available spectroscopic sample of distant tracers in the Galactic halo.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

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.

On the basis of a recently derived color-metallicity relation and stellar parameters from the Sloan Digtal Sky Survey Data Release 7 spectroscopic survey, a large sample of red horizontal branch candidates have been selected to serve as standard candles. The metallicity and kinematic distributions of these stars indicate that they mainly originate from the thick-disk and the halo populations. The typical thick disk is characterized by the first group peaking at [Fe/H] ~ -0.6, Vrot ~ 170 km/s with a vertical scale height around |Z| ~ 1.2 kpc, while stars with [Fe/H] < -0.9 are dominated by the halo population. Two sub-populations of the halo are suggested by the RHB stars peaking at [Fe/H] ~-1.3: one component with Vrot > 0 km/s (Halo I) shows a sign of metallicity gradient in the [Fe/H] versus |Z| diagram, while the other with Vrot < 0 km/s (Halo II) does not. The Halo I mainly clumps at the inner halo with R < 10 kpc and the Halo II comes both from the inner halo with R < 10 kpc and the outer halo with R > 10 kpc based on the star distribution in the R versus |Z| diagram.

We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Way's halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative "close pair distribution" (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $\rm r_{gc} < 20$ kpc.

We use a sample of blue horizontal branch (BHB) stars from the Sloan Digital Sky Survey Data Release 7 to explore the structure of the tidal tails from the Sagittarius Dwarf Galaxy. We use a method yielding BHB star candidates with up to ~70% purity from photometry alone. The resulting sample has a distance precision of roughly 5% and can probe distances in excess of 100 kpc. Using this sample, we identify a possible extension to the trailing arm at distances of 60-80 kpc from the Sun with an estimated significance of at least 3.8 sigma. Current models predict that a distant `returning' segment of the debris stream should exist, but place it substantially closer to the Sun where no debris is observed in our data. Exploiting the distance precision of our tracers, we estimate the mean line-of-sight thickness of the leading arm to be ~3 kpc, and show that the two `bifurcated' branches of the debris stream differ by only 1-2 kpc in distance. With a spectroscopic very pure BHB star subsample, we estimate the velocity dispersion in the leading arm, 37 km s^-1, which is in reasonable agreement with models of Sgr disruption. We finally present a sample of high-probability Sgr BHB stars in the leading arm of Sgr, selected to have distances and velocities consistent with Sgr membership, to allow further study.

If the stellar halos of disk galaxies are built up from the disruption of dwarf galaxies, models predict highly structured variations in the stellar populations within these halos. We test this prediction by studying the ratio of blue horizontal branch stars (BHB stars; more abundant in old, metal-poor populations) to main-sequence turn-off stars (MSTO stars; a feature of all populations) in the stellar halo of the Milky Way using data from the Sloan Digital Sky Survey. We develop and apply an improved technique to select BHB stars using ugr color information alone, yielding a sample of ~9000 g<18 candidates where ~70% of them are BHB stars. We map the BHB/MSTO ratio across ~1/4 of the sky at the distance resolution permitted by the absolute magnitude distribution of MSTO stars. We find large variations of BHB/MSTO star ratio in the stellar halo. Previously identified, stream-like halo structures have distinctive BHB/MSTO ratios, indicating different ages/metallicities. Some halo features, e.g., the low-latitude structure, appear to be almost completely devoid of BHB stars, whereas other structures appear to be rich in BHB stars. The Sagittarius tidal stream shows an apparent variation in BHB/MSTO ratio along its extent, which we interpret in terms of population gradients within the progenitor dwarf galaxy. Our detection of coherent stellar population variations between different stellar halo substructures provides yet more support to cosmologically motivated models for stellar halo growth.

We investigate the performance of some common machine learning techniques in identifying BHB stars from photometric data. To train the machine learning algorithms, we use previously published spectroscopic identifications of BHB stars from SDSS data. We investigate the performance of three different techniques, namely k nearest neighbour classification, kernel density estimation and a support vector machine (SVM). We discuss the performance of the methods in terms of both completeness and contamination. We discuss the prospect of trading off these values, achieving lower contamination at the expense of lower completeness, by adjusting probability thresholds for the classification. We also discuss the role of prior probabilities in the classification performance, and we assess via simulations the reliability of the dataset used for training. Overall it seems that no-prior gives the best completeness, but adopting a prior lowers the contamination. We find that the SVM generally delivers the lowest contamination for a given level of completeness, and so is our method of choice. Finally, we classify a large sample of SDSS DR7 photometry using the SVM trained on the spectroscopic sample. We identify 27,074 probable BHB stars out of a sample of 294,652 stars. We derive photometric parallaxes and demonstrate that our results are reasonable by comparing to known distances for a selection of globular clusters. We attach our classifications, including probabilities, as an electronic table, so that they can be used either directly as a BHB star catalogue, or as priors to a spectroscopic or other classification method. We also provide our final models so that they can be directly applied to new data.

We present some of the strategies being developed to classify and parameterize objects obtained with spectra from the Sloan Digital Sky Survey (SDSS) and the RAdial Velocity Experiment (RAVE) and present some results. We estimate stellar atmospheric parameters (effective temperature, gravity, and metallicity) from spectral and photometric data and use these to analyse Galactic populations. We demonstrate this through the selection of a sample of candidate Blue Horizontal-Branch and RR Lyrae stars selected from SDSS/SEGUE.

We derive new constraints on the mass of the Milky Way's dark matter halo, based on a set of halo stars from SDSS as kinematic tracers. Our sample comprises 2401 rigorously selected Blue Horizontal-Branch (BHB) halo stars drawn from SDSS DR-6. To interpret these distributions, we compare them to matched mock observations drawn from two different cosmological galaxy formation simulations designed to resemble the Milky Way, which we presume to have an appropriate orbital distribution of halo stars. We then determine which value of $\rm V_{cir}(r)$ brings the observed distribution into agreement with the corresponding distributions from the simulations. This procedure results in an estimate of the Milky Way's circular velocity curve to $\sim 60$ kpc, which is found to be slightly falling from the adopted value of $\rm 220 km s^{-1}$ at the Sun's location, and implies M$(<60 \rm kpc) = 4.0\pm 0.7\times 10^{11}$M$_\odot$. The radial dependence of $\rm V_{cir}(r)$, derived in statistically independent bins, is found to be consistent with the expectations from an NFW dark matter halo with the established stellar mass components at its center. If we assume an NFW halo profile of characteristic concentration holds, we can use the observations to estimate the virial mass of the Milky Way's dark matter halo, M$_{\rm vir}=1.0^{+0.3}_{-0.2} \times 10^{12}$M$_\odot$, which is lower than many previous estimates. This estimate implies that nearly 40% of the baryons within the virial radius of the Milky Way's dark matter halo reside in the stellar components of our Galaxy. A value for M$_{\rm vir}$ of only $\sim 1\times10^{12}$M$_\odot$ also (re-)opens the question of whether all of the Milky Way's satellite galaxies are on bound orbits.