SN 2014C was an unprecedented supernova (SN) that displayed a metamorphosis from Type Ib to Type IIn over $\sim$200 days. This transformation is consistent with a helium star having exploded in a cavity surrounded by a dense shell of the progenitor's stripped hydrogen envelope. For at least 5 years post-explosion, the ejecta continued to interact with an outer, extended component of circumstellar medium (CSM) that was ejected even before the dense shell. It is still unclear, however, what kind of progenitor could have undergone such a complicated mass-loss history before it produced this peculiar SN. In this paper, we report a new analysis of SN 2014C's host star cluster based on data from the Hubble Space Telescope (HST). By carefully fitting its spectral energy distribution (SED), we derive a precise cluster age of 20.0$^{+3.5}_{-2.6}$ Myr, which corresponds to the progenitor's lifetime assuming coevolution. Combined with binary stellar evolution models, we find that SN 2014C's progenitor may have been an $\sim$11-$M_\odot$ star in a relatively wide binary system. The progenitor's envelope was partially stripped by Case C or Case BC mass transfer via binary interaction, followed by a violent eruption that ejected the last hydrogen layer before terminal explosion. Thus, SN 2014C, in common with SNe 2006jc and 2015G, may be a third example that violent eruptions, with mass-loss rates matching luminous blue variable (LBV) giant eruptions, can also occur in much lower-mass massive stars if their envelopes are partially or completely stripped in interacting binaries.

Studies of young stellar objects (YSOs) in the Galaxy have found that a significant fraction exhibit photometric variability. However, no systematic investigation has been conducted on the variability of extragalactic YSOs. Here we present the first variability study of massive YSOs in a $\sim 1.5\,\mathrm{deg^2}$ region of the Large Magellanic Cloud (LMC). The aim is to investigate whether the different environmental conditions in the metal-poor LMC ($\sim$ 0.4-0.5 Z_sun) have an impact on the variability characteristics. Multi-epoch near-infrared (NIR) photometry was obtained from the VISTA Survey of the Magellanic Clouds (VMC) and our own monitoring campaign using the VISTA telescope. By applying a reduced $\chi^2$-analysis, stellar variability was identified. We found 3062 candidate variable stars from a population of 362 425 stars detected. Based on several Spitzer studies, we compiled a sample of high-reliability massive YSOs: a total of 173 massive YSOs have NIR counterparts ($K_{\mathrm{s}}\sim 18.5\,$mag) in the VMC catalogue, of which 39 display significant ($>3\sigma$) variability. They have been classified as eruptive, fader, dipper, short-term variable and long-period variable YSOs based mostly on the appearance of their $K_{\mathrm{s}}$ band light curves. The majority of YSOs are aperiodic, only five YSOs exhibit periodic lightcurves. The observed amplitudes are comparable or smaller than those for Galactic YSOs (only two Magellanic YSOs exhibit $\Delta K_{\mathrm{s}}>1\,$mag), not what would have been expected from the typically larger mass accretion rates observed in the Magellanic Clouds.

Type Ibn supernovae (SNe Ibn) are intriguing stellar explosions whose spectra exhibit narrow helium lines with little hydrogen. They trace the presence of circumstellar material (CSM) formed via pre-SN eruptions of their stripped-envelope progenitors. Early work has generally assumed that SNe Ibn come from massive Wolf-Rayet (WR) stars via single star evolution. In this paper, we report ultraviolet (UV) and optical observations of two nearby Type Ibn SNe 2006jc and 2015G conducted with the Hubble Space Telescope (HST) at late times. A point source is detected at the position of SN 2006jc, and we confirm the conclusion of Maund et al. that it is the progenitor's binary companion. Its position on the Hertzsprung-Russell (HR) diagram corresponds to a star that has evolved off the main sequence (MS); further analysis implies a low initial mass for the companion star ($M_2$ $\le$ 11.9$^{+1.2}_{-0.8}$ $M_\odot$) and a secondary-to-primary initial mass ratio very close to unity ($q$ = $M_2/M_1$ $\sim$ 1); the SN progenitor's hydrogen envelope had been stripped through binary interaction. We do not detect the binary companion of SN 2015G. For both SNe, the surrounding stellar populations have relatively old ages and argue against any massive WR stars as their progenitors. These results suggest that SNe Ibn may have lower-mass origins in interacting binaries. As a result, they also provide evidence that the giant eruptions commonly seen in massive luminous blue variables (LBVs) can also occur in much lower-mass, stripped-envelope stars just before core collapse.

Detailed studies of intermediate/low mass pre-main sequence (PMS) stars outside the Galaxy have so far been conducted only for small targeted regions harbouring known star formation complexes. The VISTA Survey of the Magellanic Clouds (VMC) provides an opportunity to study PMS populations down to solar masses on a galaxy-wide scale. Our goal is to use near-infrared data from the VMC survey to identify and characterise PMS populations down to ~1 M_sun across the Magellanic Clouds. We present our colour-magnitude diagram method, and apply it to a ~1.5 deg^2 pilot field located in the Large Magellanic Cloud. The pilot field is divided into equally-sized grid elements. We compare the stellar population in every element with the population in nearby control fields by creating K_s/(Y-K_s) Hess diagrams; the observed density excesses over the local field population are used to classify the stellar populations. Our analysis recovers all known star formation complexes in this pilot field (N44, N51, N148 and N138) and for the first time reveals their true spatial extent. In total, around 2260 PMS candidates with ages $\lesssim$ 10 Myr are found in the pilot field. PMS structures, identified as areas with a significant density excess of PMS candidates, display a power-law distribution of the number of members with a slope of -0.86+-0.12. We find a clustering of the young stellar populations along ridges and filaments where dust emission in the far-infrared (FIR) (70 micron - 500 micron) is bright. Regions with young populations lacking massive stars show a lesser degree of clustering and are usually located in the outskirts of the star formation complexes. At short FIR wavelengths (70 micron, 100 micron) we report a strong dust emission increase in regions hosting young massive stars, which is less pronounced in regions populated only by less massive ($\lesssim$ 4 M_sun) PMS stars.

VISTA observed the Small Magellanic Cloud (SMC), as part of the VISTA survey of the Magellanic Clouds system (VMC), for six years (2010-2016). The acquired multi-epoch YJKs images have allowed us to probe the stellar populations to an exceptional level of detail across an unprecedented wide area in the near-infrared. This contribution highlights the most recent VMC results obtained on the SMC focusing, in particular, on the clustering of young stellar populations, on the proper motion of stars in the main body of the galaxy and on the spatial distribution of the star formation history.

This study expands the coverage and improves the homogeneity of the distribution of MILES template stars in the parameter space, as well as extends the wavelength coverage of the template spectra to the far red beyond the Ca\,{\sc{ii}} triplet. To achieve this we have carried out a major observational campaign using the OMR long-slit spectrograph mounted on the NAOC 2.16\,m telescope and the YFOSC long-slit spectrograph mounted on the YNAO 2.4\,m telescope. The original sample is based on the MILES library, supplemented by 918 stars selected from PASTEL database. In total, 822 OMR and 1,324 YFOSC spectra, covering respectively the wavelength ranges $\lambda\lambda$3800--5180 and $\lambda\lambda$5150--9000, have been collected and reduced. The spectra have a mean resolution FWHM (full-width at half-maximum) of $\sim 3.3$\,{\AA} and are wavelength- and flux-calibrated to an accuracy of $\sim 20$\,km\,s$^{-1}$ and $\sim 5$\,per\,cent, respectively. The spectra are further corrected for systematic errors in the wavelength calibration to an accuracy of $\sim 4$\,km\,s$^{-1}$ by cross-correlating with the theoretical spectra. Almost all the spectra have an average signal to noise ratio (SNR) better than 100 per pixel. Combined with the MILES spectra, there are now 1,731, 1,542, 1,324 and 1,273 stars with spectra covering respectively $\lambda\lambda$3800--5180, $\lambda\lambda$3800--7500, $\lambda\lambda$5150--9000 and $\lambda\lambda$3800--9000. This paper describes our template star selection, the observation and data reduction, and presents the reduced spectra collected hitherto.

We recover the spatially resolved star formation history across the entire main body and Wing of the Small Magellanic Cloud (SMC), using fourteen deep tile images from the VISTA survey of the Magellanic Clouds (VMC), in the YJKs filters. The analysis is performed on 168 subregions of size 0.143 deg2, covering a total contiguous area of 23.57 deg2. We apply a colour-magnitude diagram (CMD) reconstruction method that returns the best-fitting star formation rate SFR(t), age--metallicity relation, distance and mean reddening, together with their confidence intervals, for each subregion. With respect to previous analyses, we use a far larger set of VMC data, updated stellar models, and fit the two available CMDs (Y-Ks versus Ks and J-Ks versus Ks) independently. The results allow us to derive a more complete and more reliable picture of how the mean distances, extinction values, star formation rate, and metallicities vary across the SMC, and provide a better description of the populations that form its Bar and Wing. We conclude that the SMC has formed a total mass of (5.31+-0.05)x10^8 Msun in stars over its lifetime. About two thirds of this mass is expected to be still locked in stars and stellar remnants. 50 per cent of the mass was formed prior to an age of 6.3 Gyr, and 80 per cent was formed between 8 and 3.5 Gyr ago. We also illustrate the likely distribution of stellar ages and metallicities in different parts of the CMD, to aid the interpretation of data from future astrometric and spectroscopic surveys of the SMC.

In this paper we report a clustering analysis of upper main-sequence stars in the Small Magellanic Cloud, using data from the VMC survey (the VISTA near-infrared YJKs survey of the Magellanic system). Young stellar structures are identified as surface overdensities on a range of significance levels. They are found to be organized in a hierarchical pattern, such that larger structures at lower significance levels contain smaller ones at higher significance levels. They have very irregular morphologies, with a perimeter-area dimension of 1.44 +/- 0.02 for their projected boundaries. They have a power-law mass-size relation, power-law size/mass distributions, and a lognormal surface density distribution. We derive a projected fractal dimension of 1.48 +/- 0.03 from the mass-size relation, or of 1.4 +/- 0.1 from the size distribution, reflecting significant lumpiness of the young stellar structures. These properties are remarkably similar to those of a turbulent interstellar medium (ISM), supporting a scenario of hierarchical star formation regulated by supersonic turbulence.

We use deep multi-epoch PSF photometry taken with the Visible and Infrared Survey Telescope for Astronomy (VISTA) to measure and analyze the proper motions of stars within the Galactic globular cluster 47 Tucanae (47 Tuc, NGC 104). The observations are part of the ongoing near-infrared VISTA survey of the Magellanic Cloud system (VMC). The data analyzed here corresponds to one VMC tile and covers a total sky area of 1.77 deg^2. Absolute proper motions with respect to ~9070 background galaxies are calculated from a linear regression model applied to the stellar positions in 11 epochs in the Ks filter. The data extend over a total time baseline of ~17 months. We found a median proper motion of the stars within 47 Tuc of (mu_a cos(d), mu_d) = (+5.89 +/- 0.02(stat) +/- 0.13(sys), -2.14 +/- 0.02(stat) +/- 0.08(sys)) mas/yr, based on the measurements of ~35000 individual sources between 5' and 42' from the cluster center. We compared our result to the proper motions from the newest US Naval Observatory CCD Astrograph Catalog (UCAC5) which includes data from the Gaia data release 1. Selecting cluster members (~2700 stars) we found a median proper motion of (mu_a cos(d), mu_d) = (+5.30 +/- 0.03(stat) +/- 0.70(sys), -2.70 +/- 0.03(stat) +/- 0.70(sys)) mas/yr. The values derived from the VMC data are consistent with the UCAC5 result, and are close to literature measurements obtained using HST. We combined our results with literature radial velocity measurements and reconstructed the orbit of 47 Tuc, finding that the cluster is on a low-ellipticity orbit and is confined within the inner ~7.5 kpc of the Galaxy. We show that the use of an increased time baseline combined with PSF-determined stellar centroids in crowded regions significantly improves the accuracy of the method. In future works, we will apply this method to more VMC tiles to study in detail the kinematics of the Magellanic Clouds.

Star formation is a hierarchical process, forming young stellar structures of star clusters, associations, and complexes over a wide scale range. The star-forming complex in the bar region of the Large Magellanic Cloud is investigated with upper main-sequence stars observed by the VISTA Survey of the Magellanic Clouds. The upper main-sequence stars exhibit highly non-uniform distributions. Young stellar structures inside the complex are identified from the stellar density map as density enhancements of different significance levels. We find that these structures are hierarchically organized such that larger, lower-density structures contain one or several smaller, higher-density ones. They follow power-law size and mass distributions as well as a lognormal surface density distribution. All these results support a scenario of hierarchical star formation regulated by turbulence. The temporal evolution of young stellar structures is explored by using subsamples of upper main-sequence stars with different magnitude and age ranges. While the youngest subsample, with a median age of log($\tau$/yr)~=~7.2, contains most substructure, progressively older ones are less and less substructured. The oldest subsample, with a median age of log($\tau$/yr)~=~8.0, is almost indistinguishable from a uniform distribution on spatial scales of 30--300~pc, suggesting that the young stellar structures are completely dispersed on a timescale of $\sim$100~Myr. These results are consistent with the characteristics of the 30~Doradus complex and the entire Large Magellanic Cloud, suggesting no significant environmental effects. We further point out that the fractal dimension may be method-dependent for stellar samples with significant age spreads.

We study the luminosity function of intermediate-age red-clump stars using deep, near-infrared photometric data covering $\sim$ 20 deg$^2$ located throughout the central part of the Small Magellanic Cloud (SMC), comprising the main body and the galaxy's eastern wing, based on observations obtained with the VISTA Survey of the Magellanic Clouds (VMC). We identified regions which show a foreground population ($\sim$11.8 $\pm$ 2.0 kpc in front of the main body) in the form of a distance bimodality in the red-clump distribution. The most likely explanation for the origin of this feature is tidal stripping from the SMC rather than the extended stellar haloes of the Magellanic Clouds and/or tidally stripped stars from the Large Magellanic Cloud. The homogeneous and continuous VMC data trace this feature in the direction of the Magellanic Bridge and, particularly, identify (for the first time) the inner region ($\sim$ 2 -- 2.5 kpc from the centre) from where the signatures of interactions start becoming evident. This result provides observational evidence of the formation of the Magellanic Bridge from tidally stripped material from the SMC.

We study the hierarchical stellar structures in a $\sim$1.5 deg$^2$ area covering the 30 Doradus-N158-N159-N160 star-forming complex with the VISTA Survey of the Magellanic Clouds. Based on the young upper main-sequence stars, we find that the surface densities cover a wide range of values, from log($\Sigma\cdot$pc$^2$) $\lesssim$ $-$2.0 to log($\Sigma\cdot$pc$^2$) $\gtrsim$ 0.0. Their distributions are highly non-uniform, showing groups that frequently have sub-groups inside. The sizes of the stellar groups do not exhibit characteristic values, and range continuously from several parsecs to more than 100 pc; the cumulative size distribution can be well described by a single power law, with the power-law index indicating a projected fractal dimension $D_2$ = 1.6 $\pm$ 0.3. We suggest that the phenomena revealed here support a scenario of hierarchical star formation. Comparisons with other star-forming regions and galaxies are also discussed.

Asteroseismology allows for deriving precise values of surface gravity of stars. The accurate asteroseismic determinations now available for large number of stars in the Kepler fields can be used to check and calibrate surface gravities that are currently being obtained spectroscopically for a huge numbers of stars targeted by large-scale spectroscopic surveys, such as the on-going Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Galactic survey. The LAMOST spectral surveys have obtained a large number of stellar spectra in the Kepler fields. Stellar atmospheric parameters of those stars have been determined with the LAMOST Stellar Parameter Pipeline at Peking University (LSP3), by template matching with the MILES empirical spectral library. In the current work, we compare surface gravities yielded by LSP3 with those of two asteroseismic samples - the largest sample from Huber et al. (2014) and the most accurate sample from Hekker et al. (2012, 2013). We find that LSP3 surface gravities are in good agreement with asteroseismic values of Hekker et al. (2012, 2013), with a dispersion of about 0.2 dex. Except for a few cases, asteroseismic surface gravities of Huber et al. (2014) and LSP3 spectroscopic values agree for a wide range of surface gravities. However, some patterns of differences can be identified upon close inspection. Potential ways to further improve the LSP3 spectroscopic estimation of stellar atmospheric parameters in the near future are briefly discussed. The effects of effective temperature and metallicity on asteroseismic determinations of surface gravities for giant stars are also discussed.

We present a catalog of 908 objects observed with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) in the vicinity fields of M31 and M33, targeted as globular clusters (GCs) and candidates. The targets include known GCs and candidates selected from the literature, as well as new candidates selected from the Sloan Digital Sky Survey (SDSS). Analysis shows that 356 of them are likely GCs of various degree of confidence, while the remaining ones turn out to be background galaxies and quasars, stars and HII regions in M31 or foreground Galactic stars. The 356 likely GCs include 298 bona fide GCs and 26 candidates known in the literature. Three candidates selected from the Revised Bologna Catalog of M31 GCs and candidates (RBC) and one possible cluster from Johnson et al. are confirmed to be bona fide clusters. We search for new GCs in the halo of the M31 amongst the new candidates selected from the SDSS photometry. Based on radial velocities yielded by LAMSOT spectra and visual examination of the SDSS images, we find 28 objects, 5 bona fide and 23 likely GCs. Amongst the five bona fide GCs, three have been recently discovered independently by others, the remaining 25 are our new identifications, including two bona fide ones. The new identified objects fall at projected distances ranging from 13 to 265 kpc from M31. Of the two newly discovered bona fide GCs, one is located near M33, probably a GC belonging to M33. The other bona fide GC falls on the Giant Stream with a projected distance of 78 kpc from M31. Of the 23 newly identified likely GCs, one has a projected distance of about 265 kpc from M31 and could be an intergalactic cluster.

We use 297 042 main sequence turn-off stars selected from the LSS-GAC to determine the radial and vertical gradients of stellar metallicity of the Galactic disk in the anti-center direction. We determine ages of those turn-off stars by isochrone fitting and measure the temporal variations of metallicity gradients. Our results show that the gradients, both in the radial and vertical directions, exhibit significant spatial and temporal variations. The radial gradients yielded by stars of oldest ages (>11 Gyr) are essentially zero at all heights from the disk midplane, while those given by younger stars are always negative. The vertical gradients deduced from stars of oldest ages (>11Gyr) are negative and show only very weak variations with the Galactocentric distance in the disk plane, $R$, while those yielded by younger stars show strong variations with $R$. After being essentially flat at the earliest epochs of disk formation, the radial gradients steepen as age decreases, reaching a maxima (steepest) at age 7-8 Gyr, and then they flatten again. Similar temporal trends are also found for the vertical gradients. We infer that the assemblage of the Milky Way disk may have experienced at least two distinct phases. The earlier phase is probably related to a slow, pressure-supported collapse of gas, when the gas settles down to the disk mainly in the vertical direction. In the later phase, there are significant radial flows of gas in the disk, and the rate of gas inflow near the solar neighborhood reaches a maximum around a lookback time of 7-8 Gyr. The transition of the two phases occurs around a lookback time between 8 and 11 Gyr. The two phases may be responsible for the formation of the Milky Way thick and thin disks, respectively. And, as a consequence, we recommend that stellar age is a natural, physical criterion to distinguish thin and thick disk stars. ... (abridged)

Using a sample of over 70, 000 red clump (RC) stars with $5$-$10$% distance accuracy selected from the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC), we study the radial and vertical gradients of the Galactic disk(s) mainly in the anti-center direction, covering a significant disk volume of projected Galactocentric radius $7 \leq R_{\rm GC} \leq 14$ kpc and height from the Galactic midplane $0 \leq |Z| \leq 3$ kpc. Our analysis shows that both the radial and vertical metallicity gradients are negative across much of the disk volume probed, and exhibit significant spatial variations. Near the solar circle ($7 \leq R_{\rm GC} \leq 11.5$ kpc), the radial gradient has a moderately steep, negative slope of $-0.08$ dex kpc$^{-1}$ near the midplane ($|Z| < 0.1$ kpc), and the slope flattens with increasing $|Z|$. In the outer disk ($11.5 < R_{\rm GC} \leq 14$ kpc), the radial gradients have an essentially constant, much less steep slope of $-0.01$ dex kpc $^{-1}$ at all heights above the plane, suggesting that the outer disk may have experienced an evolution path different from that of the inner disk. The vertical gradients are found to flatten largely with increasing $R_{\rm GC}$. However, the vertical gradient of the lower disk ($0 \leq |Z| \leq 1$ kpc) is found to flatten with $R_{\rm GC}$ quicker than that of the upper disk ($1 < |Z| \leq 3$ kpc). Our results should provide strong constraints on the theory of disk formation and evolution, as well as the underlying physical processes that shape the disk (e.g. gas flows, radial migration, internal and external perturbations).

We report a detailed investigation of the bulk motions of the nearby Galactic stellar disk, based on three samples selected from the LSS-GAC DR2: a global sample containing 0.57 million FGK dwarfs out to $\sim$ 2 kpc, a local subset of the global sample consisting $\sim$ 5,400 stars within 150 pc, and an anti-center sample containing $\sim$ 4,400 AFGK dwarfs and red clump stars within windows of a few degree wide centered on the Galactic anti-center. The global sample is used to construct a three-dimensional map of bulk motions of the Galactic disk from the solar vicinity out to $\sim$ 2 kpc with a spatial resolution of $\sim$ 250 pc. Typical values of the radial and vertical components of bulk motion range from $-$15 km s$^{-1}$ to 15 km s$^{-1}$, while the lag behind the circular speed dominates the azimuthal component by up to $\sim$ 15 km s$^{-1}$. The map reveals spatially coherent, kpc-scale stellar flows in the disk, with typical velocities of a few tens km s$^{-1}$. Bending- and breathing-mode perturbations are clearly visible, and vary smoothly across the disk plane. Our data also reveal higher-order perturbations, such as breaks and ripples, in the profiles of vertical motion versus height. From the local sample, we find that stars of different populations exhibit very different patterns of bulk motion. Finally, the anti-center sample reveals a number of peaks in stellar number density in the line-of-sight velocity versus distance distribution, with the nearer ones apparently related to the known moving groups. The "velocity bifurcation" reported by Liu et al. (2012) at Galactocentric radii 10--11 kpc is confirmed. However, just beyond this distance, our data also reveal a new triple-peaked structure.

We have developed and implemented an iterative algorithm of flux calibration for the LAMOST Spectroscopic Survey of the Galactic anti-center (LSS-GAC). For a given LSS-GAC plate, the spectra are first processed with a set of nominal spectral response curves (SRCs) and used to derive initial stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$ and metallicity [Fe/H]) as well as dust reddening $E(B-V)$ of all targeted stars. For each of the sixteen spectrographs, several F-type stars of good signal-to-noise ratios (SNRs) are then selected as flux standard stars for further, iterative spectral flux calibration. Comparison of spectrophotometric colours, deduced from the flux-calibrated spectra, with the photometric measurements yields average differences of 0.02$\pm$0.07 and $-$0.04$\pm$0.09\,mag for the $(g-r)$ and $(g-i)$, respectively. The relatively large negative offset in $(g-i)$ is due to the fact that we have opted not to correct for the telluric bands, most notably the atmospheric A-band in the wavelength range of $i$-band. Comparison of LSS-GAC multi-epoch observations of duplicate targets indicates that the algorithm has achieved an accuracy of about 10 per cent in relative flux calibration for the wavelength range 4000 -- 9000\,\AA. The shapes of SRC deduced for the individual LAMOST spectrographs are found to vary by up to 30 per cent for a given night, and larger for different nights, indicating that the derivation of SRCs for the individual plates is essential in order to achieve accurate flux calibration for the LAMOST spectra.

As a major component of the LAMOST Galactic surveys, the LAMOST Spectroscopic Survey of the Galactic Anti-centre (LSS-GAC) aims to survey a significant volume of the Galactic thin/thick discs and halo for a contiguous sky area of over 3,400 deg$^2$ centred on the Galactic anti-centre ($|b| \leq 30^{\circ}$, $150 \leq l \leq 210^{\circ}$), and obtain $\lambda\lambda$3700 -- 9000 low resolution ($R \sim 1,800$) spectra for a statistically complete sample of $\sim 3$\,M stars of all colours down to a limiting magnitued of $r$ $\sim$ 17.8\,mag (to 18.5\,mag for limited fields). Together with Gaia, the LSS-GAC will yield a unique dataset to advance our understanding of the structure and assemblage history of the Galaxy, in particular its disk(s). In addition to the main survey, the LSS-GAC will also target hundreds of thousands objects in the vicinity fields of M\,31 and M\,33 and survey a significant fraction (over a million) of randomly selected very bright stars (VB; $r \le 14$ mag) in the northern hemisphere. During the Pilot and the first year Regular Surveys of LAMOST, a total of 1,042,586 [750,867] spectra of a signal to noise ratio S/N(7450\AA) $\ge$ 10 [S/N(4650\AA) $\ge$ 10] have been collected. In this paper, we present a detailed description of the target selection algorithm, survey design, observations and the first data release of value-added catalogues (including radial velocities, effective temperatures, surface gravities, metallicities, values of interstellar extinction, distances, proper motions and orbital parameters) of the LSS-GAC.

We introduce the LAMOST Stellar Parameter Pipeline at Peking University --- LSP3, developed and implemented for the determinations of radial velocity $V_{\rm r}$ and stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$, metallicity [Fe/H]) for the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC). We describe the algorithms of LSP3 and examine the accuracy of parameters yielded by it. The precision and accuracy of parameters yielded are investigated by comparing results of multi-epoch observations and of candidate members of open and globular clusters, with photometric calibration, as well as with independent determinations available from a number of external databases, including the PASTEL archive, the APOGEE, SDSS and RAVE surveys, as well as those released in the LAMOST DR1. The uncertainties of LSP3 parameters are characterized and quantified as a function of the spectral signal-to-noise ratio (SNR) and stellar atmospheric parameters. We conclude that the current implementation of LSP3 has achieved an accuracy of 5.0\,km\,s$^{-1}$, 150\,K, 0.25\,dex, 0.15\,dex for the radial velocity, effective temperature, surface gravity and metallicity, respectively, for LSS-GAC spectra of FGK stars of SNRs per pixel higher than 10. The LSP3 has been applied to over a million LSS-GAC spectra collected hitherto. Stellar parameters yielded by the LSP3 will be released to the general public following the data policy of LAMOST, together with estimates of the interstellar extinction $E(B-V)$ and stellar distances, deduced by combining spectroscopic and multi-band photometric measurements using a variety of techniques.