Binaries play key roles in determining stellar parameters and exploring stellar evolution models. We build a catalog of 88 eclipsing binaries with spectroscopic information, taking advantage of observations from both the Large Sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) and the Palomar Transient Factory (PTF) surveys. A software pipeline is constructed to identify binary candidates by examining their light curves. The orbital periods of binaries are derived from the Lomb-Scargle method. The key distinguishing features of eclipsing binaries are recognized by a new filter \textit{Flat Test}. We classify the eclipsing binaries by applying Fourier analysis on the light curves. Among all the binary stars, 13 binaries are identified as eclipsing binaries for the first time. The catalog contains information: position, primary eclipsing magnitude and time, eclipsing depth, the number of photometry and radial velocity observations, largest radial velocity difference, binary type, the effective temperature of observable star $T_{\rm eff}$, and surface gravity of observable star log \emph{g}. The false-positive probability is calculated by using both a Monte Carlo simulation and real data from the SDSS Stripe 82 Standard Catalog. The binaries in the catalog are mostly with a period of less than one day. The period distribution shows a 0.22-day cut-off which is consistent with the low probability of an eclipsing binary rotating with such a period.

The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$\beta$ and Pa$\gamma$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$.

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

The radial velocity (RV) is a basic physical quantity which can be determined through Doppler shift of the spectrum of a star. The precision of RV measurement depends on the resolution of the spectrum we used and the accuracy of wavelength calibration. In this work, radial velocities of LAMOST-II medium resolution (R ~ 7500) spectra are measured for 1,594,956 spectra (each spectrum has two wavebands) through matching with templates. A set of RV standard stars are used to recalibrate the zero point of the measurement, and some reference sets with RVs derived from medium/high-resolution observations are used to evaluate the accuracy of the measurement. Comparing with reference sets, the accuracy of our measurement can get 0.0227 km s/1 with respect to radial velocities standard stars. The intrinsic precision is estimated with the multiple observations of single stars, which can achieve to 1.36 km s/1,1.08 km s/1, 0.91 km s/1 for the spectra at signal-to-noise levels of 10, 20, 50, respectively.

Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample consists of 9 binaries, where each source contains a giant star with large radial velocity variation ($\Delta V_{\rm R} > 70~{\rm km~s^{-1}}$) and periods known from light curves. We focus on the 9 sources with long periods ($T_{\rm ph} > 5$ days) and evaluate the mass $M_2$ of the optically invisible companion. Since the observed $\Delta V_{\rm R}$ from only a few repeating spectroscopic observations is a lower limit of the real amplitude, the real mass $M_2$ can be significantly higher than the current evaluation. It is likely an efficient method to place constraints on $M_2$ by combining $\Delta V_{\rm R}$ from LAMOST and $T_{\rm ph}$ from ASAS-SN, particularly by the ongoing LAMOST Medium Resolution Survey.

This paper reports on the application of the supervised machine-learning algorithm to the stellar effective temperature regression for the second $Gaia$ data release, based on the combination of the stars in four spectroscopic surveys: Large Sky Area Multi-Object Fiber Spectroscopic Telescope, Sloan Extension for Galactic Understanding and Exploration, the Apache Point Observatory Galactic Evolution Experiment and the RAdial Velocity Extension. This combination, about four million stars, enables us to construct one of the largest training sample for the regression, and further predict reliable stellar temperatures with a root-mean-squared error of 191 K. This result is more precise than that given by $Gaia$ second data release that is based on about sixty thousands stars. After a series of data cleaning processes, the input features that feed the regressor are carefully selected from the $Gaia$ parameters, including the colors, the 3D position and the proper motion. These $Gaia$ parameters is used to predict effective temperatures for 132,739,323 valid stars in the second $Gaia$ data release. We also present a new method for blind tests and a test for external regression without additional data. The machine-learning algorithm fed with the parameters only in one catalog provides us an effective approach to maximize sample size for prediction, and this methodology has a wide application prospect in future studies of astrophysics.

We investigate the chemical potential and baryon number density of the hadron-quark phase transition in neutron star matter. The hadron matter is described with relativistic mean field theory, and the quark matter is described with the Dyson-Schwinger equation approach of QCD. In order to study the first-order phase transition, we develop the sound speed interpolation scheme to construct the equation of state in the middle density region where the hadron phase and quark phase coexist. The phase transition chemical potential is constrained with the maximum mass, the tidal deformability and the radius of neutrons stars. And the most probable value of the phase transition chemical potential is found.

We propose a scheme to determine the chemical potential and baryon number density of the hadron-quark phase transition in cold dense strong interaction matter (compact star matter). The hadron matter is described with the relativistic mean field theory, and the quark matter is described with the Dyson-Schwinger equation approach of QCD. To study the first-order phase transition, we take the sound speed as the interpolation objective to construct the equation of state in the middle density region. With the maximum mass, the tidal deformability and the radius of neutron stars being taken as calibration quantities, the phase transition chemical potential is constrained to a quite small range. And the most probable value of the phase transition chemical potential is found.

We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $\Delta V_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods.

X-ray emission is an important indicator of stellar activity. In this paper, we study stellar X-ray activity using the XMM-Newton and LAMOST data for different types of stars. We provide a sample including 1259 X-ray emitting stars, of which 1090 have accurate stellar parameter estimations. Our sample size is much larger than those in previous works. We find a bimodal distribution of X-ray to optical flux ratio (log(fX/fV)) for G and K stars. We interpret that this bimodality is due to two subpopulations with different coronal heating rates. Furthermore, using the full widths at half maxima calculated from H{\alpha} and H\b{eta} lines, we show that these stars in the inactive peaks have smaller rotational velocities. This is consistent with the magnetic dynamo theory that stars with low rotational velocities have low levels of stellar activity. We also examine the correlation between log(fX/fV) and luminosity of the excess emission in the H{\alpha} line, and find a tight relation between the coronal and chromospheric activity indicators.

Radial velocity is one of key measurements in understanding the fundamental properties of stars, stellar clusters and the Galaxy. A plate of stars in the Kepler field were observed in May of 2018 with the medium-resolution spectrographs of LAMOST, aiming to test the performance of this new system which is the upgraded equipment of LAMOST after the first five-year regular survey.We present our analysis on the radial velocity measurements (RVs) derived from these data. The results show that slight and significant systematic errors exist among the RVs obtained from the spectra collected by different spectrographs and exposures, respectively. After correcting the systematic errors with different techniques, the precision of RVs reaches ~1.3, ~1.0, ~0.5 and ~0.3 km/s at S/Nr = 10, 20, 50, and 100, respectively. Comparing with the RVs of the standard stars of the APOGEE survey, our RVs are calibrated with a zero-point shift of ~7 km/s. The results indicate that the LAMOST medium-resolution spectroscopic system may provide RVs in a reasonable accuracy and precision for the selected targets.

We probed stellar X-ray activity over a wide range of stellar parameters, using {\it Chandra} and LAMOST data. We measured the X-ray-to-bolometric luminosity ratio ($R_X = L_X/L_{\rm bol}$) for 484 main-sequence stars, and found a bimodal distribution for G and K types. We interpret this bimodality as evidence of two sub-populations with different coronal temperatures, which are caused by different coronal heating rates. Using the metallicity and velocity information, we find both of the two sub-populations are mostly located in the thin disk. We find no trend of $R_X$ with stellar age for stars older than $\sim$ 4 Gyr; there is a trough in the $R_X$ vs age distribution, with the lowest range of $R_X$ appearing at ages around 2 Gyr. We then examined the correlation between $R_X$ and $R_{\rm H\alpha}$ (proxy of chromospheric activity): we find that the two quantities are well correlated, as found in many earlier studies. Finally, we selected a sample of twelve stars with X-ray flares, and studied the light-curve morphology of the flares. The variety of flare profiles and timescales observed in our sample suggests the contribution of different processes of energy release.

We present an empirical stellar spectra library created using spectra from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5. This library represents a uniform data set ranging from 3750 through 8500K in effective temperature (T_eff), from -2.5 through +1.0 dex in metallicity ([Fe/H]), and from 0 to 5.0 dex in gravity (log g). The spectra in the library have resolution R~1800, with well-calibrated fluxes and rest-framed wavelengths. Using a large number of red stars observed by LAMOST, we generated denser K type templates to fill in data missing from current empirical spectral libraries, particularly the late K type. For K giants, we calibrated the spectroscopic surface gravities against the asteroseismic surface gravities. To verify the reliability of the parameters labeled for this library, we performed an internal cross-validation by using a chi^2 minimization method to interpolate parameters of each individual spectrum using the remaining spectra in the library. We obtained precisions of 41 K, 0.11 dex, and 0.05 dex for T_eff, log g, and [Fe/H], respectively, which means the templates are labeled with correct stellar parameters. Through external comparisons, we confirmed that measurements of the stellar parameters through this library can achieve accuracies of approximately 125K in T_eff, 0.1 dex in [Fe/H] and 0.20 dex in log g without systematic offset. This empirical library is useful for stellar parameter measurements because it has large parameter coverage and full wavelength coverage from 3800 to 8900A.

Stars of spectral type Oe are very rare. To date, only 13 Oe stars have been identified within our Galaxy. In this paper, we present six new Oe stars and four new B0e stars found in LAMOST DR5. Repeated spectral observations of the same Oe stars show some emission line variability. The H$\beta$ emission of TYC 4801-17-1 shows rapid V/R variation. Phase lags in the V/R ratio of TYC 4801-17-1 spectra are also seen. We found the unusual O4.5 star RL 128 is an Oe star with variable H$\alpha$ intensity and its Ca II triplet emission appears when H$\alpha$ emission reaches maximum intensity. These newly identified early type Oe and B0e stars significantly increase the known sample.

We investigate the equation of state(EoS) and the effect of the hadron-quark phase transition of strong interaction matter in compact stars. The hadron matter is described with the relativistic mean field theory,and the quark matter is described with the Dyson-Schwinger equation approach of QCD. The complete EoS of the hybrid star matter is constructed with not only the Gibbs construction but also the 3-window interpolation. The mass-radius relation of hybrid stars is also investigated. We find that, although the EoSs of both the hadron matter with hyperon and $\Delta$-baryon and the quark matter are generally softer than that of the nucleon matter, the 3-window interpolation construction may provide an EoS stiff enough for a hybrid star with mass exceeding 2$M_{\odot}^{}$ and, in turn, solve the so called "hyperon puzzle".

We present a spectroscopic redshift catalog from the LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in the Southern Galactic Cap (SGC), which is designed to observe all sources (Galactic and extra-galactic) by using repeating observations with a limiting magnitude of $r=18.1~mag$ in two $20~deg^2$ fields. The project is mainly focusing on the completeness of LAMOST ExtraGAlactic Surveys (LEGAS) in the SGC, the deficiencies of source selection methods and the basic performance parameters of LAMOST telescope. In both fields, more than 95% of galaxies have been observed. A post-processing has been applied to LAMOST 1D spectrum to remove the majority of remaining sky background residuals. More than 10,000 spectra have been visually inspected to measure the redshift by using combinations of different emission/absorption features with uncertainty of $\sigma_{z}/(1+z)<0.001$. In total, there are 1528 redshifts (623 absorption and 905 emission line galaxies) in Field A and 1570 redshifts (569 absorption and 1001 emission line galaxies) in Field B have been measured. The results show that it is possible to derive redshift from low SNR galaxies with our post-processing and visual inspection. Our analysis also indicates that up to 1/4 of the input targets for a typical extra-galactic spectroscopic survey might be unreliable. The multi-wavelength data analysis shows that the majority of mid-infrared-detected absorption (91.3%) and emission line galaxies (93.3%) can be well separated by an empirical criterion of $W2-W3=2.4$. Meanwhile, a fainter sequence paralleled to the main population of galaxies has been witnessed both in $M_r$/$W2-W3$ and $M_*$/$W2-W3$ diagrams, which could be the population of luminous dwarf galaxies but contaminated by the edge-on/highly inclined galaxies ($\sim30\%$).

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

Sky subtraction is the key technique in data reduction of multi-fiber spectra. Knowledge of the related instrument character is necessary to determine the method adopted in sky subtraction. In this study, we described the sky subtraction method designed for LAMOST(Large sky Area Multi-Object fiber Spectroscopic Telescope) survey. The method has been intergrated into LAMOST 2D Pipeline v2.6 and applied to data of LAMOST DR3 and later. For LAMOST, sky emission line calibration is used to alleviate the position-dependent (thus time-dependent) ~4% fiber throughput uncertainty and the small wavelength instability (0.1\AA ) during observation. PCA (Principal Component Analysis) sky subtraction further reduces 25% of the sky line residual of the OH lines in the red part of the LAMOST spectra after the mater sky spectrum, which is derived from a B-spline fit of 20 sky fibers in each spectrograph, is adjusted by sky emission line and subtracted from each fiber. Further analysis shows that our wavelength calibration accuracy is about 4.5km/s, and the average sky subtraction residuals are about 3% for sky emission lines and 3% for continuum region. The relative sky subtraction residuals vary with the moon light background brightness, could reach as low as 1.5% for the sky emission line regions in the dark night. Tests on the F stars of both similar sky emission line strength and similar object continuum intensity show that the sky emission line residual of LAMOST is smaller than those of SDSS survey.

We present candidate members of the Pal 5, GD-1, Cetus Polar, and Orphan tidal stellar streams found in LAMOST DR3, SDSS DR9 and APOGEE catalogs. In LAMOST DR3, we find 20, 4, 24 high confidence candidates of tidal streams GD-1, Cetus Polar and Orphan respectively. We also list from the SDSS DR9 spectroscopic catalog 59, 118, 10 high confidence candidates of tidal streams Cetus Polar, Orphan and Pal 5, respectively. Furthermore, we find 7 high confidence candidates of the Pal 5 tidal stream in APOGEE data. Compared with SDSS, the new candidates from LAMOST DR3 are brighter, so that together, more of the color-magnitude diagram, including the giant branch can be explored. Analysis of SDSS data shows that there are 3 metallicity peaks of the Orphan stream and also shows some spatial separation. LAMOST data confirms multiple metallicities in this stream. The metallicity, given by the higher resolution APOGEE instrument, of the Pal 5 tidal stream is [Fe/H] $\sim -1.2$, higher than that given earlier by SDSS spectra. Many previously unidentified stream members are tabulated here for the first time, along with existing members, allowing future researchers to further constrain the orbits of these objects as they move within the Galaxy's dark matter potential.

We present the second release of value-added catalogues of the LAMOST Spectroscopic Survey of the Galactic Anticentre (LSS-GAC DR2). The catalogues present values of radial velocity $V_{\rm r}$, atmospheric parameters --- effective temperature $T_{\rm eff}$, surface gravity log$g$, metallicity [Fe/H], $\alpha$-element to iron (metal) abundance ratio [$\alpha$/Fe] ([$\alpha$/M]), elemental abundances [C/H] and [N/H], and absolute magnitudes ${\rm M}_V$ and ${\rm M}_{K_{\rm s}}$ deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC since September 2011 until June 2014. The catalogues also give values of interstellar reddening, distance and orbital parameters determined with a variety of techniques, as well as proper motions and multi-band photometry from the far-UV to the mid-IR collected from the literature and various surveys. Accuracies of radial velocities reach 5kms$^{-1}$ for late-type stars, and those of distance estimates range between 10 -- 30 per cent, depending on the spectral signal-to-noise ratios. Precisions of [Fe/H], [C/H] and [N/H] estimates reach 0.1dex, and those of [$\alpha$/Fe] and [$\alpha$/M] reach 0.05dex. The large number of stars, the contiguous sky coverage, the simple yet non-trivial target selection function and the robust estimates of stellar radial velocities and atmospheric parameters, distances and elemental abundances, make the catalogues a valuable data set to study the structure and evolution of the Galaxy, especially the solar-neighbourhood and the outer disk.

We investigate the nonradial oscillations of newly born neutron stars (NSs) and strange quark stars (SQSs). This is done with the relativistic nuclear field theory with hyperon degrees of freedom employed to describe the equation of state for the stellar matter in NSs, and with both the MIT bag model and the Nambu--Jona-Lasinio model adopted to construct the configurations of the SQSs. We find that the gravitational-mode ($g$-mode) eigenfrequencies of newly born SQSs are about one order of magnitude lower than those of NSs, which is independent of models implemented to describe the equation of state for the strange quark matter. Meanwhile the eigenfrequencies of the other modes of nonradial oscillations, e.g., fundamental ($f$)- and pressure ($p$)-modes, are much larger than those of the $g$-mode. In the light of the first direct observation of gravitational waves, it is promising to employ the gravitational waves to identify the QCD phase transition in high density strong interaction matter.

The task of flux calibration for LAMOST (Large sky Area Multi-Object Spectroscopic Telescope) spectra is difficult due to many factors. For example, the lack of standard stars, flat fielding for large field of view, and variation of reddening between different stars especially at low galactic latitudes etc. Poor selection, bad spectral quality, or extinction uncertainty of standard stars not only might induce errors to the calculated spectral response curve (SRC), but also might lead to failures in producing final 1D spectra. In this paper, we inspected spectra with Galactic latitude |b|>=60 degree and reliable stellar parameters, determined through the LAMOST Stellar Parameter Pipeline (LASP), to study the stability of the spectrograph. To guarantee the selected stars had been observed by each fiber, we selected 37,931 high quality exposures of 29,000 stars from LAMOST DR2, and more than 7 exposures for each fiber. We calculated the SRCs for each fiber for each exposure, and calculated the statistics of SRCs for spectrographs with both the fiber variations and time variations. The result shows that the average response curve of each spectrograph (henceforth ASPSRC) is relatively stable with statistical errors <= 10%. From the comparison between each ASPSRC and the SRCs for the same spectrograph obtained by 2D pipeline, we find that the ASPSRCs are good enough to use for the calibration. The ASPSRCs have been applied to spectra which were abandoned by LAMOST 2D pipeline due to the lack of standard stars, increasing the number of LAMOST spectra by 52,181 in DR2. Comparing those same targets with SDSS, the relative flux differences between SDSS spectra and that of LAMOST spectra with the ASPSRC method are less than 10%, which underlines that the ASPSRC method is feasible for LAMOST flux calibration.

Bolton and Schlegel presented a promising deconvolution method to extract 1D spectra from a 2D optical fiber spectral CCD image. The method could eliminate the PSF difference between fibers, extract spectra to the photo noise level, as well as improve the resolution. But the method is limited by its huge computation requirement and thus cannot be implemented in actual data reduction. In this article, we develop a practical computation method to solve the computation problem. The new computation method can deconvolve a 2D fiber spectral image of any size with actual PSFs, which may vary with positions. Our method does not require large amounts of memory and can extract a 4k multi 4k noise-free CCD image with 250 fibers in 2 hr. To make our method more practical, we further consider the influence of noise, which is thought to be an intrinsic illposed problem in deconvolution algorithms. We modify our method with a Tikhonov regularization item to depress the method induced noise. Compared with the results of traditional extraction methods, our method has the least residual and influence by cross talk and noise, even for extreme situation. Our method can convergent in 2 to 4 iterations, and the computation times are about 3.5 hr for the extreme fiber distance and about 2 hr for nonextreme cases. Finally, we apply our method to real LAMOST (Large sky Area Multi-Object fiber Spectroscopic Telescope, a.k.a. Guo Shou Jing Telescope) data. We find that the 1D spectra extracted by our method have both higher signal-to-noise ratio and resolution than the traditional methods.

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

We present a sample of metallic-line star (Am) candidates from the Large sky Area Multi-Object fiber Spectroscopic Telescope Data Release one (LAMOST DR1). According to the characteristic of under-abundance of calcium and overabundance of iron element of Am stars, we propose an empirical separation curve derived from line indices of Ca II K-line and iron lines we choose for low resolution spectra. 3537 Am candidates are ultimately selected from more than 30,000 stars which are classified as A-type or early-F stars by both LAMOST pipeline and visual inspection. Then we make some analysis on this sample and finally provide a list of these Am candidates with 10 relevant parameters. Comparing with other catalogues, Am candidates selected from LAMOST DR1 are much fainter on the whole. Obviously, our list is an important complementary to already known bright Am catalogues, and it offers valuable material for the research on this type of chemically peculiar stars.

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.

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.

Precise fiber positioning is crucial to a wide field, multi-fiber spectroscopic survey like LAMOST. Nowadays, most position error measurements are based on CCD photographic and imaging processing techniques. Those methods only work for measuring errors orthogonal to the telescope optical axis, while there also lies errors parallel to the telescope optical axis, like defocusing, and error caused by the existing deviation angle between optical axes of a fiber and the telescope. Directly measuring two latter types of position errors is difficult for individual fiber, especially during observation. Possible sources of fiber position errors are discussed in brief for LAMOST. By constructing a model of magnitude loss due to the fiber position error for a point source, we propose an indirect method to calculate the total and systematic position errors for each individual fiber from spectra data. Restrictions and applications of this method are also discussed.

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

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

Spectrograph slit is conventionally used to enhance the spectral resolution and manage how much light can be allowed to enter spectrograph. A narrower slit provides a higher resolution but sacrifices efficiency of spectrograph and results in a low signal to noise ratio (S/N) spectra product. We take GuoShouJing telescope as an example and carry out a series of experiments to study how its 2/3 slit mode affects the precision of stellar radial velocity measurement and atmosphere parameters estimate. By transforming the resolution and adding a Gaussian White Noise to the extremely high quality spectra from the Sloan Digital Sky Survey, we generate synthetic stellar spectra of various brightness with different S/Ns. Comparing the measurements on these noise added spectra with the original high quality ones, we summarize the influences of the 2/3 slit mode on the measurement accuracy of stellar radial velocity and atmospheric parameters.

With a sky brightness monitor in Xinglong station of National Astronomical Observatories of China (NAOC), we collected data from 22 dark clear nights and 90 lunar nights. We first measured the sky brightness variation with time in dark nights, found a clear correlation between the sky brightness and human activity. Then with a modified sky brightness model of moon night and data from moon night, we derived the typical value for several important parameters in the model. With these results, we calculated the sky brightness distribution under a given moon condition for Xinglong station. Furthermore, we simulated the moon night sky brightness distribution in a 5 degree field of view telescope (such as LAMOST). These simulations will be helpful to determine the magnitude limit, exposure time as well as the survey design for LAMOST at lunar night.

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.

The Andromeda galaxy was observed by the Guoshoujing Telescope (GSJT, formly named the Large Sky Area Multi-Object Fiber Spectroscopic Telescope -- LAMOST) during the 2009 commissioning phase. Due to the absence of standard stars for flux calibration, we use the photometric data of 15 intermediate bands in the Beijing-Arizona-Taipei-Connecticut (BATC) survey to calibrate the spectra. Total 59 spectra located in the bulge and disk of the galaxy are obtained. Kinematic and stellar population properties of the stellar content are derived with these spectra. We obtain the global velocity field and calculate corresponding rotation velocities outer to about 7 kpc along the major axis. These rotation velocity measurements complement those of the gas content, such as the H {\sc i} and CO. The radial velocity dispersion presents that the stars in the bulge are more dynamically thermal and the disk is more rotation-supported. The age distribution shows that the bulge was formed about 12 Gyr ago, the disk is relatively younger, and the ages of some regions along the spiral arms can reach as young as about 1 Gyr. These young stellar populations have relatively richer abundance and larger reddening. The overall average metallicity of the galaxy approximates the solar metallicity and a very weak abundance gradient is gained. The reddening map gives a picture of a dust-free bulge and a distinct dusty ring in the disk.

A number of spectroscopic surveys have been carried out or are planned to study the origin of the Milky Way. Their exploitation requires reliable automated methods and softwares to measure the fundamental parameters of the stars. Adopting the ULySS package, we have tested the effect of different resolutions and signal-to-noise ratios (SNR) on the measurement of the stellar atmospheric parameters (effective temperature Teff, surface gravity log g, and metallicity [Fe/H]). We show that ULySS is reliable to determine these parameters with medium-resolution spectra (R~2000). Then, we applied the method to measure the parameters of 771 stars selected in the commissioning database of the Guoshoujing Telescope (GSJT). The results were compared with the SDSS/SEGUE Stellar Parameter Pipeline (SSPP), and we derived precisions of 167 K, 0.34 dex, and 0.16 dex for Teff, log g and [Fe/H] respectively. Furthermore, 120 of these stars are selected to construct the primary stellar spectra template library (Version 1.0) of GSJT, and will be deployed as basic ingredients for the GSJT automated parametrization pipeline.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

We report the discovery of eight new quasars in one extragalactic field (five degree centered at RA=$08^h58^m08.2^s$, Dec=$01^o32'29.7''$) with the LAMOST commissioning observations on December 18, 2009. These quasars, with $i$ magnitudes from 16.44 to 19.34 and redshifts from 0.898 to 2.773, were not identified in the SDSS spectroscopic survey, though six of them with redshifts less than 2.5 were selected as quasar targets in SDSS. Except one source without near-IR $Y$-band data, seven of these eight new quasars meet a newly proposed quasar selection criterion involving both near-IR and optical colors. Two of them were found in the 'redshift desert' for quasars ($z$ from 2.2 to 3), indicating that the new criterion is efficient for recovering the missing quasars with similar optical colors as stars. Although LAMOST met some problems during the commissioning observations, we were still able to identify other 38 known SDSS quasars in this field, with $i$ magnitudes from 16.24 to 19.10 and redshifts from 0.297 to 4.512. Our identifications imply that a substantial fraction of quasars may be missing in the previous quasar surveys. The implication of our results to the future LAMOST quasar survey is discussed.

The redshift range from 2.2 to 3, is known as the 'redshift desert' of quasars because quasars with redshift in this range have similar optical colors as normal stars and are thus difficult to be found in optical sky surveys. A quasar candidate, SDSS J085543.40-001517.7, which was selected by a recently proposed criterion involving near-IR $Y-K$ and optical $g-z$ colors, was identified spectroscopically as a new quasar with redshift of 2.427 by the LAMOST commissioning observation in December 2009 and confirmed by the observation made with the NAOC/Xinglong 2.16m telescope in March 2010. This quasar was not targeted in the SDSS spectroscopic survey because it locates in the stellar locus of the optical color-color diagrams, while it is clearly separated from stars in the $Y-K$ vs. $g-z$ diagram. Comparing with other SDSS quasars we found this new quasar with $i$ magnitude of 16.44 is apparently the brightest one in the redshift range from 2.3 to 2.7. From the spectral properties we derived its central black hole mass as $(1.4\sim3.9) \times 10^{10} M_\odot$ and the bolometric luminosity as $3.7\times 10^{48}$ \ergs, which indicates that this new quasar is intrinsically very bright and belongs to the most luminous quasars in the universe. Our identification supports that quasars in the redshift desert can be found by the quasar selection criterion involving the near-IR colors. More missing quasars are expected to be recovered by the future LAMOST spectroscopic surveys, which is important to the study of the cosmological evolution of quasars at redshift higher than 2.2.