We present {complexity} of the Galactic outer disc {by fitting the stellar volume densities of the red giant branch stars with a two-disc component model}. {The discs are confirmed to} extend to $R\sim19$\,kpc. The radial density profile of the discs shows {two breaks at $R\sim11$ and $\sim14$\,kpc, respectively, which} separate the radial profile into three segments with different scale lengths of $2.12\pm0.26$, $1.18\pm0.08$, and $2.72$\,kpc at $R<11$, $11\leq R\leq14$, and $R>14$\,kpc, respectively. The first {break} is likely due to the sudden drop in the radial profile of the thin disc, which may be an evidence of the radial migration. {Beyond $14$\,kpc, the thick disc becomes prominent and the transition from thin to thick disc leads to the second break.} This implies that the geometrically defined thick disc is more {radially} extended than the thin disc. This is also supported by the larger scale length of the thick disc than that of the thin disc. Meanwhile, {the scale height of the thicker component increases from $0.637_{-0.036}^{+0.056}$ at $R=8$ to $1.284_{-0.079}^{+0.086}$\,kpc at $R=19$\,kpc, showing an intensive flared disc}. Moreover, rich substructures are displayed in the residuals of the stellar density. Among them, the substructures $D14+2.0$ and $O14-1.5$ show a north-south asymmetry, which can be essentially explained by southward shifting of the thick disc. However, no significant overdensity is found for the Monoceros ring. Finally, the thick disc shows a ripple-like feature with unclear origin at $9<R<10.5$\,kpc.

From the derived stellar density profile using LAMOST giant stars, we find that the Galactic disk does not show truncation or break, but smoothly transit to the halo from 19 kpc. The scale length of the outer disk is only $1.6\pm0.1$\,kpc, substantially smaller than previous results. This implies that the shapes of the inner and outer disk are different. Meanwhile, the disk flaring is not only found in older populations, but also in younger population. Moreover, the vertical oscillations of the disk are identified in a wide range or $R$ from 8 to 14 kpc. We also find that the velocity dispersion profile as a function of the Galactocentric radius is flat with scale length of $26.3\pm3.2$\,kpc. We confirm that the radial velocity profile in outer disk is significantly affected by asymmetric motion. The bar with either a slower or a faster pattern speed can induce the similar radial asymmetric motion.

We present the peculiar in-plane velocities derived from the LAMOST red clump stars, which are purified and separated by a novel approach into two groups with different ages. The samples are mostly contributed around the Galactic anti-centre direction such that we are able to map the radial profiles of the radial and azimuthal velocities in the outer disc. From the variations of the in-plane velocities with the Galactocentric radius for the younger and older populations, we find that both radial and azimuthal velocities are not axisymmetric at $8<R<14\,kpc$. The two red clump populations show that the mean radial velocity is negative within $R\sim9\,kpc$ and positive beyond. This is likely because of the perturbation induced by the rotating bar. The cross-zero radius, $R\sim9$\, kpc, essentially indicates the rough location of the outer Lindblad resonance (OLR) radius. Given the circular speed of 238\,km$\rm s^{-1}$, then the pattern speed of the bar can be approximated as $45$\,km$\rm s^{-1}\rm kpc^{-1}$. The young red clump stars show larger mean radial velocity than the old population by about 3$\,km\rm s^{-1}$ between $R\sim9$ and 12\,kpc. This is possibly because the younger population is more sensitive to the perturbation than the older one. The radial profiles of the mean azimuthal velocity for the two populations show an interesting U-shape, i.e. at $R<10.5\,kpc$, the azimuthal velocity declines with $R$ by about 10$\,km\rm s^{-1}$, while at $R>10.5$ it increases with $R$ to 240-245$\,km\rm s^{-1}$. It is not clear why the mean azimuthal velocity shows the U-shape along the Galactic anti-centre direction. Meanwhile, the azimuthal velocity for the younger population is slightly larger than the older one and the difference moderately declines with $R$. Beyond $R\sim12\,kpc$, the azimuthal velocities for the two populations are indistinguishable.

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

We present stellar density maps of the Galactic outer disc with red clump stars from the LAMOST data. These samples are separated into younger (mean age ~ 2.7 Gyr) and older (mean age ~ 4.6 Gyr) populations so that they can trace the variation of the structures with ages in the range of the Galactocentric radius R from 9 to 13.5 kpc. We show that both the scale heights for the two populations increase with R and display radial gradients of 48 +/- 6 and 40 +/- 4 pc/kpc for the older and younger populations, respectively. This is evident that the flaring occurs in the thin disc populations with a wide range of ages. Moreover, the intensity of flaring seems not significantly related to the age of the thin disc populations. On the other hand, the scale lengths of the radial surface density profiles are 4.7 +/- 0.5 kpc for the younger and 3.4 +/- 0.2 kpc for the older population, meaning that the younger disc population is more radially extended than the older one. Although the fraction of the younger population mildly increases from 28% at R ~ 9 to about 35% at R ~ 13 kpc, the older population is prominent with the fraction no less than 65% in the outer disc.

We derive the mean velocity components at various Galactocentric radii from 8 to 14 kpc using about 40,000 red clump stars observed in the LAMOST survey. We find that the vertical bulk motion for younger red clump stars are significantly larger than that for the older red clump stars. This is likely the kinematical feature of the Galactic warp around its line-of-node, which is located close to the Galactic anti-center region. It is evident that the warp are mainly contributed by the younger stars rather than the older stars. The age variation in the vertical kinematics favors a formation scenario where the Galactic warp is originated from infalling misaligned gas.

We present the peculiar in-plane velocities derived from the LAMOST red clump stars. From the variations of the in-plane velocity with the Galactocentric radius for the young and old red clump stars, we are able to identify two types of peculiar velocities: 1) both the two red clump populations show that the radial velocity is negative within $R=9.0$\,kpc and becomes positive beyond (denoted as the \emph{long-wave} mode); and 2) the young red clump stars show larger mean radial velocity than the old population by about 3\,km$\rm s^{-1}$ between $R=9$ and 12\,kpc (denoted as the \emph{short-wave} mode). We find that the elliptical disk induced by the rotating bar can well explain the \emph{long-wave} mode peculiar velocity. The axis ratio of the elliptical disk is around 0.8-0.95 and the disk keeps circular at $R=9.24\pm0.2$\,kpc, which should be the location of the outer Lindblad resonance radius (OLR). Adopting the circular speed of 238\,km$\rm s^{-1}$, the pattern speed of the bar is then derived as $48\pm3$\,km$\rm s^{-1}$kpc$^{-1}$ from the location of OLR. On the other hand, the \emph{short-wave} mode is likely the perturbation of the spiral arms as density waves.

We estimate the age for the individual stars located at the lower part of the red giant branch from the LAMOST DR2 K giant sample. Taking into account the selection effects and the volume completeness, the age--metallicity map for the stars located between 0.3 and 1.5 kpc from the Sun is obtained. A significant substructure (denoted as the \it{narrow stripe}) located from (age, [Fe/H])$\sim$(5, 0.4) to (10 Gyr, -0.4 dex) in the age--metallicity map is clearly identified. Moreover, the \it{narrow stripe} stars are found the dominate contributors to several velocity substructures, including the well-known Hercules stream. The substantially large difference between the observed guiding-center radii and the birth radii inferred from the age--metallicity relation is evident that the \it{narrow stripe} stars have been radially migrated from about R$\sim4$ kpc to the solar neighborhood. This implies that the Hercules stream may not be owe to the resonance associated with the bar, but may be the kinematic imprint of the inner disk and later moved out due to radial migration. We estimate that the traveling speed of the radial migration are roughly 1.1$\pm0.1$ kpc Gyr$^{-1}$, equivalent with about $1.1\pm0.1$ km s$^{-1}$. This is in agreement with the median $v_R$ of $2.6^{+1.8}_{-1.9}$ km s$^{-1}$ of the \it{narrow stripe}. We also obtain that about one third stars in the solar neighborhood are radially migrated from around 4 kpc. Finally, we find that the radial migration does not lead to additional disk thickening according to the distribution of $z_{max}$.

We present a sample of about 120,000 red clump candidates selected from the LAMOST DR2 catalog based on the empirical distribution model in the effective temperature vs. surface gravity plane. Although, in general, red clump stars are considered as the standard candle, they do not exactly stay in a narrow range of absolute magnitude, but may extend to more than 1 magnitude depending on their initial mass. Consequently, conventional oversimplified distance estimations with assumption of fixed luminosity may lead to systematic bias related to the initial mass or the age, which may potentially affect the study of the evolution of the Galaxy with red clump stars. We therefore employ an isochrone-based method to estimate the absolute magnitude of red clump stars from their observed surface gravities, effective temperatures, and metallicities. We verify that the estimation well removes the systematics and provide an initial mass/age independent distance estimates with accuracy less than 10%.

In this work, we select the high signal-to-noise ratio spectra of stars from the LAMOST data andmap theirMK classes to the spectral features. The equivalentwidths of the prominent spectral lines, playing the similar role as the multi-color photometry, form a clean stellar locus well ordered by MK classes. The advantage of the stellar locus in line indices is that it gives a natural and continuous classification of stars consistent with either the broadly used MK classes or the stellar astrophysical parameters. We also employ a SVM-based classification algorithm to assignMK classes to the LAMOST stellar spectra. We find that the completenesses of the classification are up to 90% for A and G type stars, while it is down to about 50% for OB and K type stars. About 40% of the OB and K type stars are mis-classified as A and G type stars, respectively. This is likely owe to the difference of the spectral features between the late B type and early A type stars or between the late G and early K type stars are very weak. The relative poor performance of the automatic MK classification with SVM suggests that the directly use of the line indices to classify stars is likely a more preferable choice.